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.

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.

Fast and accessible T2 mapping using off-resonance corrected DESPOT2 with application to 3D prostate.

PURPOSE: Most T1 and T2 mapping take long acquisitions or needs specialized sequences not widely accessible on clinical scanners. An available solution is DESPOT1/T2 (Driven equilibrium single pulse observation of T1/T2). DESPOT1/T2 uses Spoiled gradient-echo (SPGR) and balanced Steady-State Free Precession (bSSFP) sequences, offering an accessible and reliable way for 3D accelerated T1/T2 mapping. However, bSSFP is prone to off-resonance artifacts, limiting the application of DESPOT2 in regions with high susceptibility contrasts, like the prostate. Our proposal, DESPO+, employs the full bSSFP and SPGR models with a dictionary-based method to reconstruct 3D T1/T2 maps in the prostate region without off-resonance banding. METHODS: DESPO+ modifies the bSSFP acquisition of the original variable flip angle DESPOT2. DESPO+ uses variable repetition and echo times, employing a dictionary-based method of the full bSSFP and SPGR models to reconstruct T1, T2, and Proton Density (PD) simultaneously. The proposed DESPO+ method underwent testing through simulations, T1/T2 phantoms, and on fourteen healthy subjects. RESULTS: The results reveal a significant reduction in T2 map banding artifacts compared to the original DESPOT2 method. DESPO+ approach reduced T2 errors by up to seven times compared to DESPOT2 in simulations and phantom experiments. We also synthesized in-vivo T1-weighted/T2-weighted images from the acquired maps using a spin-echo model to verify the map's quality when lacking a reference. For in-vivo imaging, the synthesized images closely resemble those from the clinical MRI protocol, reducing scan time by around 50% compared to traditional spin-echo T1-weighted/T2-weighted acquisitions. CONCLUSION: DESPO+ provides an off-resonance insensitive and clinically available solution, enabling high-resolution 3D T1/T2 mapping and synthesized T1-weighted/T2-weighted images for the entire prostate, all achieved within a short scan time of 3.6 min, similar to DESPOT1/T2.

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.

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.

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.

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.

Fast T2 mapping of short-T2 tissues in knee using 3D radial dual-echo balanced steady-state free precession.

BACKGROUND: T2 mapping of short-T2 tissues in the knee (meniscus, tendon, and ligament) is needed to aid the clinical MRI knee diagnosis, which is hard to realize using traditional clinical methods. PURPOSE: To accelerate the acquisition of T2 values for short-T2 tissues in the knee by analyzing the signal equation of balanced steady-state free precession (bSSFP) sequence in MRI. METHODS: Effect of half-radial acquisition on pixel bandwidth was analyzed mathematically. A modified 3D radial dual-echo bSSFP sequence was proposed for 0.53 mm isotropic resolution knee imaging with 2 different TEs at 3 T, which alleviated the problem of off-resonance artifacts caused by traditional half-radial acquisition scheme. A novel pixel-based optimization method was proposed for efficient T2 mapping of short-T2 tissues in the knee given off-resonance values. Simulation was conducted to evaluate the sensitivity of the proposed method to other parameters. Phantom results were compared with 2D spin-echo (SE), and in vivo results were compared with SE and previously studies. RESULTS: Simulation showed that the proposed method is insensitive to T1 and B1 variations (estimation error < 1% for T1/B1 error of ±90%), avoiding the need for separated T1 and B1 scans. High isotropic resolution knee imaging was achieved using the modified dual-echo bSSFP. The total scan time was within 3.5 min, including a separate off-resonance scan for T2 measurement. Measured mean T2 values for phantoms correlated well with SE (R2 = 0.99), and no significant difference was observed (P = 0.45). In vivo meniscus T2 measurements and ligament T2 measurements agreed with the literature, while tendon T2 measurements were much lower (31.7% lower for patellar tendon, and 13.5% lower for quadriceps tendon), which might result in its bi-component property. CONCLUSIONS: The proposed method provides an efficient way for fast, robust, high-resolution imaging and T2 mapping of short-T2 tissues in the knee.

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.

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.

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.

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.

A metabolite specific 3D stack-of-spirals bSSFP sequence for improved bicarbonate imaging in hyperpolarized [1-13C]Pyruvate MRI.

13C-bicarbonate is a crucial measure of pyruvate oxidation and TCA cycle flux, but is challenging to measure due to its relatively low concentration and thus will greatly benefit from improved signal-to-noise ratio (SNR). To address this, we developed and investigated the feasibility of a 3D stack-of-spirals metabolite-specific balanced steady-state free precession (MS-bSSFP) sequence for improving the SNR and spatial resolution of dynamic 13C-bicarbonate imaging in hyperpolarized [1-13C]pyruvate studies. The bicarbonate MS-bSSFP sequence was evaluated by simulations, phantoms studies, preclinical studies on five rats, brain studies on two healthy volunteers and renal study on one renal cell carcinoma patient. The simulations and phantom results showed that the bicarbonate-specific pulse had minimal perturbation of other metabolites (<1%). In the animal studies, the MS-bSSFP sequence provided an approximately 2.6-3 × improvement in 13C-bicarbonate SNR compared to a metabolite-specific gradient echo (MS-GRE) sequence without altering the bicarbonate or pyruvate kinetics, and the shorter spiral readout in the MS-bSSFP approach reduced blurring. Using the SNR ratio between MS-bSSFP and MS-GRE, the T2 values of bicarbonate and lactate in the rat kidneys were estimated as 0.5 s and 1.1 s, respectively. The in-vivo feasibility of bicarbonate MS-bSSFP sequence was demonstrated in two human brain studies and one renal study. These studies demonstrate the potential of the sequence for in-vivo applications, laying the foundation for future studies to observe this relatively low concentration metabolite with high-quality images and improve measurements of pyruvate oxidation.

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.

High resolution single-shot myocardial imaging using bSSFP with wave encoding.

BACKGROUND: Single-shot balanced steady-state free precession (bSSFP) sequence is widely used in cardiac imaging. However, the limited scan time in one heartbeat greatly hinders its spatial resolution compared to the segmented acquisition mode. Therefore, a highly accelerated single-shot bSSFP imaging technology is needed for clinical use. PURPOSE: To develop and evaluate a wave-encoded bSSFP sequence with high acceleration rates for single-shot myocardial imaging. METHODS: The proposed Wave-bSSFP method is implemented by adding a sinusoidal wave gradient in the phase encoding direction during the readout of bSSFP sequence. Uniform undersampling is used for acceleration. Its performance was first validated via phantom studies by comparison with conventional bSSFP. Then it was evaluated in volunteer studies via anatomical imaging, T2 -prepared bSSFP, and T1 mapping in in-vivo cardiac imaging. All methods were compared with accelerated conventional bSSFP reconstructed using iterative SENSE and compressed sensing (CS), to demonstrate the advantage of wave encoding in suppressing the noise amplification and artifacts induced by acceleration. RESULTS: The proposed Wave-bSSFP method achieved a high acceleration factor of 4 for single-shot acquisitions. The proposed method showed lower average g-factors than bSSFP, and fewer blurring artifacts than CS reconstruction. The Wave-bSSFP with R = 4 achieved higher spatial and temporal resolutions compared with the conventional bSSFP with R = 2 in several applications such as T2 -prepared bSSFP and T1 mapping, and could be applied in systolic imaging. CONCLUSION: Wave encoding can be used to highly accelerate 2D bSSFP imaging with single-shot acquisitions. Compared with the conventional bSSFP sequence, the proposed Wave-bSSFP method can effectively reduce the g-factor and aliasing artifacts in cardiac imaging.

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.

Non-contrast MR angiography of pelvic arterial vasculature using the Quiescent interval slice selective (QISS) sequence.

To evaluate Quiescent Interval Slice Selective (QISS) balanced steady-state free precession (bSSFP) and QISS fast low-angle shot (FLASH) sequences for non-contrast Magnetic Resonance Angiography (MRA) of iliac arteries regarding image quality and diagnostic confidence in order to establish these sequences in daily clinical practice. A prospective study of healthy subjects (n = 10) was performed. All subjects underwent the QISS MRI protocol with bSSFP und FLASH sequences. Vessel contrast-to-background ratio (VCBR) were measured in pre-defined vessel segments. Image quality and diagnostic confidence was assessed using a Likert scale (five-point scale). Inter-reader agreement was determined using Cohen's kappa coefficient (κ). Ten healthy subjects (median age 29 years, IQR: 26.25 to 30 years) were included in this prospective study. Median MR examination time was 2:05 min (IQR 1:58 to 2:16) for QISS bSSFP and 4:11 min (IQR 3:57 to 4:32) for QISS FLASH. Both sequences revealed good VCBR in all examined vessel segments. VCBR (muscle tissue) were marginally higher for FLASH sequences (e.g., 0.82 vs. 0.78 in the right femoral artery, p = 0.035*), while bSSFP sequence showed significantly higher VCBR (fat tissue) in the majority of examined arterials vessels (e.g., 0.78 vs. 0.62 in right femoral artery, p = 0.001*). The image quality and diagnostic confidence of both sequences were rated as good to excellent. Moderate to good inter-reader agreement was found. QISS MRA using bSSFP and FLASH sequences are diagnostic for visualization of iliac arterial vasculature. The QISS bSSFP sequence might offer advantages due to the markedly shorter exam time and superior visualization of smaller vessels. The QISS FLASH sequence seems to be a robust alternative for non-contrast MRA since it is less sensitive to magnetic field inhomogeneities.

Simultaneous multi-slice T1 mapping using MOLLI with blipped CAIPIRINHA bSSFP.

BACKGROUND: This study evaluates the possibility for replacing conventional 3 slices, 3 breath-holds MOLLI cardiac T1 mapping with single breath-hold 3 simultaneous multi-slice (SMS3) T1 mapping using blipped-CAIPIRINHA SMS-bSSFP MOLLI sequence. As a major drawback, SMS-bSSFP presents unique artefacts arising from side-lobe slice excitations that are explained by imperfect RF modulation rendering and bSSFP low flip angle enhancement. Amplitude-only RF modulation (AM) is proposed to reduce these artefacts in SMS-MOLLI compared to conventional Wong multi-band RF modulation (WM). MATERIALS AND METHODS: Phantoms and ten healthy volunteers were imaged at 1.5 T using a modified blipped-CAIPIRINHA SMS-bSSFP MOLLI sequence with 3 simultaneous slices. WM-SMS3 and AM-SMS3 were compared to conventional single-slice (SMS1) MOLLI. First, SNR degradation and T1 accuracy were measured in phantoms. Second, artefacts from side-lobe excitations were evaluated in a phantom designed to reproduce fat presence near the heart. Third, the occurrence of these artefacts was observed in volunteers, and their impact on T1 quantification was compared between WM-SMS3 and AM-SMS3 with conventional MOLLI as a reference. RESULTS: In the phantom, larger slice gaps and slice thicknesses yielded higher SNR. There was no significant difference of T1 values between conventional MOLLI and SMS3-MOLLI (both WM and AM). Positive banding artefacts were identified from fat neighbouring the targeted FOV due to side-lobe excitations from WM and the unique bSSFP signal profile. AM RF pulses reduced these artefacts by 38%. In healthy volunteers, AM-SMS3-MOLLI showed similar artefact reduction compared to WM-SMS3-MOLLI (3 ± 2 vs 5 ± 3 corrupted LV segments out of 16). In-vivo native T1 values obtained from conventional MOLLI and AM-SMS3-MOLLI were equivalent in LV myocardium (SMS1-T1 = 935.5 ± 36.1 ms; AM-SMS3-T1 = 933.8 ± 50.2 ms; P = 0.436) and LV blood pool (SMS1-T1 = 1475.4 ± 35.9 ms; AM-SMS3-T1 = 1452.5 ± 70.3 ms; P = 0.515). Identically, no differences were found between SMS1 and SMS3 postcontrast T1 values in the myocardium (SMS1-T1 = 556.0 ± 19.7 ms; SMS3-T1 = 521.3 ± 28.1 ms; P = 0.626) and the blood (SMS1-T1 = 478 ± 65.1 ms; AM-SMS3-T1 = 447.8 ± 81.5; P = 0.085). CONCLUSIONS: Compared to WM RF modulation, AM SMS-bSSFP MOLLI was able to reduce side-lobe artefacts considerably, providing promising results to image the three levels of the heart in a single breath hold. However, few artefacts remained even using AM-SMS-bSSFP due to residual RF imperfections. The proposed blipped-CAIPIRINHA MOLLI T1 mapping sequence provides accurate in vivo T1 quantification in line with those obtained with a single slice acquisition.