FormalPara Key Summary Points

The CT findings of patients with mild COVID-19 have not been explored in detail.

The existing total lesion CT score only considered the area of involvement of the lesions, but not the nature of the lesions.

We analyzed 606 CT scans of 182 patients and proposed a new scoring system based on the area and nature of the lesions.

The total lesions score of mild COVID-19 pneumonia peaked 17–18 days after illness onset, and week 2 after onset may be the critical period for disease progression.

The consolidation scores may objectively reflect the severity of lung involvement and thus serve as a more accurate diagnostic measure than total lesion scores.

The typical CT features of each stage of COVID-19 pneumonia were bilateral, multiple, and ill-defined, with a GGO pattern, which were located in the posterior and peripheral regions of both lungs.

Good prognosis may be associated with increases in the proportions of the GGO and linear opacity patterns during the later stage of disease.

Digital Features

This article is published with digital features, including a summary slide, to facilitate understanding of the article. To view digital features for this article go to https://doi.org/10.6084/m9.figshare.13007915.

Introduction

The World Health Organization (WHO) reported that as of July 6, 2020, there were ten million confirmed cases of coronavirus disease 2019 (COVID-19) worldwide, including 532,340 deaths [1]. These figures are updated daily, and subsequently, the situation has largely deteriorated depending on geographical location. Among confirmed cases reported in China, 81% of 72,314 patients experienced mild pneumonia [2]. Computed tomography (CT) is instrumental for screening and effective evaluation of patients with COVID-19 [3], and understanding CT performance is of vital clinical significance for controlling the pandemic.

Nevertheless, to our knowledge, the scope of studies investigating the evolution of CT findings in COVID-19 is limited [4,5,6,7,8,9]. For example, Wang et al. [5], Zhou et al. [6], and Pan et al. [8] focused on the evolution of imaging performance within 14 or 15 days after the onset of symptoms, while more than 14 or 15 days after illness onset was classified as the last stage. However, the time from illness onset to ICU admission is 12 (IQR 8, 15) days [10]. Therefore, focusing on 14 or 15 days after illness onset as the last stage may obscure a critical change in disease course. Furthermore, certain studies [6,7,8,9] analyzed patients with mild as well as severe pneumonia, which may bias the results. For example, Liang et al. [4] analyzed only the CT findings of patients who were not discharged from the hospital during the first 3 weeks after admission, and Shi et al. included only the initial CT date of all patients [7]. Therefore, a large-scale longitudinal study on the CT manifestations of patients with mild COVID-19 pneumonia is urgently required.

Here we describe the evolution of CT appearances in 182 patients with mild COVID-19 with definitive outcomes (all patients were discharged and none became severely or critically ill or died). We aimed to compare the CT findings across the disease courses and analyze the correlation between CT scores and inpatient days, which may be helpful for understanding the nature of disease and for optimally managing patients.

Methods

The study was approved by the Research Ethics Commission of Wuhan Third Hospital (Approved number KY2020-027), and informed consent was waived in accordance with the guidelines of the Council for International Organizations of Medical Sciences. The research study was performed in accordance with the Helsinki Declaration of 1964, and its later amendments.

In this retrospective, single-center, longitudinal study we included 182 patients (Fig. 1) with confirmed COVID-19 admitted to the Wuhan Third Hospital (Tongren Hospital of Wuhan University) from January 21, 2020 to March 15, 2020. We used a published method for identifying SARS-CoV-2 infection [11]. The standards for discharge were absence of fever for more than 3 days, significant improvement in respiratory symptoms, obvious absorption of lesions in both lungs on CT, and negative SARS-CoV-2 RNA status in two consecutive nucleic acid tests (at least 1 day apart) of throat swabs [12]. In the present study, in accordance with WHO guidelines, mild COVID-19 pneumonia of adults was defined as COVID-19 pneumonia with no signs of severe pneumonia and no requirement for supplemental oxygen [13]. The inclusion criteria were as follows: (1) patients diagnosed with asymptomatic, mild, or moderate-type illness according to Chinese guidelines [12] who did not convert to the severe or critical type before discharge, (2) at least one positive CT, (3) recovery and discharge, and (4) age greater than 18 years.

Fig. 1
figure 1

Screening process of 182 patients with mild coronavirus disease 2019

CT

Chest CT images without contrast were acquired using a single inspiratory phase using two CT scanners (SOMATOM Definition AS, Siemens Healthineers, Erlangen, Germany; uCT 760, United Imaging, Shanghai, China). All patients were supine when examined. The CT protocol was as follows: 120 kVp; 120–200 mA; slice thickness 5–10 mm; matrix 512 × 512; pitch 0.625. The reconstruction section thickness was 0.5–1 mm. We collected CT scans of patients from admission to discharge.

Image Analysis

Two radiologists with 15 and 10 years of experience in chest imaging, respectively, evaluated the CT images to reach a consensus diagnosis. Definitions of CT imaging features referred to in the Fleischner Society glossary and peer-reviewed literature on viral pneumonia [7, 14,15,16]. The observers reviewed the image features as follows: lung involvement, extent of lesion involvement, predominant location, margin definition, distribution, lung segments of lesion distribution, number of lung segments and lobes involved, predominant CT pattern (Fig. 2), pure ground-glass opacity (GGO), pure consolidation, GGO and consolidation, linear opacity, adjacent pleura thickening, nodules, round cystic changes, bronchiectasis, air bronchogram sign, bronchial wall thickening, interlobular septal thickening, crazy-paving pattern, pleural effusion, thoracic lymphadenopathy, honeycomb pattern, tree-in-bud, reversed halo sign, calcification, and cavitation. Observers referred to a severity scoring system [16].

Fig. 2
figure 2

Predominant CT patterns of mild COVID-19 pneumonia. a GGO. b Consolidation. c GGO and consolidation. d Linear opacity. e Linear opacity, GGO, and consolidation. f Crazy-paving

Each lung was divided into upper, lower, and middle zones according to the horizontal lines through the carina and inferior pulmonary vein (six zones). The scoring rules for the involvement of each zone were as follows: score 0 = 0%; score 1, less than 25% involvement; score 2, 25% to less than 50%; score 3, 50% to less than 75%; and score 4, 75% or higher. The total severity score represents the sum of the scores of the six zones (possible scores, 0–24). The total lesion score represents the total score of involved lesions. We further evaluated the individual total scores for GGO, consolidation, and linear opacity of all lesions. To obtain detailed information about the dynamic evolution of CT features with disease progression, we divided the time axis into 10 periods according to the decile interval of the time from the onset of illness: stage 1 (1–5 days), stage 2 (6–8 days), stage 3 (9–11 days), stage 4 (12–14 days), stage 5 (15–16 days), stage 6 (17–18 days), stage 7 (19–21 days), stage 8 (22–25 days), stage 9 (26–30 days), and stage 10 (more than 30 days). The date of disease onset was defined as patients’ reported date of symptom onset. Moreover, studies evaluating the therapeutic effect of the COVID-19 typically required 1 week [17,18,19]. Therefore, to provide a reference for subsequent research on the evaluation of the therapeutic effects, we defined each week as one stage to further explore the dynamic changes in CT scores and imaging performance. The time axis was accordingly divided as follows: weeks 1, 2, 3, 4, and after 4 weeks. Finally, we analyzed the correlation between inpatient days and weekly CT scores.

Statistical Analysis

Statistical analyses were performed using IBM SPSS Statistics Software (version 25; IBM, New York, USA). Continuous variables are expressed as the mean and standard deviation (SD) or median and interquartile ranges (IQRs), as appropriate. Categorical variables are summarized as counts and/or percentages. The Kruskal–Wallis test was used to evaluate the differences among continuous variables of each stage, and a trend test was performed using a polynomial contrast procedure. The chi-square test was used to evaluate the differences between categorical variables in each stage, and the trend test was used with the chi-square test or curve fitting. Correlation coefficients were then calculated between inpatient days and CT scores using Spearman correlation analysis. P values less than 0.05 (two-tailed) were considered significant, and the comparison between specific pairwise comparisons used the adjusted P value.

Results

Patients’ Demographic and Clinical Characteristics

In this research, 182 patients with 606 CT scans were ultimately enrolled. Patients’ median age was 58 years (IQR 45, 66; Table 1). The most common symptoms were fever (151, 83%) and cough (130, 71%). All patients met the discharge criteria, and the lung lesions of 3 (1.6%) patients were completely absorbed upon discharge.

Table 1 Clinical characteristics of 182 patients with mild coronavirus disease 2019

CT Scores

Of the 606 CT scans (CTs), 17 were negative, and we therefore focused on the lesion scores and characteristics of the positive CTs.

In the dynamic evolution of the 10 stages, the numbers of CTs of each stage were 53, 62, 54, 70, 55, 43, 64, 72, 59, and 57, respectively. There were significant differences between different stages in scores (Fig. 3a). The differences between specific pairwise comparisons are shown in Fig. 3b. We found two peaks in the changing trends of the total lesion scores: stage 3 (9–11 days) and stage 6 (17–18 days), scores = 7.19 ± 3.66 and 8.00 ± 4.57, respectively.

Fig. 3
figure 3

Dynamic changes of CT scores and numbers of lung segments and lung lobes involved during the 10 stages described in “Methods”. a The scores for total lesion, consolidation, linear opacity, and GGO of each stage and the significant results of trend tests are shown in the table below. b Statistical analysis of pairwise comparisons of CT scores of 10 stages. c The numbers of lung segments and lung lobes involved during the 10 stages and significant results of trend tests are shown in the table below. d Pairwise comparison of the numbers of involved lung lobes and lung segments during the 10 stages. Stage 1, 1–5 days; stage 2, 6–8 days; stage 3, 9–11 days; stage 4, 12–14 days; stage 5, 15–16 days; stage 6, 17–18 days; stage 7, 19–21 days; stage 8, 22–25 days; stage 9, 26–30 days; stage 10, more than 30 days. GGO ground glass opacity; *P < 0.05; **P < 0.01

The three peaks of the GGO score occurred during stage 3 (9–11 days), stage 6 (17–18 days), and stage 8 (22–25 days), with scores of 5.56 ± 3.16, 5.86 ± 3.70, and 5.22 ± 3.93, respectively. The consolidation score peaked during stage 6 (17–18 days, score = 2.72 ± 3.07). The linear opacity score peaked during stage 6 (score = 3.81 ± 3.35). The numbers of lung segments and lobe involvement peaked during stage 6 (scores = 12.93 ± 4.61 and 4.35 ± 0.97, respectively), conforming to a quadratic function (Fig. 3c). Furthermore, there was a significant difference in the number of involved lung segments during different stages. Specific pairwise comparisons are shown in Fig. 3d.

In contrast, when a 1-week interval was used and the time axis was divided into five stages, the numbers of CTs in each stage were 103, 136, 162, 109, and 79. In the dynamic evolution of five stages, there were significant differences in total lesions (P < 0.05), consolidation (P < 0.01), and linear opacity (P < 0.01) scores (Fig. 4a). The differences between the specific pairwise comparisons are shown in Fig. 4b. The total lesion score, GGO score, consolidation score, and linear opacity score peaked during week 3 (7.3 ± 4.15), week 2 (4.98 ± 3.13), week 2 (2.54 ± 3.25), and week 3 (2.54 ± 2.46), respectively. The total lesion score and GGO score showed significantly positive correlation with the inpatient days  4 weeks after the onset of the disease (r range 0.26–0.31, P < 0.05) (Table 2). The number of lung segments and the extent of lobe involvement remained high during weeks 2–4, and the curve conformed to a quadratic function (Fig. 4c). Further, there were significant differences in the numbers of lung segments and lung lobes involvement during each stage. The specific pairwise comparisons are shown in Fig. 4d.

Fig. 4
figure 4

Dynamic changes of CT scores and the numbers of lung segments and lung lobes involved in five stages. a The total lesion, consolidation, linear opacity, and GGO scores during each stage. Significant results of trend tests are shown in the table below. b Statistical analysis of pairwise comparisons of CT scores during five stages. c The numbers of lung segments and lung lobes involved during stages and significant results of trend tests are shown in the table below. d Pairwise comparison of the numbers of involved lung lobes and lung segments during 10 stages. GGO ground glass opacity; *P < 0.05; **P < 0.01

Table 2 Correlation between inpatient days and CT scores

Imaging Features

Image features are summarized in Fig. 4. In each stage, lesions, which were mainly bilateral (86–94%), multiple (76–83%), and ill-defined (46–61%), were located in the subpleural (78–85%) and peripheral areas (78–85%) of both lungs. The percentages of the well-defined margins of the lesions conformed to a quadratic function (Fig. 5a). The lesions more frequently involved the apical posterior (54–76%), superior (68–87%), posterior basal (68–80%), and lateral basal segments (65–78%) of the left lung; and the posterior (59–75%), lateral (55–73%), superior (68–88%), posterior basal (68–83%), and lateral basal segments (65–82%) of the right lung (Fig. 5b, c). The predominant CT patterns differed significantly during each stage (P < 0.01), and GGO was the most common pattern in each stage (33–46%) (Fig. 5d). The percentages of the GGO and consolidation opacity exhibited a linearly decreasing trend (P < 0.001) (Fig. 5e).

Fig. 5
figure 5

Imaging manifestations of 182 patients during five stages. a Lung involvement, extent of lesion involvement, margin definition, predominant location, and distribution of the lesions during each stage. b Proportion of involvement of each lung segment of the left lung. c Proportion of involvement of each lung segment of the right lung. d Percentages of predominant patterns in each stage and significant results of trend tests are shown in the table below. e Percentages of pure ground glass opacity (GGO), pure consolidation, GGO and consolidation and linear opacity, and significant results of trend tests are shown in the table below. f The percentages during each stage of thickening pleura, nodules, round cystic changes, bronchiolectasis, air bronchogram sign, bronchial wall thickening, interlobular septal thickening, and crazy-paving pattern. Significant results of trend tests are shown in the table below. g Percentages during each stage of pleural effusion, lymphadenopathy, honeycomb pattern, tree-in-bud and reversed halo sign. Significant results of trend tests are shown in the table below. GGO ground glass opacity; *P < 0.05; **P < 0.01

The proportion of the pure GGO maintained a high value during each stage (72–82%), peaking during stage 2 (week 2) and stage 5 (after 4 weeks). The changing trends of the percentages of pure consolidation and linear opacity were similar to the changing trends of the consolidation and linear opacity scores, and the percentages of both peaked during week 3. Among other imaging signs, pleural thickening (56–74%), interlobular septal thickening (26–48%), and the crazy-paving pattern (15–35%) were relatively common, while lymphadenopathy (5–9%), pleural effusion (0–5%), and reversed halo signs (0–3%) were relatively infrequent (Fig. 5f, g).

Discussion

To the best of our knowledge, of all the longitudinal CT performance studies of COVID-19 pneumonia [4,5,6,7,8], the present study analyzed the largest number of patients. Through the study of 589 COVID-19 CT images of 182 patients, we found that the total lesion CT score proposed by others only considers the area of involvement of the lesions, but not the nature of the lesion [4, 6, 8, 20,21,22], which may lead to delays of imaging evaluation in clinical practice [19]. Therefore, according to the nature and area of the lesions, we first evaluated consolidation, linear opacity, and GGO scores. Second, although the total lesion score peaked during days 9–11, similar to previous reports [6, 8, 9], we found it interesting that the peak of the total lesion score was highest during days 17 and 18, because this is not reported by others. Third, the peak times of the four CT scores were basically consistent during the 10 stages. However, if the interval was extended to 1 week, the peak time of the total CT score lagged behind those of the consolidation and GGO scores, suggesting that a useful scoring system should include comprehensive assessment of the extent and histopathology of the lesion.

The peak times of the total lesion score differ among previous studies [5, 6, 8, 9]. These studies define the last stage of pneumonia as days 14 or 15 after disease onset. However, the average time from illness onset to discharge reported by other investigators is approximately 28 days [23, 24], and more than 12 days after illness onset is a critical period for disease progression [10]. Therefore, we speculate that the staging method of these studies reported the average of the true peak of the total lesion score. Clearly, in our research, the 10-stage curve confirmed our prediction that the total lesion score peaked during days 9–11 as well as on days 17–18. The latter peak corresponds to week 3 of the five-stage curve, while the previous peak was averaged. In addition, the peak time of the consolidation score was earlier than that of the total score of the five-stage change, and the curve of the consolidation score conformed to a quadratic function. Moreover, the decreased time of the total lesion score is later than the improvement time of clinical parameters [19], and greater consolidation may indicate disease progression [25]. Therefore, the consolidation score may reflect the severity of disease more promptly than the total lesion score, although the changes of the two scores may tend to be associated with the increase of the number of stages.

Among the typical features of each stage, here we show that COVID-19 pneumonia was mainly characterized by bilateral, multiple, and ill-defined lesions, with a GGO pattern, located in the posterior and peripheral regions of both lungs, similar to previous findings [5,6,7, 9]. In the predominant CT pattern of COVID-19 pneumonia, the GGO and linear opacity patterns increased rapidly from week 3 after disease onset. Histologically, GGO represents interstitial mononuclear inflammatory infiltrates and edema [26], while during the absorption stage, the absorption of exudate cells and hyaline membrane in the lung tissue forms other lesions such as consolidation and the crazy-paving pattern, gradually transforming into GGO. Evidence indicates that linear opacity is an organized feature during the late stage of disease [27]. Therefore, good prognosis may be associated with increases in the proportions of the GGO and linear opacity patterns during the later stage of disease. However, the increases in the total score and GGO score during stage 5 (after week 4) may indicate an extension of hospital stay. Furthermore, the proportions of the consolidation pattern and the consolidation score peaked during week 2 after disease onset and gradually decreased, suggesting that the events that occur during week 2 may be critical for disease progression.

In the present study, the changing trends of the proportions of consolidation, GGO, and linear opacity were not completely consistent with the changing trends of their respective scores, suggesting that previous studies [7, 27], based only on the proportion of lesions used to judge the severity of the disease, may have generated unreliable results. Thus, the proportion of lesions and the area involved should therefore be simultaneously considered. Here, pleural thickening and interlobular septal thickening were relatively common CT signs during each stage. The crazy-paving pattern sign was relatively common during the first 2 weeks after illness onset, and pleural effusion appeared during the first 3 weeks.

Obviously, CT plays an important role in controlling the spread of COVID-19. However, CT has some insurmountable shortcomings, such as high cost, exposure to radiation, and their immobile nature, which can be overcome by lung ultrasound (LUS) [28,29,30,31]. Coupled with the fact that the lesions of COVID-19 tend to occur in subpleural regions, LUS has been suggested as a potential triage and diagnostic tool for COVID-19 [32]. The latest research showed that LUS was a promising tool for early risk stratification in COVID-19, which may guide patients’ management strategies, as well as resource allocation in case of surge capacity [28, 31].

There were some limitations to our study. First, this was a single-center study, and the results require verification through a multicenter study. Second, the CT scores of this study may not be comparable with the quantitative scores determined using artificial intelligence [33], although visual scores are easier to determine. Finally, only patients with mild COVID-19 pneumonia were included, and the longitudinal CT performance of severe pneumonia must be determined using a larger sample size.

Conclusions

The total lesion score of mild COVID-19 pneumonia peaked 17–18 days after illness onset, and week 2 after onset may be the critical period for disease progression. The consolidation scores may objectively reflect the severity of lung involvement and thus serve as a more accurate diagnostic measure than total lesion scores. The mild COVID-19 pneumonia in each stage was mainly manifested as bilateral, multiple, and ill-defined lesions, with the GGO pattern, and located in the posterior and peripheral areas of both lungs. Moreover, good prognosis may be associated with increases in the proportions of the GGO and linear opacity patterns during the later stage of disease.