2471

Probing changes in lung physiology in COPD using CT, perfusion MRI and hyperpolarized xenon-129 MRI
Kun Qing1, Nicholas J. Tustison1, John P. Mugler, III1, Jaime F. Mata1, Zixuan Lin1, Li Zhao2, Da Wang3, Xue Feng1, Kai Ruppert4, Talissa A. Altes5, Joanne M. Cassani5, and Y. Michael Shim1

1University of Virginia, Charlottesville, VA, United States, 2Beth Israel Deaconess Medical Center, Boston, MA, United States, 3University of Washington, Seattle, WA, United States, 4University of Pennsylvania, Philadelphia, VA, United States, 5University of Missouri School of Medicine, Columbia, MO, United States

Synopsis

In this study, by using chest CT, Gadolinium-enhanced perfusion MRI, and hyperpolarized xenon-129 ventilation and gas uptake MRI, we assessed the quantitative changes in tissue density, pulmonary perfusion and gas uptake in patients with COPD compared to normal subjects. We found evidence for compensatory pulmonary vasoconstriction to match impairment of ventilation, and also pulmonary shunt and dead space. By incorporating a new lobar segmentation method for proton MRI, we performed statistical analysis to evaluate the regional interrelationships among different measures. We demonstrated that xenon-129 MRI has high potential to identify changes of multiple aspects of lung physiology in one acquisition.

Introduction

Clinical diagnosis and prognosis of chronic obstructive pulmonary disease (COPD) are still primarily based on irreversible airflow limitation measured by spirometry. However, multiple studies reported [1,2] heterogeneous phenotypes of COPD that are far more complicated than being simply defined by irreversible airflow limitation. In this pilot study, we exploited multiple imaging modalities, including CT, gadolinium-enhanced perfusion MRI, and hyperpolarized xenon-129 MRI [3,4], to perform a detailed investigation of changes in lung morphology and functions in COPD. Utility and strengths of xenon-129 MRI in assessing COPD were also evaluated against the other imaging modalities.

Methods

Chest CT was done based on the COPDGene protocol [5]. Lung regions were divided into five lobes based on the detection of lung fissures. The bottom 15th percentile (Perc15) of the CT Hounsfield Unit (HU) numbers was calculated for each lobe to represent the average severity of emphysema in the lobe [6]. All MR acquisitions were done using a 1.5T commercial whole-body MR scanner (Avanto; Siemens Medical Solutions, Malvern PA). For perfusion MRI, Gd-DTPA was injected at the standard dose (0.2ml/kg body weight), simultaneously with data acquisition. A Siemens TWIST sequence was used to acquire a four-dimensional data set with a real temporal resolution of 1.6-2 seconds during free breathing. First-pass dynamic Gd-enhanced images were processed as described in [7] to produce pulmonary perfusion volume (PBV), perfusion flow (PBF), and mean transit time (MTT). For Xe129 ventilation and gas uptake imaging, a combined Xe129 and proton MR acquisition during a single breath was used to enable regional quantitative analysis of lung function. Xe129 gas uptake measures, including tissue-to-gas, RBC-to-gas and RBC-to-tissue ratios were processed as described in [3], while the percentage of the lung with reduced ventilation (%VD) was calculated based on [8]. A lobar segmentation method [9] was used to divide the proton images into individual lobes, and this segmentation information was then applied to the lung function maps to calculate mean perfusion, ventilation or gas uptake measures for each lobe for statistical analysis.

Results and Discussions

Subject demographics, PFT and 6MW results are shown in Table 1. Lobar mean values of the different measures were shown in Table 2. No within-group statistical difference among lobes was observed for any of the imaging measures in the AMN or COPD group (P>0.05). However, comparing between groups, COPD subjects had significantly lower CT Perc15 (P<0.001), MRI perfusion measure PBF (P=0.047), and Xe129 gas uptake measures (P<0.001). Also the %VD in COPD subjects was significantly higher than that in the AMN group (P=0.006). Figure 1 shows representative lobar segmentation for CT and MR images/maps for a slice from subject H4. For the Xe129 MRI results, it can be seen that the ventilation (“gas”) image and the generated tissue-to-gas ratio map align with the proton images very well. Correlations of lobar mean values from different imaging measures are shown in Table 4. CT Perc15 had strong correlations with Xe129 tissue-to-gas (R=0.74, P<0.001) and RBC-to-gas ratios (R=0.62, P<0.001). This is as expected since emphysema or hyper-expansion could both cause decrease of tissue density and/or blood volume. A good example of this is Area 1 of subject C4, shown in Figure 2. With severe emphysema, Area 1 was suspected to have very little functioning tissue left as shown on CT. Both ventilation (Xe129 “gas” image) and perfusion (PBV) were very low, which would lead to virtually no gas uptake in this area. There are strong correlations between perfusion measures and the Xe129 %VD (PBV: R=-0.59, P<0.001; PBF: R=-0.62, P<0.001; MTT: R=0.59, P<0.001). This may mean perfusion variation from compensatory vasoconstriction response to ventilation changes to maintain V/Q coupling [10]. A supporting example was Area 2 from subject C4 shown in Figure 2. Lung tissue in the Area 2 appears to be intact based on CT, but ventilation and PBF were abnormally low. Mismatches between ventilation and perfusion were also observed. For example, Area 3 in subject C4 appeared normal by chest CT, and was well ventilated, but had a fissure-like perfusion defect, and abnormally low gas uptake as measured by Xe129 DP MRI. While Area 4 in subject C2 exhibited normal lung tissue without emphysema by chest CT and normal perfusion, but had poor ventilation and poor gas exchange. As expected, the Xe129 RBC-to-gas ratio, an imaging counterpart of DLCO, correlated strongly with all other imaging methods: CT Perc15 (R=0.62, P<0.001), perfusion (e.g. MTT: R=-0.52, P<0.001) and ventilation (R=-0.64, P<0.001).

Conclusion

Xe129 MRI has high potential to identify changes of multiple aspects of lung physiology in one acquisition.

Acknowledgements

This work is supported by NIH (R21HL129112, R01HL109618, R01HL132287, and R01HL132177 from the National Heart, Lung, and Blood Institute) and research funding from Siemens Medical Solutions.

References

[1]. Kim V, et. Al. Chest 2011; 140: 626-633. [2]. Vestbo J, et. Al. Am J Respir Crit Care Med. 2014; 189: 1022-1030. [3]. Qing K, et. Al. JMRI 2014; 39: 346-359. [4]. Kaushik SS, et. Al. Magn. Reson. in Med. 2016; 75: 1434-1443. [5]. Labaki WW, et. Al.. Am J Respir Crit Care Med. 2017. doi: 10.1164/rccm.201703-0451PP. [Epub ahead of print] [6]. Mohamed Hoesein FA, et. Al. Thorax 2011; 66: 782-787. [7]. Ohno Y, et. Al. JMRI 2004; 20: 353-365. [8]. Tustison NJ, et. Al. JMRI 2011; 34: 831-841. [9]. Tustison NJ, et. Al. Magn. Reson. in Med. . 2016 Jul;76(1):315-20. [10]. Stenmark KR, et. Al.. Circulation research 2006; 99: 675-691.

Figures



Figure 1. Representative lobar segmentation of Xe129 MRI, Gd-enhanced perfusion MRI, and CT images. Left: top: ventilation image (“Gas”) overlaid on proton image acquired in the same breath hold; middle: generated tissue-to-gas ratio map on top of proton image; bottom: lobar segmentation for proton image . Middle: top: Gd-enhanced perfusion MRI image of the lung; middle: generated pulmonary blood volume map (PBV) in the lung; bottom: lobar segmentation for Gd-enhanced perfusion MRI image. Right: top: chest CT with regions lower than -950 HU marked by red dots; bottom: lobar segmentation for CT was done based on fissure detection.

Figure 2. Representative CT image, PBV, and Xe129 tissue-to-gas and RBC-to-gas ratio maps from subject C4 and C2. Alterations of lung function exhibited different patterns in different areas of the lung, as shown in areas 1-4 in the figure.


Proc. Intl. Soc. Mag. Reson. Med. 26 (2018)
2471