key: cord-0070555-0k2ypfpb authors: Cena, Tiziana; Cammarota, Gianmaria; Azzolina, Danila; Barini, Michela; Bazzano, Simona; Zagaria, Domenico; Negroni, Davide; Castello, Luigi; Carriero, Alessandro; Corte, Francesco Della; Vaschetto, Rosanna title: Predictors of intubation and mortality in COVID-19 patients: a retrospective study date: 2021-11-27 journal: J Anesth Analg Crit Care DOI: 10.1186/s44158-021-00016-5 sha: 48e5d74dfccf9123fbb3b921693e62c3f94076d0 doc_id: 70555 cord_uid: 0k2ypfpb BACKGROUND: Estimating the risk of intubation and mortality among COVID-19 patients can help clinicians triage these patients and allocate resources more efficiently. Thus, here we sought to identify the risk factors associated with intubation and intra-hospital mortality in a cohort of COVID-19 patients hospitalized due to hypoxemic acute respiratory failure (ARF). RESULTS: We included retrospectively a total of 187 patients admitted to the subintensive and intensive care units of the University Hospital “Maggiore della Carità” of Novara between March 1st and April 30th, 2020. Based on these patients’ demographic characteristics, early clinical and laboratory variables, and quantitative chest computerized tomography (CT) findings, we developed two random forest (RF) models able to predict intubation and intra-hospital mortality. Variables independently associated with intubation were C-reactive protein (p < 0.001), lactate dehydrogenase level (p = 0.018) and white blood cell count (p = 0.026), while variables independently associated with mortality were age (p < 0.001), other cardiovascular diseases (p = 0.029), C-reactive protein (p = 0.002), lactate dehydrogenase level (p = 0.018), and invasive mechanical ventilation (p = 0.001). On quantitative chest CT analysis, ground glass opacity, consolidation, and fibrosis resulted significantly associated with patient intubation and mortality. The major predictors for both models were the ratio between partial pressure of arterial oxygen and fraction of inspired oxygen, age, lactate dehydrogenase, C-reactive protein, glycemia, CT quantitative parameters, lymphocyte count, and symptom onset. CONCLUSIONS: Altogether, our findings confirm previously reported demographic, clinical, hemato-chemical, and radiologic predictors of adverse outcome among COVID-19-associated hypoxemic ARF patients. The two newly developed RF models herein described show an overall good level of accuracy in predicting intra-hospital mortality and intubation in our study population. Thus, their future development and implementation may help not only identify patients at higher risk of deterioration more effectively but also rebalance the disproportion between resources and demand. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1186/s44158-021-00016-5. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes a wide spectrum of clinical manifestations, named coronavirus disease , which range from asymptomatic infections to severe interstitial pneumonia. Unfortunately, the COVID-19 pandemic has led to a large number of critically ill patients in a very short time and is currently overwhelming the healthcare systems worldwide [1] . In this scenario, delays in intensive care unit (ICU) transfers [2] [3] [4] and in patient intubations [5, 6] have been associated with increased mortality among COVID-19 patients with hypoxemic acute respiratory failure (ARF). Thus, correct triage and prompt ICU allocation of patients scheduled to receive intubation in case of oxygen therapy or noninvasive ventilation failure are crucial to achieve an effective COVID-19 pandemic response [7] . Potential prognostic factors of intra-hospital mortality and intubation in patients hospitalized due to COVID-19related hypoxemic ARF include older age [8, 9] , male gender [9] , diabetes [9] , prothrombin time and D-dimer level [10] , lymphocytopenia [11] [12] [13] , leukopenia [11] [12] [13] , Creactive protein [11] [12] [13] , lactate dehydrogenase [11] [12] [13] , sequential organ failure assessment (SOFA) score [14] , arterial oxygen pressure on inspired oxygen fraction ratio (PaO 2 /FiO 2 ) [8, 14] , and acute kidney injury [14] . The Fleischner Society Statement on Chest Imaging and COVID-19-issued on April 7th, 2020-recommends chest computerized tomography (CT) imaging for the triage of (i) patients with suspected COVID-19 presenting with moderate to severe clinical features and a high pretest probability of disease, (ii) COVID-19 patients with worsening respiratory status, and (iii) patients with functional impairment and/or hypoxemia after recovery from COVID-19 [15, 16] . Fittingly, a singlecenter study has recently shown how the compromised lung volume estimated by quantitative CT analysis is a strong predictor of the need for oxygenation support and intubation among COVID-19 patients [17] . Thus, an algorithm combining demographic and early clinical characteristics together with laboratory findings and chest CT analysis results may favor the early identification of COVID-19 patients requiring invasive mechanical ventilation (IMV) or at increased risk of mortality. Our primary aim was to retrospectively identify the risk factors associated with intubation and intra-hospital mortality in a cohort of COVID-19 patients admitted to the hospital for hypoxemic ARF by analyzing their demographic characteristics, early clinical and laboratory variables, and quantitative chest CT analysis results. As a secondary endpoint, we sought to determine the performance of a newly developed algorithm based on the aforementioned variables in predicting the probability of IMV and intra-hospital mortality in our study population. The present investigation is an observational retrospective single-center study. Ethical approval was issued by the Comitato Etico Interaziendale Novara, Italy (Chairperson Prof. G. Zulian) on May 20th, 2020 (Ethics Committee No. CE 121/20). The requirement for informed consent was waived due to the retrospective nature of the study. We analyzed 187 consecutive patients with COVID-19 pneumonia diagnosed with real-time reverse transcriptase-polymerase chain reaction (RT-PCR) nasopharyngeal swabs, subjected to chest CT images, and admitted to the subintensive and intensive care units of our hospital between March 1st and April 30th, 2020. Patients with poor-quality chest CT images were excluded. The study was reported in accordance with STROBE guidelines. Demographic information, body mass index (BMI), time from first-symptoms, comorbidity, date of admission to hospital, clinical laboratory on admission including the first PaO 2 /FiO 2 ratio, arterial oxygen saturation (SpO 2 ), blood cell counts (i.e., leukocyte and lymphocyte count), biomarkers of inflammation (i.e., lactate dehydrogenase, ferritin, C-reactive protein, procalcitonin, fibrinogen), glycemia, and troponin were collected. Furthermore, we recorded the type of oxygen assistance administeredi.e., standard oxygen therapy (low-flow oxygen nasal cannula, Venturi mask, non-rebreathing mask), noninvasive ventilation [continuous positive airway pressure (CPAP) or bilevel positive airway pressure (BiPAP)], or invasive mechanical ventilation. All data were derived from both electronic hospital records and digitization of paper documents. Criteria for intubation were cardiac or respiratory arrest; inability to protect the airway; coma or psychomotor agitation; unmanageable secretions or uncontrolled vomiting; life-threatening arrhythmias or electrocardiographic signs of ischemia; hemodynamic instability, defined as systolic arterial pressure < 90 mmHg despite adequate filling or use of vasoactive agents; intolerance to all interfaces; dyspnea during noninvasive ventilation administered as CPAP or BiPAP; respiratory rate > 30 breaths/ min; SpO 2 < 92% during CPAP or BiPAP; and acidosis. CT scan was performed within 1 day from admission. CT images were independently reviewed by two radiologists with 10 and 14 years of clinical experience: all radiologists were blinded to the clinical status of the patients. The lung parenchyma segmentation was performed through a software-based evaluation on a dedicated workstation using the open-source 3D Slicer software (Fig. 1 ). More details can be found in the Supplementary Information. A sample size was computed to ensure a predictive ability of the RF model close to 0.8 with a margin of error in the sample estimates d = 0.05. More details can be found in the Supplementary Information. Descriptive statistics were reported as median and interquartile range for continuous variables and percentages (absolute numbers) for categorical variables. Missing values were handled leaving null the estimate. The logistic regression model, odds ratio (OR) together with the 95% confidence intervals (95% CI), and p values were reported for each predictor, considering separately their association with intubation and intra-hospital mortality. The analyses were performed using R software (version 0.2) with the packages caret and rms. To identify predictors of intra-hospital intubation or mortality, a random forest (RF) algorithm was employed. The variables having less than 10% of missing values were included in the predictive tool. More details can be found in the Supplementary Information. The main characteristics of the 187 patients included in the study are summarized in Table 1 . One hundred forty patients (∼75%) were males, with a median age of 64 years, a median BMI of 28 kg/m 2 , and a median PaO 2 / FiO 2 of 258 mmHg on admission. The median symptom onset was 7 days prior to admission. The most frequent comorbidities were hypertension (51%) and diabetes (24%). Main laboratory and chest CT findings on admission are also reported in Table 1 . Forty-five patients (24%) received standard oxygen therapy, 86 patients (46%) received noninvasive ventilation, and 56 patients (30%) were admitted to ICU and required IMV. The patients' demographic, clinical, and laboratory characteristics along with chest CT analysis findings relative to both the whole study population and alive vs death subjects are listed in Table 1 . When stratifying patients according to mortality, variables independently associated with mortality were age, other cardiovascular diseases except for coronary artery disease, C-reactive protein, lactate dehydrogenase levels, and IMV. Furthermore, on quantitative chest CT examination, we found a significant positive association between GGO and consolidation and fibrosis. Conversely, female sex, PaO 2 /FiO 2 , SpO 2 , lymphocytes count, standard oxygen therapy, and evidence of well-aerated lung parenchyma on chest CT scan were all inversely associated with mortality. Table 2 enlists the patients' demographic, clinical, and laboratory characteristics along with the chest CT analysis findings relative to both the whole population and intubated vs non-intubated patients. Among our cohort of 187 patients, 29 patients were excluded from this analysis because classified as "do-not-intubate" subjects, i.e., patients deemed ineligible for intubation in case of CPAP or BiPAP failure. Variables independently associated with intubation were C-reactive protein, lactate dehydrogenase levels, and white blood cell count. On quantitative chest CT analysis, GGO, consolidation, and fibrosis resulted positively associated with intubation, while PaO 2 /FiO 2 , SpO 2 , lymphocyte count, and wellaerated parenchyma were inversely associated with intubation. Lastly, for each logistic regression and for the continuous variables the cut points maximizing the best predictive value along with their corresponding area under the curve (AUC) are shown in Tables 1 and 2 . We next evaluated the importance of the variables encompassed in the RF algorithms for intra-hospital mortality and intubation prediction. In the RF model for mortality prediction, the most important nodes (importance > 50) were C-reactive protein, age, PaO 2 /FiO 2 , glycemia, SpO 2 on admission, a well-aerated parenchyma, lactate dehydrogenase, GGO, lymphocytes, other consolidation/fibrosis, and symptom onset ( Fig. 2A) . In the RF tree model for intubation prediction, the most important nodes (importance > 50) were PaO 2 /FiO 2 , a well-aerated parenchyma, C-reactive protein, GGO, other consolidation/fibrosis, glycemia, lactate dehydrogenase, lymphocytes, and age (Fig. 2B) . The variables having more than 10% of missing values were excluded from the predictive tool (i.e., BMI, fibrinogen, procalcitonin, ferritin, and troponin). The balanced accuracy in predicting intra-hospital mortality was 0.89 (κ value = 0.72; AUC = 0.73) (Fig. 3A) , whereas the balanced accuracy in predicting intubation was 0.9 (κ value = 0.75; AUC = 0.74) (Fig. 3B) . When the quantitative CT analysis variables were removed from both RF models, the accuracy of the model predicting intra-hospital mortality dropped to 0.75, whereas that of the model predicting intubation fell to 0.69 (Fig. 4A, B ). The main findings of our investigation can be summarized as follows: (i) in our cohort of COVID-19 patients, elderly male subjects with comorbidities, such as other cardiovascular diseases, intubated for severe hypoxemic ARF associated with a major inflammatory response and a widespread pulmonary involvement on chest CT were at increased risk of intra-hospital mortality; (ii) excluding those patients classified as "do-not-intubate," subjects with severe hypoxemic ARF experiencing increased inflammatory response and poor lung aeration on chest CT were at increased risk of intubation; and (iii) our novel RF algorithms performed well in predicting intrahospital mortality and intubation in our study population. The intra-hospital mortality rate of critically ill patients admitted for COVID-19 has been reported to range from 17 to 67% [18] . Well recognized risk factors associated with low survival or poor outcome in ICU are: male gender, increasing age, comorbidities such as diabetes, hypercholesterolemia, chronic obstructive pulmonary disease, IMV at high positive end-expiratory pressure, low PaO 2 /FiO 2 on ICU admission, high SOFA score, acute kidney injury, reduced respiratory system compliance, late pulmonary infections, and cardiovascular complications [14, 19] . In our setting, intra-hospital mortality rate was 31%, in line with previous reports [18] . Risk factors for intrahospital mortality identified in our cohort confirmed all previously reported predictors [10] [11] [12] [13] . Furthermore, increased C-reactive protein and lactate dehydrogenase and reduced lymphocyte count were all associated with increased mortality in our study population, which is in good agreement with previous data showing a positive association between severity/mortality rate of COVID-19 illness and biomarkers such as C-reactive protein, lactate dehydrogenase, and lymphopenia [12, 13] . Of note, the intubation rate in our study was 35.4%, which is consistent with the IMV incidence range among COVID-19 patients (12-33.1%) [20] [21] [22] . Older age, BMI, comorbidities-i.e., hypertension, diabetes, and cardiovascular diseases-, shortness of breath, SpO 2 < 90%, and increased respiratory rate are well-known predictors of intubation in patients admitted for COVID-19 [9] . In our cohort of patients, we confirm that reduced SpO 2 and/or PaO 2 /FiO 2 are risk factors for IMV. We also show that increased C-reactive protein and lactate dehydrogenase serum concentrations and elevated total white blood cell count are predictors of intubation, which is in good agreement with previous reports demonstrating the association between the aforementioned biomarkers and illness severity [12, 13] . Lung aeration loss on chest CT scan has been previously shown to be an independent predictor of death and ICU admission in COVID-19 patients suffering from hypoxemic ARF [17, 23] . Fittingly, we found that reduced aerated lung volume and increased GGO and/or consolidation and fibrosis are indicators of poor outcome and intubation. In this regard, Colombi et al. [23] have previously demonstrated an association between mortality and exudative consolidation, which may be suggestive of concomitant bacterial infection associated with death in COVID-19 patients [24] . To date, several models have been proposed to estimate the risk of COVID-19 patients to be hospitalized or to experience a poor outcome from the infection in order to assist medical staff in triaging patients when allocating limited healthcare resources [25] . With particular regard to predictive models for mortality and progression to a more severe or critical condition, the most frequently used predictors include comorbidities, age, sex, lymphocyte count, C-reactive protein, body temperature, creatinine, and imaging features [25] . The discrimination of these models ranged from 0.68 to 0.98 for intra-hospital mortality [26] and from 0.73 to 0.99 for worsening to a more critical state [27] . Here, we propose two novel RF models developed by including all the demographic, clinical, hematochemical, and radiological variables from our cohort of COVID-19 patients having less than 10% of missing data. The predictive balanced accuracy was high for both RF models, probably because the number of nodes exceeding an importance of 50 was very high for each algorithm. Among the items included in our RF algorithms for prediction of mortality and intra-hospital intubation, blood glucose and the symptom onset duration were the two factors that, in addition to the predictors listed above, showed an importance > 50. Our findings are in keeping with recent results suggesting that hyperglycemia, even in the absence of frank diabetes, is associated with a negative outcome compared to normoglycemic individuals as well as to those with pre-existing diabetes and COVID-19 [28] . Also, the symptom onset duration was confirmed to be a poor outcome predictor, being a fever lasting more than 7 days from onset of illness associated with increased ICU admission [29] . Although our study confirms with an innovative approach (i.e., RF) the risk factors of intubation and mortality found in the recent literature, it has several limitations. First, due the retrospective nature of the present singlecenter study, our results lack of generalizability. Second, no power sample was estimated for RF model accuracy assessment. Thus, no definitive conclusions can be drawn on the precision of our algorithm for intubation and intra-hospital mortality prediction. Third, some variables were excluded during the algorithm construction due to missing data occurrence > 10% and the fact that other laboratory values, such as PaCO 2 , creatinine, and D-Dimer were not collected. Lastly, our prediction models were not validated before the present investigation. Therefore, a future prospective investigation addressing the validation of our models is clearly needed. Ethical aspects of the COVID-19 crisis: how to deal with an overwhelming shortage of acute beds Severe community-acquired pneumonia. Assessment of severity criteria The impact of a delay in intensive care unit admission for community-acquired pneumonia Outcomes of early, late, and no admission to the intensive care unit for patients hospitalized with community-acquired pneumonia Benefits and risks of success or failure of noninvasive ventilation Changing use of noninvasive ventilation in critically ill patients: trends over 15 years in francophone countries Managing ICU surge during the COVID-19 crisis: rapid guidelines Predictors of intubation in COVID-19 patients treated with out-of-ICU continuous positive airway pressure Factors Associated with intubation and prolonged intubation in hospitalized patients with COVID-19 A cohort study of 676 patients indicates D-dimer is a critical risk factor for the mortality of COVID-19 Clinical characteristics of coronavirus disease 2019 in China Hematologic, biochemical and immune biomarker abnormalities associated with severe illness and mortality in coronavirus disease 2019 (COVID-19): a metaanalysis Biomarkers and outcomes of COVID-19 hospitalisations: systematic review and meta-analysis Factors influencing liberation from mechanical ventilation in coronavirus disease 2019: multicenter observational study in fifteen Italian ICUs CT image visual quantitative evaluation and clinical classification of coronavirus disease (COVID-19) Leung AN (2020) The role of chest imaging in patient management during the COVID-19 pandemic: a multinational consensus statement from the Fleischner Society Quantitative chest CT analysis in COVID-19 to predict the need for oxygenation support and intubation Outcomes of COVID-19 patients treated with continuous positive airway pressure outside the intensive care unit Risk factors associated with mortality among patients with COVID-19 in intensive care units in Lombardy. Italy Clinical characteristics of COVID-19 in New York City Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized with COVID-19 in the Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus-infected pneumonia in Wuhan. China Qualitative and quantitative chest CT parameters as predictors of specific mortality in COVID-19 patients Bacterial coinfections with SARS-CoV-2 Prediction models for diagnosis and prognosis of covid-19: systematic review and critical appraisal Development and external validation of a prognostic multivariable model on admission for hospitalized patients with COVID-19 A model to predict SARS-CoV-2 infection based on the first threemonth surveillance data in Brazil Hyperglycemia without diabetes and new-onset diabetes are both associated with poorer outcomes in COVID-19 Fever patterns, cytokine profiles, and outcomes in COVID-19. Open Forum We thank Dr. Arsura for the language revision of the manuscript. The study was conducted on human subjects; ethical standards were followed. Informed consent was not applicable since this is a retrospective study. The research was in accordance with the 1964 Helsinki Declaration and ethical approval received from Comitato Etico Interaziendale. All authors consent for publication. The authors declare that they have no conflicts of interest. 1 Anestesia e Terapia Intensiva, Azienda Ospedaliero Universitaria "Maggiore della Carità", Novara, Italy. 2 In our cohort of COVID-19 patients, we confirmed all demographic, clinical, hemato-chemical, and radiologic predictors of adverse outcomes previously reported. In addition, our innovative RF models based on the risk factors identified in our subset showed a good level of accuracy in predicting intra-hospital mortality and intubation. Thus, this approach may help to accurately identify patients at higher risk of deterioration, which would be particularly important in case of disproportion between resources and demand. Our results await further confirmation in larger multicentric prospective studies. The online version contains supplementary material available at https://doi. org/10.1186/s44158-021-00016-5. Authors' contributions TC, GC, and RV had the idea for and designed the study, had full access to all of the data in the study, contributed to the design of the study and to the interpretation of data for the work, and drafted the paper. DA contributed to formal analysis and interpretation data and drafted the paper. MB, SB, DN, DZ, and LC contributed to the acquisition, integrity, and analysis of the data. AC and FDC contributed to manuscript revision, read, and approved the submitted version. All authors (1) agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved, (2) revised the work critically for important intellectual content, and (3) gave final approval of the version to be published.