key: cord-0830730-5925lvue
authors: Pan, Feng; Zheng, Chuansheng; Ye, Tianhe; Li, Lingli; Liu, Dehan; Li, Lin; Hesketh, Richard L.; Yang, Lian
title: Different computed tomography patterns of Coronavirus Disease 2019 (COVID-19) between survivors and non-survivors
date: 2020-07-09
journal: Sci Rep
DOI: 10.1038/s41598-020-68057-4
sha: 7966b90527495e5d26e11269414c2d2796e5ba59
doc_id: 830730
cord_uid: 5925lvue

This study aimed to compare the chest computed tomography (CT) findings between survivors and non-survivors with Coronavirus Disease 2019 (COVID-19). Between 12 January 2020 and 20 February 2020, the records of 124 consecutive patients diagnosed with COVID-19 were retrospectively reviewed and divided into survivor (83/124) and non-survivor (41/124) groups. Chest CT findings were qualitatively compared on admission and serial chest CT scans were semi-quantitively evaluated between two groups using curve estimations. On admission, significantly more bilateral (97.6% vs. 73.5%, p = 0.001) and diffuse lesions (39.0% vs. 8.4%, p < 0.001) with higher total CT score (median 10 vs. 4, p < 0.001) were observed in non-survivor group compared with survivor group. Besides, crazy-paving pattern was more predominant in non-survivor group than survivor group (39.0% vs. 12.0%, p < 0.001). From the prediction of curve estimation, in survivor group total CT score increased in the first 20 days reaching a peak of 6 points and then gradually decreased for more than other 40 days (R(2) = 0.545, p < 0.001). In non-survivor group, total CT score rapidly increased over 10 points in the first 10 days and gradually increased afterwards until ARDS occurred with following death events (R(2) = 0.711, p < 0.001). In conclusion, persistent progression with predominant crazy-paving pattern was the major manifestation of COVID-19 in non-survivors. Understanding this CT feature could help the clinical physician to predict the prognosis of the patients.

www.nature.com/scientificreports/ chain reaction (RT-PCR) test has a relatively high false-negative rate (29%) for COVID-2019 diagnosis, so chest computed tomography (CT) is recommended as the major screen modality with a higher sensitivity of 97% and faster performance [10] [11] [12] [13] . In Hubei province, the centre of the outbreak in China, the clinical diagnostic criteria were only dependent on chest CT scan, instead of the RT-PCR test before 19 February 2020 14 . However, the value of the consecutive CT scans for monitoring disease progression was still unclear. Previous studies suggested a typical time course of CT findings in survivors with COVID-19, in which initial progression was followed by recovery, the latter starting after about 2 weeks [15] [16] [17] . Case series have associated severe and critical COVID-19 with more diffuse lung involvement, development of acute respiratory distress syndrome (ARDS), and multi-organ failure 7, [18] [19] [20] . Using a case-control design, this study aims to identify the differentiating CT features and compare the temporal evolution of pulmonary involvement between recovered and died patients with COVID-19.

Patients and groups. 175 consecutive records of hospitalized patients with RT-PCR confirmed COVID- 19 were reviewed retrospectively for the period from 12 January 2020 to 20 February 2020 in this single-centre (Union Hospital, Wuhan, China). The inclusion criteria included: (1) with definite clinical outcomes (discharge or death events); (2) no comorbidities which might impair the immune or pulmonary function, such as recent chemotherapy and chronic obstructive pulmonary disease; (3) with more than three times of chest CT scans in the course for sufficient estimation of radiological patterns, unless fatal ARDS occurred resulting in impossibility to carry out the consecutive chest CT scans. Eventually, 124 patients were included and divided into two groups: survivor group (discharged patients, n = 83, including 21 patients who were preliminarily reported in the previous study 15 ) and non-survivor group (died patients, n = 41) (Fig. 1) . Clinical data (e.g. initial symptoms, past medical history, etc.) and serial chest-CT data in the follow-up (extended until 30 March 2020 in survivor www.nature.com/scientificreports/ group) were retrieved through the institutional electronic patient database. Diagnostic, isolation, grades of the disease severities (non-ARDS and ARDS), treatment, and discharge criteria were based on the published standard protocols from the continuously-updated National Health Commission of the People's Republic of China 14 .

Chest CT scan protocols. Chest CT scans were performed using two commercial multi-detector CT scanners (Philips Ingenuity Core128, Philips Medical Systems, Best, Netherlands; SOMATOM Definition AS, Siemens Healthineers, Erlangen, Germany) during a single breath-hold. The low-dose mode was set up with a tube voltage of 120 kVp and automatic tube current modulation. From the raw data, CT images were reconstructed as 1.5 mm thick axial slices and increment of 1.5 mm in transverse slice orientation with either hybrid iterative reconstruction (iDose level 5, Philips Healthcare, Best, The Netherlands) or a pulmonary B70F kernel (Siemens Healthineers, Erlangen, Germany).

Chest CT estimation. The abnormal radiological findings of CT reported using internationally standard nomenclature [21] [22] [23] . CT abnormalities included ground-glass opacity (GGO), crazy-paving pattern, and consolidation. The distribution of abnormalities was also noted as being predominantly subpleural (involving mainly the subpleural one-third of the lung), random (without predilection for subpleural or central regions), or diffuse (continuous involvement without respect to lung segments) 24 . A conventional semi-quantitative scoring system was used to evaluate the pulmonary involvement area of all these abnormalities 15, 25 . There was a score of 0-5 for each lobe on the following: 0-no involvement; 1, < 5% involvement; 2, 6-25% involvement; 3, 26-49% involvement; 4, 50-75% involvement; 5, > 75% involvement. The total CT score was the sum of the score of each lobe and ranged from 0 (no involvement) to 25 (maximum involvement). The analysis was performed using the institutional digital database system (Vue PACS, version 11.3.5.8902, Carestream Health, Oakville, Canada) by two radiologists (CZ and LY, who had 26 and 22 years of experience in thoracic radiology, respectively) and the decisions were reached in consensus. All radiologists were blinded to the groups and clinical progress of the patients to avoid information bias.

Statistical analysis. Statistical analysis was performed using IBM SPSS Statistics Software (version 24;

IBM, New York, USA). Quantitative data were presented as median with inter-quartile range (IQR) and qualitative data were presented as the percentage of the total unless otherwise specified. The comparisons of the quantitative data were statistically evaluated using the Mann-Whitney U test, according to the non-normal distribution assessed by the Shapiro-Wilk test. The comparisons of qualitative data were evaluated using the Chi-square test or Fisher's exact test. The dynamic total CT score with time from symptom onset was quantitatively assessed by using the SPSS curve estimation module 15 . A p value of < 0.05 was defined as having statistical significance.

Basic characteristics. The median age of the patients was 56 years (IQR 38-68 years) with an approximately equal male to female ratio (63:61). The median age of patients was significantly higher in non-survivor compared to non-survivors (69 years vs. 43 years, p < 0.001). The percentage of males was 38.6% and 75.6% in survivor and non-survivor groups, respectively (p < 0.001) ( Table 1) . Non-survivors were also more likely to have a history of hypertension, diabetes, and coronary heart disease than survivors (p < 0.05) ( Table 1) . Fever and cough were the most common initial symptoms (85.5% and 65.3%, respectively). Chest distress was significantly more inclined to occur in non-survivors (p < 0.001) ( Table 1 ). There was no significant difference in the period of admission from symptom onset between survivor and non-survivor groups (8 days vs. 9 days, p = 0.422) ( Table 1 ). The median survival period of non-survivor group after admission was 14 days (IQR 8-22 days) from admission, while the median hospitalized period in survivor group was 18 days (IQR 12-27 days) (p = 0.068). The survivors underwent more times of chest CT scans than non-survivors (4 vs. 2, p < 0.001) with a significantly longer duration (6 days vs. 5 days, p = 0.001) ( Table 1 ). All non-survivors aggravated to ARDS after a median of 11 days (8-14 days) from symptom onset, while only one patient aggravated to ARDS in survivor group. Multiple biochemical and haematological parameters differed significantly between the two groups such as lymphocyte count, neutrophil count, and C-reactive protein (CRP) (p < 0.05) ( Table 2 ).

Comparison of major CT findings between two groups. All 124 patients underwent a total of 436 chest CT scans with a median interval between adjacent scans of 6 days (IQR 5-12 days) ( On admission, bilateral lung involvement was more common in non-survivors than survivors (97.6% vs. 73.5%, p = 0.001) ( Table 3) . Subpleural distribution was more inclined to be observed in survivors compared with non-survivor group (69.9% vs. 43.9%, p = 0.005), while diffuse distribution was more common in nonsurvivor group compared with survivor group (39.0% vs. 8.4%) (p < 0.001) ( Table 3) . GGO (83.1%), consolidation (60.5%), and crazy-paving pattern (42.7%) were the major CT findings in both groups, while the crazy-paving pattern was more common in non-survivor group than survivor group (65.9% vs. 31.3%, p < 0.001) ( Table 3) . On admission, consolidation predominated in survivor group (37.3%), but crazy-paving pattern predominated in non-survivor group compared with survivor group (39.0% vs. 12.0%, p = 0.001) ( Table 3) . Besides, the total CT score was significantly higher in non-survivor group than survivor group (a median of 10 vs 4, p < 0.001).

Dynamic estimation of pulmonary involvement between two groups. Based Interval between adjacent chest CT scans (days) (IQR) 6 (5-12) 6 (5-13) 5 (3-9) 0.001 www.nature.com/scientificreports/ respectively; p < 0.001, each) (Fig. 2a,b ; SI Table S1 online). The optimal fitting equations were demonstrated in Fig. 2c . From the optimal fitting, in survivor group the total CT score gradually increased in the first 20 days with a peak value of 6 and then gradually decreased afterwards lasting for more than another 40 days (Fig. 2c) . The typical CT manifestation was changed from subpleural GGO to enlarged consolidation with time which was gradually absorbed afterwards leaving residual GGO and parenchymal bands (Fig. 3) . But in non-survivor group, the total CT score rapidly increased in the first 10 days and eventually approached 15 until ARDS occurred (Fig. 2c) . From the dynamic CT images, the persistently progressive pulmonary lesions from GGO with crazy-paving pattern to bilaterally extensive consolidation could be observed (Fig. 4 ).

This study compared the temporal changes in CT manifestations between survivors and non-survivors with COVID-19. It demonstrated the pulmonary involvement of subpleural GGO and sequential consolidation gradually progressed reaching the peak after 20 days since symptom onset in survivors. Afterwards, the lesions started to be absorbed lasting for more than 40 days. In contrast, non-survivors demonstrated more rapid and persistent www.nature.com/scientificreports/ progression with more extensive bilateral lesions until ARDS occurred. Crazy-paving pattern was more predominant in non-survivors on admission compared with survivors. In accordance with the previous studies, old patients (69 years, IQR 63-78 years) with more comorbidities such as hypertension, diabetes, and coronary heart disease were more inclined to develop fatal ARDS 4, 8 . Initial symptoms were similar between survivor and non-survivor groups, whilst chest distress was more common in non-survivor group. Patients in non-survivor group underwent a progressive phase which culminated in the development of ARDS after a median period of three weeks from symptom onset. As a case-controlled study, the mortality rate of ARDS caused by COVID-19 could not be evaluated, but from a previous study, it was reported mortality of 61.5% 26 .

Initial laboratory investigations on admission showed multiple haematological and biochemical abnormalities which were significantly different between survivor and non-survivor groups. This can be attributed to the systematic inflammation reaction and pulmonary vascular endothelial damage caused by a severe viral infection, similar to the systemic response seen in other types of severe pneumonia 8,24-28 . It has been postulated that COVID-19 could also damage T lymphocytes, thus, significant lymphopenia was probably a risk factor leading to the deterioration of patients' immune function and more rapid disease progression 7, 8, 26, 29 . In addition, the increased levels of CRP, lactate dehydrogenase, and d-dimer could also be indicators for development of ARDS, as reported in other types of pneumonia 23, 25, 27, 30 .

In the early stage of COVID-19, subpleural GGO was the predominant finding [15] [16] [17] 20, 31 . But in this study, patients were hospitalized after a median period of 8 and 9 days after the onset of symptoms in survivor and non-survivor groups, respectively, at which time the predominant findings in both groups corresponded with the progressive stage 15 . Thus, GGO was not the predominant finding in both groups but the consolidation and www.nature.com/scientificreports/ crazy-paving pattern. Compared with the survivors, it demonstrated the predominant CT demonstration of crazy-paving pattern in non-survivor group on admission was a major difference except for more diffuse and bilateral distributions. Pathologically, GGO may be an indicator of alveolar oedema and proteinaceous exudates 32 .

As the disease progresses, increasing alveolar oedema, exudates, and lymphocyte infiltrates fill the interstitial space leading to the radiological demonstration of diffuse "crazy-paving pattern" 22, 28, 33, 34 . Subsequent ARDS and potentially fatal respiratory failure developed as a result of diffuse alveolar oedema with loss of alveolar epithelium 22, 34 . Thus, it was speculated large area of crazy-paving pattern was probably a CT indicator of poor prognosis.

Considering the heterogeneities of the scan time among the patients, longitudinal comparisons were not appropriate. Thus, the curve estimation was used to statistically compare the temporal evolution of the disease between two groups. Being different from the static comparison of chest CT on admission using the logistic module, curve estimation could analyze the dynamic patterns of the pulmonary involvement with time 19, 35 . Thus, it could provide a more composite comprehension of the time course in COVID-19 between survivors and non-survivors. As a result, it demonstrated a gradual resolution of abnormalities after a maximal total CT score of 6 at 20 days, longer than 10 days reported in the previous report 15 . It might be ascribed to a limited sample size in the previous study, which probably underestimated the recovery duration of COVID-19. Compared with survivor groups, the total CT score in non-survivor group demonstrated a more rapid increase in the first 10 days with a higher value of more than 10 points. Although the previous study showed the feasibility of making CT score as an indicator of prognosis, it did not demonstrate the dynamic changes of CT score in the whole course 19 . This study revealed the total score persistently elevated to a higher level close to 15 points www.nature.com/scientificreports/ without any decrease in non-survivor group, until the ARDS occurred with the following death events. From one pathological study in severe acute respiratory syndrome (SARS), it found the long duration of illness resulted from the severe fibrosis and organization 28 . Considering the partial homology of SARS and COVID-19, it might explain why the lesions were rarely absorbed in non-survivors with COVID-19. This is another major difference between the two groups in the course, associated with the refractory feature of the critical COVID-19 under the present treatment protocols 14 . This study has limitations. Firstly, as a retrospective study, chest CT was used by the physician based on the clinical necessity and the status of the patient, so the heterogeneities of scanning time made it impossible to perform a conventional longitudinal comparison between two groups with regular intervals. Second, CT was not clinically feasible for patients after developing ARDS so not enough CT information was provided in the course of ARDS. Consequently, the majority of CT scans were performed in mild disease (363/436, 83.3%). To avoid data heterogeneity, the comparison of chest CT between two groups was only performed on admission due to a similar period from symptom onset and the curve estimation was used to evaluate the comprehensive trend of pulmonary involvement between two groups. Third, the multi-variate regression or propensity matching involving the CT, clinical, and laboratory parameters was not performed owing to the limited sample size and a relatively large number of parameters with significant differences between the two groups.

In summary, from comparisons between survivors and non-survivors, this study indicated that the presence of predominant crazy-paving pattern on chest CT with the high and rapidly increased CT scores may help to identify the patients at high risk of developing ARDS before clinical deterioration. A larger, prospective study is required to confirm these findings with the more accurate quantitative assessment modality of the CT images in COVID-19.

Ethical approval. This retrospective study was approved by the Ethics of Committees of Union Hospital, Tongji Medical College, Huazhong University of Science and Technology (No. 2020-0026), and followed the 1964 Helsinki Declaration and its later amendments or comparable ethical standards.

Patient and other consents. Informed consent/deceased patient permission form for this retrospective study was waived by Ethics of Committees of Union Hospital, Tongji Medical College, Huazhong University of Science and Technology because only the anonymous data was collected and analyzed to facilitate better clinical decisions and treatment.

Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia

Clinical features of patients infected with 2019 novel coronavirus in Wuhan

Situation Report 137 (Geneva: World Health Organization

Clinical characteristics of coronavirus disease 2019 in China

A Novel coronavirus from patients with pneumonia in China

Good ACE, bad ACE do battle in lung injury

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study

Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan

Clinical findings in a group of patients infected with the 2019 novel coronavirus (SARS-Cov-2) outside of Wuhan, China: Retrospective case series

Chest CT for typical 2019-nCoV pneumonia: Relationship to negative RT-PCR testing

Sensitivity of chest CT for COVID-19: Comparison to RT-PCR

Correlation of chest CT and RT-PCR testing in coronavirus disease 2019 (COVID-19) in China: A report of 1014 cases

Diagnosis of the coronavirus disease (COVID-19): rRT-PCR or CT?

National Health Commission of the People's Republic of China. Diagnosis and Treatment Protocols of Pneumonia Caused by a Novel Coronavirus

Time course of lung changes at chest CT during recovery from coronavirus disease 2019 (COVID-19)

Coronavirus disease (COVID-19): Spectrum of CT findings and temporal progression of the disease

novel coronavirus (COVID-19) pneumonia: Serial computed tomography findings

Emerging 2019 novel coronavirus (2019-nCoV) pneumonia

The clinical and chest CT features associated with severe and critical COVID-19 pneumonia

Radiological findings from 81 patients with COVID-19 pneumonia in Wuhan, China: A descriptive study

Fleischner Society: Glossary of terms for thoracic imaging

Imaging of pulmonary viral pneumonia

Radiographic and CT features of viral pneumonia

Severe acute respiratory syndrome: Temporal lung changes at thin-section CT in 30 patients

Pulmonary sequelae in convalescent patients after severe acute respiratory syndrome: Evaluation with thinsection CT

Clinical course and outcomes of critically ill patients with SARS-CoV-2 pneumonia in Wuhan, China: A singlecentered, retrospective, observational study

The association between plasma D-dimer levels and community-acquired pneumonia

Pulmonary pathology of severe acute respiratory syndrome in Toronto

T-cell immunity of SARS-CoV: Implications for vaccine development against MERS-CoV

Analysis of clinical features of 29 patients with 2019 novel coronavirus pneumonia

Novel coronavirus pneumonia outbreak in 2019: Computed tomographic findings in two cases

Pulmonary pathology of early-phase 2019 novel coronavirus (COVID-19) pneumonia in two patients with lung cancer

Crazy-paving appearance at thin-section CT: Spectrum of disease and pathologic findings

Pathological findings of COVID-19 associated with acute respiratory distress syndrome

Initial CT findings and temporal changes in patients with the novel coronavirus pneumonia (2019-nCoV): A study of 63 patients in Wuhan, China

We express our sincere gratitude to the emergency services, nurses, doctors, and a lot of medical supports from other provinces for their efforts to combat the COVID-19 outbreak in Wuhan. Besides, we are grateful to Dr Jiazheng Wang (MSC Clinical & Technical Solutions, Philips Healthcare) and Dr Dandan Zheng (MSC Clinical & Technical Solutions, Philips Healthcare) for many useful discussions through the formation and design of this work.

The datasets generated in the current study are available from the corresponding author on request.Received: 21 April 2020; Accepted: 17 June 2020

The corresponding author declares that all authors have read the manuscript and they all gave permission to submit the work in its current version. F.P. and C.Z. contributed equally to this article as the co-first authors. All authors have made substantial contributions to the conception and the design of the study, and/or acquisition of data, and/or analysis and interpretation of data, and the drafting of the article and its revision.

The authors declare no competing interests.

Supplementary information is available for this paper at https ://doi.org/10.1038/s4159 8-020-68057 -4.Correspondence and requests for materials should be addressed to L.Y.Reprints and permissions information is available at www.nature.com/reprints.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creat iveco mmons .org/licen ses/by/4.0/.