Comparison of clinical characteristics between patients with pulmonary fibrosis and with no fibrosis/ resolution of pulmonary fibrosis after 90 days from onset.
Our study included all 457 confirmed COVID-19 cases admitted to the Shenzhen Third People’s Hospital and followed up till June 20, 2020. In order to observe the persistent pulmonary consequences of COVID-19 patients, 289 confirmed COVID-19 patients who was followed up more than 90 days after onset were further divided into two groups (group A and group B) according to the progression of pulmonary fibrosis. Group A (GA) had 116 (40.14%) patients who either had no lung fibrosis, or their lung fibrosis disappeared within 90 days after onset. On the other hand, 173 (59.86%) patients who still had lung fibrosis after 90 days from onset, were categorized as Group B (GB) (Table 1).
Among the physical characteristics studied, Age and BMI were found to be two significant (p < 0.05) risk factors between the two groups. Older patients (mean age 50.68 years vs 33.06 years) and patients with higher BMI (mean BMI 24.1 kg/m2 vs 22.1 kg/m2) still showed signs of lung fibrosis even after 90 days from onset. Most symptom profiles were comparable between GA and GB. However, a higher proportion of patients in GB had fever compared to GA (~73% in GB vs ~52% in GA) that is statistically significant at p < 0.0001. In addition, severe/critical COVID patients as well as patients with pre-existing health conditions were observed to be more vulnerable to sustaining lung fibrosis even after 90 days (Table 1).
It can also be seen that, patients in the GB took significantly longer time (21 days vs 14.7 days, p < 0.0001) from onset to get virus RNA negative, required longer follow-up period from onset (117.8 days vs 67.8 days, p < 0.0001) and stayed in the hospital for longer period of time (24.4 days vs 19.2 days, p < 0.0001). These observations in the temporal patterns can be attributed to the significantly higher number (29.5% vs 4.3%, p < 0.0001) of severe/critical COVID cases in GB, since severe/critical COVID patients generally require treatment and follow-up for a longer period time. It is interesting to see that an average delay of ~2 days before hospital admission following the onset of the symptoms had a significant (p < 0.0001) detrimental effect on the patients with pneumonia, thus rendering them more vulnerable to developing and sustaining lung fibrosis for a longer period of time (Table 1)
Levels of lactic acid, white blood cells (WBC), neutrophils (N), lymphocyte (HB), alanine transaminase (ALT), aspartate aminotransferase (AST), lactate dehydrogenase (LDH), creatine kinase (CK), urea nitrogen (BUN), and D-Dimer did not differ on a statistically significant level (p > 0.05) between the two groups (Table 1). However, there were significant differences in some of the laboratory findings between the two groups. These differences (GB vs GA) include lowest oxygenation index (PaO2/FiO2), Lowest CD4 cell, L (lymphocyte), CD4 cells, Highest CRP (C-reactive protein), Highest ESR (erythrocyte sedimentation rate), ESR, Highest lactic acid, Highest PCT, CRP, Cr (creatinine), Oxygenation index (PaO2/FiO2), ALB (albumin), PLT (platelet), Highest D-Dimer, Highest IL-6 (interleukin-6), IL-6, and PCT (procalcitonin). It is seen from the lowest oxygenation index (PaO2/FiO2) that the amount of oxygen in the blood drops drastically in patients who sustain fibrosis longer (Table 1).
CT characteristics and dynamic changes in pulmonary fibrosis
We analyzed the three CT scans for every patient. Ground-glass opacities (GGO), parenchymal bands, irregular interfaces, reticulation and traction bronchiectasis were the most common CT features in all COVID-19 patients. During hospitalization and follow-up, some patients had persistent pulmonary fibrosis (Figure 1), while some patients had resolution of pulmonary fibrosis (Figure 2). Typical CT imaging of a 67-year-old man showed diffuse ground glass opacities in both lungs, and visible parenchymal band in the lower lobe of left lung on initial CT. Diffuse ground glass opacities, consolidation and irregular interfaces with a small amount of pleural effusion were observed on the first follow-up CT. For the third (88 days after symptoms onset) and latest (132 days after symptoms onset) follow-up CT after discharge, diffuse ground glass opacities, consolidation and pleural effusion were still observed in the lungs, although most lesions were resolved (Figure 1). Typical CT imaging findings of a 53-year-old woman showed multiple lesions, a mass of ground glass opacities, consolidation and irregular interfaces on initial CT. The lesions were resolved obviously on the 6th day after onset, further resolved on the 9th and 19th days after onset, and completely resolved in both lungs on 108th day after symptoms onset (Figure 2).
Some patients did not follow-up along with the time extension, a total of 457 patients was included during the 0-30 days after onset in this study. During the 31-60, 61-90, 91-120 and >120 days after onset, there were 418, 279, 207 and 79 patients included, respectively, among them, 397 (86.87%), 311 (74.40%), 222 (79.57%), 141 (68.12%) and 49 (62.03%) patients developed with pulmonary fibrosis, respectively. It is interesting to see that pulmonary fibrosis developed in COVID-19 patients could be reversed. Resolution of pulmonary fibrosis were found in 18 (4.53%), 49 (13.61%), 14 (6.31%), 30 (21.289%) and 15 (30.61%) COVID-19 patients during the 0-30, 31-60, 61-90, 91-120 and >120 days after onset, respectively (Table 2). There was a total of 397 patients had lung fibrosis, among them, 126 patients had pulmonary fibrosis reversed, the median time of resolution of pulmonary fibrosis was 70.79±37.67 days (median: 59 days, range: 8-148 days). In contrast, there were 271 patients who still had pulmonary fibrosis at the last follow-up. The overall lung fibrosis and overall resolution of fibrosis percentage were 86.87% and 31.74%, respectively (Table 2). The CT imaging features including ground-glass opacities, meshwork, parenchymal bands, irregular interface and traction bronchiectasis of the COVID-19 patients at different stages were analyzed meanwhile. Ground-glass opacities and parenchymal bands are most common CT features in these patients.
Establishment of prediction model of the persistence of pulmonary fibrosis
A total of 56 features (Table 1) were collected from each patient. After performing the ANOVA, a set of 35 statistically significant (p < 0.05) remained for developing the prediction model. However, features related to treatment measures and temporal patterns were not considered for model development since those features are determined by the physicians based on the condition of the patients and thereby are dependent variables.
Among the significant independent variables, it was observed that Age, BMI, Fever, and Highest PCT were predictive factors for sustaining fibrosis even after 90 days (Figure 3A). Among the 289 patients, 288 patients had complete data for these four parameters that were used to develop a predictive model based-on the Logistic regression method. The performance of the model was validated by 5-fold cross validation and evaluated by the ROC, accuracy, PPV, NPV, Sensitivity and Specificity. Cross validation ensures the consistency of the model’s performance while reducing model bias and variance. The confusion matrix of the five-fold cross validation is presented in Figure 4B. An average AUC (Area under the ROC Curve) of 0.84 obtained from 5-fold cross validation that affirms a good reliability of the predictive model. The accuracy, PPV, NPV, Sensitivity and Specificity of the model was 76%, 71%, 79%, 67%, and 82%, respectively (Figure 4C).
Pulmonary function of patients with pulmonary fibrosis and with no fibrosis/ resolution of pulmonary fibrosis after 90 days from onset.
In order to observe the effect of pulmonary fibrosis on lung function, 33 patients from Group A (GA) and 114 patients from Group B (GB) underwent pulmonary function testing after 90 days from onset. Six patients (18.18%) in GA and 26 patients (22.81%) in GB had a pulmonary function abnormality. Two patients (6.06%) in GA and 5 patients (4.39%) in GB were diagnosed with obstructive pulmonary disorder. Restrictive disorder was found in 1 patient (3.03%) in GA and 8 patients (7.02%) in GB. In addition, 3 patients (9.09%) in GA and 13 patients (11.40%) in GB had small airway disorder. MEF 25 decrease was found in 12 patients (36.36%) in GA and 46 patients (40.35%) patients in GB. 2 patients (6.06%) in GA and 8 patients (7.02%) in GB had a pulmonary diffusion abnormality (Table 3). Overall, 18 patients (55.55%) and 76 patients (66.67%) revealed abnormal condition in pulmonary function test. It is interesting to see that the ratio of abnormal lung function, obstructive disorder, restrictive disorder, small airway disorder, MEF 25 decrease, abnormal lung diffusion and overall abnormalities in lung function test did not differ on a statistically significant level (p > 0.05) between the two groups.