A total of 1,176 patients were admitted to the hospital with a SARS-CoV-2 infection between February 2020 and February 2022. Among these patients, 651 attended the outpatient post-COVID-19 clinic 3 months post-discharge. A total of 126 patients completed the questionnaires at both outpatient assessments and were included in the study (Figure 1).
Frequency, severity, and course of physical, cognitive, and psychological symptoms
Table 1 shows the characteristics of the included patients. The median age of the patients was 68 years, and most patients were male (67%). The most common comorbidities were hypertension (42%), obesity (33%), and chronic respiratory disease (29%). Nearly half of the patients were included in the first wave of the study (49%).
Table 1: General characteristics of the hospitalized COVID-19 patients (n=126)
Demographic and premorbid factors
|
N (%)
|
Median [IQR]
|
Age in years
|
|
68 [61-76]
|
Male, n (%)
|
84 (67)
|
|
BMI in kg/m2
|
|
27.5 [25-31]
|
Co-morbidities, n (%) present
|
|
|
Hypertension
|
53 (42)
|
|
Obesity
|
41 (33)
|
|
Chronic respiratory disease
|
36 (29)
|
|
Type 2 diabetes
|
30 (24)
|
|
Chronic cardiac disease
|
30 (24)
|
|
Autoimmune disorder
|
19 (15)
|
|
Chronic neurologic disease
|
14 (11)
|
|
Rheumatologic disorder
|
12 (10)
|
|
Chronic kidney disease
|
10 (8)
|
|
Malignant neoplasm
|
8 (6)
|
|
CCI score
|
|
3 [2-4]
|
COVID-19 illness severity
|
|
|
Hospital stay in days
|
|
7 [5-14]
|
ICU admission, n (%)
|
22 (18)
|
|
Length of ICU stay in days
|
|
13 [6-34]
|
MC admission, n (%)
|
5 (4)
|
|
Length of MC stay in days
|
|
3 [1-6]
|
Days from discharge to T1
|
|
109 [99-129]
|
Days from discharge to T2
|
|
372 [351-405]
|
Oxygen treatments during hospital stay
|
|
|
Nasal oxygen therapy, n (%)
|
115 (91)
|
|
Non-invasive ventilation, n (%)
|
7 (6)
|
|
Invasive ventilation, n (%)
|
19 (15)
|
|
Severity score, n (%)
|
|
|
Moderate
|
35 (28)
|
|
Severe
|
68 (54)
|
|
Critical
|
23 (18)
|
|
Waves, n (%)
|
|
|
First
|
62 (49)
|
|
Second
|
24 (19)
|
|
Third
|
30 (24)
|
|
Fourth
|
10 (8)
|
|
Abbreviations: BMI, body mass index; CCI, Charlson Co-morbidity Index; ICU, intensive care unit; MC, medium care
As shown in Table 2, in total, 31 and 32% of the patients reported moderate-to-severe physical symptoms on the 4DSQ at T1 and T2, respectively. Patients reported a median of three cognitive symptoms at T1 and two at T2 on the CLCE-24, with 26 and 27% reporting six or more cognitive symptoms at T1 and T2, respectively. Additionally, 22 and 12% scored higher than the cut-off on the HADS anxiety T1 and T2, respectively. For HADS depression, 17 and 18% scored higher than the cut-off at T1 and T2, respectively. In total, 15 and 14% scored above the PTSS-14 cut-off at T1 and T2, respectively.
As shown in Table 3, 68% and 57% of the patients reported physical symptoms, that is, fatigue and dyspnoea, respectively. Among pulmonary function tests, DLCO, MEP, and TLC were impaired in 42%, 29%, and 18% of patients, respectively.
Wilcoxon signed rank tests showed that the HADS anxiety score decreased in 54 patients, it remained stable in 26, and increased in 36 patients (T=-2.542, p = 0.014). There were no significant changes in the other outcomes between T1 and T2 (p > .05).
Table 2: Results of outcome measures at T1 and T2
Outcome variables
|
n
|
Median [IQR]
|
>cut-off n(%)
|
n
|
Median [IQR]
|
>cut-off n (%)
|
|
T1
|
T2
|
4DSQ
|
53
|
5.0 [1.5-13.5]
|
17 (32)
|
124
|
5.5 [2.0-12.0]
|
38 (31)
|
CLCE-24
|
121
|
3 [0-6]
|
31 (26)*
|
116
|
2 [0-6]
|
31 (27)*
|
HADS-depression
|
118
|
2 [1-6]
|
20 (17)
|
125
|
3 [1-7]
|
23 (18)
|
HADS-anxiety
|
117
|
3 [1-7]
|
26 (22)
|
125
|
2 [0-5]
|
15 (12)
|
PTSS-14
|
115
|
22 [17-34]
|
17 (15)
|
118
|
21 [16-34]
|
16 (14)
|
Abbreviations: IQR, interquartile range; 4DSQ, Four-Dimensional Symptom Questionnaire (4DSQ); CLCE-24, Checklist for Cognitive and Emotional Consequences; (HADS), HADS-depression; Hospital Anxiety and Depression Scale-depression subscale; HADS-anxiety, Hospital Anxiety and Depression Scale-anxiety Subscale; PTSS-14, Post-traumatic Stress Symptoms Checklist-14.
* ≥6 physical symptoms
Table 3: Results of predictor variables at T1
Physical variables
|
Mean±SD, Median [IQR]
|
n (%)
|
|
% pred
|
Impaired*
|
FEV1
|
93.9±21.7
|
17/126 (14)
|
FVC
|
96.9±17.1
|
11/126 (9)
|
FEV1/FVC
|
98 [89-104]
|
19/126 (15)
|
TLC
|
97 [86-108]
|
22/122 (18)
|
RV
|
90 [81-104]
|
24/121 (20)
|
DLCO
|
77.6±0.5
|
52/125 (42)
|
DLCO/VA
|
86.6±20.1
|
34/125 (27)
|
MIP
|
98 [66-122]
|
12/125 (10)
|
MEP
|
89.6±34.0
|
36/125 (29)
|
COVID-related residual pulmonary abnormalities
|
|
105/123 (85)
|
FFMI, kg/m²
|
19 [17-20]
|
|
Fatigue
|
|
85/126 (68)
|
Dyspnoea
|
|
72/126 (57)
|
Psychological factors
|
|
low ; high **
|
CERQ Self-blame
|
4 [4-7]
|
0 (0) ; 2 (4)
|
CERQ Acceptance
|
10 [6-13]
|
6 (11) ; 1 (2)
|
CERQ Rumination
|
7 [5-9]
|
1 (2) ; 1 (2)
|
CERQ Positive refocusing
|
11.9±4.6
|
0 (0) ; 10 (19)
|
CERQ Planning
|
9 [7-13]
|
11 (21) ; 0 (0)
|
CERQ Positive reappraisal
|
10.5 (4.0)
|
4 (8) ; 0 (0)
|
CERQ Putting things in perspective
|
12 [9-16]
|
0 (0) ; 1 (2)
|
CERQ Catastrophizing
|
5 [4-8]
|
0 (0) ; 5 (9)
|
CERQ Other blame
|
4 [4-5]
|
0 (0) ; 1 (2)
|
Abbreviations: IQR, interquartile range; FEV1, forced expiratory volume in one second; FVC, forced vital capacity; TLC, total lung capacity; RV, residual volume; DLCO, diffusion capacity of the lungs for carbon monoxide; DLCO/VA, DLCO per unit alveolar volume; MIP, maximal inspiratory pressure; MEP, maximal expiratory pressure; FFMI, fat free mass index; CERQ, Cognitive Emotion Regulation Questionnaire; * Impaired = below lower limit of normal (LLN) ** low = < 2 SD; high = > 2 SD
Associations with outcomes at 12 months post-discharge
Physical Symptoms
Bivariable analyses identified significant associations between physical symptoms at T2 and chronic respiratory disease (r=.389, p<.001), COVID-related residual pulmonary abnormalities, fatigue, dyspnoea, and physical symptoms at T1 (r values ranging from .329 to .771, p<.001). Additionally, cognitive symptoms, anxiety, depression, PTSS at T1 (r values from .601 to .714, p<.001), and coping strategies (i.e. acceptance, catastrophizing, and rumination) (r values from .405 to .465, p<.01) were significantly associated with physical symptoms at T2. Hierarchical regression (Table 4) showed that chronic respiratory disease, dyspnoea, and higher anxiety at T1 were significantly associated with physical symptoms at T2, explaining 65.6% of variance.
Table 4. Regression analysis with 4DSQ at T2 as dependent variable (n=100)
|
β
|
Independent variables
|
Model 1
|
Model 2
|
Model 3
|
Model 4
|
Demographic and premorbid factors
|
|
|
|
|
Chronic respiratory disease
|
0.36***
|
NE
|
0.26**
|
0.20**
|
Physical factors
|
|
|
|
|
|
NE
|
NE
|
0.30***
|
0.03
|
|
NE
|
NE
|
0.24*
|
0.15 *
|
- COVID-related residual pulmonary abnormalities
|
NE
|
NE
|
0.03
|
0.02
|
Cognitive and psychological factors
|
|
|
|
|
|
NE
|
NE
|
NE
|
0.16
|
|
NE
|
NE
|
NE
|
-0.12
|
|
NE
|
NE
|
NE
|
0.54 ***
|
|
NE
|
NE
|
NE
|
0.12
|
R²
|
0.13
|
|
0.31
|
0.66
|
Adjusted R²
|
0.12
|
|
0.28
|
0.63
|
F change
|
14.79***
|
|
8.03***
|
22.68***
|
We used the square root transformation of 4DSQ as dependent variable
Abbreviations: 4DSQ, Four-Dimensional Symptom Questionnaire (4DSQ); CLCE-24, Checklist for Cognitive and Emotional Consequences; (HADS), HADS-depression; Hospital Anxiety and Depression Scale-depression subscale; HADS-anxiety, Hospital Anxiety and Depression Scale-anxiety Subscale; PTSS-14, Post-traumatic Stress Symptoms Checklist-14.
*p<.05; ** p<.01; ***p<.001
Cognitive Symptoms
Bivariable analyses indicated that cognitive symptoms at T2 were significantly associated with BMI, CCI, chronic respiratory disease (r=.198 to .368, p<.05), wave 3 vs. wave 1 (r=.275, p<.01), dyspnoea, fatigue, physical symptoms at T1 (r=.266 to .539, p<.01), and several psychological measures (i.e. cognitive symptoms, anxiety, depression, PTSS, and rumination) (r values from .293 to .708, p<.05). Hierarchical regression (Table 5) showed that chronic respiratory disease, wave 3 vs. wave 1, and cognitive symptoms at T1 were significantly associated with more cognitive symptoms at T2, explaining 63.2% of variance.
Table 5. Regression analysis with CLCE-24 at T2 as dependent variable (n=93)
|
β
|
Independent variables
|
Model 1
|
Model 2
|
Model 3
|
Model 4
|
Demographic and premorbid factors
|
|
|
|
|
- BMI
|
0.20*
|
0.21*
|
0.16
|
0.04
|
- CCI
|
0.13
|
0.14
|
0.19*
|
0.16*
|
- Chronic respiratory disease
|
0.32**
|
0.31 **
|
0.22*
|
0.20**
|
COVID-19 illness severity
|
|
|
|
|
- Wave 3 versus wave 1
|
|
0.29**
|
0.27**
|
0.17*
|
Physical factors
|
|
|
|
|
- Fatigue
|
NE
|
NE
|
0.34***
|
0.15
|
- Dyspnoea
|
NE
|
NE
|
0.21*
|
0.13
|
Psychological factors
|
|
|
NE
|
|
- CLCE-24
|
NE
|
NE
|
NE
|
0.50***
|
- HADS-depression
|
NE
|
NE
|
NE
|
-0.00
|
- HADS-anxiety
|
NE
|
NE
|
NE
|
-0.16
|
- PTSS-14
|
NE
|
NE
|
NE
|
0.22
|
R²
|
0.17
|
0.25
|
0.43
|
0.63
|
Adjusted R²
|
0.14
|
0.22
|
0.39
|
0.59
|
F change
|
6.00***
|
10.06**
|
13.28***
|
11.28***
|
We used the log transformation of CLCE-24 as dependent variable
Abbreviations: CLCE-24, Checklist for Cognitive and Emotional Consequences; (HADS), HADS-depression; Hospital Anxiety and Depression Scale-depression subscale; HADS-anxiety, Hospital Anxiety and Depression Scale-anxiety Subscale; PTSS-14, Post-traumatic Stress Symptoms Checklist-14.
*p<.05; ** p<.01; ***p<.001
Depression
Bivariable analyses found that depressive symptoms at T2 were significantly associated with chronic respiratory disease, wave 3 vs. wave 1 (r=.274 to .347, p<.01), TLC, dyspnoea, fatigue, physical symptoms at T1 (r=-.204 to .460, p<.05), cognitive symptoms, anxiety, depression, PTSS at T1 (r values from .536 to .716, p<.001), and coping strategies (i.e. rumination, positive refocusing, putting things in perspective) (r=.277 to -.389, p<.05). Regression analysis (Table 6) revealed that chronic respiratory disease, lower TLC, depressive symptoms, and PTSS at T1 were significantly associated with more depressive symptoms at T2, explaining 62.0% of variance.
Anxiety
Bivariable analyses showed significant associations between levels of anxiety at T2 and wave 3 vs. wave 1 (r=.198, p<.05), fatigue, physical symptoms, TLC, COVID-related pulmonary abnormalities (r values from -.193 to .472, p<.001), cognitive symptoms, anxiety, depression, PTSS at T1 (r=.497 to .716, p<.001), and rumination and catastrophizing (r=.377 to .490, p<.01). Regression analysis (Table 6) identified that higher anxiety at T1 was significantly associated with higher anxiety levels at T2, explaining 56.3% of variance.
Table 6. Regression analysis with HADS-depression (n=103) and HADS-anxiety (n=106) at T2 as dependent variable
HADS-depression at T2
|
|
β
|
|
Independent variables
|
Model 1
|
Model 2
|
Model 3
|
Model 4
|
- Chronic respiratory disease
|
0.34***
|
0.35***
|
0.32***
|
0.19*
|
|
NE
|
0.23*
|
0.22*
|
0.13
|
|
NE
|
NE
|
0.19*
|
0.01
|
|
NE
|
NE
|
0.14
|
0.05
|
|
NE
|
NE
|
-0.25**
|
-0.16*
|
- CLCE-24
|
NE
|
NE
|
NE
|
-0.05
|
- HADS-depression
|
NE
|
NE
|
NE
|
0.34**
|
- HADS-anxiety
|
NE
|
NE
|
NE
|
-0.01
|
- PTSS-14
|
NE
|
NE
|
NE
|
0.38**
|
R²
|
0.12
|
0.17
|
0.32
|
0.62
|
Adjusted R²
|
0.11
|
0.15
|
0.29
|
0.58
|
F change
|
12.53**
|
5.92*
|
6.96***
|
17.24***
|
HADS-anxiety at T2
|
|
|
β
|
|
Independent variables
|
Model 1
|
Model 2
|
Model 3
|
Model 4
|
|
- Chronic respiratory disease
|
0.23*
|
0.24*
|
0.25*
|
0.12
|
|
- Wave 3 vs wave 1
|
|
0.19
|
0.17
|
0.09
|
|
- Fatigue
|
|
|
0.19
|
-0.06
|
|
- TLC
|
|
|
-0.20
|
-0.13
|
|
- COVID-related residual pulmonary abnormalities
|
|
|
-0.01
|
-0.02
|
|
- CLCE-24
|
NE
|
NE
|
NE
|
-0.07
)
|
|
- HADS-depression
|
NE
|
NE
|
NE
|
0.06
|
|
- HADS-anxiety
|
NE
|
NE
|
NE
|
0.60***
|
|
- PTSS-14
|
NE
|
NE
|
NE
|
0.10
|
|
R²
|
0.05
|
0.09
|
0.18
|
0.56
|
|
Adjusted R²
|
0.05
|
0.07
|
0.14
|
0.52
|
|
F change
|
5.53*
|
3.75
|
3.52*
|
19.05***
|
|
We used the square root transformation of HADS-depression and log transformation of HADS-anxiety and as dependent variable
Abbreviations: CLCE-24, Checklist for Cognitive and Emotional Consequences; (HADS), HADS-depression; Hospital Anxiety and Depression Scale-depression subscale; HADS-anxiety, Hospital Anxiety and Depression Scale-anxiety Subscale; PTSS-14, Post-traumatic Stress Symptoms Checklist-14.
*p<.05; ** p<.01; ***p<.001
PTSS
Bivariable analyses showed significant associations between PTSS at T2 and chronic respiratory disease (r=.241, p<.01), wave 3 vs. wave 1 (r=.193, p<.05), COVID-related residual pulmonary abnormalities, fatigue, dyspnoea, and physical symptoms at T1 (r values from .214 to .596, p<.05). Depression, anxiety, cognitive symptoms, and PTSS (r values from .581 and .762, p<.001), and coping strategies (i.e. acceptance, rumination, and catastrophizing) (r values from .318 to .552, p<.001) at T1 were also significant. Regression analysis (Table 7) indicated that COVID-related residual pulmonary abnormalities, higher anxiety, and PTSS at T1 were significantly associated with increased PTSS at T2, explaining 68.3% of variance.
Table 7. Regression analysis with PTSS-14 at T2 as dependent variable (n=98)
|
β
|
Independent variables
|
Model 1
|
Model 2
|
Model 3
|
Model 4
|
Demographic and premorbid factors
|
|
|
|
|
- Chronic respiratory disease
|
-0.26*
|
-0.25*
|
-0.18
|
-0.08
|
COVID-19 illness severity
|
|
|
|
|
- Wave 3 versus wave 1
|
NE
|
-0.18
|
-0.15
|
-0.10
|
Physical factors at T1
|
|
|
|
|
- Fatigue
|
NE
|
NE
|
-0.36***
|
-0.09
|
- Dyspnoea
|
NE
|
NE
|
-0.13
|
-0.07
|
- COVID-related residual pulmonary abnormalities
|
NE
|
NE
|
-0.11
|
-0.13*
|
Cognitive and psychological factors at T1
|
|
|
|
|
- CLCE-24
|
NE
|
NE
|
NE
|
0.01
|
- HADS-depression
|
NE
|
NE
|
NE
|
-0.02
|
- HADS-anxiety
|
NE
|
NE
|
NE
|
-0.34*
|
- PTSS-14
|
NE
|
NE
|
NE
|
-0.39**
|
R²
|
0.07
|
0.10
|
0.28
|
0.68
|
Adjusted R²
|
0.06
|
0.08
|
0.24
|
0.65
|
F change
|
6.58*
|
3.26
|
7.78***
|
27.08***
|
We used the polynomial transformation of PTSS-14 as dependent variable
Abbreviations: CLCE-24, Checklist for Cognitive and Emotional Consequences; (HADS), HADS-depression; Hospital Anxiety and Depression Scale-depression subscale; HADS-anxiety, Hospital Anxiety and Depression Scale-anxiety Subscale; PTSS-14, Post-traumatic Stress Symptoms Checklist-14.
*p<.05; ** p<.01; ***p<.001
Sensitivity Analysis
Adding a fifth predictor block (physical symptoms and coping strategies at T1) did not significantly increase variance explained in the outcomes. However, more physical symptoms at T1 were associated with more physical symptoms at T2 (β=.639, 95% CI [.028-.165], p=.008).