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.

Single Breath-Hold 3-Dimensional Magnetic Resonance Elastography Depicts Liver Fibrosis and Inflammation in Obese Patients

OBJECTIVES

Three-dimensional (3D) magnetic resonance elastography (MRE) measures liver fibrosis and inflammation but requires several breath-holds that hamper clinical acceptance. The aim of this study was to evaluate the technical and clinical feasibility of a single breath-hold 3D MRE sequence as a means of measuring liver fibrosis and inflammation in obese patients.

METHODS

From November 2020 to December 2021, subjects were prospectively enrolled and divided into 2 groups. Group 1 included healthy volunteers (n = 10) who served as controls to compare the single breath-hold 3D MRE sequence with a multiple-breath-hold 3D MRE sequence. Group 2 included liver patients (n = 10) who served as participants to evaluate the clinical feasibility of the single breath-hold 3D MRE sequence in measuring liver fibrosis and inflammation. Controls and participants were scanned at 60 Hz mechanical excitation with the single breath-hold 3D MRE sequence to retrieve the magnitude of the complex-valued shear modulus (|G*| [kPa]), the shear wave speed (Cs [m/s]), and the loss modulus (G" [kPa]). The controls were also scanned with a multiple-breath-hold 3D MRE sequence for comparison, and the participants had histopathology (Ishak scores) for correlation with Cs and G".

RESULTS

For the 10 controls, 5 were female, and the mean age and body mass index were 33.1 ± 9.5 years and 23.0 ± 2.1 kg/m 2 , respectively. For the 10 participants, 8 were female, and the mean age and body mass index were 45.1 ± 16.5 years and 33.1 ± 4.0 kg/m 2 (obese range), respectively. All participants were suspected of having nonalcoholic fatty liver disease. Bland-Altman analysis of the comparison in controls shows there are nonsignificant differences in |G*|, Cs, and G" below 6.5%, suggesting good consensus between the 2 sequences. For the participants, Cs and G" correlated significantly with Ishak fibrosis and inflammation grades, respectively ( ρ = 0.95, P < 0.001, and ρ = 0.84, P = 0.002).

CONCLUSION

The single breath-hold 3D MRE sequence may be effective in measuring liver fibrosis and inflammation in obese patients.

Reconstruction for 7T high-resolution whole-brain diffusion MRI using two-stage N/2 ghost correction and L1-SPIRiT without single-band reference

PURPOSE

To combine a new two-stage N/2 ghost correction and an adapted L1-SPIRiT method for reconstruction of 7T highly accelerated whole-brain diffusion MRI (dMRI) using only autocalibration scans (ACS) without the need of additional single-band reference (SBref) scans.

METHODS

The proposed ghost correction consisted of a 3-line reference approach in stage 1 and the reference-free entropy method in stage 2. The adapted L1-SPIRiT method was formulated within the 3D k-space framework. Its efficacy was examined by acquiring two dMRI data sets at 1.05-mm isotropic resolutions with a total acceleration of 6 or 9 (i.e., 2-fold or 3-fold slice and 3-fold in-plane acceleration). Diffusion analysis was performed to derive DTI metrics and estimate fiber orientation distribution functions (fODFs). The results were compared with those of 3D k-space GRAPPA using only ACS, all in reference to 3D k-space GRAPPA using both ACS and SBref (serving as a reference).

RESULTS

The proposed ghost correction eliminated artifacts more robustly than conventional approaches. Our adapted L1-SPIRiT method outperformed 3D k-space GRAPPA when using only ACS, improving image quality to what was achievable with 3D k-space GRAPPA using both ACS and SBref scans. The improvement in image quality further resulted in an improvement in estimation performances for DTI and fODFs.

CONCLUSION

The combination of our new ghost correction and adapted L1-SPIRiT method can reliably reconstruct 7T highly accelerated whole-brain dMRI without the need of SBref scans, increasing acquisition efficiency and reducing motion sensitivity.

Feasibility of Simultaneous Multislice Acceleration Technique in Readout-Segmented Echo-Planar Diffusion-Weighted Imaging for Assessing Rectal Cancer

BACKGROUND

Readout-segmented echo-planar imaging (rs-EPI) with simultaneous multislice (SMS) technology has been successfully applied to tumor research in many organs, but no feasibility study in rectal cancer has been reported, and the optimal acceleration of SMS with rs-EPI in rectal cancer has not been well determined yet.

OBJECTIVE

To investigate the feasibility of SMS rs-EPI of rectal cancer with different acceleration factors (AF_s) and its influence on image quality, acquisition time and apparent diffusion coefficients (ADC_s) in comparison to conventional sequences.

METHODS

All patients underwent rs-EPI and SMS rs-EPI with AF_s of 2 and 3 (2 × SMS rs-EPI and 3 × SMS rs-EPI, respectively) using a 3T scanner. Acquisition times of the three rs-EPI sequences were measured. Image qualitative parameters (5-point Likert scale), signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), geometric distortion, and apparent diffusion coefficient (ADC) values of the three sequences were compared.

RESULTS

A total of eighty-three patients were enrolled in our study. rs-EPI and 2 × SMS rs-EPI offered equivalently high overall image quality with a scan time reduction to nearly half (rs-EPI: 137 s, 2 × SM rs-EPI: 60 s). 3 × SMS rs-EPI showed significantly poorer image quality (p < 0.05). ADC values were significantly lower in 3 × SMS rs-EPI compared to rs-EPI in rectal tumors and normal tissue (tumor tissue: rs-EPI 1.19 ± 0.21 × 10^-3 mm²/s, 3 × SMS rs-EPI 1.10 ± 0.26 × 10^-3 mm²/s, p < 0.001; normal tissue: rs-EPI 1.68 ± 0.13 × 10^-3 mm²/s, 3 × SMS rs-EPI 1.54 ± 0.20 × 10^-3 mm²/s, p < 0.001).

CONCLUSIONS

SMS rs-EPI using an AF of 2 is feasible for rectal MRI resulting in substantial reductions in acquisition time while maintaining diagnostic image quality and similar ADC values to those of rs-EPI when the slice distance and number of shots are the same among three rs-EPI sequences.

A hybrid 2D/4D-MRI methodology using simultaneous multislice imaging for radiotherapy guidance

Purpose

Respiratory motion management is important in abdominothoracic radiotherapy. Fast imaging of the tumor can facilitate multileaf collimator (MLC) tracking that allows for smaller treatment margins, while repeatedly imaging the full field‐of‐view is necessary for 4D dose accumulation. This study introduces a hybrid 2D/4D‐MRI methodology that can be used for simultaneous MLC tracking and dose accumulation on a 1.5 T Unity MR‐linac (Elekta AB, Stockholm, Sweden).

Methods

We developed a hybrid 2D/4D‐MRI methodology that uses a simultaneous multislice (SMS) accelerated MRI sequence, which acquires two coronal slices simultaneously and repeatedly cycles through slice positions over the image volume. As a result, the fast 2D imaging can be used prospectively for MLC tracking and the SMS slices can be sorted retrospectively into respiratory‐correlated 4D‐MRIs for dose accumulation. Data were acquired in five healthy volunteers with an SMS‐bTFE and SMS‐TSE MRI sequence. For each sequence, a prebeam dataset and a beam‐on dataset were acquired simulating the two phases of MR‐linac treatments. Prebeam data were used to generate a 4D‐based motion model and a reference mid‐position volume, while beam‐on data were used for real‐time motion extraction and reconstruction of beam‐on 4D‐MRIs. In addition, an in‐silico computational phantom was used for validation of the hybrid 2D/4D‐MRI methodology. MLC tracking experiments were performed with the developed methodology, for which real‐time SMS data reconstruction was enabled on the scanner. A 15‐beam 8× 7.5 Gy intensity‐modulated radiotherapy plan for lung stereotactic body radiotherapy with isotropic 3 mm GTV‐to‐PTV margins was created. Dosimetry experiments were performed using a 4D motion phantom. The latency between target motion and updating the radiation beam was determined and compensated. Local gamma analyses were performed to quantify dose differences compared to a static reference delivery, and dose area histograms (DAHs) were used to quantify the GTV and PTV coverage.

Results

In‐vivo data acquisition and MLC tracking experiments were successfully performed with the developed hybrid 2D/4D‐MRI methodology. Real‐time liver–lung interface motion estimation had a Pearson's correlation of 0.996 (in‐vivo) and 0.998 (in‐silico). A median (5th–95th percentile) error of 0.0 (−0.9 to 0.7) mm and 0.0 (−0.2 to 0.2) mm was found for real‐time motion estimation for in‐vivo and in‐silico, respectively. Target motion prediction beyond the liver–lung interface had a median root mean square error of 1.6 mm (in‐vivo) and 0.5 mm (in‐silico). Beam‐on 4D MRI reconstruction required a median amount of data equal to an acquisition time of 2:21–3:17 min, which was 20% less data compared to the prebeam‐derived 4D‐MRI. System latency was reduced from 501 ± 12 ms to −1 ± 3 ms (SMS‐TSE) and from 398 ± 10 ms to −10 ± 4 ms (SMS‐bTFE) by a linear regression prediction filter. The local gamma analysis agreed within −3.8% to 3.3% (SMS‐bTFE) and −5.3% to 10% (SMS‐TSE) with a reference MRI sequence. The DAHs revealed a relative D98% GTV coverage between 97% and 100% (SMS‐bTFE) and 100% and 101% (SMS‐TSE) compared to the static reference.

Conclusions

The presented 2D/4D‐MRI methodology demonstrated the potential for accurately extracting real‐time motion for MLC tracking in abdominothoracic radiotherapy, while simultaneously reconstructing contiguous respiratory‐correlated 4D‐MRIs for dose accumulation.

A Slice-Low-Rank Plus Sparse (slice-L + S) Reconstruction Method for k-t Undersampled Multiband First-Pass Myocardial Perfusion MRI

Purpose

The synergistic use of k‐t undersampling and multiband (MB) imaging has the potential to provide extended slice coverage and high spatial resolution for first‐pass perfusion MRI. The low‐rank plus sparse (L + S) model has shown excellent performance for accelerating single‐band (SB) perfusion MRI.

Methods

A MB data consistency method employing ESPIRiT maps and through‐plane coil information was developed. This data consistency method was combined with the temporal L + S constraint to form the slice‐L + S method. Slice‐L + S was compared to SB L + S and the sequential operations of split slice‐GRAPPA and SB L + S (seq‐SG‐L + S) using synthetic data formed from multislice SB images. Prospectively k‐t undersampled MB data were also acquired and reconstructed using seq‐SG‐L + S and slice‐L + S.

Results

Using synthetic data with total acceleration rates of 6–12, slice‐L + S outperformed SB L + S and seq‐SG‐L + S (N = 7 subjects) with respect to normalized RMSE and the structural similarity index (P < 0.05 for both). For the specific case with MB factor = 3 and rate 3 undersampling, or for SB imaging with rate 9 undersampling (N = 7 subjects), the normalized RMSE values were 0.037 ± 0.007, 0.042 ± 0.005, and 0.031 ± 0.004; and the structural similarity index values were 0.88 ± 0.03, 0.85 ± 0.03, and 0.89 ± 0.02 for SB L + S, seq‐SG‐L + S, and slice‐L + S, respectively (P < 0.05 for both). For prospectively undersampled MB data, slice‐L + S provided better image quality than seq‐SG‐L + S for rate 6 (N = 7) and rate 9 acceleration (N = 7) as scored by blinded experts.

Conclusion

Slice‐L + S outperformed SB‐L + S and seq‐SG‐L + S and provides 9 slice coverage of the left ventricle with a spatial resolution of 1.5 mm × 1.5 mm with good image quality.

Model-based dynamic off-resonance correction for improved accelerated fMRI in awake behaving nonhuman primates

Purpose

To estimate dynamic off‐resonance due to vigorous body motion in accelerated fMRI of awake behaving nonhuman primates (NHPs) using the echo‐planar imaging reference navigator, in order to attenuate the effects of time‐varying off‐resonance on the reconstruction.

Methods

In NHP fMRI, the animal’s head is usually head‐posted, and the dynamic off‐resonance is mainly caused by motion in body parts that are distant from the brain and have low spatial frequency. Hence, off‐resonance at each frame can be approximated as a spatially linear perturbation of the off‐resonance at a reference frame, and is manifested as a relative linear shift in k‐space. Using GRAPPA operators, we estimated these shifts by comparing the navigator at each time frame with that at the reference frame. Estimated shifts were then used to correct the data at each frame. The proposed method was evaluated in phantom scans, simulations, and in vivo data.

Results

The proposed method is shown to successfully estimate spatially low‐order dynamic off‐resonance perturbations, including induced linear off‐resonance perturbations in phantoms, and is able to correct retrospectively corrupted data in simulations. Finally, it is shown to reduce ghosting artifacts and geometric distortions by up to 20% in simultaneous multislice in vivo acquisitions in awake‐behaving NHPs.

Conclusion

A method is proposed that does not need sequence modification or extra acquisitions and makes accelerated awake behaving NHP imaging more robust and reliable, reducing the gap between what is possible with NHP protocols and state‐of‐the‐art human imaging.

Leveraging redundancy in simultaneous multislice acquisitions to improve spike detection

PURPOSE

To improve the performance of low-level spike noise artifact detection for daily quality assurance protocols by taking advantage of redundancy in simultaneous multislice (SMS) acquisitions.

METHODS

Magnitude images were transformed into pseudo k-space images. Time series at each pseudo k-space point were detrended. A slice was determined to contain spiking artifact if it exceeded an intensity threshold and if all simultaneously acquired slices contained outliers.

RESULTS

A total of 401 112 slices were inspected. Of these, 42 showed a spike artifact, based on visual inspection of image data and k-space data. With an intensity threshold of 4.6 SDs over time for each pseudo k-space point, all slices containing artifact were correctly flagged, and only 30 slices were incorrectly flagged when using the SMS criterion. Without the SMS criterion, 12 908 slices were incorrectly flagged as containing artifact. Without the SMS criterion, sensitivity to artifact would have to be sacrificed to substantially reduce the number of incorrectly flagged slices.

CONCLUSION

This study demonstrates that the SMS criterion reduced the number of outliers reported to a manageable level while accurately identifying low-level spike artifacts. Successfully identifying low-level spikes allows early detection of hardware problems that can be fixed before the problem becomes debilitating and corrupts data. As part of a daily quality assurance protocol, the method prevents the need to retrospectively carry out time-intensive despiking and reanalysis of data.

Simultaneous Multislice Brain MRI T1 Mapping with Improved Low-Rank Modeling

To accelerate data acquisition speed in magnetic resonance imaging (MRI), multiple slices are simultaneously acquired using multiband pulses. Simultaneous multislice (SMS) imaging typically unfolds slice aliasing from the acquired collapsed slices. In this study, we extended the SMS framework to accelerated MR parameter quantification such as T1 mapping. Assuming that the slice-specific null space and signal subspace are invariant along the parameter dimension, we formulated the SMS framework as a constrained optimization problem under a joint reconstruction framework such that the noise and signal subspaces are used for slice separation and recovery, respectively. The proposed method was validated on 3T MR human brain scans. We successfully demonstrated that the proposed method outperforms competing methods in suppressing aliasing artifacts and noise at high SMS accelerations, thus leading to accurate T1 maps.

Combined simultaneous multislice bSSFP and compressed sensing for first-pass myocardial perfusion at 1.5 T with high spatial resolution and coverage

PURPOSE

To implement and evaluate a pseudorandom undersampling scheme for combined simultaneous multislice (SMS) balanced SSFP (bSSFP) and compressed-sensing (CS) reconstruction to enable myocardial perfusion imaging with high spatial resolution and coverage at 1.5 T.

METHODS

A prospective pseudorandom undersampling scheme that is compatible with SMS-bSSFP phase-cycling requirements and CS was developed. The SMS-bSSFP CS with pseudorandom and linear undersampling schemes were compared in a phantom. A high-resolution (1.4 × 1.4 mm2 ) six-slice SMS-bSSFP CS perfusion sequence was compared with a conventional (1.9 × 1.9 mm2 ) three-slice sequence in 10 patients. Qualitative assessment of image quality, perceived SNR, and number of diagnostic segments and quantitative measurements of sharpness, upslope index, and contrast ratio were performed.

RESULTS

In phantom experiments, pseudorandom undersampling resulted in residual artifact (RMS error) reduction by a factor of 7 compared with linear undersampling. In vivo, the proposed sequence demonstrated higher perceived SNR (2.9 ± 0.3 vs. 2.2 ± 0.6, P = .04), improved sharpness (0.35 ± 0.03 vs. 0.32 ± 0.05, P = .01), and a higher number of diagnostic segments (100% vs. 94%, P = .03) compared with the conventional sequence. There were no significant differences between the sequences in terms of image quality (2.5 ± 0.4 vs. 2.8 ± 0.2, P = .08), upslope index (0.11 ± 0.02 vs. 0.10 ± 0.01, P = .3), or contrast ratio (3.28 ± 0.35 vs. 3.36 ± 0.43, P = .7).

CONCLUSION

A pseudorandom k-space undersampling compatible with SMS-bSSFP and CS reconstruction has been developed and enables cardiac MR perfusion imaging with increased spatial resolution and myocardial coverage, increased number of diagnostic segments and perceived SNR, and no difference in image quality, upslope index, and contrast ratio.

Newly Developed Methods for Reducing Motion Artifacts in Pediatric Abdominal MRI: Tips and Pearls

OBJECTIVE. The purpose of this article is to review established and emerging methods for reducing motion artifacts in pediatric abdominal MRI. CONCLUSION. Clearly understanding the strengths and limitations of motion reduction methods can enable practitioners of pediatric abdominal MRI to select and combine the appropriate techniques and potentially reduce the need for sedation and anesthesia.

Generalized simultaneous multi-orientation 2D imaging

PURPOSE

Flexibility in slice prescription is critical for precise motion monitoring during MR-guided therapies. Adding more slices to improve spatial coverage during rapid 2D cine imaging often hampers temporal resolution. This work describes a framework to simultaneously acquire multiple arbitrarily oriented slices which share a common frequency encoding axis. This framework allows for higher frame rates for a given number of slices compared to conventional interleaved-slice multi-orientation cine imaging.

THEORY AND METHODS

A framework to calculate zeroth gradient moments to be played out between sequentially excited slices with multiple orientations is described here. Experiments were performed in phantom, and in vivo in the head/neck and abdomen of patients.

RESULTS

Images arbitrarily rotated relative to one another were successfully obtained in phantom and in vivo. Simultaneous multi-orientation (SMO) images were also acquired with additional in-plane acceleration to demonstrate the capability of this method to rapidly image objects moving with physiological motion.

CONCLUSIONS

The technical feasibility of the generalized SMO imaging framework was tested in this study. It shows promise for continued development for motion monitoring during MR-guided therapies.

Simultaneous multislice imaging of the heart using multiband balanced SSFP with blipped-CAIPI

Purpose

In this work, we explore the use of multiband (MB) balanced steady‐state free precession (bSSFP) with blipped‐controlled aliasing in parallel imaging (CAIPI), which avoids the issues of altered frequency response associated with RF phase cycling, and show its application to accelerating cardiac cine imaging.

Methods

Blipped and RF‐cycled CAIPI were implemented into a retrospective‐gated segmented cine multiband bSSFP sequence. The 2 methods were compared at 3T using MB2 to demonstrate the effect on frequency response. Further data (4 subjects) were acquired at both 1.5T and 3T collecting 12‐slice short axis stacks using blipped‐CAIPI with MB acceleration factors of 1–4. The impact on SNR and contrast was evaluated along with g‐factors at different accelerations.

Results

Data acquired with blipped‐CAIPI multiband bSSFP up to factor 4 yielded functional cine data with good SNR and contrast, while reliably keeping dark‐band artefacts clear of the heart at 1.5T. SAR limits the maximum MB acceleration, particularly at 3T, where minimum TR increase is problematic and leakage artefacts are more prevalent. Mean g‐factors across the heart were measured at 1.00, 1.06, and 1.12 for MB2–MB4, whereas blood‐pool SNR measures (end‐diastole) decreased by 11.8, 21.5, and 36.9%; ultimately LV‐myocardium CNR remained sufficient at 1.5T with values ranging: 15.6, 13.4, 11.9, and 9.6 (MB1–MB4).

Conclusion

Blipped‐CAIPI multiband bSSFP can be used in cardiovascular applications without affecting the frequency response because of controlled aliasing and can be readily incorporated into segmented cine acquisitions without adding any additional constraints because of phase cycling requirements. The method was used to collect full ventricular coverage within a single breath‐hold.

Non-Cartesian slice-GRAPPA and slice-SPIRiT reconstruction methods for multiband spiral cardiac MRI

PURPOSE

Spiral MRI has advantages for cardiac imaging, and multiband (MB) spiral MRI of the heart shows promise. However, current reconstruction methods for MB spiral imaging have limitations. We sought to develop improved reconstruction methods for MB spiral cardiac MRI.

METHODS

Two reconstruction methods were developed. The first is non-Cartesian slice-GRAPPA (NCSG), which uses phase demodulation and gridding operations before application of a Cartesian slice-separating kernel. The second method, slice-SPIRiT, formulates the reconstruction as a minimization problem that enforces in-plane coil consistency and consistency with the acquired MB data, and uses through-plane coil sensitivity information in the iterative solution. These methods were compared with conjugate-gradient SENSE in phantoms and volunteers. Temporal alternation of CAIPIRINHA (controlled aliasing in parallel imaging results in higher acceleration) phase and the use of a temporal filter were also investigated.

RESULTS

Phantom experiments with 3 simultaneous slices (MB = 3) showed that mean artifact power was highest for conjugate-gradient SENSE, lower for NCSG, and lowest for slice-SPIRiT. For volunteer cine imaging (MB = 3, N = 5), the artifact power was 0.182 ± 0.037, 0.148 ± 0.036, and 0.139 ± 0.034 for conjugate-gradient SENSE, NCSG, and slice-SPIRiT, respectively (P < .05, analysis of variance). Temporal alternation of CAIPIRINHA reduced artifacts for both NCSG and slice-SPIRiT.

CONCLUSION

The NCSG and slice-SPIRiT methods provide more accurate reconstructions for MB spiral cine imaging compared with conjugate-gradient SENSE. These methods hold promise for non-Cartesian MB imaging.

Accelerated whole-brain perfusion imaging using a simultaneous multislice spin-echo and gradient-echo sequence with joint virtual coil reconstruction

PURPOSE

Dynamic susceptibility contrast imaging requires high temporal sampling, which poses limits on achievable spatial coverage and resolution. Additionally, more encoding-intensive multi-echo acquisitions for quantitative imaging are desired to mitigate contrast leakage effects, which further limits spatial encoding. We present an accelerated sequence that provides whole-brain coverage at an improved spatio-temporal resolution, to allow for dynamic quantitative R2 and R2 * mapping during contrast-enhanced imaging.

METHODS

A multi-echo spin and gradient-echo sequence was implemented with simultaneous multislice acquisition. Complementary k-space sampling between repetitions and joint virtual coil reconstruction were used along with a dynamic phase-matching technique to achieve high-quality reconstruction at 9-fold acceleration, which enabled 2 × 2 × 5 mm whole-brain imaging at TR of 1.5 to 1.7 seconds. The multi-echo images from this sequence were fit to achieve quantitative R2 and R2 * maps for each repetition, and subsequently used to find perfusion measures including cerebral blood flow and cerebral blood volume.

RESULTS

Images reconstructed using joint virtual coil show improved image quality and g-factor compared with conventional reconstruction methods, resulting in improved quantitative maps with a 9-fold acceleration factor and whole-brain coverage during the dynamic perfusion acquisition.

CONCLUSION

The method presented shows the advantage of using a joint virtual coil-GRAPPA reconstruction to allow for high acceleration factors while maintaining reliable image quality for quantitative perfusion mapping, with the potential to improve tumor diagnostics and monitoring.

Brain imaging with improved acceleration and SNR at 7 Tesla obtained with 64-channel receive array

PURPOSE

Despite the clear synergy between high channel counts in a receive array and magnetic fields ≥ 7 Tesla, to date such systems have been restricted to a maximum of 32 channels. Here, we examine SNR gains at 7 Tesla in unaccelerated and accelerated images with a 64-receive channel (64Rx) RF coil.

METHODS

A 64Rx coil was built using circular loops tiled in 2 separable sections of a close-fitting form; custom designed preamplifier boards were integrated into each coil element. A 16-channel transmitter arranged in 2 rows along the z-axis was employed. The performance of the 64Rx array was experimentally compared to that of an industry-standard 32-channel receive (32Rx) array for SNR in unaccelerated images and for noise amplification under parallel imaging.

RESULTS

SNR gains were observed in the periphery but not in the center of the brain in unaccelerated imaging compared to the 32Rx coil. With either 1D or 2D undersampling of k-space, or with slice acceleration together with 1D undersampling of k-space, significant reductions in g-factor noise were observed throughout the brain, yielding effective gains in SNR in the entire brain compared to the 32Rx coil. Task-based FMRI data with 12-fold 2D (slice and phase-encode) acceleration yielded excellent quality functional maps with the 64Rx coil but was significantly beyond the capabilities of the 32Rx coil.

CONCLUSION

The results confirm the expectations from modeling studies and demonstrate that whole-brain studies with up to 16-fold, 2D acceleration would be feasible with the 64Rx coil.

Whole-heart spiral simultaneous multi-slice first-pass myocardial perfusion imaging

PURPOSE

To develop and evaluate a simultaneous multislice (SMS) spiral perfusion pulse sequence with whole-heart coverage.

METHODS

An orthogonal set of phase cycling angles following a Hadamard pattern was incorporated into a golden-angle (GA) variable density spiral perfusion sequence to perform SMS imaging at different multiband (MB) factors. Images were reconstructed using an SMS extension of L1-SPIRiT that we have termed SMS-L1-SPIRiT. The proposed sequence was evaluated in 40 subjects (10 each for MB factors of 1, 2, 3, and 4). Images were blindly graded by 2 cardiologists on a 5-point scale (5, excellent). To quantitatively evaluate the reconstruction performance against images acquired without SMS, the MB =1 data were used to retrospectively simulate data acquired at MB factors of 2 to 4.

RESULTS

Analysis of the SMS point-spread function for the desired slice showed that the proposed sampling strategy significantly canceled the main-lobe energy of the other slices and has low side-lobe energy resulting in an incoherent temporal aliasing pattern when rotated by the GA. Retrospective experiments demonstrated the SMS-L1-SPIRiT method removed aliasing from the interfering slices and showed excellent agreement with the ground-truth MB =1 images. Clinical evaluation demonstrated high-quality perfusion images with average image-quality scores of 4.3 ± 0.5 (MB =2), 4.2 ± 0.5 (MB =3), and 4.4 ± 0.4 (MB =4) with no significant quality difference in image quality between MB factors (P = 0.38).

CONCLUSION

SMS spiral perfusion at MB factors 2, 3, and 4 produces high-quality perfusion images with whole-heart coverage in a clinical setting with high sampling efficiency.

Simultaneous multislice cardiac magnetic resonance fingerprinting using low rank reconstruction

This study introduces a technique for simultaneous multislice (SMS) cardiac magnetic resonance fingerprinting (cMRF), which improves the slice coverage when quantifying myocardial T_1, T₂ , and M₀ . The single-slice cMRF pulse sequence was modified to use multiband (MB) RF pulses for SMS imaging. Different RF phase schedules were used to excite each slice, similar to POMP or CAIPIRINHA, which imparts tissues with a distinguishable and slice-specific magnetization evolution over time. Because of the high net acceleration factor (R = 48 in plane combined with the slice acceleration), images were first reconstructed with a low rank technique before matching data to a dictionary of signal timecourses generated by a Bloch equation simulation. The proposed method was tested in simulations with a numerical relaxation phantom. Phantom and in vivo cardiac scans of 10 healthy volunteers were also performed at 3 T. With single-slice acquisitions, the mean relaxation times obtained using the low rank cMRF reconstruction agree with reference values. The low rank method improves the precision in T₁ and T₂ for both single-slice and SMS cMRF, and it enables the acquisition of maps with fewer artifacts when using SMS cMRF at higher MB factors. With this technique, in vivo cardiac maps were acquired from three slices simultaneously during a breathhold lasting 16 heartbeats. SMS cMRF improves the efficiency and slice coverage of myocardial T₁ and T₂ mapping compared with both single-slice cMRF and conventional cardiac mapping sequences. Thus, this technique is a first step toward whole-heart simultaneous T₁ and T₂ quantification with cMRF.

Optimal repetition time reduction for single subject event-related functional magnetic resonance imaging

PURPOSE

Short TRs are increasingly used for fMRI as fast sequences such as simultaneous multislice excitation become available. These have been associated with apparent sensitivity improvements, although greater temporal autocorrelation at shorter TRs can inflate sensitivity measurements leading to uncertainty regarding the optimal approach.

METHODS

In volunteers (n = 10), the optimal TR was assessed at the single subject level for event-related designs (visual stimulation) with 4 frequencies of presentation at 4 TR values (412-2550 ms). T-values in the visual cortex localized in each individual were obtained and receiver operating characteristics (ROC) analysis was performed by counting voxels within and outside expected task active regions at different thresholds. This analysis was repeated using 4 different autoregressive (AR) models; SPM AR(1) and SPM AR(fast) which globally estimate autocorrelation, and fMRIstat AR(1) and AR(5) that use a local estimate.

RESULTS

The use of modest multiband factors of 2 or 3 with a reduction in TR to 1000 ± 200 ms had greater sensitivity and specificity as shown by higher T-values in visual cortex and ROC analysis. At these TRs, the ROC analysis demonstrated that a local AR model fit improved performance while high order AR models were unnecessary.

CONCLUSIONS

Modest TR reductions (to 1000 ± 200 ms) optimally improved event-related fMRI performance independent of design frequency. Autoregressive models with a local as opposed to global fit performed better, while low order autoregressive models were sufficient at the optimal TR.

Accelerated diffusion-weighted imaging for lymph node assessment in the pelvis applying simultaneous multislice acquisition: A healthy volunteer study

To evaluate the feasibility of accelerated simultaneous multislice diffusion weighted sequences (SMS-DWI) for lymph node detection in the abdominopelvic region. Sequences were evaluated regarding the number and depiction of lymph nodes detected with SMS-DWI compared with conventional diffusion weighted sequences, the most suitable SMS- acceleration factor, signal-to-noise ratio (SNR), and the overall acquisition time (TA).Eight healthy volunteers (4 men, 4 women; age range 21-39 years; median age 25 years) were examined in the pelvic region at 3T using a conventional DWI sequence and a SMS DWI sequence with different acceleration factors (AF: 2-3). Moreover, a SMS DWI sequence with AF 3 and higher slice resolution was applied. For morphological correlation of the lymph nodes and as a reference standard, an isotropic 3-dimensional T2-weighted fast-spin-echo sequence with high sampling efficiency (SPACE) was acquired. Two radiologists reviewed each DWI sequence and assessed the number of lymph nodes and the overall image quality. For each DWI sequence, SNR, SNR efficiency per time, contrast to noise (CNR), and ADC values were calculated. Values were statistically compared using a Wilcoxon test (P < .05).Overall, scan time of SMS-DWI with AF2 (AF3) decreased by 46.9% (57.2%) with respect to the conventional DWI. Compared with the SPACE sequence, the detection rate was 89.6% for conventional DWI, 69.4% for SMS-DWI with AF2, and 59.9% for SMS-DWI with AF3. The highly resolved SMS-DWI with AF3 leads to a scan time reduction of 46.9% and detection rate of 83.0%. SNR and CNR were lower in the accelerated sequences (up to 51.0%, P < .001) as compared with the conventional DWI. SNR efficiency decreased to 19.3% for AF2 and to 31.3% for AF3. In the highly resolved dataset, an SNR efficiency reduction of 51.2% was found.This study showed that lymph node detection in the abdominopelvic region with accelerated SMS-DWI sequences is feasible whereby an AF of 2 represents the best compromise between image quality, SNR, CNR, TA, and detection rate.

Disrupted Cerebrocerebellar Intrinsic Functional Connectivity in Young Adults with High-Functioning Autism Spectrum Disorder: A Data-Driven, Whole-Brain, High-Temporal Resolution Functional Magnetic Resonance Imaging Study

This study examines the resting-state functional-connectivity (RsFc) in young adults with high-functioning autism spectrum disorder (HF-ASD) using state-of-the-art fMRI data acquisition and analysis techniques. High temporal resolution fMRI using simultaneous multi-slice acquisition aided unbiased whole-brain connectome-wide multivariate pattern analysis (MVPA) techniques for assessing RsFc. MVPA revealed two clusters (Crus I/II and lobule IX) of abnormal connectivity in the cerebellum that are consistent with the notion of a triple representation of nonmotor processing in the cerebellum. Whole-brain seed-based RsFc analyses informed by these clusters showed significant under connectivity between the cerebellar and social, emotional, and language brain regions in the HF-ASD group compared to healthy controls. The results we report are coherent with existing structural, functional, and RsFc literature in autism, extend previous literature reporting cerebellar abnormalities in the neuropathology of autism, and highlight the cerebellum as a potential target for therapeutic, diagnostic, predictive, and prognostic developments in HF-ASD. The description of functional connectivity abnormalities reported in this study using whole-brain, data-driven analyses has the potential to crucially advance the development of ASD biomarkers, targets for therapeutic interventions, and neural predictors for measuring treatment response.

Simultaneous multislice acquisition without trajectory modification for hyperpolarized 13 C experiments

Purpose

To investigate the feasibility of performing large FOV hyperpolarized ¹³C metabolic imaging using simultaneous multislice excitation.

Methods

A spectral‐spatial multislice excitation pulse was constructed by cosine modulation and incorporated into a ¹³C spiral imaging sequence. Phantom and in vivo pig experiments were performed to test the feasibility of simultaneous multislice data acquisition and image reconstruction. In vivo cardiac‐gated images of hyperpolarized pyruvate, bicarbonate, and lactate were obtained at 1 × 1 × 1 cm³ resolution over a 48 × 48 × 24 cm³ FOV with 2‐fold acceleration in the slice direction. Sensitivity encoding was used for image reconstruction with both autocalibrated and numerically calculated coil sensitivities.

Results

Simultaneous multislice images obtained with 2‐fold acceleration were comparable to reference unaccelerated images. Retained SNR figures greater than 80% were achieved over the part of the image containing the heart.

Conclusion

This method is anticipated to enable large FOV imaging studies using hyperpolarized ¹³C substrates, with an aim toward whole‐body exams that have to date been out of reach.

Simultaneous Multislice Resting-State Functional Magnetic Resonance Imaging at 3 Tesla: Slice-Acceleration-Related Biases in Physiological Effects

Simultaneous multislice echo-planar imaging (SMS-EPI) can enhance the spatiotemporal resolution of resting-state functional MRI (rs-fMRI) by encoding and simultaneously imaging "groups" of slices. However, phenomena, including respiration, cardiac pulsatility, respiration volume per time (RVT), and cardiac rate variation (CRV), referred to as "physiological processes," impact SMS-EPI rs-fMRI in a manner that is yet to be well characterized. In particular, physiological noise may incur aliasing and introduce spurious signals from one slice into another within the "slice group" in rs-fMRI data, resulting in a deleterious effect on resting-state functional connectivity MRI (rs-fcMRI) maps. In the present work, we aimed to quantitatively compare the effects of physiological noise on regular EPI and SMS-EPI in terms of rs-fMRI data and resulting functional connectivity measurements. We compare SMS-EPI and regular EPI data acquired from 11 healthy young adults with matching parameters. The physiological noise characteristics were compared between the two data sets through different combinations of physiological regression steps. We observed that the physiological noise characteristics differed between SMS-EPI and regular EPI, with cardiac pulsatility contributing more to noise in regular EPI data but low-frequency heart rate variability contributing more to SMS-EPI. In addition, a significant slice-group bias was observed in the functional connectivity density maps derived from SMS-EPI data. We conclude that making appropriate corrections for physiological noise is likely more important for SMS-EPI than for regular EPI acquisitions.

High-resolution in vivo diffusion imaging of the human brain with generalized slice dithered enhanced resolution: Simultaneous multislice (gSlider-SMS)

PURPOSE

To develop an efficient acquisition for high-resolution diffusion imaging and allow in vivo whole-brain acquisitions at 600- to 700-μm isotropic resolution.

METHODS

We combine blipped-controlled aliasing in parallel imaging simultaneous multislice (SMS) with a novel slab radiofrequency (RF) encoding gSlider (generalized slice-dithered enhanced resolution) to form a signal-to-noise ratio-efficient volumetric simultaneous multislab acquisition. Here, multiple thin slabs are acquired simultaneously with controlled aliasing, and unaliased with parallel imaging. To achieve high resolution in the slice direction, the slab is volumetrically encoded using RF encoding with a scheme similar to Hadamard encoding. However, with gSlider, the RF-encoding bases are specifically designed to be highly independent and provide high image signal-to-noise ratio in each slab acquisition to enable self-navigation of the diffusion's phase corruption. Finally, the method is combined with zoomed imaging (while retaining whole-brain coverage) to facilitate low-distortion single-shot in-plane encoding with echo-planar imaging at high resolution.

RESULTS

A 10-slices-per-shot gSlider-SMS acquisition was used to acquire whole-brain data at 660 and 760 μm isotropic resolution with b-values of 1500 and 1800 s/mm2 , respectively. Data were acquired on the Connectome 3 Tesla scanner with 64-channel head coil. High-quality data with excellent contrast were achieved at these resolutions, which enable the visualization of fine-scale structures.

CONCLUSIONS

The gSlider-SMS approach provides a new, efficient way to acquire high-resolution diffusion data. Magn Reson Med 79:141-151, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Optimized amplitude modulated multiband RF pulse design

PURPOSE

Multiband pulses are characterized by highly temporally modulated waveforms. Rapid phase or frequency modulation can be extremely demanding on the performance of radiofrequency (RF) pulse generation, which can lead to errors that can be avoided if pulses are restricted to amplitude modulation (AM) only. In this work, three existing multiband pulse design techniques are modified to produce AM waveforms.

THEORY AND METHODS

Multiband refocusing pulses were designed using phase-optimization, time-shifting, and root-flipping. Each technique was constrained to produce AM pulses by exploiting conjugate symmetry in their respective frequency domain representations. Pulses were designed using the AM and unconstrained techniques for a range of multiband factors (ie, number of simultaneously excited slices), time-bandwidth products, and slice separations. Performance was compared by examining the resulting effective pulse durations. Phantom and in vivo experiments were conducted for validation.

RESULTS

Acquired data confirmed that AM pulses can produce precise results when unconstrained designs may produce artifacts. The average duration of AM pulses is longer than the unconstrained versions. Averaged across a range of parameters, the duration cost for AM pulses was 26, 38, and 20% for phase-optimizing, time-shifting and root-flipping, respectively.

CONCLUSIONS

Amplitude modulation multiband pulses can be produced for a relatively small increase in pulse duration. Magn Reson Med 78:2185-2193, 2017. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Simultaneous Multislice-Based 5-Minute Lumbar Spine MRI Protocol: Initial Experience in a Clinical Setting

BACKGROUND AND PURPOSE

Previous studies have used parallel imaging (PI) techniques to decrease spine magnetic resonance imaging (MRI) protocol acquisition times. Recently developed MRI sequences allow even faster acquisitions. Our purpose was to develop a lumbar spine MRI protocol using PI with GRAPPA (generalized autocalibrating partially parallel acquisition) and a simultaneous multislice (SMS)-based sequence and to evaluate its diagnostic performance compared to a standard lumbar spine MRI protocol.

METHODS

Ten patients were scanned in a 3-Tesla scanner (MAGNETOM Skyra, Siemens Healthcare). Each patient was imaged using a standard protocol and an optimized fast protocol acquiring the same contrasts and planes. The fast protocol included sagittal T1, T2, and fat suppressed T2 sequences accelerated with GRAPPA and an SMS-based axial T2-weighted sequence using a high-density spine coil (Siemens MR, 30 channel spine). Two blinded neuroradiologists independently assessed image quality and diagnostic accuracy for clinically relevant imaging findings.

RESULTS

The fast protocol acquisition time was 5:28 minutes, compared with 16:30 minutes for the standard protocol. Both protocols had a similar performance for definition of anatomical structures, diagnostic quality, and identification of clinically relevant imaging findings. There were more artifacts in the SMS Turbo Spin Echo (P = .014) sequence without compromising diagnostic performance. Artifacts in the remaining non-SMS sequences were similar in both protocols (P > .180). The sensitivity, specificity, and accuracy of the 5-minute protocol were 92.3%, 100.0%, and 99.6%, respectively, for the clinically relevant findings (P = 1.0, interrater agreement .57).

CONCLUSIONS

A 5-minute SMS-based MRI protocol for lumbar spine imaging is feasible and can be achieved without significant impact in the overall diagnostic quality.

Real-time multislice MRI during continuous positive airway pressure reveals upper airway response to pressure change

PURPOSE

To determine if a real-time magnetic resonance imaging (RT-MRI) method during continuous positive airway pressure (CPAP) can be used to measure neuromuscular reflex and/or passive collapsibility of the upper airway in individual obstructive sleep apnea (OSA) subjects.

MATERIALS AND METHODS

We conducted experiments on four adolescents with OSA and three healthy controls, during natural sleep and during wakefulness. Data were acquired on a clinical 3T scanner using simultaneous multislice (SMS) RT-MRI during CPAP. CPAP pressure level was alternated between therapeutic and subtherapeutic levels. Segmented airway area changes in response to rapid CPAP pressure drop and restoration were used to estimate 1) upper airway loop gain (UALG), and 2) anatomical risk factors, including fluctuation of airway area (FAA).

RESULTS

FAA significantly differed between OSA patients (2-4× larger) and healthy controls (Student's t-test, P < 0.05). UALG and FAA measurements indicate that neuromuscular reflex and passive collapsibility varied among the OSA patients, suggesting the presence of different OSA phenotypes. Measurements had high intrasubject reproducibility (intraclass correlation coefficient r > 0.7).

CONCLUSION

SMS RT-MRI during CPAP can reproducibly identify physiological traits and anatomical risk factors that are valuable in the assessment of OSA. This technique can potentially locate the most collapsible airway sites. Both UALG and FAA possess large variation among OSA patients.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2017;46:1400-1408.

Use of pattern recognition for unaliasing simultaneously acquired slices in simultaneous multislice MR fingerprinting

PURPOSE

The purpose of this study is to accelerate an MR fingerprinting (MRF) acquisition by using a simultaneous multislice method.

METHODS

A multiband radiofrequency (RF) pulse was designed to excite two slices with different flip angles and phases. The signals of two slices were driven to be as orthogonal as possible. The mixed and undersampled MRF signal was matched to two dictionaries to retrieve T₁ and T₂ maps of each slice. Quantitative results from the proposed method were validated with the gold-standard spin echo methods in a phantom. T₁ and T₂ maps of in vivo human brain from two simultaneously acquired slices were also compared to the results of fast imaging with steady-state precession based MRF method (MRF-FISP) with a single-band RF excitation.

RESULTS

The phantom results showed that the simultaneous multislice imaging MRF-FISP method quantified the relaxation properties accurately compared to the gold-standard spin echo methods. T₁ and T₂ values of in vivo brain from the proposed method also matched the results from the normal MRF-FISP acquisition.

CONCLUSION

T₁ and T₂ values can be quantified at a multiband acceleration factor of two using our proposed acquisition even in a single-channel receive coil. Further acceleration could be achieved by combining this method with parallel imaging or iterative reconstruction. Magn Reson Med 78:1870-1876, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

SMS-HSL: Simultaneous multislice aliasing separation exploiting hankel subspace learning

PURPOSE

To develop a novel, simultaneous multislice reconstruction method that exploits Hankel subspace learning (SMS-HSL) for aliasing separation in the slice direction.

METHODS

An SMS signal model with the Hankel-structured matrix was proposed. To efficiently suppress interslice leakage artifacts from a signal subspace perspective, a null space was learned from the reference data combined over all slices other than a slice of interest using singular value decomposition. Given the fact that the Hankel-structured matrix is rank-deficient while the magnitude between the reference and its estimate is similar in k-space, the SMS-HSL was reformulated as a constrained optimization problem with both low-rank and magnitude priors. SMS signals were projected onto a slice-specific subspace while undesired signals were eliminated using the null space operator. The simulations and experiments were performed with increasing multiband factors up to 6 using the SMS-HSL and the split slice-GRAPPA.

RESULTS

Compared with the split slice-GRAPPA, the SMS-HSL shows superior performance in suppressing aliasing artifacts and noises at high multiband factors even with: insufficient reference signals, a small number of coils, and a short distance between aliasing voxels.

CONCLUSION

We successfully demonstrated the effectiveness of the SMS-HSL over the split slice-GRAPPA for aliasing separation in the slice direction. Magn Reson Med 78:1392-1404, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

High-resolution gradient-recalled echo imaging at 9.4T using 16-channel parallel transmit simultaneous multislice spokes excitations with slice-by-slice flip angle homogenization

Purpose

In order to fully benefit from the improved signal‐to‐noise and contrast‐to‐noise ratios at 9.4T, the challenges of B1+ inhomogeneity and the long acquisition time of high‐resolution 2D gradient‐recalled echo (GRE) imaging were addressed.

Theory and Methods

Flip angle homogenized excitations were achieved by parallel transmission (pTx) of 3‐spoke pulses, designed by magnitude least‐squares optimization in a slice‐by‐slice fashion; the acquisition time reduction was achieved by simultaneous multislice (SMS) pulses. The slice‐specific spokes complex radiofrequency scaling factors were applied to sinc waveforms on a per‐channel basis and combined with the other pulses in an SMS slice group to form the final SMS‐pTX pulse. Optimal spokes locations were derived from simulations.

Results

Flip angle maps from presaturation TurboFLASH showed improvement of flip angle homogenization with 3‐spoke pulses over CP‐mode excitation (normalized root‐mean‐square error [NRMSE] 0.357) as well as comparable excitation homogeneity across the single‐band (NRMSE 0.119), SMS‐2 (NRMSE 0.137), and SMS‐3 (NRMSE 0.132) 3‐spoke pulses. The application of the 3‐spoke SMS‐3 pulses in a 48‐slice GRE protocol, which has an in‐plane resolution of 0.28 × 0.28 mm, resulted in a 50% reduction of scan duration (total acquisition time 6:52 min including reference scans).

Conclusion

Time‐efficient flip angle homogenized high‐resolution GRE imaging at 9.4T was accomplished by using slice‐specific SMS‐pTx spokes excitations. Magn Reson Med 78:1050–1058, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.

Two-dimensional-NGC-SENSE-GRAPPA for fast, ghosting-robust reconstruction of in-plane and slice-accelerated blipped-CAIPI echo planar imaging

Purpose

Ghosting‐robust reconstruction of blipped‐CAIPI echo planar imaging simultaneous multislice data with low computational load.

Methods

To date, Slice‐GRAPPA, with “odd–even” kernels that improve ghosting performance, has been the framework of choice for such reconstructions due to its predecessor SENSE‐GRAPPA being deemed unsuitable for blipped‐CAIPI data. Modifications to SENSE‐GRAPPA are used to restore CAIPI compatibility and to make it robust against ghosting. Two implementations are tested, one where slices and in‐plane unaliasing are dealt in the same serial manner as in Slice‐GRAPPA [referred to as one‐dimensional (1D)‐NGC‐SENSE‐GRAPPA, where NGC stands for Nyquist Ghost Corrected] and one where both are unaliased in a single step (2D‐NGC‐SENSE‐GRAPPA), which is analytically and experimentally shown to be computationally cheaper.

Results

The 1D‐NGC‐SENSE‐GRAPPA and odd‐even Slice‐GRAPPA perform identically, whereas 2D‐NGC‐SENSE‐GRAPPA shows reduced error propagation, less residual ghosting when reliable reference data were available. When the latter was not the case, error propagation was increased.

Conclusion

Unlike Slice‐GRAPPA, SENSE‐GRAPPA operates fully within the GRAPPA framework, for which improved reconstructions (e.g., iterative, nonlinear) have been developed over the past decade. It could, therefore, bring benefit to the reconstruction of SMS data as an attractive alternative to Slice‐GRAPPA. Magn Reson Med 77:998–1009, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Simultaneous multislice diffusion-weighted MRI of the liver: Analysis of different breathing schemes in comparison to standard sequences

PURPOSE

To systematically evaluate image characteristics of simultaneous-multislice (SMS)-accelerated diffusion-weighted imaging (DWI) of the liver using different breathing schemes in comparison to standard sequences.

MATERIALS AND METHODS

DWI of the liver was performed in 10 healthy volunteers and 12 patients at 1.5T using an SMS-accelerated echo planar imaging sequence performed with respiratory-triggering and free breathing (SMS-RT, SMS-FB). Standard DWI sequences served as reference (STD-RT, STD-FB). Reduction of scan time by SMS-acceleration was measured. Image characteristics of SMS-DWI and STD-DWI with both breathing schemes were analyzed quantitatively (apparent diffusion coefficient [ADC], signal-to-noise ratio [SNR]) and qualitatively (5-point Likert scale, 5 = excellent). Qualitative and quantitative parameters were compared using Friedman test and Dunn-Bonferroni post-hoc method with P-values < 0.05 considered statistically significant.

RESULTS

SMS-DWI provided diagnostic image quality in volunteers and patients both with RT and FB with a reduction of scan time of 70% (0:56 vs. 3:20 min in FB). Overall image quality did not significantly differ between FB and RT acquisition in both STD and SMS sequences (median STD-RT 5.0, STD-FB 4.5, SMS-RT: 4.75; SMS-FB: 4.5; P = 0.294). SNR in the right hepatic lobe was comparable between the four tested sequences. ADC values were significantly lower in SMS-DWI compared to STD-DWI irrespective of the breathing scheme (1.2 ± 0.2 × 10(-3) mm(2) /s vs. 1.0 ± 0.2 × 10(-3) mm(2) /s; P < 0.001).

CONCLUSION

SMS-acceleration provides considerable scan time reduction for hepatic DWI with equivalent image quality compared to the STD technique both using RT and FB. Discrepancies in ADC between STD-DWI and SMS-DWI need to be considered when transferring the SMS technique to clinical routine reading. J. MAGN. RESON. IMAGING 2016;44:865-879.

Coil compression in simultaneous multislice functional MRI with concentric ring slice-GRAPPA and SENSE

PURPOSE

Simultaneous multislice (SMS) imaging is a useful way to accelerate functional magnetic resonance imaging (fMRI). As acceleration becomes more aggressive, an increasingly larger number of receive coils are required to separate the slices, which significantly increases the computational burden. We propose a coil compression method that works with concentric ring non-Cartesian SMS imaging and should work with Cartesian SMS as well. We evaluate the method on fMRI scans of several subjects and compare it to standard coil compression methods.

METHODS

The proposed method uses a slice-separation k-space kernel to simultaneously compress coil data into a set of virtual coils. Five subjects were scanned using both non-SMS fMRI and SMS fMRI with three simultaneous slices. The SMS fMRI scans were processed using the proposed method, along with other conventional methods. Code is available at https://github.com/alcu/sms.

RESULTS

The proposed method maintained functional activation with a fewer number of virtual coils than standard SMS coil compression methods. Compression of non-SMS fMRI maintained activation with a slightly lower number of virtual coils than the proposed method, but does not have the acceleration advantages of SMS fMRI.

CONCLUSION

The proposed method is a practical way to compress and reconstruct concentric ring SMS data and improves the preservation of functional activation over standard coil compression methods in fMRI. Magn Reson Med 76:1196-1209, 2016. © 2015 Wiley Periodicals, Inc.

Design of parallel transmission pulses for simultaneous multislice with explicit control for peak power and local specific absorption rate

PURPOSE

To design parallel transmit (pTx) simultaneous multislice (SMS) spokes pulses with explicit control for peak power and local and global specific absorption rate (SAR).

METHODS

We design SMS pTx least-squares and magnitude least squares spokes pulses while constraining local SAR using the virtual observation points (VOPs) compression of SAR matrices. We evaluate our approach in simulations of a head (7T) and a body (3T) coil with eight channels arranged in two z-rows.

RESULTS

For many of our simulations, control of average power by Tikhonov regularization of the SMS pTx spokes pulse design yielded pulses that violated hardware and SAR safety limits. On the other hand, control of peak power alone yielded pulses that violated local SAR limits. Pulses optimized with control of both local SAR and peak power satisfied all constraints and therefore had the best excitation performance under limited power and SAR constraints. These results extend our previous results for single slice pTx excitations but are more pronounced because of the large power demands and SAR of SMS pulses.

CONCLUSIONS

Explicit control of local SAR and peak power is required to generate optimal SMS pTx excitations satisfying both the system's hardware limits and regulatory safety limits.

Scan time reduction for readout-segmented EPI using simultaneous multislice acceleration: Diffusion-weighted imaging at 3 and 7 Tesla

PURPOSE

Readout-segmented echo-planar imaging (rs-EPI) can provide high quality diffusion data because it is less prone to distortion and blurring artifacts than single-shot echo-planar imaging (ss-EPI), particularly at higher resolution and higher field. Readout segmentation allows shorter echo-spacing and echo train duration, resulting in reduced image distortion and blurring, respectively, in the phase-encoding direction. However, these benefits come at the expense of longer scan times because the segments are acquired in multiple repetitions times (TRs). This study shortened rs-EPI scan times by reducing the TR duration with simultaneous multislice acceleration.

METHODS

The blipped-CAIPI method for slice acceleration with reduced g-factor SNR loss was incorporated into the diffusion-weighted rs-EPI sequence. The rs- and ss-EPI sequences were compared at a range of resolutions at both 3 and 7 Tesla in terms of image fidelity and diffusion postprocessing results.

RESULTS

Slice-accelerated clinically useful trace-weighted images and tractography results are presented. Tractography analysis showed that the reduced artifacts in rs-EPI allowed better discrimination of tracts than ss-EPI.

CONCLUSION

Slice acceleration reduces rs-EPI scan times providing a practical alternative to diffusion-weighted ss-EPI with reduced distortion and high resolution. Magn Reson Med 74:136-149, 2015. © 2014 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.

Interslice leakage artifact reduction technique for simultaneous multislice acquisitions

PURPOSE

Controlled aliasing techniques for simultaneously acquired echo-planar imaging slices have been shown to significantly increase the temporal efficiency for both diffusion-weighted imaging and functional magnetic resonance imaging studies. The "slice-GRAPPA" (SG) method has been widely used to reconstruct such data. We investigate robust optimization techniques for SG to ensure image reconstruction accuracy through a reduction of leakage artifacts.

METHODS

Split SG is proposed as an alternative kernel optimization method. The performance of Split SG is compared to standard SG using data collected on a spherical phantom and in vivo on two subjects at 3 T. Slice-accelerated and nonaccelerated data were collected for a spin-echo diffusion-weighted acquisition. Signal leakage metrics and time-series SNR were used to quantify the performance of the kernel fitting approaches.

RESULTS

The Split SG optimization strategy significantly reduces leakage artifacts for both phantom and in vivo acquisitions. In addition, a significant boost in time-series SNR for in vivo diffusion-weighted acquisitions with in-plane 2× and slice 3× accelerations was observed with the Split SG approach.

CONCLUSION

By minimizing the influence of leakage artifacts during the training of SG kernels, we have significantly improved reconstruction accuracy. Our robust kernel fitting strategy should enable better reconstruction accuracy and higher slice-acceleration across many applications.

Three-dimensional Fourier encoding of simultaneously excited slices: generalized acquisition and reconstruction framework

PURPOSE

Simultaneous multislice (SMS) acquisitions have recently received much attention as a means of increasing single-shot imaging speed. SMS acquisitions combine the advantages of single-shot sampling and acceleration along the slice dimension which was previously limited to three-dimensional (3D) volumetric acquisitions. A two-dimensional description of SMS sampling and reconstruction has become established in the literature. Here, we present a more general 3D Fourier encoding and reconstruction formalism for SMS acquisitions that can easily be applied to non-Cartesian SMS acquisitions.

THEORY AND METHODS

An "SMS 3D" k-space is defined in which the field of view along the slice select direction is equal to the number of excited slices times their separation. In this picture, SMS acceleration can be viewed as an undersampling of SMS 3D k-space that can be freely distributed between the in-plane and slice directions as both are effective phase-encoding directions.

RESULTS

Use of the SMS 3D k-space picture is demonstrated in phantom and in vivo brain acquisitions including data obtained with blipped-controlled aliasing in parallel imaging sampling. SMS sensitivity encoding reconstruction is demonstrated as well as non-Cartesian SMS imaging using blipped spiral trajectories.

CONCLUSIONS

The full framework of reconstruction methods can be applied to SMS acquisitions by employing a 3D k-space approach. The blipped-controlled aliasing in parallel imaging method can be viewed as a special case of undersampling an SMS 3D k-space. The extension of SMS methods to non-Cartesian 3D sampling and reconstruction is straightforward.

Simultaneous multislice multiband parallel radiofrequency excitation with independent slice-specific transmit B1 homogenization

PURPOSE

To develop a new parallel transmit (pTx) pulse design for simultaneous multiband (MB) excitation in order to tackle simultaneously the problems of transmit B1 (B1+) inhomogeneity and total radiofrequency (RF) power, so as to allow for optimal RF excitation when using MB pulses for slice acceleration for high and ultrahigh field MRI.

METHODS

With the proposed approach, each of the bands that are simultaneously excited is subject to a band-specific set of B1 complex shim weights. The method was validated in the human brain at 7T using a 16-channel pTx system and was compared to conventional MB pulses operating in the circularly polarized (CP) mode. Further numerical simulations based on measured B1 maps were conducted.

RESULTS

The new method improved B1+ homogeneity by 60% when keeping the total RF power constant and reduced total RF power by 72% when keeping the excitation fidelity constant, as compared to the conventional CP mode.

CONCLUSION

A new pTx pulse design formalism is introduced targeting slice-specific B1+ homogenization in MB excitation while constraining total RF power. These pulses lead to significantly improved slice-wise B1+ uniformity and/or largely reduced total RF power, as compared to the conventionally employed MB pulses applied in the CP mode.

Slice accelerated diffusion-weighted imaging at ultra-high field strength

PURPOSE

Diffusion magnetic resonance imaging (dMRI) data with very high isotropic resolution can be obtained at 7T. However, for extensive brain coverage, a large number of slices is required, resulting in long acquisition times (TAs). Recording multiple slices simultaneously (SMS) promises to reduce the TA.

METHODS

A combination of zoomed and parallel imaging is used to achieve high isotropic resolution dMRI data with a low level of distortions at 7T. The blipped-CAIPI (controlled aliasing in parallel imaging) approach is used to acquire several slices simultaneously. Due to their high radiofrequency (RF) power deposition and ensuing specific absorption rate (SAR) constraints, the commonly used multiband (MB) RF pulses for SMS imaging are inefficient at 7T and entail long repetition times, counteracting the usefulness of SMS acquisitions. To address this issue, low SAR multislice Power Independent of Number of Slices RF pulses are employed.

RESULTS

In vivo dMRI results with and without SMS acceleration are presented at different isotropic spatial resolutions at ultra high field strength. The datasets are recorded at a high angular resolution to detect fiber crossings.

CONCLUSION

From the results and compared with earlier studies at these resolutions, it can be seen that scan time is significantly reduced, while image quality is preserved.