Optimized quantitative mapping of cardiopulmonary oscillations using hyperpolarized 129 Xe gas exchange MRI: Digital phantoms and clinical evaluation in CTEPH.

PURPOSE: The interaction between 129 Xe atoms and pulmonary capillary red blood cells provides cardiogenic signal oscillations that display sensitivity to precapillary and postcapillary pulmonary hypertension. Recently, such oscillations have been spatially mapped, but little is known about optimal reconstruction or sensitivity to artifacts. In this study, we use digital phantom simulations to specifically optimize keyhole reconstruction for oscillation imaging. We then use this optimized method to re-establish healthy reference values and quantitatively evaluate microvascular flow changes in patients with chronic thromboembolic pulmonary hypertension (CTEPH) before and after pulmonary thromboendarterectomy (PTE). METHODS: A six-zone digital lung phantom was designed to investigate the effects of radial views, key radius, and SNR. One-point Dixon 129 Xe gas exchange MRI images were acquired in a healthy cohort (n = 17) to generate a reference distribution and thresholds for mapping red blood cell oscillations. These thresholds were applied to 10 CTEPH participants, with 6 rescanned following PTE. RESULTS: For undersampled acquisitions, a key radius of 0.14 k max $$ 0.14{k}_{\mathrm{max}} $$ was found to optimally resolve oscillation defects while minimizing excessive heterogeneity. CTEPH participants at baseline showed higher oscillation defect + low (32 ± 14%) compared with healthy volunteers (18 ± 12%, p < 0.001). For those scanned both before and after PTE, oscillation defect + low decreased from 37 ± 13% to 23 ± 14% (p = 0.03). CONCLUSIONS: Digital phantom simulations have informed an optimized keyhole reconstruction technique for gas exchange images acquired with standard 1-point Dixon parameters. Our proposed methodology enables more robust quantitative mapping of cardiogenic oscillations, potentially facilitating effective regional quantification of microvascular flow impairment in patients with pulmonary vascular diseases such as CTEPH.

Rapid pulmonary 129Xe ventilation MRI of discharged COVID-19 patients with zigzag sampling.

PURPOSE: To demonstrate the feasibility of zigzag sampling for 3D rapid hyperpolarized 129Xe ventilation MRI in human. METHODS: Zigzag sampling in one direction was combined with gradient-recalled echo sequence (GRE-zigzag-Y) to acquire hyperpolarized 129Xe ventilation images. Image quality was compared with a balanced SSFP (bSSFP) sequence with the same spatial resolution for 12 healthy volunteers (HVs). For another 8 HVs and 9 discharged coronavirus disease 2019 subjects, isotropic resolution 129Xe ventilation images were acquired using zigzag sampling in two directions through GRE-zigzag-YZ. 129Xe ventilation defect percent (VDP) was quantified for GRE-zigzag-YZ and bSSFP acquisitions. Relationships and agreement between these VDP measurements were evaluated using Pearson correlation coefficient (r) and Bland-Altman analysis. RESULTS: For 12 HVs, GRE-zigzag-Y and bSSFP required 2.2 s and 10.5 s, respectively, to acquire 129Xe images with a spatial resolution of 3.96 × 3.96 × 10.5 mm3. Structural similarity index, mean absolute error, and Dice similarity coefficient between the two sets of images and ventilated lung regions were 0.85 ± 0.03, 0.0015 ± 0.0001, and 0.91 ± 0.02, respectively. For another 8 HVs and 9 coronavirus disease 2019 subjects, 129Xe images with a nominal spatial resolution of 2.5 × 2.5 × 2.5 mm3 were acquired within 5.5 s per subject using GRE-zigzag-YZ. VDP provided by GRE-zigzag-YZ was strongly correlated (R2 = 0.93, p < 0.0001) with that generated by bSSFP with minimal biases (bias = -0.005%, 95% limit-of-agreement = [-0.414%, 0.424%]). CONCLUSION: Zigzag sampling combined with GRE sequence provides a way for rapid 129Xe ventilation imaging.

Age, sex, and lung volume dependence of dissolved xenon-129 MRI gas exchange metrics.

PURPOSE: To characterize the dependence of Xe-MRI gas transfer metrics upon age, sex, and lung volume in a group of healthy volunteers. METHODS: Sixty-five subjects with no history of chronic lung disease were assessed with 129Xe-MRI using a four-echo 3D radial spectroscopic imaging sequence and a dose of xenon titrated according to subject height that was inhaled from a lung volume of functional residual capacity (FRC). Imaging was repeated in 34 subjects at total lung capacity (TLC). Regional maps of the fractions of dissolved xenon in red blood cells (RBC), membrane (M), and airspace (Gas) were acquired at an isotropic resolution of 2 cm, from which global averages of the ratios RBC:M, RBC:Gas, and M:Gas were computed. RESULTS: Data from 26 males and 36 females with a median age of 43 y (range: 20-69 y) were of sufficient quality to analyze. Age (p = 0.0006) and sex (p < 0.0001) were significant predictors for RBC:M, and a linear regression showed higher values and steeper decline in males: RBC:M(Males) = -0.00362 × Age + 0.60 (p = 0.01, R2 = 0.25); RBC:M(Females) = -0.00170 × Age + 0.44 (p = 0.02, R2 = 0.15). Similarly, age and sex were significant predictors for RBC:Gas but not for M:Gas. RBC:M, M:Gas and RBC:Gas were significantly lower at TLC than at FRC (plus inhaled volume), with an average 9%, 30% and 35% decrease, respectively. CONCLUSION: Expected age and sex dependence of pulmonary function concurs with 129Xe RBC:M imaging results, demonstrating that these variables must be considered when reporting Xe-MRI metrics. Xenon doses and breathing maneuvers should be controlled due to the strong dependence of Xe-MRI metrics upon lung volume.

Improving Xenon-129 lung ventilation image SNR with deep-learning based image reconstruction.

PURPOSE: To evaluate the feasibility and utility of a deep learning (DL)-based reconstruction for improving the SNR of hyperpolarized 129Xe lung ventilation MRI. METHODS: 129Xe lung ventilation MRI data acquired from patients with asthma and/or chronic obstructive pulmonary disease (COPD) were retrospectively reconstructed with a commercial DL reconstruction pipeline at five different denoising levels. Quantitative imaging metrics of lung ventilation including ventilation defect percentage (VDP) and ventilation heterogeneity index (VHI) were compared between each set of DL-reconstructed images and alternative denoising strategies including: filtering, total variation denoising and higher-order singular value decomposition. Structural similarity between the denoised and original images was assessed. In a prospective study, the feasibility of using SNR gains from DL reconstruction to allow natural-abundance xenon MRI was evaluated in healthy volunteers. RESULTS: 129Xe ventilation image SNR was improved with DL reconstruction when compared with conventionally reconstructed images. In patients with asthma and/or COPD, DL-reconstructed images exhibited a slight positive bias in ventilation defect percentage (1.3% at 75% denoising) and ventilation heterogeneity index (˜1.4) when compared with conventionally reconstructed images. Additionally, DL-reconstructed images preserved structural similarity more effectively than data denoised using alternative approaches. DL reconstruction greatly improved image SNR (greater than threefold), to a level that 129Xe ventilation imaging using natural-abundance xenon appears feasible. CONCLUSION: DL-based image reconstruction significantly improves 129Xe ventilation image SNR, preserves structural similarity, and leads to a minor bias in ventilation metrics that can be attributed to differences in the image sharpness. This tool should help facilitate cost-effective 129Xe ventilation imaging with natural-abundance xenon in the future.

A pharmacokinetic model for hyperpolarized 13C-pyruvate MRI when using metabolite-specific bSSFP sequences.

PURPOSE: Metabolite-specific balanced SSFP (MS-bSSFP) sequences are increasingly used in hyperpolarized [1-13C]Pyruvate (HP 13C) MRI studies as they improve SNR by refocusing the magnetization each TR. Currently, pharmacokinetic models used to fit conversion rate constants, kPL and kPB, and rate constant maps do not account for differences in the signal evolution of MS-bSSFP acquisitions. METHODS: In this work, a flexible MS-bSSFP model was built that can be used to fit conversion rate constants for these experiments. The model was validated in vivo using paired animal (healthy rat kidneys n = 8, transgenic adenocarcinoma of the mouse prostate n = 3) and human renal cell carcinoma (n = 3) datasets. Gradient echo (GRE) acquisitions were used with a previous GRE model to compare to the results of the proposed GRE-bSSFP model. RESULTS: Within simulations, the proposed GRE-bSSFP model fits the simulated data well, whereas a GRE model shows bias because of model mismatch. For the in vivo datasets, the estimated conversion rate constants using the proposed GRE-bSSFP model are consistent with a previous GRE model. Jointly fitting the lactate T2 with kPL resulted in less precise kPL estimates. CONCLUSION: The proposed GRE-bSSFP model provides a method to estimate conversion rate constants, kPL and kPB, for MS-bSSFP HP 13C experiments. This model may also be modified and used for other applications, for example, estimating rate constants with other hyperpolarized reagents or multi-echo bSSFP.

Particle-based MR modeling with diffusion, microstructure, and enzymatic reactions.

PURPOSE: To develop a novel particle-based in silico MR model and demonstrate applications of this model to signal mechanisms which are affected by the spatial organization of particles, including metabolic reaction kinetics, microstructural effects on diffusion, and radiofrequency (RF) refocusing effects in gradient-echo sequences. METHODS: The model was developed by integrating a forward solution of the Bloch equations with a Brownian dynamics simulator. Simulation configurations were then designed to model MR signal dynamics of interest, with a primary focus on hyperpolarized 13C MRI methods. Phantom scans and spectrophotometric assays were conducted to validate model results in vitro. RESULTS: The model accurately reproduced the reaction kinetics of enzyme-mediated conversion of pyruvate to lactate. When varying proportions of restrictive structure were added to the reaction volume, nonlinear changes in the reaction rate measured in vitro were replicated in silico. Modeling of RF refocusing effects characterized the degree of diffusion-weighted contribution from preserved residual magnetization in nonspoiled gradient-echo sequences. CONCLUSIONS: These results show accurate reproduction of a range of MR signal mechanisms, establishing the model's capability to investigate the multifactorial signal dynamics such as those underlying hyperpolarized 13C MRI data.

Hyperpolarized [1-13 C] pyruvate MRSI to detect metabolic changes in liver in a methionine and choline-deficient diet rat model of fatty liver disease.

PURPOSE: Nonalcoholic fatty liver disease is an important cause of chronic liver disease. There are limited methods for monitoring metabolic changes during progression to steatohepatitis. Hyperpolarized 13 C MRSI (HP 13 C MRSI) was used to measure metabolic changes in a rodent model of fatty liver disease. METHODS: Fifteen Wistar rats were placed on a methionine- and choline-deficient (MCD) diet for 1-18 weeks. HP 13 C MRSI, T2 -weighted imaging, and fat-fraction measurements were obtained at 3 T. Serum aspartate aminotransaminase, alanine aminotransaminase, and triglycerides were measured. Animals were sacrificed for histology and measurement of tissue lactate dehydrogenase (LDH) activity. RESULTS: Animals lost significant weight (13.6% ± 2.34%), an expected characteristic of the MCD diet. Steatosis, inflammation, and mild fibrosis were observed. Liver fat fraction was 31.7% ± 4.5% after 4 weeks and 22.2% ± 4.3% after 9 weeks. Lactate-to-pyruvate and alanine-to-pyruvate ratios decreased significantly over the study course; were negatively correlated with aspartate aminotransaminase and alanine aminotransaminase (r = -[0.39-0.61]); and were positively correlated with triglycerides (r = 0.59-0.60). Despite observed decreases in hyperpolarized lactate signal, LDH activity increased by a factor of 3 in MCD diet-fed animals. Observed decreases in lactate and alanine hyperpolarized signals on the MCD diet stand in contrast to other studies of liver injury, where lactate and alanine increased. Observed hyperpolarized metabolite changes were not explained by alterations in LDH activity, suggesting that changes may reflect co-factor depletion known to occur as a result of oxidative stress in the MCD diet. CONCLUSION: HP 13 C MRSI can noninvasively measure metabolic changes in the MCD model of chronic liver disease.

A vendor-neutral EPI sequence for hyperpolarized 13C MRI.

PURPOSE: To develop a flexible, vendor-neutral EPI sequence for hyperpolarized 13C metabolic imaging. METHODS: An open-source EPI sequence consisting of a metabolite-specific spectral-spatial RF excitation pulse and a customizable EPI readout was created using the Pulseq framework. To explore the flexibility of our sequence, we tested several versions of the sequence including a symmetric 3D readout with different spatial resolutions for each metabolite (1.0 cm3 and 1.5 cm3). A multichamber phantom constructed with a Shepp-Logan geometry, containing two chambers filled with either natural abundance 13C compounds or hyperpolarized (HP) [1-13C]pyruvate, was used to test each sequence. For experiments involving HP [1-13C]pyruvate, a single chamber was prefilled with nicotinamide adenine dinucleotide hydride and lactate dehydrogenase to facilitate the conversion of [1-13C]pyruvate to [1-13C]lactate. All experiments were performed on a Siemens Prisma 3T scanner. RESULTS: All the sequence variations localized natural-abundance 13C ethylene glycol and methanol to the appropriate compartment of the multichamber phantom. [1-13C]pyruvate was detectable in both chambers following the injection of HP [1-13C]pyruvate, whereas [1-13C]lactate was only found in the chamber containing nicotinamide adenine dinucleotide hydride and lactate dehydrogenase. The conversion rate from [1-13C]pyruvate to [1-13C]lactate (kPL) was 0.01 s-1 (95% confidence interval [0.00, 0.02]). CONCLUSION: We have developed and tested a vendor-neutral EPI sequence for imaging HP 13C agents. We have made all of our sequence creation and image reconstruction code freely available online for other investigators to use.

Distinguishing metabolic signals of liver tumors from surrounding liver cells using hyperpolarized 13 C MRI and gadoxetate.

PURPOSE: To use the hepatocyte-specific gadolinium-based contrast agent gadoxetate combined with hyperpolarized (HP) [1-13 C]pyruvate MRI to selectively suppress metabolic signals from normal hepatocytes while preserving the signals arising from tumors. METHODS: Simulations were performed to determine the expected changes in HP 13 C MR signal in liver and tumor under the influence of gadoxetate. CC531 colon cancer cells were implanted into the livers of five Wag/Rij rats. Liver and tumor metabolism were imaged at 3 T using HP [1-13 C] pyruvate chemical shift imaging before and 15 min after injection of gadoxetate. Area under the curve for pyruvate and lactate were measured from voxels containing at least 75% of normal-appearing liver or tumor. RESULTS: Numerical simulations predicted a 36% decrease in lactate-to-pyruvate (L/P) ratio in liver and 16% decrease in tumor. In vivo, baseline L/P ratio was 0.44 ± 0.25 in tumors versus 0.21 ± 0.08 in liver (p = 0.09). Following administration of gadoxetate, mean L/P ratio decreased by an average of 0.11 ± 0.06 (p < 0.01) in normal-appearing liver. In tumors, mean L/P ratio post-gadoxetate did not show a statistically significant change from baseline. Compared to baseline levels, the relative decrease in L/P ratio was significantly greater in liver than in tumors (-0.52 ± 0.16 vs. -0.19 ± 0.25, p < 0.05). CONCLUSIONS: The intracellular hepatobiliary contrast agent showed a greater effect suppressing HP 13 C MRI metabolic signals (through T1 shortening) in normal-appearing liver when compared to tumors. The combined use of HP MRI with selective gadolinium contrast agents may allow more selective imaging in HP 13 C MRI.

Functional activation of pyruvate dehydrogenase in human brain using hyperpolarized [1-13 C]pyruvate.

PURPOSE: Pyruvate, produced from either glucose, glycogen, or lactate, is the dominant precursor of cerebral oxidative metabolism. Pyruvate dehydrogenase (PDH) flux is a direct measure of cerebral mitochondrial function and metabolism. Detection of [13 C]bicarbonate in the brain from hyperpolarized [1-13 C]pyruvate using carbon-13 (13 C) MRI provides a unique opportunity for assessing PDH flux in vivo. This study is to assess changes in cerebral PDH flux in response to visual stimuli using in vivo 13 C MRS with hyperpolarized [1-13 C]pyruvate. METHODS: From seven sedentary adults in good general health, time-resolved [13 C]bicarbonate production was measured in the brain using 90° flip angles with minimal perturbation of its precursors, [1-13 C]pyruvate and [1-13 C]lactate, to test the hypothesis that the appearance of [13 C]bicarbonate signals in the brain reflects the metabolic changes associated with neuronal activation. With a separate group of healthy participants (n = 3), the likelihood of the bolus-injected [1-13 C]pyruvate being converted to [1-13 C]lactate prior to decarboxylation was investigated by measuring [13 C]bicarbonate production with and without [1-13 C]lactate saturation. RESULTS: In the course of visual stimulation, the measured [13 C]bicarbonate signal normalized to the total 13 C signal in the visual cortex increased by 17.1% ± 15.9% (p = 0.017), whereas no significant change was detected in [1-13 C]lactate. Proton BOLD fMRI confirmed the regional activation in the visual cortex with the stimuli. Lactate saturation decreased bicarbonate-to-pyruvate ratio by 44.4% ± 9.3% (p < 0.01). CONCLUSION: We demonstrated the utility of 13 C MRS with hyperpolarized [1-13 C]pyruvate for assessing the activation of cerebral PDH flux via the detection of [13 C]bicarbonate production.

First in-human evaluation of [1-13C]pyruvate in D2O for hyperpolarized MRI of the brain: A safety and feasibility study.

PURPOSE: To investigate the safety and value of hyperpolarized (HP) MRI of [1-13C]pyruvate in healthy volunteers using deuterium oxide (D2O) as a solvent. METHODS: Healthy volunteers (n = 5), were injected with HP [1-13C]pyruvate dissolved in D2O and imaged with a metabolite-specific 3D dual-echo dynamic EPI sequence at 3T at one site (Site 1). Volunteers were monitored following the procedure to assess safety. Image characteristics, including SNR, were compared to data acquired in a separate cohort using water as a solvent (n = 5) at another site (Site 2). The apparent spin-lattice relaxation time (T1) of [1-13C]pyruvate was determined both in vitro and in vivo from a mono-exponential fit to the image intensity at each time point of our dynamic data. RESULTS: All volunteers completed the study safely and reported no adverse effects. The use of D2O increased the T1 of [1-13C]pyruvate from 66.5 ± 1.6 s to 92.1 ± 5.1 s in vitro, which resulted in an increase in signal by a factor of 1.46 ± 0.03 at the time of injection (90 s after dissolution). The use of D2O also increased the apparent relaxation time of [1-13C]pyruvate by a factor of 1.4 ± 0.2 in vivo. After adjusting for inter-site SNR differences, the use of D2O was shown to increase image SNR by a factor of 2.6 ± 0.2 in humans. CONCLUSIONS: HP [1-13C]pyruvate in D2O is safe for human imaging and provides an increase in T1 and SNR that may improve image quality.

Current methods for hyperpolarized [1-13C]pyruvate MRI human studies.

MRI with hyperpolarized (HP) 13C agents, also known as HP 13C MRI, can measure processes such as localized metabolism that is altered in numerous cancers, liver, heart, kidney diseases, and more. It has been translated into human studies during the past 10 years, with recent rapid growth in studies largely based on increasing availability of HP agent preparation methods suitable for use in humans. This paper aims to capture the current successful practices for HP MRI human studies with [1-13C]pyruvate-by far the most commonly used agent, which sits at a key metabolic junction in glycolysis. The paper is divided into four major topic areas: (1) HP 13C-pyruvate preparation; (2) MRI system setup and calibrations; (3) data acquisition and image reconstruction; and (4) data analysis and quantification. In each area, we identified the key components for a successful study, summarized both published studies and current practices, and discuss evidence gaps, strengths, and limitations. This paper is the output of the "HP 13C MRI Consensus Group" as well as the ISMRM Hyperpolarized Media MR and Hyperpolarized Methods and Equipment study groups. It further aims to provide a comprehensive reference for future consensus, building as the field continues to advance human studies with this metabolic imaging modality.

3D quantitative myocardial perfusion imaging with hyperpolarized HP001(bis-1,1-(hydroxymethyl)-[1-13C]cyclopropane-d8): Application of gradient echo and balanced SSFP sequences.

PURPOSE: This study aims to show the viability of conducting three-dimensional (3D) myocardial perfusion quantification covering the entire heart using both GRE and bSSFP sequences with hyperpolarized HP001. METHODS: A GRE sequence and a bSSFP sequence, both with a stack-of-spirals readout, were designed and applied to three pigs. The images were reconstructed using 13 $$ {}^{13} $$ C coil sensitivity maps measured in a phantom experiment. Perfusion was quantified using a constrained decomposition method, and the estimated rest/stress perfusion values from 13 $$ {}^{13} $$ C GRE/bSSFP and Dynamic contrast-enhanced MRI (DCE-MRI) were individually analyzed through histograms and the mean perfusion values were compared with reference values obtained from PET( 15 $$ {}^{15} $$ O-water). The Myocardial Perfusion Reserve Index (MPRI) was estimated for 13 $$ {}^{13} $$ C GRE/bSSFP and DCE-MRI and compared with the reference values. RESULTS: Perfusion values, estimated by both DCE and 13 $$ {}^{13} $$ C MRI, were found to be lower than reference values. However, DCE-MRI's estimated perfusion values were closer to the reference values than those obtained from 13 $$ {}^{13} $$ C MRI. In the case of MPRI estimation, the 13 $$ {}^{13} $$ C estimated MPRI values (GRE/bSSFP: 2.3/2.0) more closely align with the literature value (around 3) than the DCE estimated MPRI value (1.6). CONCLUSION: This study demonstrated the feasibility of 3D whole-heart myocardial perfusion quantification using hyperpolarized HP001 with both GRE and bSSFP sequences. The 13 $$ {}^{13} $$ C perfusion measurements underestimated perfusion values compared to the 15 $$ {}^{15} $$ O PET literature value, while the 13 $$ {}^{13} $$ C estimated MPRI value aligned better with the literature. This preliminary result indicates 13 $$ {}^{13} $$ C imaging may more accurately estimate MPRI values compared to DCE-MRI.

Dynamic T 2 * relaxometry of hyperpolarized [1-13 C]pyruvate MRI in the human brain and kidneys.

PURPOSE: This study aimed to quantify T 2 * $$ {T}_2^{\ast } $$ for hyperpolarized [1-13 C]pyruvate and metabolites in the healthy human brain and renal cell carcinoma (RCC) patients at 3 T. METHODS: Dynamic T 2 * $$ {T}_2^{\ast } $$ values were measured with a metabolite-specific multi-echo spiral sequence. The dynamic T 2 * $$ {T}_2^{\ast } $$ of [1-13 C]pyruvate, [1-13 C]lactate, and 13 C-bicarbonate was estimated in regions of interest in the whole brain, sinus vein, gray matter, and white matter in healthy volunteers, as well as in kidney tumors and the contralateral healthy kidneys in a separate group of RCC patients. T 2 * $$ {T}_2^{\ast } $$ was fit using a mono-exponential function; and metabolism was quantified using pyruvate-to-lactate conversion rate maps and lactate-to-pyruvate ratio maps, which were compared with and without an estimated T 2 * $$ {T}_2^{\ast } $$ correction. RESULTS: The T 2 * $$ {T}_2^{\ast } $$ of pyruvate was shown to vary during the acquisition, whereas the T 2 * $$ {T}_2^{\ast } $$ of lactate and bicarbonate were relatively constant through time and across the organs studied. The T 2 * $$ {T}_2^{\ast } $$ of lactate was similar in gray matter (29.75 ± 1.04 ms), white matter (32.89 ± 0.9 ms), healthy kidney (34.61 ± 4.07 ms), and kidney tumor (33.01 ± 2.31 ms); and the T 2 * $$ {T}_2^{\ast } $$ of bicarbonate was different between whole-brain (108.17 ± 14.05 ms) and healthy kidney (58.45 ± 6.63 ms). The T 2 * $$ {T}_2^{\ast } $$ of pyruvate had similar trends in both brain and RCC studies, reducing from 75.56 ± 2.23 ms to 22.24 ± 1.24 ms in the brain and reducing from 122.72 ± 9.86 ms to 57.38 ± 7.65 ms in the kidneys. CONCLUSION: Multi-echo dynamic imaging can quantify T 2 * $$ {T}_2^{\ast } $$ and metabolism in a single integrated acquisition. Clear differences were observed in the T 2 * $$ {T}_2^{\ast } $$ of metabolites and in their behavior throughout the timecourse.

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.

A decay-modeled compressed sensing reconstruction approach for non-Cartesian hyperpolarized 129Xe MRI.

PURPOSE: Hyperpolarized 129Xe MRI benefits from non-Cartesian acquisitions that sample k-space efficiently and rapidly. However, their reconstructions are complex and burdened by decay processes unique to hyperpolarized gas. Currently used gridded reconstructions are prone to artifacts caused by magnetization decay and are ill-suited for undersampling. We present a compressed sensing (CS) reconstruction approach that incorporates magnetization decay in the forward model, thereby producing images with increased sharpness and contrast, even in undersampled data. METHODS: Radio-frequency, T1, and T 2 * $$ {\mathrm{T}}_2^{\ast } $$ decay processes were incorporated into the forward model and solved using iterative methods including CS. The decay-modeled reconstruction was validated in simulations and then tested in 2D/3D-spiral ventilation and 3D-radial gas-exchange MRI. Quantitative metrics including apparent-SNR and sharpness were compared between gridded, CS, and twofold undersampled CS reconstructions. Observations were validated in gas-exchange data collected from 15 healthy and 25 post-hematopoietic-stem-cell-transplant participants. RESULTS: CS reconstructions in simulations yielded images with threefold increases in accuracy. CS increased sharpness and contrast for ventilation in vivo imaging and showed greater accuracy for undersampled acquisitions. CS improved gas-exchange imaging, particularly in the dissolved-phase where apparent-SNR improved, and structure was made discernable. Finally, CS showed repeatability in important global gas-exchange metrics including median dissolved-gas signal ratio and median angle between real/imaginary components. CONCLUSION: A non-Cartesian CS reconstruction approach that incorporates hyperpolarized 129Xe decay processes is presented. This approach enables improved image sharpness, contrast, and overall image quality in addition to up-to threefold undersampling. This contribution benefits all hyperpolarized gas MRI through improved accuracy and decreased scan durations.

Compressed sensing reconstruction for high-SNR, rapid dissolved 129Xe gas exchange MRI.

PURPOSE: Three-dimensional hyperpolarized 129Xe gas exchange imaging suffers from low SNR and long breath-holds, which could be improved using compressed sensing (CS). The purpose of this work was to assess whether gas exchange ratio maps are quantitatively preserved in CS-accelerated dissolved-phase 129Xe imaging and to investigate the feasibility of CS-dissolved 129Xe imaging with reduced-cost natural abundance (NA) xenon. METHODS: 129Xe gas exchange imaging was performed at 1.5 T with a multi-echo spectroscopic imaging sequence. A CS reconstruction with an acceleration factor of 2 was compared retrospectively with conventional gridding reconstruction in a cohort of 16 healthy volunteers, 5 chronic obstructive pulmonary disease patients, and 23 patients who were hospitalized following COVID-19 infection. Metrics of comparison included normalized mean absolute error, mean gas exchange ratio, and red blood cell (RBC) image SNR. Dissolved 129Xe CS imaging with NA xenon was assessed in 4 healthy volunteers. RESULTS: CS reconstruction enabled acquisition time to be halved, and it reduced background noise. Median RBC SNR increased from 6 (2-18) to 11 (2-100) with CS, and there was strong agreement between CS and gridding mean ratio map values (R2 = 0.99). Image fidelity was maintained for gridding RBC SNR > 5, but below this, normalized mean absolute error increased nonlinearly with decreasing SNR. CS increased the mean SNR of NA 129Xe images 3-fold. CONCLUSION: CS reconstruction of dissolved 129Xe imaging improved image quality with decreased scan time, while preserving key gas exchange metrics. This will benefit patients with breathlessness and/or low gas transfer and shows promise for NA-dissolved 129Xe imaging.

Comparison of 3D UTE free-breathing lung MRI with hyperpolarized 129Xe MRI in pediatric cystic fibrosis.

PURPOSE: To compare phase-resolved functional lung (PREFUL) regional ventilation derived from a free breathing 3D UTE radial MRI acquisition to hyperpolarized 129Xe-MRI (Xe-MRI), conventional 2D multi-slice PREFUL MRI, and pulmonary function tests in pediatric cystic fibrosis (CF) lung disease. METHODS: Free-breathing 3D UTE and 2D multi-slice 1H MRI as well as Xe-MRI were acquired in 12 stable pediatric CF patients. Using PREFUL, regional ventilation (RVent) maps were calculated from the free-breathing data. Ventilation defect percentage (VDP) was determined from 3D and 2D RVent maps (2D VDPRVent and 3D VDPRVent, respectively) and Xe-MRI ventilation (VDPXe). VDP was calculated for the whole lung and for eight regions based on left/right, anterior/posterior, and superior/inferior divisions of the lung. Global and regional VDP was compared between the three methods using Bland-Altman analysis, linear mixed model-based correlation, and one-way analysis of variance and multiple comparisons tests. RESULTS: Global 3D VDPRVent, VDPXe, and 2D VDPRVent were all strongly correlated (all R2 > 0.62, p < 0.0001) and showed minimal, non-significant bias (all <2%, p > 0.05). Three dimensional and 2D VDPRVent significantly correlated to VDPXe in most of the separate lung regions (R2 = 0.18-0.74, p < 0.04), but showed lower inter-agreement. The superior/anterior lung regions showed the least agreement between all three methods (all p > 0.12). CONCLUSION: Absolute VDP assessed by 3D UTE PREFUL MRI showed good global agreement with Xe-MRI and 2D multi-slice PREFUL MRI in pediatric CF lung disease. Therefore, 3D UTE PREFUL MRI offers a sensitive and potentially more accessible alternative to Xe-MRI for regional volumetric evaluation of ventilation.

Hyperpolarized 13C and 31P MRS detects differences in cardiac energetics, metabolism, and function in obesity, and responses following treatment.

Obesity is associated with important changes in cardiac energetics and function, and an increased risk of adverse cardiovascular outcomes. Multi-nuclear MRS and MRI techniques have the potential to provide a comprehensive non-invasive assessment of cardiac metabolic perturbation in obesity. A rat model of obesity was created by high-fat diet feeding. This model was characterized using in vivo hyperpolarized [1-13C]pyruvate and [2-13C]pyruvate MRS, echocardiography and perfused heart 31P MRS. Two groups of obese rats were subsequently treated with either caloric restriction or the glucagon-like peptide-1 analogue/agonist liraglutide, prior to reassessment. The model recapitulated cardiovascular consequences of human obesity, including mild left ventricular hypertrophy, and diastolic, but not systolic, dysfunction. Hyperpolarized 13C and 31P MRS demonstrated that obesity was associated with reduced myocardial pyruvate dehydrogenase flux, altered cardiac tricarboxylic acid (TCA) cycle metabolism, and impaired myocardial energetic status (lower phosphocreatine to adenosine triphosphate ratio and impaired cardiac ΔG~ATP). Both caloric restriction and liraglutide treatment were associated with normalization of metabolic changes, alongside improvement in cardiac diastolic function. In this model of obesity, hyperpolarized 13C and 31P MRS demonstrated abnormalities in cardiac metabolism at multiple levels, including myocardial substrate selection, TCA cycle, and high-energy phosphorus metabolism. Metabolic changes were linked with impairment of diastolic function and were reversed in concert following either caloric restriction or liraglutide treatment. With hyperpolarized 13C and 31P techniques now available for human use, the findings support a role for multi-nuclear MRS in the development of new therapies for obesity.

Sampling pattern discrepancy in the application of compressed sensing hyperpolarized xenon-129 lung MRI.

Although hyperpolarized (HP) 129Xe ventilation MRI can be carried out within a breath hold, it is still challenging for many sick patients. Compressed sensing (CS) is a viable alternative to accelerate this approach. However, undersampled images with identical sampling ratios differ from one another. Twenty subjects (n = 10 healthy and n = 10 patients with asthma) were scanned using a GE MR750 3 T scanner, acquiring fully sampled 2D multi-slice HP 129Xe lung ventilation images (10 s breath hold, 128 × 80 (FE × PE-frequency encoding × phase encoding) and 16 slices). Using fully sampled data, 500 variable-density Cartesian random undersampling patterns were generated, each at eight different sampling ratios from 10% to 80%. The parallel imaging and compressed sensing (PICS) command from BART was employed to reconstruct undersampled data. The signal to noise ratio (SNR), structural similarity index measurement (SSIM) and sidelobe to peak ratio of each were subsequently compared. There was a high degree of variation in both SNR and SSIM results from each of the 500 masks of each sampling rate. As the undersampling increases, there is more variation in the quantifying metrics, for both healthy and asthmatic individuals. Our study shows that random undersampling poses a significant challenge when applied at sampling ratios less than 60%, despite fulfilling CS's incoherency criteria. Such low sampling ratios will result in a large variety of undersampling patterns. Therefore, skipped segments of k-space cannot be allowed to happen randomly at low sampling rates. By optimizing the sampling pattern, CS will reach its full potential and be able to be applied to a highly undersampled 129Xe lung dataset.