key: cord-0979279-gdnc7b8k
authors: Gamberini, Lorenzo; Mazzoli, Carlo Alberto; Prediletto, Irene; Sintonen, Harri; Scaramuzzo, Gaetano; Allegri, Davide; Colombo, Davide; Tonetti, Tommaso; Zani, Gianluca; Capozzi, Chiara; Dalpiaz, Giorgia; Agnoletti, Vanni; Cappellini, Iacopo; Melegari, Gabriele; Damiani, Federica; Fusari, Maurizio; Gordini, Giovanni; Laici, Cristiana; Lanza, Maria Concetta; Leo, Mirco; Marudi, Andrea; Papa, Raffaella; Potalivo, Antonella; Montomoli, Jonathan; Taddei, Stefania; Mazzolini, Massimiliano; Ferravante, Anna Filomena; Nicali, Roberta; Ranieri, Vito Marco; Russo, Emanuele; Volta, Carlo Alberto; Spadaro, Savino; Castelli, Andrea; Marchio, Annalisa; Bacchilega, Igor; Bernabé, Laura; Facondini, Francesca; Morini, Luca; Bissoni, Luca; Viola, Lorenzo; Meconi, Tommaso; Pavoni, Vittorio; Venni, Angelica; Pagni, Aline; Cleta, Patrizia Pompa; Cavagnino, Marco; Guzzo, Alessia; Malfatto, Anna; Adduci, Angelina; Pareschi, Silvia; Bertellini, Elisabetta; Maccieri, Jessica; Marinangeli, Elisa; Racca, Fabrizio; Verri, Marco; Falò, Giulia; Marangoni, Elisabetta; Ottaviani, Irene; Boni, Francesco; Felloni, Giulia; Baccarini, Federico Domenico; Terzitta, Marina; Maitan, Stefano; Tutino, Lorenzo; Senzi, Angelo; Consales, Guglielmo; Becherucci, Filippo; Imbriani, Michele; Orlandi, Paolo; Candini, Silvia; Golfieri, Rita; Ciccarese, Federica; Poerio, Antonio; Muratore, Francesco; Ferrari, Fabio; Mughetti, Martina; Giampalma, Emanuela; Franchini, Loredana; Neziri, Ersenad; Miceli, Marco; Minguzzi, Maria Teresa; Mellini, Lorenzo; Piciucchi, Sara; Monari, Matteo; Valli, Michele; Daniele, Federico; Ferioli, Martina; Nava, Stefano; Lazzari Agli, Luigi Arcangelo; Valentini, Ilaria; Bernardi, Eva; Balbi, Bruno; Contoli, Marco; Padovani, Marianna; Oldani, Stefano; Ravaglia, Claudia; Goti, Patrizio
title: Health-related quality of life profiles, trajectories, persistent symptoms and pulmonary function one year after ICU discharge in invasively ventilated COVID-19 patients, a prospective follow-up study
date: 2021-10-22
journal: Respir Med
DOI: 10.1016/j.rmed.2021.106665
sha: b919c46e2ac62f5dda6ae1938da04a184f428ea3
doc_id: 979279
cord_uid: gdnc7b8k

BACKGROUND: Health-related quality of life (HRQoL) impairment is often reported among COVID-19 ICU survivors, and little is known about their long-term outcomes. We evaluated the HRQoL trajectories between 3 months and 1 year after ICU discharge, the factors influencing these trajectories and the presence of clusters of HRQoL profiles in a population of COVID-19 patients who underwent invasive mechanical ventilation (IMV). Moreover, pathophysiological correlations of residual dyspnea were tested. METHODS: We followed up 178 survivors from 16 Italian ICUs up to one year after ICU discharge. HRQoL was investigated through the 15D instrument. Available pulmonary function tests (PFTs) and chest CT scans at 1 year were also collected. A linear mixed-effects model was adopted to identify factors associated with different HRQoL trajectories and a two-step cluster analysis was performed to identify HRQoL clusters. RESULTS: We found that HRQoL increased during the study period, especially for the significant increase of the physical dimensions, while the mental dimensions and dyspnea remained substantially unchanged. Four main 15D profiles were identified: full recovery (47.2%), bad recovery (5.1%) and two partial recovery clusters with mostly physical (9.6%) or mental (38.2%) dimensions affected. Gender, duration of IMV and number of comorbidities significantly influenced HRQoL trajectories. Persistent dyspnea was reported in 58.4% of patients, and weakly, but significantly, correlated with both DLCO and length of IMV. CONCLUSIONS: HRQoL impairment is frequent 1 year after ICU discharge, and the lowest recovery is found in the mental dimensions. Persistent dyspnea is often reported and weakly correlated with PFTs alterations. TRIAL REGISTRATION: NCT04411459.

As of the end of August 2021, the COVID-19 pandemic resulted in more than 4 million deaths worldwide [1] , and it considerably increased the need for healthcare resources dedicated to the acute phase of the syndrome, posing a hard challenge to the national health systems [2] .

A series of somatic and psychological consequences are described among survivors, including pulmonary function impairment, such as a reduction in forced vital capacity (FVC) and diffusion capacity of lungs for carbon monoxide (DLCO) [3] , neuropsychiatric disorders [4] , and cardiac sequelae [5] . The most frequently complained persistent symptoms are dyspnea [6, 7] and fatigue, as well as anxiety, depression, and sleep problems, which could significantly undermine the health-related quality of life (HRQoL) of these patients, and together identify the so-called long-COVID syndrome [4, 8] .

Moreover, patients who underwent ICU admission for COVID-19 related ARDS (C-ARDS) may also experience post-intensive care syndrome (PICS), which eventually overlaps other symptoms related to COVID-19 itself [9] .

Little is known about the long-term HRQoL of patients who underwent invasive mechanical ventilation (IMV) for COVID-19 pneumonia, and no studies have yet investigated if the persistence of symptoms is associated with long-term pulmonary function impairment in this particular population.

This study reports the results of the second follow-up (1 year) of a population of C-ARDS survivors who required ICU admission, intubation, and IMV [10] , and who were previously evaluated at 3 months after ICU discharge [11] . The primary objective of this study was to describe the HRQoL trajectories between 3 months and 1 year after ICU discharge. Moreover, we defined clusters of patients based on their HRQoL profiles at 1

year and investigated the factors influencing the recovery trajectories of the survivors.

Finally, a subgroup analysis was performed for patients having pulmonary function tests available between 9 and 12 months after ICU discharge to evaluate the relationship between the persistence and entity of dyspnea and pulmonary function.

This prospective multicenter observational study involved 16 Italian ICUs. Patients admitted to the participating ICUs from 22 nd February through 4 th May 2020 (the end of lockdown in Italy), who survived hospital discharge were subsequently followed up until one year after ICU discharge.

The study was approved by the Institutional Review Board (IRB) of the study coordinator centre (Maggiore Hospital, Bologna, Italy) and by each institutional review committee of the participating hospitals. Informed consent was asked at the time of the first follow-up if not obtained during ICU stay according to the approval of the local Ethics committee, and researchers analyzed anonymized data. The study was registered in ClinicalTrials.gov (NCT04411459).

Data were collected in an electronic case report form developed by YGHEA, CRO division of Ecol Studio SPA (Bologna Operational Headquarters) , and hosted by Actide Nubilaria (Novara, Italy). Collected data comprised demographic data, comorbidities, and ICU-related variables such as SAPS II and SOFA score at ICU admission and the severity of ARDS following the Berlin criteria [12] . Social related variables recorded were: marital status, occupational status, and instruction degree. Residual symptoms at 1 year evaluated were: dyspnea measured with the modified Medical Research Council (mMRC) scale [13] , arthromyalgia, palpitations, cough.

Follow-up pulmonary function tests (PFTs) data were collected for tests performed between 9 and 12 months after ICU discharge and consisted of: forced expiratory volume in the first second (FEV1), FVC, FEV1/FVC ratio, and DLCO. Follow-up CT scan data between 9 and 12 months after ICU discharge were also recorded, in particular concerning the presence of signs of fibrosis, residual ground glass, and atelectasis.

HRQoL was measured at 3 months and 1 year by telephonic interview and administration of the 15D instrument.

The 15D instrument (http://www.15d-instrument.net/15d/) is a generic, 15-dimensional multiutility instrument assessing different aspects of HRQoL (including mobility, vision, hearing, breathing, sleeping, eating, speech, excretion, usual activities, mental function, discomfort and symptoms, depression, distress, vitality, and sexual activity).

The respondents are required to answer about their state of health at the moment of the interview, and each answer is scored on a 5 points scale, with 1 being the best value and 5 being the worst [14] . The valuation system is based on an application of the multiattribute utility theory. The single index score (15D score), representing the overall HRQoL on a 0-1 scale (1=full health, 0=being dead) and the dimension level values, reflecting the goodness of the levels relative to no problems on the dimension (=1) and to being dead (=0), are calculated from the health state descriptive system using a set of population-based preference or utility weights.

Mean dimension level values are used to draw 15D profiles for groups. The minimum clinically important change or difference in the 15D score has been estimated to be ± 0.015 on the basis that people can on average feel such a difference [15] .

The 15D scores are shown to be highly reliable, sensitive and responsive to change, there is a considerable degree of agreement between health state valuations from several European countries [11] , and these latter are generalizable at least in Western-type societies [14, 16] .

Continuous variables were reported as mean and standard deviation (SD) or median and interquartile range (IQR) when appropriate, and comparisons were performed with Student's t-test or Mann-Whitney U test when J o u r n a l P r e -p r o o f appropriate. One way analysis of variance (ANOVA) and Kruskal-Wallis ANOVA were used to compare means and medians in the multiple HRQoL clusters. Categorical variables were expressed as numbers and percentages and compared using the Chi-square test or Fisher exact test. To evaluate the modifications of the 15D score and its dimensions throughout the study period, a paired-samples t-test was adopted.

Since the 15D value can be calculated only if all the 15 dimensions values are available, patients with three or more missing dimensions were excluded from the analysis while a multiple imputation technique was adopted when information regarding less than three dimensions was missing [14] . Two-Step cluster analysis, an exploratory tool designed to reveal natural groupings (clusters) within a dataset, was adopted to define the presence and number of HRQoL profile clusters one year after ICU discharge. This selection procedure allows both categorical and standardized continuous variables and it is fairly robust even when the assumptions of independence and normality of distribution of the variables are violated [17] . To determine the best number of clusters automatically, the indicators BIC (Schwarz's Bayesian Information Criterion) or AIC (Akaike's Information Criterion) are calculated for each number of clusters from a specified range. Automatic clustering was performed based on the 15 dimensions of the 15D instrument, the distance between clusters was calculated using the log-likelihood distance and both AIC and BIC were used to evaluate the best number of clusters.

A linear mixed-effects model was used to determine which factors predicted changes in HRQoL between 3 months and 1 year after ICU discharge using different demographic, ICU-related, and social variables as fixed effects and subject-level random effects while 15D score measured at 3 months and 1 year was the dependent variable. Fixed effects associated with 15D score change with a p < 0.2 in the univariate analyses were introduced in the final multivariate model. Finally, to assess eventual correlations between the 15D breathing score and mMRC scale with pulmonary function tests and ICU stay-related variables, bivariate Spearman correlation tests were performed. P-values < 0.05 were considered statistically significant, and all tests were two-sided. All data were analyzed 

Of the 470 patients initially enrolled in the study, 282 (60%) were discharged alive from the ICU and 205 (43.6%) were followed up at 3 months after ICU discharge.

Four (2%) out of the 205 patients who were initially interviewed, died between three months and 1 year after ICU discharge, while another 23 (11.2%) were lost to follow-up, therefore, the final follow-up cohort was represented by 178 patients (86.8%). The main characteristics of the study population are described in Table 1 , while supplement figure 1 shows the flow of patients throughout the study.

Globally, patients were mostly males (n=129, 72.5%), and had a median age of 64 years (IQR 55-70), the median number of comorbidities was one (IQR 0 -1) and the most frequent were hypertension (n=88, 49.4%), J o u r n a l P r e -p r o o f and diabetes (n=28, 15.7%). The median length of invasive mechanical ventilation was [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] days, and the median ICU stay was [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] days.

The multiple imputation technique was adopted only for replacing missing values in the sleeping (n = 2, 1.1%) and sexual activity (n = 37, 20.8%) dimensions. Table 2 shows the mean values for each of the 15D dimensions at 1 year. Globally, the 15D score significantly increased during the study period, with mean values increasing from 0.857 to 0.880 (p = 0.006). The change was clinically important since it is above the threshold of ± 0.015 [15] .

Most of the increase in HRQoL between 3 months and 1 year was related to physical dimensions (mobility, eating, excretion, usual activities, and sexual activity), while the psychological and breathing dimensions were substantially unchanged ( Figure 1 ).

The two-step cluster analysis identified three different clusters based on BIC and four different clusters based on AIC (see supplement tables 1 and 2). Based on the clinical relevance of the clustering process, AIC-based clustering was chosen. Figure 2 shows the 15D dimensions values among the four different clusters, 84 patients (47.2%) were classified in the "full recovery" group (15D score = 0.964 ± 0.033), and 9 patients (5.1%) were classified in the "bad recovery" group (15D score = 0.572 ± 0.112).

Finally, two different clusters of intermediate recovery were distinguished: 68 patients (38.2%) were classified in a cluster of "partial recovery with mental dimensions most affected" while 17 patients (9.6%) were grouped in a cluster of "partial recovery with physical dimensions most affected". These two intermediate clusters did not significantly differ in terms of absolute 15D score, but only in terms of 15D profiles. The partial recoveryphysical cluster indeed demonstrated significantly lower mean values in the hearing and excretion dimensions but higher values for sleep, discomfort, depression, and distress values. Patients grouped into the four clusters demonstrated different characteristics concerning age, the median number of comorbidities and prevalence of hypertension, and diabetes (Table 2) , with those in the partial and bad recovery clusters being older and affected by more comorbidities.

Most of the patients reported persistent symptoms 1 year after ICU discharge, of which dyspnea (n=104, 58.4%) and arthromyalgia (n=62, 34.8%) were the most frequent. The prevalence of symptoms and the severity of dyspnea measured with the mMRC scale demonstrated significant differences in distribution across the four clusters, with the bad recovery and partial recovery demonstrating a significantly higher prevalence of persistent symptoms and more severe dyspnea.

Univariate linear mixed-effects model analyses selected sex, age, BMI, number of comorbidities, duration of IMV, and tracheostomy for introduction in the multivariate model. The final multivariate model demonstrated that the male gender was associated with a higher increase in HRQoL from 3 to 12 months after ICU discharge while the increase in the duration of IMV and number of comorbidities were negatively associated with HRQoL change (Table 3) .

Pulmonary function tests at 1 year were available for 68 patients ( Table 4 ), 35 of these (51.5%) had a reduction in DLCO values (DLCO lower than 80% of the predicted), while 12 out of 68 (17.6%) showed a restrictive J o u r n a l P r e -p r o o f ventilatory defect (FVC lower than 80% of the predicted). Only 2 patients (0.5%) showed an obstructive ventilatory impairment (FEV1/FVC < 0.7).

Median PaO2/FiO2 ratio nadir during the first five days of IMV was significantly lower in patients showing restrictive ventilatory abnormality (81, IQR 69 -108 vs 112 IQR 83 -150, p = 0.020), while this difference was not observed about DLCO. On the other hand, DLCO impairment was more frequent in the subgroup with partial or bad recovery compared with the subgroup reporting good recovery (Table 4) , and the mean 15D score was significantly lower among those patients showing DLCO impairment (0.881± 0.105 vs 0.933± 0.075, p=0.022).

Finally, bivariate Spearman correlations showed that the 15D breathing dimension at one year after ICU discharge was significantly correlated with both DLCO and duration of IMV, even if the strength of these associations was low (ρ = 0.244 and -0.162, respectively), while the degree of dyspnea measured with mMRC scale was significantly correlated only with the duration of IMV (ρ = 0.204).

Only 37 out of 178 patients had a chest CT scan performed within 9 months and 1 year after ICU discharge, therefore data are only reported in supplement table 3, pulmonary fibrosis signs were observed in 26 (70.3%) of the available CT scans, while non-fibrosing signs such as persistent ground-glass opacities or consolidations in 15 (40.5%) patients.

Approximately half of the mechanically ventilated ICU patients subsequently develop PICS [18] , a multidimensional syndrome that concerns both the physical and psychological traits of the survivors. Signs of persistent impairment in these traits could be found up to two years after ICU discharge among ARDS patients [19, 20] . The COVID-19 pandemic produced waves of ICU survivors at risk of developing both PICS and long-COVID syndrome, moreover, these patients experienced ICU stay in a healthcare system sustaining profound structural and organizational changes in response to the pandemic.

In this cohort of COVID-19 survivors who underwent ICU admission and IMV, we found that: a) HRQoL significantly increased from 3 months to 1 year after ICU discharge; b) this increase was mainly due to an increase in physical dimensions scores, while mental and breathing dimensions scores remained substantially unchanged; c) four main clusters of HRQoL profiles could be identified, two groups at the extremities of the sample showing good and bad recovery and two distinct groups with partial recovery, a larger one showing more severe alterations in mental dimensions, and a smaller group with alterations in physical dimensions; d) factors influencing HRQoL trajectories between 3 months and 1 year were: sex, duration of IMV and number of comorbidities; e) dyspnea remains the most reported symptom in ICU survivors and it is only partially explained by pulmonary function tests.

Due to the lack of long-term information about COVID-19 sequelae, there is currently no consensus about follow-up measures after hospital discharge for COVID-19 patients [21] . Moreover, the follow-up of these J o u r n a l P r e -p r o o f patients could be furtherly complicated by the health systems overload due to the repeated pandemic waves and by the need to recover ordinary activity for non-COVID diseases [22] . A recent online survey evidenced both physical and mental sequelae one year after COVID-19 syndrome for both survivors and their relatives, with age, sex, distance from COVID-19 diagnosis, and length of hospital stay being significant predictors of HRQoL impairment [23] . This is the first literature report specifically focusing on ICU survivors who required IMV and our results strengthen the information that more than half of C-ARDS survivors report significant impairment in HRQoL and persistent symptoms, in particular dyspnea, at 1 year after ICU discharge. In fact, despite a significant increase in HRQoL from 3 months to 1 year after ICU discharge being evidenced, only 47.2% of the patients showed a complete recovery at 1 year after ICU discharge with 15D values comparable to those of the general population [11, 24] .

We have to underline that only physical functioning dimensions showed a significant improvement in this time frame, while mental dimensions remained substantially unchanged. This aspect is confirmed by previous literature about ARDS survivors demonstrating a slower recovery in psychological dimensions [25, 26] .

Therefore, a significant proportion of patients (38.2%) with partial recovery showed a long-lasting impairment in mental dimensions. Specific COVID-19 ICU policies about isolation and visits from relatives, in particular for patients admitted during the first wave, could have had a role in worsening this aspect, and this should probably demand different visiting policies and psychological support during hospital stay [27] .

Several factors influenced the slope of HRQoL trajectories. In our previous study at 3 months after ICU discharge [11] , we found that increasing age resulted negatively associated with HRQoL scores. Indeed, in the current follow-up, the "bad recovery" cluster showed a higher median age (Table 1) . Older age, however, did not significantly influence the entity of HRQoL change towards the study period, which was instead significantly affected by the number of comorbidities and duration of IMV, like in available follow-up literature for other diseases [28] . Male sex was associated with a greater increase in HRQoL over time, this is possibly related to gender-based different shapes of the HRQoL trajectories already reported in the literature [29] .

The absence of association between young age and HRQoL trajectories is another clinically relevant aspect because it underlines that even young patients, often missed at follow-up, could have significantly hampered recovery trajectories and could benefit from a longer follow-up and specific interventions. Finally, the impairment in mental dimensions has been demonstrated to be significant and long-lasting [30] , therefore, psychological advice may be considered when planning the follow-up. Moreover, a not negligible fraction of patients reported dyspnea and a variable degree of impairment in pulmonary function tests that should be furtherly discussed.

We reported a very high prevalence (58.4%) of persistent dyspnea at 1 year after ICU discharge in comparison to that reported from papers, this phenomenon could be explained by the fact that our population was made by critically ill patients (ICU survivors) and not by a mixed population of critical and non-critical COVID-19 patients [3] . Analogously, the 15-D breathing dimension score significantly contributed to the reduction in HRQoL at 3 months, and it did not significantly improve during the study period (see Figure 1 and Table 2 ).

Most of the survivors showed some kind of impairment in both DLCO (51.5%) and FVC (17.6%) values, in line with previous literature about ARDS survivors patients [31] . In particular, in our study, reduction of FVC was significantly more frequent in patients experiencing more severe forms of ARDS (Table 4 ). According to actual literature [3, 30] , respiratory functional impairment after COVID-19 pneumonia is associated with persistent pulmonary radiological abnormalities so that PFTs may therefore represent a marker of radiological sequelae of severe COVID-19.

On the other hand, in our study population DLCO impairment degree, which is the most reported respiratory functional impairment in post-COVID syndrome, was not associated with ARDS severity. Even if potentially controversial, this result confirms that COVID-19 patients requiring intensive care and IMV can fully recover [32, 33] .

Interestingly, the presence and severity of dyspnea, measured by the mMRC scale, was correlated with IMV length and not with DLCO and FVC values, while the 15D breathing dimension was correlated to both IMV length and DLCO (Table 4 ).

According to our results, pulmonary function tests alterations do not completely explain the persistence of dyspnea, so that the perceived breathing impairment complained by patients after COVID-19 may not be fully related to pulmonary function impairment [34, 35] . Other clinical and neuro-psychological aspects such as experiences of anxiety and pain, the sedative regimens, and mechanical ventilation settings during ICU stay [36, 37] , could play a role in determining the presence and entity of dyspnea among survivors.

Dyspnea is defined as a subjective perception of uncomfortable breathing and could be determined by complex and multiple mechanisms including social, psychological, and physical conditions [38] , and it is one of the major components of the post-COVID syndrome [39] . This study evaluated a population made of critically ill patients, which could be also affected by the so-called post-intensive care syndrome (PICS), a multisystemic syndrome characterized by new or worsening physical, mental and neurocognitive disorders [40, 41] , that could overlap post-COVID syndrome manifestations.

It is impossible to precisely estimate the proportion of mental and physical components potentially responsible for dyspnea as well as the relative roles of PICS and post-COVID syndrome from the available data in our population. However, the very high prevalence of psychological consequences, together with the only partial concordance between dyspnea and pulmonary functional impairment (especially DLCO), may suggest a cardinal role of neuropsychological mechanisms at the root of this symptom.

Exploring this topic in larger studies could be helpful to better understand an important aspect of both PICS and post-COVID syndrome, and to foster their management in a multidisciplinary approach.

Several limitations should be discussed: first, we adopted the telephonic interview rather than follow-up visits, this was chosen to reach the maximum number of patients taking into account the travel restrictions due to the second and third waves of the COVID-19 pandemic. Other limitations are also represented by the small proportion of patients with available PFTs at 1-year, so that the prevalence of functional impairment could be biased. Moreover, the lack of association with lung radiological data at 1-year due to data unavailability and the absence of baseline PFTs information prevented us from deriving clinical-radiological correlations and carrying out pre-post analyses.

Missing data for the sexual activity dimension exceeded 20%, despite this aspect is in line with our previous follow-up [11] and data substitution was obtained with the multiple imputation technique, cautious interpretation of these data should be warranted.

Finally, among ICU survivors significant modifications in HRQoL dimensions were reported up to ten years after discharge [42] , therefore, longer follow-up times could help in better defining the HRQoL trajectories and the effects of eventual interventions [43] . (10 -27) 14 (9 -26) 18 (13 -26) 18 (9 -32) 20 (11 -33) 0.499

Length of ICU stay -d (IQR) 23 (15 -35) 19 (14 -35) 24 (17 -36) 25 (13 - Vision -mean ± SD 0.953 ± 0.119 0.992 ± 0.040 0.942 ± 0.108 0.949 ± 0.094 0.681 ± 0.280 < 0.001

Hearing -mean ± SD 0.968 ± 0.098 1.000 ± 0.000 1.000 ± 0.000 0.745 ± 0.135 0.857 ± 0.192 < 0.001

Breathing -mean ± SD 0.746 ± 0.238 0.879 ± 0.154 0.620 ± 0.227 0.753 ± 0.223 0.438 ± 0.238 < 0.001

Sleeping -mean ± SD 0.838 ± 0.238 0.940 ± 0.135 0.716 ± 0.274 0.929 ± 0.142 0.632 ± 0.312 < 0.001

Eating -mean ± SD 0.979 ± 0.102 1.000 ± 0.000 1 .000 ± 0.000 1.000 ± 0.000 0.587 ± 0.221 < 0.001

Speech -mean ± SD 0.980 ± 0.090 0.996 ± 0.032 0.996 ± 0.036 0.948 ± 0.117 0.777 ± 0.276 < 0.001

Excretion -mean ± SD 0.974 ± 0.110 1.000 ± 0.000 1.000 ± 0.000 0.872 ± 0.191 0.720 ± 0.292 < 0.001

Usual activities -mean ± SD 0. Dyspnoea mMRC scale -(IQR) 0 (0 -2) 0 (0 -0) 1 (0 -2) Abbreviations: BMIbody mass index; ARDSacute respiratory distress syndrome; SAPSsimplified acute physiology score; SOFAsequential organ failure assessment; ICUintensive care unit.

Notes: significant p values are evidenced in bold;‫٭‬reference class: severe ARDS; * reference class: employed J o u r n a l P r e -p r o o f (25 -30) 28 (25 -31) 28 (25 -29) 0.598 28 (26 -37) 28 (25 -30) 0.292 COPD -no (%) 6 (8.8%) 2 (5.7%) 4 (12.1%) 0.667 0 (0%) 6 (10.7%) 0.581

Number of comorbidities -

Intensive Care Unit stay related variables SAPS II score (IQR) 35 (27 -43) 38 (29 -49) 34 (27 -42) 0.161 35 (27 -40) 35 (29 -44) (8 -21) 15 (9 -28) 14 (8 -18) 0.140 16 (7 -23) 14 (9 -21) 0.942

Length of ICU stay -d (IQR) 19 (14 -33) 22 (15 -42) 19 (13 -29) 0.100 32 (15 -42) 19 (14 -29) Dyspnoea mMRC scale -(IQR) 0 (0 -2) 0 (0 -0) 1 (0 -2) 1 (0 -2) 3 (2 -4) < 0.001 

COVID-19) Dashboard | WHO Coronavirus (COVID-19) Dashboard With Vaccination Data

Resource allocation in the Covid-19 health crisis: are Covid-19 preventive measures consistent with the Rule of Rescue?

3-month, 6-month, 9-month, and 12-month respiratory outcomes in patients following COVID-19-related hospitalisation: a prospective study

6-month neurological and psychiatric outcomes in 236 379 survivors of COVID-19: a retrospective cohort study using electronic health records, The Lancet Psychiatry

Cardiac sequelae after coronavirus disease 2019 recovery: a systematic review

Long-COVID': A cross-sectional study of persisting symptoms, biomarker and imaging abnormalities following hospitalisation for COVID-19

Cao, 6-month consequences of COVID-19 in patients discharged from hospital: a cohort study

Post-COVID syndrome in non-hospitalised patients with COVID-19: a longitudinal prospective cohort study

Post-intensive care syndrome after a critical COVID-19: cohort study from a J o u r n a l P r e -p r o o f Belgian follow-up clinic

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

Quality of life of COVID-19 critically ill survivors after ICU discharge: 90 days follow-up

The Berlin definition of ARDS: An expanded rationale, justification, and supplementary material

Evaluation of clinical methods for rating dyspnea

The 15D instrument of health-related quality of life: Properties and applications

Estimating the minimum important change in the 15D scores

Measuring healthrelated quality of life in women on hormone replacement therapy, Expert Rev

Using Two-Step Cluster Analysis and Latent Class Cluster Analysis to Classify the Cognitive Heterogeneity of Cross-Diagnostic Psychiatric Inpatients

Post-Intensive care syndrome and COVID-19: crisis after a crisis?

Two-year cognitive, emotional, and quality-of-life outcomes in acute respiratory distress syndrome

Physical Complications in Acute Lung Injury Survivors: A 2-Year Longitudinal Prospective Study

COVID-19 follow-up planning: what will we be missing?

The Untold Toll -The Pandemic's Effects on Patients without Covid-19

Measuring the impact of COVID-19 on the quality of life of the survivors, partners and family members: a cross-sectional international online survey

Lung function and quality of life in survivors of the acute respiratory distress syndrome (ARDS)

Health-related quality of life in intensive care survivors: Associations with social support, comorbidity, and pain interference

Impact of different visiting policies on family satisfaction in two Spanish ICUs before and during COVID-19

Hemodialysis Study Group, The influence of age on changes in health-related quality of life over three years in a cohort undergoing hemodialysis

Gender differences in the trajectories and the risk factors of depressive later life

Respiratory and Psychophysical Sequelae Among Patients With COVID-19 Four Months After Hospital Discharge

One-Year Outcomes in Survivors of the Acute Respiratory Distress Syndrome

Six months follow-up of patients with invasive mechanical ventilation due to covid-19 related ards

Short-term health-related quality of life, physical function and psychological consequences of severe COVID-19

Persistent poor health after covid-19 is not associated with respiratory complications or initial disease severity

Integrative respiratory follow-up of severe COVID-19 reveals common functional and lung imaging sequelae

Mechanisms, assessment, and management: a consensus statement

Phenotyping long COVID

Postintensive Care Syndrome: Right Care, Right Now…and Later

European Academy of Rehabilitation Medicine. Covid-19 and Post Intensive Care Syndrome: A Call for Action

Health-related quality of life in ICU survivors-10 years later

Fatigue in intensive care survivors one year after discharge, Health Qual

Worst PaO2/FiO2 5d -(IQR) 105

ARDS -acute respiratory distress syndrome; SAPS -simplified acute physiology score; SOFA -sequential organ failure assessment; ICUintensive care unit; mMRC modified Medical Research Council; CRS Respiratory system compliance; DLCOdiffusion capacity of lungs for carbon monoxide, FVC -forced vital capacity; FEV1 -forced expiratory volume in the first second; MV -mechanical ventilation Notes: significant p values are evidenced in bold, as well as the rho (ρ) values of significant bivariate Spearman correlations

Health-related quality of life profiles, trajectories, persistent symptoms and pulmonary function one year after ICU discharge in invasively ventilated COVID-19 patients, a prospective follow-up study