Influence of exposure parameters and iterative reconstruction on automatic airway segmentation and analysis on MDCT-An ex vivo phantom study

GOLD stage-specific phenotyping of emphysema and airway disease using quantitative computed tomography

Background

In chronic obstructive pulmonary disease (COPD) abnormal lung function is related to emphysema and airway obstruction, but their relative contribution in each GOLD-stage is not fully understood. In this study, we used quantitative computed tomography (QCT) parameters for phenotyping of emphysema and airway abnormalities, and to investigate the relative contribution of QCT emphysema and airway parameters to airflow limitation specifically in each GOLD stage.

Methods

Non-contrast computed tomography (CT) of 492 patients with COPD former GOLD 0 COPD and COPD stages GOLD 1-4 were evaluated using fully automated software for quantitative CT. Total lung volume (TLV), emphysema index (EI), mean lung density (MLD), and airway wall thickness (WT), total diameter (TD), lumen area (LA), and wall percentage (WP) were calculated for the entire lung, as well as for all lung lobes separately. Results from the 3rd-8th airway generation were aggregated (WT3-8, TD3-8, LA3-8, WP3-8). All subjects underwent whole-body plethysmography (FEV1%pred, VC, RV, TLC).

Results

EI was higher with increasing GOLD stages with 1.0 ± 1.8% in GOLD 0, 4.5 ± 9.9% in GOLD 1, 19.4 ± 15.8% in GOLD 2, 32.7 ± 13.4% in GOLD 3 and 41.4 ± 10.0% in GOLD 4 subjects (p < 0.001). WP3-8 showed no essential differences between GOLD 0 and GOLD 1, tended to be higher in GOLD 2 with 52.4 ± 7.2%, and was lower in GOLD 4 with 50.6 ± 5.9% (p = 0.010 - p = 0.960). In the upper lobes WP3-8 showed no significant differences between the GOLD stages (p = 0.824), while in the lower lobes the lowest WP3-8 was found in GOLD 0/1 with 49.9 ± 6.5%, while higher values were detected in GOLD 2 with 51.9 ± 6.4% and in GOLD 3/4 with 51.0 ± 6.0% (p < 0.05). In a multilinear regression analysis, the dependent variable FEV1%pred can be predicted by a combination of both the independent variables EI (p < 0.001) and WP3-8 (p < 0.001).

Conclusion

QCT parameters showed a significant increase of emphysema from GOLD 0-4 COPD. Airway changes showed a different spatial pattern with higher values of relative wall thickness in the lower lobes until GOLD 2 and subsequent lower values in GOLD3/4, whereas there were no significant differences in the upper lobes. Both, EI and WP5-8 are independently correlated with lung function decline.

Standardized airway wall thickness Pi10 from routine CT scans of COPD patients as imaging biomarker for disease severity, lung function decline, and mortality

BACKGROUND

Chest computed tomography (CT) is increasingly used for phenotyping and monitoring of patients with COPD. The aim of this work was to evaluate the association of Pi10 as a measure of standardized airway wall thickness on CT with exacerbations, mortality, and response to triple therapy.

METHODS

Patients of GOLD grades 1-4 of the COSYCONET cohort with prospective CT scans were included. Pi10 was automatically computed and analyzed for its relationship to COPD severity, comorbidities, lung function, respiratory therapy, and mortality over a 6-year period, using univariate and multivariate comparisons.

RESULTS

We included n = 433 patients (61%male). Pi10 was dependent on both GOLD grades 1-4 (p = 0.009) and GOLD groups A-D (p = 0.008); it was particularly elevated in group D, and ROC analysis yielded a cut-off of 0.26 cm. Higher Pi10 was associated to lower FEV₁ % predicted and higher RV/TLC, moreover the annual changes of lung function parameters (p < 0.05), as well as to an airway-dominated phenotype and a history of myocardial infarction (p = 0.001). These associations were confirmed in multivariate analyses. Pi10 was lower in patients receiving triple therapy, in particular in patients of GOLD groups C and D. Pi10 was also a significant predictor for mortality (p = 0.006), even after including multiple other predictors.

CONCLUSION

In summary, Pi10 was found to be predictive for the course of the disease in COPD, in particular mortality. The fact that Pi10 was lower in patients with severe COPD receiving triple therapy might hint toward additional effects of this functional therapy on airway remodeling.

REGISTRATION

ClinicalTrials.gov, Identifier: NCT01245933.

µCT to quantify muco-obstructive lung disease and effects of neutrophil elastase knockout in mice

Muco-obstructive lung diseases are characterized by airway obstruction and hyperinflation, which can be quantified by imaging. Our aim was to evaluate µCT for longitudinal quantification of muco-obstructive lung disease in β-epithelial Na⁺ channel overexpressing (Scnn1b-TG) mice and of the effects of neutrophil elastase (NE) knockout on its progression. Lungs from wild-type (WT), NE^-/-, Scnn1b-TG, and Scnn1b-TG/NE^-/- mice were scanned with 9 µm resolution at 0, 5, 14 and 60 days of age, and airway and parenchymal disease was quantified. Mucus adhesion lesions (MAL) were persistently increased in Scnn1b-TG compared to WT mice from 0 days (20.25±6.50 vs. 9.60±2.07, P<0.05), and this effect was attenuated in Scnn1b-TG/NE^-/- mice (5.33±3.67, P<0.001). Airway wall area percentage (WA%) was increased in Scnn1b-TG mice compared to WT from 14 days onward (59.2±6.3% vs. 49.8±9.0%, P<0.001) but was similar in Scnn1b-TG/NE^-/- compared to WT at 60 days (46.4±9.2% vs. 45.4±11.5%, P=0.97). Air proportion (Air%) and mean linear intercept (Lm) were persistently increased in Scnn1b-TG compared to WT from 5 days on (53.9±4.5% vs. 30.0±5.5% and 78.82±8.44µm vs. 65.66±4.15µm, respectively, P<0.001), whereas in Scnn1b-TG/NE^-/- Air% and Lm were similar to WT from birth (27.7±5.5% vs.27.2±5.9%, P =0.92 and 61.48±9.20µm vs. 61.70±6.73µm, P=0.93, respectively). Our results suggest that µCT is sensitive to detect the onset and progression of muco-obstructive lung disease and effects of genetic deletion of NE on morphology of airways and lung parenchyma in Scnn1b-TG mice, and that it may serve as a sensitive endpoint for preclinical studies of novel therapeutic interventions for muco-obstructive lung diseases.

Feasibility of bronchial wall quantification in low- and ultralow-dose third-generation dual-source CT: An ex vivo lung study

PURPOSE

To investigate image quality and bronchial wall quantification in low- and ultralow-dose third-generation dual-source computed tomography (CT).

METHODS

A lung specimen from a formerly healthy male was scanned using third-generation dual-source CT at standard-dose (51 mAs/120 kV, CTDIvol 3.41 mGy), low-dose (1/4th and 1/10th of standard dose), and ultralow-dose setting (1/20th). Low kV (70, 80, 90, and Sn100 kV) scanning was applied in each low/ultralow-dose setting, combined with adaptive mAs to keep a constant dose. Images were reconstructed at advanced modeled iterative reconstruction (ADMIRE) levels 1, 3, and 5 for each scan. Bronchial wall were semi-automatically measured from the lobar level to subsegmental level. Spearman correlation analysis was performed between bronchial wall quantification (wall thickness and wall area percentage) and protocol settings (dose, kV, and ADMIRE). ANOVA with a post hoc pairwise test was used to compare signal-to-noise ratio (SNR), noise and bronchial wall quantification values among standard- and low/ultralow-dose settings, and among ADMIRE levels.

RESULTS

Bronchial wall quantification had no correlation with dose level, kV, or ADMIRE level (|correlation coefficients| < 0.3). SNR and noise showed no statistically significant differences at different kV in the same ADMIRE level (1, 3, or 5) and in the same dose group (P > 0.05). Generally, there were no significant differences in bronchial wall quantification among the standard- and low/ultralow-dose settings, and among different ADMIRE levels (P > 0.05).

CONCLUSION

The combined use of low/ultralow-dose scanning and ADMIRE does not influence bronchial wall quantification compared to standard-dose CT. This specimen study suggests the potential that an ultralow-dose scan can be used for bronchial wall quantification.

Influence of acquisition settings and radiation exposure on CT lung densitometry-An anthropomorphic ex vivo phantom study

Objectives

To systematically evaluate the influence of acquisition settings in conjunction with raw-data based iterative image reconstruction (IR) on lung densitometry based on multi-row detector computed tomography (CT) in an anthropomorphic chest phantom.

Materials and methods

Ten porcine heart-lung explants were mounted in an ex vivo chest phantom shell, six with highly and four with low attenuating chest wall. CT (Somatom Definition Flash, Siemens Healthineers) was performed at 120kV_p and 80kV_p, each combined with current-time products of 120, 60, 30, and 12mAs, and was reconstructed with filtered back projection (FBP) and IR (Safire, Siemens Healthineers). Mean lung density (LD), air density (AD) and noise were measured by semi-automated region-of interest (ROI) analysis, with 120kV_p/120 mAs serving as the standard of reference.

Results

Using IR, noise in lung parenchyma was reduced by ~ 31% at high attenuating chest wall and by ~ 22% at low attenuating chest wall compared to FBP, respectively (p<0.05). IR induced changes in the order of ±1 HU to mean absolute LD and AD compared to corresponding FBP reconstructions which were statistically significant (p<0.05).

Conclusions

Densitometry is influenced by acquisition parameters and reconstruction algorithms to a degree that may be clinically negligible. However, in longitudinal studies and clinical research identical protocols and potentially other measures for calibration may be required.

Computed Tomography Imaging for Novel Therapies of Chronic Obstructive Pulmonary Disease

Novel therapeutic options in chronic obstructive pulmonary disease (COPD) require delicate patient selection and thus demand for expert radiologists visually and quantitatively evaluating high-resolution computed tomography (CT) with additional functional acquisitions such as paired inspiratory-expiratory scans or dynamic airway CT. The differentiation between emphysema-dominant and airway-dominant COPD phenotypes by imaging has immediate clinical value for patient management. Assessment of emphysema severity, distribution patterns, and fissure integrity are essential for stratifying patients for different surgical and endoscopic lung volume reduction procedures. This is supported by quantitative software-based postprocessing of CT data sets, which delivers objective emphysema and airway remodelling metrics. However, the significant impact of scanning and reconstruction parameters, as well as intersoftware variability still hamper comparability between sites and studies. In earlier stage COPD imaging, it is less clear as to what extent quantitative CT might impact decision making and therapy follow-up, as emphysema progression is too slow to realistically be useful as a mid-term outcome measure in an individual, and longitudinal data on airway remodelling are still very limited.

Validation of the RayStation Monte Carlo dose calculation algorithm using a realistic lung phantom

PURPOSE

Our purposes are to compare the accuracy of RaySearch's analytical pencil beam (APB) and Monte Carlo (MC) algorithms for clinical proton therapy and to present clinical validation data using a novel animal tissue lung phantom.

METHODS

We constructed a realistic lung phantom composed of a rack of lamb resting on a stack of rectangular natural cork slabs simulating lung tissue. The tumor was simulated using 70% lean ground lamb meat inserted in a spherical hole with diameter 40 ± 5 mm carved into the cork slabs. A single-field plan using an anterior beam and a two-field plan using two anterior-oblique beams were delivered to the phantom. Ion chamber array measurements were taken medial and distal to the tumor. Measured doses were compared with calculated RayStation APB and MC calculated doses.

RESULTS

Our lung phantom enabled measurements with the MatriXX PT at multiple depths in the phantom. Using the MC calculations, the 3%/3 mm gamma index pass rates, comparing measured with calculated doses, for the distal planes were 74.5% and 85.3% for the APB and 99.1% and 92% for the MC algorithms. The measured data revealed up to 46% and 30% underdosing within the distal regions of the target volume for the single and the two field plans when APB calculations are used. These discrepancies reduced to less than 18% and 7% respectively using the MC calculations.

CONCLUSIONS

RaySearch Laboratories' Monte Carlo dose calculation algorithm is superior to the pencil-beam algorithm for lung targets. Clinicians relying on the analytical pencil-beam algorithm should be aware of its pitfalls for this site and verify dose prior to delivery. We conclude that the RayStation MC algorithm is reliable and more accurate than the APB algorithm for lung targets and therefore should be used to plan proton therapy for patients with lung cancer.

Influence of Inspiratory/Expiratory CT Registration on Quantitative Air Trapping

RATIONALE AND OBJECTIVES

The aim of this study was to assess variability in quantitative air trapping (QAT) measurements derived from spatially aligned expiration CT scans.

MATERIALS AND METHODS

Sixty-four paired CT examinations, from 16 school-age cystic fibrosis subjects examined at four separate time intervals, were used in this study. For each pair, visually inspected lobe segmentation maps were generated and expiration CT data were registered to the inspiration CT frame. Measurements of QAT, the percentage of voxels on the expiration CT scan below a set threshold were calculated for each lobe and whole-lung from the registered expiration CT and compared to the true values from the unregistered data.

RESULTS

A mathematical model, which simulates the effect of variable regions of lung deformation on QAT values calculated from aligned to those from unaligned data, showed the potential for large bias. Assessment of experimental QAT measurements using Bland-Altman plots corroborated the model simulations, demonstrating biases greater than 5% when QAT was approximately 40% of lung volume. These biases were removed when calculating QAT from aligned expiration CT data using the determinant of the Jacobian matrix. We found, by Dice coefficient analysis, good agreement between aligned expiration and inspiration segmentation maps for the whole-lung and all but one lobe (Dice coefficient > 0.9), with only the lingula generating a value below 0.9 (mean and standard deviation of 0.85 ± 0.06).

CONCLUSION

The subtle and predictable variability in corrected QAT observed in this study suggests that image registration is reliable in preserving the accuracy of the quantitative metrics.

Reproducibility of an airway tapering measurement in computed tomography with application to bronchiectasis

We propose a pipeline to acquire a scalar tapering measurement from the carina to the most distal point of an individual airway visible on computed tomography (CT). We show the applicability of using tapering measurements on clinically acquired data by quantifying the reproducibility of the tapering measure. We generate a spline from the centerline of an airway to measure the area and arclength at contiguous intervals. The tapering measurement is the gradient of the linear regression between area in log space and arclength. The reproducibility of the measure was assessed by analyzing different radiation doses, voxel sizes, and reconstruction kernel on single timepoint and longitudinal CT scans and by evaluating the effect of airway bifurcations. Using 74 airways from 10 CT scans, we show a statistical difference, p = 3.4 × 10 - 4 , in tapering between healthy airways ( n = 35 ) and those affected by bronchiectasis ( n = 39 ). The difference between the mean of the two populations is 0.011    mm - 1 , and the difference between the medians of the two populations was 0.006    mm - 1 . The tapering measurement retained a 95% confidence interval of ± 0.005    mm - 1 in a simulated 25 mAs scan and retained a 95% confidence of ± 0.005    mm - 1 on simulated CTs up to 1.5 times the original voxel size. We have established an estimate of the precision of the tapering measurement and estimated the effect on precision of the simulated voxel size and CT scan dose. We recommend that the scanner calibration be undertaken with the phantoms as described, on the specific CT scanner, radiation dose, and reconstruction algorithm that are to be used in any quantitative studies.

Longitudinal airway remodeling in active and past smokers in a lung cancer screening population

OBJECTIVES

To longitudinally investigate smoking cessation-related changes of quantitative computed tomography (QCT)-based airway metrics in a group of heavy smokers.

METHODS

CT scans were acquired in a lung cancer screening population over 4 years at 12-month intervals in 284 long-term ex-smokers (ES), 405 continuously active smokers (CS), and 31 subjects who quitted smoking within 2 years after baseline CT (recent quitters, RQ). Total diameter (TD), lumen area (LA), and wall percentage (WP) of 1st-8th generation airways were computed using airway analysis software. Inter-group comparison was performed using Mann-Whitney U test or Student's t test (two groups), and ANOVA or ANOVA on ranks with Dunn's multiple comparison test (more than two groups), while Fisher's exact test or chi-squared test was used for categorical data. Multiple linear regression was used for multivariable analysis.

RESULTS

At any time, TD and LA were significantly higher in ES than CS, for example, in 5th-8th generation airways at baseline with 6.24 mm vs. 5.93 mm (p < 0.001) and 15.23 mm² vs. 13.51 mm² (p < 0.001), respectively. RQ showed higher TD (6.15 mm vs. 5.93 mm, n.s.) and significantly higher LA (14.77 mm² vs. 13.51 mm², p < 0.001) than CS after 3 years, and after 4 years. In multivariate analyses, smoking status independently predicted TD, LA, and WP at baseline, at 3 years and 4 years (p < 0.01-0.001), with stronger impact than pack years.

CONCLUSIONS

Bronchial dimensions depend on the smoking status. Smoking-induced airway remodeling can be partially reversible after smoking cessation even in long-term heavy smokers. Therefore, QCT-based airway metrics in clinical trials should consider the current smoking status besides pack years.

KEY POINTS

• Airway lumen and diameter are decreased in active smokers compared to ex-smokers, and there is a trend towards increased airway wall thickness in active smokers. • Smoking-related airway changes improve within 2 years after smoking cessation. • Smoking status is an independent predictor of airway dimensions.

Quantitative CT detects changes in airway dimensions and air-trapping after bronchial thermoplasty for severe asthma

OBJECTIVES

Bronchial thermoplasty (BT) can be considered in the treatment of severe asthma to reduce airway smooth muscle mass and bronchoconstriction. We hypothesized that BT may thus have long-term effects on airway dimensions and air-trapping detectable by quantitative computed tomography (QCT).

METHODS

Paired in- and expiratory CT and inspiratory CT were acquired in 17 patients with severe asthma before and up to two years after bronchial thermoplasty and in 11 additional conservatively treated patients with serve asthma, respectively. A fully automatic software calculated the airways metrics for wall thickness (WT), wall percentage (WP), lumen area (LA) and total diameter (TD). Furthermore, lung air-trapping was quantified by determining the quotient of mean lung attenuation in expiration vs. inspiration (E/I MLA) and relative volume change in the Hounsfield interval -950 to -856 in expiration to inspiration (RVC_856-950) in a generation- and lobe-based approach, respectively.

RESULTS

BT reduced WT for the combined analysis of the 2nd-7th airway generation significantly by 0.06 mm (p = 0.026) and WP by 2.05% (p < 0.001), whereas LA and TD did not change significantly (p = 0.147, p = 0.706). No significant changes were found in the control group. Furthermore, E/I MLA and RVC_856-950 decreased significantly after BT by 12.65% and 1.77% (p < 0.001), respectively.

CONCLUSION

BT significantly reduced airway narrowing and air-trapping in patients with severe asthma. This can be interpreted as direct therapeutic effects caused by a reduction in airway-smooth muscle mass and changes in innervation. A reduction in air-trapping indicates an influence on more peripheral airways not directly treated by the BT procedure.

Pulmonary quantitative CT imaging in focal and diffuse disease: current research and clinical applications

The frenetic development of imaging technology-both hardware and software-provides exceptional potential for investigation of the lung. In the last two decades, CT was exploited for detailed characterization of pulmonary structures and description of respiratory disease. The introduction of volumetric acquisition allowed increasingly sophisticated analysis of CT data by means of computerized algorithm, namely quantitative CT (QCT). Hundreds of thousands of CTs have been analysed for characterization of focal and diffuse disease of the lung. Several QCT metrics were developed and tested against clinical, functional and prognostic descriptors. Computer-aided detection of nodules, textural analysis of focal lesions, densitometric analysis and airway segmentation in obstructive pulmonary disease and textural analysis in interstitial lung disease are the major chapters of this discipline. The validation of QCT metrics for specific clinical and investigational needs prompted the translation of such metrics from research field to patient care. The present review summarizes the state of the art of QCT in both focal and diffuse lung disease, including a dedicated discussion about application of QCT metrics as parameters for clinical care and outcomes in clinical trials.