Respiratory support in patients with severe COVID-19 in the International Severe Acute Respiratory and Emerging Infection (ISARIC) COVID-19 study: a prospective, multinational, observational study.

BACKGROUND: Up to 30% of hospitalised patients with COVID-19 require advanced respiratory support, including high-flow nasal cannulas (HFNC), non-invasive mechanical ventilation (NIV), or invasive mechanical ventilation (IMV). We aimed to describe the clinical characteristics, outcomes and risk factors for failing non-invasive respiratory support in patients treated with severe COVID-19 during the first two years of the pandemic in high-income countries (HICs) and low middle-income countries (LMICs).
METHODS: This is a multinational, multicentre, prospective cohort study embedded in the ISARIC-WHO COVID-19 Clinical Characterisation Protocol. Patients with laboratory-confirmed SARS-CoV-2 infection who required hospital admission were recruited prospectively. Patients treated with HFNC, NIV, or IMV within the first 24 h of hospital admission were included in this study. Descriptive statistics, random forest, and logistic regression analyses were used to describe clinical characteristics and compare clinical outcomes among patients treated with the different types of advanced respiratory support.
RESULTS: A total of 66,565 patients were included in this study. Overall, 82.6% of patients were treated in HIC, and 40.6% were admitted to the hospital during the first pandemic wave. During the first 24 h after hospital admission, patients in HICs were more frequently treated with HFNC (48.0%), followed by NIV (38.6%) and IMV (13.4%). In contrast, patients admitted in lower- and middle-income countries (LMICs) were less frequently treated with HFNC (16.1%) and the majority received IMV (59.1%). The failure rate of non-invasive respiratory support (i.e. HFNC or NIV) was 15.5%, of which 71.2% were from HIC and 28.8% from LMIC. The variables most strongly associated with non-invasive ventilation failure, defined as progression to IMV, were high leukocyte counts at hospital admission (OR [95%CI]; 5.86 [4.83-7.10]), treatment in an LMIC (OR [95%CI]; 2.04 [1.97-2.11]), and tachypnoea at hospital admission (OR [95%CI]; 1.16 [1.14-1.18]). Patients who failed HFNC/NIV had a higher 28-day fatality ratio (OR [95%CI]; 1.27 [1.25-1.30]).
CONCLUSIONS: In the present international cohort, the most frequently used advanced respiratory support was the HFNC. However, IMV was used more often in LMIC. Higher leucocyte count, tachypnoea, and treatment in LMIC were risk factors for HFNC/NIV failure. HFNC/NIV failure was related to worse clinical outcomes, such as 28-day mortality. Trial registration This is a prospective observational study; therefore, no health care interventions were applied to participants, and trial registration is not applicable.

Background

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has infected over 500 million people worldwide and resulted in more than 6 million deaths (https://covid19.who.int) [1, 2]. COVID-19, the disease caused by the SARS-CoV-2, is a multisystemic disease [3]. Its most severe presentation is acute respiratory distress syndrome (ARDS), secondary to pneumonia [4-6]. Most critically ill patients with COVID-19 receive advanced respiratory support, defined as high-flow nasal cannula (HFNC), non-invasive mechanical ventilation (NIV), or invasive mechanical ventilation (IMV) [3, 7, 8]. Up to 30% of hospitalised patients with COVID-19 are treated with one of these interventions [9, 10]; however, the use and need for support have changed over time depending on COVID-19 vaccination coverage, circulating viral variants, an evolving treatment evidence base and practice variation [11, 12].

Given the high demand for respiratory support and the insufficient capacity of intensive care units (ICU) and resources during the pandemic, especially in low- and middle-income countries (LMIC), the use of less invasive alternatives emerged as an alternative to provide advanced respiratory support [13, 14]. A global survey in 2020 found that HFNC (54%) and NIV (47%) were the most frequently used types of advanced respiratory support in patients with severe COVID-19 [15]. Up to 37% of patients who received NIV support ultimately required IMV [16], with high fatality ratios, especially in Latin America [17].

The objectives of this global study are to describe the clinical characteristics and outcomes of patients treated with HFNC, NIV, and IMV during the first two years of the pandemic, to determine risk factors associated with HFNC and NIV failure, and to estimate the association of later administration of IMV on clinical outcomes. We also compare the respiratory support types used in high-income countries (HICs) with those used in LMIC.

Methods

This is a prospective observational study of hospitalised patients from five continents. The study Consortium framework is provided by the International Severe Acute Respiratory and Emerging Infection (ISARIC)—World Health Organization (WHO) Clinical Characterisation Protocol for Severe Emerging Infections [18, 19]. The protocol, case report forms, consent forms, and study information are available on the ISARIC website (https://isaric.org). This standardised protocol uses tiered data collection tailored to a range of resource settings [19]. Investigators from 69 countries collected prospective data using the ISARIC case report form (CRF) built on Research Electronic Data Capture (REDCap, version 8.11.11, Vanderbilt University, Nashville, Tenn.) [20] hosted by the University of Oxford. Other investigators collected data using locally hosted systems and submitted it to the ISARIC dataset for centralised mapping. All investigators retain full rights to their data.

This observational study required no change to clinical management and encouraged patient enrolment in other research projects. The ISARIC-WHO Clinical Characterisation Protocol was approved by the World Health Organization Ethics Review Committee (RPC571 and RPC572). Also, local ethics approval was obtained for each participating country and site according to local requirements.

Study population

We included hospitalised patients with confirmed SARS-CoV-2 infection by reverse transcription-polymerase chain reaction (RT-PCR) in a respiratory sample treated with advanced respiratory support, defined as either HFNC, NIV, or IMV [3]. Patients with no recorded demographic data or vital signs within the first 24 h of hospital admission were excluded, as were patients whose 28-day vital status was unknown.

Variables and measurement

We recorded age, sex, income classification according to the World Bank (https://data.worldbank.org/country) of the country of recruitment, comorbidities, vital signs on admission, laboratory measurements during the first 24 h of hospital admission, treatment with advanced respiratory support at any point during hospitalisation, systemic complications, and treatments used during hospitalisation. The case report form completion guide is available online (https://isaric.org).

We identified the first wave of the pandemic for each participating country and composed a dichotomous variable to evaluate the impact of being admitted during the first wave on clinical outcomes.

We stratified patients in the cohort based on the first type of respiratory support received within the first 24 h of hospital admission. High-flow nasal cannula (HFNC) was defined as respiratory support continuously applied through large-bore nasal prongs using a heated and humid gas flow at an initial flow more significant than 20 L/min (or up to 80 L per minute) and a fraction of inspired oxygen of up to 1.0. Non-invasive mechanical ventilation (NIV) was defined as any type of positive pressure therapy delivered through a fitted mask and was preferred in patients with oxygen requirements over 6–15 L/min or laboured breathing. Continuous positive pressure (CPAP) or bi-positive pressure (BiPAP) may occur and be considered NIV. Invasive mechanical ventilation (IMV) is any mechanical ventilation administered to the patient after endotracheal intubation or tracheostomy. The decision to use this modality was left to the health care providers and not per study protocol.

Patients were considered to have failed the non-invasive respiratory strategy (i.e. HFNC or NIV) if they were subsequently treated with IMV during hospitalisation.

Outcomes

The primary outcome evaluated in this study was 28-day mortality. Secondary outcomes included the rate of and risk factors for failing non-invasive respiratory support (i.e. HFNC or NIV), the association of failure with clinical outcomes, and the frequency of respiratory strategies used in HIC versus LMICs.

Statistical methods

Continuous variables were expressed as median (interquartile range), and categorical variables as counts (percentages). For the primary outcome of 28-day mortality and secondary outcome of non-invasive respiratory failure, random forest (RF) models were used to identify the factors associated with these outcomes. The RF model uses multiple randomised individual decision trees that operate as an ensemble, where each decision tree gives a predicted class. The class obtained most frequently among the decision trees becomes the RF model prediction. A total of 500 estimators were used in this model. A more detailed explanation of the RF models is presented in the supplement.

To evaluate the performance of the RF model, the area under the model's receiver operating characteristics curve (AUROC) was used; for this, a tenfold cross-validation method was used, in which the data set was divided into ten subsets, and the validation was repeated ten times. Each time, one of the subsets was used as the test cohort, and the other nine subsets were used as training subsets, then the average AUROC was calculated and reported. When used for classification, RF models perform an implicit feature selection, a general indicator of each specific feature relevance, and can be computed as the Gini importance.

Then, we fitted two multivariable logistic regression models to estimate associations with the risk of 28-day fatality ratio or non-invasive respiratory failure, respectively. Variables identified as relevant by the RF model were included as explanatory variables. Odds ratios (ORs) were presented with forest plots.

A patient treated with respiratory support might receive different strategies during hospital admission. Thus, we developed alluvia diagrams to understand how patients were treated with other respiratory methods over time, stratified by the countries’ income classification. We constructed chord diagrams to provide a graphical representation of these patients' comorbid conditions and demographics differentiated by the income classification. A significance level of < 0.001 and a confidence level of 95% was chosen to determine statistical differences. This was selected as large datasets, such as the ISARIC COVID-19 dataset, might identify minor differences as significant even when the differences are not clinically relevant. Adjusting the rejection level of the null hypothesis could control this limitation inherent to large datasets and the possibility of incurring type one error. All data processing and statistical analysis were performed using Python version 3.8 with the following data packages: Pandas version 1.2.4, Tidyverse version 1.3.0, Bioconductor version 3.12.

Results

A total of 66,565 patients were included in this study (Fig. 1). Most patients were male (63.5% [42,256/66,565]) and treated in HICs (82.6% [55,004/66,565]). Specifically, 78.2% ([52,039/66,565]) of the cohort was hospitalised and treated in Europe. Regarding the age of the patients included in the cohort, 44.0% ([29,317/66,565]) of patients were between 60 and 80 years old. During the first 24 h of hospital admission, patients were most frequently treated with HFNC (42.5% [28,256/66,565]), followed by NIV (36.2% [24,112/66,565]) and IMV (21.3% [14,197/66,565]). Demographic characteristics, physiological variables and laboratories at hospital admission are shown in Fig. 2 and Tables 1 and 2.

Fig. 1

Flow chart. This figure shows patients included in the analysis and cohort selection process

Fig. 2

Probability density of patients' basic demographics (age and sex), according to the first ventilation treatment received. A Complete cohort. B Patients from high-income countries. C Patients from low middle-income countries

Table 1

Baseline characteristics of patients, stratified by the different advance ventilatory supports

Characteristic	All n = 66,565	HFNC n = 28,256	NIV n = 24,112	IMV n = 14,197	p value
Demographics, n (%)
Female	24,309 (36.5)	11,188 (39.6)	8600 (35.7)	4521 (31.8)	< 0.001
Age 0–20 years old	558 (0.9)	188 (0.7)	156 (0.6)	214 (1.5)	< 0.001
Age 20–40 years old	4888 (7.4)	1800 (6.3)	1635 (6.8)	1453 (10.2)	< 0.001
Age 40–60 years old	19,514 (29.3)	6907 (24.5)	7429 (30.8)	5178 (36.5)	< 0.001
Age 60–80 years old	29,317 (44.0)	11,350 (40.2)	11,505 (47.7)	6462 (45.5)	< 0.001
Age 80–100 years old	12,232 (18.4)	7962 (28.1)	3382 (14.1)	888 (6.2)	< 0.001
Age ≥ 100 years old	56 (0.1)	49 (0.2)	5 (0.0)	2 (0.0)	< 0.001
Pandemic wave in which patients were admitted, n (%)
First COVID-19 wave	27,044 (40.6)	13,363 (47.3)	7888 (32.7)	5793 (40.8)	< 0.001
Continent of admission, n (%)
Africa	89 (0.1)	17 (0.1)	3 (0.0)	69 (0.5)	< 0.001
Asia	10,488 (15.8)	1520 (5.4)	2590 (10.7)	6378 (44.9)	< 0.001
Europe	52,039 (78.2)	25,586 (90.6)	20,924 (86.8)	5529 (38.9)	< 0.001
North America	2434 (3.7)	561 (2.0)	277 (1.1)	1596 (11.2)	< 0.001
Oceania	260 (0.4)	155 (0.5)	8 (0.0)	97 (0.7)	< 0.001
South America	1255 (1.9)	417 (1.5)	310 (1.3)	528 (3.7)	< 0.001
Regional income stratification, n (%)
High-income country	55,004 (82.6)	26,399 (93.4)	21,237 (88.1)	7368 (51.9)	< 0.001
Low middle-income country	11,561 (17.4)	1857 (6.6)	2875 (11.9)	6829 (48.1)	< 0.001
Chronic comorbidities, n (%)
Asthma	8097 (12.2)	3596 (12.7)	3413 (14.2)	1088 (7.7)	< 0.001
Chronic cardiac disease (not hypertension)	14,678 (22.1)	7794 (27.6)	5153 (21.4)	1731 (12.2)	< 0.001
Chronic kidney disease	7533 (11.3)	4135 (14.6)	2571 (10.7)	827 (5.8)	< 0.001
Chronic neurological disorder	4944 (7.4)	2808 (9.9)	1560 (6.5)	576 (4.1)	< 0.001
Chronic pulmonary disease (not asthma)	8856 (13.3)	4459 (15.8)	3551 (14.7)	846 (6.0)	< 0.001
Dementia	3964 (6.0)	3032 (10.7)	818 (3.4)	114 (0.8)	< 0.001
Diabetes mellitus	20,164 (30.3)	8273 (29.3)	7343 (30.5)	4548 (32.0)	< 0.001
HIV	271 (0.4)	105 (0.4)	93 (0.4)	73 (0.5)	0.08
Arterial hypertension	27,521 (41.3)	11,855 (42.0)	9874 (41.0)	5792 (40.8)	0.02
Hypothyroidism	1632 (2.5)	864 (3.1)	598 (2.5)	170 (1.2)	< 0.001
Immunosuppression	1242 (1.9)	659 (2.3)	491 (2.0)	92 (0.6)	< 0.001
Malignant neoplasm	5115 (7.7)	2811 (9.9)	1803 (7.5)	501 (3.5)	< 0.001
Malnutrition	894 (1.3)	565 (2.0)	229 (0.9)	100 (0.7)	< 0.001
Mental disorder	1042 (1.6)	541 (1.9)	418 (1.7)	83 (0.6)	< 0.001
Moderate or severe liver disease	880 (1.3)	465 (1.6)	282 (1.2)	133 (0.9)	< 0.001
Obesity	10,793 (16.2)	3960 (14.0)	4883 (20.3)	1950 (13.7)	< 0.001
Rheumatological disorder	5412 (8.1)	3033 (10.7)	1989 (8.2)	390 (2.7)	< 0.001
Smoking	15,190 (22.8)	6948 (24.6)	6521 (27.0)	1721 (12.1)	< 0.001
Solid tumour	522 (0.8)	307 (1.1)	186 (0.8)	29 (0.2)	< 0.001
Complications, n (%)
Acute kidney injury	13,353 (20.1)	5146 (18.2)	4525 (18.8)	3682 (25.9)	< 0.001
Anaemia	10,031 (15.1)	3803 (13.5)	3492 (14.5)	2736 (19.3)	< 0.001
ARDS	15,470 (23.2)	4846 (17.2)	5625 (23.3)	4999 (35.2)	< 0.001
Bacteraemia	3966 (6.0)	1191 (4.2)	1381 (5.7)	1394 (9.8)	< 0.001
Cardiac arrest	3882 (5.8)	1215 (4.3)	1275 (5.3)	1392 (9.8)	< 0.001
Cardiac arrhythmia	5989 (9.0)	2070 (7.3)	2208 (9.2)	1711 (12.1)	< 0.001
Cardiac ischemia	1175 (1.8)	471 (1.7)	426 (1.8)	278 (2.0)	0.10
Coagulation disorder	3231 (4.9)	1122 (4.0)	1346 (5.6)	763 (5.4)	< 0.001
Congestive heart failure	2188 (3.3)	1159 (4.1)	749 (3.1)	280 (2.0)	< 0.001
Gastrointestinal bleeding	1130 (1.7)	519 (1.8)	310 (1.3)	301 (2.1)	< 0.001
Liver dysfunction	5600 (8.4)	1972 (7.0)	2176 (9.0)	1452 (10.2)	< 0.001
Neurological complication	1206 (1.8)	522 (1.8)	458 (1.9)	226 (1.6)	0.08
Pleural effusion	3967 (6.0)	1858 (6.6)	1285 (5.3)	824 (5.8)	< 0.001
Pneumothorax	1590 (2.4)	458 (1.6)	671 (2.8)	461 (3.2)	< 0.001
Pulmonary embolism	1951 (2.9)	667 (2.4)	869 (3.6)	415 (2.9)	< 0.001
Stroke	918 (1.4)	358 (1.3)	302 (1.3)	258 (1.8)	< 0.001
Treatments, n (%)
Prone	15,778 (23.7)	3384 (12.0)	6628 (27.5)	5766 (40.6)	< 0.001
Vasopressors/inotropes	13,592 (20.4)	2188 (7.7)	4282 (17.8)	7122 (50.2)	< 0.001
Corticoids	40,810 (61.3)	15,586 (55.2)	17,043 (70.7)	8181 (57.6)	< 0.001
Intensive care unit	36,336 (54.6)	8302 (29.4)	14,180 (58.8)	13,854 (97.6)	< 0.001
Clinical outcomes
Hospital discharge	33,627 (50.5)	16,302 (57.7)	12,115 (50.2)	5210 (36.7)	< 0.001
28-day fatality ratio	24,021 (36.1)	9581 (33.9)	9177 (38.1)	5263 (37.1)	< 0.001
Non-invasive ventilation failure (HFNC and NIV)	10,287 (15.5)	3538 (12.5)	6749 (28.0)		< 0.001

Bold values indicate statistical significance

HFNC high-flow nasal cannula, NIV non-invasive mechanical ventilation, IMV invasive mechanical ventilation, HIV human immunodeficiency virus, ARDS acute respiratory distress syndrome

Table 2

Physiological parameters and laboratories of patients during the first 24-h hospital admission, stratified by the different advance ventilatory supports

Measure	All		HFNC		NIV		IMV		p value
	Value	n	Value	n	Value	n	Value	n
Physiological parameters on admission, median (IQR)
Diastolic blood pressure (mmHg)	74.0 (65.0–83.0)	64,667	74.0 (65.0–83.0)	27,806	75.0 (66.0–83.0)	23,713	73.0 (64.0–82.0)	13,148	< 0.001
Heart rate (beats/min)	92.0 (80.0–106.0)	64,563	91.0 (80.0–105.0)	27,564	94.0 (82.0–107.0)	23,756	93.0 (81.0–108.0)	13,243	< 0.001
Respiratory rate (breaths/min)	24.0 (20.0–28.0)	63,891	22.0 (20.0–28.0)	27,631	24.0 (20.0–30.0)	23,439	24.0 (20.0–28.0)	12,821	< 0.001
Systolic blood pressure (mmHg)	129.0 (115.0–144.0)	64,745	129.0 (114.0–144.0)	27,828	130.0 (116.0–145.0)	23,740	129.0 (113.0–142.0)	13,177	< 0.001
Temperature (C)	37.2 (36.7–38.1)	64,519	37.2 (36.6–38.1)	27,763	37.3 (36.7–38.2)	23,685	37.0 (36.7–37.8)	13,071	< 0.001
Laboratory during the first 24 h, median (IQR)
Alanine aminotransferase (U/L)	32.0 (21.0–53.0)	28,901	29.0 (19.0–48.0)	13,612	34.0 (22.0–55.0)	10,956	38.0 (25.0–62.6)	4333	< 0.001
Aspartate aminotransferase (U/L)	50.0 (33.0–77.0)	6614	47.0 (31.0–73.0)	2413	51.0 (34.0–78.0)	2085	51.0 (34.0–81.0)	2116	< 0.001
Base excess (mmol/L)	0.0 (− 2.8–2.5)	7017	0.5 (− 2.1–2.9)	2542	0.3 (− 2.2–2.6)	2379	− 1.0 (− 4.0–1.9)	2096	< 0.001
Bicarbonate (mEq/L)	23.0 (20.4–25.8)	9601	23.5 (20.9–25.9)	2254	23.1 (20.5–26.0)	2536	23.0 (20.0–25.7)	4811	< 0.001
Bilirubin (umol/L)	10.0 (7.0–14.0)	29,034	9.0 (7.0–14.0)	13,705	10.0 (7.0–14.0)	11,029	9.0 (6.0–15.0)	4300	< 0.001
C reactive protein (mg/L)	106.0 (53.0–179.0)	34,518	94.4 (48.0–161.55)	16,762	118.0 (63.0–191.0)	13,899	118.5 (45.63–213.7)	3857	< 0.001
Creatine kinase (U/L)	163.0 (76.0–427.75)	4702	151.5 (67.75–409.0)	1880	164.0 (80.0–431.0)	1557	179.0 (83.0–450.0)	1265	0.001
Creatinine (umol/L)	87.0 (68.95–120.0)	42,151	87.0 (69.0–120.0)	17,420	87.0 (69.0–116.0)	14,859	87.52 (67.19–126.41)	9872	0.20
D-Dimer (mg/L)	2.76 (0.9–375.0)	1395	1.54 (0.75–9.82)	357	6.94 (0.89–607.0)	303	3.68 (1.02–453.0)	735	< 0.001
Ferritin (ug/L)	810.0 (384.0–1523.0)	6244	735.0 (326.0–1425.7)	2542	838.0 (411.92–1535.5)	2460	942.0 (447.0–1672.75)	1242	< 0.001
Glucose (mmol/L)	7.55 (6.2–10.4)	33,299	7.0 (5.9–9.3)	12,720	7.6 (6.3–10.5)	12,543	8.49 (6.83–11.84)	8036	< 0.001
Haematocrit (%)	39.0 (34.6–42.3)	11,810	39.7 (35.1–43.0)	2389	39.0 (35.0–42.8)	3064	38.6 (34.0–42.0)	6357	< 0.001
Haemoglobin (g/L)	133.0 (117.0–146.0)	50,242	133.0 (118.0–146.0)	20,731	135.0 (121.0–148.0)	19,009	127.0 (109.0–142.0)	10,502	< 0.001
Interleukin 6 (ng/L)	67.6 (23.0–169.0)	433	43.51 (13.33–89.3)	107	49.4 (21.48–129.5)	128	131.3 (31.1–313.0)	198	< 0.001
Lactate dehydrogenase (U/L)	487.0 (349.0–684.0)	7570	445.5 (328.0–621.5)	2966	532.0 (374.0–741.0)	2837	496.0 (353.0–708.5)	1767	< 0.001
Lactic acid (mmol/L)	1.5 (1.1–2.1)	21,382	1.4 (1.05–2.0)	9154	1.5 (1.1–2.04)	8553	1.55 (1.1–2.3)	3675	< 0.001
Leukocytes (109/L)	8.07 (5.7–12.0)	50,673	7.4 (5.4–10.5)	21,036	7.8 (5.6–11.3)	19,043	10.77 (7.0–18.2)	10,594	< 0.001
Lymphocytes (109/L)	0.8 (0.58–1.2)	36,068	0.81 (0.59–1.2)	18,019	0.8 (0.58–1.11)	14,760	0.82 (0.55–1.3)	3289	< 0.001
Lymphocytes/leukocytes (%)	9.7 (5.35–15.65)	879	11.0 (6.8–16.1)	357	8.95 (5.0–16.58)	130	8.6 (4.57–15.0)	392	< 0.001
Neutrophils (109/L)	5.8 (4.0–8.63)	35,963	5.6 (3.87–8.31)	18,020	5.8 (4.0–8.43)	14,719	7.5 (4.7–11.5)	3224	< 0.001
Neutrophils/leukocytes (%)	82.0 (72.9–88.0)	697	81.9 (74.4–87.3)	293	80.4 (70.0–87.7)	101	83.2 (73.2–88.95)	303	0.26
Platelets (109/L)	199.0 (140.0–265.0)	50,263	207.0 (157.0–271.0)	20,811	202.0 (148.0–265.0)	18,903	162.0 (0.3–249.0)	10,549	< 0.001
Potassium (mmol/L)	4.1 (3.8–4.5)	35,901	4.1 (3.74–4.5)	14,512	4.1 (3.8–4.5)	12,433	4.2 (3.8–4.6)	8956	< 0.001
Procalcitonin (ug/L)	0.24 (0.12–0.7)	6234	0.2 (0.1–0.51)	2191	0.24 (0.12–0.62)	2839	0.4 (0.15–1.38)	1204	< 0.001
Prothrombin intl. (ratio)	1.1 (1.02–1.3)	2728	1.09 (1.0–1.2)	611	1.1 (1.03–1.3)	604	1.14 (1.04–1.3)	1513	< 0.001
Prothrombin time (s)	13.0 (11.3–14.5)	25,413	12.8 (11.1–14.4)	11,834	13.0 (11.4–14.5)	10,400	13.3 (11.8–14.8)	3179	< 0.001
Sodium (mmol/L)	136.0 (133.0–140.0)	38,268	137.0 (134.0–140.0)	15,894	136.0 (133.0–139.0)	13,253	137.0 (133.0–141.0)	9121	< 0.001
Troponin I (ug/L)	0.07 (0.02–6.9)	1437	0.03 (0.01–0.25)	398	0.08 (0.02–10.3)	257	0.13 (0.02–10.0)	782	< 0.001
Urea nitrogen (mmol/L)	7.7 (5.1–12.85)	46,588	7.1 (4.9–11.5)	19,388	7.3 (5.0–11.6)	17,614	10.35 (6.2–18.56)	9586	< 0.001

Bold values indicate statistical significance

HFNC high-flow nasal cannula, NIV non-invasive mechanical ventilation, IMV invasive mechanical ventilation, IQR interquartile range, PTT partial thromboplastin time

Patients’ characteristics, in-hospital treatments, and systemic complications

More than 85% of the patients had at least one comorbidity. Hypertension (41.3% [27,521/66,565]) and diabetes mellitus (30.3% [20,164/66,565]) were the most frequently reported comorbid conditions (Table 1). A total of 22.8% [15,190/66,565] of patients were current or past smokers. Complications were also common during the hospital admission (not at hospital presentation), 23.2% [15,470/66,565] developed ARDS, and 20.1% [13,353/66,565] were reported to have an acute renal injury (ARI).

During hospital admission, 61.3% (40,810/66,565) patients received corticosteroid treatment, and 54.6% [36,336/66,565] were admitted to the ICU. Vasopressor/inotrope therapy was used in a quarter of all patients (20.4% [13,592/66,565]), increasing in use according to ventilatory requirement (7.7% [2188/28,256] vs. 17.8% [4282/24,112] vs 50.2 [7122/14,197]). Approximately half of those treated with IMV received vasopressors/inotropes at some point during hospitalisation (50.2% [7122/14,197]). Almost one-quarter of the patients were placed in prone position (23.7% [15,778/66,565]), more commonly in those patients treated with IMV (12.0% [3384/28,256] vs. 27.5% [6628/24,112] vs 40.6% [5766/14,197]). A total of 15.5% [10,287/66,565] of patients failed HFNC/NIV. Moreover, 71.2% [7327/10,287] of the patients that failed HFNC/NIV were registered in HIC and 28.8% [2960/10,287] in LMIC. Finally, 28-day mortality was similar between the different advance ventilatory supports (33.9% [9581/28,256] vs. 38.1% [9177/24,112] vs. 37.1% [5263/14,197]).

Comparing respiratory support of patients admitted in HIC or LIMC

The cumulative frequency of advanced respiratory treatments was stratified by national income classification (Fig. 3). Patients admitted to the hospital in HICs were more frequently treated with HFNC (48.0% [26,399/55,004]), followed by NIV (38.6% [21,237/55,004]) and IMV (13.4% [7368/55,004]). In contrast, patients admitted in LMICs were less frequently treated with HFNC (16.1% [1857/11,561]), and the majority received IMV (59.1% [6829/11,561]) (Table 1; Fig. 3). We also found differences in distribution among the different types of respiratory support when stratified by income classification and respiratory support (Fig. 4).

Fig. 3

Cumulative frequency (net number of patients) of ventilation treatment given to patients. A Complete cohort. B Patients from high-income countries. C Patients from low middle-income countries

Fig. 4

Chord graphic with demographics and comorbidities of patients according to the type of first ventilation treatment received. A Complete cohort. B Patients from high-income countries. C Patients from low middle-income countries

Patients treated with IMV in HICs had fewer comorbid conditions and were more frequently between 40 and 70 years old. In sharp contrast, patients in LMIC who were younger than 40 years old often received IMV and were more frequently male. Also, they were mostly treated with IMV rather than non-invasive respiratory strategies (Fig. 4).

Changes in respiratory supports

Figure 5 presents the alluvia diagrams illustrating how patients progressed among respiratory support during hospital admission. Notably, patients who required more than one respiratory treatment had higher mortality than those treated with only one type of support, whether the first respiratory support was HFNC, NIV, or IMV (Fig. 5).

Fig. 5

Alluvia diagram of the patients’ transitions between ventilation treatments and clinical outcomes. The width of the links is proportional to the number of patients. A Complete cohort. B Patients from high-income countries. C Patients from low middle-income countries

Risk factors for failing HFNC or NIV as first respiratory support

The failure rate of HFNC or NIV was 15.5% [10,287/66,565]. According to the Gini importance, the variables most strongly associated with non-invasive ventilation failure (either HFNC or NIV) were age, lower platelets, and higher leukocyte count during the first 24 h of hospital admission (Fig. 6A). In the logistic regression analysis, we found that high leukocyte counts at hospital admission (OR [95% CI]; 5.86 [4.83–7.10]), treatment in an LMIC (OR [95% CI]; 2.04 [1.97–2.11]), and tachypnoea at hospital admission (OR [95% CI]; 1.16 [1.14–1.18]) were strongly associated factors with IMV treatment as rescue treatment (Fig. 6B, C).

Fig. 6

An automatised model to determine risk factors associated with non-invasive ventilation failure. A variables more strongly related to non-invasive ventilation failure according to the Gini importance. B the contribution of the variables to the output; the red values indicate a high-value contribution of the variable, and the blue values a low-value contribution. The positive values in the plot indicate a high probability of 28-day fatality, and negative values indicate a low likelihood of 28-day fatality. Panel C presents a logistic regression model, showing variables more strongly associated with the 28-day fatality ratio. The most significant variables were leucocyte count, low-/middle-income country attention, higher respiratory rate, and higher systolic blood pressure

Clinical outcomes and risk factors associated with 28-day fatality ratio

Almost half of the patients treated with HFNC [46.3%; 11,954/28,256] and 37.1% (5263/14,197) of patients treated with IMV died within 28 days. The variables identified as risk factors associated with the 28-day fatality ratio are shown in Fig. 7. Older age (OR [95% CI]; 2.42 [2.36–2.48]), cardiac arrest during hospitalisation (OR [95% CI]; 1.86 [1.81–1.92]), receiving treatment in an LMIC (OR [95% CI]; 1.56 [1.53–1.60]), and higher leukocyte counts at hospital admission (OR [95% CI]; 1.47 [1.39–1.55]) were the main adjusted risk factors associated with 28-day mortality. Notably, NIV/HFNC failure (OR [95% CI]; 1.27 [1.25–1.30]) was also highly associated with fatality. Other factors were acute kidney injury (OR [95% CI]; 1.23 [1.21–1.25]), ARDS (diagnosed during the hospital admission, not during the first 24 h) (OR [95% CI]; 1.12 [1.10–1.14]), increased heart rate at admission (OR [95% CI]; 1.15 [1.13–1.18]), increased respiratory rate at admission (OR [95% CI]; 1.15 [1.13–1.17]), chronic cardiac diseases (OR [95% CI]; 1.17 [1.14–1.19]), chronic pulmonary diseases (OR [95% CI]; 1.12 [1.10–1.14]), and diabetes mellitus (OR [95% CI]; 1.07 [1.05–1.09]). The model used to predict the 28-day fatality ratio had a good discriminatory capacity when evaluated by the AUROC (mean [SD] 0.78 [0.05], Fig. 7).

Fig. 7

An automatised model to determine risk factors associated with the 28-day fatality ratio. A The contribution of the variables to the output; the red values indicate a high-value contribution of the variable, and the blue values a low-value contribution. The positive values in the plot indicate a high probability of 28 fatalities, and negative values indicate a low likelihood of 28-day fatality. B A logistic regression model, showing variables more strongly associated with the 28-day fatality ratio. C Each cross-validation trial's receiver operative curve (ROC) for the subset of the selected variables. The blue curve represents the average of the ROC curves of each test, and the average area under the ROC is also presented. The most significant variables associated with the 28-day fatality ratio were age, cardiac arrest, low-/middle-income country attention, and leucocyte count. Also, patients that fail the non-invasive or high-flow nasal cannula are independently associated with a higher 28-day fatality ratio

Discussion

In this large, multinational, prospective cohort study, we found that patients with severe COVID-19 were mainly treated with non-invasive respiratory strategies (i.e. HNFC or NIV) in HICs; in contrast, patients with severe COVID-19 in LMICs were more frequently treated with IMV. We found that the 28-day fatality ratio was similar among patients treated with HFNC, NIV, or IMV worldwide. Notably, we found that patients treated with IMV as rescue therapy (i.e. failure of non-invasive treatments) had a higher 28-day fatality ratio than patients treated with IMV earlier in their disease course. The risk factors associated with failing the non-invasive respiratory strategies were high leukocyte counts at admission, increased heart rate at admission, and being treated in an LMIC. Notably, being admitted during the first pandemic wave did not impact clinical outcomes or respiratory treatments.

Early in the pandemic, healthcare workers identified that patients with hypoxemia could be treated with HFNC [21-23]. International guidelines also recommend non-invasive respiratory support as the first treatment, and many centres utilise HFNCs outside formal ICU settings [13]. Notably, the widespread usage of HFNC and NIV in patients with severe COVID-19 was recommended by experts and guidelines but not supported by high-quality data. Later, Ospina-Tascon et al. [12] carried out a multicentre, open randomised clinical trial and found that the early treatment with HFNC compared to conventional oxygen treatment was associated with a lower necessity of IMV (34.3 Vs 51.0, HR: 1.39; 95% CI 1.00–1.92; p = 0.04). Then, Perkins et al. [24] in the RECOVERY-RS trial found that NIV was associated with a lower requirement of tracheal intubation and lower 30-day mortality when compared to conventional oxygen therapy (absolute difference, − 8% [95% CI, − 15 to − 1%], p = 0.03). Our study found that HFNC, NIV, and IMV have similar 28-day fatality ratios, in concordance with prior literature. However, we found that HFNC was mainly used in HIC, which might be in relation to the capacity of these countries to acquire this new technology during the pandemic and the ability of these countries to expand their bed capacity to treat patients with HFNC outside of the ICU. Also, some patients or their families do not accept endotracheal intubation and prefer non-invasive strategies, though our study did not collect these data.

In contrast to HIC, the most common respiratory treatment in patients with severe COVID-19 utilised in LMIC has been IMV, as is evident in our data. Estenssoro et al. [17] described the results of a prospective observational cohort of patients admitted to 64 ICUs in Argentina. They included 1909 patients treated with IMV and found that lung-protective respiratory strategies were widely used but with a high fatality rate among patients included in the cohort (57.7%, 1101/1909). In another study in Brazil, Ranzani et al. [7] found that 23% (45,205/232,036) of patients admitted to the hospital were treated with IMV. They also found that the fatality rate among those receiving IMV was 80% during the first pandemic wave and 87% during the second wave [7, 25]. Notably, they found that 14% (5976/44,055) of the patients treated with IMV were treated outside of the ICU [25]. These results highlight that fatality rates and treatments changed during the pandemic and differed for each country. Moreover, these data align with our results, showing that IMV was frequently used in LMIC and that many patients with severe COVID-19 were treated outside of ICU [8, 26]. Notably, the impact of ICU admission on clinical outcomes was already explored in our cohort and published elsewhere [3]. We found that ICU admission was associated with better clinical outcomes independently of disease severity, treatments received, income classification, and system saturation (i.e. the number of new COVID-19 detected the day patients was admitted).

Even though non-invasive respiratory support has been proven effective in treating patients with severe hypoxemia during COVID-19, up to 30% of the patients were treated with IMV as a rescue treatment. Thus, it is essential to identify which patients might be at risk of failing under the non-invasive respiratory strategy and not to delay IMV in these patients. Rodriguez A. et al., in one of the largest prospective cohorts of patients admitted to the ICU due to severe flu infection, found that patients who failed NIV had a mortality rate three times higher than those who did not fail [27]. Also, they found that patients who failed NIV had higher mortality than those treated with IMV as initial treatment (38.4 vs 31.3, p = 0.18). In a multicentre COVID-19 study, Boscolo A. et al. found that 704 patients who failed non-invasive respiratory support had an accumulative fatality rate of 43% [28]. Our findings support that patients with severe COVID-19 who fail the initial respiratory support with non-invasive treatments have a higher mortality rate and were independently associated with 28-day fatality. Also, we found that patients with higher leukocyte counts at admission, higher respiratory rate at admission, and being in an LMIC were at higher risk of failing the non-invasive respiratory strategies. Thus, patients with these characteristics should be carefully evaluated to avoid delays in initiating IMV when appropriate.

Our study has strengths and limitations that are important to acknowledge. First, the respiratory support interventions were not according to a standardised protocol, leaving clinical teams to choose when to use HFNC, NIV, or IMV; thus, demographic or clinical characteristics may differ across the groups studied. However, we performed a robust statistical analysis using random forest analyses and logistic regression, adjusting for several confounders. This allowed us to evaluate linear and nonlinear relations in a supervised statistical approach. Second, most patients in our study were registered in Europe and HICs, which might constitute a significant selection bias. However, we had more than 11,000 patients in LMICs in Africa, South America, and Asia, including a large cohort of patients and contributing to our results' global generalisability. Third, we do not have complete data on specific respiratory parameters used during the support (i.e. peep, flows, FiO2, volumes, among many others), limiting our capacity to assess the rates of protective respiratory strategies, among other essential factors. Thus, these results cannot imply a causal association between respiratory support device treatments and clinical outcomes. Each patient should be evaluated carefully with decisions on the type of respiratory support based upon the evolving evidence base applied to their specific clinical condition and goals of care. Finally, throughout the COVID-19 pandemic, patients were treated with a large variety of medications and supportive clinical protocols; it is challenging to make conclusions about the factors associated with 28-day fatality using observational study methodologies in such a dynamic context.

Conclusions

Patients hospitalised with confirmed COVID-19 are often treated with advanced respiratory support. HFNC was the primary initial respiratory support used during the pandemic; however, this treatment was mainly used in HIC. In contrast, IMV was the primary respiratory treatment utilised in LMIC. Non-invasive respiratory treatments (i.e. HFNC and NIV) could be used as the first respiratory support in patients with severe COVID-19; however, it is crucial to identify patients at risk of failing because delaying IMV may be associated with worse clinical outcomes. Further studies are needed to confirm these associations.

Additional file 1. Supplemental methods.

Additional file 2. Conflict of interests.