If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
# These authors contributed equally to this article.
Yi Guo
Footnotes
# These authors contributed equally to this article.
Affiliations
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
# These authors contributed equally to this article.
Hao Wang
Footnotes
# These authors contributed equally to this article.
Affiliations
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
# These authors contributed equally to this article.
Xin Guan
Footnotes
# These authors contributed equally to this article.
Affiliations
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Corresponding author. Xiaodong Li, Hubei Provincial Hospital of Traditional Chinese Medicine, 4 Garden Hill, Wuhan, 430061, Hubei, China, Tel: +86-27-88920956
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Department of Occupational and Environmental Health, Key Laboratory of Environment and Health, Ministry of Education and State Key Laboratory of Environmental Health (Incubating), School of Public Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Most participants had an improved physical and psychological health status.
•
We observed sequelae symptoms in 45.2% of participants at 18.5 months.
•
We observed abnormal chest computed tomography patterns in about 30% of participants at 18.5 months.
•
The convalescents still had mild depression or anxiety disorders at the last visit.
Abstract
Objectives
This study aimed to describe the full scope of long-term outcomes and the ongoing pathophysiological alterations among COVID-19 survivors.
Methods
We established a longitudinal cohort of 208 COVID-19 convalescents and followed them at 3.3 (interquartile range [IQR]: 1.3, 4.4, visit 1), 9.2 (IQR: 9.0, 9.6, visit 2), and 18.5 (IQR: 18.2, 19.1, visit 3) months after infection, respectively. Serial changes in multiple physical and psychological outcomes were comprehensively characterized. We, in addition, explored the potential risk factors of SARS-CoV-2 antibody response and sequelae symptoms.
Results
We observed continuous improvement of sequelae symptoms, lung function, chest computed tomography (CT), 6-minute walk test, and the Borg dyspnea scale, whereas sequelae symptoms (at least one) and abnormal chest CT patterns still existed in 45.2% and about 30% of participants at 18.5 months, respectively. Anxiety and depression disorders were alleviated for the convalescents, although depression status was sustained for a longer duration.
Conclusions
Most COVID-19 convalescents had an overall improved physical and psychological health status, whereas sequelae symptoms, residual lesions on lung function, exercise impairment, and mental health disorders were still observed in a small proportion of participants at 18.5 months after infection. Implementing appropriate preventive and management strategies for the ever-growing COVID-19 population is warranted.
SARS-CoV-2, with high covertness and high transmissibility, has caused more than 0.5 billion confirmed cases and 6.3 million deaths globally as of July 13, 2022 [
]. Huge burdens for the healthcare system and the whole of society may occur due to the overwhelming COVID-19 pandemic and the rapidly growing population of post-COVID-19 patients worldwide [
]. Emerging evidence has suggested that many COVID-19 survivors suffered from a higher rate of long-term complications and limited day-to-day activities [
], less attention is paid to the recovery condition of the myocardial injury. To date, little is known about the natural history of long-term COVID-19. There is still an immediate need for studies to explore the longer health outcomes and the ongoing pathophysiology alterations among COVID-19 survivors.
In this study, we established a longitudinal cohort of COVID-19 convalescents with different disease severity with over 18.5 months of follow-up. We comprehensively characterized the serial changes of multiple indicators, including sequelae symptoms, respiratory outcomes, computed tomography (CT) scans, physical function, a biomarker of myocardial injury, SARS-CoV-2 antibody response, and mental health disorders. We further explored the potential risk factors of the SARS-CoV-2 antibody and sequelae symptoms.
Methods
Study design
In this prospective cohort study, we invited discharged patients with COVID-19 from different hospitals in multiple districts of Wuhan, China. There were 289 convalescents who agreed to participate in this study and attend designated follow-up at Hubei Provincial Hospital of Traditional Chinese Medicine on February 17, 2020. A total of 81 non-consecutive patients were lost to follow-up due to declining participation (64 participants), leaving Wuhan (15 participants), and dying before visit 1 (2 participants). There were 208 participants who were included in the final analysis and participated in three visits. The median duration of three visits was 3.3 (interquartile range [IQR]: 1.3, 4.4, visit 1), 9.2 (IQR: 9.0, 9.6, visit 2), and 18.5 (IQR: 18.2, 19.1, visit 3) months after infection, respectively. The study population and severity classification, admission and discharge criteria, and the statistical analysis were further summarized in the supplementary eText. The flowchart of the study was also shown in supplementary eFigure 1. This study was approved by the Ethics Committee of the School of Public Health, Tongji Medical College, Huazhong University of Science and Technology (approval number: 202001). All participants provided written informed consent.
Procedure
At each visit, the participants underwent a detailed questionnaire interview, physical examination, routine blood test, pulmonary function tests, high-resolution chest CT (HRCT) scan, and a 6-minute walking test (6MWT). Self-reported sequelae symptoms include cough, fatigue or muscle weakness, sleep difficulties, decreased appetite, diarrhea or vomiting, smell or taste disorder, dizziness or headache, sore throat, and chest pain.
Lung function tests were performed before and after 6MWT according to guidelines of the American Thoracic Society [
]. Parameters consisted of forced expiratory volume in one second (FEV1), forced vital capacity (FVC), FEV1/FVC ratio, and forced expiratory flow between 25% and 75% of vital capacity (FEF25-75). All measurements were performed with a spirometer (SP80B, CONTEC Medicine, Qinhuangdao, China) and expressed as predicted percentages of normal values. Measurements assessed after 6MWT were shown in the results.
Participants underwent unenhanced chest CT examinations using a 40-section CT scanner with breath holding at the end of inspiration (uCT 530, United Imaging Healthcare, Shanghai, China). Images were reconstructed at 0.55 mm slice thickness, with 513 mm × 768 mm matrix size. The collimation and rotation time was 22 mm and 0.7 seconds, respectively. Other parameters were set according to the manufacturer's standard routine. CT abnormalities were demonstrated according to the terms of international standards defined by the Fleischner Society glossary and peer-reviewed COVID-19 literature [
]. In our study, two professional radiologists and one experienced pulmonologist constituted an evaluation group, and we invited two groups with the same specialist configuration to read CT reports separately and further assess ground-glass opacity (GGO) scores and reticular pattern (RP) scores. The average score of the two groups was considered the final score for statistical analysis. Details of the GGO scoring system are: 0, no involvement; 1, involvement <25%; 2, 25-50% involvement; 3, 50-75% involvement; 4, involvement >75%. The unilateral lung is divided into three sections: upper (above the level of the carina), middle, and lower (below the level of the right lower pulmonary vein). There are six sections on both sides, including right upper, right middle, right lower, left upper, left middle, and left lower. Each lobe is scored separately, and then the total GGO score of six lobes (0-24) is calculated. Details of the RP scoring systems are: 0, no lines; 1, 1-2 short line/grid shadows; 2, multiple line/grid shadows; 3, <1 lung segment's aggregated line/grid-like shadows in the subpleural area; 4, 1-3 lung segments' large grid shadows; 5, >3 lung segments' diffuse grid shadows and/or deformed lung structure. The right upper lobe contains three segments, the left upper lobe usually contains two segments, the right middle, and the left lingual lung contain two segments and the highest score is 4 points; the lower lobe of both lungs has more than two segments and the highest score is 5 points; thus, the maximum score is 26 points. Each lobe was scored separately, and then we calculated the total RP score of the six lobes (0-26 points). All participants were scored, excluding pre-morbid lung nodules or scars.
A 6MWT was performed on a walk monitoring and analysis system (YK2020A, Wocaring Medical Equipment, Wuhan, China), following the standardized protocol of the American Thoracic Society [
American Thoracic Society Committee on Proficiency Standards for Clinical Pulmonary Function Laboratories. ATS statement: guidelines for the six-minute walk test.
]. Participants walk as far as possible indoors within 6 minutes on a flat and hard surface. Distance and predicted percentage of walking distance were calculated based on Jay's [
] method. In addition, the Borg dyspnea scale was assessed after 6MWT, as participants were asked to grade their level of breath shortness from 0 (no dyspnea at all) to 10 (excessive dyspnea) [
For biomarkers of cardiac injury, we measured cardiac troponin T (TnT) using a highly sensitive reagent of TnT (hsTnT) on the Cobas E601 immune analyzer (Roche Diagnostics), following the manufacturer's instruction. The measuring range was 3-10000 ng/l, and the intermediate precision coefficient of variation was <10%. A 99th percentile value of 14 ng/l in the general reference population had been reported previously [
]. Myocardial injury is defined as circulating cardiac Tn levels higher than the 99th percentile upper reference limit, regardless of new abnormalities of electrocardiography and echocardiography [
In addition, we assessed the SARS-CoV-2 neutralizing antibody (NAb) level by a pseudotyped virus-based assay and calculated the half-maximal inhibitory concentration (50% neutralizing titer [NT50]) according to the inhibition rate of each dilution using nonlinear regression. The details have been described in a previous study [
To further evaluate participants’ mental health, the validated Chinese versions of the Patient Health Questionnaire-9 (PHQ-9) contained a 9-item depression module ranging from 0 to 27 points and was used to evaluate post-discharge depression status [
]. Meanwhile, the Generalized Anxiety Disorder-7 (GAD-7) included a 7-item anxiety module ranging from 0 to 21 points and was applied to assess post-discharge anxiety conditions [
]. Each above item scored from 0 (not at all) to 3 (nearly every day).
Results
The demographics and comorbidities are listed in Table 1. Among 208 COVID-19 convalescents, there were 146 mild cases and 62 severe cases. The median age of participants was 58.0 years (IQR: 50.0, 64.3), and 100 (48.1%) were men. A total of 68 (32.7%) participants had college or higher education attainment. There were 101 (48.6%) participants who had more than 50,000 RMB per year of household income. The majority of participants were never-smokers (n = 183, 88.0%) and never-drinkers (n = 183, 88.0%). There were 75 (36.1%), 23 (11.1%), and 14 (6.7%) participants who had pre-existing hypertension, diabetes mellitus, and cardiovascular disease (CVD). The median body mass index (BMI) and waist circumferences were 24.3 kg/m2 and 90 cm. There were no significant differences in most variables between mild and severe cases except for BMI, with severe cases having a significantly lower median BMI of 23.4 (IQR: 22.0, 25.5) (P <0.01, Supplementary eTable 1).
Table 1Basic characteristics of COVID-19 convalescents.
Variables
Total population (N = 208)
Age, years
58 (50.0, 64.3)
Gender, n (%)
Male
100 (48.1)
Female
108 (51.9)
Education, n (%)
Middle school or lower
140 (67.3)
College or higher
68 (32.7)
Household income, n (%)
<50000 RMB/year
107 (51.4)
≥50000 RMB/year
101 (48.6)
Cigarette smoking, n (%)
Never-smoker
183 (88.0)
Ever-smoker
25 (12.0)
Alcohol consumption, n (%)
Never-drinker
183 (88.0)
Ever-drinker
25 (12.0)
Comorbidity, n (%)
Hypertension
75 (36.1)
Diabetes
23 (11.1)
Cardiovascular disease
14 (6.7)
Body mass index
24.3 (22.6, 26.5)
Waist circumference
90 (83, 97)
Duration from symptom onset to the last follow-up, months
18.5 (18.2, 19.1)
Data were expressed as median (interquartile range) or frequency (percentage).
Table 2 illustrates the longitudinal trend in sequelae symptoms, lung function, chest CT, 6MWT, hs-TnT, NT50 for serum SARS-CoV-2 NAb, and psychological condition of the COVID-19 survivors during the follow-up. The percentage of convalescents with at least one sequelae symptom decreased from 62.0% at 3.3 months to 50.0% and 45.2% at 9.2 months and 18.5 months, respectively. The most commonly reported symptoms in each visit were fatigue or muscle weakness and sleep difficulties (35.1% and 33.2%, respectively, at visit 1), with the frequency gradually declining during the follow-up. A significant decrease was observed for most sequelae symptoms between visit 1 and visit 2, whereas no significant decrease was observed for sequelae between visit 2 and visit 3. The indicators for lung function were slightly increased during the convalescence, although only a significant increase was observed for FEV1/FVC ratio. The FEV1/FVC was 83.0% at visit 1 and increased to 85.7% at visit 3. FEF25-75% was 81.2 (IQR: 71.9, 97.0) % at visit 1, while 91.8 (IQR: 75.6, 110.0) % at visit 3. A total of 95 participants attended chest CT during each visit, and the frequency of GGO decreased from 69 (72.6%) at visit 1 to 28 (29.5%) at visit 3. The CT scores of GGO and RP also experienced a drastic decrease, which was 6 and 7, respectively, at visit 1, while 1 and 3, respectively, at visit 3. The median distance of the 6MWT significantly improved, from 514.9 m (IQR: 480.2, 556.0) at 3.3 months to 535 m (IQR: 509.0, 570.0) at 18.5 months. Similarly, significant improvements were found for the percentage of predicted value and the Borg dyspnea scale after 6MWT between visit 1 and visit 3. Only four participants (4.7%) experienced dyspnea (Borg dyspnea scale ≥1) after 6MWT at 18.5 months, whereas 30 people had a Borg scale ≥ 1 (50.8%) at 3.3 months. The levels of hsTnT experienced a slight decrease over time, from 4.5 ng/l (IQR: 3.0, 7.5) at visit 1 to 4.1 ng/l (IQR: 3.0, 6.6) at visit 3, and myocardial injury frequency (hsTnT ≥14 ng/l) decreased to 3.8% at visit 3 compared with 7.0% at visit 1, although no significant difference was found between each visit. NT50 for serum SARS-CoV-2 NAb experienced a continuous decrease in the unvaccinated group, from 1153.0 (IQR: 473.5, 2095.8) at visit 1 to 281.0 (IQR: 128.2, 499.8) at visit 3. Meanwhile, the trend in the vaccinated group showed a rapid decrease between visit 1 and visit 2, with 765.5 (IQR: 343.5, 1510.5) and 287.5 (IQR: 140.2, 585.0) respectively, and then increased to 910.0 (IQR: 497.0, 1460.8) due to the COVID-19 vaccination after 1 year of pandemic. For mental health disorders, 96 and 97 participants provided information on depression and anxiety disorders during each visit, respectively. Both depression and anxiety scores showed a decreasing trend during the follow-up, with a slight difference among the three visits. The depression score gradually decreased from 5 at 3.3 months to 4 at 9.2 months and then rapidly decreased to 1 at 18.5 months. Meanwhile, the anxiety disorder scores decreased rapidly during the first two visits (median score = 4 and 1 at visit 1 and visit 2, respectively) and remained relatively stable during the last visit (median score = 0). Over 20% of our participants were suffering from at least mild depression or anxiety status at 18.5 months. The overall trends of the above parameters were similar in mild and severe convalescents (Supplementary eTable 2). Compared with the mild convalescents, severe convalescents had a relatively higher frequency of any sequelae symptoms within 18.5 months and higher CT scores of GGO and RP during the early follow-up period, as shown in Supplementary eTable 3. Meanwhile, we compared the variation of lung function, CT, 6MWT, and Borg dyspnea scale during each visit between the mild and the severe convalescents and found that the mild convalescents had a substantially better improvement in Borg dyspnea scale between visit 3 and visit 1 (Δ = -1.000 [-1.000, 0.000], P-value = 0.025) compared with the severe (Supplementary eTable 4).
Table 2Physical and psychological health status of convalescents during follow-up.
Significant difference between visit 1 and visit 3 groups (P <0.05).
1.0 (0.0, 5.0)
0.0 (0.0, 4.0)
Anxiety score ≥5, n (%)
36 (37.1)
29 (29.9)
22 (22.7)
Data were expressed as median (interquartile range) or frequency (percentage). Kruskal–Wallis test was applied for group comparisons of continuous variables, and χ² test or Fisher Exact tests were performed to analyze the categorical variables.
CT: computed tomography; FEV1, forced expiratory volume in one second; FVC, forced vital capacity; FEF25-75, forced expiratory flow between 25% and 75% of vital capacity; GGO, ground-glass opacity; RP, reticular pattern; hsTnT, highly sensitive troponin T; NT50, the half-maximal inhibitory concentration.
a Significant difference between visit 1 and visit 2 groups (P <0.05).
b Significant difference between visit 1 and visit 3 groups (P <0.05).
c Significant difference between visit 2 and visit 3 groups (P <0.05).
Table 3 provided a summary of the potential risk factors of NT50 for serum SARS-CoV-2 NAb during follow-up. We established a linear mixed model to calculate the effect size of the number of any symptoms, CT scores of GGO and RP, and depression and anxiety scores on NT50 for serum SARS-CoV-2 NAb. The number of any symptoms, CT scores of GGO, and RP were positively associated with a higher level of NT50 for serum SARS-CoV-2 NAb, and the β (SE) were 0.126 (0.045), 0.045 (0.017), and 0.045 (0.015), respectively (all P <0.01). No significant associations were observed for depression and anxiety scores with NT50 (both P >0.05). Figure 1 illustrates the potential risk factors of any sequelae symptoms. Generalized linear mixed models were established to explore the associations of any sequelae symptoms with age group, gender, education, and comorbidities, including hypertension, diabetes mellitus, and CVD, and disease severity, with multivariable adjustment conducted. Compared with participants without hypertension, participants with hypertension had an odds ratio (OR) of 0.44 (95% confidence interval [CI] 0.24-0.80) for any sequelae symptoms. Meanwhile, severe participants had an OR of 3.60 (95% CI 1.97-6.58) compared with mild participants.
Table 3Potential risk factors associated with NT50 for serum SARS-CoV-2 neutralizing antibody during follow-up.
Variables
β (SE)
P-value
Number of any symptom
0.126 (0.045)
0.005
CT scores of ground-glass opacity
0.045 (0.017)
0.009
CT scores of reticular pattern
0.045 (0.015)
0.004
Depression score
-0.002 (0.015)
0.901
Anxiety score
0.019 (0.018)
0.276
Linear mixed model of repeated measures was conducted to estimate fixed effects of any symptom, CT scores, and mental health disorders on log-transformed NT50 of SARS-CoV-2 neutralizing antibody, adjusted for age, gender, and vaccination. Random effects consisted of days of the neutralizing antibody test after symptom onset and days of the first shot of the COVID-19 vaccine since symptom onset of each participant. The number of any symptoms refers to the number of sequelae of a participant.
CT: computed tomography; NT50, the half-maximal inhibitory concentration; SE, standard error of the beta coefficient.
Generalized linear mixed models were established to explore the associations of any sequelae symptoms (categorical) with fixed effects. Each patient was included in random effect. For the age group (<60 years/≥60 years), we adjusted gender (male/female) and education (middle school or lower/college or higher). Comorbidities including hypertension (no/yes), diabetes (no/yes), and CVD (no/yes), and disease severity (mild/severe) were excluded due to the potential mediating effects. For gender, we adjusted age group, education, comorbidities, and disease severity. For education, we adjusted age group, gender, and disease severity. Comorbidities were excluded due to the potential mediating effects. For comorbidities (hypertension, diabetes mellitus, and CVD), we adjusted each for age group, gender, education, and comorbidities except itself. Disease severity (mild/severe) was excluded due to the potential mediating effects. For disease severity, we adjusted age group, gender, education, and comorbidities.
In addition, all participants in this study were hospitalized and discharged after recovery (Supplementary eTable 5). The median days of hospitalization were 27 (IQR 17.0, 35.0) for the total population and 25.0 (IQR 15.2, 32.0) and 29.0 (IQR 22.0, 37.0), respectively, for the mild and severe participants. A total of 9 (4.3%) participants were admitted to the intensive care unit (ICU). Only one participant received treatment with invasive mechanical ventilation, whereas no one underwent a tracheostomy or developed renal failure that needed renal function replacement treatment. Compared with mild participants, severe participants had a significantly higher frequency of ICU admission, more intensive treatments, and longer hospital stays. As shown in Supplementary eTable 6, the obvious improvement of CT abnormalities, 6MWT, Borg dyspnea scale, and mental health disorders was observed in both groups (hospitalization <27 days/≥27 days). Substantial improvement of sequelae symptoms was further observed in participants with a stay of ≥27 days. In the vaccinated group, the NT50 levels underwent a decrease and an increase due to the vaccination after visit 2 in both participants with hospitalization less than and more than or equal to 27 days. In Supplementary eTable 7, we observed an obvious recovery trend in any sequelae symptoms, CT abnormalities, 6MWT, Borg dyspnea scale, and mental health disorders in the non-ICU group, while no substantial improvement was observed in the ICU group during three visits.
Discussion
Due to the COVID-19 pandemic, there has been growing concern that survivors might be at higher risk of multiple complications. Large-scale, longitudinal measurements and longer-term data are urgently needed to comprehensively characterize the trend of health consequences of COVID-19 convalescents. In this study, we comprehensively assessed the sequential pathophysiology changes of COVID-19 convalescents based on a considerably long follow-up cohort in Wuhan (over 18.5 months). We observed that most convalescents had substantial improvement in their general health status during the follow-up visits. Meanwhile, some convalescents still have the sequelae symptoms, including fatigue, muscle weakness, sleep difficulties, abnormal CT patterns, and depression and anxiety disorder at 18.5 months after diagnosis. Our finding provides novel evidence for the personalized prevention and intervention of long-term outcomes among COVID-19 convalescents.
Previous studies have reported that fatigue, muscle weakness, and sleep difficulties were the most frequent post-discharge sequelae symptoms and could last at least 12 months [
]. We confirmed this finding and advanced the evidence to 18.5 months after COVID-19 diagnosis. The potential pathogenesis of fatigue and muscle weakness may include viral-induced myositis, long periods of bed rest during convalescence, and the use of systemic corticosteroid therapy [
]. We found continuous improvement in chest CT in these participants. In the meantime, the CT appearance of the abnormal radiographic pattern was still detected in around 30% of convalescents at 18.5 months. The CT abnormalities could even be sustained beyond this duration. The improvement with concurrent residual lesions of lung CT was also observed in previous studies with 1 to 2 years of follow-up [
]. Meanwhile, we observed some fluctuations in the indicators of lung function and the 6MWT. The convalescents experienced fluctuations of FEV1%, FVC%, 6MWT, and predicted distance during the three visits, with an overall recovery trend. The possible mechanisms of this phenomenon include the recovery of early injury and the increase of pulmonary ventilation from 3 to 9 months after symptom onset, and then pulmonary atelectasis and parenchymal fibrosis occurred due to the persistent effect of inflammation in the later stages of the disease course [
]. This hypothesis is supported by the CT results in our study, as 35.8% of convalescents still had reticular patterns at 18.5 months. Previous studies of acute respiratory distress syndrome also reported similar trends in lung function and 6MWTs among SARS survivors [
]. Further studies are still warranted to confirm our findings and unravel the underlying mechanisms.
Notably, the results demonstrated that both depression and anxiety disorders were alleviated for the convalescents, although at least mild depression and anxiety status were sustained in over 20% of our participants at 18.5 months. A cohort study suggested that some COVID-19 survivors had psychiatric sequelae (mood, anxiety, or psychotic disorder) and more frequent substance use at 6 months post-COVID-19, although the longitudinal assessment of the outcomes was not available in the study [
] reported that more patients had anxiety or depression at 12 months than at 6 months (26% vs 23%). In addition, another study reported that after 2 years post-symptom onset, there were still 12% of COVID-19 survivors who had anxiety or depression, compared with 5% of participants in matched non-COVID-19 controls [
The effect of tocilizumab, anakinra and prednisolone on antibody response to SARS-CoV-2 in patients with COVID-19: a prospective cohort study with multivariate analysis of factors affecting the antibody response.
] were associated with antibody immune response, and inflammation mechanisms, including the number of clusters of differentiation 4+ T cells declining and the interleukin-6 level elevation, may partially explain the associations. We included the previously mentioned three indices in multivariable analysis using a linear mixed model and observed that the number of any symptoms, CT scores of GGO, and CT scores of RP were positively associated with higher NT50 for serum SARS-CoV-2 NAb. Higher NAb titer of SARS-CoV-2 indicates a potential persistent inflammation and immune response in our participants for 18.5 months. Previous studies demonstrated that negative psychological experiences, such as less social cohesion, were associated with lower vaccine efficacy [
]. However, no significant associations were observed for depression and anxiety scores with the NT50 in our participants due to the relatively small sample size. Hence, further psychological and random behavioral intervention studies with a larger sample size are needed to explore the potential mechanism.
Hypertension was associated with lower risks of any sequelae symptoms among our COVID-19 survivors, with adjusted OR = 0.44 (95% CI: 0.24-0.80). Some studies reported that hypertension was associated with higher risks of adverse outcomes in patients with COVID-19 [
] also suggested that hypertension was not an independent predictor of the lethal outcome in patients with COVID-19. A large cohort study included 153,760 COVID-19 individuals, and over 10 million controls suggested that the incident risk of cardiovascular outcome was lower in the hypertension group (hazard ratio [HR] = 1.57, 95% CI: 1.51-1.64) than in the normotensive group (HR = 1.66, 95% CI: 1.61-1.72) [
]. Previous studies have indicated that ACE inhibitor (ACEI)/angiotensin Ⅱ receptor blocker (ARB) medication played a protective role in the COVID-19 prognosis [
Association of inpatient use of angiotensin-converting enzyme inhibitors and angiotensin II receptor blockers with mortality among patients with hypertension hospitalized with COVID-19.
]. A retrospective study involving 1128 adult patients observed that compared with ACEI/ARB non-users, the use of ACEI/ARB among hospitalized patients with COVID-19 and hypertension was associated with a lower risk of all-cause mortality (HR = 0.42, 95% CI: 0.19-0.92) [
Association of inpatient use of angiotensin-converting enzyme inhibitors and angiotensin II receptor blockers with mortality among patients with hypertension hospitalized with COVID-19.
] also reported that patients with COVID-19 with hypertension receiving ACEI/ARB therapy had a lower rate of severe diseases and better clinical outcomes compared with non-ACEI/ARB users. Despite the previously mentioned evidence, the underlying mechanisms for the observed association remain largely unclear. The findings should be interpreted cautiously, and future studies are warranted to confirm our findings and unravel the potential mechanism.
The severe participants had a higher risk for any sequelae symptoms compared with the mild participants, with adjusted OR = 3.60 (95% CI: 1.97-6.58). This is in line with Cao and his colleagues’ studies [
]. They provided consistent findings that compared with mild patients, severe or critically ill patients had a significantly poor recovery and more post-COVID symptoms at 6, 12, and 24 months after discharge [
Compared with the wild-type strain, the current dominant variant Omicron differs in transmissibility, pathogenicity, and immune escape due to the deletions and many mutations in the spike protein [
]. Patients infected with Omicron had milder clinical manifestations and decreased disease severity and mortality compared with those who were infected with the wild-type strain [
], while 45.2% of our participants (infected with wild-type strain) still experienced any sequelae symptoms at 18.5 months after symptom onset. Therefore, the generalizability of our findings should be tested in other populations infected with new variants.
A longitudinal and comprehensive profile of the disease was observed at 18.5 months of follow-up in our study. Nevertheless, this investigation has several limitations. First, participants in this study were recruited during the early period of the COVID-19 pandemic; thus, the generalizability to other new SARS-CoV-2 variants is limited [
]. Second, more than 70% of our participants were mild patients at the initial stage; further investigations that include participants of various disease severity and larger sample size are needed. Third, not all participants attended three visits, thus, selection bias could be introduced during this process. Although there were 81 non-consecutive patients lost to follow-up, we performed the additional analysis to compare the basic characteristics between included (N = 208) and declined (N = 81) participants (Supplementary eTable 8). No significant difference in basic characteristics was observed for the two groups, which confirmed the robustness of our finding. The selection bias could be minimal. Fourth, we did not assess the transfer factor for carbon monoxide (TLCO) in the pulmonary function test. We intended to conduct a comprehensive study including multiple physical and mental health indices and only assessed the key indicators for each category. The TLCO should be tested in future studies to distinguish diffusion deficit.
In summary, this longitudinal study suggested that most COVID-19 convalescents had an improved physical and psychological health status in general, whereas post-discharge sequelae symptoms, residual lesions on lung function and exercise impairment, and mental health disorders were still observed in a small proportion of our participants for 18.5 months. Clinicians and policymakers should be aware of the risk of physical and mental complications in the ever-growing COVID-19 convalescents. The development of targeted strategies for the early prevention of post-discharge sequelae symptoms, CT abnormalities, and mental health problems is warranted.
Declaration of competing interest
No potential conflict of interest was reported by the author(s).
Funding
This work was supported by the Emergency Key Program of Guangzhou Laboratory (EKPG21-30), the National Science Foundation of China (82204113 and 72061137006), and the Fellowship of China Postdoctoral Science Foundation (2020T130034ZX).
Ethical approval
This study was approved by the Ethics Committee of the School of Public Health, Tongji Medical College, Huazhong University of Science and Technology (approval number: 202001). All participants provided written informed consent.
Acknowledgments
We thank all the study participants and project staff from Tongji Medical College, Huazhong University of Science and Technology, and Hubei Provincial Hospital of Traditional Chinese Medicine for the work they have done.
Author contributions
Yi Guo, Hao Wang, Mingzhong Xiao, and Xin Guan: conceptualization, data collection, analysis, interpretation, and writing; Yanshou Lei, Tingyue Diao, Pinpin Long, Rui Zeng, Xuefeng Lai: data collection and interpretation, review; Hao Cai, Yutong You, Yuying Wen, Wenhui Li, Xi Wang, Yufei Wang, Qinlin Chen, Yuchan Yang, Yutong Qiu, Jishuai Chen, Huidan Zeng: data collection and investigation; Wei Ni, Youyun Zhao, Kani Ouyang, Jingzhi Wang: conceptualization, methodology, and investigation; Qi Wang, Li Liu, Lulu Song, Youjie Wang, Huan Guo: resources, supervision, and review; Xiaodong Li, Tangchun Wu and Yu Yuan: study design and conceptualization, methodology, data interpretation, writing, and review; All authors have read and approved the final manuscript. Yi Guo, Hao Wang, Mingzhong Xiao, and Xin Guan are joint first authors; Xiaodong Li, Tangchun Wu and Yu Yuan are joint senior authors.
The effect of tocilizumab, anakinra and prednisolone on antibody response to SARS-CoV-2 in patients with COVID-19: a prospective cohort study with multivariate analysis of factors affecting the antibody response.
Association of inpatient use of angiotensin-converting enzyme inhibitors and angiotensin II receptor blockers with mortality among patients with hypertension hospitalized with COVID-19.