Hyperpolarized 129Xenon MRI Ventilation Defect Quantification via Thresholding and Linear Binning in Multiple Pulmonary Diseases

129Xe Image Processing Pipeline: An open-source, graphical user interface application for the analysis of hyperpolarized 129Xe MRI

PURPOSE

Hyperpolarized 129Xe MRI presents opportunities to assess regional pulmonary microstructure and function. Ongoing advancements in hardware, sequences, and image processing have helped it become increasingly adopted for both research and clinical use. As the number of applications and users increase, standardization becomes crucial. To that end, this study developed an executable, open-source 129Xe image processing pipeline (XIPline) to provide a user-friendly, graphical user interface-based analysis pipeline to analyze and visualize 129Xe MR data, including scanner calibration, ventilation, diffusion-weighted, and gas exchange images.

METHODS

The customizable XIPline is designed in MATLAB to analyze data from all three major scanner platforms. Calibration data is processed to calculate optimal flip angle and determine129Xe frequency offset. Data processing includes loading, reconstructing, registering, segmenting, and post-processing images. Ventilation analysis incorporates three common algorithms to calculate ventilation defect percentage and novel techniques to assess defect distribution and ventilation texture. Diffusion analysis features ADC mapping, modified linear binning to account for ADC age-dependence, and common diffusion morphometry methods. Gas exchange processing uses a generalized linear binning for data acquired using 1-point Dixon imaging.

RESULTS

The XIPline workflow is demonstrated using analysis from representative calibration, ventilation, diffusion, and gas exchange data.

CONCLUSION

The application will reduce redundant effort when implementing new techniques across research sites by providing an open-source framework for developers. In its current form, it offers a robust and adaptable platform for 129Xe MRI analysis to ensure methodological consistency, transparency, and support for collaborative research across multiple sites and MRI manufacturers.

Quantifying Spatial Distribution of Ventilation Defects in Multiple Pulmonary Diseases With Hyperpolarized 129Xenon MRI

BACKGROUND

Hyperpolarized 129Xe MRI assesses lung ventilation, often using the ventilation defect percentage (VDP). Unlike VDP, defect distribution index (DDI) quantifies spatial clustering of defects.

PURPOSE

To quantify spatial distribution of 129Xe ventilation defects using DDI across pulmonary diseases.

STUDY TYPE

Retrospective.

SUBJECTS

Four hundred twenty-one subjects (age = 23.1 ± 17.1, female = 230), comprising healthy controls (N = 60) and subjects with obstructive conditions (asthma [N = 25], bronchiolitis obliterans syndrome [BOS, N = 18], cystic fibrosis [CF, N = 90], lymphangioleiomyomatosis [LAM, N = 50]), restrictive conditions (bleomycin-treated cancer survivors [BLEO, N = 14]; fibrotic lung diseases [FLD, N = 92]), bone marrow transplantation (BMT, N = 53), and bronchopulmonary dysplasia (BPD, N = 19).

FIELD STRENGTH/SEQUENCE

3 T, two-dimensional multi-slice gradient echo.

ASSESSMENT

Whole-lung mean DDI was extracted from DDI maps; correlated with VDP (percent of pixels <60% of whole-lung mean signal intensity) and pulmonary function tests (PFTs) including FEV1, FVC, and FEV1/FVC. DDI and DDI/VDP, a marker of defect clustering, were compared across diseases.

STATISTICAL TESTS

Pearson correlation analysis and Kruskal-Wallis tests. P < 0.0056 for disease groups, P < 0.0125 for categories.

RESULTS

DDI was significantly elevated in BMT (8.3 ± 11.5), BOS (30.1 ± 57.5), BPD (16.0 ± 46.8), CF (15.4 ± 27.2), and LAM (12.6 ± 34.2) compared to controls (1.8 ± 3.1). DDI correlated significantly with VDP in all groups (r ≥ 0.56) except BLEO, and with PFTs in CF, FLD, and LAM (r ≥ 0.56). Obstructive groups had significantly higher mean DDI (14.0 ± 32.0) than controls (1.8 ± 3.0) and restrictive groups (4.0 ± 12.0). DDI/VDP was significantly lower in the restrictive group (0.6 ± 0.6) than controls (0.8 ± 0.6) and obstructive group (1.0 ± 1.0).

DATA CONCLUSION

DDI may provide insights into the distribution of ventilation defects across diseases.

EVIDENCE LEVEL

3 TECHNICAL EFFICACY: Stage 2.

Same-Day Repeatability and 28-Day Reproducibility of Xenon MRI Ventilation in Children With Cystic Fibrosis in a Multi-Site Trial

BACKGROUND

MRI with xenon-129 gas (Xe MRI) can assess airflow obstruction and heterogeneity in lung diseases. Specifically, Xe MRI may represent a sensitive modality for future therapeutic trials of cystic fibrosis (CF) therapies. The reproducibility of Xe MRI has not yet been assessed in the context of a multi-site study.

PURPOSE

To determine the same-day repeatability and 28-day reproducibility of Xe MRI in children with CF.

STUDY TYPE

Four-center prospective, longitudinal.

POPULATION

Thirty-eight children (18 females, 47%), median interquartile range (IQR) age 12 (9-14) years old, with mild CF (forced expiratory volume in 1 second (FEV1) ≥85% predicted).

FIELD STRENGTH/SEQUENCE

3-T, two-dimensional (2D) gradient-echo (GRE) sequence.

ASSESSMENT

Xe MRI, FEV1, and nitrogen multiple-breath wash-out for lung-clearance index (LCI2.5) were performed. To assess same-day reproducibility, Xe MRI was performed twice within the first visit, and procedures were repeated at 28 days. Xe hypoventilation was quantified using ventilation-defect percentage (VDP) and reader-defect volume (RDV). For VDP, hypoventilated voxels from segmented images were identified using a threshold of <60% mean whole-lung signal and expressed as a percentage of the lung volume. For RDV, hypoventilation was identified by two trained readers and expressed as a percentage.

STATISTICAL TESTS

Inter-site comparisons were conducted using Kruskal-Wallis nonparametric tests with Dunn's multiple-comparisons tests. Differences for individuals were assessed using Wilcoxon matched-pairs tests. Bland-Altman tests were used to evaluate same-day repeatability, 28-day reproducibility, and inter-reader agreement. A P-value ≤0.05 was considered significant.

RESULTS

Median FEV1 %-predicted was 96.8% (86%-106%), and median LCI2.5 was 6.6 (6.3-7.4). Xe MRI had high same-day reproducibility (mean VDP difference 0.12%, 95% limits of agreement [-3.2, 3.4]; mean RDV difference 0.42% [-2.5, 3.3]). At 28 days, 26/31 participants (84%) fell within the same-day 95% limits of agreement.

DATA CONCLUSION

Xe MRI may offer excellent same-day and short-term reproducibility.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 2.

Analytical corrections for B1-inhomogeneity and signal decay in multi-slice 2D spiral hyperpolarized 129Xe MRI using keyhole reconstruction

PURPOSE

Hyperpolarized xenon MRI suffers from heterogeneous coil bias and magnetization decay that obscure pulmonary abnormalities. Non-physiological signal variability can be mitigated by measuring and mapping the nominal flip angle, and by rescaling the images to correct for signal bias and decay. While flip angle maps can be calculated from sequentially acquired images, scan time and breath-hold duration are doubled. Here, we exploit the low-frequency oversampling of 2D-spiral and keyhole reconstruction to measure flip angle maps from a single acquisition.

METHODS

Flip angle maps were calculated from two images generated from a single dataset using keyhole reconstructions and a Bloch-equation-based model suitable for hyperpolarized substances. Artifacts resulting from acquisition and reconstruction schemes (e.g., keyhole reconstruction radius, slice-selection profile, spiral-ordering, and oversampling) were assessed using point-spread functions. Simulated flip angle maps generated using keyhole reconstruction were compared against the paired-image approach using RMS error (RMSE). Finally, feasibility was demonstrated for in vivo xenon ventilation imaging.

RESULTS

Simulations demonstrated accurate flip angle maps and B1-inhomogeneity correction can be generated with only 1.25-fold central-oversampling and keyhole reconstruction radius = 5% (RMSE = 0.460°). These settings also generated accurate flip angle maps in a healthy control (RSME = 0.337°) and a person with cystic fibrosis (RMSE = 0.404°) in as little as 3.3 s.

CONCLUSION

Regional lung ventilation images with reduced impact of B1-inhomogeneity can be acquired rapidly by combining 2D-spiral acquisition, Bloch-equation-based modeling, and keyhole reconstruction. This approach will be especially useful for breath-hold studies where short scan durations are necessary, such as dynamic imaging and applications in children or people with severely compromised respiratory function.

Innovations in Childhood Interstitial and Diffuse Lung Disease

Children's interstitial and diffuse lung diseases (chILDs) are a heterogenous and diverse group of lung disorders presenting during childhood. Infants and children with chILD disorders present with respiratory signs and symptoms as well as diffuse lung imaging abnormalities. ChILD disorders are associated with significant health care resource utilization and high morbidity and mortality. The care of patients with chILD has been improved through multidisciplinary care, multicenter collaboration, and the establishment of patient research networks in the United Stated and abroad. This review details past and current innovations in the diagnosis and clinical care of children with chILD.

Hyperpolarized 129Xe MRI, 99mTc scintigraphy, and SPECT in lung ventilation imaging: a quantitative comparison

RATIONALE AND OBJECTIVES

The current clinical standard for functional imaging of patients with lung ailments is nuclear medicine scintigraphy and Single Photon Emission Computed Tomography (SPECT) which detect the gamma decay of inhaled radioactive tracers. Hyperpolarized (HP) Xenon-129 MRI (XeMRI) of the lungs has recently been FDA approved and provides similar functional images of the lungs with higher spatial resolution than scintigraphy and SPECT. Here we compare Technetium-99m (99mTc) diethylene-triamine-pentaacetate scintigraphy and SPECT with HP XeMRI in healthy controls, asthma, and chronic obstructive pulmonary disorder (COPD) patients.

MATERIALS AND METHODS

59 subjects, healthy, with asthma, and with COPD, underwent 99mTc scintigraphy/SPECT, standard spirometry, and HP XeMRI. XeMRI and SPECT images were registered for direct voxel-wise signal comparisons. Images were also compared using ventilation defect percentage (VDP), and a standard 6-compartment method. VDP calculated from XeMRI and SPECT images was compared to spirometry.

RESULTS

Median Pearson correlation coefficient for voxel-wise signal comparison was 0.698 (0.613-0.782) between scintigraphy and XeMRI and 0.398 (0.286-0.502) between SPECT and XeMRI. Correlation between VDP measures was r = 0.853, p < 0.05. VDP separated asthma and COPD from the control group and was significantly correlated with FEV1, FEV1/FVC, and FEF 25-75.

CONCLUSION

HP XeMRI provides equivalent information to 99mTc SPECT and standard spirometry measures. Additionally, XeMRI is non-invasive, hence it could be used for longitudinal studies for evaluating emerging treatment for lung ailments.

Simultaneous quantification of hyperpolarized xenon-129 ventilation and gas exchange with multi-breath xenon-polarization transfer contrast (XTC) MRI

PURPOSE

To demonstrate the feasibility of a multi-breath xenon-polarization transfer contrast (XTC) MR imaging approach for simultaneously evaluating regional ventilation and gas exchange parameters.

METHODS

Imaging was performed in five healthy volunteers and six chronic obstructive pulmonary disease (COPD) patients. The multi-breath XTC protocol consisted of three repeated schemes of six wash-in breaths of a xenon mixture and four normoxic wash-out breaths, with and without selective saturation of either the tissue membrane or red blood cell (RBC) resonances. Acquisitions were performed at end-exhalation while subjects maintained tidal breathing throughout the session. The no-saturation, membrane-saturation, and RBC-saturation images were fit to a per-breath gas replacement model for extracting voxelwise tidal volume (TV), functional residual capacity (FRC), and fractional ventilation (FV), as well as tissue- and RBC-gas exchange (fMem and fRBC , respectively). The sensitivity of the derived model was also evaluated via simulations.

RESULTS

With the exception of FRC, whole-lung averages for all metrics were decreased in the COPD subjects compared to the healthy cohort, significantly so for FV, fRBC , and fMem . Heterogeneity was higher overall in the COPD subjects, particularly for fRBC , fMem , and fRBC:Mem . The anterior-to-posterior gradient associated with the gravity-dependence of lung function in supine imaging was also evident for FV, fRBC , and fMem values in the healthy subjects, but noticeably absent in the COPD cohort.

CONCLUSION

Multi-breath XTC imaging generated high-resolution, co-registered maps of ventilation and gas exchange parameters acquired during tidal breathing and with low per-breath xenon doses. Clear differences between healthy and COPD subjects were apparent and consistent with spirometry.

Short-term structural and functional changes after airway clearance therapy in cystic fibrosis

BACKGROUND

Airway clearance therapy (ACT) with a high-frequency chest wall oscillation (HFCWO) vest is a common but time-consuming treatment. Its benefit to quality of life for cystic fibrosis (CF) patients is well established but has been questioned recently as new highly-effective modulator therapies begin to change the treatment landscape. 129Xe ventilation MRI has been shown to be very sensitive to lung obstruction in mild CF disease, making it an ideal tool to identify and quantify subtle, regional changes.

METHODS

20 CF patients (ages 20.7 ± 5.1 years) refrained from performing ACT before arriving for a single-day visit. Multiple-breath washout (MBW), spirometry, Xe MRI, and ultrashort echo-time (UTE) MRI were obtained twice-before and after patients performed ACT using their prescribed HFCWO vests (average 4.7 ± 0.5 h). UTE MRIs were scored for structural abnormalities, and standard functional metrics were obtained from MBW, spirometry, and Xe MRI-FEV1,pp, LCI2.5, and VDPN4, respectively.

RESULTS

Spirometry and Xe MRI detected significant improvements in lung function post-ACT. 15/20 patients showed improvements from a baseline median of 92% FEV1,pp. Similarly, 16/20 patients showed improvements in Xe MRI from a baseline median of 15.2% VDPN4. Average individual changes were +2.6% in FEV1,pp and -1.3% in VDPN4, but without spatial correlations to easily-identifiable causative structural defects (e.g. mucus plugs or bronchiectasis) on UTE MRI.

CONCLUSIONS

Lung function improved after a single instance of HFCWO-vest ACT and was detectable by spirometry and Xe MRI. The only common structural abnormalities were mucus plugs, which corresponded to ventilation defects, but ventilation defects were often present without visible abnormalities.

A single-breath-hold protocol for hyperpolarized 129 Xe ventilation and gas exchange imaging

Hyperpolarized 129 Xe MRI (Xe-MRI) is increasingly used to image the structure and function of the lungs. Because 129 Xe imaging can provide multiple contrasts (ventilation, alveolar airspace size, and gas exchange), imaging often occurs over several breath-holds, which increases the time, expense, and patient burden of scans. We propose an imaging sequence that can be used to acquire Xe-MRI gas exchange and high-quality ventilation images within a single, approximately 10 s, breath-hold. This method uses a radial one-point Dixon approach to sample dissolved 129 Xe signal, which is interleaved with a 3D spiral ("FLORET") encoding pattern for gaseous 129 Xe. Thus, ventilation images are obtained at higher nominal spatial resolution (4.2 × 4.2 × 4.2 mm3 ) compared with gas-exchange images (6.25 × 6.25 × 6.25 mm3 ), both competitive with current standards within the Xe-MRI field. Moreover, the short 10 s Xe-MRI acquisition time allows for 1 H "anatomic" images used for thoracic cavity masking to be acquired within the same breath-hold for a total scan time of about 14 s. Images were acquired using this single-breath method in 11 volunteers (N = 4 healthy, N = 7 post-acute COVID). For 11 of these participants, a separate breath-hold was used to acquire a "dedicated" ventilation scan and five had an additional "dedicated" gas exchange scan. The images acquired using the single-breath protocol were compared with those from dedicated scans using Bland-Altman analysis, intraclass correlation (ICC), structural similarity, peak signal-to-noise ratio, Dice coefficients, and average distance. Imaging markers from the single-breath protocol showed high correlation with dedicated scans (ventilation defect percent, ICC = 0.77, p = 0.01; membrane/gas, ICC = 0.97, p = 0.001; red blood cell/gas, ICC = 0.99, p < 0.001). Images showed good qualitative and quantitative regional agreement. This single-breath protocol enables the collection of essential Xe-MRI information within one breath-hold, simplifying scanning sessions and reducing costs associated with Xe-MRI.

Childhood to adulthood: Accounting for age dependence in healthy-reference distributions in 129 Xe gas-exchange MRI

PURPOSE

Xenon-129 (129 Xe) gas-exchange MRI is a pulmonary-imaging technique that provides quantitative metrics for lung structure and function and is often compared to pulmonary-function tests. Unlike such tests, it does not normalize to predictive values based on demographic variables such as age. Many sites have alluded to an age dependence in gas-exchange metrics; however, a procedure for normalizing metrics has not yet been introduced.

THEORY

We model healthy reference values for 129 Xe gas-exchange MRI against age using generalized additive models for location, scale, and shape (GAMLSS). GAMLSS takes signal data from an aggregated heathy-reference cohort and fits a distribution with flexible median, variation, skewness, and kurtosis to predict age-dependent centiles. This approach mirrors methods by the Global Lung Function Initiative for modeling pulmonary-function test data and applies it to binning methods widely used by the 129 Xe MRI community to interpret and quantify gas-exchange data.

METHODS

Ventilation, membrane-uptake, red blood cell transfer, and red blood cell:membrane gas-exchange metrics were collected on 30 healthy subjects over an age range of 5 to 68 years. A GAMLSS model was fit against age and compared against widely used linear and generalized-linear binning 129 Xe MRI analysis schemes.

RESULTS

All 4 gas-exchange metrics had significant skewness, and membrane-uptake had significant kurtosis compared to a normal distribution. Age has significant impact on distribution parameters. GAMLSS-binning produced narrower bins compared to the linear and generalized-linear binning schemes and distributed signal data closer to a normal distribution.

CONCLUSION

The proposed "proof-of-concept" GAMLSS-binning approach can improve diagnostic accuracy of 129 Xe gas-exchange MRI by providing a means of modeling voxel distribution data against age.

Bias field correction in hyperpolarized 129 Xe ventilation MRI using templates derived by RF-depolarization mapping

PURPOSE

To correct for RF inhomogeneity for in vivo 129 Xe ventilation MRI using flip-angle mapping enabled by randomized 3D radial acquisitions. To extend this RF-depolarization mapping approach to create a flip-angle map template applicable to arbitrary acquisition strategies, and to compare these approaches to conventional bias field correction.

METHODS

RF-depolarization mapping was evaluated first in digital simulations and then in 51 subjects who had undergone radial 129 Xe ventilation MRI in the supine position at 3T (views = 3600; samples/view = 128; TR/TE = 4.5/0.45 ms; flip angle = 1.5; FOV = 40 cm). The images were corrected using newly developed RF-depolarization and templated-based methods and the resulting quantitative ventilation metrics (mean, coefficient of variation, and gradient) were compared to those resulting from N4ITK correction.

RESULTS

RF-depolarization and template-based mapping methods yielded a pattern of RF-inhomogeneity consistent with the expected variation based on coil architecture. The resulting corrected images were visually similar, but meaningfully distinct from those generated using standard N4ITK correction. The N4ITK algorithm eliminated the physiologically expected anterior-posterior gradient (-0.04 ± 1.56%/cm, P < 0.001). These 2 newly introduced methods of RF-depolarization and template correction retained the physiologically expected anterior-posterior ventilation gradient in healthy subjects (2.77 ± 2.09%/cm and 2.01 ± 2.73%/cm, respectively).

CONCLUSIONS

Randomized 3D 129 Xe MRI ventilation acquisitions can inherently be corrected for bias field, and this technique can be extended to create flip angle templates capable of correcting images from a given coil regardless of acquisition strategy. These methods may be more favorable than the de facto standard N4ITK because they can remove undesirable heterogeneity caused by RF effects while retaining results from known physiology.