Abstract
Long covid follows 10-20% of first-time SARS-CoV-2 infections, but the societal burden of long covid and risk factors for the condition are not well-understood. Here, we report findings about self-reported sick leave and risk factors thereof from a hybrid survey and register study, which included 37,482 RT- PCR confirmed SARS-CoV-2 cases and 51,336 test-negative controls who were tested during the index and alpha waves. An additional 33 individuals per 1000 took substantial sick leave following acute infection compared to persons with no known history of infection, where substantial sick leave was defined as >1 month of sick leave within the period 1-9 months after the RT-PCR test date. Being female, ≥50 years, and having certain pre-existing conditions such as fibromyalgia increased risks for taking substantial sick leave. Further research exploring this heterogeneity is urgently needed and may provide important evidence for more targeted preventative strategies.
Introduction
Long covid has been termed to describe the persistence of wide-ranging, post-acute symptoms in individuals previously infected with SARS-CoV-2 (1). It has been estimated that 10-20% of people infected with SARS-CoV-2 develop long covid (2). With studies pointing to varying symptom clusters (3–6), unclear pathogenesis (1, 7), and a limited understanding of risk factors (8–11), long covid remains enigmatic. In particular, it is still unclear how many previously infected suffer from long covid to such an extent that it is severely disruptive to daily living and quality of life, and in turn also disruptive to society (12).
In Denmark, an earlier study on post-acute symptoms identified sick leave as a potential indicator of the burden of long covid symptoms, where both full- and part-time sick leave were more frequent among test-positives for SARS-CoV-2 compared to test-negatives (13). Still, the extent to which long covid symptoms translate to working ability remains understudied, and current literature on the subject has largely lacked control groups with no history of SARS-CoV-2 infection (5, 12, 14–17).
As many people have been infected and continue to be infected with SARS-CoV-2, further research is needed to understand the individual- and societal burden of long covid. Moreover, it is important to investigate the degree to which this burden is impacted by individual risk factors, such as age, sex, and preexisting health conditions. Deeper understanding of these risk factors can help guide clinical practice and target preventative public health measures.
The aim of this study was to evaluate the association between Covid-19 and post-acute sick leave and explore the possible impact of age, sex, and preexisting health conditions on this association.
Methods
Study context
Denmark’s universal SARS-CoV-2 testing strategy
In this study, it is important to consider Denmark’s national SARS-CoV-2 testing strategy and how this has facilitated the conduct of population-level studies with test-negative control groups, particularly before Omicron and the availability of home rapid antigen tests. In Denmark, universal testing for SARS- CoV-2 was implemented from the end of May 2020 and continued throughout the testing period under study from November 2020 to February 2021. Reverse transcription polymerase chain reaction (RT- PCR) tests were available and accessible for all adults free of charge and independent of the indication for acquiring a test (18). Additionally, persons admitted to hospitals were tested for SARS-CoV-2. During the testing period under study between November 2020 and February 2021 (index and alpha waves), the weekly PCR test incidence in Denmark ranged from 7,900 tests - 14,700 tests per 100,000 inhabitants (mean of 10,800 tests per 100,000 inhabitants) (19). Since February 2020, all RT-PCR test results have been registered in the Danish microbiology database (MiBa) (20).
Sick leave in Denmark
Most employees in Denmark are covered under the Act on Salaried Employees and are thus entitled to full pay during sick leave (21). Employers typically cover the first four weeks of sick leave, after which municipalities bear the costs.
Study design
In this cohort study, we merged nationwide survey- and register data. In Denmark, all residents are assigned a unique identifier (the CPR-number) in the Danish Civil Registration System, and this number is commonly used to link individual-level data from varying data sources. A summary of the characteristics of our study population and which variables we used from survey responses and register data are presented in Figure 1.
The EFTER-COVID survey
To investigate self-reported sick leave after infection with SARS-CoV-2, we used data from a nationwide Danish survey, EFTER-COVID (AFTER COVID). This survey was launched in August 2021 to investigate the general health of the Danish population during the pandemic with a particular focus on long covid. Based on RT-PCR test results recorded in MiBa, individuals were sent an invitation to participate in the study via the national digital mail system. This system (“e-Boks”) enables secure electronic communication with public authorities and is used by over 90% of all Danish residents aged ≥15 years (22). EFTER-COVID survey data were collected using Danish- or English-language web- based questionnaires created in SurveyXact (23), which could be filled out using a PC, tablet, or smartphone.
All Danish residents who had an e-Boks account and obtained a first positive RT-PCR test result for SARS-CoV-2 registered in MiBa during the period from November 2020 to February 2021 were invited to participate in the EFTER-COVID survey. Additionally, test-negative controls were randomly selected using incidence density sampling on the test date with a ratio of 2:3 between test-positive and -negative persons. This ratio was chosen to compensate for a lower expected response rate among controls compared to cases. Importantly, the test-negative controls did not have a registered positive test result in MiBa at any time point prior to receiving or responding to the questionnaire, which we could reassure due to the extensive national testing strategy. In this study, we included participants who responded to a retrospective questionnaire 9 months after their test date. Other studies which have used the EFTER- COVID survey data (13, 24, 25) include participants who responded at other points in time after testing and thus reflect different subsets of EFTER-COVID data, which contain more than 840,000 survey participants (26).
Data sources
Using the CPR-number, we enriched the EFTER-COVID questionnaire data with register-based information on age and sex, Covid-19 vaccinations from the Danish Vaccination Register (DDV), SARS- CoV-2 test results and (re)infection history registered in MiBa, and comorbidities from the Danish National Patient Register (DNPR) five years prior to each participant’s test date.
The questionnaire included questions about baseline characteristics of the participant, including height, weight, education, smoking habits, alcohol consumption, and health conditions preceding the individual’s test date. In addition, participants were asked about the amount of full- or part-time sick leave they took between their test date (indicated in the questionnaire), and the day they responded to the survey. Test-negatives were asked whether they suspected ever having had Covid-19, e.g., if they had received a seropositive test result. All questions required a response to complete the questionnaire, except for height, weight, smoking and alcohol consumption.
From the DNPR, we obtained information on in- and outpatient diagnoses coded using the 10th revision of the International Statistical Classification of Diseases and Related Problems (ICD-10), which enabled the calculation of Charlson Comorbidity Index. We also extracted information on the long covid ICD-10 code (ICD-10 code B948A) from the DNPR. A complete description of how we categorized these variables can be found in Table S1.
Exclusion criteria
Participants who did not complete the questionnaire were excluded. Furthermore, we did not include individuals who indicated that they believed they previously had SARS-CoV-2 due to receiving a seropositive result for SARS-CoV-2. Participants who were >65 years were also excluded due to retirement age, where age was calculated on the test date. See Figure S1 for a detailed flowchart of our inclusion and exclusion criteria.
Outcomes
All participants, regardless of test status, were asked whether they took sick leave around the time of their test date or at any time point after the test. Individuals who responded “yes” to taking sick leave more than four weeks after the test were then asked whether their sick leave was full-time, part-time, or both and for how long they were on full- and/or part-time sick leave.
A binary outcome variable was defined as having taken “no or ≤4 weeks of full-time sick leave >4 weeks after the test date” or “>4 weeks of full-time sick leave >4 weeks after the test date”. The latter was treated as indicative of substantial sick leave, i.e., at least one full month of sick leave was taken over an eight-month period one month after the test. An identical outcome was defined for part-time sick leave.
We chose to examine sick leave from one month after the test to capture sick leave after the acute infection period among test-positives. The self-reported sick leave was not necessarily taken consecutively such that it could capture multiple periods of sick leave due to fluctuating symptoms.
Participants who reported taking only sick leave during the week up to the test and up to four weeks after the test where included in the category “no or <4 weeks of sick leave >4 weeks after the test date”. For an illustration of how we defined our sick leave outcome, see Figure 1, Panel A.
Statistical methods
The prevalence of substantial sick leave among test-positives and -negatives were compared using risk differences (RDs), which give the difference between the risk of an outcome in the exposed group and the unexposed group. Parametric g-computation on logistic regression (27) was used to estimates RDs with 95% confidence intervals obtained using bootstrap random resampling with 1000 iterations comparing test-positive and -negative individuals with adjustments for age, sex, Charlson Comorbidity Index, and education level. RDs for the risk of substantial sick leave are expressed in percentage points. P-values for the association between test status and the survey participant characteristics were estimated using student’s t-test for continuous variables and Pearson’s Chi-squared test for categorical variables.
To investigate possible risk groups for substantial sick leave following infection with SARS-CoV-2, we conducted analyses on sub-populations defined by possible risk factors. These risk factors were defined apriori based on available variables in the survey. These included middle to older age (categorized as ≥50 years), female sex, obesity, diabetes, asthma, high blood pressure, COPD or other chronic lung disease, chronic or frequent headaches/migraines, and the following health conditions diagnosed by a medical doctor before the test: depression, anxiety, post-traumatic stress disorder, chronic fatigue syndrome, and fibromyalgia. In addition, we restricted analyses by educational level and healthcare workers. RDs with 95% confidence intervals were estimated as described above.
All statistical analyses were carried out in R version 4.1.3 (28). The R-packages “riskCommunicator” (29) was used for modelling and “forestploter” (30) for data visualization (forest plots).
Results
Overview of study population
Out of 294,035 invited, a total of 106,917 persons fully completed a questionnaire nine months after testing positive or negative for SARS-CoV-2 (response rate 36.4%). After all exclusion criteria were applied, the study population consisted of 88,818 individuals, of which 37,482 had had SARS-CoV-2 infection confirmed with a positive reverse transcription polymerase chain reaction (RT-PCR) test (Figure S1). Of all participants, 64.3% were female, and the mean age was 45 years (SD 13.8) (Table 1). Less than 1% of participants had received one or more doses of a vaccine against Covid-19. Based on self-reported height and weight, 16.6% of the population was identified as obese. Apart from obesity, the most frequent preexisting clinical characteristics were depression, high blood pressure, and anxiety (Table 1, Figure 1, Panel C).
Prevalence of substantial sick leave
The prevalence of substantial sick leave was 1.4% among test-negatives compared to 4.5% among test- positives. Among test-negatives only, the prevalence of substantial sick leave was similar across age groups. Conversely, for test-positives, the prevalence of substantial sick leave increased with age (Figure 2, Panel C). Notably, out of all test-positives who took substantial sick leave, 354 individuals (21.1%) had received a diagnosis with sequelae of SARS-CoV-2 (ICD-10 code B948A) (Figure 2, Panel A) compared to only 556 (1.6%) of the test-positives who did not have substantial sick leave.
In view of possible risk groups identified in previous literature (1, 10, 14, 31) we stratified on sex, middle to older age (>50 years), and pre-existing health conditions. Across strata, the prevalence of substantial sick leave was higher among test-positives compared to test-negatives. Among persons with pre-existing conditions, the baseline prevalence of substantial sick leave, i.e., among test-negatives, was higher than for the test-negative general population (1.4%), and the highest background prevalence was among individuals with chronic fatigue syndrome (4.1%) (Figure 3).
Risk differences for substantial sick leave
Nine months after testing for SARS-CoV-2 during the index and alpha waves, persons infected with SARS-CoV-2 had a higher risk (RD 3.3, 95% CI 3.0 to 3.5) of taking substantial sick leave after their acute infection (>1 month of sick leave within 1-9 months after the test date) compared to test-negatives with no known history of SARS-CoV-2 infection (Figure 3). Changing the definition of substantial sick leave by increasing the duration resulted in attenuation of the RD (e.g., RD = 0.5, 95% CI 0.4 to 0.6) for substantial sick leave defined as at least 6 months (Figure S3). A distributional plot underlying substantial sick leave by age (<50 years, ≥50 years) and sex is available in Figure S4.
Risk differences were greater for females (RD 3.9, 95% CI 3.6 to 4.2) than males (RD 2.1, 95% CI 1.8 to 2.5). For middle-older age (≥50 years) and all pre-existing conditions except chronic fatigue syndrome and anxiety, estimated RDs were higher than that in the general population. The largest risk differences were observed for persons with fibromyalgia (RD 10.6, 95% CI 7 to 14.6), COPD or other lung disease (RD 6.8, 95% CI 4 to 10) and diabetes (RD 5.9, 95% CI 4.2 to 7.7). Obesity, the most frequent clinical characteristic, had a larger risk difference (RD 5.7, 95% CI 4.9 to 6.4) than that in the general population (RD 3.3, 95% CI 3 to 3.5) (Figure 3). RDs by education level showed that individuals with a “higher education of 2-4 years (e.g., nurse, preschool teacher, bachelor of engineering)” had a larger risk difference (RD 3.8, 95% CI 3.4 to 4.2) than that in individuals with a “higher education of ≥5 years (e.g., master’s degree or PhD)” (RD 1.9, 95% CI 1.4 to 2.4) (Table S2). Healthcare workers (N=13,872) also had a larger risk difference (RD 4.6, 95% CI 3.9 to 5.3) than the general population.
Risk differences for part-time sick leave
Furthermore, persons infected with SARS-CoV-2 had a higher risk (RD 2.1, 95% CI 2.0 to 2.3) of substantial part-time sick leave compared to test-negatives. Similar results when stratifying on sex were observed for both part-time and full-time sick leave. However, with regards to individuals ≥50 years, a higher risk was noted for full-time sick leave than for part-time sick leave. For full-time sick leave, this age group had a higher RD (RD 4.3, 95% CI 3.9 to 4.7) than that of the full population (RD 3.3, 95% CI 3.0 to 3.5). For part-time sick leave, individuals ≥50 years had a similar risk (RD 2.3, 95% CI 2.1 to 2.6) compared to that of the full study population (RD 2.1, 95% CI 2.0 to 2.3) (Figure S2).
Discussion
Key findings
In this study, we used sick leave as a proxy for long covid and explored sick leave among a range of possible long covid risk groups. We found that individuals infected with SARS-CoV-2 during the index and alpha waves had a significantly higher risk of substantial sick leave (3.3 percentage point increase in risk), where substantial sick leave was defined as >1 month of self-reported sick leave within the period 1-9 months after the test. Notably, 20.1% of the test-positives who took substantial sick leave also had a registered ICD-10 code for long covid. Patients who receive a long covid diagnosis represent only a small proportion of individuals suffering from post-acute symptoms; as such, the fact that one in five of those experiencing substantial sick leave also have a long covid diagnosis supports the validity of using substantial sick leave as a proxy for long covid burden. We identified females, persons ≥ 50 years, obesity and a variety of pre-existing health conditions as possible risk factors, where the risk difference between testing positive and negative was greater than that of the overall study population.
Other risk factors
Other work has examined sick leave in persons previously infected with SARS-CoV-2; however, differences in risk measures, timing of measurements, national testing strategies, and definitions of sick leave vary between studies and should be carefully considered. A Danish register-based cohort study with 7,466 participants examined return-to-work following first-time infections occurring between January and May 2020, where 81.9% returned to work within four weeks of their first positive SARS- CoV-2 test (32). Although this study did not have a test-negative control group, the authors examined the cumulative incidence of return to work between patients admitted to hospitals with SARS-CoV-2 and patients admitted with influenza, which suggested that patients with SARS-CoV-2 had a reduced chance of returning to work compared to patients admitted with influenza. Additionally, a Danish cross-sectional study examined the influence of long covid on activities of daily living among 448 long covid patients reported that 56% needed sick leave and 94% were referred to rehabilitation (17). In a German register- based study with 30,950 individuals diagnosed with SARS-CoV-2, a reported 5.8% of individuals took more than four weeks of sick leave between March 2020 and February 2021 (14). These studies, together with ours, point to sick leave as a suitable indicator of the personal and societal burden of long covid. Intriguingly, what may be more telling is the heterogeneity of sick leave following SARS-CoV-2 for the purpose of identifying long covid risk groups.
Female sex has been identified as a possible risk factor for long covid (7). Some studies have indicated that females needed longer sick leave following infection with SARS-CoV-2 than their male counterparts (14, 32). Others have pointed to interactions between age and sex, e.g., less sick leave in infected females ages 20-44 compared to infected females ages 45-70 years (33). Additionally, a Swedish registry-based cohort study reported that people with recurrent sick leave were older, more often female, and more likely to have been on sick leave prior to the pandemic (16). In our study, we observed that the risk difference for substantial sick leave between SARS-CoV-2 test-positives and -negatives was greater for females (RD 3.9, 95% CI 3.6 to 4.2) than for males (RD 2.1, 95% CI 1.8 to 2.5), and that females as well as persons ≥50 years had a slightly higher prevalence of sick leave compared to the general population, irrespective of test status. Importantly, severity of SARS-CoV-2 infection is associated with increased age, pre-existing medical conditions, and male sex (34), and severity of infection has been associated with long covid (7). It is therefore unclear as to why females appear to be a risk group for substantial sick leave following acute SARS-CoV-2 infection but males are not. Further research examining age-sex interactions and the impact of multimorbidity on long covid outcomes are needed (35).
A German register-based study by Jacob et al. (14) also explored a range of possible risk factors, some of which overlapped with risk factors that we explored. As in our study, the authors also found that diabetes and high blood pressure were positively and significantly associated with long-term sick leave. However, in contrast to our study findings, the authors did not find obesity and asthma to increase the risk of long-term sick leave. Furthermore, the authors reported that anxiety and somatoform disorders were not associated with the risk of longer-term sick leave. Similarly, in our study, neither anxiety nor chronic fatigue syndrome, which is categorized as a somatoform disorder by some, increased the risk of substantial sick leave; however, new-onset chronic fatigue syndrome is a known complication of long covid (7).
Fibromyalgia, COPD/other chronic lung diseases, and obesity have also been cited as possible risk factors for long covid in other work. A US study (preprint) (N=89,843) found that long covid patients were more likely to have a history of fibromyalgia (OR 2.3, 95% CI 1.3 to 3.8) and chronic pulmonary lung disease (OR 1.9, 95% CI 1.5 to 2.6) compared to matched test-positive controls without long covid (36). Additionally, obesity has been identified as a risk factor for severe acute SARS-CoV-2 infection (1), and inpatient care following infection has been suggested to predict longer sick leave (31).
Industry
Finally, individuals’ line of work can impact the need for long covid-related sick leave. Although we were not able to examine most individuals’ professions, we did observe that individuals with a higher education of 2-4 years had a larger risk difference of taking substantial sick leave compared to individuals with a postgraduate education. We posit that these differences could be attributed to job adaptability to work-from-home and self-pacing, i.e., persons with a post-graduate education may work desk jobs whereas persons with a higher education of 2-4 years (e.g., nurse, preschool teacher, bachelor of engineering) may need to be on-site. Similar findings were reported in a Norwegian study, which estimated the industry-specific impact of the Omicron wave on sick leave compared to corresponding months from 2017-2020 (37). Persons within the food and accommodation industry had the highest increase in sick leave (4.4 percentage points increase, 95% CI 4.3-4.5), suggesting that individuals with jobs not-suited for work from home required more sick leave.
Strengths and limitations
The present study describes population-level data on post-acute sick leave as a possible consequence of infection with SARS-CoV-2. In contrast to most existing literature on the subject of long covid and working ability, this study used a test-negative control group, allowing us to consider background prevalence of sick leave for the study population and among possible risk groups, including people with various conditions which preceded their SARS-CoV-2 test. Furthermore, our sick leave outcome variable captures fluctuating illness, which is useful given the irregularity of long covid symptoms.
We consider the main limitations of this study to be its self-reporting nature and potential participation and recall bias. One shortcoming of the self-reported sick leave variable is that the duration of sick leave cannot be directly attributed to long covid symptoms, as the questionnaire asked about sick leave in general rather than sick leave associated with SARS-CoV-2 infection. A possible alternative to studying self-reported sick leave was using register-based sick leave benefits. However, the prevalence of register- based sick leave benefits was previously examined for infections occurring between January and May 2020 in Denmark (32), where the authors used data from the Danish Register for Evaluation of Marginalization (DREAM). DREAM primarily covers sick leave benefits which are granted after 30 days of sickness absence. A disadvantage of using this data is that it largely excludes shorter periods of absence (38) and thus is not well-suited for capturing fluctuating illness as a possible consequence of long covid.
Participation bias may have occurred, where individuals living with poor health or long covid symptoms may have taken more interest in participating. Alternatively, some individuals living with long covid symptoms may have felt too poorly to participate. In addition, the retrospective study design is vulnerable to recall bias, where some participants may not remember how much sick leave they took over the 9 months following their test. We sought to reduce this by multiple choice between pre-defined sick leave durations instead of using free-text. While the RD was attenuated for increasing durations of sick leave as expected, the RD for at least six months (RD 0.5, 95% CI 0.4 to 0.6) is still quite striking considering the large number of infected globally. Lastly, our results capture index- and alpha-variant infections which largely occurred before SARS-CoV-2 vaccine rollout and we cannot exclude that the absolute magnitude of our results are attenuated by vaccinations and the omicron variant.
Perspectives
Long covid is a concerning, contemporary public health issue with many unanswered questions, particularly regarding manifestations of symptom burden and risk factors. This study provides much- needed information on post-acute sick leave following SARS-CoV-2 infection in a general population and the risk of post-acute sick leave for people with comorbid conditions. The results from this study may be particularly useful to public health stakeholders in guiding evidence-based decisions concerning targeted preventative strategies. Further long covid research and funding initiatives are critical, particularly to increase knowledge of immunopathogenic mechanisms, phenotypes, risk factors, and the impact of different (sub)variants and SARS-CoV-2 vaccines.
Conclusion
Among individuals infected with SARS-CoV-2 during the index and alpha waves, an additional 33 individuals per 1000 took substantial sick leave within 1-9 months following acute infection compared to persons with no known infection. Females, older individuals, and individuals with pre-existing conditions such as fibromyalgia, COPD/chronic lung diseases, and obesity were affected. This study may be used to help inform the healthcare and research workforce of the impact of long covid on working ability and to motivate improved diagnostic and treatment options.
Funding
No specific funding was received for this work. The study was conducted as part of the governmental institution Statens Serum Institut’s advisory tasks for the Danish Ministry of Health.
Data Availability
The datasets used in this study comprise sensitive, individual-level information from completed questionnaires and national register data. According to the Danish data protection legislation, the authors are not permitted to share these sensitive data directly upon request. However, the data are available for research purposes upon request to the Danish Health Authority (register data, email: kontakt@sundhedsdata.dk) and Statens Serum Institut (questionnaire data, email: aii@ssi.dk), as well as within the framework of the Danish data protection legislation and any required permission from authorities. Data request processing can take an expected three to six months.
Contributors
This study was designed and initiated by PB, SE, AH, AK, AS, IB, and EO. The questionnaire was designed by AS, PB, NN, AK, SE.
Collection of data including programming related to this manuscript AS, JH, PB, IM, and LS. Data analysis was carried out by IB.
The first draft was written by EO, IB. Supervision of the study was carried out by AH. All authors have critically revised the manuscript.
All authors have approved the final version of the manuscript and agreed to be accountable for all aspects of the work.
EO and IB contributed equally to this manuscript. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted. EO and AH are the guarantors.
Competing interests
None.
Ethics approval
This study was performed as a surveillance study as part of the governmental institution Statens Serum Institut’s (SSI) advisory tasks for the Danish Ministry of Health. SSI’s purpose is to monitor and fight the spread of disease in accordance with section 222 of the Danish Health Act. According to Danish law, national surveillance activities carried out by SSI do not require approval from an ethics committee.
Participation in the study was voluntary. The invitation letter to participants contained information about their rights under the Danish General Data Protection Regulation (rights to access data, rectification, deletion, restriction of processing and objection). After reading this information, it was considered informed consent if participants agreed and clicked on the link to fill in the questionnaires.
Data sharing
The datasets used in this study comprise sensitive, individual-level information from completed questionnaires and national register data. According to the Danish data protection legislation, the authors are not permitted to share these sensitive data directly upon request. However, the data are available for research purposes upon request to the Danish Health Authority (register data, email: kontakt{at}sundhedsdata.dk) and Statens Serum Institut (questionnaire data, email: aii{at}ssi.dk), as well as within the framework of the Danish data protection legislation and any required permission from authorities. Data request processing can take an expected three to six months.
Supplementary Material
Acknowledgements
The authors would like to thank all EFTER-COVID participants for making this research possible by responding to the questionnaires.
We would also like to thank the members of the EFTER-COVID stakeholder group, including representatives from universities, long-term sequelae clinics, and other relevant hospital departments or health institutions, for offering their feedback on the questionnaires and contributing to stakeholder meetings.