PMC:7160614 / 3988-4125 JSONTXT 4 Projects
A diagnostic model for coronavirus disease 2019 (COVID-19) based on radiological semantic and clinical features: a multi-center study
Abstract
Objectives
Rapid and accurate diagnosis of coronavirus disease 2019 (COVID-19) is critical during the epidemic. We aim to identify differences in CT imaging and clinical manifestations between pneumonia patients with and without COVID-19, and to develop and validate a diagnostic model for COVID-19 based on radiological semantic and clinical features alone.
Methods
A consecutive cohort of 70 COVID-19 and 66 non-COVID-19 pneumonia patients were retrospectively recruited from five institutions. Patients were divided into primary (n = 98) and validation (n = 38) cohorts. The chi-square test, Student’s t test, and Kruskal-Wallis H test were performed, comparing 1745 lesions and 67 features in the two groups. Three models were constructed using radiological semantic and clinical features through multivariate logistic regression. Diagnostic efficacies of developed models were quantified by receiver operating characteristic curve. Clinical usage was evaluated by decision curve analysis and nomogram.
Results
Eighteen radiological semantic features and seventeen clinical features were identified to be significantly different. Besides ground-glass opacities (p = 0.032) and consolidation (p = 0.001) in the lung periphery, the lesion size (1–3 cm) is also significant for the diagnosis of COVID-19 (p = 0.027). Lung score presents no significant difference (p = 0.417). Three diagnostic models achieved an area under the curve value as high as 0.986 (95% CI 0.966~1.000). The clinical and radiological semantic models provided a better diagnostic performance and more considerable net benefits.
Conclusions
Based on CT imaging and clinical manifestations alone, the pneumonia patients with and without COVID-19 can be distinguished. A model composed of radiological semantic and clinical features has an excellent performance for the diagnosis of COVID-19.
Key Points
• Based on CT imaging and clinical manifestations alone, the pneumonia patients with and without COVID-19 can be distinguished.
• A diagnostic model for COVID-19 was developed and validated using radiological semantic and clinical features, which had an area under the curve value of 0.986 (95% CI 0.966~1.000) and 0.936 (95% CI 0.866~1.000) in the primary and validation cohorts, respectively.
Electronic supplementary material
The online version of this article (10.1007/s00330-020-06829-2) contains supplementary material, which is available to authorized users.
Introduction
On January 30, 2020, the World Health Organization (WHO) has declared the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreak as a global health emergency of international concern. This outbreak has infected all provinces of China and rapidly spread to the rest of the world. At the time of writing this article (March 16, 2020), there have been more than 158 countries and territories affected [1]. Whole-genome sequencing and phylogenetic analysis reveal that the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is similar to some beta coronaviruses detected in bats, but it is distinct from severe acute respiratory syndrome coronavirus (SARS-Cov) and Middle East respiratory syndrome coronavirus (MERS-CoV) [2].
Patients with COVID-19 develop pneumonia with associated symptoms of fever (98%), cough (76%), and myalgia or fatigue (44%) [3]. CT imaging plays a critical role in the diagnosis and the monitoring of disease progression [4–6]. The latest research studies described the characteristic imaging manifestations of COVID-19, including ground-glass opacities (GGO) (57 to 88%), bilateral involvement (76 to 88%), and peripheral distribution (33 to 85%) [7–10]. Other imaging features such as consolidation, cavitation, and interlobular septal thickening are also reported in some patients [11–13]. However, these imaging manifestations of COVID-19 are nonspecific and are difficult to distinguish from other pneumonia. To our knowledge, there have been no studies explicitly comparing imaging and clinical characteristics between pneumonia patients with and without COVID-19.
The current diagnostic criterion for COVID-19 is the positive result of a nucleic acid test by real-time reverse transcription polymerase chain reaction (RT-PCR) or next-generation sequencing [14]. However, false-negative results caused by unstable specimen processing are relatively high in clinical practice, which has worsened the spread of the outbreak [15–18]. Moreover, laboratory testing for SARS-CoV-2 requires a rigorous platform, which is not assembled in all hospitals. Thus, this requires specimen transfer, which may delay diagnosis for days. Early and accurate diagnosis is crucial, particularly for critically ill patients who need emergency surgery, and with pneumonia complications. To solve these problems, we hypothesize that a diagnostic model can be developed based on CT imaging and clinical manifestations alone, independent of the nucleic acid test.
In this study, we identify the differences in imaging and clinical manifestations between patients with and without COVID-19. We also develop and validate a model for COVID-19 diagnosis based on radiological semantic and clinical features.
Patients and methods
Patients
Ethical approvals by the institutional review boards were obtained for this retrospective analysis, and the need to obtain informed consent was waived.
From January 1 to February 8, 2020, seventy consecutive patients with COVID-19 admitted in 5 independent hospitals from 4 cities were enrolled in this study (mean age, 42.9 years; range, 16–69 years), including 41 men (mean age, 41.8 years; range, 16–69 years) and 29 women (mean age, 44.5 years; range, 16–66 years). All patients were confirmed with SARS-CoV-2 infection by real-time RT-PCR and next-generation sequencing. Of these patients, 24 were from Huizhou City, 25 from Shantou City, 15 from Yongzhou City, and the rest 6 from Meizhou City. At the same period, another 66 pneumonia patients without COVID-19 from Meizhou People’s Hospital were recruited as controls (mean age, 46.7 years; range, 0.3–93 years), including 43 men (mean age, 46.0 years; range, 0.3–93 years) and 23 women (mean age, 48.0 years; range, 1–86 years). All the controls were confirmed with consecutive negative RT-PCR assays. Figure E1 in the Supplementary Material shows the patient recruitment pathway for the control group, along with the inclusion and exclusion criteria.
According to previous studies [19–21], whose sample size is comparable with ours, the ratio between primary and validation cohort is 7:3. In this study, a total of 136 patients were divided into primary (n = 98) and validation (n = 38) cohorts, close to 7:3. A total of 19 COVID-19 patients from two hospitals (6 patients from Meizhou People’s Hospital and 13 patients from the First Affiliated Hospital of Shantou University Medical College) and 19 randomly selected controls from Meizhou City were incorporated into the validation cohort. The rest of the patients are incorporated in the primary cohort, including 51 COVID-19 patients from Huizhou, Yongzhou, and Shantou cities and 47 controls from Meizhou City. The primary cohort was utilized to select the most valuable features and build the predictive model, and the validation cohort was used to evaluate and validate the performance of the model.
Image and clinical data collection
The chest CT imaging data without contrast material enhancement were obtained from multiple hospitals with different CT systems, including GE CT Discovery 750 HD (General Electric Company), SCENARIA 64 CT (Hitachi Medical), Philips Ingenuity CT (PHILIPS), and Siemens SOMATOM Definition AS (Siemens). All images were reconstructed into 1-mm slices with a slice gap of 0.8 mm. Detailed acquisition parameters were summarized in the Supplementary Material (Table E1).
The clinical history, nursing records, and laboratory findings were reviewed for all patients. Clinical characteristics, including demographic information, daily body temperature, blood pressure, heart rate, clinical symptoms, and history of exposure to epidemic centers, were collected. Total white blood cell (WBC) counts, lymphocyte counts, ratio of lymphocyte, neutrophil count, ratio of neutrophil, procalcitonin (PCT), C-reactive protein level (CRP), and erythrocyte sedimentation rate (ESR) were measured. All threshold values chosen for laboratory metrics were based on the normal ranges set by each individual hospital.
Image analysis
For extraction of radiological semantic features, two senior radiologists (D.L. and X.C., more than 15 years of experience) reached a consensus, blinded to clinical and laboratory findings. The radiological semantic features included both qualitative and quantitative imaging features. The lesions in the outer third of the lung were defined as peripheral, and lesions in the inner two-thirds of the lung were defined as central [22]. The progression of COVID-19 lesions within each lung lobe was evaluated by scoring each lobe from 0 to 4 [7], corresponding to normal, 1~25% infection, 26~50% infection, 51~75% infection, and more than 75% infection, respectively. The scores were combined for all five lobes to provide a total score ranging from 0 to 20. A total of 41 radiological features (26 quantitative and 15 qualitative) were extracted for the analysis. The descriptions of radiological semantic features are listed in the Supplementary Material (Table E2). Figure 1 is one example of the evaluation of CT imaging.
Fig. 1 A 23-year-old female with a travel history to Wuhan presenting with fever. Axial noncontrast CT image shows a consolidation with ground-glass opacities in the peripheral region by the right upper lobe. Air bronchogram is found in lesion. The maximum diameter of lesion is 2.8 cm. The right upper lobe score is 1 because of the involved lung parenchyma less than 1/4
Clinical and radiological feature selection
To obtain the most valuable clinical and radiological semantic features, statistical analysis, univariate analysis, and the least absolute shrinkage and selection operator (LASSO) method were performed. In statistical analysis, the chi-square test, the Kruskal-Wallis H test, and t test were utilized to compare the radiological semantic and clinical features between COVID-19 and non-COVID-19 groups. The features with p value smaller than 0.05 were selected. Then, univariate analysis was performed for clinical and radiological candidate features to determine the COVID-19 risk factors. The features with p value smaller than 0.05 in univariate analysis were also selected. The least absolute shrinkage and selection operator (LASSO) method [23] was utilized to select the most useful features with penalty parameter tuning that was conducted by 10-fold cross-validation based on minimum criteria. Diagnostic models were then constructed by multivariate logistic regression with the selected features. The flowchart of the feature selection process for these models was presented in the Supplementary Material (Fig. E2).
Development and validation of the diagnostic model
To develop an optimal model, we evaluated 3 models by analyzing (i) the clinical features model (C model), (ii) radiological semantic features model (R model), and (iii) the combination of clinical and radiological semantic features model (CR model) by multivariate logistic regression analysis. The classification performances of the models were evaluated by the area under the receiver operating characteristic (ROC) curve. The area under the curve (AUC), accuracy, sensitivity, and specificity were also calculated. A decision curve analysis was conducted to determine the clinical usefulness of the diagnostic model by quantifying the net benefits at different threshold probabilities in the validation dataset [24]. The development of decision curve was described in the Supplementary Materials. Figure 2 depicts the flowchart of the proposed analysis pipeline described above. We also built a nomogram, which was a quantitative tool to predict the individual probability of infection by COVID-19, based on the multivariate logistic analysis of the CR model with the primary cohort. Depending on the coefficient of the predictive factors in multivariate logistic regression model, all values of each predictive factor were assigned points. A total point was obtained by summing all the points of each predictive factor. The scale also showed the relationship between the total point and the prediction probability in the nomogram. The corresponding calibration curves of the CR model in the primary cohort and validation cohort are shown in the Supplementary Material (Fig. E3).
Fig. 2 Workflow of data process and analysis in this study. Radiological semantic features, including qualitative and quantitative imaging features, are extracted from axial lung CT section. The clinical manifestation and laboratory parameters are provided by electronic case system. Statistical analysis is performed for comparing the different features between COVID-19 and non-COVID-19 patients. Univariate analysis, least absolute shrinkage, and selection operator (LASSO) are further performed to determine the COVID-19 risk factors with p < 0.05 in statistical analysis. Three models based on the selected features are established by multivariate logistic regression. These models include radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model). The performance and clinical benefits of the prediction model are assessed by the area under a receiver operating characteristic (ROC) curve and the decision curve, respectively
Statistical analysis
Statistical analysis was conducted with R software (Version: 3.6.4, http: www.r-project.org/). The reported significance levels were all two-sided, and the statistical significance level was set to 0.05. The multivariate logistic regression analysis was performed with the “stats” package. Nomogram construction was performed using the “rms” package. Decision curve analysis was performed using the “dca. R” package.
Results
Imaging and clinical manifestations between groups
The differences between patients with and without COVID-19 for all 67 features (41 imaging and 26 critical clinical features) are shown in Tables 1 and 2 and the Supplementary Materials (Tables E3 and E4). The differences between the primary cohort and validation cohort for the same features are shown in the Supplementary Materials (Tables E5 and E6). All characteristics except fatigue and white blood cell count in the CR model presented no significant difference between the primary and validation cohorts. A total of 1745 lesions were identified, with 1062 from the COVID-19 group and 683 from the non-COVID-19 group.
Table 1 Radiological semantic features of patients in COVID-19 and non-COVID-19
Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value
Number of pure GGO
Total# 1.00 (0.00, 5.05) 3.50 (0.95, 8.05) 0.018b*
Peripheral area# 1.00 (0.00, 4.05) 2.00 (0.00, 6.05) 0.032b*
Central/both peripheral and central area# 0.00 (0.00, 0.00) 0.00 (0.00, 2.00) 0.001b*
Number of mixed GGO
Total# 1.00 (0.00, 3.05) 3.00 (1.00, 9.00) 0.001b*
Peripheral area# 0.00 (0.00, 2.00) 2.50 (1.00, 6.00) < 0.001b*
Central/both peripheral and central area# 0.00 (0.00, 1.05) 0.00 (0.00, 2.00) 0.657b
Total number of consolidation
Consolidation# 1.00 (0.00, 3.00) 0.00 (0.00, 0.05) 0.001b*
Pure solid nodules# 0.00 (0.00, 0.00) 0.00 (0.00, 0.00) 0.309b
Solid nodules with GGO# 0.00 (0.00, 0.00) 0.00 (0.00, 1.00) 0.033b*
Total number of lesions
Peripheral area# 5.00 (2.00, 9.05) 7.00 (2.00, 13.00) 0.112b
Central area# 0.00 (0.00, 3.00) 0.00 (0.00, 1.05) 0.960b
Both peripheral and central area# 0.00 (0.00, 2.00) 0.00 (0.00, 2.05) 0.582b
Interlobular septal thickening 0.009a*
Negative 44 (66.67%) 31 (44.29%)
Positive 22 (33.33%) 39 (55.71%)
Crazy paving pattern < 0.001a*
Negative 60 (90.91%) 32 (45.71%)
Positive 6 (9.09%) 38 (54.29%)
Tree-in-bud sign < 0.001a*
Negative 37 (56.06%) 61 (87.14%)
Positive 29 (43.94%) 9 (12.86%)
Pleural thickening 0.030a*
Negative 46 (69.70%) 36 (51.43%)
Positive 20 (30.30%) 34 (48.57%)
Offending vessel augmentation in lesions < 0.001a*
Negative 55 (83.33%) 17 (24.29%)
Positive 11 (16.67%) 53 (75.71%)
GGO ground-glass opacities
#Results are median with interquartile range in parentheses, and the remainder results are measurements with corresponding ratio in parentheses
*Data with statistical significance. pa: chi-square test, pb: Student’s t test
Table 2 Clinical features of patients in COVID-19 and non-COVID-19
Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value
Sex
Male# 43 (65.15%) 41 (58.57%) 0.430a
Female# 23 (34.85%) 29 (41.43%)
Age (years) 46.73 ± 25.00 42.93 ± 13.32 0.275b
Vital signs
Systolic blood pressure (mmHg) 126.92 ± 23.07 127.07 ± 15.16 0.965b
Diastolic blood pressure (mmHg) 77.74 ± 15.72 80.39 ± 10.51 0.254b
Respiration rate (bpm) 25.20 ± 7.29 19.86 ± 1.90 < 0.001b*
Heart rate (bpm) 101.59 ± 20.36 86.06 ± 13.34 < 0.001b*
Temperature (°C) 37.61 ± 1.06 37.12 ± 0.83 0.003b*
Signs
Dry cough# 56 (84.85%) 48 (68.57%) 0.025a*
Fatigue# 8 (12.12%) 22 (31.43%) 0.007a*
Sore throat# 6 (9.09%) 9 (12.86%) 0.483a
Stuffy# 4 (6.06%) 2 (2.86%) 0.623a
Runny nose# 3 (4.55%) 3 (4.29%) 0.731a
White blood cell count (× 109/L) 11.48 ± 5.36 5.27 ± 2.33 < 0.001b*
White blood cell count category < 0.001c*
Low# 0 (0.00%) 2 (2.86%)
Normal# 27 (40.91%) 63 (90.00%)
High# 39 (59.09%) 5 (7.14%)
Lymphocyte count (× 109/L) 1.57 ± 1.33 1.25 ± 0.68 0.086b
Lymphocyte count category < 0.001c*
Low# 24 (36.36%) 32 (45.71%)
Normal# 35 (53.03%) 37 (52.86%)
High# 7 (10.61%) 1 (1.43%)
Neutrophil count (× 109/L) 8.97 ± 4.90 3.53 ± 2.17 < 0.001b*
Neutrophil count category < 0.001c*
Low# 3 (4.55%) 8 (11.43%)
Normal# 23 (34.85%) 59 (84.29%)
High# 40 (60.61%) 3 (4.29%)
C-reactive protein (mg/L) 69.30 ± 65.88 26.37 ± 30.97 < 0.001b*
Procalcitonin (ng/mL) 3.36 ± 8.98 0.26 ± 0.84 0.007b*
*Data with statistical significance. pa: chi-square test, pb: Student’s t test. pc: Kruskal-Wallis H test
#Results are measurements with corresponding ratio in parentheses
For imaging manifestations, 7 patients in the COVID-19 group showed normal chest CT (10%). COVID-19 patients have a greater number of pure GGO and mixed GGO than non-COVID-19 patients (p = 0.018 and p = 0.001, respectively). For pure GGO lesions, the differences are significant both in peripheral (p = 0.032) and in central areas (p = 0.001). However, the number of mixed GGO is mainly distributed at the periphery in COVID-19 patients (p < 0.001), with no statistical difference in the central area. The consolidation lesions without GGO occurred less in COVID-19 patients (p = 0.001). More lesions are between 1 and 3 cm (p = 0.027), and fewer lesions are larger than half of the lung segment (p = 0.017) in COVID-19 patients. Other significant differences between the two groups include the pleural traction sign (p = 0.019), bronchial wall thickening (p < 0.001), interlobular septal thickening (p = 0.009), crazy paving (p < 0.001), tree-in-bud (p < 0.001), pleural effusions (p < 0.001), pleural thickening (p = 0.030), and the offending vessel augmentation in lesions (p < 0.001). The lung score presents no significant difference between the COVID-19 and non-COVID-19 groups.
Comparison of clinical features between the two groups of patients with and without COVID-19 is reported in Table 2. There is no significant difference in age and sex between the two groups. Significant differences are found in common symptoms between groups, including fever (p = 0.003), dry cough (p = 0.025), and fatigue (p = 0.007). The respiration rate and heart rate also show significant differences between the two groups (both p < 0.001). Compared with non-COVID-19 pneumonia, the reduction of the WBC count is more pronounced in COVID-19 patients (p < 0.001). The ratio of lymphocyte and ratio of neutrophil also show a significant difference between COVID-19 and non-COVID-19 groups. Although lymphopenia was observed in 32 COVID-19 patients (45.71%), it is not statistically different compared with that in the non-COVID-19 group. C-creative protein (CRP) level and procalcitonin level are also significantly different between the two groups (p < 0.001 and p = 0.007, respectively). Most COVID-19 patients present normal procalcitonin level (82.86%).
Clinical and radiological feature selection
Of the features, 18 radiological features and 17 clinical features were selected to form the predictors based on the result from Tables 1 and 2. Table 3 lists the features selected by univariate analysis and LASSO.
Table 3 Selected features in C, R, and CR models
Model and individual features Coefficients
R, n = 8 (41)*
Intercept − 0.307
Total number of mixed GGO in peripheral area 0.359
Total number of consolidation − 1.262
Total number of solid nodules with ground-glass opacities 0.452
Interlobular septal thickening − 5.559
Crazy paving pattern 3.566
Tree-in-bud − 2.548
Pleural thickening 3.265
Offending vessel augmentation in lesions 5.504
C, n = 7 (26)*
Intercept 29.273
Respiration − 0.359
Heart rate − 0.054
Temperature − 0.289
White blood cell count − 0.175
Cough − 1.866
Fatigue 2.855
Lymphocyte count category − 0.028
CR, n = 10 (67)*
Intercept 45.117
Total number of mixed GGO in peripheral area 0.108
Tree-in-bud − 1.853
Offending vessel augmentation in lesions 6.000
Respiration − 0.583
Heart ratio − 0.084
Temperature − 0.536
White blood cell count − 0.471
Cough − 0.997
Fatigue − 0.228
Lymphocyte count category − 2.177
C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively
*n means corresponding selected features, and data in parentheses are total features. Coefficients: the estimate value of each feature in multivariate logistic regression model by “glm” package in R
Model development and validation
The prediction models based on (i) clinical features (C model), (ii) radiological features (R model), and (iii) the combination of clinical features and radiological features (CR model) were developed. ROC analyses for the primary and validation cohort are shown in Table 4 and Fig. 3. The CR model yielded a maximum AUC of 0.986 (95% CI 0.966~1.000) in the primary cohort with the highest accuracy and specificity, which was 0.936 (95% CI 0.866~1.000) in the validation cohort. The AUC for the C model was 0.952 (95% CI 0.988~0.915) and 0.967 (95% CI 0.919~1.000) in the primary and validation cohorts, respectively. For the R model, the AUC of the two cohorts was 0.969 (95% CI 0.940~0.997) and 0.809 (95% CI 0.669~0.948), respectively.
Table 4 Performance of the individualized prediction models
Primary cohort (n = 98) Validation cohort (n = 38)
Models AUC 95% CI Accuracy Specificity Sensitivity AUC 95% CI Accuracy Specificity Sensitivity
C model 0.952 0.915~0.988 0.888 0.894 0.882 0.967 0.919~1.000 0.868 0.859 0.842
R model 0.969 0.940~0.997 0.929 0.851 1.000 0.809 0.669~0.948 0.684 0.368 1.000
CR model 0.986 0.966~1.000 0.959 0.957 0.961 0.936 0.866~1.000 0.763 0.789 0.737
C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively. CI confidence interval
Fig. 3 ROC of the three models in primary and validation cohort curves. Comparison of receiver operating characteristic (ROC) curves among the radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model) for the diagnosis of COVID-19 in the primary (a) and validation (b) cohorts
To determine the clinical usefulness of the diagnostic model, we developed the decision curve (Fig. 4), which showed better performances for the CR model compared with that for the C model and the R model. Across the majority of the range of reasonable threshold probabilities, the decision curve analysis showed that the CR model had a higher overall benefit than the C model and R model.
Fig. 4 Decision curve analysis for each model in the primary dataset. The y-axis measures the net benefit, which is calculated by summing the benefits (true-positive findings) and subtracting the harms (false-positive findings), weighting the latter by a factor related to the relative harm of undetected metastasis compared with the harm of unnecessary treatment. The decision curve shows that if the threshold probability is over 10%, the application of the combination of clinical and radiological model (CR model) to diagnose COVID-19 adds more benefit than the clinical model (C model) and radiological model (R model)
The nomogram (Fig. 5) was developed by the CR model in the primary cohort, with the factors of the total number of mixed GGO in peripheral area (TN_Mixed_GGO_IP), tree-in-bud, offending vessel augmentation in lesions (OVAIL), respiration, heart ratio, temperature, white blood cell count, cough, fatigue and lymphocyte count category incorporated. The total points were calculated by summing the points identified on the “points” scale for each factor. By comparing the “total points” scale and the “probability” scale, the individual probability of COVID-19 infection could be obtained.
Fig. 5 Nomogram of the CR model in the primary cohort. TN_Mixed_GGO_IP represented the total number of mixed GGO in peripheral area. AVAIL represented offending vessel segmentation in lesions. N was a negative result, and P was a positive result. Norm represented normal. Note that in probability scale, 0 = non-COVID-19, 1 = COVID-19
Discussion
In this multi-center study, statistical analysis was performed in comparing imaging and clinical manifestations between pneumonia patients with and without COVID-19. Eighteen radiological semantic features and seventeen clinical features were identified to be significantly different between the two groups (p < 0.05). Three models for COVID-19 diagnosis were developed based on the refined features. The models were validated in the both primary and validation cohorts and achieved an AUC as high as 0.986. These models will play an essential role for early and easy-to-access diagnosis, especially when there are not enough RT-PCT kits or experimental platforms to test for the COVID-19 infection.
A total of 1745 lesions were evaluated for the qualitative feature, location, and size in this study. Consistent with the previous studies, the ground-glass opacities and consolidation in the lung periphery were considered to be the imaging hallmark in patients with COVID-19 infection [11, 25]. However, when we subdivided the GGO into pure GGO and mixed GGO, we found that the distribution pattern is different between these two lesions. Pure GGO show differences between groups in every location of the lungs, whereas mixed GGO only have significant differences between groups in the lung periphery. Recent studies defined four stages of lung involvement in COVID-19 [26]. Therefore, a follow-up analysis of these distributions would be significant. The lesion size in patients with COVID-19 infection was another interesting observation. Most lesions were between 1 and 3 cm, with few lesions larger than half of the lung segment, which was similar to the finding in MERS_CoV [22]. Other features similar to MERS_CoV and SARS_CoV were observed in the laboratory abnormalities, such as lymphopenia, which may be associated with the cellular immune deficiency [3, 27]. However, our results showed no significant difference in lymphopenia between the COVID-19 and non-COVID-19 patients.
To our knowledge, no diagnostic model based on imaging and clinical features alone has been proposed for the diagnosis of COVID-19. Our clinical and radiological semantic (CR) models consisted of the following features: total number of GGO with consolidation in the peripheral area, tree-in-bud, offending vessel augmentation in lesions, temperature, heart ratio, respiration, cough and fatigue, WBC count, and lymphocyte count category. The CR model outperformed the individual clinical and radiologic model. This result was in accordance with that in previous study in breast cancer, in which the model based on the combination of radiomics features and clinical features achieved a higher performance [24]. Compared with the radiomics-based model, the extraction of radiological semantic features can overcome the image discrepancy caused by different scanning parameters and/or different CT vendors. A previous study [28] also indicated that models based on semantic features determined by an experienced thoracic radiologist slightly outperformed models based on computed texture features alone.
There are a few limitations in this study. First, the sample size is relatively small because this is a retrospective analysis of a new disease and most of the cases outside of Wuhan City are imported. Second, with the multi-center retrospective design, there is a potential bias of patient selection [29], since there may be some deviations in marking semantic features among readers, though we have taken the effort to reduce this by creating pictorial examples and setting feature criteria (Supplementary Materials). Third, longitudinal CT study was not performed. Whether or not this model can be used to evaluate the follow-ups and help to guide therapy remains an open question to be further explored. Moreover, the rich high-order features of the CT image combined with radiomics or deep learning have not been studied, which may be another way to identify the patients with COVID-19. Besides, one can also focus on the role of radiological features in disease monitoring, treatment evaluation, and prognosis prediction.
In conclusion, 1745 lesions and 67 features were compared between pneumonia patients with and without COVID-19. Thirty-five features were significantly different between the two groups. A diagnostic model with AUC as high as 0.986 was developed and validated both in the primary and in the validation cohorts, which may help improve the COVID-19 diagnosis.
Electronic supplementary material
ESM 1 (DOCX 404 kb)
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Xiaofeng Chen and Yanyan Tang contributed equally to this work.
Funding information
This study has received funding by the Natural Science Foundation of China (grant numbers 81471730, 31870981) to R.W.; the Natural Science Foundation of Guangdong Province (grant number 2018A030307057) to Z.D.; and the Special Project on Prevention and Control of COVID-19 for Colleges and Universities in Guangdong Province (grant number 2020KZDZX1085) to Z.D.
Compliance with ethical standards
Guarantor
The scientific guarantor of this publication is Zhuozhi Dai.
Conflict of interest
One of the authors of this manuscript (Yuting Liao) is an employee of GE Healthcare. The remaining authors declare no relationships with any companies whose products or services may be related to the subject matter of the article.
Statistics and biometry
One of the authors, Dr. Yuting Liao, has significant statistical expertise.
Informed consent
Written informed consent was waived by the Institutional Review Board.
Ethical approval
Institutional Review Board approval was obtained.
Methodology
• retrospective
• case-control study
• multi-center study
|
Document structure show
article-title | A diagnostic model for coronavirus disease 2019 (COVID-19) based on radiological semantic and clinical features: a multi-center study |
abstract | Objectives Rapid and accurate diagnosis of coronavirus disease 2019 (COVID-19) is critical during the epidemic. We aim to identify differences in CT imaging and clinical manifestations between pneumonia patients with and without COVID-19, and to develop and validate a diagnostic model for COVID-19 based on radiological semantic and clinical features alone. Methods A consecutive cohort of 70 COVID-19 and 66 non-COVID-19 pneumonia patients were retrospectively recruited from five institutions. Patients were divided into primary (n = 98) and validation (n = 38) cohorts. The chi-square test, Student’s t test, and Kruskal-Wallis H test were performed, comparing 1745 lesions and 67 features in the two groups. Three models were constructed using radiological semantic and clinical features through multivariate logistic regression. Diagnostic efficacies of developed models were quantified by receiver operating characteristic curve. Clinical usage was evaluated by decision curve analysis and nomogram. Results Eighteen radiological semantic features and seventeen clinical features were identified to be significantly different. Besides ground-glass opacities (p = 0.032) and consolidation (p = 0.001) in the lung periphery, the lesion size (1–3 cm) is also significant for the diagnosis of COVID-19 (p = 0.027). Lung score presents no significant difference (p = 0.417). Three diagnostic models achieved an area under the curve value as high as 0.986 (95% CI 0.966~1.000). The clinical and radiological semantic models provided a better diagnostic performance and more considerable net benefits. Conclusions Based on CT imaging and clinical manifestations alone, the pneumonia patients with and without COVID-19 can be distinguished. A model composed of radiological semantic and clinical features has an excellent performance for the diagnosis of COVID-19. Key Points • Based on CT imaging and clinical manifestations alone, the pneumonia patients with and without COVID-19 can be distinguished. • A diagnostic model for COVID-19 was developed and validated using radiological semantic and clinical features, which had an area under the curve value of 0.986 (95% CI 0.966~1.000) and 0.936 (95% CI 0.866~1.000) in the primary and validation cohorts, respectively. Electronic supplementary material The online version of this article (10.1007/s00330-020-06829-2) contains supplementary material, which is available to authorized users. |
sec | Objectives Rapid and accurate diagnosis of coronavirus disease 2019 (COVID-19) is critical during the epidemic. We aim to identify differences in CT imaging and clinical manifestations between pneumonia patients with and without COVID-19, and to develop and validate a diagnostic model for COVID-19 based on radiological semantic and clinical features alone. |
title | Objectives |
p | Rapid and accurate diagnosis of coronavirus disease 2019 (COVID-19) is critical during the epidemic. We aim to identify differences in CT imaging and clinical manifestations between pneumonia patients with and without COVID-19, and to develop and validate a diagnostic model for COVID-19 based on radiological semantic and clinical features alone. |
sec | Methods A consecutive cohort of 70 COVID-19 and 66 non-COVID-19 pneumonia patients were retrospectively recruited from five institutions. Patients were divided into primary (n = 98) and validation (n = 38) cohorts. The chi-square test, Student’s t test, and Kruskal-Wallis H test were performed, comparing 1745 lesions and 67 features in the two groups. Three models were constructed using radiological semantic and clinical features through multivariate logistic regression. Diagnostic efficacies of developed models were quantified by receiver operating characteristic curve. Clinical usage was evaluated by decision curve analysis and nomogram. |
title | Methods |
p | A consecutive cohort of 70 COVID-19 and 66 non-COVID-19 pneumonia patients were retrospectively recruited from five institutions. Patients were divided into primary (n = 98) and validation (n = 38) cohorts. The chi-square test, Student’s t test, and Kruskal-Wallis H test were performed, comparing 1745 lesions and 67 features in the two groups. Three models were constructed using radiological semantic and clinical features through multivariate logistic regression. Diagnostic efficacies of developed models were quantified by receiver operating characteristic curve. Clinical usage was evaluated by decision curve analysis and nomogram. |
sec | Results Eighteen radiological semantic features and seventeen clinical features were identified to be significantly different. Besides ground-glass opacities (p = 0.032) and consolidation (p = 0.001) in the lung periphery, the lesion size (1–3 cm) is also significant for the diagnosis of COVID-19 (p = 0.027). Lung score presents no significant difference (p = 0.417). Three diagnostic models achieved an area under the curve value as high as 0.986 (95% CI 0.966~1.000). The clinical and radiological semantic models provided a better diagnostic performance and more considerable net benefits. |
title | Results |
p | Eighteen radiological semantic features and seventeen clinical features were identified to be significantly different. Besides ground-glass opacities (p = 0.032) and consolidation (p = 0.001) in the lung periphery, the lesion size (1–3 cm) is also significant for the diagnosis of COVID-19 (p = 0.027). Lung score presents no significant difference (p = 0.417). Three diagnostic models achieved an area under the curve value as high as 0.986 (95% CI 0.966~1.000). The clinical and radiological semantic models provided a better diagnostic performance and more considerable net benefits. |
sec | Conclusions Based on CT imaging and clinical manifestations alone, the pneumonia patients with and without COVID-19 can be distinguished. A model composed of radiological semantic and clinical features has an excellent performance for the diagnosis of COVID-19. |
title | Conclusions |
p | Based on CT imaging and clinical manifestations alone, the pneumonia patients with and without COVID-19 can be distinguished. A model composed of radiological semantic and clinical features has an excellent performance for the diagnosis of COVID-19. |
sec | Key Points • Based on CT imaging and clinical manifestations alone, the pneumonia patients with and without COVID-19 can be distinguished. • A diagnostic model for COVID-19 was developed and validated using radiological semantic and clinical features, which had an area under the curve value of 0.986 (95% CI 0.966~1.000) and 0.936 (95% CI 0.866~1.000) in the primary and validation cohorts, respectively. |
title | Key Points |
p | • Based on CT imaging and clinical manifestations alone, the pneumonia patients with and without COVID-19 can be distinguished. |
p | • A diagnostic model for COVID-19 was developed and validated using radiological semantic and clinical features, which had an area under the curve value of 0.986 (95% CI 0.966~1.000) and 0.936 (95% CI 0.866~1.000) in the primary and validation cohorts, respectively. |
sec | Electronic supplementary material The online version of this article (10.1007/s00330-020-06829-2) contains supplementary material, which is available to authorized users. |
title | Electronic supplementary material |
p | The online version of this article (10.1007/s00330-020-06829-2) contains supplementary material, which is available to authorized users. |
body | Introduction On January 30, 2020, the World Health Organization (WHO) has declared the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreak as a global health emergency of international concern. This outbreak has infected all provinces of China and rapidly spread to the rest of the world. At the time of writing this article (March 16, 2020), there have been more than 158 countries and territories affected [1]. Whole-genome sequencing and phylogenetic analysis reveal that the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is similar to some beta coronaviruses detected in bats, but it is distinct from severe acute respiratory syndrome coronavirus (SARS-Cov) and Middle East respiratory syndrome coronavirus (MERS-CoV) [2]. Patients with COVID-19 develop pneumonia with associated symptoms of fever (98%), cough (76%), and myalgia or fatigue (44%) [3]. CT imaging plays a critical role in the diagnosis and the monitoring of disease progression [4–6]. The latest research studies described the characteristic imaging manifestations of COVID-19, including ground-glass opacities (GGO) (57 to 88%), bilateral involvement (76 to 88%), and peripheral distribution (33 to 85%) [7–10]. Other imaging features such as consolidation, cavitation, and interlobular septal thickening are also reported in some patients [11–13]. However, these imaging manifestations of COVID-19 are nonspecific and are difficult to distinguish from other pneumonia. To our knowledge, there have been no studies explicitly comparing imaging and clinical characteristics between pneumonia patients with and without COVID-19. The current diagnostic criterion for COVID-19 is the positive result of a nucleic acid test by real-time reverse transcription polymerase chain reaction (RT-PCR) or next-generation sequencing [14]. However, false-negative results caused by unstable specimen processing are relatively high in clinical practice, which has worsened the spread of the outbreak [15–18]. Moreover, laboratory testing for SARS-CoV-2 requires a rigorous platform, which is not assembled in all hospitals. Thus, this requires specimen transfer, which may delay diagnosis for days. Early and accurate diagnosis is crucial, particularly for critically ill patients who need emergency surgery, and with pneumonia complications. To solve these problems, we hypothesize that a diagnostic model can be developed based on CT imaging and clinical manifestations alone, independent of the nucleic acid test. In this study, we identify the differences in imaging and clinical manifestations between patients with and without COVID-19. We also develop and validate a model for COVID-19 diagnosis based on radiological semantic and clinical features. Patients and methods Patients Ethical approvals by the institutional review boards were obtained for this retrospective analysis, and the need to obtain informed consent was waived. From January 1 to February 8, 2020, seventy consecutive patients with COVID-19 admitted in 5 independent hospitals from 4 cities were enrolled in this study (mean age, 42.9 years; range, 16–69 years), including 41 men (mean age, 41.8 years; range, 16–69 years) and 29 women (mean age, 44.5 years; range, 16–66 years). All patients were confirmed with SARS-CoV-2 infection by real-time RT-PCR and next-generation sequencing. Of these patients, 24 were from Huizhou City, 25 from Shantou City, 15 from Yongzhou City, and the rest 6 from Meizhou City. At the same period, another 66 pneumonia patients without COVID-19 from Meizhou People’s Hospital were recruited as controls (mean age, 46.7 years; range, 0.3–93 years), including 43 men (mean age, 46.0 years; range, 0.3–93 years) and 23 women (mean age, 48.0 years; range, 1–86 years). All the controls were confirmed with consecutive negative RT-PCR assays. Figure E1 in the Supplementary Material shows the patient recruitment pathway for the control group, along with the inclusion and exclusion criteria. According to previous studies [19–21], whose sample size is comparable with ours, the ratio between primary and validation cohort is 7:3. In this study, a total of 136 patients were divided into primary (n = 98) and validation (n = 38) cohorts, close to 7:3. A total of 19 COVID-19 patients from two hospitals (6 patients from Meizhou People’s Hospital and 13 patients from the First Affiliated Hospital of Shantou University Medical College) and 19 randomly selected controls from Meizhou City were incorporated into the validation cohort. The rest of the patients are incorporated in the primary cohort, including 51 COVID-19 patients from Huizhou, Yongzhou, and Shantou cities and 47 controls from Meizhou City. The primary cohort was utilized to select the most valuable features and build the predictive model, and the validation cohort was used to evaluate and validate the performance of the model. Image and clinical data collection The chest CT imaging data without contrast material enhancement were obtained from multiple hospitals with different CT systems, including GE CT Discovery 750 HD (General Electric Company), SCENARIA 64 CT (Hitachi Medical), Philips Ingenuity CT (PHILIPS), and Siemens SOMATOM Definition AS (Siemens). All images were reconstructed into 1-mm slices with a slice gap of 0.8 mm. Detailed acquisition parameters were summarized in the Supplementary Material (Table E1). The clinical history, nursing records, and laboratory findings were reviewed for all patients. Clinical characteristics, including demographic information, daily body temperature, blood pressure, heart rate, clinical symptoms, and history of exposure to epidemic centers, were collected. Total white blood cell (WBC) counts, lymphocyte counts, ratio of lymphocyte, neutrophil count, ratio of neutrophil, procalcitonin (PCT), C-reactive protein level (CRP), and erythrocyte sedimentation rate (ESR) were measured. All threshold values chosen for laboratory metrics were based on the normal ranges set by each individual hospital. Image analysis For extraction of radiological semantic features, two senior radiologists (D.L. and X.C., more than 15 years of experience) reached a consensus, blinded to clinical and laboratory findings. The radiological semantic features included both qualitative and quantitative imaging features. The lesions in the outer third of the lung were defined as peripheral, and lesions in the inner two-thirds of the lung were defined as central [22]. The progression of COVID-19 lesions within each lung lobe was evaluated by scoring each lobe from 0 to 4 [7], corresponding to normal, 1~25% infection, 26~50% infection, 51~75% infection, and more than 75% infection, respectively. The scores were combined for all five lobes to provide a total score ranging from 0 to 20. A total of 41 radiological features (26 quantitative and 15 qualitative) were extracted for the analysis. The descriptions of radiological semantic features are listed in the Supplementary Material (Table E2). Figure 1 is one example of the evaluation of CT imaging. Fig. 1 A 23-year-old female with a travel history to Wuhan presenting with fever. Axial noncontrast CT image shows a consolidation with ground-glass opacities in the peripheral region by the right upper lobe. Air bronchogram is found in lesion. The maximum diameter of lesion is 2.8 cm. The right upper lobe score is 1 because of the involved lung parenchyma less than 1/4 Clinical and radiological feature selection To obtain the most valuable clinical and radiological semantic features, statistical analysis, univariate analysis, and the least absolute shrinkage and selection operator (LASSO) method were performed. In statistical analysis, the chi-square test, the Kruskal-Wallis H test, and t test were utilized to compare the radiological semantic and clinical features between COVID-19 and non-COVID-19 groups. The features with p value smaller than 0.05 were selected. Then, univariate analysis was performed for clinical and radiological candidate features to determine the COVID-19 risk factors. The features with p value smaller than 0.05 in univariate analysis were also selected. The least absolute shrinkage and selection operator (LASSO) method [23] was utilized to select the most useful features with penalty parameter tuning that was conducted by 10-fold cross-validation based on minimum criteria. Diagnostic models were then constructed by multivariate logistic regression with the selected features. The flowchart of the feature selection process for these models was presented in the Supplementary Material (Fig. E2). Development and validation of the diagnostic model To develop an optimal model, we evaluated 3 models by analyzing (i) the clinical features model (C model), (ii) radiological semantic features model (R model), and (iii) the combination of clinical and radiological semantic features model (CR model) by multivariate logistic regression analysis. The classification performances of the models were evaluated by the area under the receiver operating characteristic (ROC) curve. The area under the curve (AUC), accuracy, sensitivity, and specificity were also calculated. A decision curve analysis was conducted to determine the clinical usefulness of the diagnostic model by quantifying the net benefits at different threshold probabilities in the validation dataset [24]. The development of decision curve was described in the Supplementary Materials. Figure 2 depicts the flowchart of the proposed analysis pipeline described above. We also built a nomogram, which was a quantitative tool to predict the individual probability of infection by COVID-19, based on the multivariate logistic analysis of the CR model with the primary cohort. Depending on the coefficient of the predictive factors in multivariate logistic regression model, all values of each predictive factor were assigned points. A total point was obtained by summing all the points of each predictive factor. The scale also showed the relationship between the total point and the prediction probability in the nomogram. The corresponding calibration curves of the CR model in the primary cohort and validation cohort are shown in the Supplementary Material (Fig. E3). Fig. 2 Workflow of data process and analysis in this study. Radiological semantic features, including qualitative and quantitative imaging features, are extracted from axial lung CT section. The clinical manifestation and laboratory parameters are provided by electronic case system. Statistical analysis is performed for comparing the different features between COVID-19 and non-COVID-19 patients. Univariate analysis, least absolute shrinkage, and selection operator (LASSO) are further performed to determine the COVID-19 risk factors with p < 0.05 in statistical analysis. Three models based on the selected features are established by multivariate logistic regression. These models include radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model). The performance and clinical benefits of the prediction model are assessed by the area under a receiver operating characteristic (ROC) curve and the decision curve, respectively Statistical analysis Statistical analysis was conducted with R software (Version: 3.6.4, http: www.r-project.org/). The reported significance levels were all two-sided, and the statistical significance level was set to 0.05. The multivariate logistic regression analysis was performed with the “stats” package. Nomogram construction was performed using the “rms” package. Decision curve analysis was performed using the “dca. R” package. Results Imaging and clinical manifestations between groups The differences between patients with and without COVID-19 for all 67 features (41 imaging and 26 critical clinical features) are shown in Tables 1 and 2 and the Supplementary Materials (Tables E3 and E4). The differences between the primary cohort and validation cohort for the same features are shown in the Supplementary Materials (Tables E5 and E6). All characteristics except fatigue and white blood cell count in the CR model presented no significant difference between the primary and validation cohorts. A total of 1745 lesions were identified, with 1062 from the COVID-19 group and 683 from the non-COVID-19 group. Table 1 Radiological semantic features of patients in COVID-19 and non-COVID-19 Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value Number of pure GGO Total# 1.00 (0.00, 5.05) 3.50 (0.95, 8.05) 0.018b* Peripheral area# 1.00 (0.00, 4.05) 2.00 (0.00, 6.05) 0.032b* Central/both peripheral and central area# 0.00 (0.00, 0.00) 0.00 (0.00, 2.00) 0.001b* Number of mixed GGO Total# 1.00 (0.00, 3.05) 3.00 (1.00, 9.00) 0.001b* Peripheral area# 0.00 (0.00, 2.00) 2.50 (1.00, 6.00) < 0.001b* Central/both peripheral and central area# 0.00 (0.00, 1.05) 0.00 (0.00, 2.00) 0.657b Total number of consolidation Consolidation# 1.00 (0.00, 3.00) 0.00 (0.00, 0.05) 0.001b* Pure solid nodules# 0.00 (0.00, 0.00) 0.00 (0.00, 0.00) 0.309b Solid nodules with GGO# 0.00 (0.00, 0.00) 0.00 (0.00, 1.00) 0.033b* Total number of lesions Peripheral area# 5.00 (2.00, 9.05) 7.00 (2.00, 13.00) 0.112b Central area# 0.00 (0.00, 3.00) 0.00 (0.00, 1.05) 0.960b Both peripheral and central area# 0.00 (0.00, 2.00) 0.00 (0.00, 2.05) 0.582b Interlobular septal thickening 0.009a* Negative 44 (66.67%) 31 (44.29%) Positive 22 (33.33%) 39 (55.71%) Crazy paving pattern < 0.001a* Negative 60 (90.91%) 32 (45.71%) Positive 6 (9.09%) 38 (54.29%) Tree-in-bud sign < 0.001a* Negative 37 (56.06%) 61 (87.14%) Positive 29 (43.94%) 9 (12.86%) Pleural thickening 0.030a* Negative 46 (69.70%) 36 (51.43%) Positive 20 (30.30%) 34 (48.57%) Offending vessel augmentation in lesions < 0.001a* Negative 55 (83.33%) 17 (24.29%) Positive 11 (16.67%) 53 (75.71%) GGO ground-glass opacities #Results are median with interquartile range in parentheses, and the remainder results are measurements with corresponding ratio in parentheses *Data with statistical significance. pa: chi-square test, pb: Student’s t test Table 2 Clinical features of patients in COVID-19 and non-COVID-19 Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value Sex Male# 43 (65.15%) 41 (58.57%) 0.430a Female# 23 (34.85%) 29 (41.43%) Age (years) 46.73 ± 25.00 42.93 ± 13.32 0.275b Vital signs Systolic blood pressure (mmHg) 126.92 ± 23.07 127.07 ± 15.16 0.965b Diastolic blood pressure (mmHg) 77.74 ± 15.72 80.39 ± 10.51 0.254b Respiration rate (bpm) 25.20 ± 7.29 19.86 ± 1.90 < 0.001b* Heart rate (bpm) 101.59 ± 20.36 86.06 ± 13.34 < 0.001b* Temperature (°C) 37.61 ± 1.06 37.12 ± 0.83 0.003b* Signs Dry cough# 56 (84.85%) 48 (68.57%) 0.025a* Fatigue# 8 (12.12%) 22 (31.43%) 0.007a* Sore throat# 6 (9.09%) 9 (12.86%) 0.483a Stuffy# 4 (6.06%) 2 (2.86%) 0.623a Runny nose# 3 (4.55%) 3 (4.29%) 0.731a White blood cell count (× 109/L) 11.48 ± 5.36 5.27 ± 2.33 < 0.001b* White blood cell count category < 0.001c* Low# 0 (0.00%) 2 (2.86%) Normal# 27 (40.91%) 63 (90.00%) High# 39 (59.09%) 5 (7.14%) Lymphocyte count (× 109/L) 1.57 ± 1.33 1.25 ± 0.68 0.086b Lymphocyte count category < 0.001c* Low# 24 (36.36%) 32 (45.71%) Normal# 35 (53.03%) 37 (52.86%) High# 7 (10.61%) 1 (1.43%) Neutrophil count (× 109/L) 8.97 ± 4.90 3.53 ± 2.17 < 0.001b* Neutrophil count category < 0.001c* Low# 3 (4.55%) 8 (11.43%) Normal# 23 (34.85%) 59 (84.29%) High# 40 (60.61%) 3 (4.29%) C-reactive protein (mg/L) 69.30 ± 65.88 26.37 ± 30.97 < 0.001b* Procalcitonin (ng/mL) 3.36 ± 8.98 0.26 ± 0.84 0.007b* *Data with statistical significance. pa: chi-square test, pb: Student’s t test. pc: Kruskal-Wallis H test #Results are measurements with corresponding ratio in parentheses For imaging manifestations, 7 patients in the COVID-19 group showed normal chest CT (10%). COVID-19 patients have a greater number of pure GGO and mixed GGO than non-COVID-19 patients (p = 0.018 and p = 0.001, respectively). For pure GGO lesions, the differences are significant both in peripheral (p = 0.032) and in central areas (p = 0.001). However, the number of mixed GGO is mainly distributed at the periphery in COVID-19 patients (p < 0.001), with no statistical difference in the central area. The consolidation lesions without GGO occurred less in COVID-19 patients (p = 0.001). More lesions are between 1 and 3 cm (p = 0.027), and fewer lesions are larger than half of the lung segment (p = 0.017) in COVID-19 patients. Other significant differences between the two groups include the pleural traction sign (p = 0.019), bronchial wall thickening (p < 0.001), interlobular septal thickening (p = 0.009), crazy paving (p < 0.001), tree-in-bud (p < 0.001), pleural effusions (p < 0.001), pleural thickening (p = 0.030), and the offending vessel augmentation in lesions (p < 0.001). The lung score presents no significant difference between the COVID-19 and non-COVID-19 groups. Comparison of clinical features between the two groups of patients with and without COVID-19 is reported in Table 2. There is no significant difference in age and sex between the two groups. Significant differences are found in common symptoms between groups, including fever (p = 0.003), dry cough (p = 0.025), and fatigue (p = 0.007). The respiration rate and heart rate also show significant differences between the two groups (both p < 0.001). Compared with non-COVID-19 pneumonia, the reduction of the WBC count is more pronounced in COVID-19 patients (p < 0.001). The ratio of lymphocyte and ratio of neutrophil also show a significant difference between COVID-19 and non-COVID-19 groups. Although lymphopenia was observed in 32 COVID-19 patients (45.71%), it is not statistically different compared with that in the non-COVID-19 group. C-creative protein (CRP) level and procalcitonin level are also significantly different between the two groups (p < 0.001 and p = 0.007, respectively). Most COVID-19 patients present normal procalcitonin level (82.86%). Clinical and radiological feature selection Of the features, 18 radiological features and 17 clinical features were selected to form the predictors based on the result from Tables 1 and 2. Table 3 lists the features selected by univariate analysis and LASSO. Table 3 Selected features in C, R, and CR models Model and individual features Coefficients R, n = 8 (41)* Intercept − 0.307 Total number of mixed GGO in peripheral area 0.359 Total number of consolidation − 1.262 Total number of solid nodules with ground-glass opacities 0.452 Interlobular septal thickening − 5.559 Crazy paving pattern 3.566 Tree-in-bud − 2.548 Pleural thickening 3.265 Offending vessel augmentation in lesions 5.504 C, n = 7 (26)* Intercept 29.273 Respiration − 0.359 Heart rate − 0.054 Temperature − 0.289 White blood cell count − 0.175 Cough − 1.866 Fatigue 2.855 Lymphocyte count category − 0.028 CR, n = 10 (67)* Intercept 45.117 Total number of mixed GGO in peripheral area 0.108 Tree-in-bud − 1.853 Offending vessel augmentation in lesions 6.000 Respiration − 0.583 Heart ratio − 0.084 Temperature − 0.536 White blood cell count − 0.471 Cough − 0.997 Fatigue − 0.228 Lymphocyte count category − 2.177 C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively *n means corresponding selected features, and data in parentheses are total features. Coefficients: the estimate value of each feature in multivariate logistic regression model by “glm” package in R Model development and validation The prediction models based on (i) clinical features (C model), (ii) radiological features (R model), and (iii) the combination of clinical features and radiological features (CR model) were developed. ROC analyses for the primary and validation cohort are shown in Table 4 and Fig. 3. The CR model yielded a maximum AUC of 0.986 (95% CI 0.966~1.000) in the primary cohort with the highest accuracy and specificity, which was 0.936 (95% CI 0.866~1.000) in the validation cohort. The AUC for the C model was 0.952 (95% CI 0.988~0.915) and 0.967 (95% CI 0.919~1.000) in the primary and validation cohorts, respectively. For the R model, the AUC of the two cohorts was 0.969 (95% CI 0.940~0.997) and 0.809 (95% CI 0.669~0.948), respectively. Table 4 Performance of the individualized prediction models Primary cohort (n = 98) Validation cohort (n = 38) Models AUC 95% CI Accuracy Specificity Sensitivity AUC 95% CI Accuracy Specificity Sensitivity C model 0.952 0.915~0.988 0.888 0.894 0.882 0.967 0.919~1.000 0.868 0.859 0.842 R model 0.969 0.940~0.997 0.929 0.851 1.000 0.809 0.669~0.948 0.684 0.368 1.000 CR model 0.986 0.966~1.000 0.959 0.957 0.961 0.936 0.866~1.000 0.763 0.789 0.737 C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively. CI confidence interval Fig. 3 ROC of the three models in primary and validation cohort curves. Comparison of receiver operating characteristic (ROC) curves among the radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model) for the diagnosis of COVID-19 in the primary (a) and validation (b) cohorts To determine the clinical usefulness of the diagnostic model, we developed the decision curve (Fig. 4), which showed better performances for the CR model compared with that for the C model and the R model. Across the majority of the range of reasonable threshold probabilities, the decision curve analysis showed that the CR model had a higher overall benefit than the C model and R model. Fig. 4 Decision curve analysis for each model in the primary dataset. The y-axis measures the net benefit, which is calculated by summing the benefits (true-positive findings) and subtracting the harms (false-positive findings), weighting the latter by a factor related to the relative harm of undetected metastasis compared with the harm of unnecessary treatment. The decision curve shows that if the threshold probability is over 10%, the application of the combination of clinical and radiological model (CR model) to diagnose COVID-19 adds more benefit than the clinical model (C model) and radiological model (R model) The nomogram (Fig. 5) was developed by the CR model in the primary cohort, with the factors of the total number of mixed GGO in peripheral area (TN_Mixed_GGO_IP), tree-in-bud, offending vessel augmentation in lesions (OVAIL), respiration, heart ratio, temperature, white blood cell count, cough, fatigue and lymphocyte count category incorporated. The total points were calculated by summing the points identified on the “points” scale for each factor. By comparing the “total points” scale and the “probability” scale, the individual probability of COVID-19 infection could be obtained. Fig. 5 Nomogram of the CR model in the primary cohort. TN_Mixed_GGO_IP represented the total number of mixed GGO in peripheral area. AVAIL represented offending vessel segmentation in lesions. N was a negative result, and P was a positive result. Norm represented normal. Note that in probability scale, 0 = non-COVID-19, 1 = COVID-19 Discussion In this multi-center study, statistical analysis was performed in comparing imaging and clinical manifestations between pneumonia patients with and without COVID-19. Eighteen radiological semantic features and seventeen clinical features were identified to be significantly different between the two groups (p < 0.05). Three models for COVID-19 diagnosis were developed based on the refined features. The models were validated in the both primary and validation cohorts and achieved an AUC as high as 0.986. These models will play an essential role for early and easy-to-access diagnosis, especially when there are not enough RT-PCT kits or experimental platforms to test for the COVID-19 infection. A total of 1745 lesions were evaluated for the qualitative feature, location, and size in this study. Consistent with the previous studies, the ground-glass opacities and consolidation in the lung periphery were considered to be the imaging hallmark in patients with COVID-19 infection [11, 25]. However, when we subdivided the GGO into pure GGO and mixed GGO, we found that the distribution pattern is different between these two lesions. Pure GGO show differences between groups in every location of the lungs, whereas mixed GGO only have significant differences between groups in the lung periphery. Recent studies defined four stages of lung involvement in COVID-19 [26]. Therefore, a follow-up analysis of these distributions would be significant. The lesion size in patients with COVID-19 infection was another interesting observation. Most lesions were between 1 and 3 cm, with few lesions larger than half of the lung segment, which was similar to the finding in MERS_CoV [22]. Other features similar to MERS_CoV and SARS_CoV were observed in the laboratory abnormalities, such as lymphopenia, which may be associated with the cellular immune deficiency [3, 27]. However, our results showed no significant difference in lymphopenia between the COVID-19 and non-COVID-19 patients. To our knowledge, no diagnostic model based on imaging and clinical features alone has been proposed for the diagnosis of COVID-19. Our clinical and radiological semantic (CR) models consisted of the following features: total number of GGO with consolidation in the peripheral area, tree-in-bud, offending vessel augmentation in lesions, temperature, heart ratio, respiration, cough and fatigue, WBC count, and lymphocyte count category. The CR model outperformed the individual clinical and radiologic model. This result was in accordance with that in previous study in breast cancer, in which the model based on the combination of radiomics features and clinical features achieved a higher performance [24]. Compared with the radiomics-based model, the extraction of radiological semantic features can overcome the image discrepancy caused by different scanning parameters and/or different CT vendors. A previous study [28] also indicated that models based on semantic features determined by an experienced thoracic radiologist slightly outperformed models based on computed texture features alone. There are a few limitations in this study. First, the sample size is relatively small because this is a retrospective analysis of a new disease and most of the cases outside of Wuhan City are imported. Second, with the multi-center retrospective design, there is a potential bias of patient selection [29], since there may be some deviations in marking semantic features among readers, though we have taken the effort to reduce this by creating pictorial examples and setting feature criteria (Supplementary Materials). Third, longitudinal CT study was not performed. Whether or not this model can be used to evaluate the follow-ups and help to guide therapy remains an open question to be further explored. Moreover, the rich high-order features of the CT image combined with radiomics or deep learning have not been studied, which may be another way to identify the patients with COVID-19. Besides, one can also focus on the role of radiological features in disease monitoring, treatment evaluation, and prognosis prediction. In conclusion, 1745 lesions and 67 features were compared between pneumonia patients with and without COVID-19. Thirty-five features were significantly different between the two groups. A diagnostic model with AUC as high as 0.986 was developed and validated both in the primary and in the validation cohorts, which may help improve the COVID-19 diagnosis. Electronic supplementary material ESM 1 (DOCX 404 kb) |
sec | Introduction On January 30, 2020, the World Health Organization (WHO) has declared the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreak as a global health emergency of international concern. This outbreak has infected all provinces of China and rapidly spread to the rest of the world. At the time of writing this article (March 16, 2020), there have been more than 158 countries and territories affected [1]. Whole-genome sequencing and phylogenetic analysis reveal that the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is similar to some beta coronaviruses detected in bats, but it is distinct from severe acute respiratory syndrome coronavirus (SARS-Cov) and Middle East respiratory syndrome coronavirus (MERS-CoV) [2]. Patients with COVID-19 develop pneumonia with associated symptoms of fever (98%), cough (76%), and myalgia or fatigue (44%) [3]. CT imaging plays a critical role in the diagnosis and the monitoring of disease progression [4–6]. The latest research studies described the characteristic imaging manifestations of COVID-19, including ground-glass opacities (GGO) (57 to 88%), bilateral involvement (76 to 88%), and peripheral distribution (33 to 85%) [7–10]. Other imaging features such as consolidation, cavitation, and interlobular septal thickening are also reported in some patients [11–13]. However, these imaging manifestations of COVID-19 are nonspecific and are difficult to distinguish from other pneumonia. To our knowledge, there have been no studies explicitly comparing imaging and clinical characteristics between pneumonia patients with and without COVID-19. The current diagnostic criterion for COVID-19 is the positive result of a nucleic acid test by real-time reverse transcription polymerase chain reaction (RT-PCR) or next-generation sequencing [14]. However, false-negative results caused by unstable specimen processing are relatively high in clinical practice, which has worsened the spread of the outbreak [15–18]. Moreover, laboratory testing for SARS-CoV-2 requires a rigorous platform, which is not assembled in all hospitals. Thus, this requires specimen transfer, which may delay diagnosis for days. Early and accurate diagnosis is crucial, particularly for critically ill patients who need emergency surgery, and with pneumonia complications. To solve these problems, we hypothesize that a diagnostic model can be developed based on CT imaging and clinical manifestations alone, independent of the nucleic acid test. In this study, we identify the differences in imaging and clinical manifestations between patients with and without COVID-19. We also develop and validate a model for COVID-19 diagnosis based on radiological semantic and clinical features. |
title | Introduction |
p | On January 30, 2020, the World Health Organization (WHO) has declared the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreak as a global health emergency of international concern. This outbreak has infected all provinces of China and rapidly spread to the rest of the world. At the time of writing this article (March 16, 2020), there have been more than 158 countries and territories affected [1]. Whole-genome sequencing and phylogenetic analysis reveal that the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is similar to some beta coronaviruses detected in bats, but it is distinct from severe acute respiratory syndrome coronavirus (SARS-Cov) and Middle East respiratory syndrome coronavirus (MERS-CoV) [2]. |
p | Patients with COVID-19 develop pneumonia with associated symptoms of fever (98%), cough (76%), and myalgia or fatigue (44%) [3]. CT imaging plays a critical role in the diagnosis and the monitoring of disease progression [4–6]. The latest research studies described the characteristic imaging manifestations of COVID-19, including ground-glass opacities (GGO) (57 to 88%), bilateral involvement (76 to 88%), and peripheral distribution (33 to 85%) [7–10]. Other imaging features such as consolidation, cavitation, and interlobular septal thickening are also reported in some patients [11–13]. However, these imaging manifestations of COVID-19 are nonspecific and are difficult to distinguish from other pneumonia. To our knowledge, there have been no studies explicitly comparing imaging and clinical characteristics between pneumonia patients with and without COVID-19. |
p | The current diagnostic criterion for COVID-19 is the positive result of a nucleic acid test by real-time reverse transcription polymerase chain reaction (RT-PCR) or next-generation sequencing [14]. However, false-negative results caused by unstable specimen processing are relatively high in clinical practice, which has worsened the spread of the outbreak [15–18]. Moreover, laboratory testing for SARS-CoV-2 requires a rigorous platform, which is not assembled in all hospitals. Thus, this requires specimen transfer, which may delay diagnosis for days. Early and accurate diagnosis is crucial, particularly for critically ill patients who need emergency surgery, and with pneumonia complications. To solve these problems, we hypothesize that a diagnostic model can be developed based on CT imaging and clinical manifestations alone, independent of the nucleic acid test. |
p | In this study, we identify the differences in imaging and clinical manifestations between patients with and without COVID-19. We also develop and validate a model for COVID-19 diagnosis based on radiological semantic and clinical features. |
sec | Patients and methods Patients Ethical approvals by the institutional review boards were obtained for this retrospective analysis, and the need to obtain informed consent was waived. From January 1 to February 8, 2020, seventy consecutive patients with COVID-19 admitted in 5 independent hospitals from 4 cities were enrolled in this study (mean age, 42.9 years; range, 16–69 years), including 41 men (mean age, 41.8 years; range, 16–69 years) and 29 women (mean age, 44.5 years; range, 16–66 years). All patients were confirmed with SARS-CoV-2 infection by real-time RT-PCR and next-generation sequencing. Of these patients, 24 were from Huizhou City, 25 from Shantou City, 15 from Yongzhou City, and the rest 6 from Meizhou City. At the same period, another 66 pneumonia patients without COVID-19 from Meizhou People’s Hospital were recruited as controls (mean age, 46.7 years; range, 0.3–93 years), including 43 men (mean age, 46.0 years; range, 0.3–93 years) and 23 women (mean age, 48.0 years; range, 1–86 years). All the controls were confirmed with consecutive negative RT-PCR assays. Figure E1 in the Supplementary Material shows the patient recruitment pathway for the control group, along with the inclusion and exclusion criteria. According to previous studies [19–21], whose sample size is comparable with ours, the ratio between primary and validation cohort is 7:3. In this study, a total of 136 patients were divided into primary (n = 98) and validation (n = 38) cohorts, close to 7:3. A total of 19 COVID-19 patients from two hospitals (6 patients from Meizhou People’s Hospital and 13 patients from the First Affiliated Hospital of Shantou University Medical College) and 19 randomly selected controls from Meizhou City were incorporated into the validation cohort. The rest of the patients are incorporated in the primary cohort, including 51 COVID-19 patients from Huizhou, Yongzhou, and Shantou cities and 47 controls from Meizhou City. The primary cohort was utilized to select the most valuable features and build the predictive model, and the validation cohort was used to evaluate and validate the performance of the model. Image and clinical data collection The chest CT imaging data without contrast material enhancement were obtained from multiple hospitals with different CT systems, including GE CT Discovery 750 HD (General Electric Company), SCENARIA 64 CT (Hitachi Medical), Philips Ingenuity CT (PHILIPS), and Siemens SOMATOM Definition AS (Siemens). All images were reconstructed into 1-mm slices with a slice gap of 0.8 mm. Detailed acquisition parameters were summarized in the Supplementary Material (Table E1). The clinical history, nursing records, and laboratory findings were reviewed for all patients. Clinical characteristics, including demographic information, daily body temperature, blood pressure, heart rate, clinical symptoms, and history of exposure to epidemic centers, were collected. Total white blood cell (WBC) counts, lymphocyte counts, ratio of lymphocyte, neutrophil count, ratio of neutrophil, procalcitonin (PCT), C-reactive protein level (CRP), and erythrocyte sedimentation rate (ESR) were measured. All threshold values chosen for laboratory metrics were based on the normal ranges set by each individual hospital. Image analysis For extraction of radiological semantic features, two senior radiologists (D.L. and X.C., more than 15 years of experience) reached a consensus, blinded to clinical and laboratory findings. The radiological semantic features included both qualitative and quantitative imaging features. The lesions in the outer third of the lung were defined as peripheral, and lesions in the inner two-thirds of the lung were defined as central [22]. The progression of COVID-19 lesions within each lung lobe was evaluated by scoring each lobe from 0 to 4 [7], corresponding to normal, 1~25% infection, 26~50% infection, 51~75% infection, and more than 75% infection, respectively. The scores were combined for all five lobes to provide a total score ranging from 0 to 20. A total of 41 radiological features (26 quantitative and 15 qualitative) were extracted for the analysis. The descriptions of radiological semantic features are listed in the Supplementary Material (Table E2). Figure 1 is one example of the evaluation of CT imaging. Fig. 1 A 23-year-old female with a travel history to Wuhan presenting with fever. Axial noncontrast CT image shows a consolidation with ground-glass opacities in the peripheral region by the right upper lobe. Air bronchogram is found in lesion. The maximum diameter of lesion is 2.8 cm. The right upper lobe score is 1 because of the involved lung parenchyma less than 1/4 Clinical and radiological feature selection To obtain the most valuable clinical and radiological semantic features, statistical analysis, univariate analysis, and the least absolute shrinkage and selection operator (LASSO) method were performed. In statistical analysis, the chi-square test, the Kruskal-Wallis H test, and t test were utilized to compare the radiological semantic and clinical features between COVID-19 and non-COVID-19 groups. The features with p value smaller than 0.05 were selected. Then, univariate analysis was performed for clinical and radiological candidate features to determine the COVID-19 risk factors. The features with p value smaller than 0.05 in univariate analysis were also selected. The least absolute shrinkage and selection operator (LASSO) method [23] was utilized to select the most useful features with penalty parameter tuning that was conducted by 10-fold cross-validation based on minimum criteria. Diagnostic models were then constructed by multivariate logistic regression with the selected features. The flowchart of the feature selection process for these models was presented in the Supplementary Material (Fig. E2). Development and validation of the diagnostic model To develop an optimal model, we evaluated 3 models by analyzing (i) the clinical features model (C model), (ii) radiological semantic features model (R model), and (iii) the combination of clinical and radiological semantic features model (CR model) by multivariate logistic regression analysis. The classification performances of the models were evaluated by the area under the receiver operating characteristic (ROC) curve. The area under the curve (AUC), accuracy, sensitivity, and specificity were also calculated. A decision curve analysis was conducted to determine the clinical usefulness of the diagnostic model by quantifying the net benefits at different threshold probabilities in the validation dataset [24]. The development of decision curve was described in the Supplementary Materials. Figure 2 depicts the flowchart of the proposed analysis pipeline described above. We also built a nomogram, which was a quantitative tool to predict the individual probability of infection by COVID-19, based on the multivariate logistic analysis of the CR model with the primary cohort. Depending on the coefficient of the predictive factors in multivariate logistic regression model, all values of each predictive factor were assigned points. A total point was obtained by summing all the points of each predictive factor. The scale also showed the relationship between the total point and the prediction probability in the nomogram. The corresponding calibration curves of the CR model in the primary cohort and validation cohort are shown in the Supplementary Material (Fig. E3). Fig. 2 Workflow of data process and analysis in this study. Radiological semantic features, including qualitative and quantitative imaging features, are extracted from axial lung CT section. The clinical manifestation and laboratory parameters are provided by electronic case system. Statistical analysis is performed for comparing the different features between COVID-19 and non-COVID-19 patients. Univariate analysis, least absolute shrinkage, and selection operator (LASSO) are further performed to determine the COVID-19 risk factors with p < 0.05 in statistical analysis. Three models based on the selected features are established by multivariate logistic regression. These models include radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model). The performance and clinical benefits of the prediction model are assessed by the area under a receiver operating characteristic (ROC) curve and the decision curve, respectively Statistical analysis Statistical analysis was conducted with R software (Version: 3.6.4, http: www.r-project.org/). The reported significance levels were all two-sided, and the statistical significance level was set to 0.05. The multivariate logistic regression analysis was performed with the “stats” package. Nomogram construction was performed using the “rms” package. Decision curve analysis was performed using the “dca. R” package. |
title | Patients and methods |
sec | Patients Ethical approvals by the institutional review boards were obtained for this retrospective analysis, and the need to obtain informed consent was waived. From January 1 to February 8, 2020, seventy consecutive patients with COVID-19 admitted in 5 independent hospitals from 4 cities were enrolled in this study (mean age, 42.9 years; range, 16–69 years), including 41 men (mean age, 41.8 years; range, 16–69 years) and 29 women (mean age, 44.5 years; range, 16–66 years). All patients were confirmed with SARS-CoV-2 infection by real-time RT-PCR and next-generation sequencing. Of these patients, 24 were from Huizhou City, 25 from Shantou City, 15 from Yongzhou City, and the rest 6 from Meizhou City. At the same period, another 66 pneumonia patients without COVID-19 from Meizhou People’s Hospital were recruited as controls (mean age, 46.7 years; range, 0.3–93 years), including 43 men (mean age, 46.0 years; range, 0.3–93 years) and 23 women (mean age, 48.0 years; range, 1–86 years). All the controls were confirmed with consecutive negative RT-PCR assays. Figure E1 in the Supplementary Material shows the patient recruitment pathway for the control group, along with the inclusion and exclusion criteria. According to previous studies [19–21], whose sample size is comparable with ours, the ratio between primary and validation cohort is 7:3. In this study, a total of 136 patients were divided into primary (n = 98) and validation (n = 38) cohorts, close to 7:3. A total of 19 COVID-19 patients from two hospitals (6 patients from Meizhou People’s Hospital and 13 patients from the First Affiliated Hospital of Shantou University Medical College) and 19 randomly selected controls from Meizhou City were incorporated into the validation cohort. The rest of the patients are incorporated in the primary cohort, including 51 COVID-19 patients from Huizhou, Yongzhou, and Shantou cities and 47 controls from Meizhou City. The primary cohort was utilized to select the most valuable features and build the predictive model, and the validation cohort was used to evaluate and validate the performance of the model. |
title | Patients |
p | Ethical approvals by the institutional review boards were obtained for this retrospective analysis, and the need to obtain informed consent was waived. |
p | From January 1 to February 8, 2020, seventy consecutive patients with COVID-19 admitted in 5 independent hospitals from 4 cities were enrolled in this study (mean age, 42.9 years; range, 16–69 years), including 41 men (mean age, 41.8 years; range, 16–69 years) and 29 women (mean age, 44.5 years; range, 16–66 years). All patients were confirmed with SARS-CoV-2 infection by real-time RT-PCR and next-generation sequencing. Of these patients, 24 were from Huizhou City, 25 from Shantou City, 15 from Yongzhou City, and the rest 6 from Meizhou City. At the same period, another 66 pneumonia patients without COVID-19 from Meizhou People’s Hospital were recruited as controls (mean age, 46.7 years; range, 0.3–93 years), including 43 men (mean age, 46.0 years; range, 0.3–93 years) and 23 women (mean age, 48.0 years; range, 1–86 years). All the controls were confirmed with consecutive negative RT-PCR assays. Figure E1 in the Supplementary Material shows the patient recruitment pathway for the control group, along with the inclusion and exclusion criteria. |
p | According to previous studies [19–21], whose sample size is comparable with ours, the ratio between primary and validation cohort is 7:3. In this study, a total of 136 patients were divided into primary (n = 98) and validation (n = 38) cohorts, close to 7:3. A total of 19 COVID-19 patients from two hospitals (6 patients from Meizhou People’s Hospital and 13 patients from the First Affiliated Hospital of Shantou University Medical College) and 19 randomly selected controls from Meizhou City were incorporated into the validation cohort. The rest of the patients are incorporated in the primary cohort, including 51 COVID-19 patients from Huizhou, Yongzhou, and Shantou cities and 47 controls from Meizhou City. The primary cohort was utilized to select the most valuable features and build the predictive model, and the validation cohort was used to evaluate and validate the performance of the model. |
sec | Image and clinical data collection The chest CT imaging data without contrast material enhancement were obtained from multiple hospitals with different CT systems, including GE CT Discovery 750 HD (General Electric Company), SCENARIA 64 CT (Hitachi Medical), Philips Ingenuity CT (PHILIPS), and Siemens SOMATOM Definition AS (Siemens). All images were reconstructed into 1-mm slices with a slice gap of 0.8 mm. Detailed acquisition parameters were summarized in the Supplementary Material (Table E1). The clinical history, nursing records, and laboratory findings were reviewed for all patients. Clinical characteristics, including demographic information, daily body temperature, blood pressure, heart rate, clinical symptoms, and history of exposure to epidemic centers, were collected. Total white blood cell (WBC) counts, lymphocyte counts, ratio of lymphocyte, neutrophil count, ratio of neutrophil, procalcitonin (PCT), C-reactive protein level (CRP), and erythrocyte sedimentation rate (ESR) were measured. All threshold values chosen for laboratory metrics were based on the normal ranges set by each individual hospital. |
title | Image and clinical data collection |
p | The chest CT imaging data without contrast material enhancement were obtained from multiple hospitals with different CT systems, including GE CT Discovery 750 HD (General Electric Company), SCENARIA 64 CT (Hitachi Medical), Philips Ingenuity CT (PHILIPS), and Siemens SOMATOM Definition AS (Siemens). All images were reconstructed into 1-mm slices with a slice gap of 0.8 mm. Detailed acquisition parameters were summarized in the Supplementary Material (Table E1). |
p | The clinical history, nursing records, and laboratory findings were reviewed for all patients. Clinical characteristics, including demographic information, daily body temperature, blood pressure, heart rate, clinical symptoms, and history of exposure to epidemic centers, were collected. Total white blood cell (WBC) counts, lymphocyte counts, ratio of lymphocyte, neutrophil count, ratio of neutrophil, procalcitonin (PCT), C-reactive protein level (CRP), and erythrocyte sedimentation rate (ESR) were measured. All threshold values chosen for laboratory metrics were based on the normal ranges set by each individual hospital. |
sec | Image analysis For extraction of radiological semantic features, two senior radiologists (D.L. and X.C., more than 15 years of experience) reached a consensus, blinded to clinical and laboratory findings. The radiological semantic features included both qualitative and quantitative imaging features. The lesions in the outer third of the lung were defined as peripheral, and lesions in the inner two-thirds of the lung were defined as central [22]. The progression of COVID-19 lesions within each lung lobe was evaluated by scoring each lobe from 0 to 4 [7], corresponding to normal, 1~25% infection, 26~50% infection, 51~75% infection, and more than 75% infection, respectively. The scores were combined for all five lobes to provide a total score ranging from 0 to 20. A total of 41 radiological features (26 quantitative and 15 qualitative) were extracted for the analysis. The descriptions of radiological semantic features are listed in the Supplementary Material (Table E2). Figure 1 is one example of the evaluation of CT imaging. Fig. 1 A 23-year-old female with a travel history to Wuhan presenting with fever. Axial noncontrast CT image shows a consolidation with ground-glass opacities in the peripheral region by the right upper lobe. Air bronchogram is found in lesion. The maximum diameter of lesion is 2.8 cm. The right upper lobe score is 1 because of the involved lung parenchyma less than 1/4 |
title | Image analysis |
p | For extraction of radiological semantic features, two senior radiologists (D.L. and X.C., more than 15 years of experience) reached a consensus, blinded to clinical and laboratory findings. The radiological semantic features included both qualitative and quantitative imaging features. The lesions in the outer third of the lung were defined as peripheral, and lesions in the inner two-thirds of the lung were defined as central [22]. The progression of COVID-19 lesions within each lung lobe was evaluated by scoring each lobe from 0 to 4 [7], corresponding to normal, 1~25% infection, 26~50% infection, 51~75% infection, and more than 75% infection, respectively. The scores were combined for all five lobes to provide a total score ranging from 0 to 20. A total of 41 radiological features (26 quantitative and 15 qualitative) were extracted for the analysis. The descriptions of radiological semantic features are listed in the Supplementary Material (Table E2). Figure 1 is one example of the evaluation of CT imaging. Fig. 1 A 23-year-old female with a travel history to Wuhan presenting with fever. Axial noncontrast CT image shows a consolidation with ground-glass opacities in the peripheral region by the right upper lobe. Air bronchogram is found in lesion. The maximum diameter of lesion is 2.8 cm. The right upper lobe score is 1 because of the involved lung parenchyma less than 1/4 |
figure | Fig. 1 A 23-year-old female with a travel history to Wuhan presenting with fever. Axial noncontrast CT image shows a consolidation with ground-glass opacities in the peripheral region by the right upper lobe. Air bronchogram is found in lesion. The maximum diameter of lesion is 2.8 cm. The right upper lobe score is 1 because of the involved lung parenchyma less than 1/4 |
label | Fig. 1 |
caption | A 23-year-old female with a travel history to Wuhan presenting with fever. Axial noncontrast CT image shows a consolidation with ground-glass opacities in the peripheral region by the right upper lobe. Air bronchogram is found in lesion. The maximum diameter of lesion is 2.8 cm. The right upper lobe score is 1 because of the involved lung parenchyma less than 1/4 |
p | A 23-year-old female with a travel history to Wuhan presenting with fever. Axial noncontrast CT image shows a consolidation with ground-glass opacities in the peripheral region by the right upper lobe. Air bronchogram is found in lesion. The maximum diameter of lesion is 2.8 cm. The right upper lobe score is 1 because of the involved lung parenchyma less than 1/4 |
sec | Clinical and radiological feature selection To obtain the most valuable clinical and radiological semantic features, statistical analysis, univariate analysis, and the least absolute shrinkage and selection operator (LASSO) method were performed. In statistical analysis, the chi-square test, the Kruskal-Wallis H test, and t test were utilized to compare the radiological semantic and clinical features between COVID-19 and non-COVID-19 groups. The features with p value smaller than 0.05 were selected. Then, univariate analysis was performed for clinical and radiological candidate features to determine the COVID-19 risk factors. The features with p value smaller than 0.05 in univariate analysis were also selected. The least absolute shrinkage and selection operator (LASSO) method [23] was utilized to select the most useful features with penalty parameter tuning that was conducted by 10-fold cross-validation based on minimum criteria. Diagnostic models were then constructed by multivariate logistic regression with the selected features. The flowchart of the feature selection process for these models was presented in the Supplementary Material (Fig. E2). |
title | Clinical and radiological feature selection |
p | To obtain the most valuable clinical and radiological semantic features, statistical analysis, univariate analysis, and the least absolute shrinkage and selection operator (LASSO) method were performed. In statistical analysis, the chi-square test, the Kruskal-Wallis H test, and t test were utilized to compare the radiological semantic and clinical features between COVID-19 and non-COVID-19 groups. The features with p value smaller than 0.05 were selected. Then, univariate analysis was performed for clinical and radiological candidate features to determine the COVID-19 risk factors. The features with p value smaller than 0.05 in univariate analysis were also selected. The least absolute shrinkage and selection operator (LASSO) method [23] was utilized to select the most useful features with penalty parameter tuning that was conducted by 10-fold cross-validation based on minimum criteria. Diagnostic models were then constructed by multivariate logistic regression with the selected features. The flowchart of the feature selection process for these models was presented in the Supplementary Material (Fig. E2). |
sec | Development and validation of the diagnostic model To develop an optimal model, we evaluated 3 models by analyzing (i) the clinical features model (C model), (ii) radiological semantic features model (R model), and (iii) the combination of clinical and radiological semantic features model (CR model) by multivariate logistic regression analysis. The classification performances of the models were evaluated by the area under the receiver operating characteristic (ROC) curve. The area under the curve (AUC), accuracy, sensitivity, and specificity were also calculated. A decision curve analysis was conducted to determine the clinical usefulness of the diagnostic model by quantifying the net benefits at different threshold probabilities in the validation dataset [24]. The development of decision curve was described in the Supplementary Materials. Figure 2 depicts the flowchart of the proposed analysis pipeline described above. We also built a nomogram, which was a quantitative tool to predict the individual probability of infection by COVID-19, based on the multivariate logistic analysis of the CR model with the primary cohort. Depending on the coefficient of the predictive factors in multivariate logistic regression model, all values of each predictive factor were assigned points. A total point was obtained by summing all the points of each predictive factor. The scale also showed the relationship between the total point and the prediction probability in the nomogram. The corresponding calibration curves of the CR model in the primary cohort and validation cohort are shown in the Supplementary Material (Fig. E3). Fig. 2 Workflow of data process and analysis in this study. Radiological semantic features, including qualitative and quantitative imaging features, are extracted from axial lung CT section. The clinical manifestation and laboratory parameters are provided by electronic case system. Statistical analysis is performed for comparing the different features between COVID-19 and non-COVID-19 patients. Univariate analysis, least absolute shrinkage, and selection operator (LASSO) are further performed to determine the COVID-19 risk factors with p < 0.05 in statistical analysis. Three models based on the selected features are established by multivariate logistic regression. These models include radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model). The performance and clinical benefits of the prediction model are assessed by the area under a receiver operating characteristic (ROC) curve and the decision curve, respectively |
title | Development and validation of the diagnostic model |
p | To develop an optimal model, we evaluated 3 models by analyzing (i) the clinical features model (C model), (ii) radiological semantic features model (R model), and (iii) the combination of clinical and radiological semantic features model (CR model) by multivariate logistic regression analysis. The classification performances of the models were evaluated by the area under the receiver operating characteristic (ROC) curve. The area under the curve (AUC), accuracy, sensitivity, and specificity were also calculated. A decision curve analysis was conducted to determine the clinical usefulness of the diagnostic model by quantifying the net benefits at different threshold probabilities in the validation dataset [24]. The development of decision curve was described in the Supplementary Materials. Figure 2 depicts the flowchart of the proposed analysis pipeline described above. We also built a nomogram, which was a quantitative tool to predict the individual probability of infection by COVID-19, based on the multivariate logistic analysis of the CR model with the primary cohort. Depending on the coefficient of the predictive factors in multivariate logistic regression model, all values of each predictive factor were assigned points. A total point was obtained by summing all the points of each predictive factor. The scale also showed the relationship between the total point and the prediction probability in the nomogram. The corresponding calibration curves of the CR model in the primary cohort and validation cohort are shown in the Supplementary Material (Fig. E3). Fig. 2 Workflow of data process and analysis in this study. Radiological semantic features, including qualitative and quantitative imaging features, are extracted from axial lung CT section. The clinical manifestation and laboratory parameters are provided by electronic case system. Statistical analysis is performed for comparing the different features between COVID-19 and non-COVID-19 patients. Univariate analysis, least absolute shrinkage, and selection operator (LASSO) are further performed to determine the COVID-19 risk factors with p < 0.05 in statistical analysis. Three models based on the selected features are established by multivariate logistic regression. These models include radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model). The performance and clinical benefits of the prediction model are assessed by the area under a receiver operating characteristic (ROC) curve and the decision curve, respectively |
figure | Fig. 2 Workflow of data process and analysis in this study. Radiological semantic features, including qualitative and quantitative imaging features, are extracted from axial lung CT section. The clinical manifestation and laboratory parameters are provided by electronic case system. Statistical analysis is performed for comparing the different features between COVID-19 and non-COVID-19 patients. Univariate analysis, least absolute shrinkage, and selection operator (LASSO) are further performed to determine the COVID-19 risk factors with p < 0.05 in statistical analysis. Three models based on the selected features are established by multivariate logistic regression. These models include radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model). The performance and clinical benefits of the prediction model are assessed by the area under a receiver operating characteristic (ROC) curve and the decision curve, respectively |
label | Fig. 2 |
caption | Workflow of data process and analysis in this study. Radiological semantic features, including qualitative and quantitative imaging features, are extracted from axial lung CT section. The clinical manifestation and laboratory parameters are provided by electronic case system. Statistical analysis is performed for comparing the different features between COVID-19 and non-COVID-19 patients. Univariate analysis, least absolute shrinkage, and selection operator (LASSO) are further performed to determine the COVID-19 risk factors with p < 0.05 in statistical analysis. Three models based on the selected features are established by multivariate logistic regression. These models include radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model). The performance and clinical benefits of the prediction model are assessed by the area under a receiver operating characteristic (ROC) curve and the decision curve, respectively |
p | Workflow of data process and analysis in this study. Radiological semantic features, including qualitative and quantitative imaging features, are extracted from axial lung CT section. The clinical manifestation and laboratory parameters are provided by electronic case system. Statistical analysis is performed for comparing the different features between COVID-19 and non-COVID-19 patients. Univariate analysis, least absolute shrinkage, and selection operator (LASSO) are further performed to determine the COVID-19 risk factors with p < 0.05 in statistical analysis. Three models based on the selected features are established by multivariate logistic regression. These models include radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model). The performance and clinical benefits of the prediction model are assessed by the area under a receiver operating characteristic (ROC) curve and the decision curve, respectively |
sec | Statistical analysis Statistical analysis was conducted with R software (Version: 3.6.4, http: www.r-project.org/). The reported significance levels were all two-sided, and the statistical significance level was set to 0.05. The multivariate logistic regression analysis was performed with the “stats” package. Nomogram construction was performed using the “rms” package. Decision curve analysis was performed using the “dca. R” package. |
title | Statistical analysis |
p | Statistical analysis was conducted with R software (Version: 3.6.4, http: www.r-project.org/). The reported significance levels were all two-sided, and the statistical significance level was set to 0.05. The multivariate logistic regression analysis was performed with the “stats” package. Nomogram construction was performed using the “rms” package. Decision curve analysis was performed using the “dca. R” package. |
sec | Results Imaging and clinical manifestations between groups The differences between patients with and without COVID-19 for all 67 features (41 imaging and 26 critical clinical features) are shown in Tables 1 and 2 and the Supplementary Materials (Tables E3 and E4). The differences between the primary cohort and validation cohort for the same features are shown in the Supplementary Materials (Tables E5 and E6). All characteristics except fatigue and white blood cell count in the CR model presented no significant difference between the primary and validation cohorts. A total of 1745 lesions were identified, with 1062 from the COVID-19 group and 683 from the non-COVID-19 group. Table 1 Radiological semantic features of patients in COVID-19 and non-COVID-19 Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value Number of pure GGO Total# 1.00 (0.00, 5.05) 3.50 (0.95, 8.05) 0.018b* Peripheral area# 1.00 (0.00, 4.05) 2.00 (0.00, 6.05) 0.032b* Central/both peripheral and central area# 0.00 (0.00, 0.00) 0.00 (0.00, 2.00) 0.001b* Number of mixed GGO Total# 1.00 (0.00, 3.05) 3.00 (1.00, 9.00) 0.001b* Peripheral area# 0.00 (0.00, 2.00) 2.50 (1.00, 6.00) < 0.001b* Central/both peripheral and central area# 0.00 (0.00, 1.05) 0.00 (0.00, 2.00) 0.657b Total number of consolidation Consolidation# 1.00 (0.00, 3.00) 0.00 (0.00, 0.05) 0.001b* Pure solid nodules# 0.00 (0.00, 0.00) 0.00 (0.00, 0.00) 0.309b Solid nodules with GGO# 0.00 (0.00, 0.00) 0.00 (0.00, 1.00) 0.033b* Total number of lesions Peripheral area# 5.00 (2.00, 9.05) 7.00 (2.00, 13.00) 0.112b Central area# 0.00 (0.00, 3.00) 0.00 (0.00, 1.05) 0.960b Both peripheral and central area# 0.00 (0.00, 2.00) 0.00 (0.00, 2.05) 0.582b Interlobular septal thickening 0.009a* Negative 44 (66.67%) 31 (44.29%) Positive 22 (33.33%) 39 (55.71%) Crazy paving pattern < 0.001a* Negative 60 (90.91%) 32 (45.71%) Positive 6 (9.09%) 38 (54.29%) Tree-in-bud sign < 0.001a* Negative 37 (56.06%) 61 (87.14%) Positive 29 (43.94%) 9 (12.86%) Pleural thickening 0.030a* Negative 46 (69.70%) 36 (51.43%) Positive 20 (30.30%) 34 (48.57%) Offending vessel augmentation in lesions < 0.001a* Negative 55 (83.33%) 17 (24.29%) Positive 11 (16.67%) 53 (75.71%) GGO ground-glass opacities #Results are median with interquartile range in parentheses, and the remainder results are measurements with corresponding ratio in parentheses *Data with statistical significance. pa: chi-square test, pb: Student’s t test Table 2 Clinical features of patients in COVID-19 and non-COVID-19 Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value Sex Male# 43 (65.15%) 41 (58.57%) 0.430a Female# 23 (34.85%) 29 (41.43%) Age (years) 46.73 ± 25.00 42.93 ± 13.32 0.275b Vital signs Systolic blood pressure (mmHg) 126.92 ± 23.07 127.07 ± 15.16 0.965b Diastolic blood pressure (mmHg) 77.74 ± 15.72 80.39 ± 10.51 0.254b Respiration rate (bpm) 25.20 ± 7.29 19.86 ± 1.90 < 0.001b* Heart rate (bpm) 101.59 ± 20.36 86.06 ± 13.34 < 0.001b* Temperature (°C) 37.61 ± 1.06 37.12 ± 0.83 0.003b* Signs Dry cough# 56 (84.85%) 48 (68.57%) 0.025a* Fatigue# 8 (12.12%) 22 (31.43%) 0.007a* Sore throat# 6 (9.09%) 9 (12.86%) 0.483a Stuffy# 4 (6.06%) 2 (2.86%) 0.623a Runny nose# 3 (4.55%) 3 (4.29%) 0.731a White blood cell count (× 109/L) 11.48 ± 5.36 5.27 ± 2.33 < 0.001b* White blood cell count category < 0.001c* Low# 0 (0.00%) 2 (2.86%) Normal# 27 (40.91%) 63 (90.00%) High# 39 (59.09%) 5 (7.14%) Lymphocyte count (× 109/L) 1.57 ± 1.33 1.25 ± 0.68 0.086b Lymphocyte count category < 0.001c* Low# 24 (36.36%) 32 (45.71%) Normal# 35 (53.03%) 37 (52.86%) High# 7 (10.61%) 1 (1.43%) Neutrophil count (× 109/L) 8.97 ± 4.90 3.53 ± 2.17 < 0.001b* Neutrophil count category < 0.001c* Low# 3 (4.55%) 8 (11.43%) Normal# 23 (34.85%) 59 (84.29%) High# 40 (60.61%) 3 (4.29%) C-reactive protein (mg/L) 69.30 ± 65.88 26.37 ± 30.97 < 0.001b* Procalcitonin (ng/mL) 3.36 ± 8.98 0.26 ± 0.84 0.007b* *Data with statistical significance. pa: chi-square test, pb: Student’s t test. pc: Kruskal-Wallis H test #Results are measurements with corresponding ratio in parentheses For imaging manifestations, 7 patients in the COVID-19 group showed normal chest CT (10%). COVID-19 patients have a greater number of pure GGO and mixed GGO than non-COVID-19 patients (p = 0.018 and p = 0.001, respectively). For pure GGO lesions, the differences are significant both in peripheral (p = 0.032) and in central areas (p = 0.001). However, the number of mixed GGO is mainly distributed at the periphery in COVID-19 patients (p < 0.001), with no statistical difference in the central area. The consolidation lesions without GGO occurred less in COVID-19 patients (p = 0.001). More lesions are between 1 and 3 cm (p = 0.027), and fewer lesions are larger than half of the lung segment (p = 0.017) in COVID-19 patients. Other significant differences between the two groups include the pleural traction sign (p = 0.019), bronchial wall thickening (p < 0.001), interlobular septal thickening (p = 0.009), crazy paving (p < 0.001), tree-in-bud (p < 0.001), pleural effusions (p < 0.001), pleural thickening (p = 0.030), and the offending vessel augmentation in lesions (p < 0.001). The lung score presents no significant difference between the COVID-19 and non-COVID-19 groups. Comparison of clinical features between the two groups of patients with and without COVID-19 is reported in Table 2. There is no significant difference in age and sex between the two groups. Significant differences are found in common symptoms between groups, including fever (p = 0.003), dry cough (p = 0.025), and fatigue (p = 0.007). The respiration rate and heart rate also show significant differences between the two groups (both p < 0.001). Compared with non-COVID-19 pneumonia, the reduction of the WBC count is more pronounced in COVID-19 patients (p < 0.001). The ratio of lymphocyte and ratio of neutrophil also show a significant difference between COVID-19 and non-COVID-19 groups. Although lymphopenia was observed in 32 COVID-19 patients (45.71%), it is not statistically different compared with that in the non-COVID-19 group. C-creative protein (CRP) level and procalcitonin level are also significantly different between the two groups (p < 0.001 and p = 0.007, respectively). Most COVID-19 patients present normal procalcitonin level (82.86%). Clinical and radiological feature selection Of the features, 18 radiological features and 17 clinical features were selected to form the predictors based on the result from Tables 1 and 2. Table 3 lists the features selected by univariate analysis and LASSO. Table 3 Selected features in C, R, and CR models Model and individual features Coefficients R, n = 8 (41)* Intercept − 0.307 Total number of mixed GGO in peripheral area 0.359 Total number of consolidation − 1.262 Total number of solid nodules with ground-glass opacities 0.452 Interlobular septal thickening − 5.559 Crazy paving pattern 3.566 Tree-in-bud − 2.548 Pleural thickening 3.265 Offending vessel augmentation in lesions 5.504 C, n = 7 (26)* Intercept 29.273 Respiration − 0.359 Heart rate − 0.054 Temperature − 0.289 White blood cell count − 0.175 Cough − 1.866 Fatigue 2.855 Lymphocyte count category − 0.028 CR, n = 10 (67)* Intercept 45.117 Total number of mixed GGO in peripheral area 0.108 Tree-in-bud − 1.853 Offending vessel augmentation in lesions 6.000 Respiration − 0.583 Heart ratio − 0.084 Temperature − 0.536 White blood cell count − 0.471 Cough − 0.997 Fatigue − 0.228 Lymphocyte count category − 2.177 C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively *n means corresponding selected features, and data in parentheses are total features. Coefficients: the estimate value of each feature in multivariate logistic regression model by “glm” package in R Model development and validation The prediction models based on (i) clinical features (C model), (ii) radiological features (R model), and (iii) the combination of clinical features and radiological features (CR model) were developed. ROC analyses for the primary and validation cohort are shown in Table 4 and Fig. 3. The CR model yielded a maximum AUC of 0.986 (95% CI 0.966~1.000) in the primary cohort with the highest accuracy and specificity, which was 0.936 (95% CI 0.866~1.000) in the validation cohort. The AUC for the C model was 0.952 (95% CI 0.988~0.915) and 0.967 (95% CI 0.919~1.000) in the primary and validation cohorts, respectively. For the R model, the AUC of the two cohorts was 0.969 (95% CI 0.940~0.997) and 0.809 (95% CI 0.669~0.948), respectively. Table 4 Performance of the individualized prediction models Primary cohort (n = 98) Validation cohort (n = 38) Models AUC 95% CI Accuracy Specificity Sensitivity AUC 95% CI Accuracy Specificity Sensitivity C model 0.952 0.915~0.988 0.888 0.894 0.882 0.967 0.919~1.000 0.868 0.859 0.842 R model 0.969 0.940~0.997 0.929 0.851 1.000 0.809 0.669~0.948 0.684 0.368 1.000 CR model 0.986 0.966~1.000 0.959 0.957 0.961 0.936 0.866~1.000 0.763 0.789 0.737 C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively. CI confidence interval Fig. 3 ROC of the three models in primary and validation cohort curves. Comparison of receiver operating characteristic (ROC) curves among the radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model) for the diagnosis of COVID-19 in the primary (a) and validation (b) cohorts To determine the clinical usefulness of the diagnostic model, we developed the decision curve (Fig. 4), which showed better performances for the CR model compared with that for the C model and the R model. Across the majority of the range of reasonable threshold probabilities, the decision curve analysis showed that the CR model had a higher overall benefit than the C model and R model. Fig. 4 Decision curve analysis for each model in the primary dataset. The y-axis measures the net benefit, which is calculated by summing the benefits (true-positive findings) and subtracting the harms (false-positive findings), weighting the latter by a factor related to the relative harm of undetected metastasis compared with the harm of unnecessary treatment. The decision curve shows that if the threshold probability is over 10%, the application of the combination of clinical and radiological model (CR model) to diagnose COVID-19 adds more benefit than the clinical model (C model) and radiological model (R model) The nomogram (Fig. 5) was developed by the CR model in the primary cohort, with the factors of the total number of mixed GGO in peripheral area (TN_Mixed_GGO_IP), tree-in-bud, offending vessel augmentation in lesions (OVAIL), respiration, heart ratio, temperature, white blood cell count, cough, fatigue and lymphocyte count category incorporated. The total points were calculated by summing the points identified on the “points” scale for each factor. By comparing the “total points” scale and the “probability” scale, the individual probability of COVID-19 infection could be obtained. Fig. 5 Nomogram of the CR model in the primary cohort. TN_Mixed_GGO_IP represented the total number of mixed GGO in peripheral area. AVAIL represented offending vessel segmentation in lesions. N was a negative result, and P was a positive result. Norm represented normal. Note that in probability scale, 0 = non-COVID-19, 1 = COVID-19 |
title | Results |
sec | Imaging and clinical manifestations between groups The differences between patients with and without COVID-19 for all 67 features (41 imaging and 26 critical clinical features) are shown in Tables 1 and 2 and the Supplementary Materials (Tables E3 and E4). The differences between the primary cohort and validation cohort for the same features are shown in the Supplementary Materials (Tables E5 and E6). All characteristics except fatigue and white blood cell count in the CR model presented no significant difference between the primary and validation cohorts. A total of 1745 lesions were identified, with 1062 from the COVID-19 group and 683 from the non-COVID-19 group. Table 1 Radiological semantic features of patients in COVID-19 and non-COVID-19 Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value Number of pure GGO Total# 1.00 (0.00, 5.05) 3.50 (0.95, 8.05) 0.018b* Peripheral area# 1.00 (0.00, 4.05) 2.00 (0.00, 6.05) 0.032b* Central/both peripheral and central area# 0.00 (0.00, 0.00) 0.00 (0.00, 2.00) 0.001b* Number of mixed GGO Total# 1.00 (0.00, 3.05) 3.00 (1.00, 9.00) 0.001b* Peripheral area# 0.00 (0.00, 2.00) 2.50 (1.00, 6.00) < 0.001b* Central/both peripheral and central area# 0.00 (0.00, 1.05) 0.00 (0.00, 2.00) 0.657b Total number of consolidation Consolidation# 1.00 (0.00, 3.00) 0.00 (0.00, 0.05) 0.001b* Pure solid nodules# 0.00 (0.00, 0.00) 0.00 (0.00, 0.00) 0.309b Solid nodules with GGO# 0.00 (0.00, 0.00) 0.00 (0.00, 1.00) 0.033b* Total number of lesions Peripheral area# 5.00 (2.00, 9.05) 7.00 (2.00, 13.00) 0.112b Central area# 0.00 (0.00, 3.00) 0.00 (0.00, 1.05) 0.960b Both peripheral and central area# 0.00 (0.00, 2.00) 0.00 (0.00, 2.05) 0.582b Interlobular septal thickening 0.009a* Negative 44 (66.67%) 31 (44.29%) Positive 22 (33.33%) 39 (55.71%) Crazy paving pattern < 0.001a* Negative 60 (90.91%) 32 (45.71%) Positive 6 (9.09%) 38 (54.29%) Tree-in-bud sign < 0.001a* Negative 37 (56.06%) 61 (87.14%) Positive 29 (43.94%) 9 (12.86%) Pleural thickening 0.030a* Negative 46 (69.70%) 36 (51.43%) Positive 20 (30.30%) 34 (48.57%) Offending vessel augmentation in lesions < 0.001a* Negative 55 (83.33%) 17 (24.29%) Positive 11 (16.67%) 53 (75.71%) GGO ground-glass opacities #Results are median with interquartile range in parentheses, and the remainder results are measurements with corresponding ratio in parentheses *Data with statistical significance. pa: chi-square test, pb: Student’s t test Table 2 Clinical features of patients in COVID-19 and non-COVID-19 Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value Sex Male# 43 (65.15%) 41 (58.57%) 0.430a Female# 23 (34.85%) 29 (41.43%) Age (years) 46.73 ± 25.00 42.93 ± 13.32 0.275b Vital signs Systolic blood pressure (mmHg) 126.92 ± 23.07 127.07 ± 15.16 0.965b Diastolic blood pressure (mmHg) 77.74 ± 15.72 80.39 ± 10.51 0.254b Respiration rate (bpm) 25.20 ± 7.29 19.86 ± 1.90 < 0.001b* Heart rate (bpm) 101.59 ± 20.36 86.06 ± 13.34 < 0.001b* Temperature (°C) 37.61 ± 1.06 37.12 ± 0.83 0.003b* Signs Dry cough# 56 (84.85%) 48 (68.57%) 0.025a* Fatigue# 8 (12.12%) 22 (31.43%) 0.007a* Sore throat# 6 (9.09%) 9 (12.86%) 0.483a Stuffy# 4 (6.06%) 2 (2.86%) 0.623a Runny nose# 3 (4.55%) 3 (4.29%) 0.731a White blood cell count (× 109/L) 11.48 ± 5.36 5.27 ± 2.33 < 0.001b* White blood cell count category < 0.001c* Low# 0 (0.00%) 2 (2.86%) Normal# 27 (40.91%) 63 (90.00%) High# 39 (59.09%) 5 (7.14%) Lymphocyte count (× 109/L) 1.57 ± 1.33 1.25 ± 0.68 0.086b Lymphocyte count category < 0.001c* Low# 24 (36.36%) 32 (45.71%) Normal# 35 (53.03%) 37 (52.86%) High# 7 (10.61%) 1 (1.43%) Neutrophil count (× 109/L) 8.97 ± 4.90 3.53 ± 2.17 < 0.001b* Neutrophil count category < 0.001c* Low# 3 (4.55%) 8 (11.43%) Normal# 23 (34.85%) 59 (84.29%) High# 40 (60.61%) 3 (4.29%) C-reactive protein (mg/L) 69.30 ± 65.88 26.37 ± 30.97 < 0.001b* Procalcitonin (ng/mL) 3.36 ± 8.98 0.26 ± 0.84 0.007b* *Data with statistical significance. pa: chi-square test, pb: Student’s t test. pc: Kruskal-Wallis H test #Results are measurements with corresponding ratio in parentheses For imaging manifestations, 7 patients in the COVID-19 group showed normal chest CT (10%). COVID-19 patients have a greater number of pure GGO and mixed GGO than non-COVID-19 patients (p = 0.018 and p = 0.001, respectively). For pure GGO lesions, the differences are significant both in peripheral (p = 0.032) and in central areas (p = 0.001). However, the number of mixed GGO is mainly distributed at the periphery in COVID-19 patients (p < 0.001), with no statistical difference in the central area. The consolidation lesions without GGO occurred less in COVID-19 patients (p = 0.001). More lesions are between 1 and 3 cm (p = 0.027), and fewer lesions are larger than half of the lung segment (p = 0.017) in COVID-19 patients. Other significant differences between the two groups include the pleural traction sign (p = 0.019), bronchial wall thickening (p < 0.001), interlobular septal thickening (p = 0.009), crazy paving (p < 0.001), tree-in-bud (p < 0.001), pleural effusions (p < 0.001), pleural thickening (p = 0.030), and the offending vessel augmentation in lesions (p < 0.001). The lung score presents no significant difference between the COVID-19 and non-COVID-19 groups. Comparison of clinical features between the two groups of patients with and without COVID-19 is reported in Table 2. There is no significant difference in age and sex between the two groups. Significant differences are found in common symptoms between groups, including fever (p = 0.003), dry cough (p = 0.025), and fatigue (p = 0.007). The respiration rate and heart rate also show significant differences between the two groups (both p < 0.001). Compared with non-COVID-19 pneumonia, the reduction of the WBC count is more pronounced in COVID-19 patients (p < 0.001). The ratio of lymphocyte and ratio of neutrophil also show a significant difference between COVID-19 and non-COVID-19 groups. Although lymphopenia was observed in 32 COVID-19 patients (45.71%), it is not statistically different compared with that in the non-COVID-19 group. C-creative protein (CRP) level and procalcitonin level are also significantly different between the two groups (p < 0.001 and p = 0.007, respectively). Most COVID-19 patients present normal procalcitonin level (82.86%). |
title | Imaging and clinical manifestations between groups |
p | The differences between patients with and without COVID-19 for all 67 features (41 imaging and 26 critical clinical features) are shown in Tables 1 and 2 and the Supplementary Materials (Tables E3 and E4). The differences between the primary cohort and validation cohort for the same features are shown in the Supplementary Materials (Tables E5 and E6). All characteristics except fatigue and white blood cell count in the CR model presented no significant difference between the primary and validation cohorts. A total of 1745 lesions were identified, with 1062 from the COVID-19 group and 683 from the non-COVID-19 group. Table 1 Radiological semantic features of patients in COVID-19 and non-COVID-19 Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value Number of pure GGO Total# 1.00 (0.00, 5.05) 3.50 (0.95, 8.05) 0.018b* Peripheral area# 1.00 (0.00, 4.05) 2.00 (0.00, 6.05) 0.032b* Central/both peripheral and central area# 0.00 (0.00, 0.00) 0.00 (0.00, 2.00) 0.001b* Number of mixed GGO Total# 1.00 (0.00, 3.05) 3.00 (1.00, 9.00) 0.001b* Peripheral area# 0.00 (0.00, 2.00) 2.50 (1.00, 6.00) < 0.001b* Central/both peripheral and central area# 0.00 (0.00, 1.05) 0.00 (0.00, 2.00) 0.657b Total number of consolidation Consolidation# 1.00 (0.00, 3.00) 0.00 (0.00, 0.05) 0.001b* Pure solid nodules# 0.00 (0.00, 0.00) 0.00 (0.00, 0.00) 0.309b Solid nodules with GGO# 0.00 (0.00, 0.00) 0.00 (0.00, 1.00) 0.033b* Total number of lesions Peripheral area# 5.00 (2.00, 9.05) 7.00 (2.00, 13.00) 0.112b Central area# 0.00 (0.00, 3.00) 0.00 (0.00, 1.05) 0.960b Both peripheral and central area# 0.00 (0.00, 2.00) 0.00 (0.00, 2.05) 0.582b Interlobular septal thickening 0.009a* Negative 44 (66.67%) 31 (44.29%) Positive 22 (33.33%) 39 (55.71%) Crazy paving pattern < 0.001a* Negative 60 (90.91%) 32 (45.71%) Positive 6 (9.09%) 38 (54.29%) Tree-in-bud sign < 0.001a* Negative 37 (56.06%) 61 (87.14%) Positive 29 (43.94%) 9 (12.86%) Pleural thickening 0.030a* Negative 46 (69.70%) 36 (51.43%) Positive 20 (30.30%) 34 (48.57%) Offending vessel augmentation in lesions < 0.001a* Negative 55 (83.33%) 17 (24.29%) Positive 11 (16.67%) 53 (75.71%) GGO ground-glass opacities #Results are median with interquartile range in parentheses, and the remainder results are measurements with corresponding ratio in parentheses *Data with statistical significance. pa: chi-square test, pb: Student’s t test Table 2 Clinical features of patients in COVID-19 and non-COVID-19 Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value Sex Male# 43 (65.15%) 41 (58.57%) 0.430a Female# 23 (34.85%) 29 (41.43%) Age (years) 46.73 ± 25.00 42.93 ± 13.32 0.275b Vital signs Systolic blood pressure (mmHg) 126.92 ± 23.07 127.07 ± 15.16 0.965b Diastolic blood pressure (mmHg) 77.74 ± 15.72 80.39 ± 10.51 0.254b Respiration rate (bpm) 25.20 ± 7.29 19.86 ± 1.90 < 0.001b* Heart rate (bpm) 101.59 ± 20.36 86.06 ± 13.34 < 0.001b* Temperature (°C) 37.61 ± 1.06 37.12 ± 0.83 0.003b* Signs Dry cough# 56 (84.85%) 48 (68.57%) 0.025a* Fatigue# 8 (12.12%) 22 (31.43%) 0.007a* Sore throat# 6 (9.09%) 9 (12.86%) 0.483a Stuffy# 4 (6.06%) 2 (2.86%) 0.623a Runny nose# 3 (4.55%) 3 (4.29%) 0.731a White blood cell count (× 109/L) 11.48 ± 5.36 5.27 ± 2.33 < 0.001b* White blood cell count category < 0.001c* Low# 0 (0.00%) 2 (2.86%) Normal# 27 (40.91%) 63 (90.00%) High# 39 (59.09%) 5 (7.14%) Lymphocyte count (× 109/L) 1.57 ± 1.33 1.25 ± 0.68 0.086b Lymphocyte count category < 0.001c* Low# 24 (36.36%) 32 (45.71%) Normal# 35 (53.03%) 37 (52.86%) High# 7 (10.61%) 1 (1.43%) Neutrophil count (× 109/L) 8.97 ± 4.90 3.53 ± 2.17 < 0.001b* Neutrophil count category < 0.001c* Low# 3 (4.55%) 8 (11.43%) Normal# 23 (34.85%) 59 (84.29%) High# 40 (60.61%) 3 (4.29%) C-reactive protein (mg/L) 69.30 ± 65.88 26.37 ± 30.97 < 0.001b* Procalcitonin (ng/mL) 3.36 ± 8.98 0.26 ± 0.84 0.007b* *Data with statistical significance. pa: chi-square test, pb: Student’s t test. pc: Kruskal-Wallis H test #Results are measurements with corresponding ratio in parentheses |
table-wrap | Table 1 Radiological semantic features of patients in COVID-19 and non-COVID-19 Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value Number of pure GGO Total# 1.00 (0.00, 5.05) 3.50 (0.95, 8.05) 0.018b* Peripheral area# 1.00 (0.00, 4.05) 2.00 (0.00, 6.05) 0.032b* Central/both peripheral and central area# 0.00 (0.00, 0.00) 0.00 (0.00, 2.00) 0.001b* Number of mixed GGO Total# 1.00 (0.00, 3.05) 3.00 (1.00, 9.00) 0.001b* Peripheral area# 0.00 (0.00, 2.00) 2.50 (1.00, 6.00) < 0.001b* Central/both peripheral and central area# 0.00 (0.00, 1.05) 0.00 (0.00, 2.00) 0.657b Total number of consolidation Consolidation# 1.00 (0.00, 3.00) 0.00 (0.00, 0.05) 0.001b* Pure solid nodules# 0.00 (0.00, 0.00) 0.00 (0.00, 0.00) 0.309b Solid nodules with GGO# 0.00 (0.00, 0.00) 0.00 (0.00, 1.00) 0.033b* Total number of lesions Peripheral area# 5.00 (2.00, 9.05) 7.00 (2.00, 13.00) 0.112b Central area# 0.00 (0.00, 3.00) 0.00 (0.00, 1.05) 0.960b Both peripheral and central area# 0.00 (0.00, 2.00) 0.00 (0.00, 2.05) 0.582b Interlobular septal thickening 0.009a* Negative 44 (66.67%) 31 (44.29%) Positive 22 (33.33%) 39 (55.71%) Crazy paving pattern < 0.001a* Negative 60 (90.91%) 32 (45.71%) Positive 6 (9.09%) 38 (54.29%) Tree-in-bud sign < 0.001a* Negative 37 (56.06%) 61 (87.14%) Positive 29 (43.94%) 9 (12.86%) Pleural thickening 0.030a* Negative 46 (69.70%) 36 (51.43%) Positive 20 (30.30%) 34 (48.57%) Offending vessel augmentation in lesions < 0.001a* Negative 55 (83.33%) 17 (24.29%) Positive 11 (16.67%) 53 (75.71%) GGO ground-glass opacities #Results are median with interquartile range in parentheses, and the remainder results are measurements with corresponding ratio in parentheses *Data with statistical significance. pa: chi-square test, pb: Student’s t test |
label | Table 1 |
caption | Radiological semantic features of patients in COVID-19 and non-COVID-19 |
p | Radiological semantic features of patients in COVID-19 and non-COVID-19 |
table | Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value Number of pure GGO Total# 1.00 (0.00, 5.05) 3.50 (0.95, 8.05) 0.018b* Peripheral area# 1.00 (0.00, 4.05) 2.00 (0.00, 6.05) 0.032b* Central/both peripheral and central area# 0.00 (0.00, 0.00) 0.00 (0.00, 2.00) 0.001b* Number of mixed GGO Total# 1.00 (0.00, 3.05) 3.00 (1.00, 9.00) 0.001b* Peripheral area# 0.00 (0.00, 2.00) 2.50 (1.00, 6.00) < 0.001b* Central/both peripheral and central area# 0.00 (0.00, 1.05) 0.00 (0.00, 2.00) 0.657b Total number of consolidation Consolidation# 1.00 (0.00, 3.00) 0.00 (0.00, 0.05) 0.001b* Pure solid nodules# 0.00 (0.00, 0.00) 0.00 (0.00, 0.00) 0.309b Solid nodules with GGO# 0.00 (0.00, 0.00) 0.00 (0.00, 1.00) 0.033b* Total number of lesions Peripheral area# 5.00 (2.00, 9.05) 7.00 (2.00, 13.00) 0.112b Central area# 0.00 (0.00, 3.00) 0.00 (0.00, 1.05) 0.960b Both peripheral and central area# 0.00 (0.00, 2.00) 0.00 (0.00, 2.05) 0.582b Interlobular septal thickening 0.009a* Negative 44 (66.67%) 31 (44.29%) Positive 22 (33.33%) 39 (55.71%) Crazy paving pattern < 0.001a* Negative 60 (90.91%) 32 (45.71%) Positive 6 (9.09%) 38 (54.29%) Tree-in-bud sign < 0.001a* Negative 37 (56.06%) 61 (87.14%) Positive 29 (43.94%) 9 (12.86%) Pleural thickening 0.030a* Negative 46 (69.70%) 36 (51.43%) Positive 20 (30.30%) 34 (48.57%) Offending vessel augmentation in lesions < 0.001a* Negative 55 (83.33%) 17 (24.29%) Positive 11 (16.67%) 53 (75.71%) |
tr | Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value |
th | Feature |
th | Non-COVID-19 (n = 66) |
th | COVID-19 (n = 70) |
th | p value |
tr | Number of pure GGO |
td | Number of pure GGO |
tr | Total# 1.00 (0.00, 5.05) 3.50 (0.95, 8.05) 0.018b* |
td | Total# |
td | 1.00 (0.00, 5.05) |
td | 3.50 (0.95, 8.05) |
td | 0.018b* |
tr | Peripheral area# 1.00 (0.00, 4.05) 2.00 (0.00, 6.05) 0.032b* |
td | Peripheral area# |
td | 1.00 (0.00, 4.05) |
td | 2.00 (0.00, 6.05) |
td | 0.032b* |
tr | Central/both peripheral and central area# 0.00 (0.00, 0.00) 0.00 (0.00, 2.00) 0.001b* |
td | Central/both peripheral and central area# |
td | 0.00 (0.00, 0.00) |
td | 0.00 (0.00, 2.00) |
td | 0.001b* |
tr | Number of mixed GGO |
td | Number of mixed GGO |
tr | Total# 1.00 (0.00, 3.05) 3.00 (1.00, 9.00) 0.001b* |
td | Total# |
td | 1.00 (0.00, 3.05) |
td | 3.00 (1.00, 9.00) |
td | 0.001b* |
tr | Peripheral area# 0.00 (0.00, 2.00) 2.50 (1.00, 6.00) < 0.001b* |
td | Peripheral area# |
td | 0.00 (0.00, 2.00) |
td | 2.50 (1.00, 6.00) |
td | < 0.001b* |
tr | Central/both peripheral and central area# 0.00 (0.00, 1.05) 0.00 (0.00, 2.00) 0.657b |
td | Central/both peripheral and central area# |
td | 0.00 (0.00, 1.05) |
td | 0.00 (0.00, 2.00) |
td | 0.657b |
tr | Total number of consolidation |
td | Total number of consolidation |
tr | Consolidation# 1.00 (0.00, 3.00) 0.00 (0.00, 0.05) 0.001b* |
td | Consolidation# |
td | 1.00 (0.00, 3.00) |
td | 0.00 (0.00, 0.05) |
td | 0.001b* |
tr | Pure solid nodules# 0.00 (0.00, 0.00) 0.00 (0.00, 0.00) 0.309b |
td | Pure solid nodules# |
td | 0.00 (0.00, 0.00) |
td | 0.00 (0.00, 0.00) |
td | 0.309b |
tr | Solid nodules with GGO# 0.00 (0.00, 0.00) 0.00 (0.00, 1.00) 0.033b* |
td | Solid nodules with GGO# |
td | 0.00 (0.00, 0.00) |
td | 0.00 (0.00, 1.00) |
td | 0.033b* |
tr | Total number of lesions |
td | Total number of lesions |
tr | Peripheral area# 5.00 (2.00, 9.05) 7.00 (2.00, 13.00) 0.112b |
td | Peripheral area# |
td | 5.00 (2.00, 9.05) |
td | 7.00 (2.00, 13.00) |
td | 0.112b |
tr | Central area# 0.00 (0.00, 3.00) 0.00 (0.00, 1.05) 0.960b |
td | Central area# |
td | 0.00 (0.00, 3.00) |
td | 0.00 (0.00, 1.05) |
td | 0.960b |
tr | Both peripheral and central area# 0.00 (0.00, 2.00) 0.00 (0.00, 2.05) 0.582b |
td | Both peripheral and central area# |
td | 0.00 (0.00, 2.00) |
td | 0.00 (0.00, 2.05) |
td | 0.582b |
tr | Interlobular septal thickening 0.009a* |
td | Interlobular septal thickening |
td | 0.009a* |
tr | Negative 44 (66.67%) 31 (44.29%) |
td | Negative |
td | 44 (66.67%) |
td | 31 (44.29%) |
tr | Positive 22 (33.33%) 39 (55.71%) |
td | Positive |
td | 22 (33.33%) |
td | 39 (55.71%) |
tr | Crazy paving pattern < 0.001a* |
td | Crazy paving pattern |
td | < 0.001a* |
tr | Negative 60 (90.91%) 32 (45.71%) |
td | Negative |
td | 60 (90.91%) |
td | 32 (45.71%) |
tr | Positive 6 (9.09%) 38 (54.29%) |
td | Positive |
td | 6 (9.09%) |
td | 38 (54.29%) |
tr | Tree-in-bud sign < 0.001a* |
td | Tree-in-bud sign |
td | < 0.001a* |
tr | Negative 37 (56.06%) 61 (87.14%) |
td | Negative |
td | 37 (56.06%) |
td | 61 (87.14%) |
tr | Positive 29 (43.94%) 9 (12.86%) |
td | Positive |
td | 29 (43.94%) |
td | 9 (12.86%) |
tr | Pleural thickening 0.030a* |
td | Pleural thickening |
td | 0.030a* |
tr | Negative 46 (69.70%) 36 (51.43%) |
td | Negative |
td | 46 (69.70%) |
td | 36 (51.43%) |
tr | Positive 20 (30.30%) 34 (48.57%) |
td | Positive |
td | 20 (30.30%) |
td | 34 (48.57%) |
tr | Offending vessel augmentation in lesions < 0.001a* |
td | Offending vessel augmentation in lesions |
td | < 0.001a* |
tr | Negative 55 (83.33%) 17 (24.29%) |
td | Negative |
td | 55 (83.33%) |
td | 17 (24.29%) |
tr | Positive 11 (16.67%) 53 (75.71%) |
td | Positive |
td | 11 (16.67%) |
td | 53 (75.71%) |
table-wrap-foot | GGO ground-glass opacities #Results are median with interquartile range in parentheses, and the remainder results are measurements with corresponding ratio in parentheses *Data with statistical significance. pa: chi-square test, pb: Student’s t test |
p | GGO ground-glass opacities |
p | #Results are median with interquartile range in parentheses, and the remainder results are measurements with corresponding ratio in parentheses |
p | *Data with statistical significance. pa: chi-square test, pb: Student’s t test |
table-wrap | Table 2 Clinical features of patients in COVID-19 and non-COVID-19 Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value Sex Male# 43 (65.15%) 41 (58.57%) 0.430a Female# 23 (34.85%) 29 (41.43%) Age (years) 46.73 ± 25.00 42.93 ± 13.32 0.275b Vital signs Systolic blood pressure (mmHg) 126.92 ± 23.07 127.07 ± 15.16 0.965b Diastolic blood pressure (mmHg) 77.74 ± 15.72 80.39 ± 10.51 0.254b Respiration rate (bpm) 25.20 ± 7.29 19.86 ± 1.90 < 0.001b* Heart rate (bpm) 101.59 ± 20.36 86.06 ± 13.34 < 0.001b* Temperature (°C) 37.61 ± 1.06 37.12 ± 0.83 0.003b* Signs Dry cough# 56 (84.85%) 48 (68.57%) 0.025a* Fatigue# 8 (12.12%) 22 (31.43%) 0.007a* Sore throat# 6 (9.09%) 9 (12.86%) 0.483a Stuffy# 4 (6.06%) 2 (2.86%) 0.623a Runny nose# 3 (4.55%) 3 (4.29%) 0.731a White blood cell count (× 109/L) 11.48 ± 5.36 5.27 ± 2.33 < 0.001b* White blood cell count category < 0.001c* Low# 0 (0.00%) 2 (2.86%) Normal# 27 (40.91%) 63 (90.00%) High# 39 (59.09%) 5 (7.14%) Lymphocyte count (× 109/L) 1.57 ± 1.33 1.25 ± 0.68 0.086b Lymphocyte count category < 0.001c* Low# 24 (36.36%) 32 (45.71%) Normal# 35 (53.03%) 37 (52.86%) High# 7 (10.61%) 1 (1.43%) Neutrophil count (× 109/L) 8.97 ± 4.90 3.53 ± 2.17 < 0.001b* Neutrophil count category < 0.001c* Low# 3 (4.55%) 8 (11.43%) Normal# 23 (34.85%) 59 (84.29%) High# 40 (60.61%) 3 (4.29%) C-reactive protein (mg/L) 69.30 ± 65.88 26.37 ± 30.97 < 0.001b* Procalcitonin (ng/mL) 3.36 ± 8.98 0.26 ± 0.84 0.007b* *Data with statistical significance. pa: chi-square test, pb: Student’s t test. pc: Kruskal-Wallis H test #Results are measurements with corresponding ratio in parentheses |
label | Table 2 |
caption | Clinical features of patients in COVID-19 and non-COVID-19 |
p | Clinical features of patients in COVID-19 and non-COVID-19 |
table | Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value Sex Male# 43 (65.15%) 41 (58.57%) 0.430a Female# 23 (34.85%) 29 (41.43%) Age (years) 46.73 ± 25.00 42.93 ± 13.32 0.275b Vital signs Systolic blood pressure (mmHg) 126.92 ± 23.07 127.07 ± 15.16 0.965b Diastolic blood pressure (mmHg) 77.74 ± 15.72 80.39 ± 10.51 0.254b Respiration rate (bpm) 25.20 ± 7.29 19.86 ± 1.90 < 0.001b* Heart rate (bpm) 101.59 ± 20.36 86.06 ± 13.34 < 0.001b* Temperature (°C) 37.61 ± 1.06 37.12 ± 0.83 0.003b* Signs Dry cough# 56 (84.85%) 48 (68.57%) 0.025a* Fatigue# 8 (12.12%) 22 (31.43%) 0.007a* Sore throat# 6 (9.09%) 9 (12.86%) 0.483a Stuffy# 4 (6.06%) 2 (2.86%) 0.623a Runny nose# 3 (4.55%) 3 (4.29%) 0.731a White blood cell count (× 109/L) 11.48 ± 5.36 5.27 ± 2.33 < 0.001b* White blood cell count category < 0.001c* Low# 0 (0.00%) 2 (2.86%) Normal# 27 (40.91%) 63 (90.00%) High# 39 (59.09%) 5 (7.14%) Lymphocyte count (× 109/L) 1.57 ± 1.33 1.25 ± 0.68 0.086b Lymphocyte count category < 0.001c* Low# 24 (36.36%) 32 (45.71%) Normal# 35 (53.03%) 37 (52.86%) High# 7 (10.61%) 1 (1.43%) Neutrophil count (× 109/L) 8.97 ± 4.90 3.53 ± 2.17 < 0.001b* Neutrophil count category < 0.001c* Low# 3 (4.55%) 8 (11.43%) Normal# 23 (34.85%) 59 (84.29%) High# 40 (60.61%) 3 (4.29%) C-reactive protein (mg/L) 69.30 ± 65.88 26.37 ± 30.97 < 0.001b* Procalcitonin (ng/mL) 3.36 ± 8.98 0.26 ± 0.84 0.007b* |
tr | Feature Non-COVID-19 (n = 66) COVID-19 (n = 70) p value |
th | Feature |
th | Non-COVID-19 (n = 66) |
th | COVID-19 (n = 70) |
th | p value |
tr | Sex |
td | Sex |
tr | Male# 43 (65.15%) 41 (58.57%) 0.430a |
td | Male# |
td | 43 (65.15%) |
td | 41 (58.57%) |
td | 0.430a |
tr | Female# 23 (34.85%) 29 (41.43%) |
td | Female# |
td | 23 (34.85%) |
td | 29 (41.43%) |
tr | Age (years) 46.73 ± 25.00 42.93 ± 13.32 0.275b |
td | Age (years) |
td | 46.73 ± 25.00 |
td | 42.93 ± 13.32 |
td | 0.275b |
tr | Vital signs |
td | Vital signs |
tr | Systolic blood pressure (mmHg) 126.92 ± 23.07 127.07 ± 15.16 0.965b |
td | Systolic blood pressure (mmHg) |
td | 126.92 ± 23.07 |
td | 127.07 ± 15.16 |
td | 0.965b |
tr | Diastolic blood pressure (mmHg) 77.74 ± 15.72 80.39 ± 10.51 0.254b |
td | Diastolic blood pressure (mmHg) |
td | 77.74 ± 15.72 |
td | 80.39 ± 10.51 |
td | 0.254b |
tr | Respiration rate (bpm) 25.20 ± 7.29 19.86 ± 1.90 < 0.001b* |
td | Respiration rate (bpm) |
td | 25.20 ± 7.29 |
td | 19.86 ± 1.90 |
td | < 0.001b* |
tr | Heart rate (bpm) 101.59 ± 20.36 86.06 ± 13.34 < 0.001b* |
td | Heart rate (bpm) |
td | 101.59 ± 20.36 |
td | 86.06 ± 13.34 |
td | < 0.001b* |
tr | Temperature (°C) 37.61 ± 1.06 37.12 ± 0.83 0.003b* |
td | Temperature (°C) |
td | 37.61 ± 1.06 |
td | 37.12 ± 0.83 |
td | 0.003b* |
tr | Signs |
td | Signs |
tr | Dry cough# 56 (84.85%) 48 (68.57%) 0.025a* |
td | Dry cough# |
td | 56 (84.85%) |
td | 48 (68.57%) |
td | 0.025a* |
tr | Fatigue# 8 (12.12%) 22 (31.43%) 0.007a* |
td | Fatigue# |
td | 8 (12.12%) |
td | 22 (31.43%) |
td | 0.007a* |
tr | Sore throat# 6 (9.09%) 9 (12.86%) 0.483a |
td | Sore throat# |
td | 6 (9.09%) |
td | 9 (12.86%) |
td | 0.483a |
tr | Stuffy# 4 (6.06%) 2 (2.86%) 0.623a |
td | Stuffy# |
td | 4 (6.06%) |
td | 2 (2.86%) |
td | 0.623a |
tr | Runny nose# 3 (4.55%) 3 (4.29%) 0.731a |
td | Runny nose# |
td | 3 (4.55%) |
td | 3 (4.29%) |
td | 0.731a |
tr | White blood cell count (× 109/L) 11.48 ± 5.36 5.27 ± 2.33 < 0.001b* |
td | White blood cell count (× 109/L) |
td | 11.48 ± 5.36 |
td | 5.27 ± 2.33 |
td | < 0.001b* |
tr | White blood cell count category < 0.001c* |
td | White blood cell count category |
td | < 0.001c* |
tr | Low# 0 (0.00%) 2 (2.86%) |
td | Low# |
td | 0 (0.00%) |
td | 2 (2.86%) |
tr | Normal# 27 (40.91%) 63 (90.00%) |
td | Normal# |
td | 27 (40.91%) |
td | 63 (90.00%) |
tr | High# 39 (59.09%) 5 (7.14%) |
td | High# |
td | 39 (59.09%) |
td | 5 (7.14%) |
tr | Lymphocyte count (× 109/L) 1.57 ± 1.33 1.25 ± 0.68 0.086b |
td | Lymphocyte count (× 109/L) |
td | 1.57 ± 1.33 |
td | 1.25 ± 0.68 |
td | 0.086b |
tr | Lymphocyte count category < 0.001c* |
td | Lymphocyte count category |
td | < 0.001c* |
tr | Low# 24 (36.36%) 32 (45.71%) |
td | Low# |
td | 24 (36.36%) |
td | 32 (45.71%) |
tr | Normal# 35 (53.03%) 37 (52.86%) |
td | Normal# |
td | 35 (53.03%) |
td | 37 (52.86%) |
tr | High# 7 (10.61%) 1 (1.43%) |
td | High# |
td | 7 (10.61%) |
td | 1 (1.43%) |
tr | Neutrophil count (× 109/L) 8.97 ± 4.90 3.53 ± 2.17 < 0.001b* |
td | Neutrophil count (× 109/L) |
td | 8.97 ± 4.90 |
td | 3.53 ± 2.17 |
td | < 0.001b* |
tr | Neutrophil count category < 0.001c* |
td | Neutrophil count category |
td | < 0.001c* |
tr | Low# 3 (4.55%) 8 (11.43%) |
td | Low# |
td | 3 (4.55%) |
td | 8 (11.43%) |
tr | Normal# 23 (34.85%) 59 (84.29%) |
td | Normal# |
td | 23 (34.85%) |
td | 59 (84.29%) |
tr | High# 40 (60.61%) 3 (4.29%) |
td | High# |
td | 40 (60.61%) |
td | 3 (4.29%) |
tr | C-reactive protein (mg/L) 69.30 ± 65.88 26.37 ± 30.97 < 0.001b* |
td | C-reactive protein (mg/L) |
td | 69.30 ± 65.88 |
td | 26.37 ± 30.97 |
td | < 0.001b* |
tr | Procalcitonin (ng/mL) 3.36 ± 8.98 0.26 ± 0.84 0.007b* |
td | Procalcitonin (ng/mL) |
td | 3.36 ± 8.98 |
td | 0.26 ± 0.84 |
td | 0.007b* |
table-wrap-foot | *Data with statistical significance. pa: chi-square test, pb: Student’s t test. pc: Kruskal-Wallis H test #Results are measurements with corresponding ratio in parentheses |
p | *Data with statistical significance. pa: chi-square test, pb: Student’s t test. pc: Kruskal-Wallis H test |
p | #Results are measurements with corresponding ratio in parentheses |
p | For imaging manifestations, 7 patients in the COVID-19 group showed normal chest CT (10%). COVID-19 patients have a greater number of pure GGO and mixed GGO than non-COVID-19 patients (p = 0.018 and p = 0.001, respectively). For pure GGO lesions, the differences are significant both in peripheral (p = 0.032) and in central areas (p = 0.001). However, the number of mixed GGO is mainly distributed at the periphery in COVID-19 patients (p < 0.001), with no statistical difference in the central area. The consolidation lesions without GGO occurred less in COVID-19 patients (p = 0.001). More lesions are between 1 and 3 cm (p = 0.027), and fewer lesions are larger than half of the lung segment (p = 0.017) in COVID-19 patients. Other significant differences between the two groups include the pleural traction sign (p = 0.019), bronchial wall thickening (p < 0.001), interlobular septal thickening (p = 0.009), crazy paving (p < 0.001), tree-in-bud (p < 0.001), pleural effusions (p < 0.001), pleural thickening (p = 0.030), and the offending vessel augmentation in lesions (p < 0.001). The lung score presents no significant difference between the COVID-19 and non-COVID-19 groups. |
p | Comparison of clinical features between the two groups of patients with and without COVID-19 is reported in Table 2. There is no significant difference in age and sex between the two groups. Significant differences are found in common symptoms between groups, including fever (p = 0.003), dry cough (p = 0.025), and fatigue (p = 0.007). The respiration rate and heart rate also show significant differences between the two groups (both p < 0.001). Compared with non-COVID-19 pneumonia, the reduction of the WBC count is more pronounced in COVID-19 patients (p < 0.001). The ratio of lymphocyte and ratio of neutrophil also show a significant difference between COVID-19 and non-COVID-19 groups. Although lymphopenia was observed in 32 COVID-19 patients (45.71%), it is not statistically different compared with that in the non-COVID-19 group. C-creative protein (CRP) level and procalcitonin level are also significantly different between the two groups (p < 0.001 and p = 0.007, respectively). Most COVID-19 patients present normal procalcitonin level (82.86%). |
sec | Clinical and radiological feature selection Of the features, 18 radiological features and 17 clinical features were selected to form the predictors based on the result from Tables 1 and 2. Table 3 lists the features selected by univariate analysis and LASSO. Table 3 Selected features in C, R, and CR models Model and individual features Coefficients R, n = 8 (41)* Intercept − 0.307 Total number of mixed GGO in peripheral area 0.359 Total number of consolidation − 1.262 Total number of solid nodules with ground-glass opacities 0.452 Interlobular septal thickening − 5.559 Crazy paving pattern 3.566 Tree-in-bud − 2.548 Pleural thickening 3.265 Offending vessel augmentation in lesions 5.504 C, n = 7 (26)* Intercept 29.273 Respiration − 0.359 Heart rate − 0.054 Temperature − 0.289 White blood cell count − 0.175 Cough − 1.866 Fatigue 2.855 Lymphocyte count category − 0.028 CR, n = 10 (67)* Intercept 45.117 Total number of mixed GGO in peripheral area 0.108 Tree-in-bud − 1.853 Offending vessel augmentation in lesions 6.000 Respiration − 0.583 Heart ratio − 0.084 Temperature − 0.536 White blood cell count − 0.471 Cough − 0.997 Fatigue − 0.228 Lymphocyte count category − 2.177 C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively *n means corresponding selected features, and data in parentheses are total features. Coefficients: the estimate value of each feature in multivariate logistic regression model by “glm” package in R |
title | Clinical and radiological feature selection |
p | Of the features, 18 radiological features and 17 clinical features were selected to form the predictors based on the result from Tables 1 and 2. Table 3 lists the features selected by univariate analysis and LASSO. Table 3 Selected features in C, R, and CR models Model and individual features Coefficients R, n = 8 (41)* Intercept − 0.307 Total number of mixed GGO in peripheral area 0.359 Total number of consolidation − 1.262 Total number of solid nodules with ground-glass opacities 0.452 Interlobular septal thickening − 5.559 Crazy paving pattern 3.566 Tree-in-bud − 2.548 Pleural thickening 3.265 Offending vessel augmentation in lesions 5.504 C, n = 7 (26)* Intercept 29.273 Respiration − 0.359 Heart rate − 0.054 Temperature − 0.289 White blood cell count − 0.175 Cough − 1.866 Fatigue 2.855 Lymphocyte count category − 0.028 CR, n = 10 (67)* Intercept 45.117 Total number of mixed GGO in peripheral area 0.108 Tree-in-bud − 1.853 Offending vessel augmentation in lesions 6.000 Respiration − 0.583 Heart ratio − 0.084 Temperature − 0.536 White blood cell count − 0.471 Cough − 0.997 Fatigue − 0.228 Lymphocyte count category − 2.177 C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively *n means corresponding selected features, and data in parentheses are total features. Coefficients: the estimate value of each feature in multivariate logistic regression model by “glm” package in R |
table-wrap | Table 3 Selected features in C, R, and CR models Model and individual features Coefficients R, n = 8 (41)* Intercept − 0.307 Total number of mixed GGO in peripheral area 0.359 Total number of consolidation − 1.262 Total number of solid nodules with ground-glass opacities 0.452 Interlobular septal thickening − 5.559 Crazy paving pattern 3.566 Tree-in-bud − 2.548 Pleural thickening 3.265 Offending vessel augmentation in lesions 5.504 C, n = 7 (26)* Intercept 29.273 Respiration − 0.359 Heart rate − 0.054 Temperature − 0.289 White blood cell count − 0.175 Cough − 1.866 Fatigue 2.855 Lymphocyte count category − 0.028 CR, n = 10 (67)* Intercept 45.117 Total number of mixed GGO in peripheral area 0.108 Tree-in-bud − 1.853 Offending vessel augmentation in lesions 6.000 Respiration − 0.583 Heart ratio − 0.084 Temperature − 0.536 White blood cell count − 0.471 Cough − 0.997 Fatigue − 0.228 Lymphocyte count category − 2.177 C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively *n means corresponding selected features, and data in parentheses are total features. Coefficients: the estimate value of each feature in multivariate logistic regression model by “glm” package in R |
label | Table 3 |
caption | Selected features in C, R, and CR models |
p | Selected features in C, R, and CR models |
table | Model and individual features Coefficients R, n = 8 (41)* Intercept − 0.307 Total number of mixed GGO in peripheral area 0.359 Total number of consolidation − 1.262 Total number of solid nodules with ground-glass opacities 0.452 Interlobular septal thickening − 5.559 Crazy paving pattern 3.566 Tree-in-bud − 2.548 Pleural thickening 3.265 Offending vessel augmentation in lesions 5.504 C, n = 7 (26)* Intercept 29.273 Respiration − 0.359 Heart rate − 0.054 Temperature − 0.289 White blood cell count − 0.175 Cough − 1.866 Fatigue 2.855 Lymphocyte count category − 0.028 CR, n = 10 (67)* Intercept 45.117 Total number of mixed GGO in peripheral area 0.108 Tree-in-bud − 1.853 Offending vessel augmentation in lesions 6.000 Respiration − 0.583 Heart ratio − 0.084 Temperature − 0.536 White blood cell count − 0.471 Cough − 0.997 Fatigue − 0.228 Lymphocyte count category − 2.177 |
tr | Model and individual features Coefficients |
th | Model and individual features |
th | Coefficients |
tr | R, n = 8 (41)* |
td | R, n = 8 (41)* |
tr | Intercept − 0.307 |
td | Intercept |
td | − 0.307 |
tr | Total number of mixed GGO in peripheral area 0.359 |
td | Total number of mixed GGO in peripheral area |
td | 0.359 |
tr | Total number of consolidation − 1.262 |
td | Total number of consolidation |
td | − 1.262 |
tr | Total number of solid nodules with ground-glass opacities 0.452 |
td | Total number of solid nodules with ground-glass opacities |
td | 0.452 |
tr | Interlobular septal thickening − 5.559 |
td | Interlobular septal thickening |
td | − 5.559 |
tr | Crazy paving pattern 3.566 |
td | Crazy paving pattern |
td | 3.566 |
tr | Tree-in-bud − 2.548 |
td | Tree-in-bud |
td | − 2.548 |
tr | Pleural thickening 3.265 |
td | Pleural thickening |
td | 3.265 |
tr | Offending vessel augmentation in lesions 5.504 |
td | Offending vessel augmentation in lesions |
td | 5.504 |
tr | C, n = 7 (26)* |
td | C, n = 7 (26)* |
tr | Intercept 29.273 |
td | Intercept |
td | 29.273 |
tr | Respiration − 0.359 |
td | Respiration |
td | − 0.359 |
tr | Heart rate − 0.054 |
td | Heart rate |
td | − 0.054 |
tr | Temperature − 0.289 |
td | Temperature |
td | − 0.289 |
tr | White blood cell count − 0.175 |
td | White blood cell count |
td | − 0.175 |
tr | Cough − 1.866 |
td | Cough |
td | − 1.866 |
tr | Fatigue 2.855 |
td | Fatigue |
td | 2.855 |
tr | Lymphocyte count category − 0.028 |
td | Lymphocyte count category |
td | − 0.028 |
tr | CR, n = 10 (67)* |
td | CR, n = 10 (67)* |
tr | Intercept 45.117 |
td | Intercept |
td | 45.117 |
tr | Total number of mixed GGO in peripheral area 0.108 |
td | Total number of mixed GGO in peripheral area |
td | 0.108 |
tr | Tree-in-bud − 1.853 |
td | Tree-in-bud |
td | − 1.853 |
tr | Offending vessel augmentation in lesions 6.000 |
td | Offending vessel augmentation in lesions |
td | 6.000 |
tr | Respiration − 0.583 |
td | Respiration |
td | − 0.583 |
tr | Heart ratio − 0.084 |
td | Heart ratio |
td | − 0.084 |
tr | Temperature − 0.536 |
td | Temperature |
td | − 0.536 |
tr | White blood cell count − 0.471 |
td | White blood cell count |
td | − 0.471 |
tr | Cough − 0.997 |
td | Cough |
td | − 0.997 |
tr | Fatigue − 0.228 |
td | Fatigue |
td | − 0.228 |
tr | Lymphocyte count category − 2.177 |
td | Lymphocyte count category |
td | − 2.177 |
table-wrap-foot | C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively *n means corresponding selected features, and data in parentheses are total features. Coefficients: the estimate value of each feature in multivariate logistic regression model by “glm” package in R |
p | C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively |
p | *n means corresponding selected features, and data in parentheses are total features. Coefficients: the estimate value of each feature in multivariate logistic regression model by “glm” package in R |
sec | Model development and validation The prediction models based on (i) clinical features (C model), (ii) radiological features (R model), and (iii) the combination of clinical features and radiological features (CR model) were developed. ROC analyses for the primary and validation cohort are shown in Table 4 and Fig. 3. The CR model yielded a maximum AUC of 0.986 (95% CI 0.966~1.000) in the primary cohort with the highest accuracy and specificity, which was 0.936 (95% CI 0.866~1.000) in the validation cohort. The AUC for the C model was 0.952 (95% CI 0.988~0.915) and 0.967 (95% CI 0.919~1.000) in the primary and validation cohorts, respectively. For the R model, the AUC of the two cohorts was 0.969 (95% CI 0.940~0.997) and 0.809 (95% CI 0.669~0.948), respectively. Table 4 Performance of the individualized prediction models Primary cohort (n = 98) Validation cohort (n = 38) Models AUC 95% CI Accuracy Specificity Sensitivity AUC 95% CI Accuracy Specificity Sensitivity C model 0.952 0.915~0.988 0.888 0.894 0.882 0.967 0.919~1.000 0.868 0.859 0.842 R model 0.969 0.940~0.997 0.929 0.851 1.000 0.809 0.669~0.948 0.684 0.368 1.000 CR model 0.986 0.966~1.000 0.959 0.957 0.961 0.936 0.866~1.000 0.763 0.789 0.737 C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively. CI confidence interval Fig. 3 ROC of the three models in primary and validation cohort curves. Comparison of receiver operating characteristic (ROC) curves among the radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model) for the diagnosis of COVID-19 in the primary (a) and validation (b) cohorts To determine the clinical usefulness of the diagnostic model, we developed the decision curve (Fig. 4), which showed better performances for the CR model compared with that for the C model and the R model. Across the majority of the range of reasonable threshold probabilities, the decision curve analysis showed that the CR model had a higher overall benefit than the C model and R model. Fig. 4 Decision curve analysis for each model in the primary dataset. The y-axis measures the net benefit, which is calculated by summing the benefits (true-positive findings) and subtracting the harms (false-positive findings), weighting the latter by a factor related to the relative harm of undetected metastasis compared with the harm of unnecessary treatment. The decision curve shows that if the threshold probability is over 10%, the application of the combination of clinical and radiological model (CR model) to diagnose COVID-19 adds more benefit than the clinical model (C model) and radiological model (R model) The nomogram (Fig. 5) was developed by the CR model in the primary cohort, with the factors of the total number of mixed GGO in peripheral area (TN_Mixed_GGO_IP), tree-in-bud, offending vessel augmentation in lesions (OVAIL), respiration, heart ratio, temperature, white blood cell count, cough, fatigue and lymphocyte count category incorporated. The total points were calculated by summing the points identified on the “points” scale for each factor. By comparing the “total points” scale and the “probability” scale, the individual probability of COVID-19 infection could be obtained. Fig. 5 Nomogram of the CR model in the primary cohort. TN_Mixed_GGO_IP represented the total number of mixed GGO in peripheral area. AVAIL represented offending vessel segmentation in lesions. N was a negative result, and P was a positive result. Norm represented normal. Note that in probability scale, 0 = non-COVID-19, 1 = COVID-19 |
title | Model development and validation |
p | The prediction models based on (i) clinical features (C model), (ii) radiological features (R model), and (iii) the combination of clinical features and radiological features (CR model) were developed. ROC analyses for the primary and validation cohort are shown in Table 4 and Fig. 3. The CR model yielded a maximum AUC of 0.986 (95% CI 0.966~1.000) in the primary cohort with the highest accuracy and specificity, which was 0.936 (95% CI 0.866~1.000) in the validation cohort. The AUC for the C model was 0.952 (95% CI 0.988~0.915) and 0.967 (95% CI 0.919~1.000) in the primary and validation cohorts, respectively. For the R model, the AUC of the two cohorts was 0.969 (95% CI 0.940~0.997) and 0.809 (95% CI 0.669~0.948), respectively. Table 4 Performance of the individualized prediction models Primary cohort (n = 98) Validation cohort (n = 38) Models AUC 95% CI Accuracy Specificity Sensitivity AUC 95% CI Accuracy Specificity Sensitivity C model 0.952 0.915~0.988 0.888 0.894 0.882 0.967 0.919~1.000 0.868 0.859 0.842 R model 0.969 0.940~0.997 0.929 0.851 1.000 0.809 0.669~0.948 0.684 0.368 1.000 CR model 0.986 0.966~1.000 0.959 0.957 0.961 0.936 0.866~1.000 0.763 0.789 0.737 C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively. CI confidence interval Fig. 3 ROC of the three models in primary and validation cohort curves. Comparison of receiver operating characteristic (ROC) curves among the radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model) for the diagnosis of COVID-19 in the primary (a) and validation (b) cohorts |
table-wrap | Table 4 Performance of the individualized prediction models Primary cohort (n = 98) Validation cohort (n = 38) Models AUC 95% CI Accuracy Specificity Sensitivity AUC 95% CI Accuracy Specificity Sensitivity C model 0.952 0.915~0.988 0.888 0.894 0.882 0.967 0.919~1.000 0.868 0.859 0.842 R model 0.969 0.940~0.997 0.929 0.851 1.000 0.809 0.669~0.948 0.684 0.368 1.000 CR model 0.986 0.966~1.000 0.959 0.957 0.961 0.936 0.866~1.000 0.763 0.789 0.737 C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively. CI confidence interval |
label | Table 4 |
caption | Performance of the individualized prediction models |
p | Performance of the individualized prediction models |
table | Primary cohort (n = 98) Validation cohort (n = 38) Models AUC 95% CI Accuracy Specificity Sensitivity AUC 95% CI Accuracy Specificity Sensitivity C model 0.952 0.915~0.988 0.888 0.894 0.882 0.967 0.919~1.000 0.868 0.859 0.842 R model 0.969 0.940~0.997 0.929 0.851 1.000 0.809 0.669~0.948 0.684 0.368 1.000 CR model 0.986 0.966~1.000 0.959 0.957 0.961 0.936 0.866~1.000 0.763 0.789 0.737 |
tr | Primary cohort (n = 98) Validation cohort (n = 38) |
th | Primary cohort (n = 98) |
th | Validation cohort (n = 38) |
tr | Models AUC 95% CI Accuracy Specificity Sensitivity AUC 95% CI Accuracy Specificity Sensitivity |
th | Models |
th | AUC |
th | 95% CI |
th | Accuracy |
th | Specificity |
th | Sensitivity |
th | AUC |
th | 95% CI |
th | Accuracy |
th | Specificity |
th | Sensitivity |
tr | C model 0.952 0.915~0.988 0.888 0.894 0.882 0.967 0.919~1.000 0.868 0.859 0.842 |
td | C model |
td | 0.952 |
td | 0.915~0.988 |
td | 0.888 |
td | 0.894 |
td | 0.882 |
td | 0.967 |
td | 0.919~1.000 |
td | 0.868 |
td | 0.859 |
td | 0.842 |
tr | R model 0.969 0.940~0.997 0.929 0.851 1.000 0.809 0.669~0.948 0.684 0.368 1.000 |
td | R model |
td | 0.969 |
td | 0.940~0.997 |
td | 0.929 |
td | 0.851 |
td | 1.000 |
td | 0.809 |
td | 0.669~0.948 |
td | 0.684 |
td | 0.368 |
td | 1.000 |
tr | CR model 0.986 0.966~1.000 0.959 0.957 0.961 0.936 0.866~1.000 0.763 0.789 0.737 |
td | CR model |
td | 0.986 |
td | 0.966~1.000 |
td | 0.959 |
td | 0.957 |
td | 0.961 |
td | 0.936 |
td | 0.866~1.000 |
td | 0.763 |
td | 0.789 |
td | 0.737 |
table-wrap-foot | C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively. CI confidence interval |
p | C, R, and CR indicate the predicted model based on clinical features, radiological features, and the combination of clinical features and clinical radiological features, respectively. CI confidence interval |
figure | Fig. 3 ROC of the three models in primary and validation cohort curves. Comparison of receiver operating characteristic (ROC) curves among the radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model) for the diagnosis of COVID-19 in the primary (a) and validation (b) cohorts |
label | Fig. 3 |
caption | ROC of the three models in primary and validation cohort curves. Comparison of receiver operating characteristic (ROC) curves among the radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model) for the diagnosis of COVID-19 in the primary (a) and validation (b) cohorts |
p | ROC of the three models in primary and validation cohort curves. Comparison of receiver operating characteristic (ROC) curves among the radiological mode (R model), clinical model (C model), and the combination of clinical and radiological model (CR model) for the diagnosis of COVID-19 in the primary (a) and validation (b) cohorts |
p | To determine the clinical usefulness of the diagnostic model, we developed the decision curve (Fig. 4), which showed better performances for the CR model compared with that for the C model and the R model. Across the majority of the range of reasonable threshold probabilities, the decision curve analysis showed that the CR model had a higher overall benefit than the C model and R model. Fig. 4 Decision curve analysis for each model in the primary dataset. The y-axis measures the net benefit, which is calculated by summing the benefits (true-positive findings) and subtracting the harms (false-positive findings), weighting the latter by a factor related to the relative harm of undetected metastasis compared with the harm of unnecessary treatment. The decision curve shows that if the threshold probability is over 10%, the application of the combination of clinical and radiological model (CR model) to diagnose COVID-19 adds more benefit than the clinical model (C model) and radiological model (R model) |
figure | Fig. 4 Decision curve analysis for each model in the primary dataset. The y-axis measures the net benefit, which is calculated by summing the benefits (true-positive findings) and subtracting the harms (false-positive findings), weighting the latter by a factor related to the relative harm of undetected metastasis compared with the harm of unnecessary treatment. The decision curve shows that if the threshold probability is over 10%, the application of the combination of clinical and radiological model (CR model) to diagnose COVID-19 adds more benefit than the clinical model (C model) and radiological model (R model) |
label | Fig. 4 |
caption | Decision curve analysis for each model in the primary dataset. The y-axis measures the net benefit, which is calculated by summing the benefits (true-positive findings) and subtracting the harms (false-positive findings), weighting the latter by a factor related to the relative harm of undetected metastasis compared with the harm of unnecessary treatment. The decision curve shows that if the threshold probability is over 10%, the application of the combination of clinical and radiological model (CR model) to diagnose COVID-19 adds more benefit than the clinical model (C model) and radiological model (R model) |
p | Decision curve analysis for each model in the primary dataset. The y-axis measures the net benefit, which is calculated by summing the benefits (true-positive findings) and subtracting the harms (false-positive findings), weighting the latter by a factor related to the relative harm of undetected metastasis compared with the harm of unnecessary treatment. The decision curve shows that if the threshold probability is over 10%, the application of the combination of clinical and radiological model (CR model) to diagnose COVID-19 adds more benefit than the clinical model (C model) and radiological model (R model) |
p | The nomogram (Fig. 5) was developed by the CR model in the primary cohort, with the factors of the total number of mixed GGO in peripheral area (TN_Mixed_GGO_IP), tree-in-bud, offending vessel augmentation in lesions (OVAIL), respiration, heart ratio, temperature, white blood cell count, cough, fatigue and lymphocyte count category incorporated. The total points were calculated by summing the points identified on the “points” scale for each factor. By comparing the “total points” scale and the “probability” scale, the individual probability of COVID-19 infection could be obtained. Fig. 5 Nomogram of the CR model in the primary cohort. TN_Mixed_GGO_IP represented the total number of mixed GGO in peripheral area. AVAIL represented offending vessel segmentation in lesions. N was a negative result, and P was a positive result. Norm represented normal. Note that in probability scale, 0 = non-COVID-19, 1 = COVID-19 |
figure | Fig. 5 Nomogram of the CR model in the primary cohort. TN_Mixed_GGO_IP represented the total number of mixed GGO in peripheral area. AVAIL represented offending vessel segmentation in lesions. N was a negative result, and P was a positive result. Norm represented normal. Note that in probability scale, 0 = non-COVID-19, 1 = COVID-19 |
label | Fig. 5 |
caption | Nomogram of the CR model in the primary cohort. TN_Mixed_GGO_IP represented the total number of mixed GGO in peripheral area. AVAIL represented offending vessel segmentation in lesions. N was a negative result, and P was a positive result. Norm represented normal. Note that in probability scale, 0 = non-COVID-19, 1 = COVID-19 |
p | Nomogram of the CR model in the primary cohort. TN_Mixed_GGO_IP represented the total number of mixed GGO in peripheral area. AVAIL represented offending vessel segmentation in lesions. N was a negative result, and P was a positive result. Norm represented normal. Note that in probability scale, 0 = non-COVID-19, 1 = COVID-19 |
sec | Discussion In this multi-center study, statistical analysis was performed in comparing imaging and clinical manifestations between pneumonia patients with and without COVID-19. Eighteen radiological semantic features and seventeen clinical features were identified to be significantly different between the two groups (p < 0.05). Three models for COVID-19 diagnosis were developed based on the refined features. The models were validated in the both primary and validation cohorts and achieved an AUC as high as 0.986. These models will play an essential role for early and easy-to-access diagnosis, especially when there are not enough RT-PCT kits or experimental platforms to test for the COVID-19 infection. A total of 1745 lesions were evaluated for the qualitative feature, location, and size in this study. Consistent with the previous studies, the ground-glass opacities and consolidation in the lung periphery were considered to be the imaging hallmark in patients with COVID-19 infection [11, 25]. However, when we subdivided the GGO into pure GGO and mixed GGO, we found that the distribution pattern is different between these two lesions. Pure GGO show differences between groups in every location of the lungs, whereas mixed GGO only have significant differences between groups in the lung periphery. Recent studies defined four stages of lung involvement in COVID-19 [26]. Therefore, a follow-up analysis of these distributions would be significant. The lesion size in patients with COVID-19 infection was another interesting observation. Most lesions were between 1 and 3 cm, with few lesions larger than half of the lung segment, which was similar to the finding in MERS_CoV [22]. Other features similar to MERS_CoV and SARS_CoV were observed in the laboratory abnormalities, such as lymphopenia, which may be associated with the cellular immune deficiency [3, 27]. However, our results showed no significant difference in lymphopenia between the COVID-19 and non-COVID-19 patients. To our knowledge, no diagnostic model based on imaging and clinical features alone has been proposed for the diagnosis of COVID-19. Our clinical and radiological semantic (CR) models consisted of the following features: total number of GGO with consolidation in the peripheral area, tree-in-bud, offending vessel augmentation in lesions, temperature, heart ratio, respiration, cough and fatigue, WBC count, and lymphocyte count category. The CR model outperformed the individual clinical and radiologic model. This result was in accordance with that in previous study in breast cancer, in which the model based on the combination of radiomics features and clinical features achieved a higher performance [24]. Compared with the radiomics-based model, the extraction of radiological semantic features can overcome the image discrepancy caused by different scanning parameters and/or different CT vendors. A previous study [28] also indicated that models based on semantic features determined by an experienced thoracic radiologist slightly outperformed models based on computed texture features alone. There are a few limitations in this study. First, the sample size is relatively small because this is a retrospective analysis of a new disease and most of the cases outside of Wuhan City are imported. Second, with the multi-center retrospective design, there is a potential bias of patient selection [29], since there may be some deviations in marking semantic features among readers, though we have taken the effort to reduce this by creating pictorial examples and setting feature criteria (Supplementary Materials). Third, longitudinal CT study was not performed. Whether or not this model can be used to evaluate the follow-ups and help to guide therapy remains an open question to be further explored. Moreover, the rich high-order features of the CT image combined with radiomics or deep learning have not been studied, which may be another way to identify the patients with COVID-19. Besides, one can also focus on the role of radiological features in disease monitoring, treatment evaluation, and prognosis prediction. In conclusion, 1745 lesions and 67 features were compared between pneumonia patients with and without COVID-19. Thirty-five features were significantly different between the two groups. A diagnostic model with AUC as high as 0.986 was developed and validated both in the primary and in the validation cohorts, which may help improve the COVID-19 diagnosis. |
title | Discussion |
p | In this multi-center study, statistical analysis was performed in comparing imaging and clinical manifestations between pneumonia patients with and without COVID-19. Eighteen radiological semantic features and seventeen clinical features were identified to be significantly different between the two groups (p < 0.05). Three models for COVID-19 diagnosis were developed based on the refined features. The models were validated in the both primary and validation cohorts and achieved an AUC as high as 0.986. These models will play an essential role for early and easy-to-access diagnosis, especially when there are not enough RT-PCT kits or experimental platforms to test for the COVID-19 infection. |
p | A total of 1745 lesions were evaluated for the qualitative feature, location, and size in this study. Consistent with the previous studies, the ground-glass opacities and consolidation in the lung periphery were considered to be the imaging hallmark in patients with COVID-19 infection [11, 25]. However, when we subdivided the GGO into pure GGO and mixed GGO, we found that the distribution pattern is different between these two lesions. Pure GGO show differences between groups in every location of the lungs, whereas mixed GGO only have significant differences between groups in the lung periphery. Recent studies defined four stages of lung involvement in COVID-19 [26]. Therefore, a follow-up analysis of these distributions would be significant. The lesion size in patients with COVID-19 infection was another interesting observation. Most lesions were between 1 and 3 cm, with few lesions larger than half of the lung segment, which was similar to the finding in MERS_CoV [22]. Other features similar to MERS_CoV and SARS_CoV were observed in the laboratory abnormalities, such as lymphopenia, which may be associated with the cellular immune deficiency [3, 27]. However, our results showed no significant difference in lymphopenia between the COVID-19 and non-COVID-19 patients. |
p | To our knowledge, no diagnostic model based on imaging and clinical features alone has been proposed for the diagnosis of COVID-19. Our clinical and radiological semantic (CR) models consisted of the following features: total number of GGO with consolidation in the peripheral area, tree-in-bud, offending vessel augmentation in lesions, temperature, heart ratio, respiration, cough and fatigue, WBC count, and lymphocyte count category. The CR model outperformed the individual clinical and radiologic model. This result was in accordance with that in previous study in breast cancer, in which the model based on the combination of radiomics features and clinical features achieved a higher performance [24]. Compared with the radiomics-based model, the extraction of radiological semantic features can overcome the image discrepancy caused by different scanning parameters and/or different CT vendors. A previous study [28] also indicated that models based on semantic features determined by an experienced thoracic radiologist slightly outperformed models based on computed texture features alone. |
p | There are a few limitations in this study. First, the sample size is relatively small because this is a retrospective analysis of a new disease and most of the cases outside of Wuhan City are imported. Second, with the multi-center retrospective design, there is a potential bias of patient selection [29], since there may be some deviations in marking semantic features among readers, though we have taken the effort to reduce this by creating pictorial examples and setting feature criteria (Supplementary Materials). Third, longitudinal CT study was not performed. Whether or not this model can be used to evaluate the follow-ups and help to guide therapy remains an open question to be further explored. Moreover, the rich high-order features of the CT image combined with radiomics or deep learning have not been studied, which may be another way to identify the patients with COVID-19. Besides, one can also focus on the role of radiological features in disease monitoring, treatment evaluation, and prognosis prediction. |
p | In conclusion, 1745 lesions and 67 features were compared between pneumonia patients with and without COVID-19. Thirty-five features were significantly different between the two groups. A diagnostic model with AUC as high as 0.986 was developed and validated both in the primary and in the validation cohorts, which may help improve the COVID-19 diagnosis. |
sec | Electronic supplementary material ESM 1 (DOCX 404 kb) |
title | Electronic supplementary material |
sec | ESM 1 (DOCX 404 kb) |
p | ESM 1 (DOCX 404 kb) |
label | ESM 1 |
caption | (DOCX 404 kb) |
p | (DOCX 404 kb) |
back | Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Xiaofeng Chen and Yanyan Tang contributed equally to this work. Funding information This study has received funding by the Natural Science Foundation of China (grant numbers 81471730, 31870981) to R.W.; the Natural Science Foundation of Guangdong Province (grant number 2018A030307057) to Z.D.; and the Special Project on Prevention and Control of COVID-19 for Colleges and Universities in Guangdong Province (grant number 2020KZDZX1085) to Z.D. Compliance with ethical standards Guarantor The scientific guarantor of this publication is Zhuozhi Dai. Conflict of interest One of the authors of this manuscript (Yuting Liao) is an employee of GE Healthcare. The remaining authors declare no relationships with any companies whose products or services may be related to the subject matter of the article. Statistics and biometry One of the authors, Dr. Yuting Liao, has significant statistical expertise. Informed consent Written informed consent was waived by the Institutional Review Board. Ethical approval Institutional Review Board approval was obtained. Methodology • retrospective • case-control study • multi-center study |
footnote | Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
p | Publisher’s note |
p | Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
footnote | Xiaofeng Chen and Yanyan Tang contributed equally to this work. |
p | Xiaofeng Chen and Yanyan Tang contributed equally to this work. |
notes | Funding information This study has received funding by the Natural Science Foundation of China (grant numbers 81471730, 31870981) to R.W.; the Natural Science Foundation of Guangdong Province (grant number 2018A030307057) to Z.D.; and the Special Project on Prevention and Control of COVID-19 for Colleges and Universities in Guangdong Province (grant number 2020KZDZX1085) to Z.D. |
title | Funding information |
p | This study has received funding by the Natural Science Foundation of China (grant numbers 81471730, 31870981) to R.W.; the Natural Science Foundation of Guangdong Province (grant number 2018A030307057) to Z.D.; and the Special Project on Prevention and Control of COVID-19 for Colleges and Universities in Guangdong Province (grant number 2020KZDZX1085) to Z.D. |
notes | Compliance with ethical standards Guarantor The scientific guarantor of this publication is Zhuozhi Dai. Conflict of interest One of the authors of this manuscript (Yuting Liao) is an employee of GE Healthcare. The remaining authors declare no relationships with any companies whose products or services may be related to the subject matter of the article. Statistics and biometry One of the authors, Dr. Yuting Liao, has significant statistical expertise. Informed consent Written informed consent was waived by the Institutional Review Board. Ethical approval Institutional Review Board approval was obtained. Methodology • retrospective • case-control study • multi-center study |
title | Compliance with ethical standards |
notes | Guarantor The scientific guarantor of this publication is Zhuozhi Dai. |
title | Guarantor |
p | The scientific guarantor of this publication is Zhuozhi Dai. |
notes | Conflict of interest One of the authors of this manuscript (Yuting Liao) is an employee of GE Healthcare. The remaining authors declare no relationships with any companies whose products or services may be related to the subject matter of the article. |
title | Conflict of interest |
p | One of the authors of this manuscript (Yuting Liao) is an employee of GE Healthcare. The remaining authors declare no relationships with any companies whose products or services may be related to the subject matter of the article. |
notes | Statistics and biometry One of the authors, Dr. Yuting Liao, has significant statistical expertise. |
title | Statistics and biometry |
p | One of the authors, Dr. Yuting Liao, has significant statistical expertise. |
notes | Informed consent Written informed consent was waived by the Institutional Review Board. |
title | Informed consent |
p | Written informed consent was waived by the Institutional Review Board. |
notes | Ethical approval Institutional Review Board approval was obtained. |
title | Ethical approval |
p | Institutional Review Board approval was obtained. |
notes | Methodology • retrospective • case-control study • multi-center study |
title | Methodology |
p | • retrospective |
p | • case-control study |
p | • multi-center study |
Annnotations TAB TSV DIC JSON TextAE
last updated at 2021-06-18 19:24:49 UTC
- Denotations: 1
- Blocks: 0
- Relations: 0