- Open Access
The value of high-resolution computed tomography (HRCT) to determine exercise ventilatory inefficiency and dynamic hyperinflation in adult patients with cystic fibrosis
Respiratory Researchvolume 20, Article number: 78 (2019)
In Cystic Fibrosis (CF), exercise ventilatory inefficiency and dynamic hyperinflation (DH) cause exercise limitation and induce poor exercise tolerance. High-resolution computed tomography (HRCT) of the lung can detect pulmonary abnormalities in CF patients. We aimed to identify the determinants of exercise ventilatory inefficiency and DH using HRCT-derived metrics.
Fifty-two adult CF patients were prospectively enrolled; all participants underwent cardio-pulmonary exercise test (CPET) and HRCT. Radiological impairment was evaluated by the Brody II scoring system. Slope and intercept of the minute ventilation/CO2 production (V’E/V’CO2) regression line and the ratio of inspiratory capacity/total lung capacity (IC/TLC) at rest and at peak of exercise were measured.
Four groups of patients were identified based on the combination of ventilatory efficiency (Vef) or inefficiency (Vin) and the presence/absence of DH. Compared to other groups, CF adults with Vin and DH had worse functional status and higher total (T), bronchiectasis (B) and air trapping (AT) scores at HRCT. Significant correlations were found between V’E/V’CO2 intercept and V’E/V’CO2 slope (ρ − 0.455, p = 0.001) and between V’E/V’CO2 intercept and Δ inspiratory capacity (IC) (ρ − 0.334, p = 0.015). Regression analysis identified AT score (cut-off 7.9, odds ratio-OR 3.50) as the only independent predictor of Vin and T (cut-off 53.6, OR 4.98), B (cut-off 16.1, OR 4.88), airways wall thickening (AWT) (cut-off 13, OR 3.41), and mucous plugging (MP) scores (cut-off 11.7, OR 4.18) as significant predictors of DH.
In adult CF cohort, values of HRCT metrics are determinants of Vin (AT) and DH (T, B, AWT, MP).
In Cystic Fibrosis (CF) patients, pulmonary infections lead to a progressive decline in lung function . Although spirometry and, particularly, the forced expiratory volume at 1st second (FEV1) is considered a valuable measure for lung function monitoring , variables derived from exercise testing, such as the cardio-pulmonary exercise test (CPET), are more sensitive than spirometry in detecting early structural pulmonary changes in CF patients .
In adolescents with CF, an exaggerated ventilatory response to exercise but adequate to V’CO2 (normal V’E/V’CO2 slope) was documented  and, interestingly, in adult CF patients the V’E/V’CO2 at peak of exercise was an independent determinant of exercise limitation, especially in those with a worse lung function . Moreover, in adults with CF, a high prevalence of exercise dynamic hyperinflation (DH) was found, secondary to the development of air trapping during exercise; DH was also shown to be associated to poor resting lung function, reduced exercise tolerance and increased exertional dyspnea .
High-resolution computed tomography (HRCT) of the lung is able to demonstrate a broad range of pulmonary abnormalities in CF patients with mild-to-moderate lung disease , identifying lower airway inflammation and early lung impairments . The Brody score  is a composite HRCT score including different radiological features of lung disease (Bronchiectasis-B, Airways Wall Thickening-AWT, Mucous plugging-MP, Parenchyma-P and Air Trapping-AT) and demonstrating a high inter-observer reproducibility . In CF patients, higher Brody scores indicate more severe disease ; they are also strongly associated with lung disease progression , predicting the rate of pulmonary exacerbations up to 10-year follow-up .
Interestingly, with regards to the relationship between HRCT-derived variables and the ventilatory response to exercise, the emphysema HRCT metrics were independent predictors of exercise-induced DH and V’E/V’CO2 slope in COPD patients . In CF adults, however, no studies assessed the value of the variables derived from HRCT to determine both the exercise ventilatory response to CO2 and DH. Therefore, we aimed to identify in a large cohort of adult patients with CF the determinants of exercise ventilatory response to CO2 and DH using HRCT-derived metrics, as expressed by the Brody II score.
We conducted this prospective study at the University Hospital of Parma (Italy) between June 2013 and July 2017. The main inclusion criteria were a confirmed diagnosis of CF as recommended by international guidelines , age > 18 years and a stable respiratory condition (i.e. patients free from exacerbations) for at least 4 weeks before enrolment Patients with lung transplant, pregnancy and concomitant malignancies or comorbidities (i.e. heart failure, previous rib fracture) which could preclude their safety were excluded.
The study was conducted according to the Declaration of Helsinki; all the procedures and their risks were explained to the patients, who gave their written informed consent for the participation to the study. The ethical committee of the University Hospital of Parma approved the protocol (approval number: 200084; 07 June 2013). All participants’ data were anonymously analyzed and reported. No extramural funding was used to support the study.
For all patients, anthropometric characteristics (age, gender, body mass index-BMI and fat-free mass index-FFMI), presence of diabetes, CFTR mutations (F508del/F508del, F508del/Other, or Other/Other), Pseudomonas aeruginosa colonization, long-term oxygen therapy and domiciliary non-invasive mechanical ventilation were recorded at the enrolment. BMI was calculated in kilograms per square meter. Body composition was also assessed by a foot-to-foot bioelectrical impedance analysis (BIA) method by a body composition analyzer (model SC-331S, Tanita, Tokyo, Japan). The fat-free mass (FFM) was standardized for height to calculate the FFM index (FFMI = FFM/height squared, in kilograms per square meter).
All lung function tests were performed according to the international recommendations . A flow-sensing spirometer and a body plethysmograph connected to a computer for data analysis (Vmax 22 and 6200, Sensor Medics, Yorba Linda, USA) were used for the measurements. FEV1 and forced vital capacity (FVC) were recorded; the FEV1/FVC ratio was considered as the index of airflow obstruction. Thoracic gas volume (TGV) was measured by body plethysmography with the patients painting against a closed shutter at a frequency slightly < 1 Hz and supporting their cheeks with their hands. Total lung capacity (TLC) was obtained as the sum of TGV and related inspiratory capacity (IC). To ensure reproducibility for each spirometry and lung volume variables at least three measurements were taken; the highest value was considered. The flow-sensor was calibrated before each test using a three-liter syringe. FEV1, FVC, and IC were expressed as a percentage of the predicted values .
Incremental CPET was performed according to the standardized international procedure . Patients were not pre-treated with β2 agonists before testing, but they could take their current therapies. After calibrating the oxygen and carbon dioxide analysers and flow mass sensor, patients were invited to sit on an electromagnetically braked cycle ergometer (Corival PB, Lobe Bv, Groningen, The Netherlands) and the saddle was adjusted properly to avoid the maximal extension of the knee. The exercise protocol involved an initial rest of 3 min, followed by unloaded cycling for another 3 min with an increment of 5 to 15 watts every minute, according to the patient’s anthropometry and degree of functional impairment, in order to achieve an exercise time between 8 and 12 min. Patients were invited to maintain a pedalling frequency of 60 rates per minute.
Breath-by-breath oxygen uptake (V’O2 in L/min), carbon dioxide production (V’CO2 in L/min), tidal volume (VT in L) and minute ventilation (V’E in L/min) were recorded during the test (CPX/D; Med Graphics, St Paul, MN, USA). Patients were continuously monitored with a 12-lead electrocardiogram (Welch Allyn CardioPerfect, Delft, the Netherlands) and a pulse oximeter (Pulse Oximeter 8600, Nonin Medical Inc., MPLS, Mn U.S.A.). Blood pressure was measured every two minutes. Stopping criteria consisted of symptoms, such as unsustainable dyspnea, leg fatigue or chest pain, a significant ST-segment depression at ECG, or a drop in systolic blood pressure or oxygen saturation ≤ 84% .
Workload and V’O2 at peak were recorded as the mean value of watts and V’O2 during the last 20 s of the test. V’O2 at peak was expressed as mL/kg/min. The ventilatory response during exercise was expressed as a linear regression function by plotting V’E against V’CO2 obtained every 10 s, excluding data above the ventilatory compensation point . Then, the slope and Y intercept values were obtained from the V’E/V’CO2 regression line. V’E/V’CO2 at rest, V’E/V’CO2 slope and V’E/V’CO2 intercept (L∙min− 1) were recorded. Values of V’E/V’CO2 slope ≥ 30 and < 30, derived from the median value in our cohort, were used to define patients with an exercise ventilatory inefficiency (Vin) and patients with ventilatory efficiency (Vef), respectively .
The end-tidal pressure of CO2 (PETCO2, in mm Hg) was measured as the mean of PETCO2 during the 3-min rest period and during the last 20 s of the test; moreover, the difference between PETCO2 at peak and PETCO2 at rest (Δ PETCO2) was recorded.
Changes in operational lung volumes were assessed every two minutes during exercise and at peak exercise, taking the IC measured at rest, as the baseline. After a full explanation of the procedure, satisfactory technique and reproducibility of IC maneuvers were established during an initial practice session at rest. Assuming that TLC remains constant during exercise , a change in IC (Δ IC, L) ≤ − 0.100 L at peak of exercise was considered to define patients with DH, as reported from the literature .
The cardiovascular response to exercise was expressed by the following parameters: oxygen pulse (O2 pulse), HR recovery and double product (DP) reserve. O2 pulse (in mL/beats/min) was calculated by dividing instantaneous VO2 by the HR and was recorded at rest and as the change between the peak of exercise and the resting phase (Δ O2 pulse) . HR recovery (in beats) was defined as the reduction of the HR at the peak exercise level compared to the HR after one minute of the exercise cessation . DP (mmHg*bpm) at rest and at maximal exercise was calculated by the product of systolic blood pressure and HR ; DP reserve was calculated as DP at maximal exercise minus DP at rest and was also expressed as a ratio to workload in watts.
Oxygen saturation (SpO2) was reported at rest and as the difference between values at peak of exercise and at rest (Δ SpO2). Peak dyspnea and fatigue were described as visual analogic scale (VAS) with 0–100 score and then divided by peak workload (watts) for analysis.
Patients were scanned with a 128-slice Somatom Definition Flash scanner (Siemens Medical Solutions, Forchheim, Germany). Technical parameters were as follows: volumetric acquisition, collimation 128 × 2 × 0.6 mm; 1 mm slice thickness, rotation time 0.28 msec; tube voltage 80–100 kVp; tube current modulation (CareDose) on. Images were visually scored using a window setting (− 1550 W, − 600 L). Each HRCT scan was evaluated independently by two radiologists (S.L. and P.P.), both with a 3-year imaging experience. Readers scored HRCTs using the Brody II scoring system . The total score (T) is the sum of the following sub-scores: Bronchiectasis-B, Airways Wall Thickening-AWT, Mucous Plugging-MP, Parenchyma-P (e.g. bullae, ground glass and consolidation) and Air Trapping-AT with decreased attenuation areas. Each of the five sub-scores is rated according to presence and severity, in each of the six lobes, lingula being considered as a separate lobe. All scores and sub-scores are expressed as a percentage of the maximum possible score, ranging from 0 (no pathological findings) to 100 (maximum severity).
Analyses were performed with IBM SPSS Statistics 25.0 (Armonk, New York, USA) and Shapiro-Wilk test was used to assess the normality of distribution. Data were reported as medians (1st quartile; 3rd quartile) for continuous variables with non-normal distribution or means ± standard deviation (SD) for those with normal distribution. Number of patients (%) was used for categorical variables. Comparison analysis were performed using the chi-square (X2) or the Fisher exact test for categorical variables and the t test or the nonparametric Mann-Whitney test for continuous variables. For multiple comparisons, the Bonferroni post-hoc analysis, the ANOVA, and the nonparametric Kruskal-Wallis test were used.
The interrater reliability was evaluated by a two-way random-effects model with mean of k raters and definition of consistency . The intraclass correlation coefficients (ICC) were interpreted by the following scale: ICC values < 0.50 are indicative of “poor reliability”, values between 0.50 and 0.75 indicate “moderate reliability”, values between 0.75 and 0.90 indicate “good reliability”, and values greater than 0.90 indicate “excellent reliability” .
For correlation analysis, the Pearson or Spearman correlation coefficients were used for linear or normally distributed variables and for not linear or not normally distributed variables, respectively.
Receiver operating characteristic (ROC) curves were generated to calculate the area under the curve (AUC) with 95% confidence interval (CI) and to find the best cut-off with the related sensibility and specificity. Univariate and multivariate regression logistic models (method: Stepwise) were performed to predict the probability to have Vin and DH (set as dependent variables). Odds ratios (OR) and 95% CI were calculated. A p value at two tails < 0.05 was considered significant.
Our study sample consisted of 52 adult patients with CF, normal-weight, with a wide range of lung function from normal to severe airflow obstruction and a normal-to-mildly reduced maximal exercise capacity. Variables related to spirometry and CPET demonstrated that, in comparison to patients with Vef and without DH and patients with Vin and without DH, patients with Vin and with DH had statistically significant lower values of FEV1% pred., FVC % pred., FEV1/FVC %, and lower values in V’O2 at peak ml/kg/min and workload in watts at peak without reaching statistical significance. Table 1 summarises the results about anthropometric and lung function characteristics and Table 2 CPET-related variables.
The interrater reliability between the two radiologists was “excellent” for T (ICC 0.95; 95% CI 0.91 to 0.97), B (ICC 0.93; 95% CI 0.88 to 0.96), MP (ICC 0.93; 95% CI 0.88 to 0.96), and P (ICC 0.91; 95% CI 0.84 to 0.95), while was “good” for AWT (ICC 0.86; 95% CI 0.75 to 0.92) and AT (ICC 0.81; 95% CI 0.66 to 0.88). Then, the average value of each Brody II score was considered for all analyses. Examples of HRCT metrics according to the Brody II score are provided in Fig. 1, while in Fig. 2 are illustrated for all patients the HRCT variables related to the Brody II scores. Among Vin patients (Fig. 3, bottom left), only AT differed significantly compared to those with Vef. Patients with DH, in comparison to those without DH (Fig. 3, bottom right), had higher values of T, B, AWT, and MP, respectively. Patients with Vin and with DH had higher values of T, B and AT compared to other groups (Fig. 3, top).
Significant correlations were found between V’E/V’CO2 intercept and V’E/V’CO2 slope and between V’E/V’CO2 intercept and Δ IC. No correlation, however, was found between V’E/V’CO2 slope and Δ IC (Fig. 4). Multiple comparisons between patient subgroups (Fig. 5) showed that V’E/V’CO2 intercept was significantly higher among patients with Vef and with DH compared to those with Vef and without DH and those with Vin and without DH. Moreover, V’E/V’CO2 intercept was significantly lower in patients with Vin and without DH than those with Vef and without DH and those with Vin and with DH.
Regression analysis (Table 3) showed that AT was the only independent predictor of Vin, whilst T, B, AWT and MP were significant predictors of the development of DH.
Our study in adult patients with CF allows us to highlight two main findings: 1) there is no correlation between Vin and DH and variables derived from spirometry and CPET seem to be complementary and specific to one condition rather another and 2) Brody II scores at HRCT are able to discriminate and predict Vin (AT) and DH (T, B, AWT and MP).
Ventilatory response to exercise in CF patients
There are very few data about the ventilatory response to exercise in adult CF patients. In a retrospective study on 102 CF patients with mild-to-moderate lung disease, the authors reported that VE/VCO2 at peak was not related to airflow obstruction but to the VO2 at peak . High values of VE/VCO2 slope identify greater dead space (VD) due to more extensive emphysema and/or high pulmonary vascular pressures, poorer cardiac performance, higher sympathetic drive, worse exertional hypoxemia and greater ergoreceptor stimulation . Indeed, in a sample of CF children undergoing sub-maximal constant work rate exercise testing, the arterial partial pressure of carbon dioxide (PaCO2) levels were significantly related to alveolar ventilation rather than VE; this denotes the leading role of physiologic VD in the quantitative determination of ventilatory response to exercise .
In our cohort, data about the relationship between DH and lung function impairment at rest are in line with Stevens et al. . Surprisingly, we did not observe differences in term of exercise tolerance ; however, patients with DH and Vin had a tendency to show a poorer exercise tolerance (Table 2). It is of note that in COPD patients, VE/VCO2 measurements may explain exercise intolerance, particularly in milder and in symptomatic patients with preserved FEV1 , since an increased VD may occur .
The value of HRCT metrics in Vin and DH
HRCT scores and CPET-variables have been often used separately to predict future outcomes in CF patients, such as early structural alterations  or disease progression [10, 11]. However, in adult CF patients no studies looked for determinants of two exercise-induced functional scenario (presence of Vin and DH) using HRCT scores.
Air trapping at HCRT is an important marker of distal airways disease in children with CF  and in COPD patients . Using a modified Bhalla scoring system, a retrospective study  demonstrated that HRCT abnormalities were more sensible than lung function (evaluated only by FEV1) to detect disease decline of adult CF patients; notably, patients with air trapping and moderate-to-severely impaired lung function had more rapid rate of deterioration . Moreover, in children with CF, air trapping demonstrated a significant relationship with non-invasive measures of ventilation inhomogeneity, such as the lung clearance index (LCI) . These aspects are in line with our regression analysis, in which AT score independently increases the probability to have Vin which, at least for COPD patients, may define a series of worse disease-related features . Similarly, HRCT metrics on emphysema measured by low attenuation areas (LAA%) - that conceptually may be seen as a surrogate of AT - were shown to significantly explain VE/VCO2 slope in COPD patients in a multivariate adjusted linear regression model .
MP, secondary to the basic ion-transport defect of the disease, represents a recognized HRCT hallmark of disease  and, together with perfusion abnormalities seems to play a relevant role in exacerbation . Based on our results, it may also have a potentially reversible role on DH.
The V’E/V’CO2 intercept in CF patients
The relationship between V’E and V’CO2 - when V’CO2 is equal to zero - represents the V’E/V’CO2 intercept . In normal subjects, V’E/V’CO2 intercept is a small positive value (< 3 L∙min− 1) . In COPD patients, V’E/V’CO2 intercept is related to greater DH  and it is associated with exertional dyspnea .
In our CF patients, V’E/V’CO2 intercept values lead us to some considerations. Among the whole study group, the median value was high (3.28) and, in particular, it was higher in patients with Vef (V’E/V’CO2 slope < 30) and with DH. As shown for COPD patients , a high V’E/V’CO2 intercept is a necessary compromise for the progressive mechanical respiratory constraints, increased VD and reduced V’E/V’CO2 slope. Notably, it is proved that V’E/V’CO2 intercept reflects some aspects about the VD . The increased VD and ventilation inhomogeneity  in our CF patients may explain the increased value of V’E/V’CO2 intercept in DH patients. Then, in this context it is plausible that V’E/V’CO2 intercept correlates with both V’E/V’CO2 slope and Δ IC. Because each condition (Vin and DH) influences, in different ways, the V’E/V’CO2 intercept (Table 2), the level of this parameter does not increase if the two conditions coexist (Fig. 5).
Strength and limitations of the study
Although our study has some important points of strength (originality because for the first time the exercise-induced ventilatory response was considered and the wide cohort of adult CF, in which all participant performed both HRCT and CPET) we need to report some limitations. First, our study concerns data collected in a single center in Italy, in adult CF with a wide range of lung function from normal to severe airflow obstruction and a relatively well-conserved maximal exercise capacity. The considerations derived from our analyses cannot be translated into patients with different functional and physical characteristics. Secondly, data about ATs have been collected only during the inspiration phase. Although in CF a pilot study demonstrated a very good agreement between HRCT scores from end-expiratory and end-inspiratory scans , we cannot be sure that the execution of the end-expiratory phase scans could have changed our findings.
In summary, our prospective study performed in a cohort of adult patients with CF showed the values of HRCT metrics as determinants of Vin (AT) and DH (T, B, AWT and MP). Moreover, we highlighted in CF patients the functional and clinical potential role of Vin and DH. We believe that in the future these aspects should be considered when exercise tolerance is assessed, such as for COPD patients. The comprehensive scenario for CF patients could be more interesting.
adjacent pulmonary artery outer diameter
area under the curve
bronchial lumen diameter
body mass index
bronchial wall thickening
cardio-pulmonary exercise test
- FEV1 :
forced expiratory volume at 1st second
forced vital capacity
high-resolution computed tomography
intraclass correlation coefficients
- O2 pulse:
- PETCO2 :
end-tidal pressure of CO2
receiver operating characteristic
- SpO2 :
total lung capacity
visual analogic scale
- VCO2 :
carbon dioxide production
- VE :
- Vef :
- Vin :
- VT :
Gibson RL, Burns JL, Ramsey BW. Pathophysiology and management of pulmonary infections in cystic fibrosis. Am J Respir Crit Care Med. 2003;168(8):918–51.
Tiddens HA. Detecting early structural lung damage in cystic fibrosis. Pediatr Pulmonol. 2002;34(3):228–31.
Hatziagorou E, Kampouras A, Avramidou V, et al. Exercise responses are related to structural lung damage in CF pulmonary disease. Pediatr Pulmonol. 2016;51(9):914–20.
Bongers BC, Werkman MS, Takken T, Hulzebos EH. Ventilatory response to exercise in adolescents with cystic fibrosis and mild-to-moderate airway obstruction. Springerplus. 2014;3:696.
Pastré J, Prévotat A, Tardif C, Langlois C, Duhamel A, Wallaert B. Determinants of exercise capacity in cystic fibrosis patients with mild-to-moderate lung disease. BMC Pulm Med. 2014;14:74.
Stevens D, Stephenson A, Faughnan ME, Leek E, Tullis E. Prognostic relevance of dynamic hyperinflation during cardiopulmonary exercise testing in adult patients with cystic fibrosis. J Cyst Fibros. 2013;12(6):655–61.
Brody AS, Klein JS, Molina PL, Quan J, Bean JA, Wilmott RW. High-resolution computed tomography in young patients with cystic fibrosis: distribution of abnormalities and correlation with pulmonary function tests. J Pediatr. 2004;145(1):32–8.
Davis SD, Fordham LA, Brody AS, et al. Computed tomography reflects lower airway inflammation and tracks changes in early cystic fibrosis. Am J Respir Crit Care Med. 2007;175(9):943–50.
Weber K, Paolini M, Schmitz M, et al. Cystic fibrosis in adults: short-term and long-term reproducibility of the Brody score for lung morphology in low-dose MDCT scans. Rofo. 2014;186(1):54–60.
Sanders DB, Li Z, Brody AS, Farrell PM. Chest computed tomography scores of severity are associated with future lung disease progression in children with cystic fibrosis. Am J Respir Crit Care Med. 2011;184(7):816–21.
Sanders DB, Li Z, Brody AS. Chest computed tomography predicts the frequency of pulmonary exacerbations in children with cystic fibrosis. Ann Am Thorac Soc. 2015;12(1):64–9.
Crisafulli E, Alfieri V, Silva M, et al. Relationships between emphysema and airways metrics at high-resolution computed tomography (HRCT) and ventilatory response to exercise in mild to moderate COPD patients. Respir Med. 2016;117:207–14. https://doi.org/10.1016/j.rmed.2016.06.016.
Farrell PM, White TB, Ren CL, et al. Diagnosis of Cystic Fibrosis: Consensus Guidelines from the Cystic Fibrosis Foundation. J Pediatr. 2017;181S:S4–S15.e1. https://doi.org/10.1016/j.jpeds.2016.09.064.
Miller MR, Hankinson J, Brusasco V, et al. ATS/ERS task force. Standardisation of spirometry. Eur Respir J. 2005;26(2):319–38.
Wanger J, Clausen JL, Coates A, et al. Standardisation of the measurement of lung volumes. Eur Respir J. 2005;26(3):511–22.
American Thoracic Society; American College of Chest Physicians. ATS/ACCP statement on cardiopulmonary exercise testing. Am J Respir Crit Care Med. 2003;167(2):211–77.
Puente-Maestu L, Palange P, Casaburi R, et al. Use of exercise testing in the evaluation of interventional efficacy: an official ERS statement. Eur Respir J. 2016;47(2):429–60. https://doi.org/10.1183/13993003.00745-2015.
Neder JA, Berton DC, Arbex FF, et al. Physiological and clinical relevance of exercise ventilatory efficiency in COPD. Eur Respir J. 2017;49(3). https://doi.org/10.1183/13993003.02036-2016.
Stubbing DG, Pengelly LD, Morse JL, Jones NL. Pulmonary mechanics during exercise in subjects with chronic airflow obstruction. J Appl Physiol. 1980;49:511–5.
Regnis JA, Alison JA, Henke KG, Donnelly PM, Bye PT. Changes in end-expiratory lung volume during exercise in cystic fibrosis relate to severity of lung disease. Am Rev Respir Dis. 1991;144:507–12.
Gobel FL, Norstrom LA, Nelson RR, Jorgensen CR, Wang Y. The rate-pressure product as an index of myocardial oxygen consumption during exercise in patients with angina pectoris. Circulation. 1978;57(3):549–56.
McGraw KO, Wong SP. Forming inferences about some intraclass correlation coefficients. Psychol Methods. 1996;1:30–46.
Koo TK, Li MY. A guideline of selecting and reporting Intraclass correlation coefficients for reliability research. J Chiropr Med. 2016;15(2):155–63. https://doi.org/10.1016/j.jcm.2016.02.012.
Pianosi P, Wolstein R. Carbon dioxide chemosensitivity and exercise ventilation in healthy children and in children with cystic fibrosis. Pediatr Res. 1996;40(3):508–13.
Elbehairy AF, Ciavaglia CE, Webb KA, et al. Canadian respiratory research network. Pulmonary gas exchange abnormalities in mild chronic obstructive pulmonary disease. Implications for dyspnea and exercise intolerance. Am J Respir Crit Care Med. 2015;191(12):1384–94.
Bonnel AS, Song SM, Kesavarju K, et al. Quantitative air-trapping analysis in children with mild cystic fibrosis lung disease. Pediatr Pulmonol. 2004;38(5):396–405.
Lee SM, Seo JB, Lee SM, Kim N, Oh SY, Oh YM. Optimal threshold of subtraction method for quantification of air-trapping on coregistered CT in COPD patients. Eur Radiol. 2016;26(7):2184–92.
Judge EP, Dodd JD, Masterson JB, Gallagher CG. Pulmonary abnormalities on high-resolution CT demonstrate more rapid decline than FEV1 in adults with cystic fibrosis. Chest. 2006;130(5):1424–32.
Owens CM, Aurora P, Stanojevic S, et al. London cystic fibrosis collaboration. Lung clearance index and HRCT are complementary markers of lung abnormalities in young children with CF. Thorax. 2011;66(6):481–8.
Wielpütz MO, Eichinger M, Biederer J, et al. Imaging of cystic fibrosis lung disease and clinical interpretation. Rofo. 2016 Sep;188(9):834–45.
Wielpütz MO, Puderbach M, Kopp-Schneider A, et al. Magnetic resonance imaging detects changes in structure and perfusion, and response to therapy in early cystic fibrosis lung disease. Am J Respir Crit Care Med. 2014;189:956–65.
Sun XG, Hansen JE, Garatachea N, Storer TW, Wasserman K. Ventilatory efficiency during exercise in healthy subjects. Am J Respir Crit Care Med. 2002;166(11):1443–8.
Gargiulo P, Apostolo A, Perrone-Filardi P, Sciomer S, Palange P, Agostoni P. A non invasive estimate of dead space ventilation from exercise measurements. PLoS One. 2014;9(1):e87395.
Loeve M, Lequin MH, de Bruijne M, et al. Cystic fibrosis: are volumetric ultra-low-dose expiratory CT scans sufficient for monitoring related lung disease? Radiology. 2009;253(1):223–9. https://doi.org/10.1148/radiol.2532090306.
This work was not supported by any financial source.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The ethical committee of the University Hospital of Parma approved the protocol (approval number: 200084; 07 June 2013).
Consent for publication
The authors declare that they have no competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.