Abstract

Background

The long-term pulmonary sequelae of mild coronavirus disease 2019 (COVID-19) remains unknown. In this study, we aimed to characterize lung function trajectories in individuals with mild COVID-19 from preinfection to 2 years postinfection.

Methods

We reinvited participants 2 years after infection from our matched cohort study of the Copenhagen General Population who had initially been examined 5.4 months after infection. We repeated lung tests and questionnaires. Linear mixed models were used to estimate dynamics in lung volumes in individuals with COVID-19 patients versus uninfected controls over two intervals: from pre-infection to 6 months postinfection and 6 months postinfection to 2 years postinfection.

Results

52 individuals (48.6%) attended the 2-year examination at median 1.9 years (interquartile range, 1.8–2.4) after COVID-19, all with mild infection. Individuals with COVID-19 had an adjusted excess decline in forced expiratory volume in 1 second (FEV1) of 13.0 mL per year (95% confidence interval [CI], −23.5 to −2.5; P = .02) from before infection to 6 months after infection compared to uninfected controls. From 6 to 24 months after infection, they had an excess decline of 7.5 mL per year (95% CI, −25.6–9.6; P = .40). A similar pattern was observed for forced vital capacity (FVC). Participants had a mean increase in diffusing capacity for carbon monoxide (DLco) of 3.33 (SD 7.97) between the 6- and 24-month examination.

Conclusions

Our results indicate that mild COVID-19 infection affects lung function at the time of infection with limited recovery 2 years after infection.

Pulmonary sequelae following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection have been widely reported, leading to an increased understanding of lung function trajectories at short to mid-term follow-up [1–7]. Prospective cohort studies and meta-analyses investigating lung function and lung abnormalities at 1-year follow-up report impaired diffusion capacity and restrictive pulmonary dysfunction as the most frequent pulmonary consequences of coronavirus disease 2019 (COVID-19). The prevalence of these findings varies between 8% and 31% at 1-year follow-up in survivors hospitalized with COVID-19 [8–11]. However, it remains unknown if the initial lung impairment regresses or progresses. Reports from SARS and Middle East respiratory syndrome (MERS) suggest that up to a third of admitted patients experienced pulmonary dysfunction, and in particular restrictive dysfunction, for up to 3 years following infection [12]. Therefore, there is a need for long-term follow-up studies of patients with COVID-19 to evaluate the prognosis and the potential recovery of the initial pulmonary impairment.

Furthermore, it has been established that the impact of COVID-19 on lung function is not limited to severe cases. Studies of mild COVID-19 have demonstrated an apparent decline in dynamic volumes and carbon monoxide diffusing capacity (DLco) in 16%–36% of patients evaluated from 3 to 5 months after infection [13, 14]. Our previous study on pulmonary function decline following mild COVID-19 showed an excess decline in dynamic lung volumes compared to uninfected individuals [15]. Taken together, this indicates that even mild infection causes pulmonary parenchymal damage. However, the extent and duration of this decline remains unclear. In this study, we aimed to investigate the long-term impact of mild COVID-19 on lung function and respiratory symptoms by extending our previous analyses with an additional follow-up visit 2 years after infection [15].

METHODS

Study Population and Approvals

Participants from the Copenhagen General Populations Study (CGPS) COVID-19 substudy on lung function and respiratory symptoms were invited to attend a second follow-up examination at approximately 2 years after infection. Results from the early follow-up (median 5.4 months after infection) on lung function have previously been reported [15]. The study was approved by the regional ethics committee (H-KF-01-144/0 and H-20072296) and written and oral consent were obtained from all participants. The study report follows the STROBE guidelines for cohort studies [16].

All invited individuals had performed baseline (prior to infection) lung function testing and a questionnaire on respiratory symptoms as a part of CGPS between January 2013 and November 2020 as previously described [17]. Examinations at 6 and 24 months after COVID-19 were both performed at Copenhagen University Hospital-Hvidovre Hospital, the 6-month visit was conducted from March to June of 2021 and the 24-month visit during September and October of 2022. Participants had tested positive for SARS-CoV-2 infection between 9 March 2020 and 16 January 2021.

Data Collection and Lung Function

Information on sex (male/female), age (years), body mass index (BMI; kg/m2), and smoking status (classified as current/former/never smokers) at baseline was obtained from the CGPS, along with information on cumulative tobacco consumption reported as pack years (number of cigarettes per day divided by 20 × years). The questionnaire on self-reported respiratory morbidity and symptoms was administered at all 3 visits and reviewed by the investigator together with the participant. Dyspnea was defined according to the modified Medical Research Council scale, with a score of ≥2 defined as dyspnea [18]. Sputum was defined as a daily productive cough for more than 3 months each year.

SARS-CoV-2 polymerase chain reaction (PCR) test results were obtained from microbiology departments serving the area of greater Copenhagen [15]. At 6 months after infection, information on symptoms and duration of symptoms related to COVID-19 was obtained from a specific questionnaire. Information on respiratory support were obtained from electronic medical records if participants were admitted to hospital in relation to acute infection.

Spirometry was performed at all 3 visits and diffusing capacity was measured at 6 and 24 months after infection. Baseline spirometry was conducted using a Portable EasyOne Plus ultrasonic spirometer (ndd Medical) [19]. At 6 and 24 months after infection, dynamic volumes, DLco, and total lung capacity were measured using a portable single-breath diffusing capacity device (EasyOne Pro; ndd Medical Technologies). DLCO was measured using the single-breath maneuver technique with carbon monoxide as a tracer gas. Spirometry and DLCO tests were performed according to published standards and as previously described [15, 20]. The predicted values for forced expiratory volume in 1 second (FEV1), forced vital capacity (FVC), FEV1/FVC, and lower limit of normal (LLN) were calculated based on the reference equations provided by the Global Lung Function Initiative [20].

Statistical Methods

Normality of all variables was checked using histograms and QQ plots. Continuous variables were reported as mean with standard deviation (SD) or median with interquartile range (IQR) and for comparison t test, paired t test, or nonparametric tests were used as appropriate. Categorical variables were reported as absolute values and percentages and compared using McNemar test for comparison of proportions.

Controls were individuals from the GCPS with spirometry data available prior to the pandemic and without a positive SARS-CoV-2 PCR test result. For each individual up to 5 controls from the GCPS were identified as previously described [15]. Individuals with COVID-19 were matched based on age at baseline and at 6 months after infection (±2.5 years), as well as sex and smoking status at baseline.

Our primary outcome was a comparison of the mean decline (slope) in FEV1 per year between individuals with COVID-19 and uninfected controls from baseline to 6 months after infection and from 6 to 24 months after infection. Our secondary outcome was a comparison of the mean decline (slope) in FVC per year. An additional secondary outcome was the change in respiratory symptoms in individuals with COVID-19 from baseline to 24 months after infection. We used a linear mixed model adjusted for age, sex, height, smoking status, and years at baseline. An interaction link between COVID-19 and time between examination was used to assess if individuals with COVID-19 had an increased decline. Results from the regression analysis were reported as change in mL with 95% confidence intervals (CI) and P values. A 2-sided P value of .05 was considered statistically significant. Analyses were performed using R statistical software version 4.0.1.

RESULTS

Clinical Characteristics

Of 107 individuals who attended the 6-month examination, 58 (54.2%) accepted our invitation to the attend the 24-month examination (Figure 1). Six individuals had been hospitalized in relation to infection and were excluded from all tables and analyses. Individuals who attended the 6- and 24-month examinations after infection were comparable with the original study cohort (Supplementary Table 1). Time between baseline and 6 months after infection was median 5.7 years (IQR, 3.8–8.7) and there was a median of 1.9 years (IQR, 1.8–2.4) from infection to the 24-month examination. At baseline, the 52 infected individuals and the 219 controls had comparable characteristics (Table 1). The mean age of individuals with COVID-19 at baseline was 56.2 years (SD 7.7), approximately half were female, 42.3% were never smokers, and the mean cumulative tobacco consumption at baseline was 8.9 pack-years (SD 12.1) (Table 1). Among the 52 individuals with COVID-19, 3 reported respiratory comorbidities such as asthma, and 4 (5.2%) had a spirometry result below the LLN.

Flow diagram of study participants. Abbreviations: CGPS, Copenhagen General Population Study; PCR, polymerase chain reaction; SARS-CoV-2, severe acute respiratory syndrome coronavirus 2. Figure created with Biorender.com.
Figure 1.

Flow diagram of study participants. Abbreviations: CGPS, Copenhagen General Population Study; PCR, polymerase chain reaction; SARS-CoV-2, severe acute respiratory syndrome coronavirus 2. Figure created with Biorender.com.

Table 1.

Baseline Characteristics of Individuals With COVID-19

CharacteristicIndividuals With COVID-19
(n = 52)
Uninfected Controls
(n = 219)
P Value
Age, y, mean (SD)56.2 (7.7)54.7 (7.4).24
Female sex, n (%)25 (48)114 (52.1)1.00
Ethnicity, Caucasian, n (%)52 (100)219 (100)1.00
BMI, kg/m2, mean (SD)26.6 (4.2)26.1 (4.7).53
Smoking status
 Never smokers, n (%)21 (40.4)100 (45.7).97
 Former smokers, n (%)22 (42.3)98 (44.7)
 Current smokers, n (%)5 (9.6)21 (10.4)
Cumulative smoking, pack-years, mean (SD)8.9 (12.1)21 (9.6).62
Clinical respiratory characteristics
 Self-reported asthma, n (%)a3 (5.7)11 (5.0).89
 Dyspnea, n (%)2 (3.8)7 (3.2)
 Wheezing, n (%)6 (11.5)26 (11.8).97
 Sputum, n (%)3 (7.6)22 (10.0).69
Dynamic lung volumes at baseline
 FEV1, L, mean (SD)3.1 (0.8)3.1 (0.8).84
 FEV1 predicted %, mean (SD)95.5 (15.2)95.5 (15.2).84
 FVC, L, mean (SD)4.3 (1.0)4.0 (1.0).24
 FVC predicted %, mean (SD)101.4 (13.7)97.9 (14.6).12
 FEV1/FVC, mean (SD)0.74 (0.1)0.77 (0.1).03
 FEV1/FVC < LLN, n (%)4 (5.2)11 (4.4).44
COVID-19 characteristics
 Duration of symptoms, d, median (IQR)13 (3–14)
 Time between baseline and 1st follow-up, y, median (IQR)5.7 (3.8–8.7)
 Time from positive SARS-CoV-2 RT-PCR to first follow-up, y, median (IQR)0.4 (0.3–1.0)
 Time between positive SARS-CoV-2 RT-PCR and 2nd follow-up, y, median (IQR)1.9 (1.8–2.4)
CharacteristicIndividuals With COVID-19
(n = 52)
Uninfected Controls
(n = 219)
P Value
Age, y, mean (SD)56.2 (7.7)54.7 (7.4).24
Female sex, n (%)25 (48)114 (52.1)1.00
Ethnicity, Caucasian, n (%)52 (100)219 (100)1.00
BMI, kg/m2, mean (SD)26.6 (4.2)26.1 (4.7).53
Smoking status
 Never smokers, n (%)21 (40.4)100 (45.7).97
 Former smokers, n (%)22 (42.3)98 (44.7)
 Current smokers, n (%)5 (9.6)21 (10.4)
Cumulative smoking, pack-years, mean (SD)8.9 (12.1)21 (9.6).62
Clinical respiratory characteristics
 Self-reported asthma, n (%)a3 (5.7)11 (5.0).89
 Dyspnea, n (%)2 (3.8)7 (3.2)
 Wheezing, n (%)6 (11.5)26 (11.8).97
 Sputum, n (%)3 (7.6)22 (10.0).69
Dynamic lung volumes at baseline
 FEV1, L, mean (SD)3.1 (0.8)3.1 (0.8).84
 FEV1 predicted %, mean (SD)95.5 (15.2)95.5 (15.2).84
 FVC, L, mean (SD)4.3 (1.0)4.0 (1.0).24
 FVC predicted %, mean (SD)101.4 (13.7)97.9 (14.6).12
 FEV1/FVC, mean (SD)0.74 (0.1)0.77 (0.1).03
 FEV1/FVC < LLN, n (%)4 (5.2)11 (4.4).44
COVID-19 characteristics
 Duration of symptoms, d, median (IQR)13 (3–14)
 Time between baseline and 1st follow-up, y, median (IQR)5.7 (3.8–8.7)
 Time from positive SARS-CoV-2 RT-PCR to first follow-up, y, median (IQR)0.4 (0.3–1.0)
 Time between positive SARS-CoV-2 RT-PCR and 2nd follow-up, y, median (IQR)1.9 (1.8–2.4)

Abbreviations: BMI, body mass index; COVID-19, coronavirus disease 2019; FEV1, forced expiratory volume in first second; FVC, forced vital capacity; LLN, lower limit of normal; RT-PCR, reverse transcription polymerase chain reaction; SARS-CoV-2, severe acute respiratory syndrome coronavirus 2.

aOnly for current or former smokers.

Table 1.

Baseline Characteristics of Individuals With COVID-19

CharacteristicIndividuals With COVID-19
(n = 52)
Uninfected Controls
(n = 219)
P Value
Age, y, mean (SD)56.2 (7.7)54.7 (7.4).24
Female sex, n (%)25 (48)114 (52.1)1.00
Ethnicity, Caucasian, n (%)52 (100)219 (100)1.00
BMI, kg/m2, mean (SD)26.6 (4.2)26.1 (4.7).53
Smoking status
 Never smokers, n (%)21 (40.4)100 (45.7).97
 Former smokers, n (%)22 (42.3)98 (44.7)
 Current smokers, n (%)5 (9.6)21 (10.4)
Cumulative smoking, pack-years, mean (SD)8.9 (12.1)21 (9.6).62
Clinical respiratory characteristics
 Self-reported asthma, n (%)a3 (5.7)11 (5.0).89
 Dyspnea, n (%)2 (3.8)7 (3.2)
 Wheezing, n (%)6 (11.5)26 (11.8).97
 Sputum, n (%)3 (7.6)22 (10.0).69
Dynamic lung volumes at baseline
 FEV1, L, mean (SD)3.1 (0.8)3.1 (0.8).84
 FEV1 predicted %, mean (SD)95.5 (15.2)95.5 (15.2).84
 FVC, L, mean (SD)4.3 (1.0)4.0 (1.0).24
 FVC predicted %, mean (SD)101.4 (13.7)97.9 (14.6).12
 FEV1/FVC, mean (SD)0.74 (0.1)0.77 (0.1).03
 FEV1/FVC < LLN, n (%)4 (5.2)11 (4.4).44
COVID-19 characteristics
 Duration of symptoms, d, median (IQR)13 (3–14)
 Time between baseline and 1st follow-up, y, median (IQR)5.7 (3.8–8.7)
 Time from positive SARS-CoV-2 RT-PCR to first follow-up, y, median (IQR)0.4 (0.3–1.0)
 Time between positive SARS-CoV-2 RT-PCR and 2nd follow-up, y, median (IQR)1.9 (1.8–2.4)
CharacteristicIndividuals With COVID-19
(n = 52)
Uninfected Controls
(n = 219)
P Value
Age, y, mean (SD)56.2 (7.7)54.7 (7.4).24
Female sex, n (%)25 (48)114 (52.1)1.00
Ethnicity, Caucasian, n (%)52 (100)219 (100)1.00
BMI, kg/m2, mean (SD)26.6 (4.2)26.1 (4.7).53
Smoking status
 Never smokers, n (%)21 (40.4)100 (45.7).97
 Former smokers, n (%)22 (42.3)98 (44.7)
 Current smokers, n (%)5 (9.6)21 (10.4)
Cumulative smoking, pack-years, mean (SD)8.9 (12.1)21 (9.6).62
Clinical respiratory characteristics
 Self-reported asthma, n (%)a3 (5.7)11 (5.0).89
 Dyspnea, n (%)2 (3.8)7 (3.2)
 Wheezing, n (%)6 (11.5)26 (11.8).97
 Sputum, n (%)3 (7.6)22 (10.0).69
Dynamic lung volumes at baseline
 FEV1, L, mean (SD)3.1 (0.8)3.1 (0.8).84
 FEV1 predicted %, mean (SD)95.5 (15.2)95.5 (15.2).84
 FVC, L, mean (SD)4.3 (1.0)4.0 (1.0).24
 FVC predicted %, mean (SD)101.4 (13.7)97.9 (14.6).12
 FEV1/FVC, mean (SD)0.74 (0.1)0.77 (0.1).03
 FEV1/FVC < LLN, n (%)4 (5.2)11 (4.4).44
COVID-19 characteristics
 Duration of symptoms, d, median (IQR)13 (3–14)
 Time between baseline and 1st follow-up, y, median (IQR)5.7 (3.8–8.7)
 Time from positive SARS-CoV-2 RT-PCR to first follow-up, y, median (IQR)0.4 (0.3–1.0)
 Time between positive SARS-CoV-2 RT-PCR and 2nd follow-up, y, median (IQR)1.9 (1.8–2.4)

Abbreviations: BMI, body mass index; COVID-19, coronavirus disease 2019; FEV1, forced expiratory volume in first second; FVC, forced vital capacity; LLN, lower limit of normal; RT-PCR, reverse transcription polymerase chain reaction; SARS-CoV-2, severe acute respiratory syndrome coronavirus 2.

aOnly for current or former smokers.

Between the 6- and 24-month examination, 18 individuals (34.6%) were reinfected with COVID-19, but none required hospitalization. All had completed a primary series of COVID-19 vaccination at the time of reinfection.

Change in Lung Function Volumes and Respiratory Symptoms

In individuals with COVID-19, we observed a mean annual decline of 46.1 mL per year (SD 41.0) for FEV1 and 76.7 mL per year (IQR, 41.6–97.3) for FVC from baseline to 6 months after infection. From 6 to 24 months after infection, median annual decline for FEV1 was 13.1 mL per year (IQR, −66.0 to 54.4) and for FVC it was 0.3 mL per year (IQR, −23.1 to 27.5). Table 2 presents change in lung volumes. Figure 2 shows the dynamics in respiratory symptoms from baseline to 24 months after infection. Two individuals (3.4%) reported dyspnea before infection, which increased to 7 participants (13.5%, P = .10) at 6 months after infection and to 8 (15.4%, P = .78) at 24 months after infection. Additionally, 6 (11.5%) individuals reported dyspnea for the first time at 6 months after infection and continued to report this 24 months after infection.

Respiratory symptoms during follow-up in individuals with coronavirus disease 2019 (COVID-19). Shown is the proportion of participants reporting dyspnea at each visit. Dyspnea was defined as modified Medical Research Council score ≥ 2, and sputum as daily productive cough in 3 months.
Figure 2.

Respiratory symptoms during follow-up in individuals with coronavirus disease 2019 (COVID-19). Shown is the proportion of participants reporting dyspnea at each visit. Dyspnea was defined as modified Medical Research Council score ≥ 2, and sputum as daily productive cough in 3 months.

Table 2.

Lung Function Measures and Respiratory Symptoms in Individuals With COVID-19 Who Attended All 3 Visits

Lung FunctionBaseline
(n = 52)
First Follow-up
(n = 52)
Second Follow-up
(n = 52)
P Value
FEV1, L, mean (SD)3.1 (0.8)2.9 (0.7)2.86 (0.8)
FEV1 < 80%, n (%)6 (12.1)11 (22.4)5 (13.8)
FVC, L, mean (SD)4.3 (1.0)3.8 (0.9)3.8 (1.0)
FVC < 80%, n (%)3 (6.9)8 (17.2)3 (10.3)
FEV1/FVC, mean (SD)0.7 (0.1)0.8 (0.1)0.8 (0.1)
Mean decline in FEV1, mL/y, mean/median (SD/IQR)46.1 (41.0)13.1 (−66.0 to 54.4).02
Mean decline in FVC, mL/y, median (IQR)76.7 (41.6– 97.3)0.3(−23.1 to 27.5).01
DLCO, mmol/min/kPa, mean (SD)82.3 (2.0)8.5 (2.3).02
Change in DLco, mmol/min/kPa/y, mean (SD)0.27 (0.7)
DLCO predicted, %, mean (SD)90.8 (15.9)94.9 (18.4)<.01
DLco < 80% predicted, n (%)10 (25.9)9 (22.4).56
Change in predicted DLco, %/y, mean (SD)3.9 (7.7)
DLco/Va, mmol/min/kPa/L, mean (SD)1.5 (0.2)1.4 (0.2).09
Change in DLco/Va, mmol/min/kPa/L/y, mean (SD)0.0 (0.1)
DLco/Va predicted, %, mean (SD)99.7 (15.7)98.6 (15.0).21
DLco/Va < 80% predicted, n (%)5 (13.8)5 (15.5)1.00
Change in predicted DLco/Va, %, mean (SD)−1.5 (8.0)
TLC, L, mean (SD)5.8 (1.1)6.1 (1.3)<.01
Change in TLC, L, mean (SD)0.3 (0.5)
TLC predicted, %, mean (SD)91.3 (10.1)96.5 (12.5)<.01
Change in predicted TLC, mean (SD)5.1 (7.2)
Lung FunctionBaseline
(n = 52)
First Follow-up
(n = 52)
Second Follow-up
(n = 52)
P Value
FEV1, L, mean (SD)3.1 (0.8)2.9 (0.7)2.86 (0.8)
FEV1 < 80%, n (%)6 (12.1)11 (22.4)5 (13.8)
FVC, L, mean (SD)4.3 (1.0)3.8 (0.9)3.8 (1.0)
FVC < 80%, n (%)3 (6.9)8 (17.2)3 (10.3)
FEV1/FVC, mean (SD)0.7 (0.1)0.8 (0.1)0.8 (0.1)
Mean decline in FEV1, mL/y, mean/median (SD/IQR)46.1 (41.0)13.1 (−66.0 to 54.4).02
Mean decline in FVC, mL/y, median (IQR)76.7 (41.6– 97.3)0.3(−23.1 to 27.5).01
DLCO, mmol/min/kPa, mean (SD)82.3 (2.0)8.5 (2.3).02
Change in DLco, mmol/min/kPa/y, mean (SD)0.27 (0.7)
DLCO predicted, %, mean (SD)90.8 (15.9)94.9 (18.4)<.01
DLco < 80% predicted, n (%)10 (25.9)9 (22.4).56
Change in predicted DLco, %/y, mean (SD)3.9 (7.7)
DLco/Va, mmol/min/kPa/L, mean (SD)1.5 (0.2)1.4 (0.2).09
Change in DLco/Va, mmol/min/kPa/L/y, mean (SD)0.0 (0.1)
DLco/Va predicted, %, mean (SD)99.7 (15.7)98.6 (15.0).21
DLco/Va < 80% predicted, n (%)5 (13.8)5 (15.5)1.00
Change in predicted DLco/Va, %, mean (SD)−1.5 (8.0)
TLC, L, mean (SD)5.8 (1.1)6.1 (1.3)<.01
Change in TLC, L, mean (SD)0.3 (0.5)
TLC predicted, %, mean (SD)91.3 (10.1)96.5 (12.5)<.01
Change in predicted TLC, mean (SD)5.1 (7.2)

Abbreviations: COVID-19, coronavirus disease 2019; DLco, diffusing capacity for carbon monoxide; FEV1, forced expiratory volume in first second; FVC, forced vital capacity; IQR, interquartile range; TLC, total lung capacity; Va, alveolar volume.

Table 2.

Lung Function Measures and Respiratory Symptoms in Individuals With COVID-19 Who Attended All 3 Visits

Lung FunctionBaseline
(n = 52)
First Follow-up
(n = 52)
Second Follow-up
(n = 52)
P Value
FEV1, L, mean (SD)3.1 (0.8)2.9 (0.7)2.86 (0.8)
FEV1 < 80%, n (%)6 (12.1)11 (22.4)5 (13.8)
FVC, L, mean (SD)4.3 (1.0)3.8 (0.9)3.8 (1.0)
FVC < 80%, n (%)3 (6.9)8 (17.2)3 (10.3)
FEV1/FVC, mean (SD)0.7 (0.1)0.8 (0.1)0.8 (0.1)
Mean decline in FEV1, mL/y, mean/median (SD/IQR)46.1 (41.0)13.1 (−66.0 to 54.4).02
Mean decline in FVC, mL/y, median (IQR)76.7 (41.6– 97.3)0.3(−23.1 to 27.5).01
DLCO, mmol/min/kPa, mean (SD)82.3 (2.0)8.5 (2.3).02
Change in DLco, mmol/min/kPa/y, mean (SD)0.27 (0.7)
DLCO predicted, %, mean (SD)90.8 (15.9)94.9 (18.4)<.01
DLco < 80% predicted, n (%)10 (25.9)9 (22.4).56
Change in predicted DLco, %/y, mean (SD)3.9 (7.7)
DLco/Va, mmol/min/kPa/L, mean (SD)1.5 (0.2)1.4 (0.2).09
Change in DLco/Va, mmol/min/kPa/L/y, mean (SD)0.0 (0.1)
DLco/Va predicted, %, mean (SD)99.7 (15.7)98.6 (15.0).21
DLco/Va < 80% predicted, n (%)5 (13.8)5 (15.5)1.00
Change in predicted DLco/Va, %, mean (SD)−1.5 (8.0)
TLC, L, mean (SD)5.8 (1.1)6.1 (1.3)<.01
Change in TLC, L, mean (SD)0.3 (0.5)
TLC predicted, %, mean (SD)91.3 (10.1)96.5 (12.5)<.01
Change in predicted TLC, mean (SD)5.1 (7.2)
Lung FunctionBaseline
(n = 52)
First Follow-up
(n = 52)
Second Follow-up
(n = 52)
P Value
FEV1, L, mean (SD)3.1 (0.8)2.9 (0.7)2.86 (0.8)
FEV1 < 80%, n (%)6 (12.1)11 (22.4)5 (13.8)
FVC, L, mean (SD)4.3 (1.0)3.8 (0.9)3.8 (1.0)
FVC < 80%, n (%)3 (6.9)8 (17.2)3 (10.3)
FEV1/FVC, mean (SD)0.7 (0.1)0.8 (0.1)0.8 (0.1)
Mean decline in FEV1, mL/y, mean/median (SD/IQR)46.1 (41.0)13.1 (−66.0 to 54.4).02
Mean decline in FVC, mL/y, median (IQR)76.7 (41.6– 97.3)0.3(−23.1 to 27.5).01
DLCO, mmol/min/kPa, mean (SD)82.3 (2.0)8.5 (2.3).02
Change in DLco, mmol/min/kPa/y, mean (SD)0.27 (0.7)
DLCO predicted, %, mean (SD)90.8 (15.9)94.9 (18.4)<.01
DLco < 80% predicted, n (%)10 (25.9)9 (22.4).56
Change in predicted DLco, %/y, mean (SD)3.9 (7.7)
DLco/Va, mmol/min/kPa/L, mean (SD)1.5 (0.2)1.4 (0.2).09
Change in DLco/Va, mmol/min/kPa/L/y, mean (SD)0.0 (0.1)
DLco/Va predicted, %, mean (SD)99.7 (15.7)98.6 (15.0).21
DLco/Va < 80% predicted, n (%)5 (13.8)5 (15.5)1.00
Change in predicted DLco/Va, %, mean (SD)−1.5 (8.0)
TLC, L, mean (SD)5.8 (1.1)6.1 (1.3)<.01
Change in TLC, L, mean (SD)0.3 (0.5)
TLC predicted, %, mean (SD)91.3 (10.1)96.5 (12.5)<.01
Change in predicted TLC, mean (SD)5.1 (7.2)

Abbreviations: COVID-19, coronavirus disease 2019; DLco, diffusing capacity for carbon monoxide; FEV1, forced expiratory volume in first second; FVC, forced vital capacity; IQR, interquartile range; TLC, total lung capacity; Va, alveolar volume.

Change in Diffusing Capacity for Carbon Monoxide

DLco was assessed at 6 and 24 months after infection. At 6 months after infection, the mean DLco was 90.0% of predicted (SD 15.9), which increased to 94.9% of predicted (SD 18.4) at the 24-month examination (P < .01). At 6 months after infection, 10 participants (25.9%) had diffusion impairment, which decreased to 9 (22.4%) at the 24-month examination. See Table 2 for more details on the DLco measurements.

Association of COVID-19 With Longitudinal Change in Lung Function Volumes

The longitudinal dynamics of FEV1 and FVC in individuals with COVID-19 and uninfected controls are shown in Figure 3. From our mixed model, we observed that from baseline to approximately 6 months after infection, individuals with COVID-19 had an excess decline of 13.0 mL per year (95% CI, −23.5 to −2.5; P = .02; Table 3). From 6 months after infection to 24 months, they had an excess decline of 7.5 mL per year (95% CI, −25.6–9.6; P = .40) compared to uninfected controls. In terms of FVC, individuals with COVID-19 had an excess decline of 33.1 mL per year from baseline to 6 months follow-up (95% CI, −48.5 to −17.8; P < .01) compared to 16.1 mL (95% CI, −39.1–5.0; P = .15) in the second time period. Throughout the entire study duration, from baseline to 24 months after infection, individuals with COVID-19 had an excess decline in FEV1 of 3.1 mL per year (95% CI, −32.9–5.5; P = .57) and in FVC of 19.8 mL per year (95% CI, −32.9 to −5.5; P < .01).

Dynamic lung volumes in the two study periods, from baseline to 6 months after infection and from 6 to 24 months after infection. Shown is the estimate of the slope in dynamic lung volumes from our linear mixed effects model for individuals with COVID-19 and uninfected controls. Abbreviations: COVID-19, coronavirus disease 2019; FEV1, forced expiratory volume in first second; FVC, forced vital capacity.
Figure 3.

Dynamic lung volumes in the two study periods, from baseline to 6 months after infection and from 6 to 24 months after infection. Shown is the estimate of the slope in dynamic lung volumes from our linear mixed effects model for individuals with COVID-19 and uninfected controls. Abbreviations: COVID-19, coronavirus disease 2019; FEV1, forced expiratory volume in first second; FVC, forced vital capacity.

Table 3.

Adjusted Linear Mixed Effect Analyses of the Annual Change in Lung Volumes of Individuals With COVID-19 Compared to Uninfected Controls

Change in Lung VolumesDifference, mL, (95% CI)P Value
From baseline to 6 mo after infection
 FEV1−13.0 (−23.5 to −2.5).02
 FVC−33.1 (−48.5 to −17.8)<.01
From 6 mo to 2 y after infection
 FEV1−7.5 (−25.6 to 9.6).40
 FVC−16.1 (−39.1 to 5.0).15
Change in Lung VolumesDifference, mL, (95% CI)P Value
From baseline to 6 mo after infection
 FEV1−13.0 (−23.5 to −2.5).02
 FVC−33.1 (−48.5 to −17.8)<.01
From 6 mo to 2 y after infection
 FEV1−7.5 (−25.6 to 9.6).40
 FVC−16.1 (−39.1 to 5.0).15

Each of the four analyses models the slope of FEV1 or FVC between the two examinations. The first two models, represents the slope of FEV1 and FVC from baseline to 6 months after infection. The last 2 models model the slope of FEV1 and FVC from 6 months to 2 years after infection. All models were adjusted for age, sex, height, smoking status, and cumulative smoking at baseline and COVID-19 status.

Abbreviations: COVID-19, coronavirus disease 2019; FEV1, forced expiratory volume in 1 second; FVC, forced vital capacity.

Table 3.

Adjusted Linear Mixed Effect Analyses of the Annual Change in Lung Volumes of Individuals With COVID-19 Compared to Uninfected Controls

Change in Lung VolumesDifference, mL, (95% CI)P Value
From baseline to 6 mo after infection
 FEV1−13.0 (−23.5 to −2.5).02
 FVC−33.1 (−48.5 to −17.8)<.01
From 6 mo to 2 y after infection
 FEV1−7.5 (−25.6 to 9.6).40
 FVC−16.1 (−39.1 to 5.0).15
Change in Lung VolumesDifference, mL, (95% CI)P Value
From baseline to 6 mo after infection
 FEV1−13.0 (−23.5 to −2.5).02
 FVC−33.1 (−48.5 to −17.8)<.01
From 6 mo to 2 y after infection
 FEV1−7.5 (−25.6 to 9.6).40
 FVC−16.1 (−39.1 to 5.0).15

Each of the four analyses models the slope of FEV1 or FVC between the two examinations. The first two models, represents the slope of FEV1 and FVC from baseline to 6 months after infection. The last 2 models model the slope of FEV1 and FVC from 6 months to 2 years after infection. All models were adjusted for age, sex, height, smoking status, and cumulative smoking at baseline and COVID-19 status.

Abbreviations: COVID-19, coronavirus disease 2019; FEV1, forced expiratory volume in 1 second; FVC, forced vital capacity.

From our linear mixed model, we observed that COVID-19 alone, had an impact of −214.5 mL (95% CI, −398.1 to −33.4; P = .04) on FEV1 and −300.6 mL (95% CI, −658.9 to −180.3; P < .01) on FVC from baseline to the 6-month visit. Subsequently, from 6 months postinfection to the 24-month visit, the impact on FEV1 was −145.4 mL (95% CI, −349.3 to 26.2; P = .1) and on FVC −247.1 mL (95% CI, −565.3 to −104.7; P = .02). See Supplementary Table 2 for a full list of effect estimates from our linear mixed model.

DISCUSSION

This study presents results on long-term follow-up assessment of lung function in individuals with COVID-19. By comparing longitudinal changes in individuals with mild COVID-19 and uninfected controls, our results indicate that mild COVID-19 leads to a significant excess decline in lung function during the first 6 months after infection and a slower decline thereafter, indicating lasting decline with limited recovery. However, we observed that the acceleration of the decline slowed and eventually followed the expected annual age-associated decline approximately 6 months postinfection.

A limited number of studies have investigated the long-term lung function trajectories after COVID-19 in hospitalized and nonhospitalized individuals 2 years after acute infection. The studies report conflicting results. One followed 288 hospitalized individuals at 6 months, 1 year, and 2 years [21]. Of these, 45 participants did not receive oxygen upon admission and therefore may be comparable to our study population. Interestingly, they found an unadjusted increase in FVC from admission to 6 months after infection of 10 mL, suggesting an initial recovery phase in the first 6 months. However, from the 12- to 24-month follow-up, they observed an unadjusted decline of 50 mL interpreted as an age-related decline. Another study on 144 hospitalized individuals reports a small but significant unadjusted increase in both FEV1 and FVC from 6 to 24 months across all levels of severity [22]. It is worth noting that both studies had no information on lung function volumes prior to infection, and the reported lung function volumes were not adjusted for factors that impact lung function, features that may explain the discrepancy between the studies. Moreover, the studies did not include unexposed controls, which may limit the ability to assess the effect of COVID-19. In contrast, our study included matched controls and observed an adjusted decline of 13 mL per year and 33 mL in FEV1 and FVC from prior to infection to 6 months, and seem to resume age-related decline from 6 to 24 months after infection. Our results suggest that COVID-19 results in lasting lung impairment with limited recovery.

In our study cohort, there was an increase in predicted DLCO of 3.9% over the approximately 18 months between measurements, suggesting a potential improvement, which is consistent with studies conducted in hospitalized individuals [22–25]. A 2-year follow-up study utilizing computed tomography (CT) chest scans and DLco measurements at 6, 12, and 24 months reported an increase in DLco from 80% to 84% between the 6- and 24-month examination [22]. Notably, this study observed the most significant improvement in individuals with nonfibrotic interstitial abnormalities, such as ground glass opacities, and reported gradual improvement over time up to 2 years, concluding that the observed improvement in DLco was mainly due to the resolution of these radiologic abnormalities. In this context, we have previously shown that individuals from this study cohort showed radiologic lesions at the first follow-up at 5.4 months [26]. The improvement in DLco reported here might be due resolution of these abnormalities, although we do not have chest CT scans at the 24-month examination to verify this. Conversely, the study by Zhang et al investigating DLCO trajectories reported conflicting results [21]. They observed a median absolute decrease in percent predicted DLCO of 2.2 across all levels of disease severity in hospitalized individuals from admission to the 2-year follow-up period [21]. However, it is important to note that the participants in Zhang et al's study were more severely ill, which could account for the different trajectories observed.

Despite demonstrating a mean increase in DLco from 6 to 24 months after infection, a corresponding decrease in dyspnea was not observed. At 24 months after infection, only 1 of 7 individuals who reported dyspnea for the first time 6 months after infection did not experience dyspnea. Furthermore, it is notable that participants reported dyspnea prior to infection, emphasizing both the variability in respiratory symptoms and importance of preinfection information. Previous research at mid- to long-term follow-up has demonstrated a decrease in respiratory symptoms over time simultaneously with an improvement in DLco [27, 28]. These conflicting results could be due to the limited number of participants reporting dyspnea. Furthermore, the association of DLco and respiratory symptoms remains elusive and future research with larger cohorts and more comprehensive assessments of respiratory symptoms is warranted to gain a clearer understanding of the relationship between DLco and dyspnea.

Our study has several strengths. First, it is a prospective cohort study with lung function measured at baseline and followed up for 2 years. This gave us a unique opportunity for a more detailed description of the before and long-term lung function trajectories in relation to infection. Second, matched uninfected controls were included, which made it possible to estimate whether any decline in lung function was due to natural age-related decline or related to infection. Last, the study focused on nonhospitalized individuals, making it more generalizable as most individuals experience mild infection. However, there were also several limitations to the study. In terms of generalizability, a limitation is that the entire study cohort was Caucasian, thereby constraining the applicability of the results to a more diverse population. Another limitation is that almost 50% of participants were lost to follow-up from the 6- to 24-month examination, which may have introduced nonresponse bias. However, sensitivity analysis indicated that the full cohort had the same decline from prior infection to 6 months after infection, suggesting that the sample is representative of the full cohort making this a minor concern. Another limitation is the study size, which limits the ability to perform subgroup analyses among hospitalized patients, individuals with lung disease, or those with varying smoking history and symptom duration. Such analyses could have provided deeper insights into potential diverse trajectories and the identification of factors contributing to a more rapid decline in individuals with mild infection. Additionally, there was no information on diffusion capacity of the lungs (DLco) and total lung capacity prior to infection, which made it difficult to estimate temporal changes, but this is a defect in nearly all respiratory disease studies.

In conclusion, our study found evidence of long-term lung function impairment among individuals who had recovered from mild COVID-19. The decline in dynamic lung volumes was greater in affected individuals compared to matched controls and persisted up to 2 years after the initial infection. Further research is needed to fully understand the implications of these findings, particularly as the global burden of COVID-19 continues.

Supplementary Data

Supplementary materials are available at The Journal of Infectious Diseases online (http://jid.oxfordjournals.org/). Supplementary materials consist of data provided by the author that are published to benefit the reader. The posted materials are not copyedited. The contents of all supplementary data are the sole responsibility of the authors. Questions or messages regarding errors should be addressed to the author.

Notes

Author contributions. A. R., S. A., and T. B. conceived the idea and planned the study. K. I. collected the data and performed the statistical analysis. K. I., A. R., and T. B. wrote the draft manuscript. All authors contributed to subsequent writing, critical revisions, and approved the final version of the manuscript. T. B. and B. G. N. obtained funding and provided administrative, technical, and material support.

Acknowledgments. We thank all individuals with COVID-19 for their participation.

Financial support. This work was supported by the Research Council of Rigshospitalet; AP Møller og Hustru Chastine McKinney Møllers Fond; and the Brodie Foundation.

References

1

Trinkmann
 
F
,
Müller
 
M
,
Reif
 
A
, et al.   
Residual symptoms and lower lung function in patients recovering from SARS-CoV-2 infection
.
Eur Respir J
 
2021
;
57
:
2003002
.

2

Arnold
 
DT
,
Hamilton
 
FW
,
Milne
 
A
, et al.   
Patient outcomes after hospitalisation with COVID-19 and implications for follow-up: results from a prospective UK cohort
.
Thorax
 
2021
;
76
:
399
401
.

3

Zhang
 
S
,
Bai
 
W
,
Yue
 
J
, et al.   
Eight months follow-up study on pulmonary function, lung radiographic, and related physiological characteristics in COVID-19 survivors
.
Sci Rep
 
2021
;
11
:
13854
.

4

Lewis
 
KL
,
Helgeson
 
SA
,
Tatari
 
MM
,
Mallea
 
JM
,
Baig
 
HZ
,
Patel
 
NM
.
COVID-19 and the effects on pulmonary function following infection: a retrospective analysis
.
EClinicalMedicine
 
2021
;
39
:
101079
.

5

Torres-Castro
 
R
,
Vasconcello-Castillo
 
L
,
Alsina-Restoy
 
X
, et al.   
Respiratory function in patients post-infection by COVID-19: a systematic review and meta-analysis
.
Pulmonology
 
2021
;
27
:
328
37
.

6

Fumagalli
 
A
,
Misuraca
 
C
,
Bianchi
 
A
, et al.   
Pulmonary function in patients surviving to COVID-19 pneumonia
.
Infection
 
2021
;
49
:
153
7
.

7

Guler
 
SA
,
Ebner
 
L
,
Aubry-Beigelman
 
C
, et al.   
Pulmonary function and radiological features 4 months after COVID-19: first results from the national prospective observational Swiss COVID-19 lung study
.
Eur Respir J
 
2021
;
57
:
2003690
.

8

Wu
 
X
,
Liu
 
X
,
Zhou
 
Y
, et al.   
3-month, 6-month, 9-month, and 12-month respiratory outcomes in patients following COVID-19-related hospitalisation: a prospective study
.
Lancet Respir Med
 
2021
;
9
:
747
54
.

9

Pan
 
F
,
Yang
 
L
,
Liang
 
B
, et al.   
Chest CT patterns from diagnosis to 1 year of follow-up in patients with COVID-19
.
Radiology
 
2022
;
302
:
709
19
.

10

Magdy
 
DM
,
Metwally
 
A
,
Tawab
 
DA
,
Hassan
 
SA
,
Makboul
 
M
,
Farghaly
 
S
.
Long-term COVID-19 effects on pulmonary function, exercise capacity, and health status
.
Ann Thorac Med
 
2022
;
17
:
28
36
.

11

Cho
 
JL
,
Villacreses
 
R
,
Nagpal
 
P
, et al.   
Quantitative chest CT assessment of small airways disease in post-acute SARS-CoV-2 infection
.
Radiology
 
2022
;
304
:
185
92
.

12

Hui
 
DS
,
Joynt
 
GM
,
Wong
 
KT
, et al.   
Impact of severe acute respiratory syndrome (SARS) on pulmonary function, functional capacity and quality of life in a cohort of survivors
.
Thorax
 
2005
;
60
:
401
9
.

13

Zhao
 
Y-M
,
Shang
 
Y-M
,
Song
 
W-B
, et al.   
Follow-up study of the pulmonary function and related physiological characteristics of COVID-19 survivors three months after recovery
.
EClinicalMedicine
 
2020
;
25
:
100463
.

14

Krueger
 
T
,
van den Heuvel
 
J
,
van Kampen-van den Boogaart
 
V
, et al.   
Pulmonary function three to five months after hospital discharge for COVID-19: a single centre cohort study
.
Sci Rep
 
2023
;
13
:
681
.

15

Iversen
 
KK
,
Afzal
 
S
,
Ahlström
 
MG
, et al.   
Lung function decline in relation to COVID-19 in the general population: a matched cohort study with prepandemic assessment of lung function
.
J Infect Dis
 
2022
;
225
:
1308
16
.

16

von Elm
 
E
,
Altman
 
DG
,
Egger
 
M
, et al.   
The strengthening the reporting of observational studies in epidemiology (STROBE) statement: guidelines for reporting observational studies
.
J Clin Epidemiol
 
2008
;
61
:
344
9
.

17

Çolak
 
Y
,
Afzal
 
S
,
Nordestgaard
 
BG
,
Lange
 
P
,
Vestbo
 
J
.
Importance of early COPD in young adults for development of clinical COPD: findings from the Copenhagen general population study
.
Am J Respir Crit Care Med
 
2021
;
203
:
1245
56
.

18

Stenton
 
C
.
The MRC breathlessness scale
.
Occup Med (Lond)
 
2008
;
58
:
226
7
.

19

Çolak
 
Y
,
Nordestgaard
 
BG
,
Vestbo
 
J
,
Lange
 
P
,
Afzal
 
S
.
Comparison of five major airflow limitation criteria to identify high-risk individuals with COPD: a contemporary population-based cohort
.
Thorax
 
2020
;
75
:
944
54
.

20

Quanjer
 
PH
,
Stanojevic
 
S
,
Cole
 
TJ
, et al.   
Multi-ethnic reference values for spirometry for the 3–95-yr age range: the global lung function 2012 equations
.
Eur Respir J
 
2012
;
40
:
1324
43
.

21

Zhang
 
H
,
Li
 
X
,
Huang
 
L
, et al.   
Lung-function trajectories in COVID-19 survivors after discharge: a two-year longitudinal cohort study
.
EClinicalMedicine
 
2022
;
54
:
101668
.

22

Han
 
X
,
Chen
 
L
,
Fan
 
Y
, et al.   
Longitudinal assessment of chest CT findings and pulmonary function after COVID-19 infection
.
Radiology
 
2023
;
307
:
e222888
.

23

Kimmig
 
LM
,
Rako
 
ZA
,
Ziegler
 
S
, et al.   
Long-term comprehensive cardiopulmonary phenotyping of COVID-19
.
Respir Res
 
2022
;
23
:
263
.

24

Lee
 
JH
,
Yim
 
J-J
,
Park
 
J
.
Pulmonary function and chest computed tomography abnormalities 6–12 months after recovery from COVID-19: a systematic review and meta-analysis
.
Respir Res
 
2022
;
23
:
233
.

25

Huntley
 
CC
,
Patel
 
K
,
Bil Bushra
 
S-E-S
, et al.   
Pulmonary function test and computed tomography features during follow-up after SARS, MERS and COVID-19: a systematic review and meta-analysis
.
ERJ Open Res
 
2022
;
8
:
00056-2022
.

26

Iversen
 
KK
,
Ronit
 
A
,
Kristensen
 
TS
, et al.   
Lung abnormalities in severe acute respiratory syndrome coronavirus 2 infection: an analysis of paired computed tomography scans before and after infection
.
Open Forum Infect Dis
 
2022
;
9
:
ofac596
.

27

Fortini
 
A
,
Rosso
 
A
,
Cecchini
 
P
, et al.   
One-year evolution of DLCO changes and respiratory symptoms in patients with post COVID-19 respiratory syndrome
.
Infection
 
2022
;
50
:
513
7
.

28

Faverio
 
P
,
Luppi
 
F
,
Rebora
 
P
, et al.   
One-year pulmonary impairment after severe COVID-19: a prospective, multicenter follow-up study
.
Respir Res
 
2022
;
23
:
65
.

Author notes

Potential conflicts of interest. T. B. reports grants from Pfizer, Novo Nordisk Foundation, Simonsen Foundation, and Lundbeck Foundation; grants and personal fees from GSK, Pfizer, and Gilead; and personal fees from Boehringer Ingelheim and MSD, all outside the submitted work. All other authors report no potential conflicts.

All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed.

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/pages/standard-publication-reuse-rights)

Supplementary data