Echo planar time-resolved imaging (EPTI)

Distortion-free water-fat separated diffusion-weighted imaging using spatiotemporal joint reconstruction

PURPOSE

Diffusion-weighted imaging (DWI) suffers from geometric distortion and chemical shift artifacts due to the commonly used Echo Planar Imaging (EPI) trajectory. Even with fat suppression in DWI, severe B0 and B1 variations can result in residual fat, which becomes both a source of image artifacts and a confounding factor in diffusion-weighted contrast in distinguishing benign and malignant tissues. This work presents a method for acquiring distortion-free diffusion-weighted images using spatiotemporal acquisition and joint reconstruction. Water-fat separation is performed by chemical-shift encoding.

METHODS

Spatiotemporal acquisition is employed to obtain distortion-free images at a series of echo times. Chemical-shift encoding is used for water-fat separation. Reconstruction and separation are performed jointly in the spat-spectral domain. To address the shot-to-shot motion-induced phase in DWI, an Fast Spin Echo (FSE)-based phase navigator is incorporated into the sequence to obtain distortion-free phase information. The proposed method was validated in phantoms and in vivo for the brain, head and neck, and breast.

RESULTS

The proposed method enables the acquisition of distortion-free diffusion-weighted images in the presence of B0 field inhomogenieties commonly observed in the body. Water and fat components are separated with no obvious spectral leakage artifacts. The estimated Apparent Diffusion Coefficient (ADC) is comparable to that of multishot DW-EPI.

CONCLUSION

Distortion-free, water-fat separated diffusion-weighted images in body can be obtained through the utilization of spatiotemporal acquisition and joint reconstruction methods.

Rapid submillimeter QSM and R2* mapping using interleaved multishot 3D-EPI at 7 and 3 Tesla

PURPOSE

To explore the high signal-to-noise ratio (SNR) efficiency of interleaved multishot 3D-EPI with standard image reconstruction for fast and robust high-resolution whole-brain quantitative susceptibility (QSM) andR 2 ∗ $$ {R}_2^{\ast } $$ mapping at 7 and 3T.

METHODS

Single- and multi-TE segmented 3D-EPI is combined with conventional CAIPIRINHA undersampling for up to 72-fold effective gradient echo (GRE) imaging acceleration. Across multiple averages, scan parameters are varied (e.g., dual-polarity frequency-encoding) to additionally correct forB 0 $$ {\mathrm{B}}_0 $$ -induced artifacts, geometric distortions and motion retrospectively. A comparison to established GRE protocols is made. Resolutions range from 1.4 mm isotropic (1 multi-TE average in 36 s) up to 0.4 mm isotropic (2 single-TE averages in approximately 6 min) with whole-head coverage.

RESULTS

Only 1-4 averages are needed for sufficient SNR with 3D-EPI, depending on resolution and field strength. Fast scanning and small voxels together with retrospective corrections result in substantially reduced image artifacts, which improves susceptibility andR 2 ∗ $$ {R}_2^{\ast } $$ mapping. Additionally, much finer details are obtained in susceptibility-weighted image projections through significantly reduced partial voluming.

CONCLUSION

Using interleaved multishot 3D-EPI, single-TE and multi-TE data can readily be acquired 10 times faster than with conventional, accelerated GRE imaging. Even 0.4 mm isotropic whole-head QSM within 6 min becomes feasible at 7T. At 3T, motion-robust 0.8 mm isotropic whole-brain QSM andR 2 ∗ $$ {R}_2^{\ast } $$ mapping with no apparent distortion in less than 7 min becomes clinically feasible. Stronger gradient systems may allow for even higher effective acceleration rates through larger EPI factors while maintaining optimal contrast.

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.

Reduced physiology-induced temporal instability achieved with variable-flip-angle fast low-angle excitation echo-planar technique with multishot echo planar time-resolved imaging

PURPOSE

Echo planar time-resolved imaging (EPTI) is a new imaging approach that addresses the limitations of EPI by providing high-resolution, distortion- and T2/T 2 * $$ {\mathrm{T}}_2^{\ast } $$  blurring-free imaging for functional MRI (fMRI). However, as in all multishot sequences, intershot phase variations induced by physiological processes can introduce temporal instabilities to the reconstructed time-series data. This study aims to reduce these instabilities in multishot EPTI.

THEORY AND METHODS

In conventional multishot EPTI, the time intervals between the shots comprising each slice can introduce intershot phase variations. Here, the fast low-angle excitation echo-planar technique (FLEET), in which all shots of each slice are acquired consecutively with minimal time delays, was combined with a variable flip angle (VFA) technique to improve intershot consistency and maximize signal. A recursive Shinnar-Le Roux RF pulse design algorithm was used to generate pulses for different shots to produce consistent slice profiles and signal intensities across shots. Blipped controlled aliasing in parallel imaging simultaneous multislice was also combined with the proposed VFA-FLEET EPTI to improve temporal resolution and increase spatial coverage.

RESULTS

The temporal stability of VFA-FLEET EPTI was compared with conventional EPTI at 7 T. The results demonstrated that VFA-FLEET can provide spatial-specific increase of temporal stability. We performed high-resolution task-fMRI experiments at 7 T using VFA-FLEET EPTI, and reliable BOLD responses to a visual stimulus were detected.

CONCLUSION

The intershot phase variations induced by physiological processes in multishot EPTI can manifest as specific spatial patterns of physiological noise enhancement and lead to reduced temporal stability. The VFA-FLEET technique can substantially reduce these physiology-induced instabilities in multishot EPTI acquisitions. The proposed method provides sufficient stability and sensitivity for high-resolution fMRI studies.

Spherical echo-planar time-resolved imaging (sEPTI) for rapid 3D quantitative T 2 * and susceptibility imaging

PURPOSE

To develop a 3D spherical EPTI (sEPTI) acquisition and a comprehensive reconstruction pipeline for rapid high-quality whole-brain submillimeterT 2 * $$ {\mathrm{T}}_2^{\ast } $$ and QSM quantification.

METHODS

For the sEPTI acquisition, spherical k-space coverage is utilized with variable echo-spacing and maximum kx ramp-sampling to improve efficiency and signal incoherency compared to existing EPTI approaches. For reconstruction, an iterative rank-shrinking B0 estimation and odd-even high-order phase correction algorithms were incorporated into the reconstruction to better mitigate artifacts from field imperfections. A physics-informed unrolled network was utilized to boost the SNR, where 1-mm and 0.75-mm isotropic whole-brain imaging were performed in 45 and 90 s at 3 T, respectively. These protocols were validated through simulations, phantom, and in vivo experiments. Ten healthy subjects were recruited to provide sufficient data for the unrolled network. The entire pipeline was validated on additional five healthy subjects where different EPTI sampling approaches were compared. Two additional pediatric patients with epilepsy were recruited to demonstrate the generalizability of the unrolled reconstruction.

RESULTS

sEPTI achieved 1.4× $$ \times $$ faster imaging with improved image quality and quantitative map precision compared to existing EPTI approaches. The B0 update and the phase correction provide improved reconstruction performance with lower artifacts. The unrolled network boosted the SNR, achieving high-qualityT 2 * $$ {\mathrm{T}}_2^{\ast } $$ and QSM quantification with single average data. High-quality reconstruction was also obtained in the pediatric patients using this network.

CONCLUSION

sEPTI achieved whole-brain distortion-free multi-echo imaging andT 2 * $$ {\mathrm{T}}_2^{\ast } $$ and QSM quantification at 0.75 mm in 90 s which has the potential to be useful for wide clinical applications.

QRAGE-Simultaneous multiparametric quantitative MRI of water content, T1, T2*, and magnetic susceptibility at ultrahigh field strength

PURPOSE

To introduce quantitative rapid gradient-echo (QRAGE), a novel approach for the simultaneous mapping of multiple quantitative MRI parameters, including water content, T1, T2*, and magnetic susceptibility at ultrahigh field strength.

METHODS

QRAGE leverages a newly developed multi-echo MPnRAGE sequence, facilitating the acquisition of 171 distinct contrast images across a range of inversion and TE points. To maintain a short acquisition time, we introduce MIRAGE2, a novel model-based reconstruction method that exploits prior knowledge of temporal signal evolution, represented as damped complex exponentials. MIRAGE2 minimizes local Block-Hankel and Casorati matrices. Parameter maps are derived from the reconstructed contrast images through postprocessing steps. We validate QRAGE through extensive simulations, phantom studies, and in vivo experiments, demonstrating its capability for high-precision imaging.

RESULTS

In vivo brain measurements show the promising performance of QRAGE, with test-retest SDs and deviations from reference methods of < 0.8% for water content, < 17 ms for T1, and < 0.7 ms for T2*. QRAGE achieves whole-brain coverage at a 1-mm isotropic resolution in just 7 min and 15 s, comparable to the acquisition time of an MP2RAGE scan. In addition, QRAGE generates a contrast image akin to the UNI image produced by MP2RAGE.

CONCLUSION

QRAGE is a new, successful approach for simultaneously mapping multiple MR parameters at ultrahigh field.

Ultrafast Brain MRI at 3 T for MS: Evaluation of a 51-Second Deep Learning-Enhanced T2-EPI-FLAIR Sequence

In neuroimaging, there is no equivalent alternative to magnetic resonance imaging (MRI). However, image acquisitions are generally time-consuming, which may limit utilization in some cases, e.g., in patients who cannot remain motionless for long or suffer from claustrophobia, or in the event of extensive waiting times. For multiple sclerosis (MS) patients, MRI plays a major role in drug therapy decision-making. The purpose of this study was to evaluate whether an ultrafast, T2-weighted (T2w), deep learning-enhanced (DL), echo-planar-imaging-based (EPI) fluid-attenuated inversion recovery (FLAIR) sequence (FLAIRUF) that has targeted neurological emergencies so far might even be an option to detect MS lesions of the brain compared to conventional FLAIR sequences. Therefore, 17 MS patients were enrolled prospectively in this exploratory study. Standard MRI protocols and ultrafast acquisitions were conducted at 3 tesla (T), including three-dimensional (3D)-FLAIR, turbo/fast spin-echo (TSE)-FLAIR, and FLAIRUF. Inflammatory lesions were grouped by size and location. Lesion conspicuity and image quality were rated on an ordinal five-point Likert scale, and lesion detection rates were calculated. Statistical analyses were performed to compare results. Altogether, 568 different lesions were found. Data indicated no significant differences in lesion detection (sensitivity and positive predictive value [PPV]) between FLAIRUF and axially reconstructed 3D-FLAIR (lesion size ≥3 mm × ≥2 mm) and no differences in sensitivity between FLAIRUF and TSE-FLAIR (lesion size ≥3 mm total). Lesion conspicuity in FLAIRUF was similar in all brain regions except for superior conspicuity in the occipital lobe and inferior conspicuity in the central brain regions. Further findings include location-dependent limitations of signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) as well as artifacts such as spatial distortions in FLAIRUF. In conclusion, FLAIRUF could potentially be an expedient alternative to conventional methods for brain imaging in MS patients since the acquisition can be performed in a fraction of time while maintaining good image quality.

Dynamic field mapping and distortion correction using single-shot blip-rewound EPI and joint multi-echo reconstruction

PURPOSE

To develop a method for dynamic∆ B 0 $$ \Delta {B}_0 $$ mapping and distortion correction.

METHODS

A blip-rewound EPI trajectory was developed to acquire multiple 2D EPI images in a single readout with an interleaved order, which allows a short TE difference. A joint multi-echo reconstruction was utilized to exploit the shared information between EPI images. The reconstructed images from each readout are combined to produce a final magnitude image. A∆ B 0 $$ \Delta {B}_0 $$ map is calculated from the phase of these images for distortion correction. The efficacy of the proposed method is assessed with phantom and in vivo experiments. The performance of the proposed method in the presence of subject motion is also investigated.

RESULTS

Compared to conventional multi-echo EPI, the proposed method allows dynamic∆ B 0 $$ \Delta {B}_0 $$ mapping at matched resolution with a much shorter TR. Phantom and in vivo results show that the proposed method can provide a comparable magnitude image as conventional single-shot EPI. The∆ B 0 $$ \Delta {B}_0 $$ maps calculated from the proposed method are consistent with conventional multi-echo EPI in the phantom experiment. For in vivo experiments, the proposed method provides a more accurate estimation of∆ B 0 $$ \Delta {B}_0 $$ than conventional multi-echo EPI, which is prone to phase wrapping problems due to the long TE difference. In-vivo scan with subject motion shows the proposed dynamic field mapping method can improve the temporal stability of EPI time series compared to gradient echo (GRE) based static field mapping.

CONCLUSION

The proposed method allows accurate dynamic∆ B 0 $$ \Delta {B}_0 $$ mapping for robust distortion correction without compromising spatial or temporal resolution.

Multi-Parametric Molecular Imaging of the Brain Using Optimized Multi-TE Subspace MRSI

OBJECTIVE

To develop a novel multi-TE MR spectroscopic imaging (MRSI) approach to enable label-free, simultaneous, high-resolution mapping of several molecules and their biophysical parameters in the brain.

METHODS

The proposed method uniquely integrated an augmented molecular-component-specific subspace model for multi-TE 1H-MRSI signals, an estimation-theoretic experiment optimization (nonuniform TE selection) for molecule separation and parameter estimation, a physics-driven subspace learning strategy for spatiospectral reconstruction and molecular quantification, and a new accelerated multi-TE MRSI acquisition for generating high-resolution data in clinically relevant times. Numerical studies, phantom and in vivo experiments were conducted to validate the optimized experiment design and demonstrate the imaging capability offered by the proposed method.

RESULTS

The proposed TE optimization improved estimation of metabolites, neurotransmitters and their T2's over conventional TE choices, e.g., reducing variances of neurotransmitter concentration by  ∼  40% and metabolite T2 by  ∼  60%. Simultaneous metabolite and neurotransmitter mapping of the brain can be achieved at a nominal resolution of 3.4 × 3.4 × 6.4 mm 3. High-resolution, 3D metabolite T2 mapping was made possible for the first time. The translational potential of the proposed method was demonstrated by mapping biochemical abnormality in a post-traumatic epilepsy (PTE) patient.

CONCLUSION

The feasibility for high-resolution mapping of metabolites/neurotransmitters and metabolite T2's within clinically relevant time was demonstrated. We expect our method to offer richer information for revealing and understanding metabolic alterations in neurological diseases.

SIGNIFICANCE

A novel multi-TE MRSI approach was presented that enhanced the technological capability of multi-parametric molecular imaging of the brain. The proposed method presents new technology development and application opportunities for providing richer molecular level information to uncover and comprehend metabolic changes relevant in various neurological applications.

Model-based iterative reconstruction for direct imaging with point spread function encoded echo planar MRI

BACKGROUND

Echo planar imaging (EPI) is a fast measurement technique commonly used in magnetic resonance imaging (MRI), but is highly sensitive to measurement non-idealities in reconstruction. Point spread function (PSF)-encoded EPI is a multi-shot strategy which alleviates distortion, but acquisition of encodings suitable for direct distortion-free imaging prolongs scan time. In this work, a model-based iterative reconstruction (MBIR) framework is introduced for direct imaging with PSF-EPI to improve image quality and acceleration potential.

METHODS

An MBIR platform was developed for accelerated PSF-EPI. The reconstruction utilizes a subspace representation, is regularized to promote local low-rankedness (LLR), and uses variable splitting for efficient iteration. Comparisons were made against standard reconstructions from prospectively accelerated PSF-EPI data and with retrospective subsampling. Exploring aggressive partial Fourier acceleration of the PSF-encoding dimension, additional comparisons were made against an extension of Homodyne to direct PSF-EPI in numerical experiments. A neuroradiologists' assessment was completed comparing images reconstructed with MBIR from retrospectively truncated data directly against images obtained with standard reconstructions from non-truncated datasets.

RESULTS

Image quality results were consistently superior for MBIR relative to standard and Homodyne reconstructions. As the MBIR signal model and reconstruction allow for arbitrary sampling of the PSF space, random sampling of the PSF-encoding dimension was also demonstrated, with quantitative assessments indicating best performance achieved through nonuniform PSF sampling combined with partial Fourier. With retrospective subsampling, MBIR reconstructs high-quality images from sub-minute scan datasets. MBIR was shown to be superior in a neuroradiologists' assessment with respect to three of five performance criteria, with equivalence for the remaining two.

CONCLUSIONS

A novel image reconstruction framework is introduced for direct imaging with PSF-EPI, enabling arbitrary PSF space sampling and reconstruction of diagnostic-quality images from highly accelerated PSF-encoded EPI data.

Phase-sensitive deep reconstruction method for rapid multiparametric MR fingerprinting and quantitative susceptibility mapping in the brain

INTRODUCTION

This study explores the potential of Magnetic Resonance Fingerprinting (MRF) with a novel Phase-Sensitivity Deep Reconstruction Network (PS-DRONE) for simultaneous quantification of T1, T2, Proton Density, B1+, phase and quantitative susceptibility mapping (QSM).

METHODS

Data were acquired at 3 T in vitro and in vivo using an optimized EPI-based MRF sequence. Phantom experiments were conducted using a standardized phantom for T1 and T2 maps and a custom-made agar-based gadolinium phantom for B1 and QSM maps. In vivo experiments included five healthy volunteers and one patient diagnosed with brain metastasis. PSDRONE maps were compared to reference maps obtained through standard imaging sequences.

RESULTS

Total scan time was 2 min for 32 slices and a resolution of [1 mm, 1 mm, 4.5 mm]. The reconstruction of T1, T2, Proton Density, B1+ and phase maps were reconstructed within 1 s. In the phantoms, PS-DRONE analysis presented accurate and strongly correlated T1 and T2 maps (r = 0.99) compared to the reference maps. B1 maps from PS-DRONE showed slightly higher values, though still correlated (r = 0.6) with the reference. QSM values showed a small bias but were strongly correlated (r = 0.99) with reference data. In the in vivo analysis, PS-DRONE-derived T1 and T2 values for gray and white matter matched reference values in healthy volunteers. PS-DRONE B1 and QSM maps showed strong correlations with reference values.

CONCLUSION

The PS-DRONE network enables concurrent acquisition of T1, T2, PD, B1+, phase and QSM maps, within 2 min of acquisition time and 1 s of reconstruction time.

Zero-DeepSub: Zero-shot deep subspace reconstruction for rapid multiparametric quantitative MRI using 3D-QALAS

PURPOSE

To develop and evaluate methods for (1) reconstructing 3D-quantification using an interleaved Look-Locker acquisition sequence with T2 preparation pulse (3D-QALAS) time-series images using a low-rank subspace method, which enables accurate and rapid T1 and T2 mapping, and (2) improving the fidelity of subspace QALAS by combining scan-specific deep-learning-based reconstruction and subspace modeling.

THEORY AND METHODS

A low-rank subspace method for 3D-QALAS (i.e., subspace QALAS) and zero-shot deep-learning subspace method (i.e., Zero-DeepSub) were proposed for rapid and high fidelity T1 and T2 mapping and time-resolved imaging using 3D-QALAS. Using an ISMRM/NIST system phantom, the accuracy and reproducibility of the T1 and T2 maps estimated using the proposed methods were evaluated by comparing them with reference techniques. The reconstruction performance of the proposed subspace QALAS using Zero-DeepSub was evaluated in vivo and compared with conventional QALAS at high reduction factors of up to nine-fold.

RESULTS

Phantom experiments showed that subspace QALAS had good linearity with respect to the reference methods while reducing biases and improving precision compared to conventional QALAS, especially for T2 maps. Moreover, in vivo results demonstrated that subspace QALAS had better g-factor maps and could reduce voxel blurring, noise, and artifacts compared to conventional QALAS and showed robust performance at up to nine-fold acceleration with Zero-DeepSub, which enabled whole-brain T1, T2, and PD mapping at 1 mm isotropic resolution within 2 min of scan time.

CONCLUSION

The proposed subspace QALAS along with Zero-DeepSub enabled high fidelity and rapid whole-brain multiparametric quantification and time-resolved imaging.

Romer-EPTI: rotating-view motion-robust super-resolution EPTI for SNR-efficient distortion-free in-vivo mesoscale dMRI and microstructure imaging

Purpose

To overcome the major challenges in dMRI acquisition, including low SNR, distortion/blurring, and motion vulnerability.

Methods

A novel Romer-EPTI technique is developed to provide distortion-free dMRI with significant SNR gain, high motion-robustness, sharp spatial resolution, and simultaneous multi-TE imaging. It introduces a ROtating-view Motion-robust supEr-Resolution technique (Romer) combined with a distortion/blurring-free EPTI encoding. Romer enhances SNR by a simultaneous multi-thick-slice acquisition with rotating-view encoding, while providing high motion-robustness through a motion-aware super-resolution reconstruction, which also incorporates slice-profile and real-value diffusion, to resolve high-isotropic-resolution volumes. The in-plane encoding is performed using distortion/blurring-free EPTI, which further improves effective spatial resolution and motion robustness by preventing not only T2/T2*-blurring but also additional blurring resulting from combining encoded volumes with inconsistent geometries caused by dynamic distortions. Self-navigation was incorporated to enable efficient phase correction. Additional developments include strategies to address slab-boundary artifacts, achieve minimal TE for SNR gain at 7T, and achieve high robustness to strong phase variations at high b-values.

Results

Using Romer-EPTI, we demonstrate distortion-free whole-brain mesoscale in-vivo dMRI at both 3T (500-μm-iso) and 7T (485-μm-iso) for the first time, with high SNR efficiency (e.g., 25 × ), and high image quality free from distortion and slab-boundary artifacts with minimal blurring. Motion experiments demonstrate Romer-EPTI's high motion-robustness and ability to recover sharp images in the presence of motion. Romer-EPTI also demonstrates significant SNR gain and robustness in high b-value (b=5000s/mm2) and time-dependent dMRI.

Conclusion

Romer-EPTI significantly improves SNR, motion-robustness, and image quality, providing a highly efficient acquisition for high-resolution dMRI and microstructure imaging.

Neuroimaging of Parkinson's disease by quantitative susceptibility mapping

Parkinson's disease (PD) is a common neurodegenerative disease, and apart from a few rare genetic causes, its pathogenesis remains largely unclear. Recent scientific interest has been captured by the involvement of iron biochemistry and the disruption of iron homeostasis, particularly within the brain regions specifically affected in PD. The advent of Quantitative Susceptibility Mapping (QSM) has enabled non-invasive quantification of brain iron in vivo by MRI, which has contributed to the understanding of iron-associated pathogenesis and has the potential for the development of iron-based biomarkers in PD. This review elucidates the biochemical underpinnings of brain iron accumulation, details advancements in iron-sensitive MRI technologies, and discusses the role of QSM as a biomarker of iron deposition in PD. Despite considerable progress, several challenges impede its clinical application after a decade of QSM studies. The initiation of multi-site research is warranted for developing robust, interpretable, and disease-specific biomarkers for monitoring PD disease progression.

Correcting synthetic MRI contrast-weighted images using deep learning

Synthetic magnetic resonance imaging (MRI) offers a scanning paradigm where a fast multi-contrast sequence can be used to estimate underlying quantitative tissue parameter maps, which are then used to synthesize any desirable clinical contrast by retrospectively changing scan parameters in silico. Two benefits of this approach are the reduced exam time and the ability to generate arbitrary contrasts offline. However, synthetically generated contrasts are known to deviate from the contrast of experimental scans. The reason for contrast mismatch is the necessary exclusion of some unmodeled physical effects such as partial voluming, diffusion, flow, susceptibility, magnetization transfer, and more. The inclusion of these effects in signal encoding would improve the synthetic images, but would make the quantitative imaging protocol impractical due to long scan times. Therefore, in this work, we propose a novel deep learning approach that generates a multiplicative correction term to capture unmodeled effects and correct the synthetic contrast images to better match experimental contrasts for arbitrary scan parameters. The physics inspired deep learning model implicitly accounts for some unmodeled physical effects occurring during the scan. As a proof of principle, we validate our approach on synthesizing arbitrary inversion recovery fast spin-echo scans using a commercially available 2D multi-contrast sequence. We observe that the proposed correction visually and numerically reduces the mismatch with experimentally collected contrasts compared to conventional synthetic MRI. Finally, we show results of a preliminary reader study and find that the proposed method statistically significantly improves in contrast and SNR as compared to synthetic MR images.

Single-shot Echo Planar Time-resolved Imaging for multi-echo functional MRI and distortion-free diffusion imaging

Purpose

To develop EPTI, a multi-shot distortion-free multi-echo imaging technique, into a single-shot acquisition to achieve improved robustness to motion and physiological noise, increased temporal resolution, and high SNR efficiency for dynamic imaging applications.

Methods

A new spatiotemporal encoding was developed to achieve single-shot EPTI by enhancing spatiotemporal correlation in k-t space. The proposed single-shot encoding improves reconstruction conditioning and sampling efficiency, with additional optimization under various accelerations to achieve optimized performance. To achieve high SNR efficiency, continuous readout with minimized deadtime was employed that begins immediately after excitation and extends for an SNR-optimized length. Moreover, k-t partial Fourier and simultaneous multi-slice acquisition were integrated to further accelerate the acquisition and achieve high spatial and temporal resolution.

Results

We demonstrated that ss-EPTI achieves higher tSNR efficiency than multi-shot EPTI, and provides distortion-free imaging with densely-sampled multi-echo images at resolutions ~1.25-3 mm at 3T and 7T-with high SNR efficiency and with comparable temporal resolutions to ss-EPI. The ability of ss-EPTI to eliminate dynamic distortions common in EPI also further improves temporal stability. For fMRI, ss-EPTI also provides early-TE images (e.g., 2.9ms) to recover signal-intensity and functional-sensitivity dropout in challenging regions. The multi-echo images provide TE-dependent information about functional fluctuations, successfully distinguishing noise-components from BOLD signals and further improving tSNR. For diffusion MRI, ss-EPTI provides high-quality distortion-free diffusion images and multi-echo diffusion metrics.

Conclusion

ss-EPTI provides distortion-free imaging with high image quality, rich multi-echo information, and enhanced efficiency within comparable temporal resolution to ss-EPI, offering a robust and efficient acquisition for dynamic imaging.

mcLARO: Multi-contrast learned acquisition and reconstruction optimization for simultaneous quantitative multi-parametric mapping

PURPOSE

To develop a method for rapid sub-millimeter T1 , T2 ,T 2 * $$ {\mathrm{T}}_2^{\ast } $$ , and QSM mapping in a single scan using multi-contrast learned acquisition and reconstruction optimization (mcLARO).

METHODS

A pulse sequence was developed by interleaving inversion recovery and T2 magnetization preparations and single-echo and multi-echo gradient echo acquisitions, which sensitized k-space data to T1 , T2 ,T 2 * $$ {\mathrm{T}}_2^{\ast } $$ , and magnetic susceptibility. The proposed mcLARO optimized both the multi-contrast k-space under-sampling pattern and image reconstruction based on image feature fusion using a deep learning framework. The proposed mcLARO method withR = 8 $$ R=8 $$ under-sampling was validated in a retrospective ablation study and compared with other deep learning reconstruction methods, including MoDL and Wave-MoDL, using fully sampled data as reference. Various under-sampling ratios in mcLARO were investigated. mcLARO was also evaluated in a prospective study using separately acquired conventionally sampled quantitative maps as reference standard.

RESULTS

The retrospective ablation study showed improved image sharpness of mcLARO compared to the baseline network without the multi-contrast sampling pattern optimization or image feature fusion module. The under-sampling ratio experiment showed that higher under-sampling ratios resulted in blurrier images and lower precision of quantitative values. The prospective study showed that small or negligible bias and narrow 95% limits of agreement on regional T1 , T2 ,T 2 * $$ {\mathrm{T}}_2^{\ast } $$ , and QSM values by mcLARO (5:39 mins) compared to reference scans (40:03 mins in total). mcLARO outperformed MoDL and Wave-MoDL in terms of image sharpness, noise suppression, and artifact removal.

CONCLUSION

mcLARO enabled fast sub-millimeter T1 , T2 ,T 2 * $$ {\mathrm{T}}_2^{\ast } $$ , and QSM mapping in a single scan.

Ultrafast Brain MRI Protocol at 1.5 T Using Deep Learning and Multi-shot EPI

RATIONALE AND OBJECTIVES

To evaluate clinical feasibility and image quality of a comprehensive ultrafast brain MRI protocol with multi-shot echo planar imaging and deep learning-enhanced reconstruction at 1.5T.

MATERIALS AND METHODS

Thirty consecutive patients who underwent clinically indicated MRI at a 1.5 T scanner were prospectively included. A conventional MRI (c-MRI) protocol, including T1-, T2-, T2*-, T2-FLAIR, and diffusion-weighted images (DWI)-weighted sequences were acquired. In addition, ultrafast brain imaging with deep learning-enhanced reconstruction and multi-shot EPI (DLe-MRI) was performed. Subjective image quality was evaluated by three readers using a 4-point Likert scale. To assess interrater agreement, Fleiss' kappa (ϰ) was determined. For objective image analysis, relative signal intensity levels for grey matter, white matter, and cerebrospinal fluid were calculated.

RESULTS

Time of acquisition (TA) of c-MRI protocols added up to 13:55 minutes, whereas the TA of DLe-MRI-based protocol added up to 3:04 minutes, resulting in a time reduction of 78%. All DLe-MRI acquisitions yielded diagnostic image quality with good absolute values for subjective image quality. C-MRI demonstrated slight advantages for DWI in overall subjective image quality (c-MRI: 3.93 [+/- 0.25] vs DLe-MRI: 3.87 [+/- 0.37], P = .04) and diagnostic confidence (c-MRI: 3.93 [+/- 0.25] vs DLe-MRI: 3.83 [+/- 3.83], P = .01). For most evaluated quality scores, moderate interobserver agreement was found. Objective image evaluation revealed comparable results for both techniques.

CONCLUSION

DLe-MRI is feasible and allows for highly accelerated comprehensive brain MRI within 3minutes at 1.5 T with good image quality. This technique may potentially strengthen the role of MRI in neurological emergencies.

Cramér-Rao Bound Optimized Subspace Reconstruction in Quantitative MRI

We extend the traditional framework for estimating subspace bases that maximize the preserved signal energy to additionally preserve the Cramér-Rao bound (CRB) of the biophysical parameters and, ultimately, improve accuracy and precision in the quantitative maps. To this end, we introduce an approximate compressed CRB based on orthogonalized versions of the signal's derivatives with respect to the model parameters. This approximation permits singular value decomposition (SVD)-based minimization of both the CRB and signal losses during compression. Compared to the traditional SVD approach, the proposed method better preserves the CRB across all biophysical parameters with negligible cost to the preserved signal energy, leading to reduced bias and variance of the parameter estimates in simulation. In vivo, improved accuracy and precision are observed in two quantitative neuroimaging applications, permitting the use of smaller basis sizes in subspace reconstruction and offering significant computational savings.

High-fidelity direct contrast synthesis from magnetic resonance fingerprinting

PURPOSE

This work was aimed at proposing a supervised learning-based method that directly synthesizes contrast-weighted images from the Magnetic Resonance Fingerprinting (MRF) data without performing quantitative mapping and spin-dynamics simulations.

METHODS

To implement our direct contrast synthesis (DCS) method, we deploy a conditional generative adversarial network (GAN) framework with a multi-branch U-Net as the generator and a multilayer CNN (PatchGAN) as the discriminator. We refer to our proposed approach as N-DCSNet. The input MRF data are used to directly synthesize T1-weighted, T2-weighted, and fluid-attenuated inversion recovery (FLAIR) images through supervised training on paired MRF and target spin echo-based contrast-weighted scans. The performance of our proposed method is demonstrated on in vivo MRF scans from healthy volunteers. Quantitative metrics, including normalized root mean square error (nRMSE), peak signal-to-noise ratio (PSNR), structural similarity (SSIM), learned perceptual image patch similarity (LPIPS), and Fréchet inception distance (FID), were used to evaluate the performance of the proposed method and compare it with others.

RESULTS

In-vivo experiments demonstrated excellent image quality with respect to that of simulation-based contrast synthesis and previous DCS methods, both visually and according to quantitative metrics. We also demonstrate cases in which our trained model is able to mitigate the in-flow and spiral off-resonance artifacts typically seen in MRF reconstructions, and thus more faithfully represent conventional spin echo-based contrast-weighted images.

CONCLUSION

We present N-DCSNet to directly synthesize high-fidelity multicontrast MR images from a single MRF acquisition. This method can significantly decrease examination time. By directly training a network to generate contrast-weighted images, our method does not require any model-based simulation and therefore can avoid reconstruction errors due to dictionary matching and contrast simulation (code available at:https://github.com/mikgroup/DCSNet).

Latent signal models: Learning compact representations of signal evolution for improved time-resolved, multi-contrast MRI

PURPOSE

To improve time-resolved reconstructions by training auto-encoders to learn compact representations of Bloch-simulated signal evolution and inserting the decoder into the forward model.

METHODS

Building on model-based nonlinear and linear subspace techniques, we train auto-encoders on dictionaries of simulated signal evolution to learn compact, nonlinear, latent representations. The proposed latent signal model framework inserts the decoder portion of the auto-encoder into the forward model and directly reconstructs the latent representation. Latent signal models essentially serve as a proxy for fast and feasible differentiation through the Bloch equations used to simulate signal. This work performs experiments in the context of T2 -shuffling, gradient echo EPTI, and MPRAGE-shuffling. We compare how efficiently auto-encoders represent signal evolution in comparison to linear subspaces. Simulation and in vivo experiments then evaluate if reducing degrees of freedom by incorporating our proxy for the Bloch equations, the decoder portion of the auto-encoder, into the forward model improves reconstructions in comparison to subspace constraints.

RESULTS

An auto-encoder with 1 real latent variable represents single-tissue fast spin echo, EPTI, and MPRAGE signal evolution to within 0.15% normalized RMS error, enabling reconstruction problems with 3 degrees of freedom per voxel (real latent variable + complex scaling) in comparison to linear models with 4-8 degrees of freedom per voxel. In simulated/in vivo T2 -shuffling and in vivo EPTI experiments, the proposed framework achieves consistent quantitative normalized RMS error improvement over linear approaches. From qualitative evaluation, the proposed approach yields images with reduced blurring and noise amplification in MPRAGE-shuffling experiments.

CONCLUSION

Directly solving for nonlinear latent representations of signal evolution improves time-resolved MRI reconstructions.

High-fidelity mesoscale in-vivo diffusion MRI through gSlider-BUDA and circular EPI with S-LORAKS reconstruction

PURPOSE

To develop a high-fidelity diffusion MRI acquisition and reconstruction framework with reduced echo-train-length for less T2* image blurring compared to typical highly accelerated echo-planar imaging (EPI) acquisitions at sub-millimeter isotropic resolution.

METHODS

We first proposed a circular-EPI trajectory with partial Fourier sampling on both the readout and phase-encoding directions to minimize the echo-train-length and echo time. We then utilized this trajectory in an interleaved two-shot EPI acquisition with reversed phase-encoding polarity, to aid in the correction of off-resonance-induced image distortions and provide complementary k-space coverage in the missing partial Fourier regions. Using model-based reconstruction with structured low-rank constraint and smooth phase prior, we corrected the shot-to-shot phase variations across the two shots and recover the missing k-space data. Finally, we combined the proposed acquisition/reconstruction framework with an SNR-efficient RF-encoded simultaneous multi-slab technique, termed gSlider, to achieve high-fidelity 720 µm and 500 µm isotropic resolution in-vivo diffusion MRI.

RESULTS

Both simulation and in-vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide distortion-corrected diffusion imaging at the mesoscale with markedly reduced T2*-blurring. The in-vivo results of 720 µm and 500 µm datasets show high-fidelity diffusion images with reduced image blurring and echo time using the proposed approaches.

CONCLUSIONS

The proposed method provides high-quality distortion-corrected diffusion-weighted images with ∼40% reduction in the echo-train-length and T2* blurring at 500µm-isotropic-resolution compared to standard multi-shot EPI.

Pseudo multishot echo-planar imaging for geometric distortion improvement

Conventional echo-planar imaging (EPI) uses a radiofrequency pulse for excitation and a prolonged echo train to sample k space, while off-resonance and T₂ * decay effects caused by magnetic susceptibility variation accumulate within each echo, leading to geometric distortion. Multishot EPI methods, which divide k space into segments, can shorten the effective echo spacing and reduce the distortion on EPI images. But multiple shots cost longer scan time and render susceptibility to motion. In this study, we propose a new "multishot" EPI method termed pseudo multishot EPI (pmsEPI), in which phase-encoding lines are segmented as in multishot EPI but are collected within a single shot. With the magnetization divided into different pathways via interleaved excitation instead of refocusing in a single long echo train, the total phase error accumulation is reduced in each segmented acquisition, thereby improving distortion of the resultant EPI image. The performance of the pmsEPI method is demonstrated by phantom and in vivo brain experiments on a 3-T scanner. The experimental results show that the distortion displacements of pmsEPI acquisition compared with conventional EPI decrease by 50% with two pseudo shots and 66% with three pseudo shots, validating the ability of the method to obtain images with reduced distortion in a single shot, although magnetization splitting may induce more than 40% SNR loss and minor artifacts. Specifically, the ability of pmsEPI in diffusion-weighted imaging with different trajectory options is highlighted, and the flexibility is demonstrated in a single-shot blip up and down acquisition.

3D-EPI blip-up/down acquisition (BUDA) with CAIPI and joint Hankel structured low-rank reconstruction for rapid distortion-free high-resolution T 2 * mapping

PURPOSE

This work aims to develop a novel distortion-free 3D-EPI acquisition and image reconstruction technique for fast and robust, high-resolution, whole-brain imaging as well as quantitativeT 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping.

METHODS

3D Blip-up and -down acquisition (3D-BUDA) sequence is designed for both single- and multi-echo 3D gradient recalled echo (GRE)-EPI imaging using multiple shots with blip-up and -down readouts to encode B0 field map information. Complementary k-space coverage is achieved using controlled aliasing in parallel imaging (CAIPI) sampling across the shots. For image reconstruction, an iterative hard-thresholding algorithm is employed to minimize the cost function that combines field map information informed parallel imaging with the structured low-rank constraint for multi-shot 3D-BUDA data. Extending 3D-BUDA to multi-echo imaging permitsT 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping. For this, we propose constructing a joint Hankel matrix along both echo and shot dimensions to improve the reconstruction.

RESULTS

Experimental results on in vivo multi-echo data demonstrate that, by performing joint reconstruction along with both echo and shot dimensions, reconstruction accuracy is improved compared to standard 3D-BUDA reconstruction. CAIPI sampling is further shown to enhance image quality. ForT 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping, parameter values from 3D-Joint-CAIPI-BUDA and reference multi-echo GRE are within limits of agreement as quantified by Bland-Altman analysis.

CONCLUSIONS

The proposed technique enables rapid 3D distortion-free high-resolution imaging andT 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping. Specifically, 3D-BUDA enables 1-mm isotropic whole-brain imaging in 22 s at 3T and 9 s on a 7T scanner. The combination of multi-echo 3D-BUDA with CAIPI acquisition and joint reconstruction enables distortion-free whole-brainT 2 * $$ {\mathrm{T}}_2^{\ast } $$ mapping in 47 s at 1.1 × 1.1 × 1.0 mm3 resolution.

APART-QSM: an improved sub-voxel quantitative susceptibility mapping for susceptibility source separation using an iterative data fitting method

The brain tissue phase contrast in MRI sequences reflects the spatial distributions of multiple substances, such as iron, myelin, calcium, and proteins. These substances with paramagnetic and diamagnetic susceptibilities often colocalize in one voxel in brain regions. Both opposing susceptibilities play vital roles in brain development and neurodegenerative diseases. Conventional QSM methods only provide voxel-averaged susceptibility value and cannot disentangle intravoxel susceptibilities with opposite signs. Advanced susceptibility imaging methods have been recently developed to distinguish the contributions of opposing susceptibility sources for QSM. The basic concept of separating paramagnetic and diamagnetic susceptibility proportions is to include the relaxation rate R₂^* with R₂^' in QSM. The magnitude decay kernel, describing the proportionality coefficient between R₂^' and susceptibility, is an essential reconstruction coefficient for QSM separation methods. In this study, we proposed a more comprehensive complex signal model that describes the relationship between 3D GRE signal and the contributions of paramagnetic and diamagnetic susceptibility to the frequency shift and R₂^* relaxation. The algorithm is implemented as a constrained minimization problem in which the voxel-wise magnitude decay kernel and sub-voxel susceptibilities are determined alternately in each iteration until convergence. The calculated voxel-wise magnitude decay kernel could realistically model the relationship between the R₂^' relaxation and the volume susceptibility. Thus, the proposed method effectively prevents the errors of the magnitude decay kernel from propagating to the final susceptibility separation reconstruction. Phantom studies, ex vivo macaque brain experiments, and in vivo human brain imaging studies were conducted to evaluate the ability of the proposed method to distinguish paramagnetic and diamagnetic susceptibility sources. The results demonstrate that the proposed method provides state-of-the-art performances for quantifying brain iron and myelin compared to previous QSM separation methods. Our results show that the proposed method has the potential to simultaneously quantify whole brain iron and myelin during brain development and aging. The proposed model was also deployed with multiple-orientation complex GRE data input measurements, resulting in high-quality QSM separation maps with more faithful tissue delineation between brain structures compared to those reconstructed by single-orientation QSM separation methods.

A vision of 14 T MR for fundamental and clinical science

OBJECTIVE

We outline our vision for a 14 Tesla MR system. This comprises a novel whole-body magnet design utilizing high temperature superconductor; a console and associated electronic equipment; an optimized radiofrequency coil setup for proton measurement in the brain, which also has a local shim capability; and a high-performance gradient set.

RESEARCH FIELDS

The 14 Tesla system can be considered a 'mesocope': a device capable of measuring on biologically relevant scales. In neuroscience the increased spatial resolution will anatomically resolve all layers of the cortex, cerebellum, subcortical structures, and inner nuclei. Spectroscopic imaging will simultaneously measure excitatory and inhibitory activity, characterizing the excitation/inhibition balance of neural circuits. In medical research (including brain disorders) we will visualize fine-grained patterns of structural abnormalities and relate these changes to functional and molecular changes. The significantly increased spectral resolution will make it possible to detect (dynamic changes in) individual metabolites associated with pathological pathways including molecular interactions and dynamic disease processes.

CONCLUSIONS

The 14 Tesla system will offer new perspectives in neuroscience and fundamental research. We anticipate that this initiative will usher in a new era of ultra-high-field MR.

An untrained deep learning method for reconstructing dynamic MR images from accelerated model-based data

PURPOSE

To implement physics-based regularization as a stopping condition in tuning an untrained deep neural network for reconstructing MR images from accelerated data.

METHODS

The ConvDecoder (CD) neural network was trained with a physics-based regularization term incorporating the spoiled gradient echo equation that describes variable-flip angle data. Fully-sampled variable-flip angle k-space data were retrospectively accelerated by factors of R = {8, 12, 18, 36} and reconstructed with CD, CD with the proposed regularization (CD + r), locally low-rank (LR) reconstruction, and compressed sensing with L1-wavelet regularization (L1). Final images from CD + r training were evaluated at the "argmin" of the regularization loss; whereas the CD, LR, and L1 reconstructions were chosen optimally based on ground truth data. The performance measures used were the normalized RMS error, the concordance correlation coefficient, and the structural similarity index.

RESULTS

The CD + r reconstructions, chosen using the stopping condition, yielded structural similarity indexs that were similar to the CD (p = 0.47) and LR structural similarity indexs (p = 0.95) across R and that were significantly higher than the L1 structural similarity indexs (p = 0.04). The concordance correlation coefficient values for the CD + r T1 maps across all R and subjects were greater than those corresponding to the L1 (p = 0.15) and LR (p = 0.13) T1 maps, respectively. For R ≥ 12 (≤4.2 min scan time), L1 and LR T1 maps exhibit a loss of spatially refined details compared to CD + r.

CONCLUSION

The use of an untrained neural network together with a physics-based regularization loss shows promise as a measure for determining the optimal stopping point in training without relying on fully-sampled ground truth data.

Deep Learning Initialized Compressed Sensing (Deli-CS) in Volumetric Spatio-Temporal Subspace Reconstruction

Introduction

Spatio-temporal MRI methods enable whole-brain multi-parametric mapping at ultra-fast acquisition times through efficient k-space encoding, but can have very long reconstruction times, which limit their integration into clinical practice. Deep learning (DL) is a promising approach to accelerate reconstruction, but can be computationally intensive to train and deploy due to the large dimensionality of spatio-temporal MRI. DL methods also need large training data sets and can produce results that don't match the acquired data if data consistency is not enforced. The aim of this project is to reduce reconstruction time using DL whilst simultaneously limiting the risk of deep learning induced hallucinations, all with modest hardware requirements.

Methods

Deep Learning Initialized Compressed Sensing (Deli-CS) is proposed to reduce the reconstruction time of iterative reconstructions by "kick-starting" the iterative reconstruction with a DL generated starting point. The proposed framework is applied to volumetric multi-axis spiral projection MRF that achieves whole-brain T1 and T2 mapping at 1-mm isotropic resolution for a 2-minute acquisition. First, the traditional reconstruction is optimized from over two hours to less than 40 minutes while using more than 90% less RAM and only 4.7 GB GPU memory, by using a memory-efficient GPU implementation. The Deli-CS framework is then implemented and evaluated against the above reconstruction.

Results

Deli-CS achieves comparable reconstruction quality with 50% fewer iterations bringing the full reconstruction time to 20 minutes.

Conclusion

Deli-CS reduces the reconstruction time of subspace reconstruction of volumetric spatio-temporal acquisitions by providing a warm start to the iterative reconstruction algorithm.

Improved MR temperature imaging at 0.5 T using view-sharing accelerated multiecho thermometry for MR-guided laser interstitial thermal therapy

The aim of the current study was to improve temperature-monitoring precision using multiecho proton resonance frequency shift-based thermometry with view-sharing acceleration for MR-guided laser interstitial thermal therapy (MRgLITT) on a 0.5-T low-field MR system. Both precision and speed of the temperature measurement for clinical MRgLITT treatments suffer at low field, due to reduced image signal-to-noise ratio (SNR), decreased temperature-induced phase changes, and limited RF receiver channels. In this work, a bipolar multiecho gradient-recalled echo sequence with a temperature-to-noise ratio optimal weighted echo combination is applied to improve the temperature precision. A view-sharing-based approach is utilized to accelerate signal acquisitions while preserving image SNRs. The method was evaluated using ex vivo (pork and pig brain) LITT heating experiments and in vivo (human brain) nonheating experiments on a high-performance 0.5-T scanner. In terms of results, (1) after echo combination, multiecho thermometry (i.e., ~7.5-40.5 ms, 7 TEs) provides ~1.5-1.9 times higher temperature precision than the no echo combination case (i.e., TE7 = 40.5 ms) within the same readout bandwidth. Additionally, echo registration is necessary for the bipolar multiecho sequence; (2) for a threefold acceleration, the view-sharing approach with variable-density subsampling shows around 1.8 times lower temperature errors than the GRAPPA method. Particularly for view-sharing, variable-density subsampling performs better than Interleave subsampling; and (3) ex vivo heating and in vivo nonheating experiments demonstrated that the temperature accuracy was less than 0.5 ° C $$ {}^{{}^{\circ}}\mathrm{C} $$ and that the temperature precision was less than 0.6 ° C $$ {}^{{}^{\circ}}\mathrm{C} $$ using the proposed 0.5-T thermometry. It was concluded that view-sharing accelerated multiecho thermometry is a practical temperature measurement approach for MRgLITT at 0.5 T.

High-Resolution 3D Magnetic Resonance Fingerprinting With a Graph Convolutional Network

Magnetic resonance fingerprinting (MRF) is a novel quantitative imaging framework for rapid and simultaneous quantification of multiple tissue properties. 3D MRF allows higher through-plane resolution, but the acquisition process is slow when whole-brain coverage is needed. Existing methods for acceleration mainly rely on GRAPPA for k-space interpolation in the partition-encoding direction, limiting the acceleration factor to 2 or 3. In this work, we replace GRAPPA with a deep learning approach for accurate tissue quantification with greater acceleration. Specifically, a graph convolution network (GCN) is developed to cater to the non-Cartesian spiral sampling trajectories typical in MRF acquisition. The GCN maintains high quantification accuracy with up to 6-fold acceleration and allows 1mm isotropic resolution whole-brain 3D MRF data to be acquired in 3min and submillimeter 3D MRF (0.8mm) in 5min, greatly improving the feasibility of MRF in clinical settings.

Improving robustness of 3D multi-shot EPI by structured low-rank reconstruction of segmented CAIPI sampling for fMRI at 7T

Three-dimensional (3D) encoding methods are increasingly being explored as alternatives to two-dimensional (2D) multi-slice acquisitions in fMRI, particularly in cases where high isotropic resolution is needed. 3D multi-shot EPI acquisition, as the workhorse of 3D fMRI imaging, is susceptible to physiological fluctuations which can induce inter-shot phase variations, and thus reducing the achievable tSNR, negating some of the benefit of 3D encoding. This issue can be particularly problematic at ultra-high fields like 7T, which have more severe off-resonance effects. In this work, we aim to improve the temporal stability of 3D multi-shot EPI at 7T by improving its robustness to inter-shot phase variations. We presented a 3D segmented CAIPI sampling trajectory ("seg-CAIPI") and an improved reconstruction method based on Hankel structured low-rank matrix recovery. Simulation and in-vivo results demonstrate that the combination of the seg-CAIPI sampling scheme and the proposed structured low-rank reconstruction is a promising way to effectively reduce the unwanted temporal variance induced by inter-shot physiological fluctuations, and thus improve the robustness of 3D multi-shot EPI for fMRI.

Accelerated Simultaneous T2 and T2* Mapping of Multiple Sclerosis Lesions Using Compressed Sensing Reconstruction of Radial RARE-EPI MRI

(1) Background: Radial RARE-EPI MRI facilitates simultaneous T₂ and T₂* mapping (2in1-RARE-EPI). With modest undersampling (R = 2), the speed gain of 2in1-RARE-EPI relative to Multi-Spin-Echo and Multi-Gradient-Recalled-Echo references is limited. Further reduction in scan time is crucial for clinical studies investigating T₂ and T₂* as imaging biomarkers. We demonstrate the feasibility of further acceleration, utilizing compressed sensing (CS) reconstruction of highly undersampled 2in1-RARE-EPI. (2) Methods: Two-fold radially-undersampled 2in1-RARE-EPI data from phantoms, healthy volunteers (n = 3), and multiple sclerosis patients (n = 4) were used as references, and undersampled (R_extra = 1-12, effective undersampling R_eff = 2-24). For each echo time, images were reconstructed using CS-reconstruction. For T₂ (RARE module) and T₂* mapping (EPI module), a linear least-square fit was applied to the images. T₂ and T₂* from CS-reconstruction of undersampled data were benchmarked against values from CS-reconstruction of the reference data. (3) Results: We demonstrate accelerated simultaneous T₂ and T₂* mapping using undersampled 2in1-RARE-EPI with CS-reconstruction is feasible. For R_extra = 6 (TA = 01:39 min), the overall MAPE was ≤8% (T₂*) and ≤4% (T₂); for R_extra = 12 (TA = 01:06 min), the overall MAPE was <13% (T₂*) and <5% (T₂). (4) Conclusion: Substantial reductions in scan time are achievable for simultaneous T₂ and T₂* mapping of the brain using highly undersampled 2in1-RARE-EPI with CS-reconstruction.

Rapid high-fidelity T 2 * mapping using single-shot overlapping-echo acquisition and deep learning reconstruction

PURPOSE

To develop and evaluate a single-shot quantitative MRI technique called GRE-MOLED (gradient-echo multiple overlapping-echo detachment) for rapid T 2 * $$ {T}_2^{\ast } $$ mapping.

METHODS

In GRE-MOLED, multiple echoes with different TEs are generated and captured in a single shot of the k-space through MOLED encoding and EPI readout. A deep neural network, trained by synthetic data, was employed for end-to-end parametric mapping from overlapping-echo signals. GRE-MOLED uses pure GRE acquisition with a single echo train to deliver T 2 * $$ {T}_2^{\ast } $$ maps less than 90 ms per slice. The self-registered B₀ information modulated in image phase was utilized for distortion-corrected parametric mapping. The proposed method was evaluated in phantoms, healthy volunteers, and task-based FMRI experiments.

RESULTS

The quantitative results of GRE-MOLED T 2 * $$ {T}_2^{\ast } $$ mapping demonstrated good agreement with those obtained from the multi-echo GRE method (Pearson's correlation coefficient = 0.991 and 0.973 for phantom and in vivo brains, respectively). High intrasubject repeatability (coefficient of variation <1.0%) were also achieved in scan-rescan test. Enabled by deep learning reconstruction, GRE-MOLED showed excellent robustness to geometric distortion, noise, and random subject motion. Compared to the conventional FMRI approach, GRE-MOLED also achieved a higher temporal SNR and BOLD sensitivity in task-based FMRI.

CONCLUSION

GRE-MOLED is a new real-time technique for T 2 * $$ {T}_2^{\ast } $$ quantification with high efficiency and quality, and it has the potential to be a better quantitative BOLD detection method.

Spatiotemporal encoding MRI using subspace-constrained sampling and locally-low-rank regularization: Applications to diffusion weighted and diffusion kurtosis imaging of human brain and prostate

The benefits of performing locally low-rank (LLR) reconstructions on subsampled diffusion weighted and diffusion kurtosis imaging data employing spatiotemporal encoding (SPEN) methods, is investigated. SPEN allows for self-referenced correction of motion-induced phase errors in case of interleaved diffusion-oriented acquisitions, and allows one to overcome distortions otherwise observed along EPI's phase-encoded dimension. In combination with LLR-based reconstructions of the pooled imaging data and with a joint subsampling of b-weighted and interleaved images, additional improvements in terms of sensitivity as well as shortened acquisition times are demonstrated, without noticeable penalties. Details on how the LLR-regularized, subspace-constrained image reconstructions were adapted to SPEN are given; the improvements introduced by adopting these reconstruction frameworks for the accelerated acquisition of diffusivity and of kurtosis imaging data in both relatively homogeneous regions like the human brain and in more challenging regions like the human prostate, are presented and discussed within the context of similar efforts in the field.

A joint linear reconstruction for multishot diffusion weighted non-Carr-Purcell-Meiboom-Gill fast spin echo with full signal

PURPOSE

Diffusion weighted Fast Spin Echo (DW-FSE) is a promising approach for distortionless DW imaging that is robust to system imperfections such as eddy currents and off-resonance. Due to non-Carr-Purcell-Meiboom-Gill (CPMG) magnetization, most DW-FSE sequences discard a large fraction of the signal (2 - 2 × $$ \sqrt{2}-2\times $$ ), reducing signal-to-noise ratio (SNR) efficiency compared to DW-EPI. The full FSE signal can be preserved by quadratically incrementing the transmit phase of the refocusing pulses, but this method of resolving non-CPMG magnetization has only been applied to single-shot DW-FSE due to challenges associated with image reconstruction. We present a joint linear reconstruction for multishot quadratic phase increment data that addresses these challenges and corrects ghosting from both shot-to-shot phase and intrashot signal oscillations. Multishot imaging reduces T2 blur and joint reconstruction of shots improves g-factor performance. A thorough analysis on the condition number of the proposed linear system is described.

METHODS

A joint multishot reconstruction is derived from the non-CPMG signal model. Multishot quadratic phase increment DW-FSE was tested in a standardized diffusion phantom and compared to single-shot DW-FSE and DW-EPI in vivo in the brain, cervical spine, and prostate. The pseudo multiple replica technique was applied to generate g-factor and SNR maps.

RESULTS

The proposed joint shot reconstruction eliminates ghosting from shot-to-shot phase and intrashot oscillations. g-factor performance is improved compared to previously proposed reconstructions, permitting efficient multishot imaging. apparent diffusion coefficient estimates in phantom experiments and in vivo are comparable to those obtained with conventional methods.

CONCLUSION

Multi-shot non-CPMG DW-FSE data with full signal can be jointly reconstructed using a linear model.

Single-shot multi-parametric mapping based on multiple overlapping-echo detachment (MOLED) imaging

Multi-parametric quantitative magnetic resonance imaging (mqMRI) allows the characterization of multiple tissue properties non-invasively and has shown great potential to enhance the sensitivity of MRI measurements. However, real-time mqMRI during dynamic physiological processes or general motions remains challenging. To overcome this bottleneck, we propose a novel mqMRI technique based on multiple overlapping-echo detachment (MOLED) imaging, termed MQMOLED, to enable mqMRI in a single shot. In the data acquisition of MQMOLED, multiple MR echo signals with different multi-parametric weightings and phase modulations are generated and acquired in the same k-space. The k-space data is Fourier transformed and fed into a well-trained neural network for the reconstruction of multi-parametric maps. We demonstrated the accuracy and repeatability of MQMOLED in simultaneous mapping apparent proton density (APD) and any two parameters among T₂, T₂*, and apparent diffusion coefficient (ADC) in 130-170 ms. The abundant information delivered by the multiple overlapping-echo signals in MQMOLED makes the technique potentially robust to system imperfections, such as inhomogeneity of static magnetic field or radiofrequency field. Benefitting from the single-shot feature, MQMOLED exhibits a strong motion tolerance to the continuous movements of subjects. For the first time, it captured the synchronous changes of ADC, T₂, and T₁-weighted APD in contrast-enhanced perfusion imaging on patients with brain tumors, providing additional information about vascular density to the hemodynamic parametric maps. We expect that MQMOLED would promote the development of mqMRI technology and greatly benefit the applications of mqMRI, including therapeutics and analysis of metabolic/functional processes.

Motion-corrected 3D-EPTI with efficient 4D navigator acquisition for fast and robust whole-brain quantitative imaging

PURPOSE

To develop a motion estimation and correction method for motion-robust three-dimensional (3D) quantitative imaging with 3D-echo-planar time-resolved imaging.

THEORY AND METHODS

The 3D-echo-planar time-resolved imaging technique was designed with additional four-dimensional navigator acquisition (x-y-z-echoes) to achieve fast and motion-robust quantitative imaging of the human brain. The four-dimensional-navigator is inserted into the relaxation-recovery deadtime of the sequence in every pulse TR (∼2 s) to avoid extra scan time, and to provide continuous tracking of the 3D head motion and B0 -inhomogeneity changes. By using an optimized spatiotemporal encoding combined with a partial-Fourier scheme, the navigator acquires a large central k-t data block for accurate motion estimation using only four small-flip-angle excitations and readouts, resulting in negligible signal-recovery reduction to the 3D-echo-planar time-resolved imaging acquisition. By incorporating the estimated motion and B0 -inhomogeneity changes into the reconstruction, multi-contrast images can be recovered with reduced motion artifacts.

RESULTS

Simulation shows the cost to the SNR efficiency from the added navigator acquisitions is <1%. Both simulation and in vivo retrospective experiments were conducted, that demonstrate the four-dimensional navigator provided accurate estimation of the 3D motion and B0 -inhomogeneity changes, allowing effective reduction of image artifacts in quantitative maps. Finally, in vivo prospective undersampling acquisition was performed with and without head motion, in which the motion corrupted data after correction show close image quality and consistent quantifications to the motion-free scan, providing reliable quantitative measurements even with head motion.

CONCLUSION

The proposed four-dimensional navigator acquisition provides reliable tracking of the head motion and B0 change with negligible SNR cost, equips the 3D-echo-planar time-resolved imaging technique for motion-robust and efficient quantitative imaging.

Mapping oxidative metabolism in the human brain with calibrated fMRI in health and disease

Conventional functional MRI (fMRI) with blood-oxygenation level dependent (BOLD) contrast is an important tool for mapping human brain activity non-invasively. Recent interest in quantitative fMRI has renewed the importance of oxidative neuroenergetics as reflected by cerebral metabolic rate of oxygen consumption (CMR_O2) to support brain function. Dynamic CMR_O2 mapping by calibrated fMRI require multi-modal measurements of BOLD signal along with cerebral blood flow (CBF) and/or volume (CBV). In human subjects this "calibration" is typically performed using a gas mixture containing small amounts of carbon dioxide and/or oxygen-enriched medical air, which are thought to produce changes in CBF (and CBV) and BOLD signal with minimal or no CMR_O2 changes. However non-human studies have demonstrated that the "calibration" can also be achieved without gases, revealing good agreement between CMR_O2 changes and underlying neuronal activity (e.g., multi-unit activity and local field potential). Given the simpler set-up of gas-free calibrated fMRI, there is evidence of recent clinical applications for this less intrusive direction. This up-to-date review emphasizes technological advances for such translational gas-free calibrated fMRI experiments, also covering historical progression of the calibrated fMRI field that is impacting neurological and neurodegenerative investigations of the human brain.

Rapid Geometry-Corrected Echo-Planar Diffusion Imaging at Ultrahigh Field: Fusing View Angle Tilting and Point-Spread Function Mapping

PURPOSE

Severe geometric distortions induced by tissue susceptibility, water-fat chemical shift, and eddy currents pose a substantial obstacle in single-shot EPI, especially for high-resolution imaging at ultrahigh field. View angle tilting (VAT)-EPI can mitigate in-plane distortion. However, the accompanied strong image blurring prevented its widespread applications. On the other hand, point-spread function mapping (PSF)-EPI can correct distortion and blurring accurately but requires prolonged scan time. We present fused VAT-PSF-EPI and possibilities for acceleration.

METHODS

MR signal equations were explicitly derived to quantify image blurring in VAT-EPI and the maximum acceleration capacity in VAT-PSF-EPI. To validate the theoretical prediction, phantom measurements with varying in-plane parallel imaging factors, slice thicknesses, and RF pulses were conducted at 7 Tesla. In addition, in vivo human brain scans were acquired with T₂ and diffusion weighting to assess distortion and blurring correction.

RESULTS

VAT can effectively suppress distortion, and the introduced image blurring is corrected through PSF encoding. Up to fourfold acceleration (only 5 shots) in VAT-PSF-EPI was achieved compared with standard PSF-EPI without VAT. VAT-induced signal loss was mitigated by adjusting the sequence parameters and EPI resolution. In vivo T₂ -weighted EPI data with 1.4 mm³ resolution demonstrate immunity to water-fat chemical shift-induced distortion. Very high-spatial resolution diffusion-weighted EPI (0.7 × 0.7 × 2.8 mm³ and 1.2 mm³ ) demonstrates the immunity to eddy current-induced distortion.

CONCLUSION

VAT-PSF-EPI is a novel spin-echo EPI-based sequence for fast high-resolution diffusion imaging at ultrahigh field.

A Geometric View of Signal Sensitivity Metrics in multi-echo fMRI

In multi-echo fMRI (ME-fMRI), two metrics have been widely used to measure the performance of various acquisition and analysis approaches. These are temporal SNR (tSNR) and differential contrast-to-noise ratio (dCNR). A key step in ME-fMRI is the weighted combination of the data from multiple echoes, and prior work has examined the dependence of tSNR and dCNR on the choice of weights. However, most studies have focused on only one of these two metrics, and the relationship between the two metrics has not been examined. In this work, we present a geometric view that offers greater insight into the relation between the two metrics and their weight dependence. We identify three major regimes: (1) a tSNR robust regime in which tSNR is robust to the weight selection with most weight variants achieving close to optimal performance, whereas dCNR shows a pronounced dependence on the weights with most variants achieving suboptimal performance; (2) a dCNR robust regime in which dCNR is robust to the weight selection with most weight variants achieving close to optimal performance, while tSNR exhibits a strong dependence on the weights with most variants achieving significantly lower than optimal performance; and (3) a within-type robust regime in which both tSNR and dCNR achieve nearly optimal performance when the form of the weights are variants of their respective optimal weights and exhibit a moderate decrease in performance for other weight variants. Insight into the behavior observed in the different regimes is gained by considering spherical representations of the weight dependence of the components used to form each metric. For multi-echo acquisitions, dCNR is shown to be more directly related than tSNR to measures of CNR and signal-to-noise ratio (SNR) for both task-based and resting-state fMRI scans, respectively.

Accelerating multi-echo chemical shift encoded water-fat MRI using model-guided deep learning

PURPOSE

To accelerate chemical shift encoded (CSE) water-fat imaging by applying a model-guided deep learning water-fat separation (MGDL-WF) framework to the undersampled k-space data.

METHODS

A model-guided deep learning water-fat separation framework is proposed for the acceleration using Cartesian/radial undersampling data. The proposed MGDL-WF combines the power of CSE water-fat imaging model and data-driven deep learning by jointly using a multi-peak fat model and a modified residual U-net network. The model is used to guide the image reconstruction, and the network is used to capture the artifacts induced by the undersampling. A data consistency layer is used in MGDL-WF to ensure the output images to be consistent with the k-space measurements. A Gauss-Newton iteration algorithm is adapted for the gradient updating of the networks.

RESULTS

Compared with the compressed sensing water-fat separation (CS-WF) algorithm/2-step procedure algorithm, the MGDL-WF increased peak signal-to-noise ratio (PSNR) by 5.31/5.23, 6.11/4.54, and 4.75 dB/1.88 dB with Cartesian sampling, and by 4.13/6.53, 2.90/4.68, and 1.68 dB/3.48 dB with radial sampling, at acceleration rates (R) of 4, 6, and 8, respectively. By using MGDL-WF, radial sampling increased the PSNR by 2.07 dB at R = 8, compared with Cartesian sampling.

CONCLUSIONS

The proposed MGDL-WF enables exploiting features of the water images and fat images from the undersampled multi-echo data, leading to improved performance in the accelerated CSE water-fat imaging. By using MGDL-WF, radial sampling can further improve the image quality with comparable scan time in comparison with Cartesian sampling.

Synthetic MR: Physical principles, clinical implementation, and new developments

Current clinical MR imaging practices rely on the qualitative assessment of images for diagnosis and treatment planning. While contrast in MR images is dependent on the spin parameters of the imaged tissue, pixel values on MR images are relative and are not scaled to represent any tissue properties. Synthetic MR is a fully featured imaging workflow consisting of efficient multiparameter mapping acquisition, synthetic image generation, and volume quantitation of brain tissues. As the application becomes more widely available on multiple vendors and scanner platforms, it has also gained widespread adoption as clinicians begin to recognize the benefits of rapid quantitation. This review will provide details about the sequence with a focus on the physical principles behind its relaxometry mechanisms. It will present an overview of the products in their current form and some potential issues when implementing it in the clinic. It will conclude by highlighting some recent advances of the technique, including a 3D mapping method and its associated applications. This article is protected by copyright. All rights reserved.

Blip up-down acquisition for spin- and gradient-echo imaging (BUDA-SAGE) with self-supervised denoising enables efficient T2 , T2 *, para- and dia-magnetic susceptibility mapping

PURPOSE

To rapidly obtain high resolution T₂ , T₂ *, and quantitative susceptibility mapping (QSM) source separation maps with whole-brain coverage and high geometric fidelity.

METHODS

We propose Blip Up-Down Acquisition for Spin And Gradient Echo imaging (BUDA-SAGE), an efficient EPI sequence for quantitative mapping. The acquisition includes multiple T₂ *-, T₂ '-, and T₂ -weighted contrasts. We alternate the phase-encoding polarities across the interleaved shots in this multi-shot navigator-free acquisition. A field map estimated from interim reconstructions was incorporated into the joint multi-shot EPI reconstruction with a structured low rank constraint to eliminate distortion. A self-supervised neural network (NN), MR-Self2Self (MR-S2S), was used to perform denoising to boost SNR. Using Slider encoding allowed us to reach 1 mm isotropic resolution by performing super-resolution reconstruction on volumes acquired with 2 mm slice thickness. Quantitative T₂ (=1/R₂ ) and T₂ * (=1/R₂ *) maps were obtained using Bloch dictionary matching on the reconstructed echoes. QSM was estimated using nonlinear dipole inversion on the gradient echoes. Starting from the estimated R₂ /R₂ * maps, R₂ ' information was derived and used in source separation QSM reconstruction, which provided additional para- and dia-magnetic susceptibility maps.

RESULTS

In vivo results demonstrate the ability of BUDA-SAGE to provide whole-brain, distortion-free, high-resolution, multi-contrast images and quantitative T₂ /T₂ * maps, as well as yielding para- and dia-magnetic susceptibility maps. Estimated quantitative maps showed comparable values to conventional mapping methods in phantom and in vivo measurements.

CONCLUSION

BUDA-SAGE acquisition with self-supervised denoising and Slider encoding enables rapid, distortion-free, whole-brain T₂ /T₂ * mapping at 1 mm isotropic resolution under 90 s.

3D Echo Planar Time-resolved Imaging (3D-EPTI) for ultrafast multi-parametric quantitative MRI

Multi-parametric quantitative MRI has shown great potential to improve the sensitivity and specificity of clinical diagnosis and to enhance our understanding of complex brain processes, but suffers from long scan time especially at high spatial resolution. To address this longstanding challenge, we introduce a novel approach, termed 3D Echo Planar Time-resolved Imaging (3D-EPTI), which significantly increases the acceleration capacity of MRI sampling, and provides high acquisition efficiency for multi-parametric MRI. This is achieved by exploiting the spatiotemporal correlation of MRI data at multiple timescales through new encoding strategies within and between efficient continuous readouts. Specifically, an optimized spatiotemporal CAIPI encoding within the readouts combined with a radial-block sampling strategy across the readouts enables an acceleration rate of 800 fold in the k-t space. A subspace reconstruction was employed to resolve thousands of high-quality multi-contrast images. We have demonstrated the ability of 3D-EPTI to provide robust and repeatable whole-brain simultaneous T1, T2, T2*, PD and B1+ mapping at high isotropic resolution within minutes (e.g., 1-mm isotropic resolution in 3 minutes), and to enable submillimeter multi-parametric imaging to study detailed brain structures.

MULTI-parametric MR imaging with fLEXible design (MULTIPLEX)

PURPOSE

To introduce a gradient echo (GRE) -based method, namely MULTIPLEX, for single-scan 3D multi-parametric MRI with high resolution, signal-to-noise ratio (SNR), accuracy, efficiency, and acquisition flexibility.

THEORY

With a comprehensive design with dual-repetition time (TR), dual flip angle (FA), multi-echo, and optional flow modulation features, the MULTIPLEX signals contain information on radiofrequency (RF) B_1t fields, proton density, T₁ , susceptibility and blood flows, facilitating multiple qualitative images and parametric maps.

METHODS

MULTIPLEX was evaluated on system phantom and human brains, via visual inspection for image contrasts and quality or quantitative evaluation via simulation, phantom scans and literature comparison. Region-of-interest (ROI) analysis was performed on parametric maps of the system phantom and brain scans, extracting the mean and SD of the T₁ , T2∗ , proton density (PD), and/or quantitative susceptibility mapping (QSM) values for comparison with reference values or literature.

RESULTS

One MULTIPLEX scan offers multiple sets of images, including but not limited to: composited PDW/T₁ W/ T2∗ W, aT₁ W, SWI, MRA (optional), B_1t map, T₁ map, T2∗ / R2∗ maps, PD map, and QSM. The quantitative error of phantom T₁ , T2∗ and PD mapping were <5%, and those in brain scans were in good agreement with literature. MULTIPLEX scan times for high resolution (0.68 × 0.68 × 2 mm³ ) whole brain coverage were about 7.5 min, while processing times were <1 min. With flow modulation, additional MRA images can be obtained without affecting the quality or accuracy of other images.

CONCLUSION

The proposed MUTLIPLEX method possesses great potential for multi-parametric MR imaging.

Functional connectome of brainstem nuclei involved in autonomic, limbic, pain and sensory processing in living humans from 7 Tesla resting state fMRI

Despite remarkable advances in mapping the functional connectivity of the cortex, the functional connectivity of subcortical regions is understudied in living humans. This is the case for brainstem nuclei that control vital processes, such as autonomic, limbic, nociceptive and sensory functions. This is because of the lack of precise brainstem nuclei localization, of adequate sensitivity and resolution in the deepest brain regions, as well as of optimized processing for the brainstem. To close the gap between the cortex and the brainstem, on 20 healthy subjects, we computed a correlation-based functional connectome of 15 brainstem nuclei involved in autonomic, limbic, nociceptive, and sensory function (superior and inferior colliculi, ventral tegmental area-parabrachial pigmented nucleus complex, microcellular tegmental nucleus-prabigeminal nucleus complex, lateral and medial parabrachial nuclei, vestibular and superior olivary complex, superior and inferior medullary reticular formation, viscerosensory motor nucleus, raphe magnus, pallidus, and obscurus, and parvicellular reticular nucleus - alpha part) with the rest of the brain. Specifically, we exploited 1.1mm isotropic resolution 7 Tesla resting-state fMRI, ad-hoc coregistration and physiological noise correction strategies, and a recently developed probabilistic template of brainstem nuclei. Further, we used 2.5mm isotropic resolution resting-state fMRI data acquired on a 3 Tesla scanner to assess the translatability of our results to conventional datasets. We report highly consistent correlation coefficients across subjects, confirming available literature on autonomic, limbic, nociceptive and sensory pathways, as well as high interconnectivity within the central autonomic network and the vestibular network. Interestingly, our results showed evidence of vestibulo-autonomic interactions in line with previous work. Comparison of 7 Tesla and 3 Tesla findings showed high translatability of results to conventional settings for brainstem-cortical connectivity and good yet weaker translatability for brainstem-brainstem connectivity. The brainstem functional connectome might bring new insight in the understanding of autonomic, limbic, nociceptive and sensory function in health and disease.

Functional connectome of arousal and motor brainstem nuclei in living humans by 7 Tesla resting-state fMRI

Brainstem nuclei play a pivotal role in many functions, such as arousal and motor control. Nevertheless, the connectivity of arousal and motor brainstem nuclei is understudied in living humans due to the limited sensitivity and spatial resolution of conventional imaging, and to the lack of atlases of these deep tiny regions of the brain. For a holistic comprehension of sleep, arousal and associated motor processes, we investigated in 20 healthy subjects the resting-state functional connectivity of 18 arousal and motor brainstem nuclei in living humans. To do so, we used high spatial-resolution 7 Tesla resting-state fMRI, as well as a recently developed in-vivo probabilistic atlas of these nuclei in stereotactic space. Further, we verified the translatability of our brainstem connectome approach to conventional (e.g. 3 Tesla) fMRI. Arousal brainstem nuclei displayed high interconnectivity, as well as connectivity to the thalamus, hypothalamus, basal forebrain and frontal cortex, in line with animal studies and as expected for arousal regions. Motor brainstem nuclei showed expected connectivity to the cerebellum, basal ganglia and motor cortex, as well as high interconnectivity. Comparison of 3 Tesla to 7 Tesla connectivity results indicated good translatability of our brainstem connectome approach to conventional fMRI, especially for cortical and subcortical (non-brainstem) targets and to a lesser extent for brainstem targets. The functional connectome of 18 arousal and motor brainstem nuclei with the rest of the brain might provide a better understanding of arousal, sleep and accompanying motor function in living humans in health and disease.

Distortion-Free Diffusion Imaging Using Self-Navigated Cartesian Echo-Planar Time Resolved Acquisition and Joint Magnitude and Phase Constrained Reconstruction

Echo-planar time resolved imaging (EPTI) is an effective approach for acquiring high-quality distortion-free images with a multi-shot EPI (ms-EPI) readout. As with traditional ms-EPI acquisitions, inter-shot phase variations present a main challenge when incorporating EPTI into a diffusion-prepared pulse sequence. The aim of this study is to develop a self-navigated Cartesian EPTI-based (scEPTI) acquisition together with a magnitude and phase constrained reconstruction for distortion-free diffusion imaging. A self-navigated Cartesian EPTI-based diffusion-prepared pulse sequence is designed. The different phase components in EPTI diffusion signal are analyzed and an approach to synthesize a fully phase-matched navigator for the inter-shot phase correction is demonstrated. Lastly, EPTI contains richer magnitude and phase information than conventional ms-EPI, such as the magnitude and phase correlations along the temporal dimension. The potential of these magnitude and phase correlations to enhance the reconstruction is explored. The reconstruction results with and without phase matching and with and without phase or magnitude constraints are compared. Compared with reconstruction without phase matching, the proposed phase matching method can improve the accuracy of inter-shot phase correction and reduce signal corruption in the final diffusion images. Magnitude constraints further improve image quality by suppressing the background noise and thereby increasing SNR, while phase constraints can mitigate possible image blurring from adding magnitude constraints. The high-quality distortion-free diffusion images and simultaneous diffusion-relaxometry imaging capacity provided by the proposed EPTI design represent a highly valuable tool for both clinical and neuroscientific assessments of tissue microstructure.

Simultaneous pure T2 and varying T2'-weighted BOLD fMRI using Echo Planar Time-resolved Imaging for mapping cortical-depth dependent responses

Spin-echo (SE) BOLD fMRI has high microvascular specificity, and thus provides a more reliable means to localize neural activity compared to conventional gradient-echo BOLD fMRI. However, the most common SE BOLD acquisition method, SE-EPI, is known to suffer from T2' contrast contamination with undesirable draining vein bias. To address this, in this study, we extended a recently developed distortion/blurring-free multi-shot EPI technique, Echo-Planar Time-resolved Imaging (EPTI), to cortical-depth dependent SE-fMRI at 7T to test whether it could provide purer SE BOLD contrast with minimal T2' contamination for improved neuronal specificity. From the same acquisition, the time-resolved feature of EPTI also provides a series of asymmetric SE (ASE) images with varying T2' weightings, and enables extraction of data equivalent to conventional SE EPI with different echo train lengths (ETLs). This allows us to systematically examine how T2'-contribution affects different SE acquisition strategies using a single dataset. A low-rank spatiotemporal subspace reconstruction was implemented for the SE-EPTI acquisition, which incorporates corrections for both shot-to-shot phase variations and dynamic B0 drifts. SE-EPTI was used in a visual task fMRI experiment to demonstrate that i) the pure SE image provided by EPTI results in the highest microvascular specificity; ii) the ASE EPTI series, with a graded introduction of T2' weightings at time points farther away from the pure SE, show a gradual sensitivity increase along with increasing draining vein bias; iii) the longer ETL seen in conventional SE EPI acquisitions will induce more draining vein bias. Consistent results were observed across multiple subjects, demonstrating the robustness of the proposed technique for SE-BOLD fMRI with high specificity.