Introduction

The coronavirus disease 2019 (COVID-19) rapidly invaded the world affecting millions of people and becoming a global health emergency.

Sparse autopsies on patients with COVID-19 have found interstitial pneumonia with diffuse alveolar damage, pulmonary arterioles thrombosis [1], and right ventricle (RV) dilation [2], suggesting that increased RV afterload due to endothelial injury [3] with lung vessel micro-thrombosis [4] might be a pathological driver in critical COVID-19 illness.

In two recent studies [5, 6] on 120 and 115 COVID-19 patients undergoing echocardiography, non-survivors displayed larger right-heart chambers, reduced RV function, and elevated pulmonary artery systolic pressure compared to survivors. In light of the possible impact of RV afterload on patients’ outcome and treatment [7,8,9], more robust data coming from larger populations are needed to strengthen the current evidence. Data from echocardiography are scarce, because of technical [5] and practical limitations, mainly related to limited resources in a scenario of health emergency and to the need of reducing healthcare workers’ exposure.

Chest CT has been widely adopted for COVID-19 pneumonia diagnosis [10], monitoring, and prognostication [11]. Beyond the extraction of parameters describing the features of lung involvement and its severity, the main pulmonary artery diameter (MPAD) can be easily measured from a chest CT scan. Enlarged MPAD on CT has been studied as a method for the screening and diagnosis of pulmonary hypertension (PH) [12,13,14,15]. Reference values for healthy patients have been previously established [16] and showed excellent sensitivity for excluding PH [13]. In the pre-COVID-19 era, enlargement of MPAD resulted to be associated with PH severity at RV catheterization and to higher rate of mortality [13]. Hence, we hypothesized that the measurement of MPAD, a known marker of pulmonary hypertension, could be useful for risk stratification of COVID-19 patients.

The aim of the present study was to evaluate the prognostic value of MPAD enlargement at the admission chest CT on overall survival in COVID-19 patients considering a multivariable setting, and to validate the final multivariable model.

Materials and methods

Study setting and participants

This is a multicenter retrospective cohort study. Study participants were recruited among all adult patients (age 18 years or older) with confirmed RT-PCR for SARS-CoV-2, submitted to chest CT within 72 h from admission in seven tertiary level hospitals located in Northern Italy ((1) IRCCS San Raffaele Scientific Institute, Milan (OSR); (2) Guglielmo da Saliceto Hospital, Piacenza (PCZ); (3) Ospedale Maggiore, Bologna (BOL); (4) ASST Bolognini Hospital Bergamo Est (BGE); (5) Ospedale di Cremona (CRE); (6) San Gerardo Hospital, Monza (SGM); (7) Centro Cardiologico Monzino IRCCS, Milano (CCM)), between March 1 and April 20, 2020. Missing data on comorbidities or follow-up were considered exclusion criteria. The study was approved by the local ethics committees and written informed consent was obtained.

Clinical data were collected by each center according to a centralized electronic case report form (CRF). Clinical data at the admission consisted of the following: demographic characteristics (sex and age), comorbidities (hypertension, diabetes, chronic lung disease, cardiovascular disease), and laboratory tests: white blood cell count (WBC), creatinine, C-reactive protein (CRP), lactate dehydrogenase (LDH), troponin I, interleukin-6, and D-dimer. Outcome data were oro-tracheal intubation and death.

All chest CT images were collected and analyzed in a single core-lab (Experimental Imaging Center, IRCCS San Raffaele Scientific Institute, Milan, Italy) blinded to clinical data.

The derivation cohort consisted of patients from the three hospitals (OSR, PCZ, BOL) that completed the CRF and shared chest CT images for centralized analysis within May 20, 2020, established as the first deadline for data collection. Patients from the remaining four participating hospitals (BGE, CRE, SGM, CCM), who provided data and CT images within the end of June, were used for the external validation. Prior to the analysis, data were cross-checked with medical charts and verified by data managers and clinicians for accuracy.

A total of 1469 consecutive patients fulfilled the inclusion criteria: 68 of the 829 patients of the derivation cohort and 7 of the 640 patients of the validation cohort were excluded for lacking data on comorbidities or follow-up (Fig. 1). The median follow-up was 51 days for the derivation cohort and 63 days for the validation cohort.

Fig. 1
figure 1

Enrollment flowchart

Complete demographic and clinical features of both cohorts are presented in Table 1.

Table 1 Patient characteristics

Outcomes

The primary outcome was overall survival (OS) measured from admission to the emergency department. Secondary outcomes were in-hospital mortality and time-to-orotracheal intubation from admission to the emergency department.

Chest CT scan

All chest CT examinations were performed on multidetector scanners with at least 16 detector rows. The list of scanners and acquisition protocols are reported in Supplementary Materials (ESM 1). Briefly, all volumetric chest scans were sent to the core-lab where they were reformatted at a 2.5-mm slice thickness without gap or overlap. Images were reconstructed with sharp kernel for lung parenchyma evaluation and a soft kernel for mediastinum evaluation, and they were visualized using a standard window (lung: width 1400 HU, center − 450 HU; mediastinum: width 350 HU, center 40 HU).

Chest CT image analysis

The analysis was performed by a radiologist with 9 years of experience in cardiothoracic imaging, blinded to all clinical data in order to reduce the possible bias in the analysis and to increase homogeneity and reliability of data. Parameters of lung involvement (residual respiratory lung reserve, pneumonia extension, and features) and pulmonary artery metrics were extracted (Fig. 2).

Fig. 2
figure 2

Imaging parameters extracted from chest CT in COVID-19 patients. For each patient, the diameter of the main pulmonary artery (red line) was measured at the level of its bifurcation and the diameter of right and left pulmonary arteries at the level of their origin (green and blue lines, respectively) on the mediastinal window reconstruction (top left of the panel). Using dedicated software, well-aerated lung volume was automatically extracted (violet parenchyma in 3D lung volume rendering on the top right of the panel). Pneumonia involvement was scored from 0 to 4 (score 0: 0%; score 1: 1–25%; score 2: 26–50%; score 3: 51–75% score 4: > 75%); well-aerated lung parenchyma is displayed in red, and pneumonia in green for the right lung and in blue for the left lung (bottom of the panel)

Residual respiratory lung reserve was measured as well-aerated lung volume, quantified using automatic software (IntelliSpace Portal 8.0, Philips Medical Systems) based on a previously defined HU threshold value [17]. Pneumonia extension was scored in 0% (absent, score 0), 1–25% (minimal, score 1); 26–50% (mild, score 2); 51–75% (moderate, score 3); and > 75% (severe, score 4). Qualitative features of pneumonia were scored as follows: score 0 for absent pneumonia, score 1 for prevalent ground-glass opacities (GGOs), score 2 for prevalent consolidation, and score 3 for GGOs and consolidation equally represented. The metrics of pulmonary arteries included the measurement of the MPAD at the level of its bifurcation, of the left (LPAD) and of the right pulmonary artery diameter (RPAD) at their origin. MPAD was classified according to the severity classification system for the diagnosis and prognosis of PH [13] in four classes: normal (≤ 27 mm for females and ≤ 29 mm for males, according to the Framingham sex-specific normative values [16]); mild enlargement (from > 27 to < 31 mm for females and > 29 to < 31 mm for males); moderate enlargement (from ≥ 31 to 34 mm for both sexes), and severe enlargement (> 34 mm for both sexes). To assess inter-observer reliability, a second reader with 5 years of experience in cardiothoracic imaging, blinded to the measurements of the first reader, assessed MPAD and pneumonia score.

Statistical analysis

Comparisons between numerical variables were performed with the Mann-Whitney test, and comparisons between categorical variables with Fisher’s exact test. P values were adjusted with Bonferroni’s correction to account for multiple testing. Follow-up data were censored at 40 days from hospital admission. OS was estimated with the Kaplan-Meier estimator and groups were compared with the log-rank test. The cumulative incidence of orotracheal intubation was estimated using the competing risk approach (with death without intubation as a competing event) and groups were compared with Gray’s test.

Cox’s regression analysis was employed for evaluating the role of MPAD classes [13, 16] in predicting the OS, un-adjusting and adjusting for demographics and comorbidities (age, sex, hypertension, diabetes, history of cardiovascular disease, and chronic lung disease). Multivariable Cox’s regression analysis was employed for evaluating the simultaneous role of CT parameters (MPAD, pneumonia extension, well-aerated lung volume), laboratory tests (WBC, CRP, creatinine), demographics, and comorbidities in predicting OS. The final model was obtained with a backward variable selection with a removal significant level of 0.05. Model calibration was performed as detailed in Supplementary Materials (ESM 1). Model discrimination was evaluated with the optimism corrected c-index. To assess the role of gender in categorized MPAD with respect to the outcome, the final Cox model was estimated adding an interaction term between the two variables.

The adjusted Cox regression analysis, the calibration curve, and the c-index of the final multivariable Cox’s model were computed on the validation cohort using CT data from two independent readers.

Inter-observer reliability was assessed with Rothery’s non-parametric intra-class correlation coefficient (ICC) [18] for MPAD and Cohen’s kappa coefficient for the pneumonia score.

Missing data were not imputed; thus, each analysis considered only complete cases for the variables used in the analysis. All tests were 2-sided and the significance level was set at 0.05. Confidence intervals (CIs) were computed at a 95% level. All statistical analyses were performed using R 3.5.0 (http://www.R-project.org/). All details are reported in Supplementary Materials (ESM 1).

Results

Derivation cohort: clinical and lung CT characteristics

The median age was 69 (interquartile range IQR, 58–77) years and 537 (70.6%) patients were male (Table 1). Hypertension was the most frequent comorbidity (484 patients, 63.6%) and median oxygen saturation was 92% (IQR, 88–95%).

At chest CT, 752 patients (98.8%) had pneumonia, with prevalent GGOs (432 patients, 56.8%). Pneumonia extension greater than 50% was recorded in 139 cases (18.3%).

Median MPAD was 27 [25–29] mm. Twenty-one (2.8%) patients were not hospitalized after emergency department evaluation, while 459 (60.3%) patients were discharged within 40 days from admission with a median hospitalization length of 11 (IQR, 7–17) days. The cumulative incidence of orotracheal intubation at 40 days was 21.28%. At 40 days, non-survivors were 182/761 patients (23.9%), with death occurring at a median time from an admission of 7 (IQR, 4–12) days. The characteristics of survivors and non-survivors are reported in Table 2. OS at 40 days was 76.08% (95% CI: 73.11–79.17%).

Table 2 Clinical and CT features of survivors vs non-survivors in the derivation cohort

Derivation cohort: pulmonary enlargement and mortality

Out of 761 patients, 86 (11.3%) had mild enlargement, 98 (12.9%) moderate enlargement, and 36 (4.7%) had severe MPAD enlargement.

The cumulative incidence of orotracheal intubation was not significantly different among the classes of MPAD, despite a numerically higher incidence in patients with a severe enlargement (at 40 days [95% CI]: for normal 21.44% [18.09–24.99%], for mild enlargement 17.44% [10.27–26.19%], for moderate enlargement 21.43% [13.91–30.02%], and for severe enlargement 27.78% [14.34–42.96%]; Gray’s test p = 0.609).

The mortality rate was 25.6% (22/86), 38.8% (38/98), and 47.2% (17/36) in patients with mild, moderate, and severe MPAD enlargement, respectively, while it was 19.4% (105/541) in patients with normal MPAD (global Fisher’s exact test p < 0.001). Indeed, the classes of MPAD were associated with a different risk of death (40-day overall survival [95% CI]: for normal 80.56% [77.32–83.99%], for mild enlargement 74.40%[65.72–84.23%], for moderate enlargement 61.22% [52.30–71.67%], and for severe enlargement 52.78% [38.75–71.89%]; log-rank test p < 0.001) (Fig. 3a). At Cox’s regression analysis, moderate to severe MPAD enlargement (≥ 31 mm) resulted to predict a significantly reduced overall survival with respect to normal MPAD (hazard ratio, HR [95% CI]: mild enlargement vs normal 1.387 [0.876–2.196], p = 0.163; moderate enlargement vs normal 2.311 [1.594–3.350], p < 0.001; severe enlargement vs normal 3.001 [1.797–5.010], p < 0.001). Similar results were obtained also at unadjusted Cox’s regression analysis (Table S1). Therefore, in the subsequent analyses, the MPAD was dichotomized into moderate-severe MPAD enlargement (≥ 31 mm) vs normal-mild MPAD enlargement (< 31 mm). The percentage of moderate-severe MPAD resulted to be similar between genders (M: 18.6% vs F: 15.2%, p = 1) (Table S2). The new categorization achieved HR = 1.741 (95% CI: 1.253–2.418, p < 0.001) at adjusted Cox’s regression (at unadjusted analysis: HR = 2.367 [1.725; 3.249], p < 0.001) (Fig. 3b).

Fig. 3
figure 3

Overall survival in COVID-19 patients according to four-tier (a) and two-tier (b) classification of main pulmonary artery enlargement. The Kaplan-Meier curves show the overall survival in COVID-19 patients according to the four-tier classification system of MPAD classification (a) [13]. The four classes of MPAD were associated with a different risk of death (40-day overall survival [95% CI]: 80.56% [77.32–83.99%] for normal MPAD, 74.40% [65.72–84.23%] for mild enlargement, 61.22% [52.30–71.67%] for moderate enlargement, 52.78% [38.75–71.89%] for severe enlargement; log-rank test p < 0.001). In b are reported the Kaplan-Meier curves (solid lines) and the adjusted Kaplan-Meier curves (dashed lines) with inverse probability weighting based on propensity score (detailed in Supplementary Materials ESM 1) showing the overall survival in COVID-19 patients according to two-tier MPAD enlargement (moderate-severe vs normal-mild). The two classes of MPAD enlargement were associated with a different risk of death (40-day OS [95% CI]: 79.74% [76.65–82.95%] for normal-mild enlargement MPAD and 58.96% [51.19–67.90%] for moderate-severe enlargement; log-rank test p < 0.001; adjusted log-rank test with inverse probability weighting based on propensity score p = 0.002)

Moderate-severe MPAD enlargement remained associated with a higher risk of death at multiple Cox’s regression analysis also when pneumonia extension, well-aerated lung volume, laboratory tests, demographics, and comorbidities were considered (HR [95% CI]: 1.592 [1.154–2.196], p = 0.005) (Table 3). Figure 4 reports two exemplifying cases with normally and moderately enlarged MPAD, respectively. Besides MPAD, the other variables associated with higher risk of mortality in the model were old age (HR [95% CI]: 1.063 [1.052–1.075], p < 0.001), male gender (1.586 [1.079–2.329], p = 0.019), higher creatinine (2.468 [1.824–3.339], p < 0.001), and greater pneumonia score (1.397 [1.084–1.799], p = 0.010), while greater well-aerated lung volume was a protective factor (0.699 [0.581–0.841], p < 0.001). The final model showed a good calibration (Figure S1) and discrimination (c-index [95% CI]: 0.826 [95% CI: 0.796–0.851]).

Table 3 Cox’s multivariable regression model for predicting overall survival estimated in the derivation cohort (c-index [95% CI] = 0.826 [0.796–0.851]; baseline survival at 28 days = 99.92%)
Fig. 4
figure 4

Exemplifying cases: a survivor (top) and non-survivor (bottom) patients. CT 3D volume rendering of lungs in antero-posterior view (a, e) and lateral view (b, f), CT axial images with mediastinal window (c, g) for main pulmonary artery diameter (MPAD) measurement (red lines in c, g), and 3D volumetric reconstruction of pulmonary arteries (d, h) in a survivor (top) and a non-survivor (bottom) COVID-19 patient, both males of 74 and 76 years, respectively, both with history of hypertension and treated ischemic cardiomyopathy, both suffering from fever (> 37.5 °C) and caught from 7 days. At the admission, WBC/mm3 were 4.36 and 3.8, CRP 0.93 and 1.1 mg/dL, LDH 295 and 378 U/L, and creatinine 0.93 and 1.1 mg/dL in survivor and non-survivor, respectively. The oxygen saturation was 90% and 94% for survivor and non-survivor, respectively, with higher well-aerated lung volume in the survivor patient (3882 mL vs 2382 mL, violet parenchyma on 3D lung volume rendering in a, b, e, f) with pneumonia (bright parenchyma in a, b, e, f), involving from 25 to 50% of lung volume in both cases. Pulmonary artery diameter was normal (27 mm) in the survivor (c, d) while it was enlarged (32 mm) in the non-survivor (g, h) who died 9 days after hospital admission

An interaction term between MPAD and gender was included in the final multiple Cox regression model and resulted not to be significant (p = 0.724), thus supporting that MPAD has a similar role in both sexes.

Validation of enlarged MPAD as a predictor of mortality

In the validation cohort, the median age was 69 years (IQR, 59–77 years) and 421 (66.5%) patients were males (Table 1). Hypertension was the most frequent comorbidity (285 patients, 45%) and median oxygen saturation was 92% (IQR, 88–96%).

Most of the patients (626, 98.9%) had pneumonia at chest CT, with prevalent GGOs (328 patients, 51.8%). Thirty-one percent of cases (196) had lung involvement greater than 50%.

Mild, moderate, and severe MPAD enlargement was found in 58 (9.2%), 101 (15.9%), and 33 (5.2%) patients, respectively.

Inter-reader agreement was excellent for both MPAD (ICC [95% CI]: 0.937 [0.931–1]) and pneumonia score (Cohen’s kappa [95% CI]: 0.886 [0.856–0.916]).

The mortality rate of the validation cohort was 22.1% (140/633 patients). OS at 40 days was 77.69% (95% CI = 74.50–81.03%).

At Cox’s regression analysis adjusted with respect to demographics and comorbidities, the risk of mortality associated with moderate-severe MPAD enlargement was slightly higher than in the derivation cohort (HR [95% CI]: 2.524 [1.784–3.571], p < 0.001). The discrimination of the model obtained in the derivation cohort was confirmed, using CT measurement from both the first and second reader (c-index [95% CI]: 0.789 [0.758–0.823] and 0.790 [0.753–0.825], respectively).

Discussion

The main finding of our study is that moderate-severe enlargement of the pulmonary artery diameter (MPAD ≥ 31 mm) at chest CT performed within 72 h from the admission is an independent predictor of death in COVID-19 patients. The multivariable model, coming from the analysis of 761 consecutive COVID-19 patients admitted to three third-level hospitals in Northern Italy, was confirmed on an external validation cohort including 633 patients from four other third-level hospitals in the same Italian territory.

Enlargement of MPAD was associated with a higher rate of in-hospital mortality: 19.4% (105/541) in patients with normal MPAD and 25.6% (22/86), 38.8% (38/98), and 47.2% (17/36) in patients with mild, moderate, and severe MPAD enlargement (global Fisher’s exact test p < 0.001).

Enlarged pulmonary artery on chest CT scan was previously demonstrated to be associated with pulmonary hypertension in other clinical settings before COVID-19 outbreak [12, 13, 19]; it is considered a consequence of elevated pressure, reflecting disease severity [13] and duration [20].

MPAD ≥ 31 mm was previously found to have excellent specificity (from 80 to 98%) and positive predictive value (from 85 to 98%) for the diagnosis of PH [13], with a 2–3-fold increased risk of mortality than in patients with normal MPAD [13].

Our results, based on a single time-point CT examination, do not allow to assert if the enlargement of MPAD is acute or not. However, the statistical adjustment for demographics and comorbidities suggested that the observed phenomenon is more likely an acute complication of COVID-19 pneumonia. Moreover, a previous CT study [21] on a limited sample size (44 COVID-19 patients) showed an MPAD enlargement with respect to a previous CT examination in non-survivors, further supporting the hypothesis that MPAD enlargement occurred in response to increased peripheral resistance for COVID-19-associated lung endothelial injury [3, 22, 23].

In the setting of non-COVID-19-related acute respiratory distress syndrome [24], PH resulted in highly prevalent (almost 50%) and entails higher risk of cardiac failure, shock [25, 26], and mortality [26, 27]. The etiology of acute PH is considered to be multifactorial, being related to the refractory hypoxemia, vasoconstriction, pulmonary edema, and microvascular thrombosis [28].

In the COVID-19 setting, where angiotensin-converting enzyme 2 on epithelial cells was recognized as a functional receptor for SARS-CoV-2, low levels of angiotensin and decreased downregulation of angiotensin II may contribute to determine acute PH by increasing pulmonary vasoconstriction [29]. Moreover, pulmonary small-vessel thrombosis could be the main cause of PH in COVID-19 pneumonia. In the autopsy series of COVID-19 patients [1, 2, 30], hallmarks of classic ARDS were not prominent [30], while a pauci-inflammatory thrombogenic vasculopathy [30], with small-vessel thrombosis [1, 2, 23], was reported as a common finding. In our cohort, similar to previous studies [4, 31], higher levels of IL-6, CRP, LDH, D-dimer, and troponin were found in non-survivors compared to survivors, further supporting the hypothesis of microvascular COVID-19 lung thrombo-inflammatory syndrome [4].

The measurement of the MPAD on chest CT is easy, fast (few seconds), highly reproducible, does not require any dedicated software, and could be a useful quantitative marker to be integrated into models for risk stratification. Moreover, the detection of enlarged MPAD on chest CT, being a known sign of PH, may have a potential impact on patient management and treatment. In fact, patients with PH may rapidly deteriorate with RV dysfunction and worsening oxygenation for ventilation-perfusion mismatch. In patients with severe pneumonia and dilated MPAD, alveolar ventilation could be optimized to limit hypoxic vasoconstriction, and prone ventilation preferred for minimizing positive end-expiratory pressure and hemodynamic impact [7]. Close monitoring of pulmonary vascular resistance should be considered to optimize inotropic support and pulmonary vasodilatation [7, 32], decreasing RV afterload.

The limitation of the present study was the lack of data about pulmonary artery pressure estimated by right heart catheterization or echocardiography. Despite echocardiography being commonly used to investigate pulmonary hypertension, it was not routinely performed in all COVID-19 patients. Moreover, its diagnostic accuracy is affected by body habitus, heart rate, acoustic window, operator’s experience [33], and notably by the presence of significant lung disease [5].

In our study, no differences in orotracheal intubation rate were found between patients with normal or enlarged MPAD, despite numerically higher rate in severely enlarged MPAD. This result might be influenced by bias of treatment in overwhelmed hospitals closely related to the scarcity of resources.

Non-survivor rates in our derivation and validation cohorts (23.9% and 22.1%, respectively) are consistent with the mortality rate recorded in Italy during the COVID-19 outbreak [34], suggesting that our cohort could be considered a representative sample of the Italian COVID-19 population.

In conclusion, in the present study, on a large cohort of COVID-19 patients from seven third-level hospitals from Northern Italy, we found that enlarged MPAD measured on chest CT at admission is an independent predictor of mortality. The MPAD can be easily measured on chest CT performed for pneumonia extension evaluation, offering additional prognostic information, and helping clinicians in patients’ management.