0432

Imaging Regional Capillary Cardio-Pulmonary Blood Flow Dynamics using Hyperpolarized 129Xe MRI and Keyhole Reconstruction
Peter James Niedbalski1, Elianna A Bier2,3, Ziyi Wang2,3, Matthew M Willmering1, Bastiaan Driehuys2,3,4, and Zackary I Cleveland1,5
1Center for Pulmonary Imaging Research, Cincinnati Children's Hospital Medical Center, Cincinnati, OH, United States, 2Department of Biomedical Engineering, Duke University, Durham, NC, United States, 3Center for In Vivo Microscopy, Duke University Medical Center, Durham, NC, United States, 4Department of Radiology, Duke University Medical Center, Durham, NC, United States, 5Department of Pediatrics, University of Cincinnati Medical Center, Cincinnati, OH, United States

Synopsis

Hyperpolarized 129Xe MRI offers the ability to analyze pulmonary gas transfer by imaging 129Xe dissolved in red blood cells (RBCs) separately from 129Xe in other tissues. A notable feature of the dissolved 129Xe signal is the presence of small cardiogenic oscillations in the 129Xe RBC signal, which have been used to characterized global abnormalities in pulmonary microvascular hemodynamics. Here, we demonstrate that these cardiogenic oscillations can be mapped 3-dimensionally to image capillary bed hemodynamics. Our approach uses keyhole reconstruction of standard 129Xe gas exchange MR acquisitions. Metrics obtained from these maps distinguished healthy from disease cohorts and predicted disease progression.

Introduction

Cardiovascular abnormalities stemming from dysfunction in the pulmonary microvasculature are a common comorbidity or defining characteristic of a variety of lung diseases.1 While CT and MRI metrics have been used to quantify defects in larger pulmonary vasculature, sensitivity to the smaller vessels of the lungs (e.g., the gas exchange capillaries) is lacking in current imaging techniques.2 In recent years, hyperpolarized (HP) 129Xe MRI has emerged as a method to probe regional lung function, and of particular interest, MRI of HP 129Xe dissolved in red blood cells (RBCs) and other tissues has been used to image regional gas exchange using 1-Point Dixon imaging.3 These images have been shown to be sensitive to impaired gas transfer in a variety of lung diseases.4-6 Additional metrics derived from whole-lung 129Xe dynamic spectroscopy have revealed global cardiogenic oscillations in the 129Xe signal dissolved in RBCs that appear to reflect pulmonary vascular disease.7,8 Notably, these same cardiogenic oscillations are also present in the raw data of HP 129Xe gas exchange imaging (Figure 1). Herein, we show that these cardiogenic oscillations can be mapped regionally using keyhole image reconstruction9, providing images of local hemodynamics within the pulmonary capillary bed.

Methods

Healthy (N=12), idiopathic pulmonary fibrosis (IPF) (N=12), pulmonary arterial hypertension (PAH) (N=10), and left heart failure (LHF) (N=6) subjects underwent HP 129Xe 1-point Dixon gas exchange MRI after inhaling 1L of xenon. Some subjects (N=7 Healthy, N=7 IPF, N=10 PAH, N=5 LHF) also underwent whole-lung 129Xe dynamic spectroscopy for analysis of cardiogenic oscillations. Of the subjects who were imaged, healthy (N=8) and IPF (N=12) were imaged on a 1.5T scanner (GE 15M4 EXCITE), while all others were imaged on a 3T scanner (SIEMENS MAGNETOM Trio). Of the 12 IPF subjects, N=8 were imaged a second time ~12 months following the initial scan.

The first point on each dissolved phase radial projection (k0) was phase-shifted such that the mean real to imaginary signal ratio matched the spectroscopically-acquired RBC/Barrier ratio, thus separating RBC and Barrier signals into these channels.3 The RBC k0 signal exhibited oscillations, but the Barrier did not (Figure 1). RBC signal at k0 was then temporally smoothed, detrended, and oscillation amplitudes were normalized to bin projections into high (>0.6×max) and low (<-0.6×max) bins (Figure 1). These binned projections were used to define “keys” for keyhole reconstruction. The key radius was chosen as the largest k-space radius at which k-space was at least 50% sampled. Images were reconstructed using unsmoothed, detrended data for the “high” key, the “low” key, and the complete, non-keyhole dataset.

After reconstruction, images were separated into RBC and barrier components again using the RBC/Barrier ratio. Oscillation amplitude maps were defined voxel by voxel by
(High Key – Low Key)/Non-Keyholed Image×100%
for total dissolved and Barrier images, and by
(High Key – Low Key)/mean(Non-Keyholed Image)×100%
for RBC images (Figure 2) to accommodate the lower SNR in these images. Voxels with SNR<1.5 were removed. RBC oscillation amplitude images were then binned based on the total mean and standard deviation (SD) oscillation for all healthy subjects (Figure 3). The percent of lung voxels binned into “Low” (<healthy mean – 1SD) and “High” (>healthy mean + 1SD) were compared between subjects.

Results

Oscillation amplitude was minimal in the total dissolved phase and barrier images for all subjects. However, RBC images displayed significant oscillation amplitude (Figure 2). Whole lung mean oscillation amplitudes correlated reasonably well with the amplitude measured using whole lung spectroscopy (R2=0.60). Healthy subjects displayed a relatively homogeneous pattern of oscillation, with many subjects exhibiting slightly elevated oscillation amplitude in the posterior portions of the lungs (Figure 3). IPF subjects consistently displayed elevated oscillation amplitude (Figure 4), with the percent of the lungs binned into “high” oscillation significantly higher than healthy (p<0.001) and PAH (p=0.03) subject groups. LHF subjects likewise had a larger percentage binned into the “high” region than healthy subjects (p=0.03). In contrast, PAH and LHF subjects had a greater percent of the lungs binned into “low” oscillation than healthy (PAH: p=0.002, LHF: p=0.03) and IPF (PAH: p=0.001, LHF: p=0.02) (Figure 5). In longitudinal imaging of IPF subjects, neither the percent of lung binned into low nor the high regions changed significantly over time. However, the percent of lung outside of the normal region (%low + %high) increased significantly (p=0.048) (Figure 5).

Discussion and Conclusion

By combining keyhole reconstruction and hyperpolarized 129Xe gas exchange imaging, we are able to image cardiopulmonary blood flow dynamics in the capillary bed of the lungs. Metrics derived from this method correlated well with global spectroscopy, and distinguished between the four subject populations studied based on the percentages of the lungs exhibiting low or high oscillation amplitude as compared to the healthy subject cohort. Furthermore, oscillation amplitude mapping over 1 year in a small cohort of IPF subjects who were imaged twice show similar oscillation patterns and generally increasing percentages of lung showing abnormal oscillations. As this technique involves only post-processing, cardiogenic oscillation imaging can readily be derived from existing data. These oscillation data can thus be included in a standard gas exchange acquisition and analysis pipeline to provide metrics of regional microvasculature function for both existing and future data.

Acknowledgements

This work was supported by the NIH (T32HL007752, R01HL143011, and R01HL126771) and the Cincinnati Children’s Research Foundation.

References

1. Gaskill CF, Carrier EJ, Kropski JA, et al. Disruption of lineage specification in adult pulmonary mesenchymal progenitor cells promotes microvascular dysfunction. The Journal of Clinical Investigation. 2017;127(6):2262-2276.

2. Johns CS, Kiely DG, Swift AJ. Novel imaging techniques in pulmonary hypertension. Current Opinion in Cardiology. 2018;33(6):587-593.

3. Kaushik SS, Robertson SH, Freeman MS, et al. Single-breath clinical imaging of hyperpolarized Xe-129 in the airspaces, barrier, and red blood cells using an interleaved 3D radial 1-Point Dixon acquisition. Magn Reson Med. 2016;75(4):1434-1443.

4. Wang JM, Robertson SH, Wang Z, et al. Using hyperpolarized (129)Xe MRI to quantify regional gas transfer in idiopathic pulmonary fibrosis. Thorax. 2018;73(1):21-28.

5. Wang Z, Bier EA, Swaminathan A, et al. Diverse cardiopulmonary diseases are associated with distinct xenon MRI signatures. European Respiratory Journal. 2019; DOI: 10.1183/13993003.00831-2019.

6. Qing K, Mugler JP, Altes TA, et al. Assessment of lung function in asthma and COPD using hyperpolarized 129Xe chemical shift saturation recovery spectroscopy and dissolved-phase MRI. NMR Biomed. 2014;27(12):1490-1501.

7. Bier EA, Robertson SH, Schrank GM, et al. A protocol for quantifying cardiogenic oscillations in dynamic (129) Xe gas exchange spectroscopy: The effects of idiopathic pulmonary fibrosis. NMR Biomed. 2019;32(1):e4029-e4029.

8. Ruppert K, Altes TA, Mata JF, Ruset IC, Hersman FW, Mugler III JP. Detecting pulmonary capillary blood pulsations using hyperpolarized xenon-129 chemical shift saturation recovery (CSSR) MR spectroscopy. Magn Reson Med. 2016;75(4):1771-1780.

9. Lethmate R, Wajer FTAW, Crémillieux Y, van Ormondt D, Graveron-Demilly D. Dynamic MR-Imaging with radial scanning, a post-acquisition keyhole approach. EURASIP Journal on Advances in Signal Processing. 2003;2003(5):794878.

Figures

Figure 1. k0 signal intensity for (a) Barrier and (b) RBC signal with insets showing Fourier transform of data. (c) Correlation of oscillation frequency and recorded heart rate at the time of scan. (d-f) Workflow for selecting “high” and “low” oscillation bins. Data are smoothed and fit to a biexponential decay (d), detrended (e), oscillations normalized (f), and high and low regions selected with thresholds of ±0.6 times the maximum. (g-i) Keyhole reconstruction: Images are reconstructed corresponding to all data (g), high binned projections (h), and low binned projections (i).

Figure 2. Representative axial and coronal key images of HP 129Xe dissolved in RBCs from a healthy subject and associated oscillation intensity maps. Oscillation amplitude is very near to zero in the barrier images, as expected based on whole lung spectroscopy, which shows minimal oscillation of the barrier signal. In the total dissolved signal images, oscillation amplitude is similarly low. Significant oscillations are observed in the RBC images only, which matches expectation based on whole lung spectroscopy and raw k0 data.

Figure 3. (a) Oscillation amplitude histogram from all lung voxels for all 8 healthy subjects imaged on a 1.5 T scanner. Each color corresponds to one standard deviation away from the healthy subject cohort mean. (b,c) Sample binned oscillation images for healthy subjects imaged on both scanner platforms. Increased oscillation amplitude in the posterior of lungs, possibly indicative of gravitational gradients, is visible in several subjects.

Figure 4. Axial and coronal oscillation amplitude maps for representative healthy subjects and Idiopathic Pulmonary Fibrosis (IPF), Pulmonary Arterial Hypertension (PAH), and Left Heart Failure (LHF) patients. Histograms showing bin counts of percent RBC oscillations are shown next to images with the percent of voxels with oscillation amplitude less than (Healthy Mean – 1SD) (L) and greater than (Healthy Mean + 1SD) (H) listed.

Figure 5. Box plots comparing the percent of the lung binned into low (a) and high (b) oscillation bins among the four subject populations. PAH and LHF subjects had a greater percentage of the lung with low oscillations compared to healthy while IPF and LHF subjects had a greater percentage with high oscillations. (c) Slope plot showing an increase in lung percentage with abnormal (> 1SD from healthy cohort mean) oscillation in IPF subjects after ~12 months of disease progression. (d) Representative binned oscillation images for IPF subjects showing disease progression over ~12 months.

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
0432