Controlled aliasing in parallel imaging results in higher acceleration (CAIPIRINHA) for multi-slice imaging

Multiphoton simultaneous multislice imaging

PURPOSE

To develop multiphoton excitation techniques for simultaneous multislice (SMS) imaging and evaluate their performance and specific absorption rate (SAR) benefit. To improve multiphoton SMS reconstruction quality with a novel CAIPIRINHA (controlled aliasing in parallel imaging results in higher acceleration) design.

THEORY AND METHODS

When a conventional single-slice RF field is applied together with an oscillating gradient field, the two can combine to generate multiphoton excitation at multiple discrete spatial locations. Because the conventional RF is reused at multiple spatial locations, multiphoton excitation offers reduced SAR for SMS applications. CAIPIRINHA shifts are often used to improve parallel-imaging acceleration. Interestingly, CAIPIRINHA-type shifts can be obtained for multiphoton SMS by updating the oscillating gradient phase at every phase encode. In this work, both a gradient-echo and a spin-echo sequence with multiphoton CAIPIRINHA-SMS excitation pulses are implemented for in vivo human imaging at 3 T.

RESULTS

For three slices, multiphoton SMS provides a 51% reduction in SAR compared with conventional superposition SMS, whereas for five slices, SAR is reduced by 66%. Multiphoton SMS outperforms PINS (power independent of number of slices) and MultiPINS in terms of SAR reduction especially when the pulse duration is short, slices are thin, and/or the slice spacing is large. A custom CAIPIRINHA phase-encoding design for multiphoton SMS significantly improves reconstruction quality.

CONCLUSION

Multiphoton SMS excitation can be obtained by combining conventional single-slice RF pulses with an oscillating gradient and offers significant SAR benefits compared with conventional superposition SMS. A novel CAIPIRINHA design allows higher multiband factors for multiphoton SMS imaging.

Deep learning-based rapid image reconstruction and motion correction for high-resolution cartesian first-pass myocardial perfusion imaging at 3T

PURPOSE

To develop and evaluate a deep learning (DL) -based rapid image reconstruction and motion correction technique for high-resolution Cartesian first-pass myocardial perfusion imaging at 3T with whole-heart coverage for both single-slice (SS) and simultaneous multi-slice (SMS) acquisitions.

METHODS

3D physics-driven unrolled network architectures were utilized for the reconstruction of high-resolution Cartesian perfusion imaging. The SS and SMS multiband (MB) = 2 networks were trained from 135 slices from 20 subjects. Structural similarity index (SSIM), peak SNR (PSNR), and normalized RMS error (NRMSE) were assessed, and prospective images were blindly graded by two experienced cardiologists (5, excellent; 1, poor). For respiratory motion correction, a 2D U-Net based motion corrected network was proposed, and the temporal fidelity and second-order derivative were calculated to assess the performance of the motion correction.

RESULTS

Excellent performance was demonstrated in the proposed technique with high SSIM and PSNR, and low NRMSE. Image quality scores were (4.3 [4.3, 4.4], 4.5 [4.4, 4.6], 4.3 [4.3, 4.4], and 4.5 [4.3, 4.5]) for SS DL and SS L1-SENSE, MB = 2 DL and MB = 2 SMS-L1-SENSE, respectively, showing no statistically significant difference (p > 0.05 for SS and SMS) between (SMS)-L1-SENSE and the proposed DL technique. The network inference time was around 4 s per dynamic perfusion series with 40 frames while the time of (SMS)-L1-SENSE with GPU acceleration was approximately 30 min.

CONCLUSION

The proposed DL-based image reconstruction and motion correction technique enabled rapid and high-quality reconstruction for SS and SMS MB = 2 high-resolution Cartesian first-pass perfusion imaging at 3T.

Improved abdominal T1 weighted imaging at 0.55T

PURPOSE

Breath-held fat-suppressed volumetric T1-weighted MRI is an important and widely-used technique for evaluating the abdomen. Both fat-saturation and Dixon-based fat-suppression methods are used at conventional field strengths; however, both have challenges at lower field strengths (<1.5T) due to insufficient fat suppression and/or inadequate resolution. Specifically, at lower field strengths, fat saturation often fails due to the short T1 of lipid; and Cartesian Dixon imaging provides poor spatial resolution due to the need for a long ∆TE, due to the smaller ∆f between water and lipid. The purpose of this work is to demonstrate a new approach capable of simultaneously achieving excellent fat suppression and high spatial resolution on a 0.55T whole-body system.

METHODS

We applied 3D stack-of-spirals Dixon imaging at 0.55T, with compensation of concomitant field phase during reconstruction. The spiral readouts make efficient use of the requisite ∆TE. We compared this with 3D Cartesian Dixon imaging. Experiments were performed in 2 healthy and 10 elevated liver fat volunteers.

RESULTS

Stack-of-spirals Dixon imaging at 0.55T makes excellent use of the required ∆TE, provided high SNR efficiency and finer spatial resolution (1.7 × 1.7 × 5 mm3) compared Cartesian Dixon (3.5 × 3.5 × 5 mm3), within a 17-s breath-hold. We observed successful fat suppression, and improved definition of structures such as the liver, kidneys, and bowel.

CONCLUSION

We demonstrate that high-resolution single breath-hold volumetric abdominal T1-weighted imaging is feasible at 0.55T using spiral sampling and concomitant field correction. This is an attractive alternative to existing Cartesian-based methods, as it simultaneously provides high-resolution and excellent fat-suppression.

Fast reconstruction of SMS bSSFP myocardial perfusion images using noise map estimation network (NoiseMapNet): a head-to-head comparison with parallel imaging and iterative reconstruction

Background

Simultaneous multi-slice (SMS) bSSFP imaging enables stress myocardial perfusion imaging with high spatial resolution and increased spatial coverage. Standard parallel imaging techniques (e.g., TGRAPPA) can be used for image reconstruction but result in high noise level. Alternatively, iterative reconstruction techniques based on temporal regularization (ITER) improve image quality but are associated with reduced temporal signal fidelity and long computation time limiting their online use. The aim is to develop an image reconstruction technique for SMS-bSSFP myocardial perfusion imaging combining parallel imaging and image-based denoising using a novel noise map estimation network (NoiseMapNet), which preserves both sharpness and temporal signal profiles and that has low computational cost.

Methods

The proposed reconstruction of SMS images consists of a standard temporal parallel imaging reconstruction (TGRAPPA) with motion correction (MOCO) followed by image denoising using NoiseMapNet. NoiseMapNet is a deep learning network based on a 2D Unet architecture and aims to predict a noise map from an input noisy image, which is then subtracted from the noisy image to generate the denoised image. This approach was evaluated in 17 patients who underwent stress perfusion imaging using a SMS-bSSFP sequence. Images were reconstructed with (a) TGRAPPA with MOCO (thereafter referred to as TGRAPPA), (b) iterative reconstruction with integrated motion compensation (ITER), and (c) proposed NoiseMapNet-based reconstruction. Normalized mean squared error (NMSE) with respect to TGRAPPA, myocardial sharpness, image quality, perceived SNR (pSNR), and number of diagnostic segments were evaluated.

Results

NMSE of NoiseMapNet was lower than using ITER for both myocardium (0.045 ± 0.021 vs. 0.172 ± 0.041, p < 0.001) and left ventricular blood pool (0.025 ± 0.014 vs. 0.069 ± 0.020, p < 0.001). There were no significant differences between all methods for myocardial sharpness (p = 0.77) and number of diagnostic segments (p = 0.36). ITER led to higher image quality than NoiseMapNet/TGRAPPA (2.7 ± 0.4 vs. 1.8 ± 0.4/1.3 ± 0.6, p < 0.001) and higher pSNR than NoiseMapNet/TGRAPPA (3.0 ± 0.0 vs. 2.0 ± 0.0/1.3 ± 0.6, p < 0.001). Importantly, NoiseMapNet yielded higher pSNR (p < 0.001) and image quality (p < 0.008) than TGRAPPA. Computation time of NoiseMapNet was only 20s for one entire dataset.

Conclusion

NoiseMapNet-based reconstruction enables fast SMS image reconstruction for stress myocardial perfusion imaging while preserving sharpness and temporal signal profiles.

Deep Learning-based Image Enhancement Techniques for Fast MRI in Neuroimaging

Despite its superior soft tissue contrast and non-invasive nature, MRI requires long scan times due to its intrinsic signal acquisition principles, a main drawback which technological advancements in MRI have been focused on. In particular, scan time reduction is a natural requirement in neuroimaging due to detailed structures requiring high resolution imaging and often volumetric (3D) acquisitions, and numerous studies have recently attempted to harness deep learning (DL) technology in enabling scan time reduction and image quality improvement. Various DL-based image reconstruction products allow for additional scan time reduction on top of existing accelerated acquisition methods without compromising the image quality.

Non-Cartesian 3D-SPARKLING vs Cartesian 3D-EPI encoding schemes for functional Magnetic Resonance Imaging at 7 Tesla

The quest for higher spatial and/or temporal resolution in functional MRI (fMRI) while preserving a sufficient temporal signal-to-noise ratio (tSNR) has generated a tremendous amount of methodological contributions in the last decade ranging from Cartesian vs. non-Cartesian readouts, 2D vs. 3D acquisition strategies, parallel imaging and/or compressed sensing (CS) accelerations and simultaneous multi-slice acquisitions to cite a few. In this paper, we investigate the use of a finely tuned version of 3D-SPARKLING. This is a non-Cartesian CS-based acquisition technique for high spatial resolution whole-brain fMRI. We compare it to state-of-the-art Cartesian 3D-EPI during both a retinotopic mapping paradigm and resting-state acquisitions at 1mm3 (isotropic spatial resolution). This study involves six healthy volunteers and both acquisition sequences were run on each individual in a randomly-balanced order across subjects. The performances of both acquisition techniques are compared to each other in regards to tSNR, sensitivity to the BOLD effect and spatial specificity. Our findings reveal that 3D-SPARKLING has a higher tSNR than 3D-EPI, an improved sensitivity to detect the BOLD contrast in the gray matter, and an improved spatial specificity. Compared to 3D-EPI, 3D-SPARKLING yields, on average, 7% more activated voxels in the gray matter relative to the total number of activated voxels.

Quantification of whole-organ individual and bilateral renal metabolic rate of oxygen

PURPOSE

Renal metabolic rate of oxygen (rMRO2 ) is a potentially important biomarker of kidney function. The key parameters for rMRO2 quantification include blood flow rate (BFR) and venous oxygen saturation (SvO2 ) in a draining vessel. Previous approaches to quantify renal metabolism have focused on the single organ. Here, both kidneys are considered as one unit to quantify bilateral rMRO2 . A pulse sequence to facilitate bilateral rMRO2 quantification is introduced.

METHODS

To quantify bilateral rMRO2 , measurements of BFR and SvO2 are made along the inferior vena cava (IVC) at suprarenal and infrarenal locations. From the continuity equation, these four parameters can be related to derive an expression for bilateral rMRO2 . The recently reported K-MOTIVE pulse sequence was implemented at four locations: left kidney, right kidney, suprarenal IVC, and infrarenal IVC. A dual-band variant of K-MOTIVE (db-K-MOTIVE) was developed by incorporating simultaneous-multi-slice imaging principles. The sequence simultaneously measures BFR and SvO2 at suprarenal and infrarenal locations in a single pass of 21 s, yielding bilateral rMRO2 .

RESULTS

SvO2 and BFR are higher in suprarenal versus infrarenal IVC, and the renal veins are highly oxygenated (SvO2  >90%). Bilateral rMRO2 quantified in 10 healthy subjects (8 M, 30 ± 8 y) was found to be 291 ± 247 and 349 ± 300 (μmolO2 /min)/100 g, derived from K-MOTIVE and db-K-MOTIVE, respectively. In comparison, total rMRO2 from combining left and right was 329 ± 273 (μmolO2 /min)/100 g.

CONCLUSION

The present work demonstrates that bilateral rMRO2 quantification is feasible with fair reproducibility and physiological plausibility. The indirect method is a promising approach to compute bilateral rMRO2 when individual rMRO2 quantification is difficult.

The effect of and correction for through-slice dephasing on 2D gradient-echo double angle B 1 + mapping

PURPOSE

To show thatB 0 $$ {\mathrm{B}}_0 $$ variations through slice and slice profile effects are two major confounders affecting 2D dual angleB 1 + $$ {\mathrm{B}}_1^{+} $$ maps using gradient-echo signals and thus need to be corrected to obtain accurateB 1 + $$ {\mathrm{B}}_1^{+} $$ maps.

METHODS

The 2D gradient-echo transverse complex signal was Bloch-simulated and integrated across the slice dimension including nonlinear variations inB 0 $$ {\mathrm{B}}_0 $$ inhomogeneities through slice. A nonlinear least squares fit was used to find theB 1 + $$ {\mathrm{B}}_1^{+} $$ factor corresponding to the best match between the two gradient-echo signals experimental ratio and the Bloch-simulated ratio. The correction was validated in phantom and in vivo at 3T.

RESULTS

For our RF excitation pulse, the error in theB 1 + $$ {\mathrm{B}}_1^{+} $$ factor scales by approximately 3.8% for every 10 Hz/cm variation inB 0 $$ {\mathrm{B}}_0 $$ along the slice direction. Higher accuracy phantomB 1 + $$ {\mathrm{B}}_1^{+} $$ maps were obtained after applying the proposed correction; the root mean squareB 1 + $$ {\mathrm{B}}_1^{+} $$ error relative to the gold standardB 1 + $$ {\mathrm{B}}_1^{+} $$ decreased from 6.4% to 2.6%. In vivo whole-liverT 1 $$ {\mathrm{T}}_1 $$ maps using the correctedB 1 + $$ {\mathrm{B}}_1^{+} $$ map registered a significant decrease inT 1 $$ {\mathrm{T}}_1 $$ gradient through slice.

CONCLUSION

B 0 $$ {\mathrm{B}}_0 $$ inhomogeneities varying through slice were seen to have an impact on the accuracy of 2D double angleB 1 + $$ {\mathrm{B}}_1^{+} $$ maps using gradient-echo sequences. Consideration of this confounder is crucial for research relying on accurate knowledge of the true excitation flip angles, as is the case ofT 1 $$ {\mathrm{T}}_1 $$ mapping using a spoiled gradient recalled echo sequence.

Accelerated 2D Cartesian MRI with an 8-channel local B0 coil array combined with parallel imaging

PURPOSE

In MRI, the magnetization of nuclear spins is spatially encoded with linear gradients and radiofrequency receivers sensitivity profiles to produce images, which inherently leads to a long scan time. Cartesian MRI, as widely adopted for clinical scans, can be accelerated with parallel imaging and rapid magnetic field modulation during signal readout. Here, by using an 8-channel localB 0 $$ {\mathrm{B}}_0 $$ coil array, the modulation scheme optimized for sampling efficiency is investigated to speed up 2D Cartesian scans.

THEORY AND METHODS

An 8-channel localB 0 $$ {\mathrm{B}}_0 $$ coil array is made to carry sinusoidal currents during signal readout to accelerate 2D Cartesian scans. An MRI sampling theory based on reproducing kernel Hilbert space is exploited to visualize the efficiency of nonlinear encoding in arbitrary sampling duration. A field calibration method using current monitors for localB 0 $$ {\mathrm{B}}_0 $$ coils and the ESPIRiT algorithm is proposed to facilitate image reconstruction. Image acceleration with various modulation field shapes, aliasing control, and distinct modulation frequencies are scrutinized to find an optimized modulation scheme. A safety evaluation is conducted. In vivo 2D Cartesian scans are accelerated by the localB 0 $$ {\mathrm{B}}_0 $$ coils.

RESULTS

For 2D Cartesian MRI, the optimal modulation field by this localB 0 $$ {\mathrm{B}}_0 $$ array converges to a nearly linear gradient field. With the field calibration technique, it accelerates the in vivo scans (i.e., proved safe) by threefold and eightfold free of visible artifacts, without and with SENSE, respectively.

CONCLUSION

The nonlinear encoding analysis tool, the field calibration method, the safety evaluation procedures, and the in vivo reconstructed scans make significant steps to push MRI speed further with the localB 0 $$ {\mathrm{B}}_0 $$ coil array.

Magnetic resonance coherence pathway unraveling

Efficiently acquiring multi-contrast magnetic resonance imaging data is crucial for patient comfort and clinical throughput. Developing scan acceleration methods tailored for specific applications drastically improves the value of an MRI examination. Here, we propose a novel method to control the aliasing of simultaneously acquired images of multiple spin echo coherence pathways with the goal of producing high quality multi-contrast images from a single acquisition. Modulating the radiofrequency phase of several pulses applied in brief succession also uniquely modulates the phase of spin echo coherence pathways. A method, termed magnetic resonance coherence pathway unraveling (MR-CPU), to control the aliasing of simultaneously acquired coherence pathway images is developed here along with parallel imaging-based reconstruction methods to separate them. MR-CPU was validated in phantom experiments and tested in vivo. High levels of correlation between reference pathway images and MR-CPU-derived coherence pathway images were found from the phantom experiments. Minimal artifacts arising from the separation of the overlapped coherence pathway images were observed in vivo. MR-CPU provides a novel mechanism through which to acquire and separate multiple overlapped coherence pathway images, thus adding to the diagnostic potential of an MRI exam without the penalty of additional scan time.

Real-time automatic image-based slice tracking of gadolinium-filled balloon wedge catheter during MR-guided cardiac catheterization: A proof-of-concept study

PURPOSE

MR-guided cardiac catheterization procedures currently use passive tracking approaches to follow a gadolinium-filled catheter balloon during catheter navigation. This requires frequent manual tracking and repositioning of the imaging slice during navigation. In this study, a novel framework for automatic real-time catheter tracking during MR-guided cardiac catheterization is presented.

METHODS

The proposed framework includes two imaging modes (Calibration and Runtime). The sequence starts in Calibration mode, in which the 3D catheter coordinates are determined using a stack of 10-20 contiguous saturated slices combined with real-time image processing. The sequence then automatically switches to Runtime mode, where three contiguous slices (acquired with partial saturation), initially centered on the catheter balloon using the Calibration feedback, are acquired continuously. The 3D catheter balloon coordinates are estimated in real time from each Runtime slice stack using image processing. Each Runtime stack is repositioned to maintain the catheter balloon in the central slice based on the prior Runtime feedback. The sequence switches back to Calibration mode if the catheter is not detected. This framework was evaluated in a heart phantom and 3 patients undergoing MR-guided cardiac catheterization. Catheter detection accuracy and rate of catheter visibility were evaluated.

RESULTS

The automatic detection accuracy for the catheter balloon during the Calibration/Runtime mode was 100%/95% in phantom and 100%/97 ± 3% in patients. During Runtime, the catheter was visible in 82% and 98 ± 2% of the real-time measurements in the phantom and patients, respectively.

CONCLUSION

The proposed framework enabled real-time continuous automatic tracking of a gadolinium-filled catheter balloon during MR-guided cardiac catheterization.

Accelerated Musculoskeletal Magnetic Resonance Imaging

With a substantial growth in the use of musculoskeletal MRI, there has been a growing need to improve MRI workflow, and faster imaging has been suggested as one of the solutions for a more efficient examination process. Consequently, there have been considerable advances in accelerated MRI scanning methods. This article aims to review the basic principles and applications of accelerated musculoskeletal MRI techniques including widely used conventional acceleration methods, more advanced deep learning-based techniques, and new approaches to reduce scan time. Specifically, conventional accelerated MRI techniques, including parallel imaging, compressed sensing, and simultaneous multislice imaging, and deep learning-based accelerated MRI techniques, including undersampled MR image reconstruction, super-resolution imaging, artifact correction, and generation of unacquired contrast images, are discussed. Finally, new approaches to reduce scan time, including synthetic MRI, novel sequences, and new coil setups and designs, are also reviewed. We believe that a deep understanding of these fast MRI techniques and proper use of combined acceleration methods will synergistically improve scan time and MRI workflow in daily practice. EVIDENCE LEVEL: 3 TECHNICAL EFFICACY: Stage 1.

Cardiac MR Fingerprinting: Overview, Technical Developments, and Applications

Cardiovascular magnetic resonance (CMR) is an established imaging modality with proven utility in assessing cardiovascular diseases. The ability of CMR to characterize myocardial tissue using T1 - and T2 -weighted imaging, parametric mapping, and late gadolinium enhancement has allowed for the non-invasive identification of specific pathologies not previously possible with modalities like echocardiography. However, CMR examinations are lengthy and technically complex, requiring multiple pulse sequences and different anatomical planes to comprehensively assess myocardial structure, function, and tissue composition. To increase the overall impact of this modality, there is a need to simplify and shorten CMR exams to improve access and efficiency, while also providing reproducible quantitative measurements. Multiparametric MRI techniques that measure multiple tissue properties offer one potential solution to this problem. This review provides an in-depth look at one such multiparametric approach, cardiac magnetic resonance fingerprinting (MRF). The article is structured as follows. First, a brief review of single-parametric and (non-Fingerprinting) multiparametric CMR mapping techniques is presented. Second, a general overview of cardiac MRF is provided covering pulse sequence implementation, dictionary generation, fast k-space sampling methods, and pattern recognition. Third, recent technical advances in cardiac MRF are covered spanning a variety of topics, including simultaneous multislice and 3D sampling, motion correction algorithms, cine MRF, synthetic multicontrast imaging, extensions to measure additional clinically important tissue properties (proton density fat fraction, T2 *, and T1ρ ), and deep learning methods for image reconstruction and parameter estimation. The last section will discuss potential clinical applications, concluding with a perspective on how multiparametric techniques like MRF may enable streamlined CMR protocols. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 1.

Three-dimensional echo-shifted EPI with simultaneous blip-up and blip-down acquisitions for correcting geometric distortion

PURPOSE

EPI with blip-up/down acquisition (BUDA) can provide high-quality images with minimal distortions by using two readout trains with opposing phase-encoding gradients. Because of the need for two separate acquisitions, BUDA doubles the scan time and degrades the temporal resolution when compared to single-shot EPI, presenting a major challenge for many applications, particularly fMRI. This study aims at overcoming this challenge by developing an echo-shifted EPI BUDA (esEPI-BUDA) technique to acquire both blip-up and blip-down datasets in a single shot.

METHODS

A 3D esEPI-BUDA pulse sequence was designed by using an echo-shifting strategy to produce two EPI readout trains. These readout trains produced a pair of k-space datasets whose k-space trajectories were interleaved with opposite phase-encoding gradient directions. The two k-space datasets were separately reconstructed using a 3D SENSE algorithm, from which time-resolved B0 -field maps were derived using TOPUP in FSL and then input into a forward model of joint parallel imaging reconstruction to correct for geometric distortion. In addition, Hankel structured low-rank constraint was incorporated into the reconstruction framework to improve image quality by mitigating the phase errors between the two interleaved k-space datasets.

RESULTS

The 3D esEPI-BUDA technique was demonstrated in a phantom and an fMRI study on healthy human subjects. Geometric distortions were effectively corrected in both phantom and human brain images. In the fMRI study, the visual activation volumes and their BOLD responses were comparable to those from conventional 3D echo-planar images.

CONCLUSION

The improved imaging efficiency and dynamic distortion correction capability afforded by 3D esEPI-BUDA are expected to benefit many EPI applications.

3D CEST MRI with an unevenly segmented RF irradiation scheme: A feasibility study in brain tumor imaging

PURPOSE

To integrate 3D CEST EPI with an unevenly segmented RF irradiation module and preliminarily demonstrate it in the clinical setting.

METHODS

A CEST MRI with unevenly segmented RF saturation was implemented, including a long primary RF saturation to induce the steady-state CEST effect, maintained with repetitive short secondary RF irradiation between readouts. This configuration reduces relaxation-induced blur artifacts during acquisition, allowing fast 3D spatial coverage. Numerical simulations were performed to select parameters such as flip angle (FA), short RF saturation duration (Ts2), and the number of readout segments. The sequence was validated experimentally with data from a phantom, healthy volunteers, and a brain tumor patient.

RESULTS

Based on the numerical simulation and l-carnosine gel phantom experiment, FA, Ts2, and the number of segments were set to 20°, 0.3 s, and the range from 4 to 8, respectively. The proposed method minimized signal modulation in the human brain images in the kz direction during the acquisition and provided the blur artifacts-free CEST contrast over the whole volume. Additionally, the CEST contrast in the tumor tissue region is higher than in the contralateral normal tissue region.

CONCLUSIONS

It is feasible to implement a highly accelerated 3D EPI CEST imaging with unevenly segmented RF irradiation.

Super-resolution QSM in little or no additional time for imaging (NATIve) using 2D EPI imaging in 3 orthogonal planes

Quantitative Susceptibility Mapping has the potential to provide additional insights into neurological diseases but is typically based on a quite long (5-10 min) 3D gradient-echo scan which is highly sensitive to motion. We propose an ultra-fast acquisition based on three orthogonal (sagittal, coronal and axial) 2D simultaneous multi-slice EPI scans with 1 mm in-plane resolution and 3 mm thick slices. Images in each orientation are corrected for susceptibility-related distortions and co-registered with an iterative non-linear Minimum Deformation Averaging (Volgenmodel) approach to generate a high SNR, super-resolution data set with an isotropic resolution of close to 1 mm. The net acquisition time is 3 times the volume acquisition time of EPI or about 12 s, but the three volumes could also replace "dummy scans" in fMRI, making it feasible to acquire QSM in little or No Additional Time for Imaging (NATIve). NATIve QSM values agreed well with reference 3D GRE QSM in the basal ganglia in healthy subjects. In patients with multiple sclerosis, there was also a good agreement between the susceptibility values within lesions and control ROIs and all lesions which could be seen on 3D GRE QSMs could also be visualized on NATIve QSMs. The approach is faster than conventional 3D GRE by a factor of 25-50 and faster than 3D EPI by a factor of 3-5. As a 2D technique, NATIve QSM was shown to be much more robust to motion than the 3D GRE and 3D EPI, opening up the possibility of studying neurological diseases involving iron accumulation and demyelination in patients who find it difficult to lie still for long enough to acquire QSM data with conventional methods.

A 128-channel receive array for cortical brain imaging at 7 T

PURPOSE

A 128-channel receive-only array for brain imaging at 7 T was simulated, designed, constructed, and tested within a high-performance head gradient designed for high-resolution functional imaging.

METHODS

The coil used a tight-fitting helmet geometry populated with 128 loop elements and preamplifiers to fit into a 39 cm diameter space inside a built-in gradient. The signal-to-noise ratio (SNR) and parallel imaging performance (1/g) were measured in vivo and simulated using electromagnetic modeling. The histogram of 1/g factors was analyzed to assess the range of performance. The array's performance was compared to the industry-standard 32-channel receive array and a 64-channel research array.

RESULTS

It was possible to construct the 128-channel array with body noise-dominated loops producing an average noise correlation of 5.4%. Measurements showed increased sensitivity compared with the 32-channel and 64-channel array through a combination of higher intrinsic SNR and g-factor improvements. For unaccelerated imaging, the 128-channel array showed SNR gains of 17.6% and 9.3% compared to the 32-channel and 64-channel array, respectively, at the center of the brain and 42% and 18% higher SNR in the peripheral brain regions including the cortex. For R = 5 accelerated imaging, these gains were 44.2% and 24.3% at the brain center and 86.7% and 48.7% in the cortex. The 1/g-factor histograms show both an improved mean and a tighter distribution by increasing the channel count, with both effects becoming more pronounced at higher accelerations.

CONCLUSION

The experimental results confirm that increasing the channel count to 128 channels is beneficial for 7T brain imaging, both for increasing SNR in peripheral brain regions and for accelerated imaging.

Volumetric measurements of weak current-induced magnetic fields in the human brain at high resolution

PURPOSE

Clinical use of transcranial electrical stimulation (TES) requires accurate knowledge of the injected current distribution in the brain. MR current density imaging (MRCDI) uses measurements of the TES-induced magnetic fields to provide this information. However, sufficient sensitivity and image quality in humans in vivo has only been documented for single-slice imaging.

METHODS

A recently developed, optimally spoiled, acquisition-weighted, gradient echo-based 2D-MRCDI method has now been advanced for volume coverage with densely or sparsely distributed slices: The 3D rectilinear sampling (3D-DENSE) and simultaneous multislice acquisition (SMS-SPARSE) were optimized and verified by cable-loop experiments and tested with 1-mA TES experiments for two common electrode montages.

RESULTS

Comparisons between the volumetric methods against the 2D-MRCDI showed that relatively long acquisition times of 3D-DENSE using a single slab with six slices hindered the expected sensitivity improvement in the current-induced field measurements but improved sensitivity by 61% in the Laplacian of the field, on which some MRCDI reconstruction methods rely. Also, SMS-SPARSE acquisition of three slices, with a factor 2 CAIPIRINHA (controlled aliasing in parallel imaging results in higher acceleration) acceleration, performed best against the 2D-MRCDI with sensitivity improvements for the∆ B z , c $$ \Delta {B}_{z,c} $$ and Laplacian noise floors of 56% and 78% (baseline without current flow) as well as 43% and 55% (current injection into head). SMS-SPARSE reached a sensitivity of 67 pT for three distant slices at 2 × 2 × 3 mm3 resolution in 10 min of total scan time, and consistently improved image quality.

CONCLUSION

Volumetric MRCDI measurements with high sensitivity and image quality are well suited to characterize the TES field distribution in the human brain.

Global effects of aging on the hemodynamic response function in the human brain

In functional magnetic resonance imaging, the hemodynamic response function (HRF) is a transient, stereotypical response to local changes in cerebral hemodynamics and oxygen metabolism due to briefly (< 4 s) evoked neural activity. Accordingly, the HRF is often used as an impulse response with the assumption of linearity in data analysis. In cognitive aging studies, it has been very common to interpret differences in brain activation as age-related changes in neural activity. Contrary to this assumption, however, evidence has accrued that normal aging may also significantly affect the vasculature, thereby affecting cerebral hemodynamics and metabolism, confounding interpretation of fMRI aging studies. In this study, use was made of a multisensory stimulus to evoke the HRF in ~ 87% of cerebral cortex in cognitively intact adults with ages ranging from 22-75 years. The stimulus evokes both positive and negative HRFs, which were characterized using model-free parameters in native-space coordinates. Results showed significant age trends in HRF parameter distributions in terms of both amplitudes (e.g., peak amplitude and CNR) and temporal dynamics (e.g., full-width-at-half-maximum). This work sets the stage for using HRF methods as a biomarker for age-related pathology.

Optimal flip angles for in vivo liver 3D T1 mapping and B1+ mapping at 3T

PURPOSE

The spoiled gradient recalled echo (SPGR) sequence with variable flip angles (FAs) enables whole liverT 1 $$ {T}_1 $$ mapping at high spatial resolutions but is strongly affected byB 1 + $$ {B}_1^{+} $$ inhomogeneities. The aim of this work was to study how the precision of acquiredT 1 $$ {T}_1 $$ maps is affected by theT 1 $$ {T}_1 $$ andB 1 + $$ {B}_1^{+} $$ ranges observed in the liver at 3T, as well as how noise propagates from the acquired signals into the resultingT 1 $$ {T}_1 $$ map.

THEORY

TheT 1 $$ {T}_1 $$ variance was estimated through the Fisher information matrix with a total noise variance including, for the first time, theB 1 + $$ {B}_1^{+} $$ map noise as well as contributions from the SPGR noise.

METHODS

Simulations were used to find the optimal FAs for both theB 1 + $$ {B}_1^{+} $$ mapping andT 1 $$ {T}_1 $$ mapping. The simulations results were validated in 10 volunteers.

RESULTS

Four optimized SPGR FAs of 2°, 2°, 15°, and 15° (TR = 4.1 ms) andB 1 + $$ {B}_1^{+} $$ map FAs of 65° and 130° achieved aT 1 $$ {T}_1 $$ coefficient of variation of 6.2 ± 1.7% across 10 volunteers and validated our theoretical model. Four optimal FAs outperformed five uniformly spaced FAs, saving the patient one breath-hold. For the liverB 1 + $$ {B}_1^{+} $$ andT 1 $$ {T}_1 $$ parameter space at 3T, a higher return inT 1 $$ {T}_1 $$ precision was obtained by investing FAs in the SPGR acquisition rather than in theB 1 + $$ {B}_1^{+} $$ map.

CONCLUSION

A novel framework was developed and validated to calculate the SPGRT 1 $$ {T}_1 $$ variance. This framework efficiently identifies optimal FA values and determines the total number of SPGR andB 1 + $$ {B}_1^{+} $$ measurements needed to achieve a desiredT 1 $$ {T}_1 $$ precision.

Atlas-Based Adaptive Hadamard-Encoded MR Spectroscopic Imaging at 3T

BACKGROUND

This study aimed to develop a time-efficient method of acquiring simultaneous, dual-slice MR spectroscopic imaging (MRSI) for the evaluation of brain metabolism.

METHODS

Adaptive Hadamard-encoded pulses were developed and integrated with atlas-based automatic prescription. The excitation profiles were evaluated via simulation, phantom and volunteer experiments. The feasibility of γ-aminobutyric acid (GABA)-edited dual-slice MRSI was also assessed.

RESULTS

The signal between slices in the dual-band MRSI was less than 1% of the slice profiles. Data from a homemade phantom containing separate, interfacing compartments of creatine and acetate solutions demonstrated ~0.4% acetate signal contamination relative to the amplitude in the excited creatine compartment. The normalized signal-to-noise ratios from atlas-based acquisitions in volunteers were found to be comparable between dual-slice, Hadamard-encoded MRSI and 3D acquisitions. The mean and standard deviation of the coefficients of variation for NAA/Cho from the repeated volunteer scans were 8.2% ± 0.8% and 10.1% ± 3.7% in the top and bottom slices, respectively. GABA-edited, dual-slice MRSI demonstrated simultaneous detection of signals from GABA and coedited macromolecules (GABA+) from both superior grey and deep grey regions of volunteers.

CONCLUSION

This study demonstrated a fully automated dual-slice MRSI acquisition using atlas-based automatic prescription and adaptive Hadamard-encoded pulses.

Improved Simultaneous Multi-slice imaging with Composition of k-space Interpolations (SMS-COOKIE) for myocardial T1 mapping

The aim of this study is to develop and evaluate a regularized Simultaneous Multi-Slice (SMS) reconstruction method for improved Cardiac Magnetic Resonance Imaging (CMR). The proposed reconstruction method, SMS with COmpOsition of k-space IntErpolations (SMS-COOKIE) combines the advantages of Iterative Self-consistent Parallel Imaging Reconstruction (SPIRiT) and split slice-Generalized Autocalibrating Partially Parallel Acquisitions (GRAPPA), while allowing regularization for further noise reduction. The proposed SMS-COOKIE was implemented with and without regularization, and validated using a Saturation Pulse-Prepared Heart rate Independent inversion REcovery (SAPPHIRE) myocardial T1 mapping sequence. The performance of the proposed reconstruction method was compared to ReadOut (RO)-SENSE-GRAPPA and split slice-GRAPPA, on both retrospectively and prospectively three-fold SMS-accelerated data with an additional two-fold in-plane acceleration. All SMS reconstruction methods yielded similar T1 values compared to single band imaging. SMS-COOKIE showed lower spatial variability in myocardial T1 with significant improvement over RO-SENSE-GRAPPA and split slice-GRAPPA (P < 10-4). The proposed method with additional locally low rank (LLR) regularization reduced the spatial variability, again with significant improvement over RO-SENSE-GRAPPA and split slice-GRAPPA (P < 10-4). In conclusion, improved reconstruction quality was achieved with the proposed SMS-COOKIE, which also provided lower spatial variability with significant improvement over split slice-GRAPPA.

1.5-T MR relaxometry in quantifying splenic and pancreatic iron: retrospective comparison of a commercial 3D-Dixon sequence and an established 2D multi-gradient echo sequence

OBJECTIVES

To compare the quantitative measurement of splenic and pancreatic iron content using a commercial 3D-Dixon sequence (qDixon) versus an established fat-saturated R2* relaxometry method (ME-GRE).

METHODS

We analyzed splenic and pancreatic iron levels in 143 MR examinations (1.5 T) using the qDixon and a ME-GRE sequence (108 patients: 65 males, 43 females, mean age 61.31 years). Splenic and pancreatic R2* values were compared between both methods using Bland-Altman plots, concordance correlation coefficients (CCC), and linear regression analyses. Iron overload (R2* > 50 1/s) was defined for both organs and compared using contingency tables, overall agreement, and Gwet's AC1 coefficient.

RESULTS

Of all analyzable examinations, the median splenic R2* using the qDixon sequence was 25.75 1/s (range: 5.6-433) and for the ME-GRE sequence 35.35 1/s (range: 10.9-400.8) respectively. Concerning the pancreas, a median R2* of 29.93 1/s (range: 14-111.45) for the qDixon and 31.25 1/s (range: 14-97) for the ME-GRE sequence was found. Bland-Altman analysis showed a mean R2* difference of 2.12 1/s with a CCC of 0.934 for the spleen and of 0.29 1/s with a CCC of 0.714 for the pancreas. Linear regression for the spleen/pancreas resulted in a correlation coefficient of 0.94 (p < 0.001)/0.725 (p < 0.001). Concerning iron overload, the proportion of overall agreement between the two methods was 91.43% for the spleen and 93.18% for the pancreas.

CONCLUSIONS

Our data show good concordance between R2* values obtained with a commercial qDixon sequence and a validated ME-GRE relaxometry method. The 3D-qDixon sequence, originally intended for liver assessment, seems to be a reliable tool for non-invasive evaluation of iron content also in the spleen and the pancreas.

KEY POINTS

• A 3D chemical shift imaging sequence and 2D multi-gradient echo sequence show good conformity quantifying splenic and pancreatic R2* values. • The 3D chemical shift imaging sequence allows a reliable analysis also of splenic and pancreatic iron status. • In addition to the liver, the analysis of the spleen and pancreas is often helpful for further differential diagnostic clarification and patient guidance regarding the iron status.

Acquisition sequences and reconstruction methods for fast chemical exchange saturation transfer imaging

Chemical exchange saturation transfer (CEST) imaging is an emerging molecular magnetic resonance imaging (MRI) technique that has been developed and employed in numerous diseases. Based on the unique saturation transfer principle, a family of CEST-detectable biomolecules in vivo have been found capable of providing valuable diagnostic information. However, CEST MRI needs a relatively long scan time due to the common long saturation labeling module and typical acquisition of multiple frequency offsets and signal averages, limiting its widespread clinical applications. So far, a plethora of imaging schemes and techniques has been developed to accelerate CEST MRI. In this review, the key acquisition and reconstruction methods for fast CEST imaging are summarized from a practical and systematic point of view. The first acquisition sequence section describes the major development of saturation schemes, readout patterns, ultrafast z-spectroscopy, and saturation-editing techniques for rapid CEST imaging. The second reconstruction method section lists the important advances of parallel imaging, compressed sensing, sparsity in the z-spectrum, and algorithms beyond the Fourier transform for speeding up CEST MRI.

M4Raw: A multi-contrast, multi-repetition, multi-channel MRI k-space dataset for low-field MRI research

Recently, low-field magnetic resonance imaging (MRI) has gained renewed interest to promote MRI accessibility and affordability worldwide. The presented M4Raw dataset aims to facilitate methodology development and reproducible research in this field. The dataset comprises multi-channel brain k-space data collected from 183 healthy volunteers using a 0.3 Tesla whole-body MRI system, and includes T1-weighted, T2-weighted, and fluid attenuated inversion recovery (FLAIR) images with in-plane resolution of ~1.2 mm and through-plane resolution of 5 mm. Importantly, each contrast contains multiple repetitions, which can be used individually or to form multi-repetition averaged images. After excluding motion-corrupted data, the partitioned training and validation subsets contain 1024 and 240 volumes, respectively. To demonstrate the potential utility of this dataset, we trained deep learning models for image denoising and parallel imaging tasks and compared their performance with traditional reconstruction methods. This M4Raw dataset will be valuable for the development of advanced data-driven methods specifically for low-field MRI. It can also serve as a benchmark dataset for general MRI reconstruction algorithms.

Evaluation of simultaneous multi-slice readout-segmented diffusion-weighted MRI acquisition in human optic nerve measurements

Diffusion-weighted magnetic resonance imaging (dMRI) is the only available method to measure the tissue properties of white matter tracts in living human brains and has opened avenues for neuroscientific and clinical studies on human white matter. However, dMRI using conventional simultaneous multi-slice (SMS) single-shot echo planar imaging (ssEPI) still presents challenges in the analyses of some specific white matter tracts, such as the optic nerve, which are heavily affected by susceptibility-induced artifacts. In this study, we evaluated dMRI data acquired by using SMS readout-segmented EPI (rsEPI), which aims to reduce susceptibility-induced artifacts by dividing the acquisition space into multiple segments along the readout direction to reduce echo spacing. To this end, we acquired dMRI data from 11 healthy volunteers by using SMS ssEPI and SMS rsEPI, and then compared the dMRI data of the human optic nerve between the SMS ssEPI and SMS rsEPI datasets by visual inspection of the datasets and statistical comparisons of fractional anisotropy (FA) values. In comparison with the SMS ssEPI data, the SMS rsEPI data showed smaller susceptibility-induced distortion and exhibited a significantly higher FA along the optic nerve. In summary, this study demonstrates that despite its prolonged acquisition time, SMS rsEPI is a promising approach for measuring the tissue properties of the optic nerve in living humans and will be useful for future neuroscientific and clinical investigations of this pathway.

Free-Breathing Liver Fat, R₂* and B₀ Field Mapping Using Multi-Echo Radial FLASH and Regularized Model-Based Reconstruction

This work introduced a stack-of-radial multi-echo asymmetric-echo MRI sequence for free-breathing liver volumetric acquisition. Regularized model-based reconstruction was implemented in Berkeley Advanced Reconstruction Toolbox (BART) to jointly estimate all physical parameter maps (water, fat, R2∗ , and B0 field inhomogeneity maps) and coil sensitivity maps from self-gated k -space data. Specifically, locally low rank and temporal total variation regularization were employed directly on physical parameter maps. The proposed free-breathing radial technique was tested on a water/fat & iron phantom, a young volunteer, and obesity/diabetes/hepatic steatosis patients. Quantitative fat fraction and R2∗ accuracy were confirmed by comparing our technique with the reference breath-hold Cartesian scan. The multi-echo radial sampling sequence achieves fast k -space coverage and is robust to motion. Moreover, the proposed motion-resolved model-based reconstruction allows for free-breathing liver fat and R2∗ quantification in multiple motion states. Overall, our proposed technique offers a convenient tool for non-invasive liver assessment with no breath holding requirement.

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.

Simultaneous Multislice Accelerated TSE for Improved Spatiotemporal Resolution and Diagnostic Accuracy in Magnetic Resonance Neurography: A Feasibility Study

OBJECTIVES

This study aims to evaluate the utility of simultaneous multislice (SMS) acceleration for routine magnetic resonance neurography (MRN) at 3 T.

MATERIALS AND METHODS

Patients with multiple sclerosis underwent MRN of the sciatic nerve consisting of a standard fat-saturated T2-weighted turbo spin echo (TSE) sequence using integrated parallel acquisition technique (PAT2) acceleration and 2 T2 TSE sequences using a combination of PAT-SMS acceleration (1) to reduce scan time (PAT2-SMS2; SMS-TSE FAST ) and (2) for time neutral increase of in-plane resolution (PAT1-SMS2; SMS-TSE HR ). Acquisition times were 5:29 minutes for the standard T2 TSE, 3:12 minutes for the SMS-TSE FAST , and 5:24 minutes for the SMS-TSE HR . Six qualitative imaging parameters were analyzed by 2 blinded readers using a 5-point Likert scale and T2 nerve lesions were quantified, respectively. Qualitative and quantitative image parameters were compared, and both interrater and intrarater reproducibility were statistically assessed. In addition, signal-to-noise ratio/contrast-to-noise ratio (CNR) was obtained in healthy controls using the exact same imaging protocol.

RESULTS

A total of 15 patients with MS (mean age ± standard deviation, 38.1 ± 11 years) and 10 healthy controls (mean age, 29.1 ± 7 years) were enrolled in this study. CNR analysis was highly reliable (intraclass correlation coefficient, 0.755-0.948) and revealed a significant CNR decrease for the sciatic nerve for both SMS protocols compared with standard T2 TSE (SMS-TSE FAST /SMS-TSE HR , -39%/-55%; P ≤ 0.01). Intrarater and interrater reliability of qualitative image review was good to excellent (κ: 0.672-0.971/0.617-0.883). Compared with the standard T2 TSE sequence, both SMS methods were shown to be superior in reducing pulsatile flow artifacts ( P < 0.01). Ratings for muscle border sharpness, detailed muscle structures, nerve border sharpness, and nerve fascicular structure did not differ significantly between the standard T2 TSE and the SMS-TSE FAST ( P > 0.05) and were significantly better for the SMS-TSE HR than for standard T2 TSE ( P < 0.001). Muscle signal homogeneity was mildly inferior for both SMS-TSE FAST ( P > 0.05) and SMS-TSE HR ( P < 0.001). A significantly higher number of T2 nerve lesions were detected by SMS-TSE HR ( P ≤ 0.01) compared with the standard T2 TSE and SMS-TSE FAST , whereas no significant difference was observed between the standard T2 TSE and SMS-TSE FAST .

CONCLUSIONS

Implementation of SMS offers either to substantially reduce acquisition time by over 40% without significantly impeding image quality compared with the standard T2 TSE or to increase in-plane resolution for a high-resolution approach and improved depiction of T2 nerve lesions while keeping acquisition times constant. This addresses the specific needs of MRN by providing different imaging approaches for 2D clinical MRN.

Quantitative myocardial perfusion with a hybrid 2D simultaneous multi-slice sequence

PURPOSE

To evaluate a novel 2D simultaneous multi-slice (SMS) myocardial perfusion acquisition and compare directly to a published quantitative 3D stack-of-stars (SoS) acquisition.

METHODS

A hybrid saturation recovery radial 2D SMS sequence following a single saturation was created for the quantification of myocardial blood flow (MBF). This sequence acquired three slices simultaneously and generated an arterial input function (AIF) using the first 24 rays. Validation was done in a novel way by alternating heartbeats between the hybrid 2D SMS and the 3D SoS acquisitions. Initial studies were done to study the effects of using only every other beat for the 2D SMS in two subjects, and for the 3D SoS in four subjects. The proposed alternating acquisitions were then performed in ten dog studies at rest, four dog studies at adenosine stress, and two human resting studies. Quantitative MBF analysis was performed for 2D SMS and 3D SoS separately, using a compartment model.

RESULTS

Acquiring every-other-beat data resulted in 6 ± 5% ("ideal") and 11 ± 8% ("practical") perfusion changes for both 2D SMS and 3D SoS methods. For alternating acquisitions, 2D SMS and 3D SoS quantitative perfusion values were comparable for both the twelve rest studies (2D SMS: 0.69 ± 0.16 vs 3D: 0.69 ± 0.15 ml/g/min, p = 0.55) and the four stress studies (2D SMS: 1.28 ± 0.22 vs 3D: 1.30 ± 0.24 ml/g/min, p = 0.61).

CONCLUSION

Every-other-beat acquisition changed estimated perfusion values relatively little for both sequences. The quantitative hybrid radial 2D SMS myocardial first-pass perfusion imaging sequence gave results similar to 3D perfusion when compared directly with an alternating beat acquisition.

Simultaneous multislice imaging with slice-specific z-shim

PURPOSE

To implement slice-specific z-shim in simultaneous multislice (SMS) imaging in order to minimize signal losses in slice-accelerated T₂ *-weighted acquisitions, such as for spinal cord functional neuroimaging.

METHODS

The RF envelopes of the individual slice bands are temporally shifted on the plateau of the slice-selection gradient pulse before being combined to the multiband RF envelope. Thus, optimum z-shims can be realized for each slice of an SMS excitation, which is in contrast to conventional z-shimming. EPI with 2-fold SMS acceleration was performed on a 3T whole-body MR system in phantoms and the cervical spinal cord of healthy volunteers (i) without z-shim, (ii) with conventional z-shim using the average value of the slices of the SMS excitation, and (iii) with optimal, slice-specific z-shims for each slice using envelope shifts.

RESULTS

Phantom experiments demonstrate the equivalence of the envelope shift and conventional z-shimming for non-SMS excitations. With SMS, the best image quality is obtained with "mixed" z-shim, where only the z-shim differences of the slices of an SMS excitation are implemented by an envelope shift while their mean z-shim is applied conventionally with a gradient pulse after the echoes acquired for N/2 ghost correction. In phantoms and in vivo, this setup outperforms the approaches without slice-specific z-shim with respect to signal amplitude and temporal SNR at the expense of slight TE differences (<1 ms) between the slices.

CONCLUSION

With RF envelope shifts, slice-specific z-shims can be combined with SMS imaging, which could improve slice-accelerated functional neuroimaging in the spinal cord.

Magnetic resonance imaging at 9.4 T: the Maastricht journey

The 9.4 T scanner in Maastricht is a whole-body magnet with head gradients and parallel RF transmit capability. At the time of the design, it was conceptualized to be one of the best fMRI scanners in the world, but it has also been used for anatomical and diffusion imaging. 9.4 T offers increases in sensitivity and contrast, but the technical ultra-high field (UHF) challenges, such as field inhomogeneities and constraints set by RF power deposition, are exacerbated compared to 7 T. This article reviews some of the 9.4 T work done in Maastricht. Functional imaging experiments included blood oxygenation level-dependent (BOLD) and blood-volume weighted (VASO) fMRI using different readouts. BOLD benefits from shorter T₂* at 9.4 T while VASO from longer T₁. We show examples of both ex vivo and in vivo anatomical imaging. For many applications, pTx and optimized coils are essential to harness the full potential of 9.4 T. Our experience shows that, while considerable effort was required compared to our 7 T scanner, we could obtain high-quality anatomical and functional data, which illustrates the potential of MR acquisitions at even higher field strengths. The practical challenges of working with a relatively unique system are also discussed.

MR-self Noise2Noise: self-supervised deep learning-based image quality improvement of submillimeter resolution 3D MR images

OBJECTIVES

The study aimed to develop a deep neural network (DNN)-based noise reduction and image quality improvement by only using routine clinical scans and evaluate its performance in 3D high-resolution MRI.

METHODS

This retrospective study included T1-weighted magnetization-prepared rapid gradient-echo (MP-RAGE) images from 185 clinical scans: 135 for DNN training, 11 for DNN validation, 20 for qualitative evaluation, and 19 for quantitative evaluation. Additionally, 18 vessel wall imaging (VWI) data were included to evaluate generalization. In each scan of the DNN training set, two noise-independent images were generated from the k-space data, resulting in an input-label pair. 2.5D U-net architecture was utilized for the DNN model. Qualitative evaluation between conventional MP-RAGE and DNN-based MP-RAGE was performed by two radiologists in image quality, fine structure delineation, and lesion conspicuity. Quantitative evaluation was performed with full sampled data as a reference by measuring quantitative error metrics and volumetry at 7 different simulated noise levels. DNN application on VWI was evaluated by two radiologists in image quality.

RESULTS

Our DNN-based MP-RAGE outperformed conventional MP-RAGE in all image quality parameters (average scores = 3.7 vs. 4.9, p < 0.001). In the quantitative evaluation, DNN showed better error metrics (p < 0.001) and comparable (p > 0.09) or better (p < 0.02) volumetry results than conventional MP-RAGE. DNN application to VWI also revealed improved image quality (3.5 vs. 4.6, p < 0.001).

CONCLUSION

The proposed DNN model successfully denoises 3D MR image and improves its image quality by using routine clinical scans only.

KEY POINTS

• Our deep learning framework successfully improved conventional 3D high-resolution MRI in all image quality parameters, fine structure delineation, and lesion conspicuity. • Compared to conventional MRI, the proposed deep neural network-based MRI revealed better quantitative error metrics and comparable or better volumetry results. • Deep neural network application to 3D MRI whose pulse sequences and parameters were different from the training data showed improvement in image quality, revealing the potential to generalize on various clinical MRI.

Multi-readout DWI with a reduced FOV for studying the coupling between diffusion and T 2 * relaxation in the prostate

PURPOSE

To develop a DWI sequence with multiple readout echo-trains in a single shot (multi-readout DWI) over a reduced FOV, and to demonstrate its ability to achieve high data acquisition efficiency in the study of coupling between diffusion and relaxation in the human prostate.

METHODS

The proposed multi-readout DWI sequence plays out multiple EPI readout echo-trains after a Stejskal-Tanner diffusion preparation module. Each EPI readout echo-train corresponded to a distinct effective TE. To maintain a high spatial resolution with a relatively short echo-train for each readout, a 2D RF pulse was used to limit the FOV. Experiments were performed on the prostate of six healthy subjects to acquire a set of images with three b values (0, 500, and 1000 s/mm² ) and three TEs (63.0, 78.8, and 94.6 ms), producing three ADC maps at different TEs and three T 2 * $$ {T}_2^{\ast } $$ maps at different b values.

RESULTS

Multi-readout DWI enabled a threefold acceleration without compromising the spatial resolution when compared with a conventional single-readout sequence. Images with three b values and three TEs were obtained in 3 min 40 s with an adequate SNR (≥ 26.9). The ADC values (1.45 ± 0.13, 1.52 ± 0.14, and 1.58 ± 0.15  μm 2 / ms $$ {\upmu \mathrm{m}}^2/\mathrm{ms} $$ ; P < 0.01) exhibited an increasing trend as TEs increased (63.0 ms, 78.8 ms, and 94.6 ms), whereas T 2 * $$ {T}_2^{\ast } $$ values (74.78 ± 13.21, 63.21 ± 7.84, and 56.61 ± 5.05 ms; P < 0.01) decreases as the b values increased (0, 500, and 1000 s/mm² ).

CONCLUSION

The multi-readout DWI sequence over a reduced FOV provides a time-efficient technique to study the coupling between diffusion and relaxation times.

Alteration of brain function and systemic inflammatory tone in older adults by decreasing the dietary palmitic acid intake

Prior studies in younger adults showed that reducing the normally high intake of the saturated fatty acid, palmitic acid (PA), in the North American diet by replacing it with the monounsaturated fatty acid, oleic acid (OA), decreased blood concentrations and secretion by peripheral blood mononuclear cells (PBMCs) of interleukin (IL)-1β and IL-6 and changed brain activation in regions of the working memory network. We examined the effects of these fatty acid manipulations in the diet of older adults. Ten subjects, aged 65-75 years, participated in a randomized, cross-over trial comparing 1-week high PA versus low PA/high OA diets. We evaluated functional magnetic resonance imaging (fMRI) using an N-back test of working memory and a resting state scan, cytokine secretion by lipopolysaccharide (LPS)-stimulated PBMCs, and plasma cytokine concentrations. During the low PA compared to the high PA diet, we observed increased activation for the 2-back minus 0-back conditions in the right dorsolateral prefrontal cortex (Broadman Area (BA) 9; p < 0.005), but the effect of diet on working memory performance was not significant (p = 0.09). We observed increased connectivity between anterior regions of the salience network during the low PA/high OA diet (p < 0.001). The concentrations of IL-1β (p = 0.026), IL-8 (p = 0.013), and IL-6 (p = 0.009) in conditioned media from LPS-stimulated PBMCs were lower during the low PA/high OA diet. This study suggests that lowering the dietary intake of PA down-regulated pro-inflammatory cytokine secretion and altered working memory, task-based activation and resting state functional connectivity in older adults.

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.

Simultaneous multi-segment (SMSeg) EPI over multiple focal regions

Objective.This study aimed at developing a simultaneous multi-segment (SMSeg) imaging technique using a two-dimensional (2D) RF pulse in conjunction with echo planar imaging (EPI) to image multiple focal regions.Approach.The SMSeg technique leveraged periodic replicates of the excitation profile of a 2D RF pulse to simultaneously excite multiple focal regions at different locations. These locations were controlled by rotating and scaling transmit k-space trajectories. The resulting multiple isolated focal regions were projected into a composite 'slice' for display. GRAPPA-based parallel imaging was incorporated into SMSeg by taking advantage of coil sensitivity variations in both the phase-encoded and slice-selection directions. The SMSeg technique was implemented at 3 T in a single-shot gradient-echo EPI sequence and demonstrated in a phantom and human brains for both anatomic imaging and functional imaging.Main results.In both the phantom and the human brain, SMSeg images from three focal regions were simultaneously acquired. SMSeg imaging enabled up to a six-fold acceleration in parallel imaging without causing appreciable residual aliasing artifacts when compared with a conventional gradient-echo EPI sequence with the same acceleration factor. In the functional imaging experiment, BOLD activations associated with a visuomotor task were simultaneously detected in two non-coplanar segments (each with a size of 240 × 30 mm²), corresponding to visual and motor cortices, respectively.Significance.Our study has demonstrated that SMSeg imaging can be a viable method for studying multiple focal regions simultaneously.

Clinical Evaluation of Scout Accelerated Motion Estimation and Reduction Technique for 3D MR Imaging in the Inpatient and Emergency Department Settings

BACKGROUND AND PURPOSE

A scout accelerated motion estimation and reduction (SAMER) framework has been developed for efficient retrospective motion correction. The goal of this study was to perform an initial evaluation of SAMER in a series of clinical brain MR imaging examinations.

MATERIALS AND METHODS

Ninety-seven patients who underwent MR imaging in the inpatient and emergency department settings were included in the study. SAMER motion correction was retrospectively applied to an accelerated T1-weighted MPRAGE sequence that was included in brain MR imaging examinations performed with and without contrast. Two blinded neuroradiologists graded images with and without SAMER motion correction on a 5-tier motion severity scale (none = 1, minimal = 2, mild = 3, moderate = 4, severe = 5).

RESULTS

The median SAMER reconstruction time was 1 minute 47 seconds. SAMER motion correction significantly improved overall motion grades across all examinations (P < .005). Motion artifacts were reduced in 28% of cases, unchanged in 64% of cases, and increased in 8% of cases. SAMER improved motion grades in 100% of moderate motion cases and 75% of severe motion cases. Sixty-nine percent of nondiagnostic motion cases (grades 4 and 5) were considered diagnostic after SAMER motion correction. For cases with minimal or no motion, SAMER had negligible impact on the overall motion grade. For cases with mild, moderate, and severe motion, SAMER improved the motion grade by an average of 0.3 (SD, 0.5), 1.1 (SD, 0.3), and 1.1 (SD, 0.8) grades, respectively.

CONCLUSIONS

SAMER improved the diagnostic image quality of clinical brain MR imaging examinations with motion artifacts. The improvement was most pronounced for cases with moderate or severe motion.

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

BACKGROUND

This study evaluates the possibility for replacing conventional 3 slices, 3 breath-holds MOLLI cardiac T1 mapping with single breath-hold 3 simultaneous multi-slice (SMS3) T1 mapping using blipped-CAIPIRINHA SMS-bSSFP MOLLI sequence. As a major drawback, SMS-bSSFP presents unique artefacts arising from side-lobe slice excitations that are explained by imperfect RF modulation rendering and bSSFP low flip angle enhancement. Amplitude-only RF modulation (AM) is proposed to reduce these artefacts in SMS-MOLLI compared to conventional Wong multi-band RF modulation (WM).

MATERIALS AND METHODS

Phantoms and ten healthy volunteers were imaged at 1.5 T using a modified blipped-CAIPIRINHA SMS-bSSFP MOLLI sequence with 3 simultaneous slices. WM-SMS3 and AM-SMS3 were compared to conventional single-slice (SMS1) MOLLI. First, SNR degradation and T1 accuracy were measured in phantoms. Second, artefacts from side-lobe excitations were evaluated in a phantom designed to reproduce fat presence near the heart. Third, the occurrence of these artefacts was observed in volunteers, and their impact on T1 quantification was compared between WM-SMS3 and AM-SMS3 with conventional MOLLI as a reference.

RESULTS

In the phantom, larger slice gaps and slice thicknesses yielded higher SNR. There was no significant difference of T1 values between conventional MOLLI and SMS3-MOLLI (both WM and AM). Positive banding artefacts were identified from fat neighbouring the targeted FOV due to side-lobe excitations from WM and the unique bSSFP signal profile. AM RF pulses reduced these artefacts by 38%. In healthy volunteers, AM-SMS3-MOLLI showed similar artefact reduction compared to WM-SMS3-MOLLI (3 ± 2 vs 5 ± 3 corrupted LV segments out of 16). In-vivo native T1 values obtained from conventional MOLLI and AM-SMS3-MOLLI were equivalent in LV myocardium (SMS1-T1 = 935.5 ± 36.1 ms; AM-SMS3-T1 = 933.8 ± 50.2 ms; P = 0.436) and LV blood pool (SMS1-T1 = 1475.4 ± 35.9 ms; AM-SMS3-T1 = 1452.5 ± 70.3 ms; P = 0.515). Identically, no differences were found between SMS1 and SMS3 postcontrast T1 values in the myocardium (SMS1-T1 = 556.0 ± 19.7 ms; SMS3-T1 = 521.3 ± 28.1 ms; P = 0.626) and the blood (SMS1-T1 = 478 ± 65.1 ms; AM-SMS3-T1 = 447.8 ± 81.5; P = 0.085).

CONCLUSIONS

Compared to WM RF modulation, AM SMS-bSSFP MOLLI was able to reduce side-lobe artefacts considerably, providing promising results to image the three levels of the heart in a single breath hold. However, few artefacts remained even using AM-SMS-bSSFP due to residual RF imperfections. The proposed blipped-CAIPIRINHA MOLLI T1 mapping sequence provides accurate in vivo T1 quantification in line with those obtained with a single slice acquisition.

Signal intensity informed multi-coil encoding operator for physics-guided deep learning reconstruction of highly accelerated myocardial perfusion CMR

PURPOSE

To develop a physics-guided deep learning (PG-DL) reconstruction strategy based on a signal intensity informed multi-coil (SIIM) encoding operator for highly-accelerated simultaneous multislice (SMS) myocardial perfusion cardiac MRI (CMR).

METHODS

First-pass perfusion CMR acquires highly-accelerated images with dynamically varying signal intensity/SNR following the administration of a gadolinium-based contrast agent. Thus, using PG-DL reconstruction with a conventional multi-coil encoding operator leads to analogous signal intensity variations across different time-frames at the network output, creating difficulties in generalization for varying SNR levels. We propose to use a SIIM encoding operator to capture the signal intensity/SNR variations across time-frames in a reformulated encoding operator. This leads to a more uniform/flat contrast at the output of the PG-DL network, facilitating generalizability across time-frames. PG-DL reconstruction with the proposed SIIM encoding operator is compared to PG-DL with conventional encoding operator, split slice-GRAPPA, locally low-rank (LLR) regularized reconstruction, low-rank plus sparse (L + S) reconstruction, and regularized ROCK-SPIRiT.

RESULTS

Results on highly accelerated free-breathing first pass myocardial perfusion CMR at three-fold SMS and four-fold in-plane acceleration show that the proposed method improves upon the reconstruction methods use for comparison. Substantial noise reduction is achieved compared to split slice-GRAPPA, and aliasing artifacts reduction compared to LLR regularized reconstruction, L + S reconstruction and PG-DL with conventional encoding. Furthermore, a qualitative reader study indicated that proposed method outperformed all methods.

CONCLUSION

PG-DL reconstruction with the proposed SIIM encoding operator improves generalization across different time-frames /SNRs in highly accelerated perfusion CMR.

Myocardial First-Pass Perfusion With Increased Anatomic Coverage at 3 T Using Autocalibrated Multiband Imaging

BACKGROUND

Myocardial first-pass perfusion (FPP) imaging is a useful cardiac MRI method for the diagnosis of coronary artery disease. However, conventional 2D multislice FPP acquisitions usually have gaps between myocardium slices, which limits the overall assessment of myocardial ischemia.

PURPOSE

To increase the anatomic coverage of myocardial FPP imaging at 3 T by implementing both autocalibrated multiband (MB) acquisition and k-t space acceleration with compress sensing (CS) reconstruction, without the need for additional reference scans.

STUDY TYPE

Phantom and prospective human studies.

PHANTOM/SUBJECTS

A T1MES (T1 Mapping and ECV Standardization in cardiovascular magnetic resonance) phantom and 20 subjects (12 healthy subjects and 8 patients, 10 males, age 42 ± 16 years).

FIELD STRENGTH/SEQUENCE

A 3 T/saturation recovery prepared gradient echo sequence with contrast administration.

ASSESSMENT

Phantom experiments were performed to compare the performance of autocalibrated MB-FPP with k-t acceleration using slice-GRAPPA and CS reconstructions. In vivo experiments were performed to compare the performance of conventional FPP (2.5× acceleration) with autocalibrated MB + CS-FPP (6× acceleration). In phantom experiments, the error maps were calculated. In in vivo experiments, the contrast ratio (CR) and blurring were quantitatively measured, while image quality, perceived signal-to-noise ratio (SNR), and artifact level were qualitatively graded by three cardiologists on a 4-point scale.

STATISTICAL TESTS

Wilcoxon signed-rank test, paired t-test. A P value <0.05 was considered statistically significant.

RESULTS

In phantom experiments, residual artifact was reduced using the MB + CS-FPP reconstruction method compared with using the MB + slice-GRAPPA reconstruction method. In in vivo experiments, the proposed autocalibrated MB + CS-FPP method demonstrated significantly higher CR (3.52 ± 0.78 vs 2.91 ± 0.81) and had significantly better perceived SNR (2.69 ± 0.29 vs 2.48 ± 0.31) compared to the conventional sequence. Compared with conventional FPP, MB + CS-FPP doubled the spatial coverage (MB + CS-FPP vs conventional FPP) without compromising the image quality (2.69 ± 0.26 vs 2.60 ± 0.30) or increasing the artifact level (2.60 ± 0.26 vs 2.52 ± 0.31).

CONCLUSION

Autocalibrated MB + CS-FPP improved the myocardial coverage and achieved comparable image quality with the same spatial resolution and scan time as conventional FPP and is a promising technique for clinical myocardial perfusion imaging.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 3.

Framework for Accurate Classification of Self-Reported Stress From Multisession Functional MRI Data of Veterans With Posttraumatic Stress

Background: Posttraumatic stress disorder (PTSD) is a significant burden among combat Veterans returning from the wars in Iraq and Afghanistan. While empirically supported treatments have demonstrated reductions in PTSD symptomatology, there remains a need to improve treatment effectiveness. Functional magnetic resonance imaging (fMRI) neurofeedback has emerged as a possible treatment to ameliorate PTSD symptom severity. Virtual reality (VR) approaches have also shown promise in increasing treatment compliance and outcomes. To facilitate fMRI neurofeedback-associated therapies, it would be advantageous to accurately classify internal brain stress levels while Veterans are exposed to trauma-associated VR imagery. Methods: Across 2 sessions, we used fMRI to collect neural responses to trauma-associated VR-like stimuli among male combat Veterans with PTSD symptoms (N = 8). Veterans reported their self-perceived stress level on a scale from 1 to 8 every 15 s throughout the fMRI sessions. In our proposed framework, we precisely sample the fMRI data on cortical gray matter, blurring the data along the gray-matter manifold to reduce noise and dimensionality while preserving maximum neural information. Then, we independently applied 3 machine learning (ML) algorithms to this fMRI data collected across 2 sessions, separately for each Veteran, to build individualized ML models that predicted their internal brain states (self-reported stress responses). Results: We accurately classified the 8-class self-reported stress responses with a mean (± standard error) root mean square error of 0.6 (± 0.1) across all Veterans using the best ML approach. Conclusions: The findings demonstrate the predictive ability of ML algorithms applied to whole-brain cortical fMRI data collected during individual Veteran sessions. The framework we have developed to preprocess whole-brain cortical fMRI data and train ML models across sessions would provide a valuable tool to enable individualized real-time fMRI neurofeedback during VR-like exposure therapy for PTSD.

Wave-controlled aliasing in parallel imaging (Wave-CAIPI): Accelerating speed for the MRI-based diagnosis of enhancing intracranial lesions compared to magnetization-prepared gradient echo

PURPOSE

We aimed to validate the diagnostic performance of accelerated post-contrast magnetization-prepared rapid gradient-echo (MPRAGE) using wave-controlled aliasing in parallel imaging (Wave-CAIPI) for enhancing intracranial lesions, compared with conventional MPRAGE.

METHODS

A total of 233 consecutive patients who underwent post-contrast Wave-CAIPI and conventional MPRAGE (scan time: 2 min 39 s vs. 4 min 30 s) were retrospectively evaluated. Two radiologists independently assessed whole images for the presence and diagnosis of enhancing lesions. The diagnostic performance for non-enhancing lesions, quantitative parameters (diameter of enhancing lesions, signal-to-noise ratio [SNR], contrast-to-noise ratio [CNR], and contrast rate), qualitative parameters (grey-white matter differentiation and conspicuity of enhancing lesions), and image qualities (overall image quality and motion artifacts) were also surveyed. The weighted kappa and percent agreement were used to evaluate the diagnostic agreement between the two sequences.

RESULTS

Wave-CAIPI MPRAGE achieved significantly high agreement for the detection (98.7%[460/466], κ = 0.965) and diagnosis (97.8%[455/466], κ = 0.955) of enhancing intracranial lesions with conventional MPRAGE in pooled analysis. Detection and diagnosis of non-enhancing lesions (97.6% and 96.9% agreement), and diameter of enhancing lesions (P>0.05) also demonstrated high agreements between two sequences. Although Wave-CAIPI MPRAGE show lower SNR (P<0.01) than conventional MRAGE, it fulfilled comparable CNR (P = 0.486) and higher contrast rate (P<0.01). The qualitative parameters show similar value (P>0.05). The overall image quality was slightly poor, however, motion artifacts were better in Wave-CAIPI MPRAGE (both P = 0.005).

CONCLUSION

Wave-CAIPI MPRAGE provides reliable diagnostic performance for enhancing intracranial lesions within half of the scan time compared with conventional MPRAGE.

Diagnosis of intracranial lesions using accelerated 3D T1 MPRAGE with wave-CAIPI technique: comparison with conventional 3D T1 MPRAGE

We aimed to evaluate the agreement in the diagnosis of intracranial lesions between conventional pre-contrast 3D T1 magnetization-prepared rapid gradient echo (MPRAGE) and wave-CAIPI (wave-controlled aliasing in parallel imaging) MPRAGE. Institutional review board approval was obtained and informed consent was waived for this retrospective study. We included 149 consecutive patients who had undergone brain MR with both conventional MPRAGE (scan time: 5 min 42 s) and wave-CAIPI MPRAGE (scan time: 2 min 44 s) from February to June 2018. All images were independently reviewed by two radiologists for the diagnosis of intracranial lesion and scored image quality using visual analysis. One technician measured signal-to-noise ratio. The agreement for diagnosis of intracranial lesion was calculated, and the intra- and interobserver agreements were analyzed by using kappa value. For the diagnosis of intracranial lesion, the conventional and wave-CAIPI MPRAGE demonstrated 99.7% of agreement (297 of 298) in the pooled analysis with very good agreement (k = 0.994). Intra- and inter-observer agreement showed very good (k > 0.9 in all) and good (k > 0.75) agreement, respectively. In the quantitative analysis, the signal-to-noise ratio had no difference (P > 0.05 for all). The overall image quality was poorer in images of wave-CAIPI MPRAGE (P < 0.001), but motion artifact had no difference between two sequences (P = 0.06). Compared to conventional MPRAGE, pre-contrast 3D T1 wave-CAIPI MPRAGE achieved higher agreement for the diagnosis of intracranial lesions and reduced the scan time by approximately 50%.

Wave-Encoded Model-Based Deep Learning for Highly Accelerated Imaging with Joint Reconstruction

A recently introduced model-based deep learning (MoDL) technique successfully incorporates convolutional neural network (CNN)-based regularizers into physics-based parallel imaging reconstruction using a small number of network parameters. Wave-controlled aliasing in parallel imaging (CAIPI) is an emerging parallel imaging method that accelerates imaging acquisition by employing sinusoidal gradients in the phase- and slice/partition-encoding directions during the readout to take better advantage of 3D coil sensitivity profiles. We propose wave-encoded MoDL (wave-MoDL) combining the wave-encoding strategy with unrolled network constraints for highly accelerated 3D imaging while enforcing data consistency. We extend wave-MoDL to reconstruct multicontrast data with CAIPI sampling patterns to leverage similarity between multiple images to improve the reconstruction quality. We further exploit this to enable rapid quantitative imaging using an interleaved look-locker acquisition sequence with T₂ preparation pulse (3D-QALAS). Wave-MoDL enables a 40 s MPRAGE acquisition at 1 mm resolution at 16-fold acceleration. For quantitative imaging, wave-MoDL permits a 1:50 min acquisition for T₁, T₂, and proton density mapping at 1 mm resolution at 12-fold acceleration, from which contrast-weighted images can be synthesized as well. In conclusion, wave-MoDL allows rapid MR acquisition and high-fidelity image reconstruction and may facilitate clinical and neuroscientific applications by incorporating unrolled neural networks into wave-CAIPI reconstruction.

Pediatric magnetic resonance imaging: faster is better

Magnetic resonance imaging (MRI) has emerged as the preferred imaging modality for evaluating a wide range of pediatric medical conditions. Nevertheless, the long acquisition times associated with this technique can limit its widespread use in young children, resulting in motion-degraded or non-diagnostic studies. As a result, sedation or general anesthesia is often necessary to obtain diagnostic images, which has implications for the safety profile of MRI, the cost of the exam and the radiology department's clinical workflow. Over the last decade, several techniques have been developed to increase the speed of MRI, including parallel imaging, single-shot acquisition, controlled aliasing techniques, compressed sensing and artificial-intelligence-based reconstructions. These are advantageous because shorter examinations decrease the need for sedation and the severity of motion artifacts, increase scanner throughput, and improve system efficiency. In this review we discuss a framework for image acceleration in children that includes the synergistic use of state-of-the-art MRI hardware and optimized pulse sequences. The discussion is framed within the context of pediatric radiology and incorporates the authors' experience in deploying these techniques in routine clinical practice.

Generalized self-calibrating simultaneous multi-slice MR image reconstruction from 3D Fourier encoding perspective

This work introduces a novel, k-space based one-step solution for simultaneous multi-slice MR image reconstruction from 3D Fourier encoding perspective. With undersampled SMS imaging, image reconstruction suffers from both inter-slice leakages and in-plane aliasing artifacts. Aliasing separation becomes further challenging in the presence of discrepancies between calibration and imaging. To address them, in this work a measured SMS 3D k-space with additional calibrating signals is decomposed into SMS imaging and self-calibrating data sets. Extended controlled aliasing is performed by upsampling the measured data in the k_z-direction. A slice-specific null space operator is then learned using extended self-calibration exploiting target slices and additional in-plane-shifted images. Inter-slice leakages and in-plane aliasing artifacts are jointly resolved in a single step by solving a constrained optimization problem in which null space reconstruction consistency is balanced with a Hankel-structured low rank prior while data fidelity in 3D Fourier space is enforced. Retrospective and prospective studies are performed to validate the effectiveness of the proposed method in various regions including knee and L-spine.

Open and reproducible neuroimaging: From study inception to publication

Empirical observations of how labs conduct research indicate that the adoption rate of open practices for transparent, reproducible, and collaborative science remains in its infancy. This is at odds with the overwhelming evidence for the necessity of these practices and their benefits for individual researchers, scientific progress, and society in general. To date, information required for implementing open science practices throughout the different steps of a research project is scattered among many different sources. Even experienced researchers in the topic find it hard to navigate the ecosystem of tools and to make sustainable choices. Here, we provide an integrated overview of community-developed resources that can support collaborative, open, reproducible, replicable, robust and generalizable neuroimaging throughout the entire research cycle from inception to publication and across different neuroimaging modalities. We review tools and practices supporting study inception and planning, data acquisition, research data management, data processing and analysis, and research dissemination. An online version of this resource can be found at https://oreoni.github.io. We believe it will prove helpful for researchers and institutions to make a successful and sustainable move towards open and reproducible science and to eventually take an active role in its future development.

Customized whole brain-covering 3D GRASE in multi-delay pseudo-continuous arterial spin labeling for duplex distinct hemodynamic mapping contrasts of brain tissues and circulation pathways

Gradient and spin echo (GRASE) is widely employed in arterial spin labeling (ASL) as an efficient readout sequence. Hemodynamic parameter mappings of perfusion, such as cerebral blood flow (CBF) and arterial transit time (ATT), can be derived via multi-delay ASL acquisitions. Multi-delay ASL perfusion imaging inevitably suffers limited signal-to-noise ratio (SNR) since a motion-sensitized vessel suppressing module has to be employed to highlight perfusion signals. The present work reveals that in multi-delay ASL, manipulation of GRASE sequence on either planar imaging echo echo train for adjusted spatial resolutions or FSE echo train for modulated extent of T 2-blurring can significantly alter the mapping contrasts among tissues and among cerebral lobes under different pathways of blood circulation, and meanwhile regulates SNR. Four separate multi-delay ASL scans with different echo train designs in 3D whole brain covering GRASE were carried out for healthy subjects to evaluate the variations in regard to the parameter quantifications and SNR. Based on the quantification mappings, the GRASE acquisition with moderate spatial resolution (3.5 × 3.5 × 4 mm3) and segmented k z scheme was recognized for the first time to be recommended for more unambiguous CBF and ATT contrasts between GM and WM in conjunction with more enhanced ATT contrast between anterior and posterior cerebral circulations, with reasonably good SNR. The technical proposal is of great value for the cutting-edge research of a variety of neurological diseases of global concerns.

Highly accelerated EPI with wave encoding and multi-shot simultaneous multislice imaging

PURPOSE

To introduce wave-encoded acquisition and reconstruction techniques for highly accelerated EPI with reduced g-factor penalty and image artifacts.

THEORY AND METHODS

Wave-EPI involves application of sinusoidal gradients during the EPI readout, which spreads the aliasing in all spatial directions, thereby taking better advantage of 3D coil sensitivity profiles. The amount of voxel spreading that can be achieved by the wave gradients during the short EPI readout period is constrained by the slew rate of the gradient coils and peripheral nerve stimulation monitor. We propose to use a "half-cycle" sinusoidal gradient to increase the amount of voxel spreading that can be achieved while respecting the slew and stimulation constraints. Extending wave-EPI to multi-shot acquisition minimizes geometric distortion and voxel blurring at high in-plane resolutions, while structured low-rank regularization mitigates shot-to-shot phase variations. To address gradient imperfections, we propose to use different point spread functions for the k-space lines with positive and negative polarities, which are calibrated with a FLEET-based reference scan.

RESULTS

Wave-EPI enabled whole-brain single-shot gradient-echo (GE) and multi-shot spin-echo (SE) EPI acquisitions at high acceleration factors at 3T and was combined with g-Slider encoding to boost the SNR level in 1 mm isotropic diffusion imaging. Relative to blipped-CAIPI, wave-EPI reduced average and maximum g-factors by up to 1.21- and 1.37-fold at R_in  × R_sms  = 3 × 3, respectively.

CONCLUSION

Wave-EPI allows highly accelerated single- and multi-shot EPI with reduced g-factor and artifacts and may facilitate clinical and neuroscientific applications of EPI by improving the spatial and temporal resolution in functional and diffusion imaging.