K-t PCA accelerated in-plane balanced steady-state free precession phase-contrast (PC-SSFP) for all-in-one diastolic function evaluation.

PURPOSE: Diastolic function evaluation requires estimates of early and late diastolic mitral filling velocities (E and A) and of mitral annulus tissue velocity (e'). We aimed to develop an MRI method for simultaneous all-in-one diastolic function evaluation in a single scan by generating a 2D phase-contrast (PC) sequence with balanced steady-state free precession (bSSFP) contrast (PC-SSFP). E and A could then be measured with PC, and e' estimated by valve tracking on the magnitude images, using an established deep learning framework. METHODS: Our PC-SSFP used in-plane flow-encoding, with zeroth and first moment nulling over each TR. For further acceleration, different k-t principal component analysis (PCA) methods were investigated with both retrospective and prospective undersampling. PC-SSFP was compared to separate balanced SSFP cine and PC-gradient echo acquisitions in phantoms and in 10 healthy subjects. RESULTS: Phantom experiments showed that PC-SSFP measured accurate velocities compared to PC-gradient echo (r = 0.98 for a range of pixel-wise velocities -80 cm/s to 80 cm/s). In subjects, PC-SSFP generated high SNR and myocardium-blood contrast, and excellent agreement for E (limits of agreement [LOA] 0.8 ± 2.4 cm/s, r = 0.98), A (LOA 2.5 ± 4.1 cm/s, r = 0.97), and e' (LOA 0.3 ± 2.6 cm/s, r = 1.00), versus the standard methods. The best k-t PCA approach processed the complex difference data and substituted in raw k-space data. With prospective k-t PCA acceleration, higher frame rates were achieved (50 vs. 25 frames per second without k-t PCA), yielding a 13% higher e'. CONCLUSION: The proposed PC-SSFP method achieved all-in-one diastolic function evaluation.

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.

Robust T2 estimation with balanced steady state free precession.

PURPOSE: To develop a novel signal representation for balanced steady state free precession (bSSFP) displaying its T2 independence on B1 and on magnetization transfer (MT) effects. METHODS: A signal model for bSSFP is developed that shows only an explicit dependence (up to a scaling factor) on E2 (and, therefore, T2) and a novel parameter c (with implicit dependence on the flip angle and E1). Moreover, it is shown that MT effects, entering the bSSFP signal via a binary spin bath model, can be captured by a redefinition of T1 and, therefore, leading to modification of E1, resulting in the same signal model. Various sets of phase-cycled bSSFP brain scans (different flip angles, different TR, different RF pulse durations, and different number of phase cycles) were recorded at 3 T. The parameters T2 (E2) and c were estimated using a variable projection (VARPRO) method and Monte-Carlo simulations were performed to assess T2 estimation precision. RESULTS: Initial experiments confirmed the expected independence of T2 on various protocol settings, such as TR, the flip angle, B1 field inhomogeneity, and the RF pulse duration. Any variation (within the explored range) appears to directly affect the estimation of the parameter c only-in agreement with theory. CONCLUSION: BSSFP theory predicts an extraordinary feature that all MT and B1-related variational aspects do not enter T2 estimation, making it a potentially robust methodology for T2 quantification, pending validation against existing standards.

Submillimeter balanced SSFP BOLD-functional MRI accelerated with 3D stack-of-spirals at 9.4 T.

PURPOSE: This work aims to improve the speed of balanced SSFP (bSSFP) acquisition with segmented 3D stack-of-spirals for functional brain studies at ultrahigh field. METHODS: Functional experiments were performed with an accelerated 3D stack-of-spirals sequence with water excitation for fat suppression. The resulting data were reconstructed using an iterative algorithm with corrections for system imperfections such as trajectory deviations and B0 inhomogeneity. In the first set of experiments, we evaluated the signal change and stability with respect to echo and TR for a full-field checkerboard stimulus. To demonstrate the high spatio-temporal resolution of the developed method, the results of three optimized protocols at submillimeter resolution (0.6-mm isotropic and 0.8-mm isotropic) and at 1.2 mm isotropic resolution for whole-brain coverage were shown. RESULTS: Water excitation and the model-based iterative reconstruction improved image quality. The BOLD-related signal changes increased with longer TE and longer TR. We observed an increase in thermal noise performance at lower TE and higher TR. However, signal stability deteriorates at higher TE and TR. Therefore, optimized protocols used shorter TE and moderately long TR to maximize the sensitivity and speed. Reproducible activations were detected along the gray-matter gyri in the submillimeter protocols with a median signal change of approximately 4% across subjects. CONCLUSIONS: Three-dimensional stack-of-spirals enables passband balanced SSFP functional imaging at a much higher spatial and temporal scale, compared with conventional spoiled gradient-echo train sequences.

Intra-scan RF power amplifier drift correction.

PURPOSE: The drift in radiofrequency (RF) power amplifiers (RFPAs) is assessed and several contributing factors are investigated. Two approaches for prospective correction of drift are proposed and their effectiveness is evaluated. METHODS: RFPA drift assessment encompasses both intra-pulse and inter-pulse drift analyses. Scan protocols with varying flip angle (FA), RF length, and pulse repetition time (TR) are used to gauge the influence of these parameters on drift. Directional couplers (DICOs) monitor the forward waveforms of the RFPA outputs. DICOs data is stored for evaluation, allowing calculation of correction factors to adjust RFPAs' transmit voltage. Two correction methods, predictive and run-time, are employed: predictive correction necessitates a calibration scan, while run-time correction calculates factors during the ongoing scan. RESULTS: RFPA drift is indeed influenced by the RF duty-cycle, and in the cases examined with a maximum duty-cycle of 66%, the potential drift is approximately 41% or 15%, depending on the specific RFPA revision. Notably, in low transmit voltage scenarios, FA has minimal impact on RFPA drift. The application of predictive and run-time drift correction techniques effectively reduces the average drift from 10.0% to less than 1%, resulting in enhanced MR signal stability. CONCLUSION: Utilizing DICO recordings and implementing a feedback mechanism enable the prospective correction of RFPA drift. Having a calibration scan, predictive correction can be utilized with fewer complexity; for enhanced performance, a run-time approach can be employed.

Vessel-specific quantification of cerebral venous oxygenation with velocity-encoding preparation and rapid acquisition.

PURPOSE: Non-invasive measurement of cerebral venous oxygenation (Yv) is of critical importance in brain diseases. The present work proposed a fast method to quantify regional Yv map for both large and small veins. METHODS: A new sequence was developed, referred to as TRU-VERA (T2 relaxation under velocity encoding and rapid acquisition, which isolates blood spins from static tissue with velocity-encoding preparation, modulates the T2 weighting of venous signal with T2-preparation and utilizes a bSSFP readout to achieve fast acquisition with high resolution. The sequence was first optimized to achieve best sensitivity for both large and small veins, and then validated with TRUST (T2 relaxation under spin tagging), TRUPC (T2 relaxation under phase contrast), and accelerated TRUPC MRI. Regional difference of Yv was evaluated, and test-retest reproducibility was examined. RESULTS: Optimal Venc was determined to be 3 cm/s, while recovery time and balanced SSFP flip angle within reasonable range had minimal effect on SNR efficiency. Venous T2 measured with TRU-VERA was highly correlated with T2 from TRUST (R2 = 0.90), and a conversion equation was established for further calibration to Yv. TRU-VERA sequences showed consistent Yv estimation with TRUPC (R2 = 0.64) and accelerated TRUPC (R2 = 0.79). Coefficient of variation was 0.84% for large veins and 2.49% for small veins, suggesting an excellent test-retest reproducibility. CONCLUSION: The proposed TRU-VERA sequence is a promising method for vessel-specific oxygenation assessment.

Getting the phase consistent: The importance of phase description in balanced steady-state free precession MRI of multi-compartment systems.

PURPOSE: Determine the correct mathematical phase description for balanced steady-state free precession (bSSFP) signals in multi-compartment systems. THEORY AND METHODS: Based on published bSSFP signal models, different phase descriptions can be formulated: one predicting the presence and the other predicting the absence of destructive interference effects in multi-compartment systems. Numerical simulations of bSSFP signals of water and acetone were performed to evaluate the predictions of these different phase descriptions. For experimental validation, bSSFP profiles were measured at 3T using phase-cycled bSSFP acquisitions performed in a phantom containing mixtures of water and acetone, which replicates a system with two signal components. Localized single voxel MRS was performed at 7T to determine the relative chemical shift of the acetone-water mixtures. RESULTS: Based on the choice of phase description, the simulated bSSFP profiles of water-acetone mixtures varied significantly, either displaying or lacking destructive interference effects, as predicted theoretically. In phantom experiments, destructive interference was consistently observed in the measured bSSFP profiles of water-acetone mixtures, supporting the theoretical description that predicts such interference effects. The connection between the choice of phase description and predicted observation enables unambiguous experimental identification of the correct phase description for multi-compartment bSSFP profiles, which is consistent with the Bloch equations. CONCLUSION: The study emphasizes that consistent phase descriptions are crucial for accurately describing multi-compartment bSSFP signals, as incorrect phase descriptions result in erroneous predictions.

Sequence comparison of spoiled gradient echo and balanced steady-state free precession for pulmonary free-breathing proton MRI in patients and healthy volunteers: Correspondence, repeatability, and validation with dynamic contrast-enhanced MRI.

Phase-resolved functional lung (PREFUL) MRI is a proton-based, contrast agent-free technique derived from the Fourier decomposition approach to measure regional ventilation and perfusion dynamics during free-breathing. Besides the necessity of extensive PREFUL postprocessing, the utilized MRI sequence must fulfill specific requirements. This study investigates the impact of sequence selection on PREFUL-MRI-derived functional parameters by comparing the standard spoiled gradient echo (SPGRE) sequence with a lung-optimized balanced steady-state free precession (bSSFP) sequence, thereby facilitating PREFULs clinical application in pulmonary disease assessment. This study comprised a prospective dataset of healthy volunteers and a retrospective dataset of patients with suspected chronic thromboembolic pulmonary hypertension. Both cohorts underwent PREFUL-MRI with both sequences to assess the correspondence of PREFUL ventilation and perfusion parameters (A). Additionally, healthy subjects were scanned a second time to evaluate repeatability (B), whereas patients received dynamic contrast-enhanced (DCE)-MRI, considered the perfusion gold standard for comparison with PREFUL-MRI (C). Signal-to-noise ratio (SNR), calculated from the unprocessed images, was compared alongside median differences of PREFUL-MRI-derived parameters using a paired Wilcoxon signed rank test. Further evaluations included calculation of the Pearson correlation, intraclass-correlation coefficient for repeatability assessment, and spatial overlap (SO) for regional comparison of PREFUL-MRI and DCE-MRI. bSSFP showed a clear SNR advantage over SPGRE (median: 23 vs. 9, p < 0.001). (A) Despite significant differences, parameter values were strongly correlated (r ≥ 0.75). After thresholding, binary maps showed high healthy overlap across both cohorts (SOHealthy > 86%) and high defect overlap in the patient cohort (SODefect ≥ 48%). (B) bSSFP demonstrated slightly higher repeatability across most parameters. (C) Both sequences demonstrated comparable correspondence to DCE-MRI, with SPGRE excelling in absolute quantification and bSSFP in spatial agreement. Although bSSFP showed superior SNR results, both sequences displayed spatial defect concordance and highly correlated PREFUL parameters with deviations regarding repeatability and alignment with DCE-MRI.

Unfolding coil localized errors from an imperfect slice profile using a structured autocalibration matrix: An application to reduce outflow effects in cine bSSFP imaging.

PURPOSE: Balanced steady-state free precession (bSSFP) imaging is susceptible to outflow effects where excited spins leaving the slice as part of the blood stream are misprojected back onto the imaging plane. Previous work proposed using slice-encoding steps to localize these outflow effects from corrupting the target slice, at the expense of prolonged scan time. This present study extends this idea by proposing a means of significantly reducing most of the outflowing signal from the imaged slice using a coil localization method that acquires a slice-encoded calibration scan in addition to the 2D data, without being nearly as time-demanding as our previous method. This coil localization method is titled UNfolding Coil Localized Errors from an imperfect slice profile using a Structured Autocalibration Matrix (UNCLE SAM). METHODS: Retrospective and prospective evaluations were carried out. Both featured a 2D acquisition and a separate slice-encoded calibration of the center in-plane k $$ k $$ -space lines across all desired slice-encoding steps. RESULTS: Retrospective results featured a slice-by-slice comparison of the slice-encoded images with UNCLE SAM. UNCLE SAM's subtraction from the slice-encoded image was compared with a subtraction from the flow-corrupted 2D image, to demonstrate UNCLE SAM's capability to unfold outflowing spins. UNCLE SAM's comparison with slice encoding showed that UNCLE SAM was able to unfold up to 74% of what slice encoding achieved. Prospective results showed significant reduction in outflow effects with only a marginal increase in scan time from the 2D acquisition. CONCLUSIONS: We developed a method that effectively unfolds most outflowing spins from corrupting the target slice and does not require the explicit use of slice-encoding gradients. This development offers a method to reduce most outflow effects from the target slice within a clinically feasible scan duration compared with the fully sampled slice-encoding technique.

Incorporation of view sharing and KWIC filtering into GRASP-Pro improves spatial resolution of single-shot, multi-TI, late gadolinium enhancement MRI.

While single-shot late gadolinium enhancement (LGE) is useful for imaging patients with arrhythmia and/or dyspnea, it produces low spatial resolution. One approach to improve spatial resolution is to accelerate data acquisition using compressed sensing (CS). Our previous work described a single-shot, multi-inversion time (TI) LGE pulse sequence using radial k-space sampling and CS, but over-regularization resulted in significant image blurring that muted the benefits of data acceleration. The purpose of the present study was to improve the spatial resolution of the single-shot, multi-TI LGE pulse sequence by incorporating view sharing (VS) and k-space weighted contrast (KWIC) filtering into a GRASP-Pro reconstruction. In 24 patients (mean age = 61 ± 16 years; 9/15 females/males), we compared the performance of our improved multi-TI LGE and standard multi-TI LGE, where clinical standard LGE was used as a reference. Two clinical raters independently graded multi-TI images and clinical LGE images visually on a five-point Likert scale (1, nondiagnostic; 3, clinically acceptable; 5, best) for three categories: the conspicuity of myocardium or scar, artifact, and noise. The summed visual score (SVS) was defined as the sum of the three scores. Myocardial scar volume was quantified using the full-width at half-maximum method. The SVS was not significantly different between clinical breath-holding LGE (median 13.5, IQR 1.3) and multi-TI LGE (median 12.5, IQR 1.6) (P = 0.068). The myocardial scar volumes measured from clinical standard LGE and multi-TI LGE were strongly correlated (coefficient of determination, R2 = 0.99) and in good agreement (mean difference = 0.11%, lower limit of the agreement = -2.13%, upper limit of the agreement = 2.34%). The inter-rater agreement in myocardial scar volume quantification was strong (intraclass correlation coefficient = 0.79). The incorporation of VS and KWIC into GRASP-Pro improved spatial resolution. Our improved 25-fold accelerated, single-shot LGE sequence produces clinically acceptable image quality, multi-TI reconstruction, and accurate myocardial scar volume quantification.

SIMPLEX: Multiple phase-cycled bSSFP quantitative magnetization transfer imaging with physic-guided simulation learning of neural network.

Most quantitative magnetization transfer (qMT) imaging methods require acquiring additional quantitative maps (such as T1) for data fitting. A method based on multiple phase-cycled bSSFP was recently proposed to enable high-resolution 3D qMT imaging based on least square fitting without any extra acquisition, and thus has high potential for simplifying the qMT procedure. However, the quantification of qMT parameters with this method was suboptimal, limiting its potential for clinical application despite its simpler protocol and higher spatial resolution. To improve the fitting of qMT data obtained with multiple phase-cycled bSSFP, we propose SIMulation-based Physics-guided Learning of neural network for qMT parameters EXtraction, or SIMPLEX. In contrast to previous deep learning supervised approaches for quantitative MR that require the acquisition of input data and corresponding ground truth for training, we leveraged the MR signal model to generate training samples without expensive data curation. The network was trained exclusively with simulation data by predicting the simulation parameters. The same network was applied directly to in-vivo data without additional training. The approach was verified with both simulation and in-vivo data. SIMPLEX showed a decrease in fitting mean squared error for all simulation data compared to the existing least-square fitting method. The in-vivo experiment revealed that the network performed well with the real in vivo data unseen during training. For all experiments, we observed that SIMPLEX consistently improved the quantification quality of the qMT parameters whilst being more robust to noise compared to the prior technique. The proposed SIMPLEX will expedite the routine clinical application of qMT by providing qMT parameters (exchange rate, pool fraction) as well as T1, T2, and ΔB0 maps simultaneously with high spatial resolution, better reliability, and reduced processing time.

Motion resilience of the balanced steady-state free precession geometric solution.

PURPOSE: Many MRI sequences are sensitive to motion and its associated artifacts. The linearized geometric solution (LGS), a balanced steady-state free precession (bSSFP) off-resonance signal demodulation technique, is evaluated with respect to motion artifact resilience. THEORY AND METHODS: The mechanism and extent of LGS motion artifact resilience is examined in simulated, flow phantom, and in vivo clinical imaging. Motion artifact correction capabilities are decoupled from susceptibility artifact correction when feasible to permit controlled analysis of motion artifact correction when comparing the LGS with standard and phase-cycle-averaged (complex sum) bSSFP imaging. RESULTS: Simulations reveal that the LGS demonstrates motion artifact reduction capabilities similar to standard clinical bSSFP imaging techniques, with slightly greater resilience in high SNR regions and for shorter-duration motion. Flow phantom experiments assert that the LGS reduces shorter-duration motion artifact error by ∼24%-65% relative to the complex sum, whereas reconstructions exhibit similar error reduction for constant motion. In vivo analysis demonstrates that in the internal auditory canal/orbits, the LGS was deemed to have less artifact in 24%/49% and similar artifact in 76%/51% of radiological assessments relative to the complex sum, and the LGS had less artifact in 97%/81% and similar artifact in 3%/16% of assessments relative to standard bSSFP. Only 2 of 63 assessments deemed the LGS inferior to either complex sum or standard bSSFP in terms of artifact reduction. CONCLUSION: The LGS provides sufficient bSSFP motion artifact resilience to permit robust elimination of susceptibility artifacts, inspiring its use in a wide variety of applications.

Spectrally selective bSSFP using off-resonant RF excitations permits high spatiotemporal resolution 3D metabolic imaging of hyperpolarized [1-13 C]Pyruvate-to-[1-13 C]lactate conversion.

PURPOSE: To develop a high spatiotemporal resolution 3D dynamic pulse sequence for preclinical imaging of hyperpolarized [1-13 C]pyruvate-to-[1-13 C]lactate metabolism at 7T. METHODS: A standard 3D balanced SSFP (bSSFP) sequence was modified to enable alternating-frequency excitations. RF pulses with 2.33 ms duration and 900 Hz FWHM were placed off-resonance of the target metabolites, [1-13 C]pyruvate (by approximately -245 Hz) and [1-13 C]lactate (by approximately 735 Hz), to selectively excite those resonances. Relatively broad bandwidth (compared to those metabolites' chemical shift offset) permits a short TR of 6.29 ms, enabling higher spatiotemporal resolution. Bloch equation simulations of the bSSFP response profile guided the sequence parameter selection to minimize spectral contamination between metabolites and preserve magnetization over time. RESULTS: Bloch equation simulations, phantom studies, and in vivo studies demonstrated that the two target resonances could be cleanly imaged without substantial bSSFP banding artifacts and with little spectral contamination between lactate and pyruvate and from pyruvate hydrate. High spatiotemporal resolution 3D images were acquired of in vivo pyruvate-lactate metabolism in healthy wild-type and endogenous pancreatic tumor-bearing mice, with 1.212 s acquisition time per single-metabolite image and (1.75 mm)3 isotropic voxels with full mouse abdomen 56 × 28 × 21 mm3 FOV and fully-sampled k-space. Kidney and tumor lactate/pyruvate ratios of two consecutive measurements in one animal, 1 h apart, were consistent. CONCLUSION: Spectrally selective bSSFP using off-resonant RF excitations can provide high spatio-temporal resolution 3D dynamic images of pyruvate-lactate metabolic conversion.

SPARCQ: A new approach for fat fraction mapping using asymmetries in the phase-cycled balanced SSFP signal profile.

PURPOSE: To develop SPARCQ (Signal Profile Asymmetries for Rapid Compartment Quantification), a novel approach to quantify fat fraction (FF) using asymmetries in the phase-cycled balanced SSFP (bSSFP) profile. METHODS: SPARCQ uses phase-cycling to obtain bSSFP frequency profiles, which display asymmetries in the presence of fat and water at certain TRs. For each voxel, the measured signal profile is decomposed into a weighted sum of simulated profiles via multi-compartment dictionary matching. Each dictionary entry represents a single-compartment bSSFP profile with a specific off-resonance frequency and relaxation time ratio. Using the results of dictionary matching, the fractions of the different off-resonance components are extracted for each voxel, generating quantitative maps of water and FF and banding-artifact-free images for the entire image volume. SPARCQ was validated using simulations, experiments in a water-fat phantom and in knees of healthy volunteers. Experimental results were compared with reference proton density FFs obtained with 1 H-MRS (phantoms) and with multiecho gradient-echo MRI (phantoms and volunteers). SPARCQ repeatability was evaluated in six scan-rescan experiments. RESULTS: Simulations showed that FF quantification is accurate and robust for SNRs greater than 20. Phantom experiments demonstrated good agreement between SPARCQ and gold standard FFs. In volunteers, banding-artifact-free quantitative maps and water-fat-separated images obtained with SPARCQ and ME-GRE demonstrated the expected contrast between fatty and non-fatty tissues. The coefficient of repeatability of SPARCQ FF was 0.0512. CONCLUSION: SPARCQ demonstrates potential for fat quantification using asymmetries in bSSFP profiles and may be a promising alternative to conventional FF quantification techniques.

Deuterium imaging of the Warburg effect at sub-millimolar concentrations by joint processing of the kinetic and spectral dimensions.

Deuterium metabolic imaging (DMI) is a promising molecular MRI approach, which follows the administration of deuterated substrates and their metabolization. [6,6'-2 H2 ]-glucose for instance is preferentially converted in tumors to [3,3'-2 H2 ]-lactate as a result of the Warburg effect, providing a distinct resonance whose mapping using time-resolved spectroscopic imaging can diagnose cancer. The MR detection of low-concentration metabolites such as lactate, however, is challenging. It has been recently shown that multi-echo balanced steady-state free precession (ME-bSSFP) increases the signal-to-noise ratio (SNR) of these experiments approximately threefold over regular chemical shift imaging; the present study examines how DMI's sensitivity can be increased further by advanced processing methods. Some of these, such as compressed sensing multiplicative denoising and block-matching/3D filtering, can be applied to any spectroscopic/imaging methods. Sensitivity-enhancing approaches were also specifically tailored to ME-bSSFP DMI, by relying on priors related to the resonances' positions and to features of the metabolic kinetics. Two new methods are thus proposed that use these constraints for enhancing the sensitivity of both the spectral images and the metabolic kinetics. The ability of these methods to improve DMI is evidenced in pancreatic cancer studies carried at 15.2 T, where suitable implementations of the proposals imparted eightfold or more SNR improvement over the original ME-bSSFP data, at no informational cost. Comparisons with other propositions in the literature are briefly discussed.

Quantification of Prostate Cancer Metabolism Using 3D Multiecho bSSFP and Hyperpolarized [1-13 C] Pyruvate: Metabolism Differs Between Tumors of the Same Gleason Grade.

BACKGROUND: Three-dimensional (3D) multiecho balanced steady-state free precession (ME-bSSFP) has previously been demonstrated in preclinical hyperpolarized (HP) 13 C-MRI in vivo experiments, and it may be suitable for clinical metabolic imaging of prostate cancer (PCa). PURPOSE: To validate a signal simulation framework for the use of sequence parameter optimization. To demonstrate the feasibility of ME-bSSFP for HP 13 C-MRI in patients. To evaluate the metabolism in PCa measured by ME-bSSFP. STUDY TYPE: Retrospective single-center cohort study. PHANTOMS/POPULATION: Phantoms containing aqueous solutions of [1-13 C] lactate (2.3 M) and [13 C] urea (8 M). Eight patients (mean age 67 ± 6 years) with biopsy-confirmed Gleason 3 + 4 (n = 7) and 4 + 3 (n = 1) PCa. FIELD STRENGTH/SEQUENCES: 1 H MRI at 3 T with T2 -weighted turbo spin-echo sequence used for spatial localization and spoiled dual gradient-echo sequence used for B0 -field measurement. ME-bSSFP sequence for 13 C MR spectroscopic imaging with retrospective multipoint IDEAL metabolite separation. ASSESSMENT: The primary endpoint was the analysis of pyruvate-to-lactate conversion in PCa and healthy prostate regions of interest (ROIs) using model-free area under the curve (AUC) ratios and a one-directional kinetic model (kP ). The secondary objectives were to investigate the correlation between simulated and experimental ME-bSSFP metabolite signals for HP 13 C-MRI parameter optimization. STATISTICAL TESTS: Pearson correlation coefficients with 95% confidence intervals and paired t-tests. The level of statistical significance was set at P < 0.05. RESULTS: Strong correlations between simulated and empirical ME-bSSFP signals were found (r > 0.96). Therefore, the simulation framework was used for sequence optimization. Whole prostate metabolic HP 13 C-MRI, observing the conversion of pyruvate into lactate, with a temporal resolution of 6 seconds was demonstrated using ME-bSSFP. Both assessed metrics resulted in significant differences between PCa (mean ± SD) (AUC = 0.33 ± 012, kP  = 0.038 ± 0.014) and healthy (AUC = 0.15 ± 0.10, kP  = 0.011 ± 0.007) ROIs. DATA CONCLUSION: Metabolic HP 13 C-MRI in the prostate using ME-bSSFP allows for differentiation between aggressive PCa and healthy tissue. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 1.

Dynamic cardiac MRI with high spatiotemporal resolution using accelerated spiral-out and spiral-in/out bSSFP pulse sequences at 1.5 T.

OBJECTIVE: To develop two spiral-based bSSFP pulse sequences combined with L + S reconstruction for accelerated ungated, free-breathing dynamic cardiac imaging at 1.5 T. MATERIALS AND METHODS: Tiny golden angle rotated spiral-out and spiral-in/out bSSFP sequences combined with view-sharing (VS), compressed sensing (CS), and low-rank plus sparse (L + S) reconstruction were evaluated and compared via simulation and in vivo dynamic cardiac imaging studies. The proposed methods were then validated against the standard cine, in terms of quantitative image assessment and qualitative quality rating. RESULTS: The L + S method yielded the least residual artifacts and the best image sharpness among the three methods. Both spiral cine techniques showed clinically diagnostic images (score > 3). Compared to standard cine, there were significant differences in global image quality and edge sharpness for spiral cine techniques, while there was significant difference in image contrast for the spiral-out cine but no significant difference for the spiral-in/out cine. There was good agreement in left ventricular ejection fraction for both the spiral-out cine (- 1.6 [Formula: see text] 3.1%) and spiral-in/out cine (- 1.5 [Formula: see text] 2.8%) against standard cine. DISCUSSION: Compared to the time-consuming standard cine (~ 5 min) which requires ECG-gating and breath-holds, the proposed spiral bSSFP sequences achieved ungated, free-breathing cardiac movies at a similar spatial (1.5 × 1.5 × 8 mm3) and temporal resolution (36 ms) per slice for whole heart coverage (10-15 slices) within 45 s, suggesting the clinical potential for improved patient comfort or for imaging patients with arrhythmias or who cannot hold their breath.

Pure balanced steady-state free precession imaging (pure bSSFP).

PURPOSE: To show that for tissues the conspicuous asymmetries in the frequency response function of bSSFP can be mitigated by using a short enough TR. THEORY AND METHODS: Configuration theory indicates that bSSFP becomes apparently "pure" (i.e., exhibiting a symmetric profile) in the limit of TR → 0 . To this end, the frequency profile of bSSFP was measured as a function of the TR using a manganese-doped aqueous probe, as well as brain tissue that was shown to exhibit a pronounced asymmetry due to its microstructure. The frequency response function was sampled using N = 72 (phantom) and N = 36 (in vivo) equally distributed linear RF phase increments in the interval [ 0 , 2 π ) . Imaging was performed with 2.0 mm isotropic resolution over a TR range of 1.5-8 ms at 3 and 1.5 T. RESULTS: As expected, pure substances showed a symmetric TR-independent frequency profile, whereas brain tissue revealed a pronounced asymmetry. The observed asymmetry for the tissue, however, decreases with decreasing TR and gives strong evidence that the frequency response function of bSSFP becomes symmetric in the limit of TR → 0 , in agreement with theory. The limit of apparently pure bSSFP imaging can thus be achieved for a TR ∼ 1.5 ms at 1.5 T, whereas at 3 T, tissues still show some residual asymmetry. CONCLUSION: In the limit of short enough TR, tissues become apparently pure for bSSFP. This limit can be reached for brain tissue at 1.5 T with TR ∼ 1-2 ms at clinically relevant resolutions.

Development of specialized magnetic resonance acquisition techniques for human hyperpolarized [13 C,15 N2 ]urea + [1-13 C]pyruvate simultaneous perfusion and metabolic imaging.

PURPOSE: This study aimed to develop and demonstrate the in vivo feasibility of a 3D stack-of-spiral balanced steady-state free precession(3D-bSSFP) urea sequence, interleaved with a metabolite-specific gradient echo (GRE) sequence for pyruvate and metabolic products, for improving the SNR and spatial resolution of the first hyperpolarized 13 C-MRI human study with injection of co-hyperpolarized [1-13 C]pyruvate and [13 C,15 N2 ]urea. METHODS: A metabolite-specific bSSFP urea imaging sequence was designed using a urea-specific excitation pulse, optimized TR, and 3D stack-of-spiral readouts. Simulations and phantom studies were performed to validate the spectral response of the sequence. The image quality of urea data acquired by the 3D-bSSFP sequence and the 2D-GRE sequence was evaluated with 2 identical injections of co-hyperpolarized [1-13 C]pyruvate and [13 C,15 N2 ]urea formula in a rat. Subsequently, the feasibility of the acquisition strategy was validated in a prostate cancer patient. RESULTS: Simulations and phantom studies demonstrated that 3D-bSSFP sequence achieved urea-only excitation, while minimally perturbing other metabolites (<1%). An animal study demonstrated that compared to GRE, bSSFP sequence provided an ∼2.5-fold improvement in SNR without perturbing urea or pyruvate kinetics, and bSSFP approach with a shorter spiral readout reduced blurring artifacts caused by J-coupling of [13 C,15 N2 ]urea. The human study demonstrated the in vivo feasibility and data quality of the acquisition strategy. CONCLUSION: The 3D-bSSFP urea sequence with a stack-of-spiral acquisition demonstrated significantly increased SNR and image quality for [13 C,15 N2 ]urea in co-hyperpolarized [1-13 C]pyruvate and [13 C,15 N2 ]urea imaging studies. This work lays the foundation for future human studies to achieve high-quality and high-SNR metabolism and perfusion images.

High-resolution neural network-driven mapping of multiple diffusion metrics leveraging asymmetries in the balanced steady-state free precession frequency profile.

We propose to utilize the rich information content about microstructural tissue properties entangled in asymmetric balanced steady-state free precession (bSSFP) profiles to estimate multiple diffusion metrics simultaneously by neural network (NN) parameter quantification. A 12-point bSSFP phase-cycling scheme with high-resolution whole-brain coverage is employed at 3 and 9.4 T for NN input. Low-resolution target diffusion data are derived based on diffusion-weighted spin-echo echo-planar-imaging (SE-EPI) scans, that is, mean, axial, and radial diffusivity (MD, AD, and RD), fractional anisotropy (FA), as well as the spherical coordinates (azimuth Φ and inclination Ï´) of the principal diffusion eigenvector. A feedforward NN is trained with incorporated probabilistic uncertainty estimation. The NN predictions yielded highly reliable results in white matter (WM) and gray matter structures for MD. The quantification of FA, AD, and RD was overall in good agreement with the reference but the dependence of these parameters on WM anisotropy was somewhat biased (e.g. in corpus callosum). The inclination Ï´ was well predicted for anisotropic WM structures, while the azimuth Φ was overall poorly predicted. The findings were highly consistent across both field strengths. Application of the optimized NN to high-resolution input data provided whole-brain maps with rich structural details. In conclusion, the proposed NN-driven approach showed potential to provide distortion-free high-resolution whole-brain maps of multiple diffusion metrics at high to ultrahigh field strengths in clinically relevant scan times.