Skip to main content

Subphenotypes in patients with acute respiratory distress syndrome treated with high-flow oxygen



Acute respiratory distress syndrome (ARDS) subphenotypes differ in outcomes and treatment responses. Subphenotypes in high-flow nasal oxygen (HFNO)-treated ARDS patients have not been investigated.


To identify biological subphenotypes in HFNO-treated ARDS patients.


Secondary analysis of a prospective multicenter observational study including ARDS patients supported with HFNO. Plasma inflammation markers (interleukin [IL]-6, IL-8, and IL-33 and soluble suppression of tumorigenicity-2 [sST2]) and lung epithelial (receptor for advanced glycation end products [RAGE] and surfactant protein D [SP-D]) and endothelial (angiopoietin-2 [Ang-2]) injury were measured. These biomarkers and bicarbonate were used in K-means cluster analysis to identify subphenotypes. Logistic regression was performed on biomarker combinations to predict clustering. We chose the model with the best AUROC and the lowest number of variables. This model was used to describe the HAIS (High-flow ARDS Inflammatory Subphenotype) score.


Among 41 HFNO patients, two subphenotypes were identified. Hyperinflammatory subphenotype (n = 17) showed higher biomarker levels than hypoinflammatory (n = 24). Despite similar baseline characteristics, the hyperinflammatory subphenotype had higher 60-day mortality (47 vs 8.3% p = 0.014) and longer ICU length of stay (22.0 days [18.0–30.0] vs 39.5 [25.5–60.0], p = 0.034). The HAIS score, based on IL-8 and sST2, accurately distinguished subphenotypes (AUROC 0.96 [95%CI: 0.90–1.00]). A HAIS score ≥ 7.45 was predictor of hyperinflammatory subphenotype.


ARDS patients treated with HFNO exhibit two biological subphenotypes that have similar clinical characteristics, but hyperinflammatory patients have worse outcomes. The HAIS score may identify patients with hyperinflammatory subphenotype and might be used for enrichment strategies in future clinical trials.


Acute respiratory distress syndrome (ARDS) is a common condition in ICU, carrying a high mortality burden. Recent research and the latest guidelines on ARDS management underlined the importance of identifying clinical and biological features to classify ARDS patients into subphenotypes that might have different outcomes and respond differently to specific therapies [1]. Identification of subphenotypes in ARDS mechanically ventilated patients has mainly relied on inflammation-related biomarkers. Consistently across studies, two subphenotypes have been reported, termed hypoinflammatory and hyperinflammatory, that have different outcomes and distinct responses to clinical interventions [2,3,4,5,6,7]. The stability of the subphenotypes over time was also shown in the same individual [8]. More recently, these subphenotypes were also identified in a large cohort of patients with COVID-19-associated ARDS [9]. Differentiating hypoinflammatory and hyperinflammatory ARDS enables prognosis enrichment and enhances the likelihood of treatment response [10].

Since the use of high-flow nasal oxygen (HFNO) in adult patients with hypoxemic acute respiratory failure was first described [11], its use has continuously increased [12,13,14]. Traditionally, patients treated with HFNO with bilateral opacities in the chest X-ray and a PaO2/FiO2 ratio ≤ 300 mmHg were not considered ARDS patients as they were not on positive pressure ventilation. However, some studies including non-intubated patients who meet all other ARDS criteria showed that they had similar biological characteristics [15] and clinical outcomes [16] to mechanically ventilated ARDS patients. Similarly, the results of a recent study in COVID-19 patients showed that more than 90% of the patients initially treated with HFNO remained with a PaO2/FiO2 ratio ≤ 300 after intubation [17]. All this evidence led to the addition of HFNO in a recently published more global definition of ARDS [18]. The inclusion of HFNO patients would allow us to identify patients with ARDS in earlier stages of lung injury and, consequently, to test treatments prior to intubation and mechanical ventilation (MV). It will also broaden the applicability of the ARDS definition by focusing on all patients with clinically significant lung injury requiring high levels of oxygen support and including those patients treated in resource-limited settings.

Patients treated with HFNO may present better outcomes than MV ARDS patients [17, 19]. These differences in outcomes may be, at least partially, explained by the presence of different subphenotypes. However, the characterization of subphenotypes in patients with ARDS treated with HFNO has not been explored. We hypothesized that, as happen in MV ARDS patients, ARDS patients treated with HFNO might present different subphenotypes with different outcomes. Therefore, we aimed to identify inflammatory subphenotypes in a prospective cohort of ARDS patients treated with HFNO.


We conducted a post hoc analysis of a multicenter prospective observational study comparing plasma biomarkers in mechanically ventilated and non-intubated ARDS patients [15]. Clinical data from the original cohort have been published elsewhere [15]. Briefly, patients were enrolled between 2014 and 2016 across 3 general ICU in tertiary hospitals in Spain.

In the present study inclusion criteria were non-intubated patients, treated with HFNO for hypoxemia defined by PaO2/FiO2 ≤ 300 or pulse oximetry (SpO2)/FiO2 ≤ 315, with bilateral radiographic opacities not fully explained by cardiac failure. FiO2 was set to ensure a SpO2 > 92%. The flow was adjusted according to patient tolerance. Need for MV was left at the discretion of the attending physician. Ethics Committee of each hospital approved the study. Written consent from patients or their relatives was obtained.

Baseline-recorded data included demographic characteristics, comorbidities, and etiology of ARDS. General respiratory and hemodynamic variables were collected at inclusion. Severity was assessed with the Acute Physiology and Chronic Health Evaluation (APACHE) II score within 24 h of ICU admission.

Sequential Organ Failure Assessment (SOFA) at ARDS onset while treated with HFNO was calculated. Acute renal failure was defined as a serum level of creatinine of 1.2 mg/dL or higher, and shock was defined by use of vasopressors. Follow-up for survival was performed during the 60 days after inclusion.

Blood samples were collected within 24 h of ARDS onset while being treated with HFNO. Plasma biomarkers of lung epithelial (receptor for advanced glycation end products [RAGE] and surfactant protein D [SP-D]) and endothelial (angiopoietin-2 [Ang-2]) injury as well as inflammation markers (interleukin [IL]-6, IL-8, and IL-33 and soluble suppression of tumorigenicity-2 [sST2]) levels were measured using enzyme-linked immunosorbent assay kits.

A detailed description of the statistical analysis is provided in the additional material (Additional file 1). Briefly, missing data were imputed using multiple imputation. sST2, IL-6, IL-8, angiopoietin-2, and s-RAGE were log10. Cluster analysis by K-means was used on the following variables: IL-33, sST2 (log), IL-6 (log), IL-8 (log), SP-D, angiopoietin-2 (log), RAGE (log), and bicarbonate level. A heatmap was used to represent the contribution of each biomarker to cluster identification. Values range from −1 to 1. A value close to zero has little or no impact on cluster attribution, while being close to 1 or −1 means a high impact on cluster attribution. The sign of the value indicates how biomarker’s variation influences cluster attribution (positive meaning that an increase in the value is correlated with cluster attribution). NbClust package was used to determine the optimal number of cluster [20]. Stability index was calculated by bootstrapping. After cluster attribution, patients were compared for initial characteristics and outcomes. Continuous variables were expressed as median [inter-quartile] and compared with Wilcoxon test. Categorical variables were reported as frequency (percentage) and compared with Chi-square.

To reduce the number of dimensions of our model, logistic regression was performed on all possible variable combinations to predict cluster belonging. We chose the model with the best area under the curve (AUC) of the receiving operative characteristics (ROC) and the lowest number of variables. A score to predict cluster belonging was designed by multiplying the variable's value by its corresponding odds ratio in the chosen model. The score threshold was then determined to achieve the optimal sensitivity and specificity [21].

Statistical analyses were conducted using R software (version 4.2.2). Results were considered statistically significant for a p value < 0.05.


Study population

During a 3-year period, 170 ARDS patients were enrolled in the original cohort [15]. Among those patients, 127 were under mechanical ventilation at inclusion and, therefore, they were not included in the present study. Forty-three patients were initially treated with HFNO, of whom 41 had measurement of biomarkers at inclusion. The baseline characteristics of the non-intubated ARDS population treated with HFNO included are presented in Table 1.

Table 1 Baseline characteristics

Identification of subphenotypes

The correlation matrix did not identify any couple of variables with a correlation coefficient R ≥ 0,50. The optimal number of clusters was defined as two. Using unsupervised clustering analysis and considering the plasma concentration of the previously mentioned biomarkers, we identified 2 clusters of patients within the ARDS patients treated with HFNO (see Additional file 2). The two clusters were well separated with no overlap (see Additional file 3). The contribution of each biomarker to cluster identification was depicted using a heatmap (Fig. 1). IL8 and sST2 are the two biomarkers that differ more between clusters. The stability of our model was evaluated by calculating the stability index by bootstrapping 200 times the samples. The stability index achieved a value of 0.63.

Fig. 1
figure 1

Heat map representing the strength of each variable to influence cluster belonging. Values range from −1 to 1. A value close to zero has little or no impact on cluster attribution, whereas a value close to 1 or −1 indicates a high impact on cluster attribution. The sign of the value indicates how the biomarker variation influences cluster attribution (when positive, an increase in the value is correlated with cluster attribution). IL: interleukin; sST2: soluble suppression of tumorigenicity-2; SP-D: surfactant protein D; RAGE: receptor for advanced glycation end products; ANG-2: angiopoietin; HCO3: bicarbonate

Clinical and biological characteristics and outcomes of the two subphenotypes

We defined cluster 2 as the hyperinflammatory subphenotype as they had significantly increased levels of IL-33, sST2, IL-8, SP-D, and ANG-2 compared to cluster 1 (Fig. 2 and see Additional file 4).

Fig. 2
figure 2

Comparison of biomarkers levels between clusters. IL: interleukin; sST2: soluble suppression of tumorigenicity-2; SP-D: surfactant protein D; RAGE: receptor for advanced glycation end products

No differences were found in IL-6, RAGE, and bicarbonate (HCO3). Cluster 1 was consequently defined as the hypoinflammatory subphenotype. No differences in baseline clinical characteristics were observed between both subphenotypes (Table 1).

The need for MV in the hyperinflammatory subphenotype was 64.7 versus 33.3% in the hypoinflammatory (p = 0.096). Patients belonging to the hyperinflammatory subphenotype were more frequently intubated at day 7 compared to the hypoinflammatory subphenotype (7 (30.4%) versus 9 (52.9%), p = 0.049). The ICU length of stay was longer in the hyperinflammatory subphenotype (22.0 days [18.0, 30.0] versus 39.5 [25.5, 60.0], p = 0.034). The hyperinflammatory patients had higher 60-day mortality compared to hypoinflammatory patients (47.1% vs. 9.7%; p = 0.020) (Table 2).

Table 2 Clinical outcomes during ICU course

Subphenotype prediction with reduced model

The ROC curves illustrating the capacity of each biomarker to predict cluster appurtenance are presented in Fig. 3A. IL-8 alone had an area under the ROC curve of 0.92 (0.83–1.00) for cluster prediction. Using a dimension reduction model by logistic regression analysis, IL-8 and sST2 were able to accurately predict the two subphenotypes (area under the ROC curve of 0.96 [95%-CI: 0.90–1.00]) (Fig. 3B).

Fig. 3
figure 3

Diagnostic accuracy for subphenotype identification of different biomarkers. A ROC curves describing the ability of each variable to predict cluster appurtenance. B ROC curve of the IL-8 + sST2 model

A simple score based on the plasma concentration of IL-8 and sST2, named HAIS score (High-flow ARDS Inflammatory Subphenotype), was able to determine the subphenotype of each patient (HAIS score = 1.48* log10([IL-8]) + 1.52* log10([ST-2])). A HAIS score ≥ 7.45 indicates that the patient belongs to subphenotype 2 with a sensitivity of 0.87 and a specificity of 0.94. It has a positive and negative predictive value of 0.95 and 0.84, respectively.


This is the first study that analyzed the presence of subphenotypes in non-intubated ARDS patients treated with HFNO. Based on inflammatory biomarkers, we identified two different subphenotypes that, despite presenting similar clinical characteristics at baseline, had divergent outcomes. Indeed, patients of hyperinflammatory subphenotype had higher mortality and longer ICU length of stay. Moreover, a score based on 2 inflammatory biomarkers (sT2 and IL-8) can accurately identify the subphenotype of HFNO patients who met ARDS criteria.

ARDS is characterized by a pulmonary and systemic inflammatory response which injures lung epithelium and endothelium resulting in protein-rich pulmonary edema [22]. Traditionally, the use of biomarkers has provided valuable insights into ARDS pathophysiology and has been used for risk stratification. Higher plasma concentrations of RAGE and SP-D, which are markers of epithelial injury, have been observed in ARDS patients [23, 24]. Similarly, elevated plasma Ang-2 predicts ARDS development in MV patients [25]. More recent studies have also shown the importance of the IL33/ST2 axis in the generation and modulation of lung injury [26, 27]. Finally, it has also been shown how these and other biomarkers, such as IL6 or IL8, are associated with outcomes, in both MV [28] and HFNO [15] ARDS patients. Interestingly, HFNO and MV ARDS patients have similar plasma levels of inflammatory biomarkers, suggesting that both have the same underlying pathophysiological alterations [15].

Subphenotypes in ARDS patients were first described by Calfee et al. [2] using latent class analysis of two retrospective cohorts of ARDS patients previously included in two randomized controlled trials. Subsequently, the same subphenotypes were observed in other secondary analyses of randomized controlled trials [3,4,5] and also in unselected populations of ARDS patients [6, 7, 9]. All studies agreed on the description of two subphenotypes that have different outcomes. As differences between subphenotypes were mainly driven by plasma levels of different biomarkers, they were named as the hyperinflammatory subphenotype, which had a higher mortality, and the hypoinflammatory subphenotype, which had a lower mortality rate. The subphenotype was typically evaluated on day one of the ARDS diagnosis. However, it has also been shown that ARDS subphenotypes are stable over the first three days of enrollment in the trial [8]. Finally, subphenotypes responded differently to various therapeutic interventions. In this sense, patients with the hyperinflammatory subphenotype have a lower mortality rate when treated with higher PEEP [2], a conservative fluid strategy [3], or simvastatin [4], suggesting that phenotyping ARDS patients may be a useful enrichment strategy to increase the likelihood of positive results in clinical trials with MV ARDS patients.

In contrast, most studies on HFNO have focused on identifying predictors of intubation based on patients’ clinical characteristics. Indeed, no studies have attempted to identify HFNO patients with a higher risk of death. To our knowledge, it is the first time that subphenotypes are identified across non-intubated ARDS patients. Indeed, another prospective observational study (NCT04009330) is currently recruiting patients to study ARDS subphenotypes in classical MV ARDS and in non-intubated ARDS. Importantly, this is the first time that we have been able to predict mortality in patients with HFNO. It should also be noted that our results are consistent with the subphenotypes observed in MV ARDS and their association with mortality. The identification of subphenotypes may therefore enhance enrichment strategies in clinical trials involving HFNO patients aiming to reduce the heterogeneity of treatment effect and to increase the likelihood of benefit of any given treatment.

The present study has some limitations. First, it is a secondary analysis that included a relatively small sample size. However, for k-means clustering analysis, the power to detect clustering primarily depends on cluster separation rather than on sample size [29]. Indeed, 20 observations per subgroup resulted in sufficient power to detect the presence of subgroups with k-means, also providing near-perfect accuracy for detecting the true number of clusters, and very high classification accuracy of individual observation’s group membership. Moreover, it is a unique cohort of HFNO patients in whom a panel of biomarker concentration has been determined. Second, our results need to be prospectively validated in an external cohort and the stability over time of these subphenotypes in HFNO patients should also be assessed. However, the results presented are consistent with previously reported data on MV ARDS patients [2,3,4,5,6,7, 9]. Third, the number of biomarkers compared to the number of patients exposed to the risk of overfitting. Nevertheless, this risk is compensated by the fact that each biomarker we used has been highly validated in MV ARDS studies, giving information on prognosis, or allowing the identification of subphenotypes. Fourth, the post hoc design of our study does not allow to study the correlation between biological subphenotypes and other clinical features such as radiological phenotype (diffuse or focal) or the delay between the first respiratory symptoms and non-intubated ARDS diagnosis. Likewise, as SpO2 was not recorded, we could not calculate the ROX index. However, PaO2/FiO2 ratio and RR were similar between groups. Therefore, it is unlikely that significant differences existed in the ROX index between groups.

In conclusion, using cluster analysis on inflammatory biomarkers we identified two subphenotypes in non-intubated ARDS patients under HFNO. These subphenotypes were indistinguishable according to the baseline clinical characteristics. However, they had radically different outcomes. A score based on sST2 and IL-8 serum concentrations can identify these subphenotypes with excellent accuracy and might be used for enrichment strategies in future clinical trials involving HFNO patients.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.





Acute Physiology and Chronic Health Evaluation


Acute respiratory distress syndrome


Area under the curve of the receiving operative characteristics

HAIS score:

High-flow ARDS Inflammatory Subphenotype score

HCO3 :



High-flow nasal oxygen


Intensive care unit








Mechanical ventilation

PaO2/FiO2 :

Partial pressure of oxygen in arterial blood to the fraction of inspiratory oxygen concentration ratio


Positive end-expiratory pressure


Receptor for advanced glycation end products


Sequential Organ Failure Assessment


Surfactant protein D


Peripheral oxygen saturation


Soluble suppression of tumorigenicity-2


  1. Grasselli G, Calfee CS, Camporota L, Poole D, Amato MBP, Antonelli M, Arabi YM, Baroncelli F, Beitler JR, Bellani G, Bellingan G, Blackwood B, Bos LDJ, Brochard L, Brodie D, Burns KEA, Combes A, D'Arrigo S, De Backer D, Demoule A, Einav S, Fan E, Ferguson ND, Frat JP, Gattinoni L, Guérin C, Herridge MS, Hodgson C, Hough CL, Jaber S, Juffermans NP, Karagiannidis C, Kesecioglu J, Kwizera A, Laffey JG, Mancebo J, Matthay MA, McAuley DF, Mercat A, Meyer NJ, Moss M, Munshi L, Myatra SN, Ng Gong M, Papazian L, Patel BK, Pellegrini M, Perner A, Pesenti A, Piquilloud L, Qiu H, Ranieri MV, Riviello E, Slutsky AS, Stapleton RD, Summers C, Thompson TB, Valente Barbas CS, Villar J, Ware LB, Weiss B, Zampieri FG, Azoulay E, Cecconi M. ESICM guidelines on acute respiratory distress syndrome: definition, phenotyping and respiratory support strategies. Intensive Care Med 2023.

  2. Calfee CS, Delucchi K, Parsons PE, Thompson BT, Ware LB, Matthay MA. Subphenotypes in acute respiratory distress syndrome: latent class analysis of data from two randomised controlled trials. Lancet Respir Med. 2014;2:611–20.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Famous KR, Delucchi K, Ware LB, Kangelaris KN, Liu KD, Thompson BT, Calfee CS. Acute respiratory distress syndrome subphenotypes respond differently to randomized fluid management strategy. Am J Respir Crit Care Med. 2017;195:331–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Calfee CS, Delucchi KL, Sinha P, Matthay MA, Hackett J, Shankar-Hari M, McDowell C, Laffey JG, O’Kane CM, McAuley DF. Acute respiratory distress syndrome subphenotypes and differential response to simvastatin: secondary analysis of a randomised controlled trial. Lancet Respir Med. 2018;6:691–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Sinha P, Delucchi KL, Thompson BT, McAuley DF, Matthay MA, Calfee CS. Latent class analysis of ARDS subphenotypes: a secondary analysis of the statins for acutely injured lungs from sepsis (SAILS) study. Intensive Care Med. 2018;44:1859–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Sinha P, Delucchi KL, Chen Y, Zhuo H, Abbott J, Wang C, Wickersham N, McNeil JB, Jauregui A, Ke S, Vessel K, Gomez A, Hendrickson CM, Kangelaris KN, Sarma A, Leligdowicz A, Liu KD, Matthay MA, Ware LB, Calfee CS. Latent class analysis-derived subphenotypes are generalisable to observational cohorts of acute respiratory distress syndrome: a prospective study. Thorax. 2022;77:13–21.

    Article  PubMed  Google Scholar 

  7. Bos LD, Schouten LR, van Vught LA, Wiewel MA, Ong DSY, Cremer O, Artigas A, Martin-Loeches I, Hoogendijk AJ, van der Poll T, Horn J, Juffermans N, Calfee CS, Schultz MJ. Identification and validation of distinct biological phenotypes in patients with acute respiratory distress syndrome by cluster analysis. Thorax. 2017;72:876–83.

    Article  CAS  PubMed  Google Scholar 

  8. Delucchi K, Famous KR, Ware LB, Parsons PE, Thompson BT, Calfee CS. Stability of ARDS subphenotypes over time in two randomised controlled trials. Thorax. 2018;73:439–45.

    Article  PubMed  Google Scholar 

  9. Sinha P, Furfaro D, Cummings MJ, Abrams D, Delucchi K, Maddali MV, He J, Thompson A, Murn M, Fountain J, Rosen A, Robbins-Juarez SY, Adan MA, Satish T, Madhavan M, Gupta A, Lyashchenko AK, Agerstrand C, Yip NH, Burkart KM, Beitler JR, Baldwin MR, Calfee CS, Brodie D, O’Donnell MR. Latent Class Analysis reveals COVID-19-related acute respiratory distress syndrome subgroups with differential responses to corticosteroids. Am J Respir Crit Care Med. 2021;204:1274–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Shankar-Hari M, Fan E, Ferguson ND. Acute respiratory distress syndrome (ARDS) phenotyping. Intensive Care Med. 2019;45:516–9.

    Article  CAS  PubMed  Google Scholar 

  11. Roca O, Riera J, Torres F, Masclans JR. High-flow oxygen therapy in acute respiratory failure. Respir Care. 2010;55:408–13.

    PubMed  Google Scholar 

  12. Rochwerg B, Granton D, Wang DX, Helviz Y, Einav S, Frat JP, Mekontso-Dessap A, Schreiber A, Azoulay E, Mercat A, Demoule A, Lemiale V, Pesenti A, Riviello ED, Mauri T, Mancebo J, Brochard L, Burns K. High flow nasal cannula compared with conventional oxygen therapy for acute hypoxemic respiratory failure: a systematic review and meta-analysis. Intensive Care Med. 2019;45:563–72.

    Article  CAS  PubMed  Google Scholar 

  13. Ferreyro BL, Angriman F, Munshi L, Del Sorbo L, Ferguson ND, Rochwerg B, Ryu MJ, Saskin R, Wunsch H, da Costa BR, Scales DC. Association of noninvasive oxygenation strategies with all-cause mortality in adults with acute hypoxemic respiratory failure: a systematic review and meta-analysis. JAMA. 2020;324:57–67.

    Article  CAS  PubMed  Google Scholar 

  14. Rochwerg B, Einav S, Chaudhuri D, Mancebo J, Mauri T, Helviz Y, Goligher EC, Jaber S, Ricard JD, Rittayamai N, Roca O, Antonelli M, Maggiore SM, Demoule A, Hodgson CL, Mercat A, Wilcox ME, Granton D, Wang D, Azoulay E, Ouanes-Besbes L, Cinnella G, Rauseo M, Carvalho C, Dessap-Mekontso A, Fraser J, Frat JP, Gomersall C, Grasselli G, Hernandez G, Jog S, Pesenti A, Riviello ED, Slutsky AS, Stapleton RD, Talmor D, Thille AW, Brochard L, Burns KEA. The role for high flow nasal cannula as a respiratory support strategy in adults: a clinical practice guideline. Intensive Care Med. 2020;46:2226–37.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Garcia-de-Acilu M, Marin-Corral J, Vazquez A, Ruano L, Magret M, Ferrer R, Masclans JR, Roca O. Hypoxemic patients with bilateral infiltrates treated with high-flow nasal cannula present a similar pattern of biomarkers of inflammation and injury to acute respiratory distress syndrome patients. Crit Care Med. 2017;45:1845–53.

    Article  PubMed  Google Scholar 

  16. Kangelaris KN, Ware LB, Wang CY, Janz DR, Zhuo H, Matthay MA, Calfee CS. Timing of intubation and clinical outcomes in adults with acute respiratory distress syndrome. Crit Care Med. 2016;44:120–9.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Ranieri VM, Tonetti T, Navalesi P, Nava S, Antonelli M, Pesenti A, Grasselli G, Grieco DL, Menga LS, Pisani L, Boscolo A, Sella N, Pasin L, Mega C, Pizzilli G, Dell’Olio A, Dongilli R, Rucci P, Slutsky AS. High-flow nasal oxygen for severe hypoxemia: oxygenation response and outcome in patients with COVID-19. Am J Respir Crit Care Med. 2022;205:431–9.

    Article  CAS  PubMed  Google Scholar 

  18. Matthay MA, Arabi Y, Arroliga AC, Bernard G, Bersten AD, Brochard LJ, Calfee CS, Combes A, Daniel BM, Ferguson ND, Gong MN, Gotts JE, Herridge MS, Laffey JG, Liu KD, Machado FR, Martin TR, McAuley DF, Mercat A, Moss M, Mularski RA, Pesenti A, Qiu H, Ramakrishnan N, Ranieri M, Riviello ED, Rubin E, Slutsky A, Thompson BT, Twagirumugabe T, Ware LB, Wick KD. A new global definition of acute respiratory distress syndrome. Am J Respir Crit Care Med 2023.

  19. van der Ven F, Valk CMA, Blok S, Brouwer MG, Go DM, Lokhorst A, Swart P, van Meenen DMP, Paulus F, Schultz MJ. Broadening the Berlin definition of ARDS to patients receiving high-flow nasal oxygen: an observational study in patients with acute hypoxemic respiratory failure due to COVID-19. Ann Intensive Care. 2023;13:64.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Charrad M, Ghazzali N, Boiteau V, Niknafs A. NbClust: An R package for determining the relevant number of clusters in a data set. J Stat Softw. 2014;61:1–36.

    Article  Google Scholar 

  21. Youden WJ. Index for rating diagnostic tests. Cancer. 1950;3:32–5.

    Article  CAS  PubMed  Google Scholar 

  22. Matthay MA, Zemans RL, Zimmerman GA, Arabi YM, Beitler JR, Mercat A, Herridge M, Randolph AG, Calfee CS. Acute respiratory distress syndrome. Nat Rev Dis Primers. 2019;5:18.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Ware LB, Koyama T, Zhao Z, Janz DR, Wickersham N, Bernard GR, May AK, Calfee CS, Matthay MA. Biomarkers of lung epithelial injury and inflammation distinguish severe sepsis patients with acute respiratory distress syndrome. Crit Care. 2013;17:R253.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Jabaudon M, Futier E, Roszyk L, Chalus E, Guerin R, Petit A, Mrozek S, Perbet S, Cayot-Constantin S, Chartier C, Sapin V, Bazin JE, Constantin JM. Soluble form of the receptor for advanced glycation end products is a marker of acute lung injury but not of severe sepsis in critically ill patients. Crit Care Med. 2011;39:480–8.

    Article  CAS  PubMed  Google Scholar 

  25. Agrawal A, Matthay MA, Kangelaris KN, Stein J, Chu JC, Imp BM, Cortez A, Abbott J, Liu KD, Calfee CS. Plasma angiopoietin-2 predicts the onset of acute lung injury in critically ill patients. Am J Respir Crit Care Med. 2013;187:736–42.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Zou L, Dang W, Tao Y, Zhao H, Yang B, Xu X, Li Y. The IL-33/ST2 axis promotes acute respiratory distress syndrome by natural killer T cells. Shock. 2023;59:902–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Martínez-González I, Roca O, Masclans JR, Moreno R, Salcedo MT, Baekelandt V, Cruz MJ, Rello J, Aran JM. Human mesenchymal stem cells overexpressing the IL-33 antagonist soluble IL-1 receptor-like-1 attenuate endotoxin-induced acute lung injury. Am J Respir Cell Mol Biol. 2013;49:552–62.

    Article  PubMed  Google Scholar 

  28. Parsons PE, Eisner MD, Thompson BT, Matthay MA, Ancukiewicz M, Bernard GR, Wheeler AP. Lower tidal volume ventilation and plasma cytokine markers of inflammation in patients with acute lung injury. Crit Care Med. 2005;33:1–6 (discussion 230-232).

    Article  CAS  PubMed  Google Scholar 

  29. Dalmaijer ES, Nord CL, Astle DE. Statistical power for cluster analysis. BMC Bioinformatics 2022; 23.

Download references


Not applicable.


The princeps study was supported, in part, by a grant from Instituto de Salud Carlos III-Fondo Europeo de Desarrollo Regional (FEDER) (PI14/01420) and a grant from Federación Española del Enfermo Crítico (FEEC) 2015. No specific support was done for this re-analysis.

Author information

Authors and Affiliations



P-LB; B-GC; OR; and MG wrote the first draft. P-LB and JC did the statistical analysis. MS, AP, MM, J-RM, OR, and MG did the princeps study, took care of the patients, and collected the data. All the authors contributed to the interpretation of the results and reviewed and approved the final manuscript.

Corresponding author

Correspondence to Oriol Roca.

Ethics declarations

Ethics approval and consent to participate

This is a post hoc analysis of a multicenter prospective observational study previously published. The primary study was approved by the Ethics Committee of each participating hospital, and written consent from patients or their relatives was obtained.

Consent for publication

Not applicable.

Competing interests

O.R reports receiving a research grant from Hamilton Medical AG and Fisher&Paykel Healthcare Ltd, speaker fees from Hamilton Medical AG, Fisher&Paykel Healthcare Ltd, Aerogen Ltd and Ambu, and non-financial research support from Timpel; all outside the submitted work. B-G.C. received speaker fees from Baxter and was a member of an advisory board for Roche Diagnostics. J-R.M received Grants, travels, and non-financial support from Fisher & Paykel.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Detailed statistical analysis.

Additional file 2.

Variables included in cluster analysis.

Additional file 3.

Cluster plot.

Additional file 4.

Plasma concentration of the different biomarkers.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Blot, PL., Chousterman, B.G., Santafè, M. et al. Subphenotypes in patients with acute respiratory distress syndrome treated with high-flow oxygen. Crit Care 27, 419 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: