key: cord-0735881-0q0k5s78
authors: Marti, Sergi; Carsin, Anne-Elie; Sampol, Júlia; Pallero, Mercedes; Aldas, Irene; Marin, Toni; Lujan, Manel; Lalmolda, Cristina; Sabater, Gladis; Bonnin-Vilaplana, Marc; Peñacoba, Patricia; Martinez-Llorens, Juana; Tárrega, Julia; Bernadich, Óscar; Córdoba-Izquierdo, Ana; Lozano, Lourdes; Mendez, Susana; Vélez-Segovia, Eduardo; Prina, Elena; Eizaguirre, Saioa; Balañá-Corberó, Ana; Ferrer, Jaume; Garcia-Aymerich, Judith
title: Higher mortality and intubation rate in COVID-19 patients treated with noninvasive ventilation compared with high-flow oxygen or CPAP
date: 2022-04-20
journal: Sci Rep
DOI: 10.1038/s41598-022-10475-7
sha: ee94305699f3d15908392a219b85adfa8c986066
doc_id: 735881
cord_uid: 0q0k5s78

The effectiveness of noninvasive respiratory support in severe COVID-19 patients is still controversial. We aimed to compare the outcome of patients with COVID-19 pneumonia and hypoxemic respiratory failure treated with high-flow oxygen administered via nasal cannula (HFNC), continuous positive airway pressure (CPAP) or noninvasive ventilation (NIV), initiated outside the intensive care unit (ICU) in 10 university hospitals in Catalonia, Spain. We recruited 367 consecutive patients aged ≥ 18 years who were treated with HFNC (155, 42.2%), CPAP (133, 36.2%) or NIV (79, 21.5%). The main outcome was intubation or death at 28 days after respiratory support initiation. After adjusting for relevant covariates and taking patients treated with HFNC as reference, treatment with NIV showed a higher risk of intubation or death (hazard ratio 2.01; 95% confidence interval 1.32–3.08), while treatment with CPAP did not show differences (0.97; 0.63–1.50). In the context of the pandemic and outside the intensive care unit setting, noninvasive ventilation for the treatment of moderate to severe hypoxemic acute respiratory failure secondary to COVID-19 resulted in higher mortality or intubation rate at 28 days than high-flow oxygen or CPAP. This finding may help physicians to choose the best noninvasive respiratory support treatment in these patients. Clinicaltrials.gov identifier: NCT04668196.

Study design. A multicentre, retrospective cohort study of COVID-19 patients followed from NIRS initiation up to 28 days or death, whichever occurred first.

We included a consecutive sample of patients aged at least 18 years who had initiated NIRS treatment for HARF related to COVID-19 pneumonia outside the ICU at any of the 10 participating university hospitals, during the first pandemic surge, between 1 March and 30 April 2020. All participating hospitals belong to the National Health System of Catalonia, Spain, and attend a population of around 4.3 million inhabitants. COVID-19 diagnosis was confirmed through reverse-transcriptase-polymerase-chain-reaction assays performed on nasopharyngeal swab specimens. From a total of 419 candidate patients, we excluded those with: (1) respiratory failure not related to COVID-19 (e.g., cardiogenic pulmonary edema as primary cause of respiratory failure); (2) rejection or early intolerance to any NIRS treatment; (3) pregnancy; (4) nosocomial infection; and (5) PaCO 2 above 45 mm Hg. A total of 367 patients were finally included in the study (Fig. 1) , which was approved by the research ethics committee at each participating hospital (study coordinator centre, Hospital Vall d'Hebron, Barcelona; protocol No. PR(AG)265/2020). Research was performed in accordance with the Declaration of Helsinki. The requirement of informed consent was waived due to the retrospective nature of the study. Treatment strategies. The NIRS treatments evaluated were high-flow oxygen administered via nasal cannula (HFNC), continuous positive airway pressure (CPAP), and noninvasive ventilation (NIV). According to current Spanish recommendations 8 , criteria for initiating respiratory support were moderate to severe dyspnoea, respiratory rate > 30 bpm, or PaO2/FiO2 < 200 mmHg, screened either at hospital admission or ward admission. The decision regarding the choice of treatment was taken by the pulmonologist in charge of the patient's care, with HFNC usually as the first step after the failure of conventional oxygen therapy 8 , and taking into account the availability of NIRS devices at each centre.

In the HFNC group, heated and humidified oxygen was applied through nasal prongs, at an initial flow rate of 50-60 lpm if tolerated. CPAP was initially set at 8-10 cm H 2 O and then adjusted according to tolerance and clinical response. In the NIV group, a pressure support ventilator mode was adjusted; a high positive end-expiratory pressure (PEEP) and a low support pressure were used to set a tidal volume < 9 ml/kg of predicted body weight 8 .

Patients' characteristics. Study data were collected and managed using REDCap electronic data capture tools hosted at ISGlobal (Institut de Salut Global, Barcelona) 23 . We obtained patients' data from electronic medical records using a modified version of the standardized International Severe Acute Respiratory and Emerging Infection Consortium (ISARIC) COVID-19 case report forms 24 , including: (i) demographics (age, sex, ethnicity); (ii) smoking status; (iii) chronic conditions (cardiac disease, respiratory disease, kidney disease, neoplasm, dementia, obesity, neurological conditions, liver disease, diabetes, and a modified Charlson comorbidity index) 25 ; (iv) symptoms at admission and physical signs at NIRS initiation (days since the onset of COVID-19 www.nature.com/scientificreports/ symptoms, temperature, heart rate, systolic and diastolic blood pressure, respiratory rate, and Quick Sequential Organ Failure Assessment (qSOFA) score) 26 ; (v) arterial blood gases at NIRS initiation (PaO 2 /F I O 2 ratio calculated for patients with available PaO 2 , and imputed from SpO 2 for the 33% of patients without PaO 2 ) 27 ; (vi) laboratory blood parameters at NIRS initiation; (vii) chest X-ray findings (unilateral or bilateral pneumonia); and (viii) treatment received during admission (highest level of care received outside ICU, ICU admission, NIRS as ceiling of treatment, awake prone positioning, and drug treatments).

Study outcomes. The primary outcome was treatment failure, defined as endotracheal intubation or death within 28 days of NIRS initiation. Secondary outcomes were 28-day mortality, endotracheal intubation at day 28, in-hospital mortality, and duration of hospital stay.

With an expected frequency of 50% for intubation or death in patients with HARF and treated by NIRS 28 , 300 patients were needed in order to detect a significant difference greater than 20% between the types of NIRS evaluated in the present study, with an alpha risk of 0.05 and a statistical power of 80%. Characteristics of the patients at baseline according to NIRS treatment were described by mean and standard deviation, median and 25th and 75th percentiles (P25 and P75) and by absolute and relative frequencies, and compared using Chi2, Anova and Kruskal Wallis tests. Given the small number of missing information and that missing were considered at random, we conducted a complete case approach.

Kaplan-Meier curves described the crude event-free rate in each NIRS group and were compared by means of the log-rank test. Multivariable Cox proportional-hazards regression models were used to estimate the hazard ratios (HR) for patients treated with NIV and CPAP as compared to HFNC (the reference group), adjusting for age, sex, and variables found to be significantly different between treatments at baseline (hospital, date of admission and sleep apnea). D-dimer levels and respiratory rate at baseline were also significantly associated with treatment, but since they had missing values for 82 and 41 patients respectively, these variables were only included in a sensitivity analysis.

To account for the potential effect modification, analyses were stratified according to hypoxemia severity (moderate-severe: PaO 2 /F I O 2 < 150 mm Hg; mild-moderate: PaO 2 /F I O 2 ≥ 150 mm Hg) 4 . To assess the potential impact of NIRS treatment settings, we compared outcomes within NIRS-group according to: flow in the HFNC group (> 50 vs. ≤ 50 L/min), pressure in the CPAP group (> 10 vs. ≤ 10 cm H 2 O), and PEEP in the NIV group (> 10 vs. ≤ 10 cm H 2 O).

Sensitivity analyses included: (1) repeating models excluding patients who changed their initial NIRS treatment during the course of the hospitalization to another NIRS treatment (crossover, n = 44); (2) excluding patients with missing measured PaO 2 /F I O 2 (n = 123); (3) excluding patients receiving NIRS as ceiling of treatment (n = 140); and (4) additionally adjusting models for, one at a time, D-dimer levels, respiratory rate, systemic corticosteroid use and Charlson index.

All analyses were performed using StataCorp. 2019. Stata Statistical Software: Release 16. College Station, TX: StataCorp LLC.

All authors have approved the submission and provide consent to publish.

Patients' characteristics. Among the 367 patients included in the study, 155 were treated with HFNC (42.2%), 133 with CPAP (36.2%), and 79 with NIV (21.5%).

Most patients were male (72%), and the mean age was 67.5 years (SD 11.2). Chronic conditions were frequent (35% of the sample had a Charlson comorbidity index ≥ 2) and did not differ between NIRS treatment groups, except for sleep apnea (more common in the NIV-treated group, Table 1 and Table S1 ). At the initiation of NIRS, patients had moderate to severe hypoxemia (median PaO 2 /F I O 2 125.5 mm Hg, P25-P75: 81-174). Clinical severity and laboratory values were well balanced between the groups ( Table 2 and Table S2 ), except for respiratory rate (higher in patients treated with NIV).

Treatments. The NIRS treatments applied were not equally distributed among participating hospitals, although HFNC or CPAP were the first NIRS treatment choice at all centers (Table S1 ). Differences were also found in the NIRS treatments applied according to the date of admission: HFNC was the most frequent treatment early in the period (before 23 March), while CPAP was the most frequent choice in the second and the third periods (Table 1 , p = 0.008). Noninvasive respiratory support treatments were applied as ceiling of treatment in 140 patients (38%) ( Table 3) . Table S3 shows the NIRS settings.

Primary and secondary outcomes. The cumulative percentage of patients who had received intubation or who had died by day 28 (primary outcome) was 45.8% in the HFNC group, 36.8% in the CPAP group, and 60.8% in the NIV group (Fig. 2a) . After adjustment, and taking patients treated with HFNC as reference, patients who underwent NIV had a higher risk of intubation or death at 28 days (HR 2.01, 95% CI 1.32-3.08), while those treated with CPAP did not present differences (HR 0.97, 95% CI 0.63-1.50) ( Table 4) .

As for secondary outcomes, patients treated with NIV had a significantly higher risk of endotracheal intubation, 28-day mortality, and in-hospital mortality than patients treated with HFNC, while no differences were observed between CPAP and HFNC (Fig. 2b, www.nature.com/scientificreports/ patients were treated with NIRS as ceiling of treatment and 26 (11.5%) were treated with NIRS not regarded as ceiling of treatment. Days between NIRS initiation and intubation (median (P25-P75) 3 (1-5), 3.5 (2-7), and 3 (3) (4) (5) , for HFNC, CPAP, and NIV groups respectively; p = 0.341) and the length of hospital stay did not differ between groups (Table 4) . Outcomes by hospital are listed in Table S4 .

Stratified and sensitivity analyses. In patients with mild-moderate hypoxaemia, CPAP, but not NIV, treatment was associated with reduced outcome risk compared to HFNC (Table S5 ). The analyses excluding patients with missing PaO 2 /F I O 2 or receiving NIRS as ceiling of treatment showed similar associations to those observed in the main analysis (Tables S6 and S7 , respectively). No differences were found when we performed within NIRS-group comparisons according to settings applied (Table S8) . www.nature.com/scientificreports/ During the follow-up period, 44 patients (12%) switched to another NIRS treatment: eight (5%) in the HFNC group (treated subsequently with NIV), 28 (21%) in the CPAP group (13 switched to HFNC, and 15 to NIV), and eight (10%) in the NIV group (seven treated with HFNC, and one with CPAP). Excluding these patients showed no relevant changes in the associations observed (Table S9) . Additional adjustment for D-dimer, respiratory rate, Charlson index, or treatment with systemic corticosteroids produced very similar results (Table S10 ). 

Charlson index † , med (P25-P75) 1 (0-2) 1 (0-2) 1 (0-2) 1 (0-2) 0.

This study shows that noninvasive ventilation initiated outside the ICU for the treatment of hypoxemic acute respiratory failure secondary to COVID-19 resulted in higher mortality or intubation rate at 28 days (i.e., treatment failure) than high-flow oxygen or CPAP. These results were robust to a number of stratified and sensitivity analyses.

Most previous data on the effectiveness of NIRS treatments in severe COVID-19 patients came from studies which had limited sample sizes and were not designed to compare the different techniques [13] [14] [15] 17, 18 . In the only available study (also observational) comparing NIV, HFNC and CPAP outside the ICU 16 , conducted in Italy, the authors did not find differences between treatments in mortality or intubation at 30 days. The discrepancy between these results and ours may be due to differences in the characteristics of the patients included. First, in the Italian study, the mean PaO 2 /F I O 2 ratio was 152 mm Hg, suggesting a less severe respiratory failure than in our patients (125 mm Hg). In the stratified analysis of our cohort, planned a priori, patients with a PaO 2 /F I O 2 ratio above 150 responded similarly to HFNC and NIV treatments, suggesting that the severity of the hypoxemia might predict the success of NIV, as previously reported in non-COVID patients 4, 28, 29 . Second, the Italian study did not provide data on PaCO 2 , meaning that the improvements with NIV might have been attributable to the inclusion of some patients with hypercapnic respiratory failure, who were excluded in our study.

Recently, the effectiveness of CPAP or HFNC compared with conventional oxygen therapy was assessed in the RECOVERY-RS multicentric randomized clinical trial, in 1,273 COVID-19 patients with HARF who were deemed suitable for tracheal intubation if treatment escalation was required 20 . In this study, the requirement of intubation or mortality within 30 days (primary outcome) was significantly lower with CPAP (36%) than with conventional oxygen therapy (45%; absolute difference, − 8% [95% CI, − 15% to − 1%], p = 0.03). This improvement was mostly driven by a reduction in the need of intubation, but no differences in mortality were seen (16.7% vs 19.2%, respectively). No significant differences in the main outcome were found between HFNC (44%) vs conventional oxygen therapy (45%; absolute difference, − 1% [95% CI, − 8% to 6%], p = 0.83). The main difference in respect to our study was the better outcomes of CPAP compared with HFNC. However, the RECOVERY-RS study may have been underpowered for the comparison of HFNC vs conventional oxygen therapy due to early study termination and the number of crossovers among groups (11.5% of HFNC and 23.6% of conventional oxygen treated patients).

There are several possible explanations for the poor outcome of COVID-19 patients undergoing NIV in our study. First, NIV has been reported to produce overdistension, compounded by the respiratory effort itself 30 , which could result in ventilation-induced lung injury due to the excessive increases in tidal volumes 28, 31 . This risk would be avoided in CPAP and HFNC because they improve oxygenation without changing tidal volume 32, 33 . Unfortunately, tidal volume measurements during NIV were not available in our study to support or reject this hypothesis. Second, patient-ventilator asynchronies might have arisen in NIV-treated patients making more difficult their management outside the ICU setting and thereby explaining, at least partially, their worse outcomes. Third, a bench study has recently reported that some approaches to minimize aerosol dispersion can modify ventilator performance 34 . In short, the addition of intentional leaks, as in our study, led to a lower maximal pressure without a significant impact on the work of breathing and without increasing patient-ventilator asynchronies 34 . Table 3 . Inpatient characteristics and treatments according to non-invasive respiratory support group. P-value from Chi 2 test. *Anticoagulation had 1 missing value. † Systemic corticosteroids included prednisone (n = 6), methylprednisolone (n = 223), dexamethasone (n = 21) and hydrocortisone (n = 1). www.nature.com/scientificreports/ Then, in the present work, we believe that the availability of trained pulmonologists to adjust ventilator settings may have overcome this aspect. Fourth, non-responders to NIV could have suffered a delay in intubation, but in our study the time to intubation was similar in the three NIRS groups, thus making this explanation less likely. Fifth, we cannot exclude the possibility that NIV implied a more complicated clinical course than HFNC or CPAP. Patients undergoing NIV may require some degree of sedation to tolerate the technique, but unfortunately we have no data on this regard. Furthermore, NIV and CPAP may impair expectoration which could contribute to www.nature.com/scientificreports/ bacterial infections, although this hypothesis remains unknown with the present data. Finally, we cannot rule out the possibility that NIV was tolerated worse than HFNC or CPAP, which would have reduced adherence and lowered the effectiveness of the therapy. However, the number of patients abandoning their original treatment was nearly twice as high in the CPAP group than in the NIV group.

In the treatment of HARF with CPAP or NIV the interface via which these treatments are applied should be considered, since better outcomes have been reported with a helmet interface than with face masks in non-COVID patients 6, 35 , possibly due to a greater tolerance of the helmet and a more effective delivery of PEEP 36 . As noted above, a single randomized study has evaluated helmet NIV against HFNC in COVID- 19 19 , and, in spite of the lower intubation rate in the helmet NIV group, no differences in 28-day mortality were registered. More studies are needed to define the place of treatment with helmet CPAP or NIV in respiratory failure due to COVID-19, together with other NIRS strategies. In our study, CPAP and NIV treatments were applied via oronasal and full face masks, reflecting the fact that most hospitals in our country have little experience with the helmet interface.

Our study supports several guidelines 37,38 that favor HFNC and CPAP over NIV for the treatment of HARF in COVID-19 patients, but to our knowledge no previous data have been published in support of this recommendation. Furthermore, our results suggest that the severity of the hypoxemic respiratory failure might help physicians to decide which specific NIRS technique could be better for a patient. However, the retrospective design of our study does not allow establishing a causative link between NIV and the worse clinical outcomes observed. Obviously, reaching a definitive conclusion on this point will require further studies with better phenotypic characterization of patients, and considering additional factors implicated in the response to therapies such as the interface used or the monitoring of the inspiratory effort.

This study has some limitations. First, the observational design could have resulted in residual confounding by selection bias. However, the inclusion of patients was consecutive and the collection of variables was really comprehensive. Moreover, NIRS treatment groups exhibited only minor differences which were accounted for in the multivariable and sensitivity analyses thus minimizing the selection bias risk. Although treatment received and outcomes differed by hospital, this fact was taken into account through adjustment. Second, we must be cautious before extrapolating our results to other nonemergency situations. Our study was carried out during the first wave of the pandemics when the healthcare system was overwhelmed and many patients were treated outside ICU facilities. As mentioned above, NIV might have better outcomes in a more controlled setting allowing an optimal critical care 39 . However, the scarcity of critical care resources has remained along the different pandemic surges until now and this scenario is unfortunately frequent in other health care systems around the world. Thus, we believe that our results may be useful for a great number of physicians treating COVID-19 patients around the world. Third, crossovers could have been responsible for differences observed between NIRS treatments but their proportion was small (12%) and our results did not change when these patients were excluded. Fourth, it could be argued that changes in treatment strategies over the timeframe of the study may have led to differential effects of the NIRS. Nevertheless, we do not think it may have influenced our results, because analyses were adjusted for relevant treatments such as systemic corticosteroids 40 and included the time period as a covariate. And finally, due to the shortage of critical care ventilators at the height of the pandemic, some patients were treated with home devices with limited FiO2 delivery capability and, therefore, could have been undertreated 41, 42 . www.nature.com/scientificreports/ However, as more home devices were used in the CPAP group (81.6% vs. 38% in the NIV group; Table S3 ), and better outcomes were recorded in the CPAP-treated patients, our result do not support this concern.

The main strength of this study is, in our opinion, its real-life design that allows obtaining the effectiveness of these techniques in the clinical setting. In the current situation with few available data from randomized control trials regarding the best choice to treat COVID-19 patients with noninvasive respiratory support, data from reallife studies like ours may be appropriate 43 . These data are complementary and still useful later on by including some patients usually excluded from randomized studies; patients with do-not-intubate orders are an example and, obviously, they represent a challenge for the physician responsible to decide the best therapeutic strategy.

In conclusion, the present real-life study shows that, in the context of the pandemic and outside the intensive care unit setting, noninvasive ventilation for the treatment of hypoxemic acute respiratory failure secondary to COVID-19 resulted in higher treatment failure than high-flow oxygen or CPAP. These findings may be relevant for many physicians elsewhere since the successive pandemic surges result in overwhelmed health care systems, leading to the need for severe COVID-19 patients to be treated out of critical care settings.

All data generated or analyzed during this study are included in this published article and its supplementary information files. 

Critical care utilization for the COVID-19 outbreak in Lombardy, Italy: Early experience and forecast during an emergency response

Features of 20 133 UK patients in hospital with covid-19 using the ISARIC WHO Clinical Characterization Protocol: Prospective observational cohort study

Treatment of acute hypoxemic nonhypercapnic respiratory insufficiency with continuous positive airway pressure delivered by a face mask: A randomized controlled trial

Noninvasive ventilation of patients with acute respiratory distress syndrome. Insights from the LUNG SAFE study

Official ERS/ATS clinical practice guidelines: Noninvasive ventilation for acute respiratory failure

Association of noninvasive oxygenation strategies with all-cause mortality in adults with acute hypoxemic respiratory failure: A systematic review and meta-analysis

Preliminary findings on control of dispersion of aerosols and droplets during high-velocity nasal insufflation therapy using a simple surgical mask: Implications for the high-flow nasal cannula

Clinical consensus recommendations regarding non-invasive respiratory support in the adult patient with acute respiratory failure secondary to SARS-CoV-2 infection

Specialty Guides for Patient Management During the Coronavirus Pandemic. Guidance for the Role and Use of Non-invasive Respiratory Support in Adult Patients with COVID-19

Surviving sepsis campaign: Guidelines on the management of critically ill adults with coronavirus disease 2019 (COVID-19)

Early consensus management for non-ICU acute respiratory failure SARS-CoV-2 emergency in Italy: From ward to trenches

High-flow, noninvasive ventilation and awake (nonintubation) proning in patients with coronavirus disease 2019 with respiratory failure

Continuous positive airway pressure to avoid intubation in SARS-CoV-2 pneumonia: A two-period retrospective case-control study

High-flow nasal cannula oxygen therapy to treat patients with hypoxemic acute respiratory failure consequent to SARS-CoV-2 infection

Helmet CPAP treatment in patients with COVID-19 pneumonia: A multicentre cohort study

Feasibility and clinical impact of out-of-ICU noninvasive respiratory support in patients with COVID-19-related pneumonia

High-flow nasal cannula in critically III patients with severe COVID-19

Continuous positive airway pressure in COVID-19 patients with moderate-to-severe respiratory failure

Effect of helmet noninvasive ventilation vs. high-flow nasal oxygen on days free of respiratory support in patients with COVID-19 and moderate to severe hypoxemic respiratory failure: The HENIVOT randomized clinical trial

Effect of noninvasive respiratory strategies on intubation or mortality among patients with acute hypoxemic respiratory failure and COVID-19 The RECOVERY-RS randomized clinical trial

Oxygen therapy for acutely ill medical patients: A clinical practice guideline

Ethical recommendations for a difficult decision-making in intensive care units due to the exceptional situation of crisis by the COVID-19 pandemia: A rapid review & consensus of experts

The REDCap consortium: Building an international community of software platform partners

The International Severe Acute Respiratory and Emerging Infection Consortium (ISARIC

A new method of classifying prognostic comorbidity in longitudinal studies: Development and validation

The third international consensus definitions for sepsis and septic shock (Sepsis-3)

Nonlinear imputation of PaO2/FiO2 from SpO2/FiO2 among patients with acute respiratory distress syndrome

High-flow oxygen through nasal cannula in acute hypoxemic respiratory failure

Non-invasive ventilation for acute hypoxemic respiratory failure: Intubation rate and risk factors

Mechanical ventilation to minimize progression of lung injury in acute respiratory failure

Failure of noninvasive ventilation for de novo acute hypoxemic respiratory failure: Role of tidal volume

Physiologic effects of noninvasive ventilation during acute lung injury

Physiologic effects of high-flow nasal cannula in acute hypoxemic respiratory failure

Recommended approaches to minimize aerosol dispersion of SARS-CoV-2 during noninvasive ventilatory support can cause ventilator performance deterioration: A benchmark comparative study

Effect of noninvasive ventilation delivered helmet vs. face mask on the rate of endotracheal intubation in patients with acute respiratory distress syndrome: A randomized clinical trial

Patient self-inflicted lung injury and positive end-expiratory pressure for safe spontaneous breathing

Management of hospitalised adults with coronavirus disease 2019 (COVID-19): A European Respiratory Society living guideline

Expert consensus statements for the management of COVID-19-related acute respiratory failure using Delphi method

Effect of prone position on respiratory parameters, intubation and death rate in COVID-19 patients: Systematic review and meta-analysis

Dexamethasone in hospitalized patients with Covid-19

Inspired oxygen fraction achieved with a portable ventilator: Determinant factors

Oxygen supplementation in noninvasive home mechanical ventilation: The crucial roles of CO2 exhalation systems and leakages

Noninvasive strategies in COVID-19: Epistemology, randomised trials, guidelines, physiology

The authors wish to thank Barcelona Research Network (BRN) for their logistical and administrative support and to Rosa Llòria for her assistance and technical help in the edition of the paper. 

ISGlobal acknowledges support from the Spanish Ministry of Science and Innovation through the "Centro de Excelencia Severo Ochoa 2019-2023" Program (CEX2018-000806-S), and from the Generalitat de Catalunya through the CERCA Program.

The authors declare no competing interests.

The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-022-10475-7.Correspondence and requests for materials should be addressed to S.M.Reprints and permissions information is available at www.nature.com/reprints.Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.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 http:// creat iveco mmons. org/ licen ses/ by/4. 0/.