key: cord-0926796-z8xb31w3
authors: Dalpiaz, Giorgia; Gamberini, Lorenzo; Carnevale, Aldo; Spadaro, Savino; Mazzoli, Carlo Alberto; Piciucchi, Sara; Allegri, Davide; Capozzi, Chiara; Neziri, Ersenad; Bartolucci, Maurizio; Muratore, Francesco; Coppola, Francesca; Poerio, Antonio; Giampalma, Emanuela; Baldini, Luca; Tonetti, Tommaso; Cappellini, Iacopo; Colombo, Davide; Zani, Gianluca; Mellini, Lorenzo; Agnoletti, Vanni; Damiani, Federica; Gordini, Giovanni; Laici, Cristiana; Gola, Giuliano; Potalivo, Antonella; Montomoli, Jonathan; Ranieri, Vito Marco; Russo, Emanuele; Taddei, Stefania; Volta, Carlo Alberto; Scaramuzzo, Gaetano
title: Clinical implications of microvascular CT scan signs in COVID-19 patients requiring invasive mechanical ventilation
date: 2022-01-16
journal: Radiol Med
DOI: 10.1007/s11547-021-01444-7
sha: d4b47276b92d11c3f3bc58bd5d502bab3e24b805
doc_id: 926796
cord_uid: z8xb31w3

PURPOSE: COVID-19-related acute respiratory distress syndrome (ARDS) is characterized by the presence of signs of microvascular involvement at the CT scan, such as the vascular tree in bud (TIB) and the vascular enlargement pattern (VEP). Recent evidence suggests that TIB could be associated with an increased duration of invasive mechanical ventilation (IMV) and intensive care unit (ICU) stay. The primary objective of this study was to evaluate whether microvascular involvement signs could have a prognostic significance concerning liberation from IMV. MATERIAL AND METHODS: All the COVID-19 patients requiring IMV admitted to 16 Italian ICUs and having a lung CT scan recorded within 3 days from intubation were enrolled in this secondary analysis. Radiologic, clinical and biochemical data were collected. RESULTS: A total of 139 patients affected by COVID-19 related ARDS were enrolled. After grouping based on TIB or VEP detection, we found no differences in terms of duration of IMV and mortality. Extension of VEP and TIB was significantly correlated with ground-glass opacities (GGOs) and crazy paving pattern extension. A parenchymal extent over 50% of GGO and crazy paving pattern was more frequently observed among non-survivors, while a VEP and TIB extent involving 3 or more lobes was significantly more frequent in non-responders to prone positioning. CONCLUSIONS: The presence of early CT scan signs of microvascular involvement in COVID-19 patients does not appear to be associated with differences in duration of IMV and mortality. However, patients with a high extension of VEP and TIB may have a reduced oxygenation response to prone positioning. Trial Registration: NCT04411459 SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11547-021-01444-7.

COVID-19 is a multisystem syndrome due to cellular infection via the ACE 2 receptors by SARS-CoV-2. This syndrome most commonly involves the respiratory system, with 15% of cases showing severe respiratory failure and 5% requiring intensive care unit (ICU) admission for acute respiratory distress syndrome (ARDS) or multiorgan failure [1] .

Thrombotic complications are frequently observed in COVID-19 patients [2] and, especially for COVID-19 related ARDS, the observed incidence of pulmonary embolism (PE) and alveolar-capillary microthrombosis are significantly higher than for other ARDS causes [3, 4] .

A combination of conventional mechanisms of clot formation and immunothrombosis (particularly for smaller clots) seems to underlie these manifestations, and more recent insights about the interactions between SARS-CoV-2 and the vascular endothelium led to the definition of a COVID-19-associated coagulopathy.

This syndrome is characterized by mild thrombocytopenia, mildly prolonged prothrombin time, increased fibrinogen and raised D-dimer, all of which are more pronounced as disease severity increases [5, 6] .

From the histopathological point of view, other than the classical signs of ARDS such as diffuse alveolar damage and hyaline membrane formation, some distinctive pulmonary microvascular abnormalities occur in COVID-19, including intravascular fibrin deposition, perivascular monocyte infiltration, angiogenesis and microthrombi formation [4] .

Endothelial inflammation and microthrombosis can justify a series of CT scan abnormalities, such as vascular enlargement pattern (VEP) and vascular tree in bud (TIB), which have been increasingly described in COVID-19 CT scans [7, 8] .

The VEP is characterized by dilatation of pulmonary vessels around and within the lung abnormalities on CT images [8, 9] , and different descriptive approaches for VEP were adopted in the available literature, such as vascular thickening [10] , subsegmental vessels greater than 3 mm [11] or dilated small vessels in the lesion [12] .

It appears to be related to the vasodilatation and damage characterized by endothelial injury and enhanced permeability due to the massive accumulation of pro-inflammatory factors in the lung [13] . An increased representation of the VEP is observed along with the increase in the extent of ground-glass opacities (GGOs) and with the clinical severity of COVID-19 symptoms [8] .

On the other hand, TIB is the expression of pulmonary microthrombosis that enhances the peripheral pulmonary vessels on CT imaging, therefore producing the TIB pattern [14] . Vascular TIB identification was correlated with an increased duration of mechanical ventilation (MV) and length of hospital stay in a cohort of 39 critically ill invasively ventilated COVID-19 patients [7] .

It is not clear if these radiological signs represent only an expression of the extension of the COVID-19-related lung inflammation or if they could have an independent role in determining the severity of the clinical manifestations and the duration of invasive mechanical ventilation (IMV). Moreover, no evidence exists that patients responding differently to some of the COVID-19 related ARDS treatments, such as prone positioning, have different patterns of vascular involvement [15] .

Therefore, the primary objective of this study was to assess whether the presence of microvascular damage signs (VEP and TIB) in the chest CT scans acquired within 3 days before or after endotracheal intubation could be independently associated with prolonged IMV in critically ill COVID-19 patients.

Secondary objectives were to find eventual influences of the extension of these signs on the response to prone positioning and ICU outcome measures.

This is a secondary analysis of patients enrolled in a previous prospective study [16] initially conducted in 15 ICUs from Italian hospitals between February 22 and May 4, 2020. The data from another ICU obtained after the publication of the first paper were included in the present report.

The study was approved by the Institutional Review Board of the study coordinator centre (Maggiore Hospital, Bologna, Italy, approval number: 273/2020/OSS/AUSLBO) and by each institutional review committee of the participating hospitals. Informed consent was partially waived according to the approval of the local Ethics committee and analysis was conducted on anonymized individual data. The study was registered in ClinicalTrials.gov (NCT04411459).

The primary outcome measure was the number of ventilatorfree days (VFDs) at day 28 (see below for the full definition) between two groups based on the presence/absence of vascular TIB [7] or VEP.

Secondary outcome measures were the differences in VFD at day 28, PaO 2 /FiO 2 ratio change after pronation and ICU mortality between groups of patients showing different extensions of the microvascular and parenchymal signs. In particular, the microvascular involvement of 3 or more lobes for VES and TIB and the parenchymal involvement of over 50% for the parenchymal signs were used to define groups of patients with extended radiologic involvement.

Vascular TIB was defined as the presence of dilated peripheral subpleural pulmonary small vessels with a branching aspect resembling a budding tree [7, 17] .

VEP was defined as an ectasia of segmental and subsegmental vessels with a diameter larger than expected for the point within the vascular tree, characterized by: (a) vessel diameter larger than that in adjacent portions of the nondiseased lung, (b) vessel diameter larger than that in comparable regions of the non-diseased contralateral lung, or (c) focal dilatation or non-tapering of vessels as they course towards the lung periphery [18] .

Detailed images of TIB and VEP are available in Fig. 1 .

For the other signs, the definitions were referred to the Glossary of Terms for Thoracic Imaging of the Fleischner Society [19] .

Mechanical ventilation was considered invasive if delivered through an endotracheal tube or a tracheostomy cannula. The duration of IMV was defined as the time elapsed from intubation to successful extubation or successful permanent disconnection from mechanical ventilation for tracheostomized patients. This latter was considered effective if sustained 24 h per day. VFDs at 28 days were defined as 28 minus the number of days elapsed from the last successful extubation in intubated patients, whether or not NIV was used after extubation. VFDs were defined as zero in patients who died during the 28 days regardless of their extubation status. In tracheostomized patients, intermittent disconnections were not counted and VFD were defined as 28 minus the number of days from the last successful sustained disconnection from mechanical ventilation [20] .

To be included in the present study, patients should be tested positive for SARS-CoV-2 infection (confirmed by real-time reverse transcription-polymerase chain reaction assays), older than 18 years, undergo IMV, fulfil the criteria for ARDS according to the Berlin definition [21] , and have at least one CT scan available performed within three days before or after orotracheal intubation.

The single non-inclusion criterion was the use of noninvasive ventilation (NIV) during the entire clinical course.

The sample size was determined considering 5 VFD as a clinically important difference on a baseline mean of 11 ± 7.5 VFD [16] , and considering an ∝ -error of 0.05 and a β-error of 0.8 with a 1:1 ratio between the two groups, a sample size of at least 126 patients was required.

Baseline data and patient's history were collected using an electronic case report form developed by YGHEA, CRO division of Ecol Studio SPA (Bologna Operational Headquarters, Bologna, Italy) and hosted by Actide Nubilaria (Novara, Italy). Collected data included demographic data, underlying comorbidities, laboratory data, chest radiological reports, respiratory parameters before the intubation and during the first five days of IMV (static compliance-C RS , PaO 2 /FiO 2 ratio), and clinical outcomes.

Radiological data recorded consisted in the presence, extent and distribution of VEP, TIB and peripheral micronodules, scored on the basis of the number of lobes involved; GGO, crazy paving and consolidations, evaluated for percentage of lung involvement on a scale of 0 to 4 as follows: 0 = 0% involvement, 1 = less than 25%, 2 = 25% to less than 50%, 3 = 50% to less than 75%, 4 = 75% or greater involvement. The concomitant presence of pleural effusion and pulmonary artery dilatation (diameter greater than 29 mm) were also recorded.

Concerning the extent of the lesions, VEP, TIB and peripheral nodules were considered extended if they involved 3 or more lobes, while GGO, crazy paving and consolidations were considered extended if they involved more than 50% of the pulmonary parenchyma.

Concerning pronation-related data, the PaO 2 /FiO 2 was recorded within 3 h before prone positioning (T1) and from 1 to 3 h after re-supination (T2), all data were referred to the first pronation and both measurements were taken, therefore, in the supine position, and the duration of the PP session was recorded. A PaO 2 /FiO 2 increase of at least 49% was used as a cut-off value to identify responders to pronation [15] .

Chest CT studies were interpreted by the means of a consensus assessment from two radiologists with at least 5 years of experience for each centre. If consensus was not reached, a senior radiologist was consulted to address the final decision.

CT scan acquisitions were performed using different scanners in different institutions: GE VCT Lightspeed 64, GE VCT Lightspeed 16 (GE Healthcare, Buckinghamshire, UK), Philips CT Ingenuity 64, Philips Brilliance 64 (Philips Medical Systems, Best, The Netherlands), Toshiba Aquilion S16 TSX-101A and Toshiba Aquilion S32 TSX-101A (Toshiba Medical Systems Europe, Zoetermeer, The Netherlands).

The scans included non-contrast acquisition and were obtained in the supine position from lung apices to bases at full-suspended inspiration, with ≤ 1.25-mm section thicknesses, using standard acquisition parameters adjusted to patients' biometrics (100-120 effective mAs and 120-140 kVp). Images were reconstructed using both sharp kernels and visualized at window settings optimized for lung parenchyma (window width, 1500 HU; window level, -500 HU), and medium-soft kernels with soft-tissue window settings (window width, 300 HU; window level, 40 HU). Maximum intensity projection (MIP) images were used to facilitate the recognition of the TIB and VEP.

Data were analysed using Stata/CI 16 (StataCorp LLC, Texas, USA) and GraphPad Prism version 8.4.3 for Windows (GraphPad Software, San Diego, California, USA).

Continuous variables were tested for normality with the Shapiro-Wilk test and were expressed as mean and standard deviation (SD) or median and interquartile range [IQR] based on their distribution. Comparisons between continuous variables were made with the Student's t-test or the Mann-Whitney U test when appropriate. Categorical variables were expressed as numbers and percentages and compared using the Chi-square test or Fisher's exact test when needed.

The log-rank test was used to compare the cumulative incidences of liberation from mechanical ventilation up to day 28 between the groups of patients with and without TIB and VEP. Kaplan-Meier curves were also built for the same groups.

Finally, bivariate Spearman correlation coefficients were calculated taking into account the extension of the different radiological signs.

During the study period, 632 patients were assessed for eligibility. After excluding those who underwent exclusively NIV or with missing data, the original study enrolled a total of 470 patients. Of these, 139 had a CT scan performed within ± 3 days from orotracheal intubation and were therefore included in the final analysis; supplement Fig. 1 shows the flow of the patients throughout the study.

Globally, the TIB sign was observed in 75 out of the 139 patients (54%), while the VEP in 108 out of 139 (77.7%); Table 1 shows the main demographic and clinical variables of the patients included, while Table 2 describes their clinical course.

The TIB sign was associated with a significantly lower BMI and a small, but statistically significant, delay from hospital to ICU admission and orotracheal intubation; these differences were not observed concerning the presence of the VEP. Patients' demographics were not significantly different with regard to age, comorbidities, clinical severity at ICU admission (SAPS II and SOFA score) as well as for the clinical severity of COVID-19-related ARDS in terms of worst PaO 2 /FiO 2 ratio and lowest respiratory system compliance (C RS ) observed during the first five days of IMV for both the TIB and VEP comparisons.

Supplement Table 1 shows the initial laboratory and ventilatory data. The nearest D-dimer concentration from the CT scan as well as the peak concentration of IL-6 and LDH was not significantly different between the two groups; however, the prevalence of patients with LDH > 250 U/L [22] was significantly higher in the VEP group (99% vs 85.7%, p = 0.008). There were no significant differences in terms of responders to prone positioning either in the TIB (58.7% vs 55.8%, p = 0.832) and VEP groups (50.8% vs 73.1%, p = 0.063). Table 3 shows the extent and distribution of the main radiologic signs found in the CT scans included, while Fig. 2 represents the correlation matrix of the radiologic variables taken into account. In most cases, both TIB and VEP involved two or three lobes; the extents of the two signs were significantly intercorrelated and also correlated with the percentage of pulmonary parenchyma affected by GGO. The magnitude of the association was weak to moderate (rho coefficients range: 0.32-0.35, p < 0.01), TIB extent was also significantly correlated with the extent of peripheral micronodules (rho: 0.26, p < 0.01).

Vascular TIB was mostly observed within the GGO areas, while VEP was observed inside the GGO areas in 61.1% and both inside and outside the GGO areas in 34.3% of the cases.

Interestingly, pulmonary artery enlargement was significantly more frequently observed in patients with TIB sign, while this observation was not confirmed for VEP.

Duration of mechanical ventilation was similar between patients with and without TIB in CT scan (15 [8-23] versus 16 [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] , respectively, p = 0.297) as well as 28d-VFD (0 [0-15] vs 0 [0-12], p = 0.979). Similar results were found when dividing the two groups based on the presence of the VEP. Figure 3 shows the Kaplan-Meier graph for the cumulative incidence of liberation from mechanical ventilation within the first 28 days after orotracheal intubation based on the presence of the TIB sign or the VEP and the relative logrank tests; no differences were found in the two functions.

Finally, ICU cardiovascular complications and the need for renal replacement therapy were not significantly different in either of the two grouping modalities (VEP and TIB presence), as well as ICU mortality (Table 2 ). Figure 4 shows the distribution of extended radiological signs among different subgroups of patients.

Extended VEP and TIB signs were significantly more frequent in poor responders to pronation, while consolidations were more frequent in responders. Ground glass and crazy paving pattern extension over 50% of pulmonary parenchyma were significantly more frequent in non-survivors. No differences were observed among the radiologic signs concerning weaning from MV within day 28 and PaO 2 /FiO 2 ratio before orotracheal intubation.

Microvascular involvement is a key determinant of COVID-19 syndrome, and microvascular signs such as VEP and TIB are part of the typical radiological aspects of this syndrome and have an established role in its diagnosis [8, 12, 23] . 

Time from hospital admission to intubation-d (IQR) However, the clinical implications of these signs are still poorly investigated.

The main findings of this study were: (a) the evidence of VEP or TIB on CT scans performed within 3 days from orotracheal intubation was not associated with different clinical courses in terms of duration of MV or mortality; (b) no differences were either found about the median concentration of biomarkers related to disease severity and thrombosis; however, VEP was significantly associated with a slightly higher prevalence of LDH > 250 U/L and TIB was associated with an increased prevalence of pulmonary artery enlargement; (c) the extent of microvascular signs was significantly correlated with that of GGO, peripheral micronodules and crazy paving pattern, while no correlation was found with the extent of consolidations; (d) the involvement of 3 or more lobes by VEP and TIB was more frequent among nonresponders to prone positioning.

Recent literature found a significant association between the evidence of VEP or TIB in CT scan imaging and disease severity in terms of symptoms and biomarkers, impaired coagulation and an increase in the duration of MV and hospital stay [7, 23] .

Among COVID-19 patients showing vascular changes in CT scan, a significantly increased proportion of elevated C-reactive protein and LDH was observed [23] , and the elevation of these biomarkers was significantly associated with ICU admission and mortality in the available literature [1, 24] .

Our results showed no differences in median values of LDH; however, VEP was associated with a significantly increased proportion of patients with peak LDH values above 250 U/L during the first 7 days of hospitalization [22] .

On the other hand, TIB was associated with a significantly increased prevalence of pulmonary artery dilatation, a radiologic sign recently recognized as a predictor of clinical severity and mortality in COVID-19 patients [25, 26] .

However, despite the higher prevalence of negative prognostic factors in the VEP and TIB populations, we were not able to demonstrate any significant difference in terms of liberation from mechanical ventilation and mortality. Moreover, the number of lobes involved by VEP and TIB was correlated with a higher extent of the GGO (Fig. 2) , which has been already associated with worse outcomes in several studies [27] [28] [29] , and this aspect was also confirmed in our study where we found a significantly higher prevalence of extended GGO among non-survivors (Fig. 4) .

The pulmonary epithelial and endothelial damage due to SARS-CoV-2, the subsequent release of inflammatory cytokines and the proinflammatory/anti-inflammatory unbalance seem to be a shared mechanism between the development of the microvascular abnormalities, oedema and exudate formation [30] , of which GGOs are the radiologic expression [23, 31] .

Therefore, based on our results, VEP and TIB extent could be a partial epiphenomenon of GGO extent, without a prognostic significance in terms of duration of IMV or mortality in critically ill patients.

Interestingly, extended microvascular signs were significantly more frequent among non-responders to prone positioning, this aspect is in line with our previous research [15] , which found that the PaO 2 /FiO 2 increase observed after prone positioning in sustained responders could not be fully explained by differences in driving pressure. Possibly, disease phenotypes showing an extensive microvascular involvement could blunt the oxygenation response to prone positioning, while phenotypes showing extended consolidations rely more on parenchymal recruitment and ventilation/perfusion ratio redistribution for their oxygenation response. Unfortunately, prone positioning was performed in only 83 (59.7%) of the patients enrolled in the present study; therefore, we were limited in performing subgroup analyses concerning the anteroposterior distribution of vascular lesions or extent of consolidations. The mixed population of patients treated and not treated with prone positioning in this study could have also played a role concerning the observed liberation from MV.

Our results could have also been partially influenced by the selection criteria adopted for patients' enrolment because, since this study focused exclusively on critically ill COVID-19 patients, a higher prevalence of negative prognostic risk factors is expected, and the relative eventual prognostic value of vascular signs could have been smoothened in favour of other factors that notoriously play a role in critically ill patients such as the number of concomitant organ failures [16, 32] , hospital-acquired superinfections [33, 34] and other ICU complications [2, 35] .

Another limitation of this study is related to the different experimental treatments that were adopted for patients undergoing IMV during the study period and that could have partially influenced our results.

Moreover, we could not get more detailed lung perfusion images due to the lack of systematic acquisition of pulmonary CT angiograms in the patients included and the unavailability of dual-energy computed tomography (DECT) machines in the centres participating in this study, differently from previous works [7] . Finally, it should be noted that consolidations imply the increase in parenchymal attenuation that obscures the margins of vessels, and this may have influenced the prevalence of microvascular signs in patients demonstrating more consolidated lung areas.

Future research should further define if vascular signs are only a radiologic epiphenomenon of a common shared pathogenetic mechanism with GGO or have an independent clinical meaning and their eventual relation with COVID-19-related coagulopathy. 

Vascular TIB and VEP are frequently observed in the CT scans of critically ill COVID-19 patients requiring IMV. Their extent is significantly correlated with GGO and crazy paving pattern areas. Despite previous studies suggesting a potential prognostic role of these signs concerning the duration of IMV and survival, and the detection of a higher prevalence of prognostic factors associated with prolonged ventilation in our population, we were not able to confirm this hypothesis.

Patients not responding to prone positioning had a higher extent of microvascular involvement signs. More focused studies on the potential relation between vascular signs and oxygenation response to prone positioning could better clarify if these signs could help in choosing which patients may benefit from prone positioning. 

Radiology Collaborators (to be indexed and searchable into PubMed). Maggiore Hospital Carlo Alberto Pizzardi

Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: summary of a report of 72314 cases from the chinese center for disease control and prevention

Venous and arterial thromboembolic complications in COVID-19 patients admitted to an academic hospital in

High risk of thrombosis in patients with severe SARS-CoV-2 infection: a multicenter prospective cohort study

Pulmonary vascular endothelialitis, thrombosis, and angiogenesis in Covid-19

COVID-19, immunothrombosis and venous thromboembolism: biological mechanisms

Pathophysiology of COVID-19-associated acute respiratory distress syndrome: a multicentre prospective observational study

Pulmonary angiopathy in severe COVID-19: physiologic, imaging, and hematologic observations

CT features of coronavirus disease 2019 (COVID-19) with an emphasis on the vascular enlargement pattern

Pulmonary vascular enlargement on thoracic CT for diagnosis and differential diagnosis of COVID-19: a systematic review and meta-analysis

Performance of radiologists in differentiating COVID-19 from non-COVID-19 viral pneumonia at chest CT

Chest CT features of COVID-19 in

CT features of coronavirus disease 2019 (COVID-19) pneumonia in 62 patients in Wuhan, China

Overlapping and discrete aspects of the pathology and pathogenesis of the emerging human pathogenic coronaviruses SARS-CoV, MERS-CoV, and 2019-nCoV

Computed tomography vascular tree-inbud: a novel prognostic imaging biomarker in COVID-19?

Sustained oxygenation improvement after first prone positioning is associated with liberation from mechanical ventilation and mortality in critically ill COVID-19 patients: a cohort study

Factors influencing liberation from mechanical ventilation in coronavirus disease 2019: multicenter observational study in fifteen Italian ICUs

Tree-in-bud pattern in neoplastic pulmonary emboli

Pulmonary vascular manifestations of COVID-19 pneumonia

Fleischner society: glossary of terms for thoracic imaging

Reappraisal of ventilator-free days in critical care research

Acute respiratory distress syndrome: the Berlin definition

Prognostic factors for COVID-19 pneumonia progression to severe symptoms based on earlier clinical features: a retrospective analysis

A new imaging sign in COVID-19 pneumonia: vascular changes and their correlation with clinical severity of the disease

Elevated levels of IL-6 and CRP predict the need for mechanical ventilation in COVID-19

Chest CT-derived pulmonary artery enlargement at the admission predicts overall survival in COVID-19 patients: insight from 1461 consecutive patients in Italy

Pulmonary artery trunk enlargement on admission as a predictor of mortality in in-hospital patients with COVID-19

Assigning computed tomography involvement score in COVID-19 patients: prognosis prediction and impact on management

COVID-19 chest computed tomography to stratify severity and disease extension by artificial neural network computer-aided diagnosis

Vascular mechanisms and manifestations of COVID-19

COVID-19 and multiorgan failure: a narrative review on potential mechanisms

Analysis of critical care severity of illness scoring systems in patients with coronavirus disease 2019: a retrospective analysis of three UK ICUs

Nosocomial infections associated to COVID-19 in the intensive care unit: clinical characteristics and outcome

Ventilator-associated pneumonia in patients with covid-19: a systematic review and meta-analysis

Acute complications and mortality in hospitalized patients with coronavirus disease 2019: a systematic review and meta-analysis

LG, AC, GS and CAM drafted the article. GD, GS and CAV substantively revised the article. The ICU-RER COVID-19 collaboration was involved in data collection (see acknowledgments). All the authors have approved the submitted version and have agreed both to be personally accountable for the author's own contributions and to ensure that questions related to the accuracy or integrity of any part of the work, even ones in which the author was not personally involved.Funding This research did not receive any specific grant from funding agencies in the public, commercial or not-for-profit sectors.

The dataset related to this publication is available at https:// doi. org/ 10. 17632/ cmctw 928j4.1.

The authors have nothing to disclose.

The study was approved by the Institutional Review Board (IRB) of the study coordinator centre (Maggiore Hospital, Bologna, Italy, Approval Number: 273/2020/OSS/ AUSLBO) and by each institutional review committee of the participating hospitals. Informed consent was waived for unconscious patients while it was acquired for conscious patients or after liberation from mechanical ventilation or at the time of follow-up. The researchers analysed anonymized individual data.