Dante P. I. Capaldi

Globally optimal co-segmentation of three-dimensional pulmonary 1 H and hyperpolarized 3 He MRI with spatial consistence prior

Pulmonary imaging using hyperpolarized (3)He/(129)Xe gas is emerging as a new way to understand the regional nature of pulmonary ventilation abnormalities in obstructive lung diseases. However, the quantitative information derived is completely dependent on robust methods to segment both functional and structural/anatomical data. Here, we propose an approach to jointly segment the lung cavity from (1)H and (3)He pulmonary magnetic resonance images (MRI) by constraining the spatial consistency of the two segmentation regions, which simultaneously employs the image features from both modalities. We formulated the proposed co-segmentation problem as a coupled continuous min-cut model and showed that this combinatorial optimization problem can be solved globally and exactly by means of convex relaxation. In particular, we introduced a dual coupled continuous max-flow model to study the convex relaxed coupled continuous min-cut model under a primal and dual perspective. This gave rise to an efficient duality-based convex optimization algorithm. We implemented the proposed algorithm in parallel using general-purpose programming on graphics processing unit (GPGPU), which substantially increased its computational efficiency. Our experiments explored a clinical dataset of 25 subjects with chronic obstructive pulmonary disease (COPD) across a wide range of disease severity. The results showed that the proposed co-segmentation approach yielded superior performance compared to single-channel image segmentation in terms of precision, accuracy and robustness.

Pulmonary Imaging Biomarkers of Gas Trapping and Emphysema in COPD: (3)He MR Imaging and CT Parametric Response Maps.

PURPOSE To directly compare magnetic resonance (MR) imaging and computed tomography (CT) parametric response map (PRM) measurements of gas trapping and emphysema in ex-smokers both with and without chronic obstructive pulmonary disease (COPD). MATERIALS AND METHODS Participants provided written informed consent to a protocol that was approved by a local research ethics board and Health Canada and was compliant with the HIPAA (Institutional Review Board Reg. #00000940). The prospectively planned study was performed from March 2014 to December 2014 and included 58 ex-smokers (mean age, 73 years ± 9) with (n = 32; mean age, 74 years ± 7) and without (n = 26; mean age, 70 years ± 11) COPD. MR imaging (at functional residual capacity plus 1 L), CT (at full inspiration and expiration), and spirometry or plethysmography were performed during a 2-hour visit to generate ventilation defect percent (VDP), apparent diffusion coefficient (ADC), and PRM gas trapping and emphysema measurements. The relationships between pulmonary function and imaging measurements were determined with analysis of variance (ANOVA), Holm-Bonferroni corrected Pearson correlations, multivariate regression modeling, and the spatial overlap coefficient (SOC). RESULTS VDP, ADC, and PRM gas trapping and emphysema (ANOVA, P < .001) measurements were significantly different in healthy ex-smokers than they were in ex-smokers with COPD. In all ex-smokers, VDP was correlated with PRM gas trapping (r = 0.58, P < .001) and with PRM emphysema (r = 0.68, P < .001). VDP was also significantly correlated with PRM in ex-smokers with COPD (gas trapping: r = 0.47 and P = .03; emphysema: r = 0.62 and P < .001) but not in healthy ex-smokers. In a multivariate model that predicted PRM gas trapping, the forced expiratory volume in 1 second normalized to the forced vital capacity (standardized coefficients [βS] = -0.69, P = .001) and airway wall area percent (βS = -0.22, P = .02) were significant predictors. PRM emphysema was predicted by the diffusing capacity for carbon monoxide (βS = -0.29, P = .03) and VDP (βS = 0.41, P = .001). Helium 3 ADC values were significantly elevated in PRM gas-trapping regions (P < .001). The spatial relationship for ventilation defects was significantly greater with PRM gas trapping than with PRM emphysema in patients with mild (for gas trapping, SOC = 36% ± 28; for emphysema, SOC = 1% ± 2; P = .001) and moderate (for gas trapping, SOC = 34% ± 28; for emphysema, SOC = 7% ± 15; P = .006) COPD. For severe COPD, the spatial relationship for ventilation defects with PRM emphysema (SOC = 64% ± 30) was significantly greater than that for PRM gas trapping (SOC = 36% ± 18; P = .01). CONCLUSION In all ex-smokers, ADC values were significantly elevated in regions of PRM gas trapping, and VDP was quantitatively and spatially related to both PRM gas trapping and PRM emphysema. In patients with mild to moderate COPD, VDP was related to PRM gas trapping, whereas in patients with severe COPD, VDP correlated with both PRM gas trapping and PRM emphysema.

Regional Heterogeneity of Chronic Obstructive Pulmonary Disease Phenotypes: Pulmonary (3)He Magnetic Resonance Imaging and Computed Tomography.

Abstract Pulmonary ventilation may be visualized and measured using hyperpolarized 3He magnetic resonance imaging (MRI) while emphysema and its distribution can be quantified using thoracic computed tomography (CT). Our objective was to phenotype ex-smokers with COPD based on the apical-to-basal distribution of ventilation abnormalities and emphysema to better understand how these phenotypes change regionally as COPD progresses. We evaluated 100 COPD ex-smokers who provided written informed consent and underwent spirometry, CT and 3He MRI. 3He MRI ventilation imaging was used to quantify the ventilation defect percent (VDP) for whole-lung and individual lung lobes. Regional VDP was used to generate the apical-lung (AL)-to-basal-lung (BL) difference (ΔVDP); a positive ΔVDP indicated AL-predominant and negative ΔVDP indicated BL-predominant ventilation defects. Emphysema was quantified using the relative-area-of-the-lung ≤−950HU (RA950) of the CT density histogram for whole-lung and individual lung lobes. The AL-to-BL RA950 difference (ΔRA950) was generated with a positive ΔRA950 indicating AL-predominant emphysema and a negative ΔRA950 indicating BL-predominant emphysema. Seventy-two ex-smokers reported BL-predominant MRI ventilation defects and 71 reported AL-predominant CT emphysema. BL-predominant ventilation defects (AL/BL: GOLD I = 18%/82%, GOLD II = 24%/76%) and AL-predominant emphysema (AL/BL: GOLD I = 84%/16%, GOLD II = 72%/28%) were the major phenotypes in mild-moderate COPD. In severe COPD there was a more uniform distribution for ventilation defects (AL/BL: GOLD III = 40%/60%, GOLD IV = 43%/57%) and emphysema (AL/BL: GOLD III = 64%/36%, GOLD IV = 43%/57%). Basal-lung ventilation defects predominated in mild-moderate GOLD grades, and a more homogeneous distribution of ventilation defects was observed in more advanced grade COPD; these differences suggest that over time, regional ventilation abnormalities become more homogenously distributed during disease progression.

Anatomical pulmonary magnetic resonance imaging segmentation for regional structure-function measurements of asthma

PURPOSE Pulmonary magnetic-resonance-imaging (MRI) and x-ray computed-tomography have provided strong evidence of spatially and temporally persistent lung structure-function abnormalities in asthmatics. This has generated a shift in their understanding of lung disease and supports the use of imaging biomarkers as intermediate endpoints of asthma severity and control. In particular, pulmonary (1)H MRI can be used to provide quantitative lung structure-function measurements longitudinally and in response to treatment. However, to translate such biomarkers of asthma, robust methods are required to segment the lung from pulmonary (1)H MRI. Therefore, their objective was to develop a pulmonary (1)H MRI segmentation algorithm to provide regional measurements with the precision and speed required to support clinical studies. METHODS The authors developed a method to segment the left and right lung from (1)H MRI acquired in 20 asthmatics including five well-controlled and 15 severe poorly controlled participants who provided written informed consent to a study protocol approved by Health Canada. Same-day spirometry and plethysmography measurements of lung function and volume were acquired as well as (1)H MRI using a whole-body radiofrequency coil and fast spoiled gradient-recalled echo sequence at a fixed lung volume (functional residual capacity + 1 l). We incorporated the left-to-right lung volume proportion prior based on the Potts model and derived a volume-proportion preserved Potts model, which was approximated through convex relaxation and further represented by a dual volume-proportion preserved max-flow model. The max-flow model led to a linear problem with convex and linear equality constraints that implicitly encoded the proportion prior. To implement the algorithm, (1)H MRI was resampled into ∼3 × 3 × 3 mm(3) isotropic voxel space. Two observers placed seeds on each lung and on the background of 20 pulmonary (1)H MR images in a randomized dataset, on five occasions, five consecutive days in a row. Segmentation accuracy was evaluated using the Dice-similarity-coefficient (DSC) of the segmented thoracic cavity with comparison to five-rounds of manual segmentation by an expert observer. The authors also evaluated the root-mean-squared-error (RMSE) of the Euclidean distance between lung surfaces, the absolute, and percent volume error. Reproducibility was measured using the coefficient of variation (CoV) and intraclass correlation coefficient (ICC) for two observers who repeated segmentation measurements five-times. RESULTS For five well-controlled asthmatics, forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) was 83% ± 7% and FEV1 was 86 ± 9%pred. For 15 severe, poorly controlled asthmatics, FEV1/FV C = 66% ± 17% and FEV1 = 72 ± 27%pred. The DSC for algorithm and manual segmentation was 91% ± 3%, 92% ± 2% and 91% ± 2% for the left, right, and whole lung, respectively. RMSE was 4.0 ± 1.0 mm for each of the left, right, and whole lung. The absolute (percent) volume errors were 0.1 l (∼6%) for each of right and left lung and ∼0.2 l (∼6%) for whole lung. Intra- and inter-CoV (ICC) were <0.5% (>0.91%) for DSC and <4.5% (>0.93%) for RMSE. While segmentation required 10 s including ∼6 s for user interaction, the smallest detectable difference was 0.24 l for algorithm measurements which was similar to manual measurements. CONCLUSIONS This lung segmentation approach provided the necessary and sufficient precision and accuracy required for research and clinical studies.

Second-order Texture Measurements of 3He Ventilation MRI:

RATIONALE AND OBJECTIVES (3)He magnetic resonance imaging (MRI) can be used to quantify functional responses to asthma therapy and provocation. Ventilation imaging offers quantitative information beyond ventilation defects that have not yet been exploited. Therefore, our objective was to evaluate hyperpolarized (3)He MRI ventilation defect percent (VDP) and compare this and pulmonary function measurements to ventilation image texture features and their changes post-bronchodilator administration in patients with asthma. MATERIALS AND METHODS Volunteers with a diagnosis of asthma provided written informed consent to an ethics board-approved protocol and underwent pulmonary function tests and MRI before and after salbutamol inhalation. MR images were analyzed using VDP, and their texture was evaluated via gray-level run-length matrices. These texture classifiers were compared to VDP in responders to bronchodilation based on VDP (VDP responders) and forced expiratory volume in 1 s (FEV1) (FEV1 responders). RESULTS In total, 47 patients with asthma (18 males 39 ± 13 years, FEV1 = 79 ± 21%) reported significantly improved FEV1, FEV1/forced vital capacity (FVC), residual volume (RV)/total lung capacity (TLC) (all P = .0001) and VDP (P = .01) post-salbutamol. Post-salbutamol, VDP responders and nonresponders to salbutamol were significantly different for coarse-texture features including long-run emphasis (LRE) and long-run, low gray-level emphasis (LRLGE, both P < .05) and for FEV1 responders to salbutamol, there was significantly different long-run, high gray-level emphasis (LRHGE, P = .04). There were significant relationships for VDP with LRE (R = .50, P = .0003), LRLGE (R = .34, P = .02), and LRHGE (R = .56, P = .0001). Receiver operating characteristic curves showed VDP with the strongest performance (AUC = .92), followed by coarse-texture classifier LRHGE (AUC = .83), FEV1 (AUC = .80), LRE (AUC = .66), FVC (AUC = .58), and LRLGE (AUC = .42). CONCLUSIONS In patients with asthma, differences in ventilation patchiness post-salbutamol can be quantified using coarse-texture classifiers that are significantly different in bronchodilator responders.

Ultrashort echo time MRI biomarkers of asthma.

To develop and assess ultrashort echo‐time (UTE) magnetic resonance imaging (MRI) biomarkers of lung function in asthma patients.

Magnetic resonance imaging biomarkers of chronic obstructive pulmonary disease prior to radiation therapy for non-small cell lung cancer

Highlights • Three imaging phenotypes of COPD and ventilation heterogeneity.• We examine relationships for non-tumour lobe ventilation voids and clinical tests.• Smoking history and airflow obstruction were diagnostics for imaging phenotypes.

Thoracic CT‐MRI coregistration for regional pulmonary structure–function measurements of obstructive lung disease

Purpose Recent pulmonary imaging research has revealed that in patients with chronic obstructive pulmonary disease (COPD) and asthma, structural and functional abnormalities are spatially heterogeneous. This novel information may help optimize treatment in individual patients, monitor interventional efficacy, and develop new treatments. Moreover, by automating the measurement of regional biomarkers for the 19 different anatomical lung segments, there is an opportunity to embed imaging biomarkers into clinically acceptable clinical workflows and improve lung disease clinical care. Therefore, to exploit the regional structure–function information provided by thoracic imaging, and as a first step toward this goal, our objective was to develop a fully automated registration pipeline for thoracic x‐ray computed tomography (CT) and inhaled gas functional magnetic resonance imaging (MRI) whole lung and segmental structure–function biomarkers. Methods Thirty‐five patients including 15 severe, poorly controlled asthmatics and 20 COPD patients [classified according to the global initiative for chronic obstructive lung disease (GOLD) criteria)] provided written informed consent to a study protocol approved by Health Canada and underwent pulmonary function tests, MRI, and CT during a single 2‐hour visit. Using this diverse patient dataset, we developed and evaluated a joint deformable registration approach to simultaneously coregister CT with both 1H and 3He MRI by enforcing the similarity of the deformation fields from the two individual registrations. We derived a simpler model that was equivalent to the original challenging optimization problem through variational analysis and the simpler model gave rise to an efficient numerical solver that was parallelized on a graphics processing unit. The coregistered CT‐3He MRI and whole lung/segmental lung masks were used to generate whole lung and segmental 3He MRI ventilation defect percent (VDP). To estimate fiducial localization reproducibility, a single observer manually identified 109 pairs of CT and 3He MRI fiducials for 35 patient images on five separate occasions and determined the fiducial localization error (FLE). CT‐3He MRI registration accuracy was evaluated using the target registration error (TRE). Whole lung VDP generated using the algorithm was compared with VDP generated using a previously validated semiautomated approach and computational efficiency was evaluated using run time. Results In 35 patients including 15 with severe asthma and 20 with COPD, mean forced expiratory volume in 1 s (FEV1) was 63±24%pred and FEV1/forced vital capacity (FVC) was 54 ± 17%. FLE was 0.16 mm and 0.34 mm for 3He MRI and CT, respectively. TRE was 4.5 ± 2.0 mm, 4.0 ± 1.7 mm, 4.8 ± 2.3 mm for asthma, COPD GOLD II, and GOLD III groups, respectively, with a mean of 4.4 ± 2.0 mm for the entire dataset. TRE was significantly improved for joint CT‐1H/3He MRI registration compared with CT‐1H MRI rigid registration (P < 0.0001). Whole lung VDP generated using the pipeline was not significantly different (P = 0.37) compared to a semiautomated method with which it was strongly correlated (r = 0.93, P < 0.0001). The fully automated pipeline required 11 ± 0.4 min to generate whole lung and segmental VDP. Conclusions For a diverse group of patients with COPD and asthma, whole lung and segmental VDP was measured using an automated lung image analysis pipeline which provides a way to incorporate lung functional biomarkers into clinical research and patient care.

Free-breathing Functional Pulmonary MRI: Response to Bronchodilator and Bronchoprovocation in Severe Asthma

RATIONALE AND OBJECTIVES Ventilation heterogeneity is a hallmark feature of asthma. Our objective was to evaluate ventilation heterogeneity in patients with severe asthma, both pre- and post-salbutamol, as well as post-methacholine (MCh) challenge using the lung clearance index, free-breathing pulmonary 1H magnetic resonance imaging (FDMRI), and inhaled-gas MRI ventilation defect percent (VDP). MATERIALS AND METHODS Sixteen severe asthmatics (49 ± 10 years) provided written informed consent to an ethics board-approved protocol. Spirometry, plethysmography, and multiple breath nitrogen washout to measure the lung clearance index were performed during a single visit within 15 minutes of MRI. Inhaled-gas MRI and FDMRI were performed pre- and post-bronchodilator to generate VDP. For asthmatics with forced expiratory volume in 1 second (FEV1) >70%predicted, MRI was also performed before and after MCh challenge. Wilcoxon signed-rank tests, Spearman correlations, and a repeated-measures analysis of variance were performed. RESULTS Hyperpolarized 3He (P = .02) and FDMRI (P = .02) VDP significantly improved post-salbutamol and for four asthmatics who could perform MCh (n = 4). 3He and FDMRI VDP significantly increased at the provocative concentration of MCh, resulting in a 20% decrease in FEV1 (PC20) and decreased post-bronchodilator (P = .02), with a significant difference between methods (P = .01). FDMRI VDP was moderately correlated with 3He VDP (ρ = .61, P = .01), but underestimated VDP relative to 3He VDP (-6 ± 9%). Whereas 3He MRI VDP was significantly correlated with the lung clearance index, FDMRI was not (ρ = .49, P = .06). CONCLUSIONS FDMRI VDP generated in free-breathing asthmatic patients was correlated with static inspiratory breath-hold 3He MRI VDP but underestimated VDP relative to 3He MRI VDP. Although less sensitive to salbutamol and MCh, FDMRI VDP may be considered for asthma patient evaluations at centers without inhaled-gas MRI.

TH‐CD‐202‐09: Free‐Breathing Proton MRI Functional Lung Avoidance Maps to Guide Radiation Therapy

PURPOSE Pulmonary functional MRI using inhaled gas contrast agents was previously investigated as a way to identify well-functioning lung in patients with NSCLC who are clinical candidates for radiotherapy. Hyperpolarized noble-gas (3 He and 129 Xe) MRI has also been optimized to measure functional lung information, but for a number of reasons, the clinical translation of this approach to guide radiotherapy planning has been limited. As an alternative, free-breathing pulmonary 1H MRI using clinically available MRI systems and pulse sequences provides a non-contrast-enhanced method to generate both ventilation and perfusion maps. Free-breathing 1 H MRI exploits non-rigid registration and Fourier decomposition of MRI signal intensity differences (Bauman et al., MRM, 2009) that may be generated during normal tidal breathing. Here, our objective was to generate free-breathing 1 H MRI ventilation and lung function avoidance maps in patients with NSCLC as a way to guide radiation therapy planning. METHODS Stage IIIA/IIIB NSCLC patients (n=8, 68±9yr) provided written informed consent to a randomized controlled clinical trial (https://clinicaltrials.gov/ct2/show/NCT02002052) that aimed to compare outcomes related to image-guided versus conventional radiation therapy planning. Hyperpolarized 3 He/129 Xe and dynamic free tidal-breathing 1 H MRI were acquired as previously described (Capaldi et al., Acad Radiol, 2015). Non-rigid registration was performed using the modality-independent-neighbourhood-descriptor (MIND) deformable approach (Heinrich et al., Med Image Anal, 2012). Ventilation-defect-percent (3 He:VDPHe , 129 Xe:VDPXe , Free-breathing-1 H:VDPFB ) and the corresponding ventilation maps were compared using Pearson correlation coefficients (r) and the Dice similarity coefficient (DSC). RESULTS VDPFB was significantly related to VDPHe (r=.71; p=.04) and VDPXe (r=.80; p=.01) and there were also strong spatial relationships (DSCHe /DSCXe =89±3%/77±11%). CONCLUSION In this proof of concept study in NSCLC patients, free-breathing 1 H MRI ventilation defects were quantitatively and spatially related to inhaled-noble-gas MRI ventilation defects. Free-breathing 1 H MRI measures lung function/ventilation that can be used to optimize radiotherapy planning in NSCLC patients.