Abstract
Rationale Post-coronavirus disease 2019 (COVID-19) survivors frequently have dyspnoea that can lead to exercise intolerance and lower quality of life. Despite recent advances, the pathophysiological mechanisms of exercise intolerance in the post-COVID-19 patients remain incompletely characterised. The objectives of the present study were to clarify the mechanisms of exercise intolerance in post-COVID-19 survivors after hospitalisation.
Methods This prospective study evaluated consecutive patients previously hospitalised due to moderate-to-severe/critical COVID-19. Within mean±sd 90±10 days of onset of acute COVID-19 symptoms, patients underwent a comprehensive cardiopulmonary assessment, including cardiopulmonary exercise testing with earlobe arterialised capillary blood gas analysis.
Measurements and main results 87 patients were evaluated; mean±sd peak oxygen consumption was 19.5±5.0 mL·kg−1·min−1, and the tertiles were ≤17.0, 17.1–22.2 and ≥22.3 mL·kg−1·min−1. Hospitalisation severity was similar among the three groups; however, at the follow-up visit, patients with peak oxygen consumption ≤17.0 mL·kg−1·min−1 reported a greater sensation of dyspnoea, along with indices of impaired pulmonary function, and abnormal ventilatory, gas-exchange and metabolic responses during exercise compared to patients with peak oxygen consumption >17 mL·kg−1·min−1. By multivariate logistic regression analysis (receiver operating characteristic curve analysis) adjusted for age, sex and prior pulmonary embolism, a peak dead space fraction of tidal volume ≥29 and a resting forced vital capacity ≤80% predicted were independent predictors of reduced peak oxygen consumption.
Conclusions Exercise intolerance in the post-COVID-19 survivors was related to a high dead space fraction of tidal volume at peak exercise and a decreased resting forced vital capacity, suggesting that both pulmonary microcirculation injury and ventilatory impairment could influence aerobic capacity in this patient population.
Abstract
Post-COVID-19 survivors may have exercise intolerance; in this study, this was related to high VD/VT at exercise and decreased FVC % pred, suggesting that pulmonary microcirculatory injury and ventilatory impairment influence aerobic capacity https://bit.ly/41AYvwl
Introduction
In March 2020, coronavirus disease 2019 (COVID-19) was characterised by the World Health Organization as a pandemic infection and has been considered an international public health emergency for the past 2 years. A few months after the pandemic's start, Brazil had the second highest number of confirmed COVID-19 cases worldwide. In April 2021, Brazil had become the epicentre of the COVID-19 pandemic, with >4000 deaths per day [1].
COVID-19 infection may be asymptomatic in the acute phase, but clinical presentation might also range from mild respiratory symptoms to severe respiratory failure with associated acute respiratory distress syndrome (ARDS). Additionally, clinical presentation might include extrapulmonary symptoms [2]. After hospitalisation, patients may remain symptomatic and this could be related to cardiac/lung sequalae and/or post-COVID-19 syndrome [3].
The post-COVID-19 syndrome is defined by the presence of persistent symptoms 12 weeks after the onset of COVID-19, and is not attributable to other known causes [3]. Among the most frequent signs and symptoms reported in post-COVID-19 syndrome are fatigue, muscle weakness, dyspnoea, hypoxaemia, depression, anxiety and sleep and cognitive disorders, along with exercise intolerance [3–5], the last of which might lead to a significant decrease in functional capacity and quality of life. Different hypotheses for mechanisms of exercise intolerance after COVID-19 infection have been explored so far, and physical deconditioning has been described as one of the most likely driving forces of symptoms [6, 7], despite the complexity of COVID-19 and the potential for multiorgan involvement.
In this context, recent findings suggest that exercise limitation in post-COVID-19 survivors in more severe patients may be related to 1) central cardiocirculatory disorder due to chronic myocardial inflammation and/or pulmonary microvascular injury [8]; 2) ventilatory inefficiency [9, 10] due to increased dead space (VD) as a fraction of tidal volume (VT), possibly related to endothelial and/or microvascular dysfunction [11]; 3) reduced peripheral muscle oxygen extraction [11, 12]. In patients with mild post-COVID-19 syndrome, dysfunctional breathing was a relevant mechanism of exercise intolerance [12]. Nevertheless, despite these recent advances, the pathophysiological mechanisms of exercise intolerance in post-COVID-19 survivors remain incompletely characterised. In the current study, we aimed to clarify the mechanisms of exercise intolerance associated with reduced aerobic capacity after moderate-to-severe/critical COVID-19 hospitalisation.
Materials and methods
Study design and participants
The current study is part of an observational prospective Brazilian initiative to evaluate clinical symptoms and respiratory, radiological and metabolomic function in patients who were hospitalised due to COVID-19 (FENIX Study; Brazilian Clinical Trials Registry ReBEC identifier RBR-8j9kqy).
The current report presents data from consecutive adult patients from the post-COVID-19 outpatient clinic of the Federal University of São Paulo. All included patients had the first medical visit after hospital discharge between August 2020 and May 2021 and had the following characteristics at the time of COVID-19 hospitalisation (inclusion criteria): 1) confirmed diagnosis of COVID-19 by reverse transcription PCR; 2) received supplemental oxygen (O2) support; and 3) had acute lung parenchymal involvement confirmed by chest computed tomography (CT) scan.
Patients were invited to participate in the study in their first clinical outpatient evaluation after hospital discharge. Those patients who fulfilled the study inclusion criteria and signed an informed consent form had their clinical information recorded, and within 90±10 days after the onset of COVID-19 acute symptoms, performed a comprehensive cardiopulmonary assessment, including a cardiopulmonary exercise testing (CPET) with earlobe arterialised capillary blood gas analysis. All other tests, pulmonary lung function, echocardiogram and high-resolution chest CT (HRCT), were performed within 10 days of CPET (figure 1).
Study protocol and patient inclusion. COVID-19: coronavirus disease 2019; RT: reverse transcription; CT: computed tomography; HRCT: high-resolution CT; CPET: cardiopulmonary exercise testing.
Patients in palliative cancer care, with psychiatric disturbances, musculoskeletal impairment to perform the exercise and uncontrolled known cardiovascular, endocrine–metabolic or renal diseases were excluded from the study. Patients who could not complete the study follow-up visit were also excluded (supplementary figure E1).
The methodological description of pulmonary function tests and modified Medical Council Research (mMRC) dyspnoea scale are described in the supplementary material [13–15].
CPET
Patients performed symptom-limited, ramp-incremental cycle ergometer CPET using a computer-based exercise system with breath-by-breath analysis of metabolic, ventilatory and cardiovascular variables (ULTIMA CardioO2; Med Graphics, Saint Paul, MN, USA). The work rate was individually selected to provide an incremental phase of 7–12 min (5–20 W·min−1) and started after a 2-min unloading warm-up period. The measures obtained was described elsewhere [16] and are included in the supplementary material. Earlobe arterialised capillary blood gas samples (Heparinated 200-I microtubes; Radiometer, Copenhagen, Denmark), were drawn at rest and at peak exercise after applying vasodilator capsaicin cream (Moment 0.075%; Apsen Pharmaceutical, São Paulo, Brazil). The blood analyses were performed immediately (ABL800; Radiometer, Brønshøj, Denmark) to obtain lactate and gas exchange variables (arterial oxygen partial pressure, arterial carbon dioxide partial pressure (PaCO2) and arterial oxygen saturation). Measures of alveolar–arterial O2 gradient, arterial end-expiratory carbon dioxide gradient (PaETCO2) and VD/VT (Enghoff modification of the Bohr equation) were then calculated [16].
Data analysis
In the study design, there were not enough studies for sample calculation; for this sample, the confidence interval was used for a population proportion (95% CI) considering a third of the patients with reduced peak oxygen uptake (V′O2peak). Descriptive statistics are present as mean±sd, median and interquartile range of frequencies. Patients were categorised according to V′O2peak tertiles: ≤17.0 mL·kg−1·min−1, 17.1–22.2 mL·kg−1·min−1 or ≥22.3 mL·kg−1·min−1. Comparisons between more than two groups were performed with one-way ANOVA with Bonferroni or Kruskal–Wallis post hoc analysis, according to the data distribution. Correlation analyses were performed using Pearson's or Spearman's coefficients to identify variables significantly associated with V′O2peak mL·kg−1·min−1. Receiver operating characteristic (ROC) curves were drawn for variables that had a high correlation with V′O2peak while accounting for the presence or absence of a V′O2peak ≤17.0 mL·kg−1·min−1. The thresholds for each ROC curve were obtained from the points with the greatest sum of sensitivity and specificity. After dichotomising the variables of interest according to ROC thresholds, univariate logistic regression was performed to explore potential V′O2peak ≤17.0 mL·kg−1·min−1 predictors. Noncollinear variables (r≥0.6) from the univariate analysis from different pathophysiological domains (i.e. symptoms, lung function, ventilatory, gas-exchange or metabolic responses to exercise) were included in multivariate logistic regression models adjusted for age, sex and prior pulmonary embolism to estimate the probability of having a V′O2peak ≤17.0 mL·kg−1·min−1, a second model was analysed with adjustment for age, sex and the presence of any comorbidity (supplementary table E4). The accepted statistical significance value was <0.050. Graphs were created with GraphPad Prism (version 9.3.0 for Windows; GraphPad Software), and statistical analyses were performed using SPSS for Windows (version 21.0; IBM, Armonk, NY, USA).
Results
96 patients were eligible to participate in this study. Nine patients were excluded. Patient exclusion occurred due to acute arthritis (n=1), severe thrombocytopenia (n=1), acute deep vein thrombosis (n=1), uncontrolled systemic arterial hypertension (n=1), acute metabolic acidosis (n=1) and inability to perform the study follow-up visit (n=4). Therefore, the study sample was composed of 87 patients.
Of the 87 included patients, 54% were admitted to the intensive care unit (ICU) and 49% had ≥50% ground-glass opacities on chest CT scan. The mean age was 53±13 years; 62% were male; and 63% had two or more comorbidities (table 1). Systemic hypertension, previous smoking history and obesity were the most common comorbidities among the patients studied (supplementary table E1). Detailed information regarding patients' comorbidities, medications of continuous use and COVID-19-related acute symptoms are provided in the supplementary table E1).
Baseline characteristics of coronavirus disease 2019 patients
The mean V′O2peak for the entire study sample was 19.5±5.0 mL·kg−1·min−1, corresponding to 93±21% of predicted V′O2 (30% had V′O2peak ≤80% pred). V′O2peak tertiles were ≤17.0, 17.1–22.2 and ≥22.3 mL·kg−1·min−1. Patients with V′O2peak ≤17.0 mL·kg−1·min−1 had similar hospitalisation severity as patients with V′O2peak 17.1–22.2 and ≥22.3 mL·kg−1·min−1, including days in ICU, need for mechanical ventilation and radiological severity on chest CT at admission. However, at the study follow-up visit (90±10 days after the onset of COVID-19), patients with V′O2peak ≤17.0 mL·kg−1·min−1 reported a greater sensation of dyspnoea (mMRC ≥1) compared to the other two groups (table 2). Additionally, patients with V′O2peak ≤17.0 mL·kg−1·min−1 had lower forced vital capacity (FVC), total lung capacity (TLC), diffusing capacity of the lung for carbon monoxide (DLCO) and residual volume compared to the other groups (table 2). The persistence of lung parenchymal involvement on HRCT and cardiac function by echocardiogram at the follow-up visit was similar between groups (table 2).
Coronavirus disease 2019 patients’ characteristics during hospitalisation and lung function tests, chest computed tomography and echocardiogram according to peak oxygen uptake (V′O2peak) tertiles
CPET findings are presented in table 3. Patients with V′O2peak ≤17.0 mL·kg−1·min−1 achieved lower peak work rate (WR), peak heart rate and lower ΔV′O2/ΔWR. At the anaerobic threshold, patients with V′O2peak ≤17.0 mL·kg−1·min−1 had higher minute ventilation (V′E)/carbon dioxide production (V′CO2) and lower PETCO2 and no difference on V′O2 (table 3). Additionally, patients with V′O2peak ≤17.0 mL·kg−1·min−1 had higher ΔV′E/ΔV′CO2 at respiratory compensation point (RCP), peak respiratory rate/VT, peak VD/VT, peak PaETCO2 and associated with a lower peak arterial oxygen content (CaO2) and higher level of lactate/WR and a greater sensation of dyspnoea and fatigue in proportion to WR compared to patients with V′O2peak 17.1–22.2 and ≥22.3 mL·kg−1·min−1 (figure 2).
Cardiopulmonary exercise testing (CPET) responses and blood gas analysis of coronavirus disease 2019 patients at rest and at peak exercise according to peak oxygen uptake (V′O2peak; mL·kg−1·min−1) tertiles
Comparison of peak oxygen uptake (V′O2peak) (mL·kg−1·min−1) in cardiopulmonary exercise testing responses after 3 months of symptoms in survivors of coronavirus disease 2019. a) Ventilatory equivalents for carbon dioxide at respiratory compensation point (RCP); b) respiratory rate (RR) of tidal volume (VT) at peak exercise; c) dead space volume (VD) fraction of VT at peak exercise; d) arterial to end-tidal carbon dioxide difference at peak exercise (PaETCO2peak); e) relationship of V′O2 and arterial oxygen content at peak exercise (CaO2peak); f) lactate by work rate (WR) at peak exercise; g) dyspnoea Borg scale by WR at peak exercise; h) fatigue Borg scale by WR at peak exercise. V′E: minute ventilation; V′CO2: carbon dioxide production. p-values calculated by ANOVA or Kruskal–Wallis.
There was a positive correlation between V′O2peak and FVC, DLCO, V′E/maximal voluntary ventilation (MVV) and peak CaO2. There was a negative correlation between V′O2peak, several comorbidities, dyspnoea (mMRC), ΔV′E/ΔV′CO2RCP, peak respiratory rate/VT, peak VD/VT, peak PaETCO2 and peak lactate/WR. No correlation was found between V′O2peak and days of hospitalisation or days in ICU (supplementary table E2).
The ROC curve analyses to identify the presence of a V′O2peak ≤17.0 mL·kg−1·min−1 showed a statistically significant area under the curve for symptoms (mMRC), FVC, DLCO, peak respiratory rate/VT, peak V′E/MVV, peak VD/VT, ΔV′E/ΔV′CO2RCP, PaETCO2, peak WR, peak CaO2, peak lactate and peak lactate/WR (supplementary table E3).
The univariate logistic regression analysis to predict a V′O2peak ≤17.0 mL·kg−1·min−1, including relevant variables from different pathophysiological domains (i.e. symptoms, lung function, ventilatory, gas-exchange or metabolic responses to exercise) is presented in table 4. Among noncollinear variables, the multivariate logistic regression model adjusted for age, sex and presence of pulmonary embolism identified that a FVC ≤80% pred and a peak VD/VT ≥29 were independent predictors of a V′O2peak ≤17.0 mL·kg−1·min−1 (table 4). A second multivariate logistic regression model was performed, with adjustment for age, sex and the presence of any comorbidity and FVC % pred and VD/VT remained as predictors of V′O2peak (supplementary table E4). Of note, VD/VT had a negative correlation with DLCO % pred (r=0.64, p<0.01); a positive correlation with peak VD (r=0.62, p<0.001); and a positive correlation with PaETCO2 (r=0.88, p<0.001). Interestingly, FVC and VD/VT were not significantly correlated (r=0.14, p=0.292).
Univariate and multivariate logistic analysis adjusted for sex, age and prior pulmonary embolism for peak oxygen uptake (V′O2peak) ≤17.0 mL·kg−1·min−1 according to persistence of symptoms, lung function and cardiopulmonary exercise testing (CPET) variables
Discussion
The present observational study showed that exercise intolerance in post-COVID-19 survivors with a relatively short hospital stay (15±10 days) was related to high VD/VT at peak exercise and low FVC % pred after 90±10 days of acute infection. This finding suggests that both pulmonary microcirculation injury and pulmonary ventilatory impairment might play a role in influencing aerobic capacity in the post-COVID-19 survivors.
VD/VT is related to the physiological dead space ratio, divided into anatomical dead space (i.e. airways that do not participate in gas exchange), and alveolar dead space. A high VD/VT results from areas of normal ventilation and low perfusion that contribute to ventilation–perfusion (V′/Q′) mismatch. A low VD/VT results from areas of low ventilation and normal perfusion. Both high and low VD/VT can be present in the same disease [17]. It is important to note that VD/VT is expected to reach a level <0.20 after the anaerobic threshold in physiological conditions due to the increased perfusion of areas of the lungs with high V′/Q′ ratios at rest and a relatively greater increase in VT tidal volume than anatomical dead space, the abnormal response is dependent on severity of pulmonary lesions [16]. In our sample, VD/VT decreased during exercise in all three groups, but peak VD/VT progressively increased from the subgroup V′O2peak >22.2 mL·kg−1·min−1 to the subgroup V′O2peak ≤17.0 mL·kg−1·min−1. Additionally, despite reducing during exercise, VD/VT did not reach physiological values in all three groups. A high VD/VT might be related to ventilatory inefficiency (high V′E/V′CO2), and dyspnoea sensation, being associated or not with enhanced chemosensitivity and a decreased carbon dioxide set point [17].
Our results show that a high VD/VT at peak exercise (≥0.29) is an independent predictor of a V′O2peak ≤17.0 mL·min·kg−1 (table 4). In addition to the high VD/VT, a high peak exercise PaETCO2 (figure 2) might corroborate the presence of V′/Q′ inequality in the studied population. Some studies in post-COVID-19 patients showed an increase in VD/VT; however, they did not link its association to patients’ exercise intolerance [11, 18]. Baratto et al. [11] showed that exercise hyperventilation after COVID-19 acute infection was related to enhanced chemoreflex sensitivity rather than increased VD/VT. Conversely, others have demonstrated that a reduced V′O2peak was associated with a mild increase of V′E/V′CO2 and have suggested that the observed hyperventilation could be related to increased chemoreflex sensitivity secondary to deconditioning, dysfunctional breathing or even dysautonomia [6, 9, 19, 20]. Acute COVID-19 lung lesions have been related to diffuse alveolar damage, interstitial fibrosis and endothelial vascular injuries, which result in areas of shunt (low V′/Q′) and/or dead space (high V′/Q′). Along the lines, studies comparing ARDS in COVID-19 versus non-COVID-19 patients showed that COVID-19 ARDS patients have a higher dead space ventilation compared to non-COVID-19-ARDS, despite a similar pulmonary compliance [21]. The aforementioned lung insults can potentially cause transitory or persistent lung sequelae [22–24]. In our study, VD/VT had a negative correlation with low DLCO, a positive correlation with VD and a positive correlation with PaETCO2. Similar findings have been shown for cardiocirculatory diseases such as left heart failure and pulmonary arterial hypertension [25–28]. It is important to note that a low DLCO was found in the long term after severe acute respiratory syndrome coronavirus 1 (SARS-CoV-1) and SARS-CoV-2 patients, even in those with normal lung parenchyma on HRCT [2, 29–33]. Furthermore, during acute COVID-19 infection, dual-energy thoracic CT studies showed the presence of pulmonary perfusion heterogeneity along with pulmonary ischaemic areas in the absence of visible pulmonary arterial thrombosis and in areas not related to ground-glass opacities or any parenchymal lesions, which may reflect the presence of microvascular injury [34]. Based on this and our study findings, we speculate that chronic lung microvascular injury might be a pathophysiological mechanism leading to high VD/VT during exercise in post-COVID-19 patients. This hypothesis is supported by the multivariate regression (table 4), where history of pulmonary embolism was not a determining factor for the increased VD/VT. The same occurs when the regression is adjusted for the presence of any comorbidity (supplementary table E4), suggesting that VD/VT might be elevated due to microcirculation injury. Of note, this microvascular involvement had no repercussions on the findings of resting echocardiogram in our patients.
FVC % pred was also identified as an independent predictor of a V′O2peak ≤17.0 mL·kg−1·min−1; however, FVC and VD/VT were not significantly correlated. A low FVC has been reported in post-COVID-19 patients as far as 1 year after the acute infection, and similar results have been demonstrated in SARS-CoV-1 survivors [29, 33]. Considering that a low FVC might be related to the ARDS severity, it might indicate the development of restrictive ventilatory impairment secondary to lung interstitial sequelae [35]. This is in line with a tachypnoea pattern, proven by high respiratory rate/VT. Nonetheless, we did not identify significant differences in TLC and acute parenchymal lung involvement on HRCT according to V′O2peak severity (table 2).
In addition to a potential interstitial lung disease development impacting FVC, we should also consider pulmonary neuromuscular dysfunction as a possible cause of reduced FVC. Inspiratory muscle weakness and decreases in peripheral muscle strength have been described in post-COVID-19 patients, and were associated with reduced aerobic capacity [35–37]. However, our results did not identify a significant difference in maximal inspiratory pressure according to V′O2peak severity (table 2).
Interestingly, lactate/WR was higher according to V′O2 tertiles (figure 2), despite the similar anaerobic threshold (table 3). This finding has been demonstrated previously in patients with oxidative myopathy [38]. It suggests that the mechanisms of lactate clearance fail to keep pace with lactate production in post-COVID-19 patients, and/or there is an impairment in O2 utilisation at higher levels of exercise [39]. In our study, the elevated lactate/WR observed in patients with V′O2peak ≤17.0 mL·kg−1·min−1 might be a consequence of a mildly reduced O2 delivery (low CaO2) and/or an imbalance in O2 muscle utilisation due to a decrease in oxidative fibres secondary to prolonged hospitalisation, neuromuscular drug toxicity, direct viral mitochondrial injury by immediate viral effect and/or systemic inflammation [3, 40]. As a result, the aforementioned mechanisms will stimulate a rapid respiratory rate and increase the neural perception of dyspnoea, but further studies are required to investigate this hypothesis in post-COVID-19 patients.
Our study has some limitations that should be considered. We did not include a healthy-control group; nonetheless, patients with V′O2peak >22.2 mL·kg−1·min−1 had a more preserved aerobic capacity and therefore could be considered from an exercise physiology perspective as a control for the subgroup with V′O2peak ≤17.0 mL·kg−1·min−1. Despite not having a healthy-control group, our exercise findings are similar to Skjørten et al. [6]. Along these lines, it is important to note that all patients included in the subgroup V′O2peak 22.2 mL·kg−1·min−1 had a V′O2peak >80% pred, and that most patients with a V′O2 ≤80% pred were included in the subgroup V′O2peak ≤17.0 mL·kg−1·min−1. In physiological terms, the V′O2 in absolute value decreases with ageing, and more so in females than males. In our study, age was different across V′O2peak subgroups. It is known that age and sex might influence some ventilatory responses due to lower VTpeak and less efficient ventilation during exercise (without abnormally high VD/VT), probably related to increased airway resistance and mechanical constraint with a reduced compliance of the lungs. This phenomenon is more pronounced in older females but, in general, with little impact on exercise capacity. Of note, sex per se does not affect gas exchange, but ageing could indeed change the PaCO2 equilibrium [41, 42]. Considering this and aiming to minimise the possible effects of age and sex on exercise physiological responses and in the study findings, the multivariate model was adjusted for age and sex. We did not perform exercise haemodynamics, single-photon emission lung CT or dual-energy CT thoracic angiography, and therefore we can only speculate on the association between high VD/VT during exercise and the hypothesis of pulmonary microvascular dysfunction. Additionally, we did not perform comprehensive muscle-related studies, and therefore we are not able to undoubtedly confirm muscle weakness as a potential cause for a reduced V′O2peak. Finally, the control of breathing during exercise is complex, multifactorial and not completely understood. The current study could not explain or phenotype the pathophysiological mechanisms of exercise intolerance in post-COVID-19 patients.
In summary, the current study demonstrates that a high VD/VT at peak exercise and a low resting FVC are associated with a reduced V′O2peak in moderate-to-severe/critical post-COVID-19 patients. The high peak exercise VD/VT might suggest the role of pulmonary microvascular dysfunction on dyspnoea and exercise intolerance in the post-COVID-19 survivors. The low FVC suggests that pulmonary ventilatory dysfunction might be an additional factor influencing aerobic capacity in this patient population. Further studies are needed to confirm whether post-COVID-19 survivors will develop pulmonary vascular disease and/or clinically relevant interstitial pulmonary disease in the long term.
Supplementary material
Supplementary Material
Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.
Supplementary material 00538-2022.supplement
Acknowledgements
The authors thank all participating investigators of the SEFICE (Pulmonary Function and Clinical Exercise Physiology Sector) from the Hospital Sao Paulo – UNIFESP/EPM for their contribution to the collected data and review of the article.
Footnotes
Provenance: Submitted article, peer reviewed.
Support statement: This study was supported by the Sao Paulo Research Foundation (Fapesp) (protocol number 2020/08996-1) and M.L. Lafetá receives a PhD bursary from CAPES (Coordination for Improvement of Higher Education Personnel) (process number 88887.508806/2020-00). Funding information for this article has been deposited with the Crossref Funder Registry.
Author contributions: M.L. Lafetá, S.E. Tanni, R.K.F. Oliveira and E.V.M. Ferreira were responsible for the data analysis and content of the article. M.L. Lafetá, A.L.P. Albuquerque, S.E. Tanni and E.V.M. Ferreira were responsible for drafting the manuscript. M.L. Lafetá, V.C. Souza, T.C.F. Menezes, P.A. Sperandio, J.P. Carlstron and E.V.M. Ferreira were responsible for data collection. C.G.Y. Verrastro, M. Izbicki, F.J. Mancuso and E.V.M. Ferreira were responsible for image analysis. L.E. Nery, A.L.P. Albuquerque, P.A. Sperandio, R.F.K. Oliveira and E.V.M. Ferreira revised the draft for important intellectual content. M.L. Lafetá, A.L.P. Albuquerque, L.E. Nery, P.A. Sperandio, R.K.F. Oliveira and E.V.M. Ferreira gave final approval of the content of the manuscript.
Conflict of interest: M.L. Lafetá has nothing to disclose.
Conflict of interest: V.C. Souza has nothing to disclose.
Conflict of interest: T.C.F. Menezes has nothing to disclose.
Conflict of interest: C.G.Y. Verrastro has nothing to disclose.
Conflict of interest: F.J. Mancuso has nothing to disclose.
Conflict of interest: A.L.P. Albuquerque has nothing to disclose.
Conflict of interest: S.E. Tanni is President of the Sao Paulo Thoracic Society outside the submitted work.
Conflict of interest: M. Izbicki has nothing to disclose.
Conflict of interest: J.P. Carlstron has nothing to disclose.
Conflict of interest: L.E. Nery has nothing to disclose.
Conflict of interest: R.K.F. Oliveira report grants from the National Council for Scientific and Technological Development (CNPq), Brazil (grant 313284/2021-0), and personal fees from Janssen Brazil, outside the submitted work.
Conflict of interest: P.A. Sperandio has nothing to disclose.
Conflict of interest: E.V.M. Ferreira reports speaker fees from Janssen, and personal fees from Aché, AstraZeneca, Bayer, Boeringer, GSK, Novo Nordisk, Jassen-Cilag, J&J and Zambon, outside the submitted work.
- Received October 21, 2022.
- Accepted February 26, 2023.
- Copyright ©The authors 2023
This version is distributed under the terms of the Creative Commons Attribution Non-Commercial Licence 4.0. For commercial reproduction rights and permissions contact permissions{at}ersnet.org