In vivo lung morphometry with hyperpolarized 3He diffusion MRI: theoretical background

Age-Trajectories of Higher-Order Diffusion Properties of Major Brain Metabolites in Cerebral and Cerebellar Gray Matter Using In Vivo Diffusion-Weighted MR Spectroscopy at 3T

Healthy brain aging involves changes in both brain structure and function, including alterations in cellular composition and microstructure across brain regions. Unlike diffusion-weighted MRI (dMRI), diffusion-weighted MR spectroscopy (dMRS) can assess cell-type specific microstructural changes, providing indirect information on both cell composition and microstructure through the quantification and interpretation of metabolites' diffusion properties. This work investigates age-related changes in the higher-order diffusion properties of total N-Acetyl-aspartate (neuronal biomarker), total choline (glial biomarker), and total creatine (both neuronal and glial biomarker) beyond the classical apparent diffusion coefficient in cerebral and cerebellar gray matter of healthy human brain. Twenty-five subjects were recruited and scanned using a diffusion-weighted semi-LASER sequence in two brain regions-of-interest (ROI) at 3T: posterior-cingulate (PCC) and cerebellar cortices. Metabolites' diffusion was characterized by quantifying metrics from both Gaussian and non-Gaussian signal representations and biophysical models. All studied metabolites exhibited lower apparent diffusivities and higher apparent kurtosis values in the cerebellum compared to the PCC, likely stemming from the higher microstructural complexity of cellular composition in the cerebellum. Multivariate regression analysis (accounting for ROI tissue composition as a covariate) showed slight decrease (or no change) of all metabolites' diffusivities and slight increase of all metabolites' kurtosis with age, none of which statistically significant (p > 0.05). The proposed age-trajectories provide benchmarks for identifying anomalies in the diffusion properties of major brain metabolites which could be related to pathological mechanisms altering both the brain microstructure and cellular composition.

Machine Learning and CT Texture Features in Ex-smokers with no CT Evidence of Emphysema and Mildly Abnormal Diffusing Capacity

RATIONALE AND OBJECTIVES

Ex-smokers without spirometry or CT evidence of chronic obstructive pulmonary disease (COPD) but with mildly abnormal diffusing capacity of the lungs for carbon monoxide (DLCO) are at higher risk of developing COPD. It remains difficult to make clinical management decisions for such ex-smokers without other objective assessments consistent with COPD. Hence, our objective was to develop a machine-learning and CT texture-analysis pipeline to dichotomize ex-smokers with normal and abnormal DLCO (DLCO≥75%pred and DLCO<75%pred).

MATERIALS AND METHODS

In this retrospective study, 71 ex-smokers (50-85yrs) without COPD underwent spirometry, plethysmography, thoracic CT, and 3He MRI to generate ventilation defect percent (VDP) and apparent diffusion coefficients (ADC). PyRadiomics was utilized to extract 496 CT texture-features; Boruta and principal component analysis were used for feature selection and various models were investigated for classification. Machine-learning classifiers were evaluated using area under the receiver operator characteristic curve (AUC), sensitivity, specificity, and F1-measure.

RESULTS

Of 71 ex-smokers without COPD, 29 with mildly abnormal DLCO had significantly different MRI ADC (p < .001), residual-volume to total-lung-capacity ratio (p = .003), St. George's Respiratory Questionnaire (p = .029), and six-minute-walk distance (6MWD) (p < .001), but similar relative area of the lung < -950 Hounsfield-units (RA950) (p = .9) compared to 42 ex-smokers with normal DLCO. Logistic-regression machine-learning mixed-model trained on selected texture-features achieved the best classification accuracy of 87%. All clinical and imaging measurements were outperformed by high-high-pass filter high-gray-level-run-emphasis texture-feature (AUC=0.81), which correlated with DLCO (ρ = -0.29, p = .02), MRI ADC (ρ = 0.23, p = .048), and 6MWD (ρ = -0.25, p = .02).

CONCLUSION

In ex-smokers with no CT evidence of emphysema, machine-learning models exclusively trained on CT texture-features accurately classified ex-smokers with abnormal diffusing capacity, outperforming conventional quantitative CT measurements.

Chest MRI and CT Predictors of 10-Year All-Cause Mortality in COPD

Pulmonary imaging measurements using magnetic resonance imaging (MRI) and computed tomography (CT) have the potential to deepen our understanding of chronic obstructive pulmonary disease (COPD) by measuring airway and parenchymal pathologic information that cannot be provided by spirometry. Currently, MRI and CT measurements are not included in mortality risk predictions, diagnosis, or COPD staging. We evaluated baseline pulmonary function, MRI and CT measurements alongside imaging texture-features to predict 10-year all-cause mortality in ex-smokers with (n = 93; 31 females; 70 ± 9years) and without (n = 69; 29 females, 69 ± 9years) COPD. CT airway and vessel measurements, helium-3 (3He) MRI ventilation defect percent (VDP) and apparent diffusion coefficients (ADC) were quantified. MRI and CT texture-features were extracted using PyRadiomics (version2.2.0). Associations between 10-year all-cause mortality and all clinical and imaging measurements were evaluated using multivariable regression model odds-ratios. Machine-learning predictive models for 10-year all-cause mortality were evaluated using area-under-receiver-operator-characteristic-curve (AUC), sensitivity and specificity analyses. DLCO (%pred) (HR = 0.955, 95%CI: 0.934-0.976, p < 0.001), MRI ADC (HR = 1.843, 95%CI: 1.260-2.871, p < 0.001), and CT informational-measure-of-correlation (HR = 3.546, 95% CI: 1.660-7.573, p = 0.001) were the strongest predictors of 10-year mortality. A machine-learning model trained on clinical, imaging, and imaging textures was the best predictive model (AUC = 0.82, sensitivity = 83%, specificity = 84%) and outperformed the solely clinical model (AUC = 0.76, sensitivity = 77%, specificity = 79%). In ex-smokers, regardless of COPD status, addition of CT and MR imaging texture measurements to clinical models provided unique prognostic information of mortality risk that can allow for better clinical management.Clinical Trial Registration: www.clinicaltrials.gov NCT02279329.

Preclinical MRI Using Hyperpolarized 129Xe

Although critical for development of novel therapies, understanding altered lung function in disease models is challenging because the transport and diffusion of gases over short distances, on which proper function relies, is not readily visualized. In this review we summarize progress introducing hyperpolarized ¹²⁹Xe imaging as a method to follow these processes in vivo. The work is organized in sections highlighting methods to observe the gas replacement effects of breathing (Gas Dynamics during the Breathing Cycle) and gas diffusion throughout the parenchymal airspaces (3). We then describe the spectral signatures indicative of gas dissolution and uptake (4), and how these features can be used to follow the gas as it enters the tissue and capillary bed, is taken up by hemoglobin in the red blood cells (5), re-enters the gas phase prior to exhalation (6), or is carried via the vasculature to other organs and body structures (7). We conclude with a discussion of practical imaging and spectroscopy techniques that deliver quantifiable metrics despite the small size, rapid motion and decay of signal and coherence characteristic of the magnetically inhomogeneous lung in preclinical models (8).

Mapping Alveolar Oxygen Partial Pressure in COPD Using Hyperpolarized Helium-3: The Multi-Ethnic Study of Atherosclerosis (MESA) COPD Study

Chronic obstructive pulmonary disease (COPD) and emphysema are characterized by functional and structural damage which increases the spaces for gaseous diffusion and impairs oxygen exchange. Here we explore the potential for hyperpolarized (HP) 3He MRI to characterize lung structure and function in a large-scale population-based study. Participants (n = 54) from the Multi-Ethnic Study of Atherosclerosis (MESA) COPD Study, a nested case-control study of COPD among participants with 10+ packyears underwent HP 3He MRI measuring pAO2, apparent diffusion coefficient (ADC), and ventilation. HP MRI measures were compared to full-lung CT and pulmonary function testing. High ADC values (>0.4 cm2/s) correlated with emphysema and heterogeneity in pAO2 measurements. Strong correlations were found between the heterogeneity of global pAO2 as summarized by its standard deviation (SD) (p < 0.0002) and non-physiologic pAO2 values (p < 0.0001) with percent emphysema on CT. A regional study revealed a strong association between pAO2 SD and visual emphysema severity (p < 0.003) and an association with the paraseptal emphysema subtype (p < 0.04) after adjustment for demographics and smoking status. HP noble gas pAO2 heterogeneity and the fraction of non-physiological pAO2 results increase in mild to moderate COPD. Measurements of pAO2 are sensitive to regional emphysematous damage detected by CT and may be used to probe pulmonary emphysema subtypes. HP noble gas lung MRI provides non-invasive information about COPD severity and lung function without ionizing radiation.

Rat Brain Global Ischemia-Induced Diffusion Changes Revisited: Biophysical Modeling of the Water and NAA MR "Diffusion Signal"

PURPOSE

To assess changes in intracellular diffusion as a mechanism for the reduction in water ADC that accompanies brain injury. Using NAA as a marker of neuronal cytoplasmic diffusion, NAA diffusion was measured before and after global ischemia (immediately postmortem) in the female Sprague-Dawley rat.

METHODS

Diffusion-weighted PRESS spectra, with diffusion encoding in a single direction, were acquired from large voxels of rat brain gray matter in vivo and postischemia employing either pairs of pulsed half-sine-shaped gradients (in vivo and postischemia, b_max  = 19 ms/μm² ) or sinusoidal oscillating gradients (in vivo only) with frequencies of 99.2-250 Hz. A 2D randomly oriented cylinder (neurite) model gave estimates of longitudinal and transverse diffusivities (D_L and D_T , respectively). In this model, D_L represents the "free" diffusivity of NAA, whereas D_T reflects highly restricted diffusion. Using oscillating gradients, the frequency dependence of D_T [D_T (ω)] gave estimates of the cylinder (axon/dendrite) radius.

RESULTS

A 10% decrease in D_L,NAA followed global ischemia, dropping from 0.391 ± 0.012 μm² /ms to 0.350 ± 0.009 μm² /ms. Modeling D_T,NAA (ω) provided an estimate of the neurite radius of 1.0 ± 0.6 μm.

CONCLUSION

Whereas the increase in apparent intraneuronal viscosity suggested by changes in D_L,NAA may contribute to the overall reduction in water ADC associated with brain injury, it is not sufficient to be the sole explanation. Estimates of neurite radius based on D_T (ω) were consistent with literature values.

Improving hyperpolarized 129 Xe ADC mapping in pediatric and adult lungs with uncertainty propagation

RATIONALE

Hyperpolarized (HP) ¹²⁹ Xe-MRI provides non-invasive methods to quantify lung function and structure, with the ¹²⁹ Xe apparent diffusion coefficient (ADC) being a well validated measure of alveolar airspace size. However, the experimental factors that impact the precision and accuracy of HP ¹²⁹ Xe ADC measurements have not been rigorously investigated. Here, we introduce an analytical model to predict the experimental uncertainty of ¹²⁹ Xe ADC estimates. Additionally, we report ADC dependence on age in healthy pediatric volunteers.

METHODS

An analytical expression for ADC uncertainty was derived from the Stejskal-Tanner equation and simplified Bloch equations appropriate for HP media. Parameters in the model were maximum b-value (b_max ), number of b-values (N_b ), number of phase encoding lines (N_ph ), flip angle and the ADC itself. This model was validated by simulations and phantom experiments, and five fitting methods for calculating ADC were investigated. To examine the lower range for ¹²⁹ Xe ADC, 32 healthy subjects (age 6-40 years) underwent diffusion-weighted ¹²⁹ Xe MRI.

RESULTS

The analytical model provides a lower bound on ADC uncertainty and predicts that decreased signal-to-noise ratio yields increases in relative uncertainty (ϵADC) . As such, experimental parameters that impact non-equilibrium ¹²⁹ Xe magnetization necessarily impact the resulting ϵADC . The values of diffusion encoding parameters (N_b and b_max ) that minimize ϵADC strongly depend on the underlying ADC value, resulting in a global minimum for ϵADC . Bayesian fitting outperformed other methods (error < 5%) for estimating ADC. The whole-lung mean ¹²⁹ Xe ADC of healthy subjects increased with age at a rate of 1.75 × 10^-4  cm² /s/yr (p = 0.001).

CONCLUSIONS

HP ¹²⁹ Xe diffusion MRI can be improved by minimizing the uncertainty of ADC measurements via uncertainty propagation. Doing so will improve experimental accuracy when measuring lung microstructure in vivo and should allow improved monitoring of regional disease progression and assessment of therapy response in a range of lung diseases.

Accelerate gas diffusion-weighted MRI for lung morphometry with deep learning

OBJECTIVES

Multiple b-value gas diffusion-weighted MRI (DW-MRI) enables non-invasive and quantitative assessment of lung morphometry, but its long acquisition time is not well-tolerated by patients. We aimed to accelerate multiple b-value gas DW-MRI for lung morphometry using deep learning.

METHODS

A deep cascade of residual dense network (DC-RDN) was developed to reconstruct high-quality DW images from highly undersampled k-space data. Hyperpolarized ¹²⁹Xe lung ventilation images were acquired from 101 participants and were retrospectively collected to generate synthetic DW-MRI data to train the DC-RDN. Afterwards, the performance of the DC-RDN was evaluated on retrospectively and prospectively undersampled multiple b-value ¹²⁹Xe MRI datasets.

RESULTS

Each slice with size of 64 × 64 × 5 could be reconstructed within 7.2 ms. For the retrospective test data, the DC-RDN showed significant improvement on all quantitative metrics compared with the conventional reconstruction methods (p < 0.05). The apparent diffusion coefficient (ADC) and morphometry parameters were not significantly different between the fully sampled and DC-RDN reconstructed images (p > 0.05). For the prospectively accelerated acquisition, the required breath-holding time was reduced from 17.8 to 4.7 s with an acceleration factor of 4. Meanwhile, the prospectively reconstructed results showed good agreement with the fully sampled images, with a mean difference of -0.72% and -0.74% regarding global mean ADC and mean linear intercept (L_m) values.

CONCLUSIONS

DC-RDN is effective in accelerating multiple b-value gas DW-MRI while maintaining accurate estimation of lung microstructural morphometry, facilitating the clinical potential of studying lung diseases with hyperpolarized DW-MRI.

KEY POINTS

• The deep cascade of residual dense network allowed fast and high-quality reconstruction of multiple b-value gas diffusion-weighted MRI at an acceleration factor of 4. • The apparent diffusion coefficient and morphometry parameters were not significantly different between the fully sampled images and the reconstructed results (p > 0.05). • The required breath-holding time was reduced from 17.8 to 4.7 s and each slice with size of 64 × 64 × 5 could be reconstructed within 7.2 ms.

The relationship between diffusion heterogeneity and microstructural changes in high-grade gliomas using Monte Carlo simulations

PURPOSE

Diffusion-weighted imaging (DWI) may aid accurate tumor grading. Decreased diffusivity and increased diffusion heterogeneity measures have been observed in high-grade gliomas using the non-monoexponential models for DWI. However, DWI measures concerning tissue characteristics in terms of pathophysiological and structural changes are yet to be established. Thus, this study aims to investigate the relationship between the diffusion measurements and microstructural changes in the presence of high-grade gliomas using a three-dimensional Monte Carlo simulation with systematic changes of microstructural parameters.

METHODS

Water diffusion was simulated in a microenvironment along with changes associated with the presence of high-grade gliomas, including increases in cell density, nuclear volume, extracellular volume (VF_ex), and extracellular tortuosity (λ_ex), and changes in membrane permeability (P_mem). DWI signals were simulated using a pulsed gradient spin-echo sequence. The sequence parameters, including the maximum gradient strength and diffusion time, were set to be comparable to those of clinical scanners and advanced human MRI systems. The DWI signals were fitted using the gamma distribution and diffusional kurtosis models with b-values up to 6000 and 2500 s/mm², respectively.

RESULTS

The diffusivity measures (apparent diffusion coefficients (ADC), D_gamma of the gamma distribution model and D_app of the diffusional kurtosis model) decreased with increases in cell density and λ_ex, and a decrease in P_mem. These diffusivity measures increased with increases in nuclear volume and VF_ex. The diffusion heterogeneity measures (σ_gamma of the gamma distribution model and K_app of the diffusional kurtosis model) increased with increases in cell density or nuclear volume at the low P_mem, and a decrease in P_mem. Increased σ_gamma was also associated with an increase in VF_ex.

CONCLUSION

Among simulated microstructural changes, only increases in cell density at low P_mem or decreases in P_mem corresponded to both the decreased diffusivity and increased diffusion heterogeneity measures. The results suggest that increases in cell density at low P_mem or decreases in P_mem may be associated with the diffusion changes observed in high-grade gliomas.

Validating in vivo hyperpolarized 129 Xe diffusion MRI and diffusion morphometry in the mouse lung

PURPOSE

Diffusion and lung morphometry imaging using hyperpolarized gases are promising tools to quantify pulmonary microstructure noninvasively in humans and in animal models. These techniques assume the motion encoded is exclusively diffusive gas displacement, but the impact of cardiac motion on measurements has never been explored. Furthermore, although diffusion morphometry has been validated against histology in humans and mice using 3 He, it has never been validated in mice for 129 Xe. Here, we examine the effect of cardiac motion on diffusion imaging and validate 129 Xe diffusion morphometry in mice.

THEORY AND METHODS

Mice were imaged using gradient-echo-based diffusion imaging, and apparent diffusion-coefficient (ADC) maps were generated with and without cardiac gating. Diffusion-weighted images were fit to a previously developed theoretical model using Bayesian probability theory, producing morphometric parameters that were compared with conventional histology.

RESULTS

Cardiac gating had no significant impact on ADC measurements (dual-gating: ADC = 0.020 cm2 /s, single-gating: ADC = 0.020 cm2 /s; P = .38). Diffusion-morphometry-generated maps of ADC (mean, 0.0165 ± 0.0001 cm2 /s) and acinar dimensions (alveolar sleeve depth [h] = 44 µm, acinar duct radii [R] = 99 µm, mean linear intercept [Lm ] = 74 µm) that agreed well with conventional histology (h = 45 µm, R = 108 µm, Lm = 63 µm).

CONCLUSION

Cardiac motion has negligible impact on 129 Xe ADC measurements in mice, arguing its impact will be similarly minimal in humans, where relative cardiac motion is reduced. Hyperpolarized 129 Xe diffusion morphometry accurately and noninvasively maps the dimensions of lung microstructure, suggesting it can quantify the pulmonary microstructure in mouse models of lung disease.

Mapping kidney tubule diameter ex vivo by diffusion MRI

Tubular pathologies are a common feature of kidney disease. Current metrics to assess kidney health, in vivo or in transplant, are generally based on urinary or serum biomarkers and pathological findings from kidney biopsies. Biopsies, usually taken from the kidney cortex, are invasive and prone to sampling error. Tools to directly and noninvasively measure tubular pathology could provide a new approach to assess kidney health. This study used diffusion magnetic resonance imaging (dMRI) as a noninvasive tool to measure the size of the tubular lumen in ex vivo, perfused kidneys. We first used Monte Carlo simulations to demonstrate that dMRI is sensitive to restricted tissue water diffusion at the scale of the kidney tubule. We applied dMRI and biophysical modeling to examine the distribution of tubular diameters in ex vivo, fixed kidneys from mice, rats, and a human donor. The biophysical model to fit the dMRI signal was based on a superposition of freely diffusing water and water diffusing inside infinitely long cylinders of different diameters. Tubular diameters measured by dMRI were within 10% of those measured by histology within the same tissue. Finally, we applied dMRI to investigate kidney pathology in a mouse model of folic-acid-induced acute kidney injury. dMRI detected heterogeneity in the distribution of tubules within the kidney cortex of mice with acute kidney injury compared with control mice. We conclude that dMRI can be used to measure the distribution of tubule diameters in the kidney cortex ex vivo and that dMRI may provide a new noninvasive biomarker of tubular pathology.NEW & NOTEWORTHY Tubular pathologies are a common feature of kidney disease. Current metrics to assess kidney health, in vivo or in transplant, are generally based on urinary or serum biomarkers and pathological findings from kidney biopsies. Diffusion MRI can be used to measure the distribution of tubule diameters in the kidney cortex ex vivo and may provide a new noninvasive biomarker of tubular pathology.

Finite element simulations of hyperpolarized gas DWI in micro-CT meshes of acinar airways: validating the cylinder and stretched exponential models of lung microstructural length scales

PURPOSE

This work assesses the accuracy of the stretched exponential (SEM) and cylinder models of lung microstructural length scales that can be derived from hyperpolarized gas DWI. This was achieved by simulating ³ He and ¹²⁹ Xe DWI signals within two micro-CT-derived realistic acinar airspace meshes that represent healthy and idiopathic pulmonary fibrosis lungs.

METHODS

The healthy and idiopathic pulmonary fibrosis acinar airway meshes were derived from segmentations of 3D micro-CT images of excised human lungs and meshed for finite element simulations of the Bloch-Torrey equations. ³ He and ¹²⁹ Xe multiple b value DWI experiments across a range of diffusion times (³ He Δ = 1.6 ms; ¹²⁹ Xe Δ = 5 to 20 ms) were simulated in each mesh. Global SEM mean diffusive length scale and cylinder model mean chord length value was derived from each finite element simulation and compared against each mesh's mean linear intercept length, calculated from intercept length measurements within micro-CT segmentation masks.

RESULTS

The SEM-derived mean diffusive length scale was within ±10% of the mean linear intercept length for simulations with both ³ He (Δ = 1.6 ms) and ¹²⁹ Xe (Δ = 7 to 13 ms) in the healthy mesh, and with ¹²⁹ Xe (Δ = 13 to 20 ms) for the idiopathic pulmonary fibrosis mesh, whereas for the cylinder model-derived mean chord length the closest agreement with mean linear intercept length (11.7% and 22.6% difference) was at ¹²⁹ Xe Δ = 20 ms for both healthy and IPF meshes, respectively.

CONCLUSION

This work validates the use of the SEM for accurate estimation of acinar dimensions and indicates that the SEM is relatively robust across a range of experimental conditions and acinar length scales.

In vivo methods and applications of xenon-129 magnetic resonance

Hyperpolarised gas lung MRI using xenon-129 can provide detailed 3D images of the ventilated lung airspaces, and can be applied to quantify lung microstructure and detailed aspects of lung function such as gas exchange. It is sensitive to functional and structural changes in early lung disease and can be used in longitudinal studies of disease progression and therapy response. The ability of 129Xe to dissolve into the blood stream and its chemical shift sensitivity to its local environment allow monitoring of gas exchange in the lungs, perfusion of the brain and kidneys, and blood oxygenation. This article reviews the methods and applications of in vivo129Xe MR in humans, with a focus on the physics of polarisation by optical pumping, radiofrequency coil and pulse sequence design, and the in vivo applications of 129Xe MRI and MRS to examine lung ventilation, microstructure and gas exchange, blood oxygenation, and perfusion of the brain and kidneys.

Lung microstructure in adolescent idiopathic scoliosis before and after posterior spinal fusion

Adolescent idiopathic scoliosis (AIS) is associated with decreased respiratory quality of life and impaired diaphragm function. Recent hyperpolarized helium (HHe) MRI studies show alveolarization continues throughout adolescence, and mechanical forces are known to impact alveolarization. We therefore hypothesized that patients with AIS would have alterations in alveolar size, alveolar number, or alveolar septal dimensions compared to adolescents without AIS, and that posterior spinal fusion (PSF) might reverse these differences. We conducted a prospective observational trial using HHe MRI to test for changes in alveolar microstructure in control and AIS subjects at baseline and one year. After obtaining written informed consent from subjects’ legal guardians and assent from the subjects, we performed HHe and proton MRI in 14 AIS and 16 control subjects aged 8–21 years. The mean age of control subjects (12.9 years) was significantly less than AIS (14.9 years, p = 0.003). At baseline, there were no significant differences in alveolar size, number, or alveolar duct morphometry between AIS and control subjects or between the concave (compressed) and convex (expanded) lungs of AIS subjects. At one year after PSF AIS subjects had an increase in alveolar density in the formerly convex lung (p = 0.05), likely reflecting a change in thoracic anatomy, but there were no other significant changes in lung microstructure. Modeling of alveolar size over time demonstrated similar rates of alveolar growth in control and AIS subjects in both right and left lungs pre- and post-PSF. Although this study suffered from poor age-matching, we found no evidence that AIS or PSF impacts lung microstructure.

Trial registration: Clinical trial registration number NCT03539770.

Application of a stretched-exponential model for morphometric analysis of accelerated diffusion-weighted 129Xe MRI of the rat lung

OBJECTIVE

Diffusion-weighted, hyperpolarized ¹²⁹Xe MRI is useful for the characterization of microstructural changes in the lung. A stretched exponential model was proposed for morphometric extraction of the mean chord length (L_m) from diffusion-weighted data. The stretched exponential model enables accelerated mapping of L_m in a single-breathhold using compressed sensing. Our purpose was to compare L_m maps obtained from stretched-exponential model analysis of accelerated versus unaccelerated diffusion-weighted ¹²⁹Xe MRI data obtained from healthy/injured rat lungs.

MATERIAL AND METHODS

L_m maps were generated using a stretched-exponential model analysis of previously acquired fully sampled diffusion-weighted ¹²⁹Xe rat data (b values = 0 … 110 s/cm²) and compared to L_m maps generated from retrospectively undersampled data simulating acceleration factors of 7/10. The data included four control rats and five rats receiving whole-lung irradiation to mimic radiation-induced lung injury. Mean L_m obtained from the accelerated/unaccelerated maps were compared to histological mean linear intercept.

RESULTS

Accelerated L_m estimates were similar to unaccelerated L_m estimates in all rats, and similar to those previously reported (< 12% different). L_m was significantly reduced (p < 0.001) in the irradiated rat cohort (90 ± 20 µm/90 ± 20 µm) compared to the control rats (110 ± 20 µm/100 ± 15 µm) and agreed well with histological mean linear intercept.

DISCUSSION

Accelerated mapping of L_m using a stretched-exponential model analysis is feasible, accurate and agrees with histological mean linear intercept. Acceleration reduces scan time, thus should be considered for the characterization of lung microstructural changes in humans where breath-hold duration is short.

Accelerated 129 Xe MRI morphometry of terminal airspace enlargement: Feasibility in volunteers and those with alpha-1 antitrypsin deficiency

PURPOSE

Multi-b diffusion-weighted hyperpolarized inhaled-gas MRI provides imaging biomarkers of terminal airspace enlargement including ADC and mean linear intercept (L_m ), but clinical translation has been limited because image acquisition requires relatively long or multiple breath-holds that are not well-tolerated by patients. Therefore, we aimed to accelerate single breath-hold 3D multi-b diffusion-weighted ¹²⁹ Xe MRI, using k-space undersampling in imaging direction using a different undersampling pattern for different b-values combined with the stretched exponential model to generate maps of ventilation, apparent transverse relaxation time constant ( T 2 ∗ ), ADC, and L_m values in a single, short breath-hold; accelerated and non-accelerated measurements were directly compared.

METHODS

We evaluated multi-b (0, 12, 20, 30, and 45.5 s/cm² ) diffusion-weighted ¹²⁹ Xe T 2 ∗ /ADC/morphometry estimates using acceleration factor (AF = 1 and 7) and multi-breath sampling in 3 volunteers (HV), and 6 participants with alpha-1 antitrypsin deficiency (AATD).

RESULTS

For the HV subgroup, mean differences of 5%, 2%, and 8% were observed between fully sampled and undersampled k-space for ADC, L_m , and T 2 ∗ values, respectively. For the AATD subgroup, mean differences were 9%, 6%, and 12% between fully sampled and undersampled k-space for ADC, L_m and T 2 ∗ values, respectively. Although mean differences of 1% and 4.5% were observed between accelerated and multi-breath sampled ADC and L_m values, respectively, mean ADC/L_m estimates were not significantly different from corresponding mean ADC^M /L_m ^M or mean ADC^A /L_m ^A estimates (all P > 0.60 , ^A = undersampled and ^M = multi-breath sampled).

CONCLUSIONS

Accelerated multi-b diffusion-weighted ¹²⁹ Xe MRI is feasible at AF = 7 for generating pulmonary ADC and L_m in AATD and normal lung.

Alveolar dynamics during mechanical ventilation in the healthy and injured lung

Mechanical ventilation is a life-saving therapy in patients with acute respiratory distress syndrome (ARDS). However, mechanical ventilation itself causes severe co-morbidities in that it can trigger ventilator-associated lung injury (VALI) in humans or ventilator-induced lung injury (VILI) in experimental animal models. Therefore, optimization of ventilation strategies is paramount for the effective therapy of critical care patients. A major problem in the stratification of critical care patients for personalized ventilation settings, but even more so for our overall understanding of VILI, lies in our limited insight into the effects of mechanical ventilation at the actual site of injury, i.e., the alveolar unit. Unfortunately, global lung mechanics provide for a poor surrogate of alveolar dynamics and methods for the in-depth analysis of alveolar dynamics on the level of individual alveoli are sparse and afflicted by important limitations. With alveolar dynamics in the intact lung remaining largely a "black box," our insight into the mechanisms of VALI and VILI and the effectiveness of optimized ventilation strategies is confined to indirect parameters and endpoints of lung injury and mortality.In the present review, we discuss emerging concepts of alveolar dynamics including alveolar expansion/contraction, stability/instability, and opening/collapse. Many of these concepts remain still controversial, in part due to limitations of the different methodologies applied. We therefore preface our review with an overview of existing technologies and approaches for the analysis of alveolar dynamics, highlighting their individual strengths and limitations which may provide for a better appreciation of the sometimes diverging findings and interpretations. Joint efforts combining key technologies in identical models to overcome the limitations inherent to individual methodologies are needed not only to provide conclusive insights into lung physiology and alveolar dynamics, but ultimately to guide critical care patient therapy.

MRI and CT lung biomarkers: Towards an in vivo understanding of lung biomechanics

BACKGROUND

The biomechanical properties of the lung are necessarily dependent on its structure and function, both of which are complex and change over time and space. This makes in vivo evaluation of lung biomechanics and a deep understanding of lung biomarkers, very challenging. In patients and animal models of lung disease, in vivo evaluations of lung structure and function are typically made at the mouth and include spirometry, multiple-breath gas washout tests and the forced oscillation technique. These techniques, and the biomarkers they provide, incorporate the properties of the whole organ system including the parenchyma, large and small airways, mouth, diaphragm and intercostal muscles. Unfortunately, these well-established measurements mask regional differences, limiting their ability to probe the lung's gross and micro-biomechanical properties which vary widely throughout the organ and its subcompartments. Pulmonary imaging has the advantage in providing regional, non-invasive measurements of healthy and diseased lung, in vivo. Here we summarize well-established and emerging lung imaging tools and biomarkers and how they may be used to generate lung biomechanical measurements.

METHODS

We review well-established and emerging lung anatomical, microstructural and functional imaging biomarkers generated using synchrotron x-ray tomographic-microscopy (SRXTM), micro-x-ray computed-tomography (micro-CT), clinical CT as well as magnetic resonance imaging (MRI).

FINDINGS

Pulmonary imaging provides measurements of lung structure, function and biomechanics with high spatial and temporal resolution. Imaging biomarkers that reflect the biomechanical properties of the lung are now being validated to provide a deeper understanding of the lung that cannot be achieved using measurements made at the mouth.

Advanced pulmonary MRI to quantify alveolar and acinar duct abnormalities: Current status and future clinical applications

There are serious clinical gaps in our understanding of chronic lung disease that require novel, sensitive, and noninvasive in vivo measurements of the lung parenchyma to measure disease pathogenesis and progressive changes over time as well as response to treatment. Until recently, our knowledge and appreciation of the tissue changes that accompany lung disease has depended on ex vivo biopsy and concomitant histological and stereological measurements. These measurements have revealed the underlying pathologies that drive lung disease and have provided important observations about airway occlusion, obliteration of the terminal bronchioles and airspace enlargement, or fibrosis and their roles in disease initiation and progression. ex vivo tissue stereology and histology are the established gold standards and, more recently, micro-computed tomography (CT) measurements of ex vivo tissue samples has also been employed to reveal new mechanistic findings about the progression of obstructive lung disease in patients. While these approaches have provided important understandings using ex vivo analysis of excised samples, recently developed hyperpolarized noble gas MRI methods provide an opportunity to noninvasively measure acinar duct and terminal airway dimensions and geometry in vivo, and, without radiation burden. Therefore, in this review we summarize emerging pulmonary MRI morphometry methods that provide noninvasive in vivo measurements of the lung in patients with bronchopulmonary dysplasia and chronic obstructive pulmonary disease, among others. We discuss new findings, future research directions, as well as clinical opportunities to address current gaps in patient care and for testing of new therapies. Level of Evidence: 5 Technical Efficacy: Stage 5 J. Magn. Reson. Imaging 2019;50:28-40.

3He diffusion MRI in human lungs

Hyperpolarized 3He gas allows the air spaces of the lungs to be imaged via MRI. Imaging of restricted diffusion is addressed here, which allows the microstructure of the lung to be characterized through the physical restrictions to gas diffusion presented by airway and alveolar walls in the lung. Measurements of the apparent diffusion coefficient (ADC) of 3He at time scales of milliseconds and seconds are compared; measurement of acinar airway sizes by determination of the microscopic anisotropy of diffusion is discussed. This is where Dr. JJH Ackerman's influence was greatest in aiding the formation of the Washington University 3He group, involving early a combination of physicists, radiologists, and surgeons, as the first applications of 3He ADC were to COPD and its destruction/modification of lung microstructure via emphysema. The sensitivity of the method to early COPD is demonstrated, as is its validation by direct comparison to histology. More recently the method has been used broadly in adult and pediatric obstructive lung diseases, from severe asthma to cystic fibrosis to bronchopulmonary dysplasia, a result of premature birth. These applications of the technique are discussed briefly.

Pulmonary 3He Magnetic Resonance Imaging Biomarkers of Regional Airspace Enlargement in Alpha-1 Antitrypsin Deficiency

RATIONALE AND OBJECTIVES

Thoracic x-ray computed tomography (CT) and hyperpolarized ³He magnetic resonance imaging (MRI) provide quantitative measurements of airspace enlargement in patients with emphysema. For patients with panlobular emphysema due to alpha-1 antitrypsin deficiency (AATD), sensitive biomarkers of disease progression and response to therapy have been difficult to develop and exploit, especially those biomarkers that correlate with outcomes like quality of life. Here, our objective was to generate and compare CT and diffusion-weighted inhaled-gas MRI measurements of emphysema including apparent diffusion coefficient (ADC) and MRI-derived mean linear intercept (L_m) in patients with AATD, chronic obstructive pulmonary disease (COPD) ex-smokers, and elderly never-smokers.

MATERIALS AND METHODS

We enrolled patients with AATD (n = 8; 57 ± 7 years), ex-smokers with COPD (n = 8; 77 ± 6 years), and a control group of never-smokers (n = 5; 64 ± 2 years) who underwent thoracic CT, MRI, spirometry, plethysmography, the St. George's Respiratory Questionnaire, and the 6-minute walk test during a single 2-hour visit. MRI-derived ADC, L_m, surface-to-volume ratio, and ventilation defect percent were generated for the apical, basal, and whole lung as was CT lung area ≤-950 Hounsfield units (RA₉₅₀), low attenuating clusters, and airway count.

RESULTS

In patients with AATD, there was a significantly different MRI-derived ADC (P = .03), L_m (P < .0001), and surface-to-volume ratio (P < .0001), but not diffusing capacity of carbon monoxide, residual volume or total lung capacity, or CT RA₉₅₀ (P > .05) compared to COPD ex-smokers with a significantly different St. George's Respiratory Questionnaire.

CONCLUSIONS

In this proof-of-concept demonstration, we evaluated CT and MRI lung emphysema measurements and observed significantly worse MRI biomarkers of emphysema in patients with AATD compared to patients with COPD, although CT RA₉₅₀ and diffusing capacity of carbon monoxide were not significantly different, underscoring the sensitivity of MRI measurements of AATD emphysema.

3D diffusion-weighted 129 Xe MRI for whole lung morphometry

Purpose

To obtain whole lung morphometry measurements from ¹²⁹Xe in a single breath‐hold with 3D multiple b‐value ¹²⁹Xe diffusion‐weighted MRI (DW‐MRI) with an empirically optimized diffusion time and compressed sensing for scan acceleration.

Methods

Prospective three‐fold undersampled 3D multiple b‐value hyperpolarized ¹²⁹Xe DW‐MRI datasets were acquired, and the diffusion time (Δ) was iterated so as to provide diffusive length scale (Lm_D) estimates from the stretched exponential model (SEM) that are comparable to those from ³He. The empirically optimized ¹²⁹Xe diffusion time was then implemented with a four‐fold undersampling scheme and was prospectively benchmarked against ³He measurements in a cohort of five healthy volunteers, six ex‐smokers, and two chronic obstructive pulmonary disease patients using both SEM‐derived Lm_D and cylinder model (CM)‐derived mean chord length (Lm).

Results

Good agreement between the mean ¹²⁹Xe and ³He Lm_D (mean difference, 2.2%) and Lm (mean difference, 1.1%) values was obtained in all subjects at an empirically optimized ¹²⁹Xe Δ = 8.5 ms.

Conclusion

Compressed sensing has facilitated single‐breath 3D multiple b‐value ¹²⁹Xe DW‐MRI acquisitions, and results at ¹²⁹Xe Δ = 8.5 ms indicate that ¹²⁹Xe provides a viable alternative to ³He for whole lung morphometry mapping with either the SEM or CM. Magn Reson Med 79:2986–2995, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Diffusion lung imaging with hyperpolarized gas MRI

Lung imaging using conventional 1 H MRI presents great challenges because of the low density of lung tissue, lung motion and very fast lung tissue transverse relaxation (typical T2 * is about 1-2 ms). MRI with hyperpolarized gases (3 He and 129 Xe) provides a valuable alternative because of the very strong signal originating from inhaled gas residing in the lung airspaces and relatively slow gas T2 * relaxation (typical T2 * is about 20-30 ms). However, in vivo human experiments should be performed very rapidly - usually during a single breath-hold. In this review, we describe the recent developments in diffusion lung MRI with hyperpolarized gases. We show that a combination of the results of modeling of gas diffusion in lung airspaces and diffusion measurements with variable diffusion-sensitizing gradients allows the extraction of quantitative information on the lung microstructure at the alveolar level. From an MRI scan of less than 15 s, this approach, called in vivo lung morphometry, allows the provision of quantitative values and spatial distributions of the same physiological parameters as measured by means of 'standard' invasive stereology (mean linear intercept, surface-to-volume ratio, density of alveoli, etc.). In addition, the approach makes it possible to evaluate some advanced Weibel parameters characterizing lung microstructure: average radii of alveolar sacs and ducts, as well as the depth of their alveolar sleeves. Such measurements, providing in vivo information on the integrity of pulmonary acinar airways and their changes in different diseases, are of great importance and interest to a broad range of physiologists and clinicians. We also discuss a new type of experiment based on the in vivo lung morphometry technique combined with quantitative computed tomography measurements, as well as with gradient echo MRI measurements of hyperpolarized gas transverse relaxation in the lung airspaces. Such experiments provide additional information on the blood vessel volume fraction, specific gas volume and length of the acinar airways, and allow the evaluation of lung parenchymal and non-parenchymal tissue. Copyright © 2015 John Wiley & Sons, Ltd.

Fundamentals of diffusion MRI physics

Diffusion MRI is commonly considered the "engine" for probing the cellular structure of living biological tissues. The difficulty of this task is threefold. First, in structurally heterogeneous media, diffusion is related to structure in quite a complicated way. The challenge of finding diffusion metrics for a given structure is equivalent to other problems in physics that have been known for over a century. Second, in most cases the MRI signal is related to diffusion in an indirect way dependent on the measurement technique used. Third, finding the cellular structure given the MRI signal is an ill-posed inverse problem. This paper reviews well-established knowledge that forms the basis for responding to the first two challenges. The inverse problem is briefly discussed and the reader is warned about a number of pitfalls on the way.

Pulmonary MRI morphometry modeling of airspace enlargement in chronic obstructive pulmonary disease and alpha-1 antitrypsin deficiency

PURPOSE

We generated lung morphometry measurements using single-breath diffusion-weighted MRI and three different acinar duct models in healthy participants and patients with emphysema stemming from chronic obstructive lung disease (COPD) and alpha-1 antitrypsin deficiency (AATD).

METHODS

Single-breath-inhaled ³ He MRI with five diffusion sensitizations (b-value = 0, 1.6, 3.2, 4.8, and 6.4 s/cm² ) was used, and signal intensities were fit using a cylindrical and single-compartment acinar-duct model to estimate MRI-derived mean linear intercept (L_m ) and surface-to-volume ratio (S/V). A stretched exponential model was also developed to estimate the mean airway length and L_m .

RESULTS

We evaluated 42 participants, including 15 elderly never-smokers (69 ± 5 years), 12 ex-smokers without COPD (67 ± 11 years), 9 COPD ex-smokers (80 ± 6 years), and 6 AATD patients (59 ± 6 years). In the never- and ex-smokers, the diffusing capacity of the lung for carbon monoxide (DL_CO ) and computed tomography relative area of less than -950 Hounsfield units (RA₉₅₀ ) were normal, but these were abnormal in the COPD and AATD patients, which is reflective of emphysema. Although cylindrical and stretched-exponential-model estimates of L_m and S/V were not significantly different, the single-compartment-model estimates were significantly different (P < 0.05) for the never- and ex-smoker subgroups. All models estimated significantly worse L_m and S/V in the AATD and COPD subgroups compared with the never- and ex-smokers without emphysema.

CONCLUSIONS

Differences in airspace enlargement may be estimated using L_m and S/V, generated using MRI and a stretched-exponential or cylindrical model of the acinar ducts. Magn Reson Med 79:439-448, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Detection of the mild emphysema by quantification of lung respiratory airways with hyperpolarized xenon diffusion MRI

PURPOSE

To demonstrate the feasibility to quantify the lung respiratory airway in vivo with hyperpolarized xenon diffusion magnetic resonance imaging (MRI), which is able to detect mild emphysema in the rat model.

MATERIALS AND METHODS

The lung respiratory airways were quantified in vivo using hyperpolarized xenon diffusion MRI (7T) with eight b values (5, 10, 15, 20, 25, 30, 35, 40 s/cm² ) in five control rats and five mild emphysematous rats, which were induced by elastase. The morphological results from histology were acquired and used for comparison.

RESULTS

The parameters D_L (longitudinal diffusion coefficient), r (internal radius), h (alveolar sleeve depth), Lm (mean linear intercept), and S/V (surface area to lung volume ratio) derived from the hyperpolarized xenon diffusion MRI in the emphysematous group showed significant differences from those in the control group (P < 0.05). Additionally, these parameters correlated well with the Lm obtained by the traditional histological sections (Pearson's correlation coefficients >0.8).

CONCLUSION

The lung respiratory airways can be quantified by hyperpolarized xenon diffusion MRI, showing the potential for mild emphysema diagnosis. Also, the study suggested that the hyperpolarized xenon D_L is more sensitive than D_T (transverse diffusion coefficient) to detect mild emphysema.

LEVEL OF EVIDENCE

1 J. Magn. Reson. Imaging 2017;45:879-888.

Whole lung morphometry with 3D multiple b-value hyperpolarized gas MRI and compressed sensing

Purpose

To demonstrate three‐dimensional (3D) multiple b‐value diffusion‐weighted (DW) MRI of hyperpolarized ³He gas for whole lung morphometry with compressed sensing (CS).

Methods

A fully‐sampled, two b‐value, 3D hyperpolarized ³He DW‐MRI dataset was acquired from the lungs of a healthy volunteer and retrospectively undersampled in the k_y and k_z phase‐encoding directions for CS simulations. Optimal k‐space undersampling patterns were determined by minimizing the mean absolute error between reconstructed and fully‐sampled ³He apparent diffusion coefficient (ADC) maps. Prospective three‐fold, undersampled, 3D multiple b‐value ³He DW‐MRI datasets were acquired from five healthy volunteers and one chronic obstructive pulmonary disease (COPD) patient, and the mean values of maps of ADC and mean alveolar dimension (Lm_D) were validated against two‐dimensional (2D) and 3D fully‐sampled ³He DW‐MRI experiments.

Results

Reconstructed undersampled datasets showed no visual artifacts and good preservation of the main image features and quantitative information. A good agreement between fully‐sampled and prospective undersampled datasets was found, with a mean difference of +3.4% and +5.1% observed in mean global ADC and Lm_D values, respectively. These differences were within the standard deviation range and consistent with values reported from healthy and COPD lungs.

Conclusions

Accelerated CS acquisition has facilitated 3D multiple b‐value ³He DW‐MRI scans in a single breath‐hold, enabling whole lung morphometry mapping. Magn Reson Med 77:1916–1925, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Fluorine-19 MRI Contrast Agents for Cell Tracking and Lung Imaging

Fluorine-19 (¹⁹F)-based contrast agents for magnetic resonance imaging stand to revolutionize imaging-based research and clinical trials in several fields of medical intervention. First, their use in characterizing in vivo cell behavior may help bring cellular therapy closer to clinical acceptance. Second, their use in lung imaging provides novel noninvasive interrogation of the ventilated airspaces without the need for complicated, hard-to-distribute hardware. This article reviews the current state of ¹⁹F-based cell tracking and lung imaging using magnetic resonance imaging and describes the link between the methods across these fields and how they may mutually benefit from solutions to mutual problems encountered when imaging ¹⁹F-containing compounds, as well as hardware and software advancements.

Hyperpolarized 83Kr magnetic resonance imaging of alveolar degradation in a rat model of emphysema

Hyperpolarized ⁸³Kr surface quadrupolar relaxation (SQUARE) generates MRI contrast that was previously shown to correlate with surface-to-volume ratios in porous model surface systems. The underlying physics of SQUARE contrast is conceptually different from any other current MRI methodology as the method uses the nuclear electric properties of the spin I = 9/2 isotope ⁸³Kr. To explore the usage of this non-radioactive isotope for pulmonary pathophysiology, MRI SQUARE contrast was acquired in excised rat lungs obtained from an elastase-induced model of emphysema. A significant ⁸³Kr T₁ relaxation time increase in the SQUARE contrast was found in the elastase-treated lungs compared with the baseline data from control lungs. The SQUARE contrast suggests a reduction in pulmonary surface-to-volume ratio in the emphysema model that was validated by histology. The finding supports usage of ⁸³Kr SQUARE as a new biomarker for surface-to-volume ratio changes in emphysema.

Probing lung microstructure with hyperpolarized noble gas diffusion MRI: theoretical models and experimental results

The introduction of hyperpolarized gases ((3)He and (129)Xe) has opened the door to applications for which gaseous agents are uniquely suited-lung MRI. One of the pulmonary applications, diffusion MRI, relies on measuring Brownian motion of inhaled hyperpolarized gas atoms diffusing in lung airspaces. In this article we provide an overview of the theoretical ideas behind hyperpolarized gas diffusion MRI and the results obtained over the decade-long research. We describe a simple technique based on measuring gas apparent diffusion coefficient (ADC) and an advanced technique, in vivo lung morphometry, that quantifies lung microstructure both in terms of Weibel parameters (acinar airways radii and alveolar depth) and standard metrics (mean linear intercept, surface-to-volume ratio, and alveolar density) that are widely used by lung researchers but were previously available only from invasive lung biopsy. This technique has the ability to provide unique three-dimensional tomographic information on lung microstructure from a less than 15 s MRI scan with results that are in good agreement with direct histological measurements. These safe and sensitive diffusion measurements improve our understanding of lung structure and functioning in health and disease, providing a platform for monitoring the efficacy of therapeutic interventions in clinical trials.

Microstructural Analysis of Peripheral Lung Tissue through CPMG Inter-Echo Time R2 Dispersion

Since changes in lung microstructure are important indicators for (early stage) lung pathology, there is a need for quantifiable information of diagnostically challenging cases in a clinical setting, e.g. to evaluate early emphysematous changes in peripheral lung tissue. Considering alveoli as spherical air-spaces surrounded by a thin film of lung tissue allows deriving an expression for Carr-Purcell-Meiboom-Gill transverse relaxation rates R₂ with a dependence on inter-echo time, local air-tissue volume fraction, diffusion coefficient and alveolar diameter, within a weak field approximation. The model relaxation rate exhibits the same hyperbolic tangent dependency as seen in the Luz-Meiboom model and limiting cases agree with Brooks et al. and Jensen et al. In addition, the model is tested against experimental data for passively deflated rat lungs: the resulting mean alveolar radius of R_A = 31.46 ± 13.15 μm is very close to the literature value (∼34 μm). Also, modeled radii obtained from relaxometer measurements of ageing hydrogel foam (that mimics peripheral lung tissue) are in good agreement with those obtained from μCT images of the same foam (mean relative error: 0.06 ± 0.01). The model’s ability to determine the alveolar radius and/or air volume fraction will be useful in quantifying peripheral lung microstructure.

Noninvasive quantification of alveolar morphometry in elderly never- and ex-smokers

Diffusion-weighted magnetic resonance imaging (MRI) provides a way to generate in vivo lung images with contrast sensitive to the molecular displacement of inhaled gas at subcellular length scales. Here, we aimed to evaluate hyperpolarized ³He MRI estimates of the alveolar dimensions in 38 healthy elderly never-smokers (73 ± 6 years, 15 males) and 21 elderly ex-smokers (70 ± 10 years, 14 males) with (n = 8, 77 ± 6 years) and without emphysema (n = 13, 65 ± 10 years). The ex-smoker and never-smoker subgroups were significantly different for FEV₁/FVC (P = 0.0001) and DL_CO (P = 0.009); while ex-smokers with emphysema reported significantly diminished FEV₁/FVC (P = 0.02) and a trend toward lower DL_CO (P = 0.05) than ex-smokers without emphysema. MRI apparent diffusion coefficients (ADC) and CT measurements of emphysema (relative area–CT density histogram, RA₉₅₀) were significantly different (P = 0.001 and P = 0.007) for never-smoker and ex-smoker subgroups. In never-smokers, the MRI estimate of mean linear intercept (260 ± 27 μm) was significantly elevated as compared to the results previously reported in younger never-smokers (210 ± 30 μm), and trended smaller than in the age-matched ex-smokers (320 ± 72 μm, P = 0.06) evaluated here. Never-smokers also reported significantly smaller internal (220 ± 24 μm, P = 0.01) acinar radius but greater alveolar sheath thickness (120 ± 4 μm, P < 0.0001) than ex-smokers. Never-smokers were also significantly different than ex-smokers without emphysema for alveolar sheath thickness but not ADC, while ex-smokers with emphysema reported significantly different ADC but not alveolar sheath thickness compared to ex-smokers without CT evidence of emphysema. Differences in alveolar measurements in never- and ex-smokers demonstrate the sensitivity of MRI measurements to the different effects of smoking and aging on acinar morphometry.

CT-Definable Subtypes of Chronic Obstructive Pulmonary Disease: A Statement of the Fleischner Society

The purpose of this statement is to describe and define the phenotypic abnormalities that can be identified on visual and quantitative evaluation of computed tomographic (CT) images in subjects with chronic obstructive pulmonary disease (COPD), with the goal of contributing to a personalized approach to the treatment of patients with COPD. Quantitative CT is useful for identifying and sequentially evaluating the extent of emphysematous lung destruction, changes in airway walls, and expiratory air trapping. However, visual assessment of CT scans remains important to describe patterns of altered lung structure in COPD. The classification system proposed and illustrated in this article provides a structured approach to visual and quantitative assessment of COPD. Emphysema is classified as centrilobular (subclassified as trace, mild, moderate, confluent, and advanced destructive emphysema), panlobular, and paraseptal (subclassified as mild or substantial). Additional important visual features include airway wall thickening, inflammatory small airways disease, tracheal abnormalities, interstitial lung abnormalities, pulmonary arterial enlargement, and bronchiectasis.

Probing lung microstructure with hyperpolarized 3He gradient echo MRI

In this paper we demonstrate that gradient echo MRI with hyperpolarized (3)He gas can be used for simultaneously extracting in vivo information about lung ventilation properties, alveolar geometrical parameters, and blood vessel network structure. This new approach is based on multi-gradient-echo experimental measurements of hyperpolarized (3)He gas MRI signal from human lungs and a proposed theoretical model of this signal. Based on computer simulations of (3)He atoms diffusing in the acinar airway tree in the presence of an inhomogeneous magnetic field induced by the susceptibility differences between lung tissue (alveolar septa, blood vessels) and lung airspaces, we derive analytical expressions relating the time-dependent MR signal to the geometrical parameters of acinar airways and the blood vessel network. Data obtained on eight healthy volunteers are in good agreement with literature values. This information is complementary to the information obtained by means of the in vivo lung morphometry technique with hyperpolarized 3He diffusion MRI previously developed by our group, and opens new opportunities to study lung microstructure in health and disease.

Commentary on "The influence of lung airways branching structure and diffusion time on measurements and models of short-range 3He gas MR diffusion"

In a recently published paper by Parra-Robles and Wild, the authors challenge the in vivo lung morphometry technique (based on hyperpolarized gas diffusion MRI) developed by our Washington University research group. In this Commentary we demonstrate that the main conclusion of Parra-Robles and Wild, that our MRI-based lung morphometry technique "produces inaccurate estimates of the airway dimensions", does not have any scientific basis and is not in agreement with the considerable body of peer-reviewed scientific reports as well as with Parra-Robles and Wild's own data. On the contrary, our technique has a strong theoretical background, is validated, and provides accurate 3D tomographic information on lung microstructural parameters previously available only from invasive biopsy specimens. This technique has already produced a number of results related to lung morphology and function that were not previously available. In our Commentary we also discuss a number of other incorrect statements in and shortcomings of Parra-Robles and Wild's paper.

Response to Commentary on "The influence of lung airways branching structure and diffusion time on measurements and models of short-range 3He gas MR diffusion"

Our extensive investigation of the cylinder model theory through numerical modelling and purpose-designed experiments has demonstrated that it does produce inaccurate estimates of airway dimensions at all diffusion times currently used. This is due to a variety of effects: incomplete treatment of non-Gaussian effects, finite airway size, branching geometry, background susceptibility gradients and diffusion time dependence of the (3)He MR diffusion behaviour in acinar airways. The cylinder model is a good starting point for the development of a lung morphometry technique from (3)He diffusion MR but its limitations need to be understood and documented in the interest of reliable clinical interpretation.

Longitudinal, noninvasive monitoring of compensatory lung growth in mice after pneumonectomy via (3)He and (1)H magnetic resonance imaging

In rodents and some other mammals, partial pneumonectomy (PNX) of adult lungs results in rapid compensatory lung growth. In the past, quantification of compensatory lung growth and realveolarization could only be accomplished after killing the animal, removal of lungs, and histologic analysis of lungs at single time points. Hyperpolarized (3)He diffusion magnetic resonance imaging (MRI) allows in vivo morphometry of human lungs; our group has adapted this technique for application to mouse lungs. Through imaging, we can obtain maps of lung microstructural parameters that allow quantification of morphometric and physiologic measurements. In this study, we employed our (3)He MRI technique to image in vivo morphometry at baseline and to serially assess compensatory growth after left PNX of mice. (1)H and hyperpolarized (3)He diffusion MRI were performed at baseline (pre-PNX), 3-days, and 30-days after PNX. Compared with the individual mouse's own baseline, MRI was able to detect and serially quantify changes in lung volume, alveolar surface area, alveolar number, and regional changes in alveolar size that occurred during the course of post-PNX lung growth. These results are consistent with morphometry measurements reported in the literature for mouse post-PNX compensatory lung growth. In addition, we were also able to serially assess and quantify changes in the physiologic parameter of lung compliance during the course of compensatory lung growth; this was consistent with flexiVent data. With these techniques, we now have a noninvasive, in vivo method to serially assess the effectiveness of therapeutic interventions on post-PNX lung growth in the same mouse.

Imaging the interaction of atelectasis and overdistension in surfactant-depleted lungs

OBJECTIVE

Atelectasis and surfactant depletion may contribute to greater distension-and thereby injury-of aerated lung regions; recruitment of atelectatic lung may protect these regions by attenuating such overdistension. However, the effects of atelectasis (and recruitment) on aerated airspaces remain elusive. We tested the hypothesis that during mechanical ventilation, surfactant depletion increases the dimensions of aerated airspaces and that lung recruitment reverses these changes.

DESIGN

Prospective imaging study in an animal model.

SETTING

Research imaging facility.

SUBJECTS

Twenty-seven healthy Sprague Dawley rats.

INTERVENTIONS

Surfactant depletion was obtained by saline lavage in anesthetized, ventilated rats. Alveolar recruitment was accomplished using positive end-expiratory pressure and exogenous surfactant administration.

MEASUREMENTS AND MAIN RESULTS

Airspace dimensions were estimated by measuring the apparent diffusion coefficient of He, using diffusion-weighted hyperpolarized gas magnetic resonance imaging. Atelectasis was demonstrated using computerized tomography and by measuring oxygenation. Saline lavage increased atelectasis (increase in nonaerated tissue from 1.2% to 13.8% of imaged area, p < 0.001), and produced a concomitant increase in mean apparent diffusion coefficient (~33%, p < 0.001) vs. baseline; the heterogeneity of the computerized tomography signal and the variance of apparent diffusion coefficient were also increased. Application of positive end-expiratory pressure and surfactant reduced the mean apparent diffusion coefficient (~23%, p < 0.001), and its variance, in parallel to alveolar recruitment (i.e., less computerized tomography densities and heterogeneity, increased oxygenation).

CONCLUSIONS

Overdistension of aerated lung occurs during atelectasis is detectable using clinically relevant magnetic resonance imaging technology, and could be a key factor in the generation of lung injury during mechanical ventilation. Lung recruitment by higher positive end-expiratory pressure and surfactant administration reduces airspace distension.

Pulmonary 3He magnetic resonance imaging of childhood asthma

BACKGROUND

Magnetic resonance imaging (MRI) with (3)He does not require ionizing radiation and has been shown to detect regional abnormalities in lung ventilation and structure in adults with asthma, but the method has not been extended to children with asthma. Measurements of regional lung ventilation and microstructure in subjects with childhood asthma could advance our understanding of disease mechanisms.

OBJECTIVE

We sought to determine whether (3)He MRI in children can identify abnormalities related to the diagnosis of asthma or prior history of respiratory illness.

METHODS

Forty-four children aged 9 to 10 years were recruited from a birth cohort at increased risk of asthma and allergic diseases. For each subject, a time-resolved 3-dimensional image series and a 3-dimensional diffusion-weighted image were acquired in separate breathing maneuvers. The numbers and sizes of ventilation defects were scored, and regional maps and statistics of average (3)He diffusion lengths were calculated.

RESULTS

Children with mild-to-moderate asthma had lower average root-mean-square diffusion length (X(rms)) values (P = .004), increased regional SD of diffusion length values (P = .03), and higher defect scores (P = .03) than those without asthma. Children with histories of wheezing illness with rhinovirus infection before the third birthday had lower X(rms) values (P = .01) and higher defect scores (P = .05).

CONCLUSION

MRI with (3)He detected more and larger regions of ventilation defect and a greater degree of restricted gas diffusion in children with asthma compared with those seen in children without asthma. These measures are consistent with regional obstruction and smaller and more regionally variable dimensions of the peripheral airways and alveolar spaces.

The influence of lung airways branching structure and diffusion time on measurements and models of short-range 3He gas MR diffusion

Hyperpolarized (3)He diffusion experiments have been shown to be sensitive to changes in acinar structure due to emphysematous lung disease. Extracting quantitative information about lung microstructure from the diffusion signal is complicated due its dependence on a number of factors including diffusion time and the complex branching acinar geometry. A theoretical model (cylinder model) has been proposed as a means of estimating acinar airway dimensions from measured diffusivities. This model assumes that the effects of acinar branching geometry and finite airway length upon (3)He diffusion behaviour are negligible. In this work, we use finite element simulations of diffusion in a model of branching alveolar ducts to investigate in detail the effects of acinar branching structure and finite airway length on short-range (3)He diffusion measurements. The results show that branching effects have a significant influence upon (3)He diffusivity, even at short diffusion times. The expressions of the cylinder model theory do not account for significant dependences upon diffusion time, branching geometry and airway length, as a consequence of the oversimplified geometrical model used. The effect of diffusion time on (3)He ADC was also investigated through experiments with healthy human volunteers. The results demonstrate that the cylinder model can produce inaccurate estimates of the airway dimensions as a consequence of incompletely accounting for the diffusion-time dependence in the model equations and confirmed the predicted limitations of the cylinder model for reliable lung morphometry measurements. The results and models presented in this work may help in the development of a more realistic theoretical framework for 'in vivo lung morphometry' using (3)He diffusion MR.

Lung morphometry with hyperpolarized 129Xe: theoretical background

The (3) He lung morphometry technique, based on MRI measurements of hyperpolarized (3) He gas diffusion in lung airspaces, provides unique information on the lung microstructure at the alveolar level. In vivo 3D tomographic images of standard morphological parameters (airspace chord length, lung parenchyma surface-to-volume ratio, and number of alveoli per unit volume) can be generated from a rather short (several seconds) MRI scan. The technique is based on a theory of gas diffusion in lung acinar airways and experimental measurements of diffusion-attenuated MRI signal. The present work aims at developing the theoretical background of a similar technique based on hyperpolarized (129) Xe gas. As the diffusion coefficient and gyromagnetic ratio of (129) Xe gas are substantially different from those of (3) He gas, the specific details of the theory and experimental measurements with (129) Xe should be amended. We establish phenomenological relationships between acinar airway geometrical parameters and the diffusion-attenuated MR signal for human and small animal lungs, both normal lungs and lungs with mild emphysema. Optimal diffusion times are shown to be about 5 ms for human and 1.3 ms for small animals. The expected uncertainties in measuring main morphometrical parameters of the lungs are estimated in the framework of Bayesian probability theory.

Morphometric changes in the human pulmonary acinus during inflation

Despite decades of research into the mechanisms of lung inflation and deflation, there is little consensus about whether lung inflation occurs due to the recruitment of new alveoli or by changes in the size and/or shape of alveoli and alveolar ducts. In this study we use in vivo (3)He lung morphometry via MRI to measure the average alveolar depth and alveolar duct radius at three levels of inspiration in five healthy human subjects and calculate the average alveolar volume, surface area, and the total number of alveoli at each level of inflation. Our results indicate that during a 143 ± 18% increase in lung gas volume, the average alveolar depth decreases 21 ±5%, the average alveolar duct radius increases 7 ± 3%, and the total number of alveoli increases by 96 ± 9% (results are means ± SD between subjects; P < 0.001, P < 0.01, and P < 0.00001, respectively, via paired t-tests). Thus our results indicate that in healthy human subjects the lung inflates primarily by alveolar recruitment and, to a lesser extent, by anisotropic expansion of alveolar ducts.