Abstract
Introduction
Cystic fibrosis (CF) is an autosomal recessive disease caused by mutations in the
Functional respiratory imaging (FRI) is a quantifiable method to measure biological responses to therapeutic interventions using high-resolution computed tomography (HRCT) and computational fluid dynamics (CFD). A 3D reconstruction of the respiratory system combined with flow simulations can be used to determine parameters such as patient-specific airflow distribution, airway resistance and volume, lobar volume, and air trapping. This functional and structural evaluation of the lungs allows for an in-depth description of respiratory health and post-treatment effects. 9 Previous studies have demonstrated the validity of the lobar ventilation models and CFD simulations by comparing the results to 3He magnetic resonance imaging (MRI) ventilation measurements and combined single photon emission computed tomography and CT (SPECT/CT) in patients with asthma.10,11 Besides a detailed report of characteristics on a lobar level, this method has shown to be a more sensitive technique than conventional lung function tests in various obstructive airway diseases. In chronic obstructive pulmonary disease (COPD) and asthma, several studies assessing the effects of bronchodilators and anti-inflammatory compounds showed significant improvements in FRI parameters, while these changes were less pronounced in the results of spirometry and body plethysmography.12–14 In the CF population, FRI has demonstrated its value by establishing patient-specific airway models to predict the deposition of inhaled antibiotics. 15 Given these previous findings, we hypothesize that FRI parameters related to hyperinflation, airway volume and airway wall thickness could provide valuable insights into the effects of LUM/IVA as imaging biomarkers on lung structure and functionality. This study aims to offer additional and novel information about the mode of action of the CFTR-modulating drug by assessing lung structure and function across multiple FRI parameters.
Methods
Study design
An open-label, single-arm study was conducted to investigate the effects of LUM/IVA over an intervention period of 12 weeks. The investigational product, Orkambi®, was administered as two tablets (lumacaftor 200 mg/ivacaftor 125 mg) twice a day in compliance with the prescribing information and the European Medicines Agency (EMA) approval. Treatment was started after the completion of all baseline measurements. Participants visited the hospital every 4 weeks during the 12-week period at which following study assessments were carried out: physical examination and spirometry. Other assessments were only performed at baseline and at 12 weeks: CT chest imaging, nitrogen multiple breath washout (N2MBW) test, Six-Minute Walk Test (6MWT), sweat chloride test and the Cystic Fibrosis Questionnaire–Revised (CFQ-R). Adverse events and adverse drug reactions were monitored throughout the study. The study was approved by the Ethics committee of the Antwerp University Hospital (EudraCT 2018-001573-24).
Participants
Participants aged 12 years or older were eligible for inclusion if they had a confirmed diagnosis of CF, homozygous for the
Study assessments
Functional respiratory imaging
HRCT scans were taken with a GE VCT LightSpeed 64-slice scanner (GE Healthcare, Chicago, Illinois, USA) at two breathing levels, total lung capacity (TLC) and functional residual capacity (FRC) monitored with a pneumotachograph. A description of specific CT settings can be found in the Supplementary Material. Images were imported into a medical image processing software package, Mimics (Materialize, Leuven, Belgium), for segmentation and 3D reconstruction of the airways and lung lobes. All segmentation and postprocessing steps were automated, except for the identification of the fissures to separate the lung lobes, which were manually adjusted. The computed models were used to determine structural parameters, such as lobar volume, airway volume, and air trapping. After segmentation and postprocessing, the models were used for CFD simulations by solving Reynolds-averaged Navier-Stokes equations, to calculate regional airway resistance and to predict the deposition pattern of inhaled drugs. Pulmonary vasculature was segmented using algorithms based on shape analysis to recognize tubular structures. The cross-sectional area of each identified blood vessel was determined to compute pulmonary blood distribution. Although Orkambi® is not delivered by aerosol, simulations were performed to evaluate whether the treatment would affect the lung deposition of one of the most frequently inhaled drugs by CF patients, dornase alfa. All subjects were coupled with the Pari eFlow Nebulizer (Starnberg, Germany) and simulated using a tidal breathing profile. The particle characteristics of dornase alfa were taken from literature. 16 Details of the analysis can be found in the Supplementary Material and previously published work.11,17 The following FRI parameters were evaluated in this study: air trapping, lobar volume, airway volume, airway wall volume, airway resistance, pulmonary blood distribution, ventilation/perfusion matching, internal airflow distribution and aerosol deposition.
Secondary outcomes
Structural abnormalities on inspiratory and expiratory chest CTs were visually scored using the CF-CT scoring system, a validated upgraded version of the Brody II score.18,19 The scans were scored by two independent observers blinded to subject ID and study visit, of whom one experienced chest radiologist and one certified researcher trained in CF-CT scoring. The presence and extent of the following abnormalities were evaluated: bronchiectasis, mucus plugging, parenchymal abnormalities and air trapping. All scores were expressed as a percentage of the maximum score. Data recorded during the monthly physical examination included body weight and length, peripheral oxygen saturation (SpO2), pulse rate and blood pressure. Pulmonary function was measured by spirometry and N2MBW. Exercise tolerance was evaluated by the Six-Minute Walk Distance (6MWD). During the exercise test, pulse oximetry was used to monitor SpO2. Also, the Borg score for dyspnoea and fatigue were questioned before and immediately after the 6MWT. Health-related quality of life (QoL) was evaluated using the CFQ-R. Sweat chloride concentration was measured by pilocarpine iontophoresis. A more detailed description of the study assessments can be found in the Supplementary Material.
Statistical analyses
An initial sample size of 14 subjects was set to reach a power of 80% with a significance level of 0.05. Since no pilot data were available, the calculation was based on the change in airway volume (Cohen’s effect size of 0.82) reported in a previous FRI study including patients with COPD. 12
The analysis of the FRI parameters was based on a linear mixed-effects model, with visit, lobe and their interaction as fixed effects and lung lobe within each subject as a random effect. The heterogeneity across lobes (within subject) was modelled using an unstructured variance-covariance matrix, with independence assumed between subjects. For spirometry, anthropometrics and vital signs a similar approach was taken, excluding the lobar component. Changes in CF-CT scores, exercise tolerance, CFQ-R domain scores, lung clearance index and sweat chloride concentration were explored with the paired samples
Results
Participants
Twelve participants were enrolled in the study, of whom one discontinued treatment after 10 weeks due to respiratory complaints, including cough and increased sputum. Therefore, the latter subject was excluded in further analyses determining the effects of LUM/IVA. Besides the discontinuation of the study at the patient’s own initiative, no serious adverse events were reported throughout the study. Baseline characteristics are presented in Table 1. The intended sample size of 14 subjects was not reached, since no other participants could be motivated to enter the trial. The main reason was the relatively short intervention period in combination with the fact that patients were not able to continue taking the drug after the trial, as Orkambi® was not reimbursed in Belgium at that time.
Baseline characteristics (
BMI, body mass index.
Data are presented as mean ± standard deviation or median [range].
FRI parameters
The FRI analysis showed a decrease in air trapping (16.12 ± 19.13% to 13.14 ± 17.55%,

Mean change and individual changes of total (a) air trapping (%), (b) lobar volume at FRC (L), and (c) airway wall volume (mL/L).

Treatment response in subject 5. The average change per lung lobe from baseline in (a) air trapping, (b) lobar volume at FRC, and (c) airway wall volume.
Results FRI (
AT, air trapping; BV5, blood volume of vessels smaller than 5 mm2; BV5%, BV5 in relation to total blood volume; BV5_10, blood volume of vessels between 5 and 10 mm2; BV5_10%, BV5_10 in relation to total blood volume; BV10, blood volume of vessels larger than 10 mm2; BV10%, BV10 in relation to total blood volume; C/P ratio, central-to-peripheral deposition; FRI, functional respiratory imaging; IAD, internal airflow distribution; iVlobe, lobar volume; PLD, peripheral lung deposition; siRaw, specific airway resistance; siVaw, specific airway volume; siVaww, specific airway wall volume; TLD, total lung deposition; V/Q, ventilation/perfusion.
Data are presented as mean ± standard deviation.
Statistically significant (
Secondary outcomes
The interobserver variability analysis of the CF-CT scoring method showed a moderate to excellent agreement between the two observers (Table 3). None of the CF-CT components nor the total score changed significantly after treatment.
Interobserver variability of the CF-CT scoring method presented as ICC [95% CI].
CF, cystic fibrosis; CI, confidence interval; CT, computed tomography; ICC, intraclass correlation coefficient.
ppFEV1 increased by 2.0 ± 4.0% (
The median 6MWD increased from 650 m [512; 790] to 725 m [525; 839] (
The CFQ-R questionnaire showed a trend towards an improvement of the physical domain score as the majority of the patients scored higher after treatment (88% (17; 100) to 88% (58; 100),
Body mass index (BMI) increased steadily during 12 weeks of treatment (22.2 ±–3.7 kg/m2 to 22.9 ± 4.0 kg/m2,
Finally, sweat chloride concentration tended to decrease after 12 weeks (91.6 ± 14.5 mmol/L to 81.6 ± 14.7 mmol/L,
Results secondary outcomes (
BMI, body mass index; CF, cystic fibrosis; CFQ-R, Cystic Fibrosis Questionnaire–Revised; CT, computed tomography; 6MWD, Six-Minute Walk Distance.
Data are presented as mean ± standard deviation or median [range].
Statistically significant (
Correlations FRI and secondary outcomes
Significant correlations were found at baseline between air trapping and ppFEV1 (

Baseline correlations between air trapping measured by FRI and ppFEV1, 6MWD, CFQ-R physical and CF-CT total score.
Between secondary outcomes, ppFEV1 was correlated with CF-CT total score (
Discussion
This study aimed to gain more insight into the short-term effects of LUM/IVA by evaluating the patient’s lung structure and function across multiple FRI parameters. The results showed a decrease in lung hyperinflation indicated by a significant decrease in air trapping and in lobar volume after normal expiration. Furthermore, airway wall volume decreased and the proportion of pulmonary blood volume shifted from larger to smaller blood vessels. The larger proportion of blood in smaller blood vessels would indicate an average decrease in pulmonary vascular resistance downstream that resulted from an improved regional ventilation. The absolute volume of BV5, however, did not significantly change. We hypothesize that both the redistribution of pulmonary blood and the decrease in airway wall volume could be related to a reduction in mucus impaction, due to the similar attenuation properties of mucus, airway walls and blood vessels in the CT images causing bias. Consequently, the calculation of the total airway wall thickness includes mucus adhering to the airway wall. Second, mucus plugs obstructing small airways could have been recognized by the applied algorithms as larger blood vessels due to the corresponding shape and attenuation values. This could explain the significant correlation between the CF-CT subscore for mucus plugging and the absolute values of BV5_10 and BV10 at baseline. Although these assumptions indicate improved mucociliary clearance, the CF-CT score for mucus plugging did not show any differences after treatment. This could be related to the insensitivity of the CF-CT score to detect subtle changes, as the smallest increment of the score represents an increase in the extent of one-third of a lobe.
21
The same reasoning can be applied to explain the lack of change in other CF-CT scores. More particularly, differences in air trapping and airway wall volume determined by FRI were not reflected in CF-CT scores for air trapping and bronchial wall thickening, respectively. In contrast to our study, a recent retrospective study evaluating structural changes using the modified Brody score after the initiation of LUM/IVA therapy showed a significant improvement in the CT total score, which was mostly related to a decrease in mucus plugging.
22
The higher number of participants and longer treatment period compared with our study could explain these significant differences. In addition, Sheikh
Another important finding of this study is that the changes in several FRI parameters were also not reflected in any of the spirometric parameters in this small group of patients. The limitations of spirometry to measure therapeutic responses have been widely recognized: they are considered too insensitive to detect treatment effects and the heterogeneity of the pathology cannot be captured.
25
Nevertheless, FEV1 is still considered one of the most important primary endpoints by the EMA for research in CF, since FEV1 is an established marker for disease progression and a strong predictor of survival.26,27 Over the past decade, alternative endpoints for clinical research in CF have been proposed, of which the LCI, CT scores and magnetic resonance imaging (MRI) are regarded to be the most promising.28,29 LCI, CT scores as well as MRI outcomes have shown to be more sensitive to detect early lung disease, which is an important advantage since the progression of CF lung disease is slowing down and the FEV1 remains within the normal range at older age.25,30,31 Two previous observational studies showed that LUM/IVA therapy was associated with a decrease in LCI, while no significant changes were observed in ppFEV1.32,33 Unfortunately, a recent real-world study by Reix
Although FRI as described in this study is a new tool in the field of CF research, previous studies in other obstructive lung diseases, such as COPD and asthma, have shown that FRI is able to overcome several of the shortcomings or challenges associated with other methods as described earlier.12,39 Although FRI is based on CT imaging, this imaging modality currently remains the gold standard to evaluate structural abnormalities of the respiratory system, due to its higher spatial resolution compared with MRI. FRI biomarkers obtain regional information of both structural and functional characteristics using objective image processing techniques. Therefore, FRI has the potential to provide valuable contributions to future interventional studies targeting CF lung disease. Similar to studies in COPD and asthma patients, our results confirm that FRI is more sensitive than conventional outcomes, such as spirometry and CF-CT scores, to pick up local changes in the respiratory system. Consequently, a smaller sample size is required and the intervention period of a clinical trial can be reduced. This is an important benefit since disease progression in CF is slowing down as mentioned earlier, which would imply larger clinical studies of a longer duration to demonstrate significant changes in conventional outcomes. 40
In this study, important correlations were found between FRI and other respiratory outcomes, including physical functioning, spirometry, and structural abnormalities on chest CT. This indicates that improvements in specific FRI parameters could be related to an improved clinical outcome in the long term. Correlations between CT findings and lung function in CF have been reported in previous research, showing a higher correlation between LCI and CT than FEV1 and CT.41,42 Therefore, it would have been interesting to compare the LCI with FRI parameters, but unfortunately the number of successful MBW tests was insufficient to perform such analysis.
Although the result of ppFEV1 in this study was not significant, the mean change from baseline was comparable to the findings of the pivotal clinical trials. 6 Changes in most secondary outcomes were also similar to those reported in previous research, including nutritional status, 7 exercise tolerance, 43 QoL6,44 and arterial blood pressure. 7
The study is subject to several limitations. First, the obtained results could not be compared to a control group. A high variability was found between subjects across several FRI parameters as well as secondary outcomes, which can be related to the varying levels of disease severity and wide age range representing a heterogeneous group of patients. No previous FRI studies have been performed in the CF population to determine which changes can be expected over time and what treatment effects can be considered clinically important.
Second, only 12 participants were enrolled instead of the calculated sample size of 14. Although this is relatively small sample size, it should be stressed that all FRI parameters were evaluated at a lobar level, which implies five data points per subject. This partially explains why significant results could be detected by FRI in contrast to secondary outcomes. In retrospect, it would have been more appropriate for a study including CF patients to perform the sample size calculation using an FRI parameter other than airway volume, such as air trapping. Due to the presence of structural airway damage, in particular bronchiectasis, LUM/IVA is not expected to affect airway volume considerably during a short intervention period. The choice of airway volume to estimate the effect size was based on a previous study including COPD patients, since no other FRI studies with CF patients have been performed to date.
Over the past few years, new CFTR modulators have been proposed for the treatment of CF patients homozygous for
Conclusion
This study is the first to perform FRI to evaluate the effects of CFTR-modulating drugs on lung functionality. The analysis showed that LUM/IVA decreased lung hyperinflation in combination with a potential decrease in mucus impaction, which can be related to an improved mucociliary transport. These FRI parameters, reflecting regional and distal lung structures, appear to be more sensitive to changes caused by LUM/IVA than conventional outcomes, such as spirometry and subjective CT scores. The objective and quantitative FRI methodology could play an important role in acquiring profound knowledge about the local effects in the respiratory system of new CFTR modulators in future research.
Supplemental Material
sj-docx-1-tar-10.1177_17534666211046774 – Supplemental material for The short-term effects of ORKAMBI (lumacaftor/ivacaftor) on regional and distal lung structures using functional respiratory imaging
Supplemental material, sj-docx-1-tar-10.1177_17534666211046774 for The short-term effects of ORKAMBI (lumacaftor/ivacaftor) on regional and distal lung structures using functional respiratory imaging by Eline Lauwers, Dennis Belmans, Benjamin Mignot, Kris Ides, Kim Van Hoorenbeeck, Annemiek Snoeckx, Cedric Van Holsbeke, Vicky Nowé, Eva Van Braeckel, Wilfried De Backer, Jan De Backer and Stijn Verhulst in Therapeutic Advances in Respiratory Disease
Footnotes
Conflict of interest statement
Funding
Supplemental material
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
