Quantitative T2 mapping using accelerated 3D stack-of-spiral gradient echo readout

Whole-cerebrum guanidino and amide CEST mapping at 3 T by a 3D stack-of-spirals gradient echo acquisition

PURPOSE

To develop a 3D, high-sensitivity CEST mapping technique based on the 3D stack-of-spirals (SOS) gradient echo readout, the proposed approach was compared with conventional acquisition techniques and evaluated for its efficacy in concurrently mapping of guanidino (Guan) and amide CEST in human brain at 3 T, leveraging the polynomial Lorentzian line-shape fitting (PLOF) method.

METHODS

Saturation time and recovery delay were optimized to achieve maximum CEST time efficiency. The 3DSOS method was compared with segmented 3D EPI (3DEPI), turbo spin echo, and gradient- and spin-echo techniques. Image quality, temporal SNR (tSNR), and test-retest reliability were assessed. Maps of Guan and amide CEST derived from 3DSOS were demonstrated on a low-grade glioma patient.

RESULTS

The optimized recovery delay/saturation time was determined to be 1.4/2 s for Guan and amide CEST. In addition to nearly doubling the slice number, the gradient echo techniques also outperformed spin echo sequences in tSNR: 3DEPI (193.8 ± 6.6), 3DSOS (173.9 ± 5.6), and GRASE (141.0 ± 2.7). 3DSOS, compared with 3DEPI, demonstrated comparable GuanCEST signal in gray matter (GM) (3DSOS: [2.14%-2.59%] vs. 3DEPI: [2.15%-2.61%]), and white matter (WM) (3DSOS: [1.49%-2.11%] vs. 3DEPI: [1.64%-2.09%]). 3DSOS also achieves significantly higher amideCEST in both GM (3DSOS: [2.29%-3.00%] vs. 3DEPI: [2.06%-2.92%]) and WM (3DSOS: [2.23%-2.66%] vs. 3DEPI: [1.95%-2.57%]). 3DSOS outperforms 3DEPI in terms of scan-rescan reliability (correlation coefficient: 3DSOS: 0.58-0.96 vs. 3DEPI: -0.02 to 0.75) and robustness to motion as well.

CONCLUSION

The 3DSOS CEST technique shows promise for whole-cerebrum CEST imaging, offering uniform contrast and robustness against motion artifacts.

Accelerated 3D multi-echo spin-echo sequence with a subspace constrained reconstruction for whole mouse brain T 2 mapping

PURPOSE

To accelerate whole-brain quantitativeT 2 $$ {\mathrm{T}}_2 $$ mapping in preclinical imaging setting.

METHODS

A three-dimensional (3D) multi-echo spin echo sequence was highly undersampled with a variable density Poisson distribution to reduce the acquisition time. Advanced iterative reconstruction based on linear subspace constraints was employed to recover high-quality raw images. Different subspaces, generated using exponential or extended-phase graph (EPG) simulations or from low-resolution calibration images, were compared. The subspace dimension was investigated in terms ofT 2 $$ {\mathrm{T}}_2 $$ precision. The method was validated on a phantom containing a wide range ofT 2 $$ {\mathrm{T}}_2 $$ and was then applied to monitor metastasis growth in the mouse brain at 4.7T. Image quality andT 2 $$ {\mathrm{T}}_2 $$ estimation were assessed for 3 acceleration factors (6/8/10).

RESULTS

The EPG-based dictionary gave robust estimations of a large range ofT 2 $$ {\mathrm{T}}_2 $$ . A subspace dimension of 6 was the best compromise betweenT 2 $$ {\mathrm{T}}_2 $$ precision and image quality. Combining the subspace constrained reconstruction with a highly undersampled dataset enabled the acquisition of whole-brainT 2 $$ {\mathrm{T}}_2 $$ maps, the detection and the monitoring of metastasis growth of less than 500μ m 3 $$ \mu {\mathrm{m}}^3 $$ .

CONCLUSION

Subspace-based reconstruction is suitable for 3DT 2 $$ {\mathrm{T}}_2 $$ mapping. This method can be used to reach an acceleration factor up to 8, corresponding to an acquisition time of 25 min for an isotropic 3D acquisition of 156μ $$ \mu $$ m on the mouse brain, used here for monitoring metastases growth.

Integrated Short-TE and Hadamard-edited Multi-Sequence (ISTHMUS) for advanced MRS

BACKGROUND

To examine data quality and reproducibility using ISTHMUS, which has been implemented as the standardized MR spectroscopy sequence for the multi-site Healthy Brain and Child Development (HBCD) study.

METHODS

ISTHMUS is the consecutive acquisition of short-TE PRESS (32 transients) and long-TE HERCULES (224 transients) data with dual-TE water reference scans. Voxels were positioned in the centrum semiovale, dorsal anterior cingulate cortex, posterior cingulate cortex and bilateral thalamus regions. After acquisition, ISTHMUS data were separated into the PRESS and HERCULES portions for analysis and modeled separately using Osprey. In vivo experiments were performed in 10 healthy volunteers (6 female; 29.5±6.6 years). Each volunteer underwent two scans on the same day. Differences in metabolite measurements were examined. T2 correction based on the dual-TE water integrals were compared with: 1) T2 correction based on the default white matter and gray matter T2 reference values in Osprey and 2) shorter WM and GM T2 values from recent literature.

RESULTS

No significant difference in linewidth was observed between PRESS and HERCULES. Bilateral thalamus spectra had produced significantly higher (p<0.001) linewidth compared to the other three regions. Linewidth measurements were similar between scans, with scan-to-scan differences under 1 Hz for most subjects. Paired t-tests indicated a significant difference only in PRESS NAAG between the two thalamus scans (p=0.002). T2 correction based on shorter T2 values showed better agreement to the dual-TE water integral ratio.

CONCLUSIONS

ISTHMUS facilitated data acquisition and post-processing and reduced operator workload to eliminate potential human error.

Integrated Short-TE and Hadamard-edited Multi-Sequence (ISTHMUS) for Advanced MRS

Background

To examine data quality and reproducibility using ISTHMUS, which has been implemented as the standardized MR spectroscopy sequence for the multi-site Healthy Brain and Child Development (HBCD) study.

Methods

ISTHMUS is the consecutive acquisition of short-TE PRESS (32 transients) and long-TE HERCULES (224 transients) data with dual-TE water reference scans. Voxels were positioned in the centrum semiovale, dorsal anterior cingulate cortex, posterior cingulate cortex and bilateral thalamus regions. After acquisition, ISTHMUS data were separated into the PRESS and HERCULES portions for analysis and modeled separately using Osprey. In vivo experiments were performed in 10 healthy volunteers (6 female; 29.5±6.6 years). Each volunteer underwent two scans on the same day. Differences in metabolite measurements were examined. T2 correction based on the dual-TE water integrals were compared with: 1) T2 correction based the default white matter and gray matter T2 reference values in Osprey; 2) shorter WM and GM T2 values from recent literature; and 3) reduced CSF fractions.

Results

No significant difference in linewidth was observed between PRESS and HERCULES. Bilateral thalamus spectra had produced significantly higher (p<0.001) linewidth compared to the other three regions. Linewidth measurements were similar between scans, with scan-to-scan differences under 1 Hz for most subjects. Paired t-tests indicated a significant difference only in PRESS NAAG between the two thalamus scans (p=0.002). T2 correction based on shorter T2 values showed better agreement to the dual-TE water integral ratio.

Conclusions

ISTHMUS facilitated and standardized acquisition and post-processing and reduced operator workload to eliminate potential human error.

Technical note: Accuracy and precision of T2 and T2* with a gradient-echo spin-echo (GESE) sequence for cardiac imaging

BACKGROUND

The use of a gradient echo spin echo (GESE) method to obtain rapid T2 and T2* estimation in the heart has been proposed. The effect of acquisition parameter settings on T2 and T2* bias and precision have not been investigated in depth.

PURPOSE

To understand factors impacting the quantification of T2 and T2* values with a gradient echo spin echo (GESE) method using echo planar imaging (EPI) readouts in a reduced field of view acquisition.

METHODS

The GESE method is implemented with a reduced field-of-view using an outer volume suppression (OVS) technique to minimize the time for multi-echo EPI readouts. The number of EPI readouts (images) for the GESE is optimized using Cramer-Rao Lower Bound (CRLB) and Monte Carlo simulations with a nonlinear least-square (NLLS) estimator. The SNR requirements were studied using the latter simulation method for a selected range of T2 and T2* values and T2/T2* ratios. Two healthy control subjects were imaged with the proposed GESE sequence and evaluated with the NLLS estimation method. In addition, the proposed OVS method was compared with a saturation bands OVS method in one subject. Clinical T2 and T2* mappings were used as the reference.

RESULTS

The optimal number of EPI readouts is five and the performance is slightly better when the refocusing pulse is placed between the 2nd and 3rd readouts. The SNR requirement for achieving a target bias < 1 ms and standard deviation (SD) < 5 ms is more demanding when T2/T2* ratio increases. The minimum SNR requirement in the GESE acquisition should vary from 6 to 20 depending on specific myocardial T2 and T2* values at 3T. The T2 and T2* estimates using the proposed OVS method and the saturation bands OVS method are both similar to the reference.

CONCLUSION

The GESE sequence with five EPI readouts is a feasible and efficient technique that can estimate T2 and T2* values in the septal myocardium within a heartbeat when the SNR requirement can be satisfied.

Evaluation of 3D stack-of-spiral turbo FLASH acquisitions for pseudo-continuous and velocity-selective ASL-derived brain perfusion mapping

PURPOSE

The most-used 3D acquisitions for ASL are gradient and spin echo (GRASE)- and stack-of-spiral (SOS)-based fast spin echo, which require multiple shots. Alternatively, turbo FLASH (TFL) allows longer echo trains, and SOS-TFL has the potential to reduce the number of shots to even single-shot, thus improving the temporal resolution. Here we compare the performance of 3D SOS-TFL and 3D GRASE for ASL at 3T.

METHODS

The 3D SOS-TFL readout was optimized with respect to fat suppression and excitation flip angles for pseudo-continuous ASL- and velocity-selective (VS)ASL-derived cerebral blood flow (CBF) mapping as well as for VSASL-derived cerebral blood volume (CBV) mapping. Results were compared with 3D GRASE readout on healthy volunteers in terms of perfusion quantification and temporal SNR (tSNR) efficiency. CBF and CBV mapping derived from 3D SOS-TFL-based ASL was demonstrated on one stroke patient, and the potential for single-shot acquisitions was exemplified.

RESULTS

SOS-TFL with a 15° flip angle resulted in adequate tSNR efficiency with negligible image blurring. Selective water excitation was necessary to eliminate fat-induced artifacts. For pseudo-continuous ASL- and VSASL-based CBF and CBV mapping, compared to the employed four-shot 3D GRASE with an acceleration factor of 2, the fully sampled 3D SOS-TFL delivered comparable performance (with a similar scan time) using three shots, which could be further undersampled to achieve single-shot acquisition with higher tSNR efficiency. SOS-TFL had reduced CSF contamination for VSASL-CBF.

CONCLUSION

3D SOS-TFL acquisition was found to be a viable substitute for 3D GRASE for ASL with sufficient tSNR efficiency, minimal relaxation-induced blurring, reduced CSF contamination, and the potential of single-shot, especially for VSASL.

Efficient phase-cycling strategy for high-resolution 3D gradient-echo quantitative parameter mapping

Magnetization-prepared (MP) gradient-echo (GRE) sequences suffer from signal contaminations from T₁ recovery during the readout train, which can be eliminated by paired RF phase cycling (PC) at the cost of doubling the scan time. The objective of this study was to develop and validate a novel unpaired PC strategy to eliminate the time penalty for high-resolution quantitative parameter mapping in 3D MP-GRE sequences. Based on the observation that the contaminating T₁ recovery signal along the GRE readout train is independent of magnetization preparation, its impact can be eliminated using a novel curve-fitting approach with complex-valued data without needing paired PC acquisitions. Four new unpaired PC schemes were compared with two traditional paired PC schemes in both phantom and in vivo human knee studies at 3 T using a MP angle-modulated partitioned k-space spoiled gradient-echo snapshots (MAPSS) T_1ρ mapping sequence. In the phantom study, all methods resulted in consistent T_1ρ measurements (∆T_1ρ  < 0.5%) at the center slice when B₀ /B₁ values were uniform. Results were not consistent when off-center slices with nonideal B₀ /B₁ were included. Two unpaired PC schemes had comparable or significantly improved quantitative accuracy and scan-rescan reproducibility compared with the paired PC schemes. There was no significant T_1ρ quantitative variability increase or spatial fidelity loss using the new unpaired PC schemes. Unpaired PC schemes also had different T_1ρ spectral responses at different B₀ frequency offsets, which can potentially be exploited to reduce sensitivity to B₀ field inhomogeneities. The human knee study results were consistent with the phantom study findings. In conclusion, an unpaired PC strategy potentially allows more accurate quantitative parameter mapping with halved scan time compared with the paired PC approach to eliminate signal contaminations from T₁ recovery. It therefore offers additional flexibility in SNR optimization, spatial resolution improvement, and choice of imaging sampling points to obtain more accurate quantitative parameter mapping.

Ultrafast B1 mapping with RF-prepared 3D FLASH acquisition: Correcting the bias due to T1 -induced k-space filtering effect

Purpose

The traditional radiofrequency (RF)‐prepared B₁ mapping technique consists of one scan with an RF preparation module for flip angle‐encoding and a second scan without this module for normalizing. To reduce the T₁‐induced k‐space filtering effect, this method is limited to 2D FLASH acquisition with a two‐parameter method. A novel 3D RF‐prepared three‐parameter method for ultrafast B₁‐mapping is proposed to correct the T₁‐induced quantification bias.

Theory

The point spread function analysis of FLASH shows that the prepared longitudinal magnetization before the FLASH acquisition and the image signal obeys a linear (not proportional) relationship. The intercept of the linear function causes the quantification bias and can be captured by a third saturated scan.

Methods

Using the 2D double‐angle method (DAM) as the reference, a 3D RF‐prepared three‐parameter protocol with 9 s duration was compared with the two‐parameter method, as well as the saturated DAM (SDAM) method, the dual refocusing echo acquisition mode (DREAM) method, and the actual flip‐angle imaging (AFI) method, for B₁ mapping of brain, breast, and abdomen with different orientations and shim settings at 3T.

Results

The 3D RF‐prepared three‐parameter method with complex‐subtraction delivered consistently lower RMS error, error mean, error standard deviation, and higher concordance correlation coefficients values than the two‐parameter method, the three‐parameter method with magnitude‐subtraction, the multi‐slice DREAM and the 3D AFI, and were close to the results of 2D or multi‐slice SDAM.

Conclusion

The proposed ultrafast 3D RF‐prepared three‐parameter method with complex‐subtraction was demonstrated with high accuracy for B₁ mapping of brain, breast, and abdomen.

A revisit of the k-space filtering effects of magnetization-prepared 3D FLASH and balanced SSFP acquisitions: Analytical characterization of the point spread functions

PURPOSE

3D FLASH and balanced SSFP (bSSFP) are increasingly used in quantitative MRI after contrast preparation. The acquired k-space data are modulated by T₁ relaxation (or additional T₂ for bSSFP). Three separate sequence parameters including the number of phase-encoding steps per shot (N), flip angle (FA), and TR have made the transient state of rapid gradient echo (GRE) imaging difficult for analysis and optimization. Here we aim to analytically characterize the k-space filtering effect of magnetization-prepared FLASH and bSSFP with the point spread functions (PSF).

METHODS

The amplitude effect is characterized with the peak magnitude of the PSF, i.e. PSF(0), which, due to their approaching from transient state to steady-state for the GRE acquisitions, obeys a linear (with a slope and an intercept, not proportional) relationship with the prepared longitudinal magnetization (M_prep). The blurring effect is characterized by the FWHM of the PSF. The magnetization-prepared acquisition-dependent image contrast efficiency is characterized with the relative contrast-to-noise ratio (CNR) per unit time (ruCNR).

RESULTS

The slope of PSF(0) characterizes the relative contrast between different M_prep levels. The intercept of PSF(0) could lead to quantification bias for magnetization-prepared imaging. FLASH and bSSFP experience very little blurring effect, which is to the contrary of conventional fast spin echo (FSE). Analytical selections of N, FA, and TR are provided to optimize ruCNR for different scenarios.

CONCLUSIONS

PSFs of the FLASH and bSSFP acquisitions are analytically derived and numerically validated, and compared with the FSE acquisition, thus providing a useful tool for optimizing magnetization-prepared GRE acquisitions.

T2 quantification in brain using 3D fast spin-echo imaging with long echo trains

Purpose

Three‐dimensional fast spin‐echo (FSE) sequences commonly use very long echo trains (>64 echoes) and severely reduced refocusing angles. They are increasingly used in brain exams due to high, isotropic resolution and reasonable scan time when using long trains and short interecho spacing. In this study, T₂ quantification in 3D FSE is investigated to achieve increased resolution when comparing with established 2D (proton‐density dual‐echo and multi‐echo spin‐echo) methods.

Methods

The FSE sequence design was explored to use long echo trains while minimizing T₂ fitting error and maintaining typical proton density and T₂‐weighted contrasts. Constant and variable flip angle trains were investigated using extended phase graph and Bloch equation simulations. Optimized parameters were analyzed in phantom experiments and validated in vivo in comparison to 2D methods for eight regions of interest in brain, including deep gray‐matter structures and white‐matter tracts.

Results

Phantom and healthy in vivo brain T₂ measurements showed that optimized variable echo‐train 3D FSE performs similarly to previous 2D methods, while achieving three‐fold‐higher slice resolution, evident visually in the 3D T₂ maps. Optimization resulted in better T₂ fitting and compared well with standard multi‐echo spin echo (within the 8‐ms confidence limits defined based on Bland‐Altman analysis).

Conclusion

T₂ mapping using 3D FSE with long echo trains and variable refocusing angles provides T₂ accuracy in agreement with 2D methods with additional high‐resolution benefits, allowing isotropic views while avoiding incidental magnetization transfer effects. Consequently, optimized 3D sequences should be considered when choosing T₂ mapping methods for high anatomic detail.

Spiral 3D time-of-flight MR angiography for rapid non-contrast carotid artery imaging: Clinical feasibility and protocol optimization

PURPOSE

To assess the clinical feasibility of spiral 3D Time-Of-Flight (TOF) MR Angiography (MRA) sequence variants for rapid non-contrast carotid artery imaging.

METHODS

Nine different 3D TOF MRA sequences were acquired in nine healthy volunteers on a standard clinical 1.5 T scanner. Three cartesian sequences (fully sampled (10:15 min), accelerated with SENSE (05:08 min), accelerated with Compressed SENSE (03:32 min)) and six different spiral sequences were acquired (spiral acquisition windows ranging from 10 to 5 ms (01:32 min-03:05 min)). Three readers graded the images qualitatively in terms of overall image quality, vessel sharpness, inhomogeneous intraluminal signal, background noise, visualization of large and small vessels and overall impression of the number of visible vessels. Cross-sectional areas of the vessel lumen were measured and vessel sharpness was quantified.

RESULTS

The SENSE and Compressed SENSE accelerated cartesian sequences and the Spiral 6 ms and 5 ms sequences were deemed comparable to the fully sampled cartesian sequence in most qualitative categories (p > 0.05) based on exact binomial tests. The Spiral 6 ms and 5 ms sequences achieved a scan time reduction of 75.3% and 69.9% respectively compared to the fully sampled cartesian sequence. The spiral sequences (generally) exhibited improved subjective vessel sharpness (p < 0.01-p = 0.13) but increased background noise (p = 0.03-p = 0.25). Cross-sectional area measurements were similar between all sequences (Krippendorff's alpha: 0.955-0.982). Quantitative vessel sharpness was increased for all spiral sequences compared to all cartesian sequences (p = 0.004).

CONCLUSIONS

Spiral 3D TOF MRA sequences with a spiral acquisition window of 5 ms or 6 ms may be used for accurate, rapid, clinical non-contrast carotid artery imaging.

T2-Weighted Whole-Brain Intracranial Vessel Wall Imaging at 3 Tesla With Cerebrospinal Fluid Suppression

Background

T₂-weighted (T₂w) intracranial vessel wall imaging (IVWI) provides good contrast to differentiate intracranial vasculopathies and discriminate various important plaque components. However, the strong cerebrospinal fluid (CSF) signal in T₂w images interferes with depicting the intracranial vessel wall. In this study, we propose a T₂-prepared sequence for whole-brain IVWI at 3T with CSF suppression.

Methods

A preparation module that combines T₂ preparation and inversion recovery (T₂IR) was used to suppress the CSF signal and was incorporated into the commercial three-dimensional (3D) turbo spin echo sequence-Sampling Perfection with Application optimized Contrast using different flip angle Evolution (SPACE). This new technique (hereafter called T₂IR-SPACE) was evaluated on nine healthy volunteers and compared with two other commonly used 3D T₂-weighted sequences: T₂w-SPACE and FLAIR-SPACE (FLAIR: fluid-attenuated inversion recovery). The signal-to-noise ratios (SNRs) of the vessel wall (VW) and CSF and contrast-to-noise ratios (CNRs) between them were measured and compared among these three T₂-weighted sequences. Subjective wall visualization of the three T₂-weighted sequences was scored blindly and independently by two radiologists using a four-point scale followed by inter-rater reproducibility analysis. A pilot study of four stroke patients was performed to preliminarily evaluate the diagnostic value of this new sequence, which was compared with two conventional T₂-weighted sequences.

Results

T₂IR-SPACE had the highest CNR (11.01 ± 6.75) compared with FLAIR-SPACE (4.49 ± 3.15; p < 0.001) and T₂w-SPACE (-56.16 ± 18.58; p < 0.001). The subjective wall visualization score of T₂IR-SPACE was higher than those of FLAIR-SPACE and T₂w-SPACE (T₂IR-SPACE: 2.35 ± 0.59; FLAIR-SPACE: 0.52 ± 0.54; T₂w-SPACE: 1.67 ± 0.58); the two radiologists' scores showed excellent agreement (ICC = 0.883).

Conclusion

The T₂IR preparation module markedly suppressed the CSF signal without much SNR loss of the other tissues (i.e., vessel wall, white matter, and gray matter) compared with the IR pulse. Our results suggest that T₂IR-SPACE is a potential alternative T₂-weighted sequence for assessing intracranial vascular diseases.

A novel spectrally selective fat saturation pulse design with robustness to B0 and B1 inhomogeneities: A demonstration on 3D T1-weighted breast MRI at 3 T

PURPOSE

Spectrally selective fat saturation (FatSat) sequence is commonly used to suppress signal from adipose tissue. Conventional SINC-shaped pulses are sensitive to B₀ off-resonance and B₁⁺ offset. Uniform fat saturation with large spatial coverage is especially challenging for the body and breast MRI. The aim of this study is to develop spectrally selective FatSat pulses that offer more immunity to B₀/B₁⁺ field inhomogeneities than SINC pulses and evaluate them in bilateral breast imaging at 3 T.

MATERIALS AND METHODS

Optimized composite pulses (OCP) were designed based on the optimal control theory with robustness to a targeted B₀/ B₁⁺ conditions. OCP pulses also allows flexible flip angles to meet different requirements. Comparisons with the vendor-provided SINC pulses were conducted by numerical simulation and in vivo scans using a 3D T₁-weighted (T₁w) gradient-echo (GRE) sequence with coverage of the whole-breast.

RESULTS

Simulation revealed that OCP pulses yielded almost half of the transition band and much less sensitivity to B₁⁺ inhomogeneity compared to SINC pulses with B₀ off-resonance within ±200 Hz and B₁⁺ scale error within ±0.3 (P < 0.001). Across five normal subjects, OCP FatSat pulses produced 25-41% lower residual fat signals (P < 0.05) with 27-36% less spatial variation (P < 0.05) than SINC.

CONCLUSION

In contrast to conventional SINC-shaped pulses, the newly designed OCP FatSat pulses mitigated challenges of wide range of B₀/ B₁⁺ field inhomogeneities and achieved more uniform fat suppression in bilateral breast T₁w imaging at 3 T.