Fmed 08 651556
Fmed 08 651556
Fmed 08 651556
yielded an area under curve (AUC) of 0.963 and 0.943, with sensitivity of 89.7 and 86.2%,
specificity of 89.7 and 89.7%, and accuracy of 89.7 and 87.1%, respectively.
Conclusions: Combination of distribution pattern and category of pulmonary opacity
on chest CT with clinical features facilitates the differentiation of COVID-19 pneumonia
from H1N1 pneumonia.
Keywords: coronavirus disease 2019, influenza A (H1N1), computed tomography, multivariate analysis, differential
diagnosis
INTRODUCTION features of the two diseases and further establish a model based
on clinicoradiologic features so as to provide some guidance for
Coronavirus disease 2019 (COVID-19), caused by the novel their early identification.
severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2,
previously known as 2019-nCoV), has become a global health
concern that threaten human life and public health security.
METHODS
As of 3 January, 2021, more than 83 million cases and more Study Population
than 1.8 million deaths have been reported worldwide according A search of the medical records in two hospitals’ information
to World Health Organization (WHO) statistics (1). When system was conducted, and 303 COVID-19 patients from January
assessing COVID-19, it is noteworthy that influenza viruses 20 to February 13, 2020, and 224 patients with H1N1 pneumonia
occur in the same season. Influenza A (H1N1) is the most from May 24, 2009, to January 29, 2010, were identified. All
common influenza, and it caused a worldwide pandemic in 2009– patients were diagnosed according to the diagnostic criteria of
2010 with more than 18,449 deaths (2) and has now become the National Health Commission of China (22, 23) and were
an annual seasonal influenza, leading to a large number of confirmed by real-time Polymerase Chain Reaction (RT-PCR)
hospitalizations and deaths (3, 4). Because of the differences in detection of virus nuclear acid. Our exclusion criteria were (a)
therapy, prognosis, and protective measure between COVID- patients without CT scans; (b) patients with normal CT imaging;
19 and H1N1, it is important for clinicians and radiologists to and (c) COVID-19 patients with positive influenza A nuclear acid
identify these two respiroviral infections. tests. Finally, 291 patients with COVID-19 pneumonia and 97
Both COVID-19 and H1N1 pneumonias share similar clinical patients with H1N1 pneumonia were enrolled. All patients were
manifestations, such as mild to moderate flulike syndromes, and randomly grouped in seven-to-three ratio into a primary cohort
are often cured by symptomatic treatments (5, 6). But a few (n = 272) and a validation cohort (n = 116), respectively, by using
patients develop severe or even lethal acute respiratory distress computer-generated random numbers. The flow chart of patient
syndrome (ARDS), particularly in patients with comorbidities selection, grouping, and disease subtypes is shown in Figure 1.
(7–11), and their laboratory exams often display lymphopenia,
abnormalities in liver function, and myocardial zymogram (7, 8, Ethics Statement
10–14). This work was approved by the Ethics Committee of Shanghai
Significantly, the typical findings on computed tomography Public Health Clinical Center, China, and Wuhan Union
(CT) seem different in both pneumonias: bilateral, multifocal Red Cross Hospital, China, and informed consent for this
ground-glass opacities (GGOs) with or without consolidations retrospective study was waived (YJ-2020-S035-01).
or intralobular lines, in a predominant peripheral distribution
usually present in COVID-19 pneumonia (15–17); nodules and Clinical Data Collection
bud signs (18, 19), bronchiectasis (18, 20), and pleural effusion The clinical data on admission were retrospectively collected.
(20) are common in H1N1 pneumonia but rare in COVID-19 Particular attention was paid to the demographics, comorbidities,
pneumonia (16). Unlike the predominant peripheral distribution coinfection, symptoms, and laboratory findings. The
in COVID-19 pneumonia, H1N1 pneumonia presents as a comorbidities or underlying medical conditions mainly
predominant peribronchovascular (18, 21) or peripheral (18, 21) included chronic pulmonary, cardiac, renal and hepatic diseases,
or mixed distribution (18, 19). The differences in CT findings diabetes, cerebrovascular disease, hyperlipemia, malignancy,
in category and the distribution patterns of pulmonary opacity and immunosuppression.
suggest their significance in terms of a differential diagnosis of
both pneumonias. Imaging Acquiring
Although it is easy to identify these typical lesions, CT All patients underwent CT examinations at full inspiration from
manifestations of COVID-19 and H1N1 pneumonia are very the thoracic inlet to the costophrenic angle level. CT scans
diverse. To date, knowledge regarding the comprehensive were performed with one of two scanners (Hitachi Scenaria 64,
identification of both pneumonias still remains limited Hitachi Medical Systems, Tokyo, Japan; or Siemens Sensation
and cannot meet the urgent clinical needs. Therefore, this 16, Siemens Medical Systems, Forchheim, Germany) using
retrospective study aimed to assess initial CT and clinical automatic exposure control with the following parameters: tube
FIGURE 1 | Flow chart shows patient selection, grouping and disease subtypes. COVID-19, coronavirus disease 2019; RT-PCR, real-time Polymerase Chain Reaction.
voltage, 120 or 140 kV; tube current, 150–250 mA; detector Fleischner Society (24). The images were analyzed independently
width, 64 × 0.625 mm or 16 × 0.75 mm; pitch, 1.57 or 1; rotation by two radiologists (Weiya Shi and Fei Shan, with 12 and 19
time, 0.35 or 0.5 s; field of view (FOV), 350 mm; and matrix, 512 years of experience in chest radiology, respectively). In cases of
× 512. The reconstruction kernel used was lung smooth with disagreement, the results were determined by consensus.
a thickness of 1 mm and an interval of 0.8 mm. The following
windows were used: a mediastinal window with a window width Statistical Analysis
of 350 Hounsfield unit (Hu) and a window level of 40 Hu, and a The continuous data are expressed as the median and
lung window with a width of 1,200 Hu and a level of −600 Hu. interquartile range (IQR, 25th and 75th percentiles), because
a majority of them did not follow a normal distribution.
Image Interpretation The Fisher exact test and the chi-square test were used
CT images were assessed for the presence and distribution of for categorical variables, and the Wilcoxon rank sum test
pulmonary opacities, including pure GGOs, which manifested as was used for continuous variables when comparing the
a hazy opacity without obscuring the underlying vessels; GGO clinicoradiologic features of COVID-19 pneumonia with those of
with interlobular septal thickening or reticulation, which was H1N1 pneumonia.
defined as a crazy-paving sign; GGOs with consolidation, which In order to evaluate the interobserver agreement (IA) between
was defined as an area of opacification obscuring the underlying the two radiologists, the Cohen’s Kappa test was used for
vessels in GGO; consolidation; centrilobular nodule or tree-in- categorical variables and the intraclass correlation coefficient
bud sign, which was regarded present when centrilobular nodules (ICC) for continuous variables. The kappa coefficient (k) between
or nodular branching structures resembled a budding tree. 0.00 and 0.20; 0.21 and 0.40; 0.41 and 0.60; 0.61 and 0.80; and
The distribution pattern of pulmonary opacities was 0.81 and 1.00 indicated slight, fair, moderate, substantial, and
assessed as being in a predominant peripheral (outer third almost perfect agreement, respectively. The ICC values between
of the lungs), peribronchovascular, both of peripheral and 0.00 and 0.25; 0.26 and 0.40; 0.41 and 0.60; 0.61 and 0.75;
peribronchovascular, or diffuse distribution, or lacking a 0.75 and 1.00, indicated poor, low, fair, good, and excellent
specific distribution (Figure 2). The laterality (unilateral and agreement, respectively.
bilateral) and predominant involved pulmonary lobes (upper, The clinicoradiologic characteristics found to be significant
middle/lingula, lower, or diffuse) were also assessed. in univariate analysis were inputted into the least absolute
Bronchial wall thickening or bronchiectasis, focal pulmonary shrinkage and selection operator (LASSO) logistic regression
fibrosis (including reticulation and liner opacity), pleural analysis to identify the optimal subset of clinicoradiologic
effusion, and mediastinal lymphadenopathy (>1 cm in short- features in order to develop a model for identification. Receiver
axis diameter) were noted. The number of pulmonary segments operating characteristic (ROC) curves were plotted for assessing
involved was counted. All the terms were defined according to the the performance of the model in the primary and validation
FIGURE 2 | The distribution patterns of pneumonia due to COVID-19 or H1N1. (A) A 74-year-old male with COVID-19 pneumonia presents the onset symptom of
fever (37.5◦ C), and the CT shows GGO with consolidation and crazy-paving sign mainly along subpleural lungs, namely, as a peripheral distribution pattern. (B) A
29-year-old male with H1N1 pneumonia presents the onset symptoms of fever (39.5◦ C), cough, and shortness of breath, and the CT shows consolidation and small
centrilobular nodules mainly along bronchovascular bundles, namely, as a peribronchovascular distribution pattern. (C) A 59-year-old female with H1N1 pneumonia
presents the onset symptoms of fever (38.5◦ C), cough, and shortness of breath, and the CT shows consolidation with GGO along bronchovascular bundles and
subpleural lungs, namely, as a distribution pattern of both peripheral and peribronchovascular. (D) A 71-year-old male with COVID-19 pneumonia presents the onset
symptoms of fever (38◦ C) and dyspnea, and CT shows diffuse GGO with consolidation and crazy-paving sign in both lungs, namely, as a diffuse distribution pattern.
(E) A 50-year-old female with COVID-19 pneumonia presents the onset symptoms of fever (39◦ C) and cough, and CT shows very small non-rounded GGO located in
the middle lobe of the right lung lacking a specific distribution. COVID-19, coronavirus disease 2019; GGOs, ground-glass opacities.
TABLE 1 | Baseline clinical and laboratory characteristics of the patients with COVID-19 pneumonia vs. those with H1N1 pneumonia.
Continuous variables are presented as median (interquartile range), and dichotomous variables are presented as numbers of patients, with percentages in parentheses. *P < 0.05;
P-values are from Fisher’s exact test or Wilcoxon rank-sum test when comparing characteristics of COVID-19 pneumonia with those of H1N1 pneumonia; COVID-19, coronavirus
disease 2019; AST, aspartate aminotransferase.
cohorts. The cut-off values were defined based on the maximal 0.001, p < 0.001, p = 0.032, 0.011, 0.001, 0.010, 0.017, and
Youden index. 0.023, respectively).
Statistical analysis was performed with R version 3.6.1 (R
Project for Statistical Computing, Vienna, Austria). A two-tailed
CT Characteristics
α < 0.05 was considered statistically significant.
The IA between the two radiologists was almost perfect for all CT
findings except the tree-in-bud sign, for which it was substantial
RESULTS (k = 0.789).
The baseline CT characteristics of the 388 cases were shown in
Clinical and Laboratory Features Table 2. In terms of lesions’ distribution pattern, 90.0% (262/291)
The baseline clinical and laboratory characteristics of the 388 of COVID-19 pneumonia had a peripheral distribution pattern in
cases were shown in Table 1. Compared with H1N1 pneumonia contrast to only 20.6% (20/97) in H1N1 pneumonia (p < 0.001).
(97 patients), patients with COVID-19 pneumonia (291 patients) In the H1N1 pneumonia, a peribronchovascular distribution
were older (51.0 vs. 31.0 years, p < 0.001) and had a lower pattern (52/97, 53.6%) was the most common, though this is rare
proportion of men (156/291, 53.6% vs. 66/97, 68.0%, p = 0.013). in COVID-19 pneumonia (9/291, 3.1%, p < 0.001). The H1N1
They had lower fever incidence (239/291, 82.1% vs. 93/97, 95.9%, pneumonia was more likely to exist in a distribution pattern
p = 0.001) and lower body temperatures (38.0 vs. 38.8◦ C, p < of both peripheral and peribronchovascular (17/97, 17.5% vs.
0.001). H1N1 patients were more likely to have coinfection than 14/291, 4.8%, p < 0.001) or be lacking a specific distribution
COVID-19 patients (16/97, 16.5% vs. 13/291, 4.5%, p < 0.001). (6/97, 6.2% vs. 3/291, 1.0%, p = 0.003). The incidence of diffuse
Among them, all were bacterial infections except for 2 H1N1 distribution pattern was 1.0% (3/291) and 2.1% (2/97) in COVID-
cases with bacterial and fungal coinfections. 19 pneumonia and H1N1 pneumonia, respectively, and had no
The laboratory exams displayed lower serum levels of significant difference.
alanine aminotransferase, aspartate aminotransferase (AST), With respect to other CT characteristics, the COVID-
lactate dehydrogenase, total bilirubin, and c-reactive protein, 19 pneumonia was more likely to present a crazy-paving
higher level of CD4+ T counts and blood urea nitrogen, as well sign (206/291, 70.8% vs. 39/97, 40.2%, p < 0.001), GGO
as more frequency of increased serum creatinine in COVID- with consolidation (216/291, 74.2% vs. 49/97, 50.5%, p <
19 pneumonia in contrast to those in H1N1 pneumonia (p < 0.001), bilateral involvement (241/291, 82.8% vs. 71/97, 73.2%,
TABLE 2 | Baseline CT characteristics of the patients with COVID-19 pneumonia vs. those with H1N1 pneumonia.
Characteristics COVID-19 (n = 291) H1N1 (n = 97) P-value Agreement (k or ICC; 95% CI)
Continuous variables are presented as median (interquartile range), and dichotomous variables are presented as numbers of patients, with percentages in parentheses. *P < 0.05;
P-values are from Fisher’s exact test or Wilcoxon rank-sum test when comparing characteristics of COVID-19 pneumonia with those of H1N1 pneumonia. For categorical variables
and continuous variables, the interobserver agreement is assessed with the kappa coefficient (k) and the intraclass correlation coefficient (ICC), respectively. The k and ICC values are
reported with 95% confidence intervals (95% CI). COVID-19, coronavirus disease 2019; GGOs, ground-glass opacities.
p = 0.039), and focal pulmonary fibrosis (103/291, 35.4% vs. and a peribronchovascular distribution pattern were inversely
21/97, 21.6%, p = 0.012); H1N1 pneumonia, however, was more associated with COVID-19 pneumonia.
likely to present consolidation (79/97, 81.4% vs. 91/291, 31.3%, For the primary and validation cohorts, the LASSO model
p < 0.001), centrilobular nodule or tree-in-bud sign (43/97, containing above eight features yielded an area under curve
44.3% vs. 3/291, 1.0%, p < 0.001), predominant middle/lingula (AUC) of 0.963 (95% CI: 0.942–0.984) and 0.943 (95% CI: 0.900–
involvement (12/97, 12.4% vs. 12/291, 4.1%, p = 0.003), bronchial 0.986), with sensitivity of 89.7 and 86.2%, specificity of 89.7 and
wall thickening or bronchiectasis (28/97, 28.9% vs. 18/291, 6.2%, 89.7%, accuracy of 89.7 and 87.1%, positive predictive value of
p < 0.001), and pleural effusion (27/97, 27.8% vs. 14/291, 4.8%, 96.3 and 96.1%, and negative predictive value of 74.4 and 68.4%,
p < 0.001). respectively (Table 4 and Figure 4).
FIGURE 3 | The selection of clinico-radiological features using LASSO logistic regression. (A) Optimal feature selection according to AUC value. (B) LASSO coefficient
profiles of the features. Vertical line is drawn at the selected value using 10-fold cross-validation, where optimal λ results in eight non-zero coefficients. LASSO, least
absolute shrinkage and selection operator; AUC, area under curve.
aimed to accurately identify these two diseases. In our study, (AUC, 0.943; sensitivity, 86.2%; specificity, 89.7%; accuracy,
CT manifestations of consolidation, centrilobular nodule or 87.1%) in the validation cohort, which provides guidance for
tree-in-bud sign, bronchial wall thickening or bronchiectasis, clinical diagnosis.
peripheral distribution pattern and peribronchovascular Several previous studies have found that COVID-19 typically
distribution pattern, together with the clinical features such presents GGO with or without consolidation in a peripheral
as age, body temperature and AST were identified the optimal distribution (15–17), which was endorsed by the Society of
features subset for differentiating COVID-19 pneumonia from Thoracic Radiology, the American College of Radiology, and
H1N1 pneumonia. Our model had high diagnostic efficiency RSNA (16) as guidance for COVID-19 diagnosis. Our study was
consistent with previous reports that the COVID-19 pneumonia is, coinfection (mainly bacterial infection) was more frequent in
mainly presented a peripheral distribution pattern (262/291, H1N1 pneumonia (16/97, 16.5%) than in COVID-19 pneumonia
90.0%); in contrast, H1N1 pneumonia most commonly presented (13/291, 4.5%).
a peribronchovascular distribution pattern (52/97, 53.6%). The Besides radiological features, clinical features such as age,
differences in CT imaging between these two pneumonias body temperature, and AST should be also taken into
may result from their distinct pathological changes in lungs. consideration for differentiation. Our cohort showed an older
The pathological findings of COVID-19 pneumonia include median age of the COVID-19 patients; however, it should be
exudative diffuse alveolar damage with alveolar and interstitial viewed cautiously, the statistics were based on the early stage
edema, alveolar fibrinous exudate with hyaline membranes, of the outbreak in Shanghai and Hubei province, China. With
and reactive pneumocytes (26), whereas H1N1 pneumonia, in the global spread, it has been found that the youth are also a
addition to diffuse alveolar damage, is usually accompanied susceptible population (30).
by necrotizing bronchiolitis and alveolar hemorrhage (27). Although the model had high diagnostic efficiency, it should
These pathologic lesion-dependent distribution patterns, which be noted that it had a negative predictive value of 68.4% in
are very conspicuous at the first glance of the images, are the validation cohort due to the misdiagnosis of 12 cases of
valuable indicators for differentiating COVID-19 pneumonia COVID-19 pneumonia as H1N1 pneumonia. The misdiagnosed
from H1N1 pneumonia. patients were relatively young (median, 38.5 years; IQR, 32.5–
The bronchiolitis causes central lobular nodules or tree-bud 58.5), with high fever (median, 39.0◦ C; IQR, 38.3–39.7) and
signs, and bronchial wall thickening or bronchiectasis; therefore, moderately elevated AST (median, 30.0 U/L; IQR, 22.5–48.0),
it is quite understandable that these two signs were more and 75.0% (9/12) of them presented consolidation and 58.3%
common in H1N1 pneumonia than in COVID-19 pneumonia (7/12) a non-peripheral distribution pattern (two cases with
(both p < 0.001), consisting with other studies (16, 18, 21). a peribronchovascular distribution pattern, four cases with a
We also found that consolidation was more frequent in H1N1 distribution pattern of both peripheral and peribronchovascular,
pneumonia (79/97, 81.4%) than in COVID-19 pneumonia and one case with a diffuse distribution pattern). Radiologists
(91/291, 31.3%), which is consistent with previous studies should pay more attention to these atypical clinicoradiologic
and is possibly associated with pathologic basis (8), disease manifestations of COVID-19 in young individuals so as to
progression or more severe disease (28), bacterial coinfection avoid misdiagnosis.
(29). The latter phenomenon was also found in our study, that However, our study had limitations. Firstly, there was an
imbalance between the sample sizes of COVID-19 pneumonia
and H1N1 pneumonia, and the proportion of severe and critically
ill H1N1 patients was greater than those of the COVID-
TABLE 3 | The coefficients of the elected features by LASSO logistic regression 19 cohort, which may have led to statistical disequilibrium.
analysis. Secondly, there is a bias in the laboratory tests because there
were several laboratory changes in COVID-19 patients compared
Feature selected Coefficient
to H1N1 patients, using only the tests common to both groups.
Age, years 0.013
Thirdly, patients comorbid with chronic pulmonary disease
Body temperature, ◦ C −0.523
were not excluded, which may lead to a bias in this study,
AST, U/L −2.998e-03
although the proportion of these patients was very low and
Consolidation (yes vs. no) −0.300
there was no statistically significant difference between both
Centrilobular nodule or tree-in-bud sign (yes vs. no) −0.984
groups. Fourthly, coinfection is common in patients of both
groups, especially in patients with N1H1. Due to the complex
Bronchial wall thickening or bronchiectasis (yes vs. no) −0.207
nature of the clinical situation, we believe that differential
Peripheral distribution pattern (yes vs. no) 1.505
diagnosis is also necessary for these patients to obtain effective
Peribronchovascular distribution pattern (yes vs. no) −0.712
subsequent treatment. Therefore, we did not exclude these
LASSO, least absolute shrinkage and selection operator; AST, aspartate patients when constructing our analysis model, but this may lead
aminotransferase. to a bias.
TABLE 4 | AUC values of the LASSO regression model for differentiating COVID-19 pneumonia from H1N1 pneumonia in the primary and validation cohorts.
LASSO regression model AUC value Sensitivity Specificity Accuracy PPV NPV
(95% CI) (%) (%) (%) (%) (%)
COVID-19, coronavirus disease 2019; AUC, area under curve; LASSO, least absolute shrinkage and selection operator; PPV, positive predictive value; NPV, negative predictive value;
95% CI, 95% confidence intervals.
FIGURE 4 | The performance of the LASSO logistic regression model for differentiating COVID-19 pneumonia from influenza A (H1N1) pneumonia. It was presented
by ROC curves in (A) the primary cohort (AUC, 0.963; 95% CI: 0.942–0.984) and (B) the validation cohort (AUC, 0.943; 95% CI: 0.900–0.986).
In conclusion, CT characteristics, including the distribution to participate in this study in accordance with the national
pattern and category of pulmonary opacity, combined legislation and the institutional requirements.
with clinical features, can help the early differentiation
of COVID-19 pneumonia from H1N1 pneumonia. CT AUTHOR CONTRIBUTIONS
manifestations of peripheral distribution patterns, together
with older age, low-grade fever, and slightly elevated AST, W-YS and FS participated in the study design and
indicate COVID-19 pneumonia; however, CT presentations conceptualization. W-YS, X-LZ, H-LZ, SZ, and S-PH participated
of peribronchovascular distribution patterns, centrilobular in the acquisition of data. W-YS, FS, and Y-HT participated in
nodule or tree-in-bud sign, consolidation, and bronchial analysis and interpretation of data. W-YS and Y-HT participated
wall thickening or bronchiectasis, together with younger in drafting of the manuscript and participated in the statistical
age, hyperpyrexia and higher level of AST, suggest analysis. FS and T-FL participated in critical revision of the
H1N1 pneumonia. manuscript for important intellectual content. Z-YZ, Y-XS, and
NX participated in administrative, technical, material support,
DATA AVAILABILITY STATEMENT and participated in study supervision. All authors contributed to
the article and approved the submitted version.
The original contributions presented in the study are included
in the article/supplementary material, further inquiries can be FUNDING
directed to the corresponding author/s.
This research was funded by the Novel Coronavirus Special
ETHICS STATEMENT Research Foundation of the Shanghai Municipal Science and
Technology Commission (Grant Number: 20441900600).
The studies involving human participants were reviewed and
approved by the Ethics Committee of Shanghai Public Health ACKNOWLEDGMENTS
Clinical Center, China, and Wuhan Union Red Cross Hospital,
China, and informed consent for this retrospective study was We thank all the doctors, nurses, disease control workers, and
waived (YJ-2020-S035-01). Written informed consent from researchers who have fought bravely and ceaseless against the
the participants’ legal guardian/next of kin was not required virus on the frontline during the COVID-19 epidemic.