Comparison of Echo-Planar Imaging and Compressed Sensing in the Estimation of Flow Metrics from Aortic 4D Flow MR Imaging: A Healthy Volunteer Study

PURPOSE

Prolonged scanning of time-resolved 3D phase-contrast MRI (4D flow MRI) limits its routine use in clinical practice. An echo-planar imaging (EPI)-based sequence and compressed sensing can reduce the scan duration. We aimed to determine the impact of EPI for 4D flow MRI on the scan duration, image quality, and quantitative flow metrics.

METHODS

This was a prospective study of 15 healthy volunteers (all male, mean age 33 ± 5 years). Conventional sensitivity encoding (SENSE), EPI with SENSE (EPI), and compressed SENSE (CS) (reduction factors: 6 and 12, respectively) were scanned.Scan duration, qualitative indexes of image quality, and quantitative flow parameters of net flow volume, maximum flow velocity, wall shear stress (WSS), and energy loss (EL) in the ascending aorta were assessed. Two-dimensional phase-contrast cine MRI (2D-PC) was considered the gold standard of net flow volume and maximum flow velocity.

RESULTS

Compared to SENSE, EPI and CS12 shortened scan durations by 71% and 73% (EPI, 4 min 39 sec; CS6, 7 min 29 sec; CS12, 4 min 14 sec; and SENSE, 15 min 51 sec). Visual image quality was significantly better for EPI than for SENSE and CS (P < 0.001). The net flow volumes obtained with SENSE, EPI, and CS12 and those obtained with 2D-PC were correlated well (r = 0.950, 0.871, and 0.850, respectively). However, the maximum velocity obtained with EPI was significantly underestimated (P < 0.010). The average WSS was significantly higher with EPI than with SENSE, CS6, and CS12 (P < 0.001, P = 0.040, and P = 0.012, respectively). The EL was significantly lower with EPI than with CS6 and CS12 (P = 0.002 and P = 0.007, respectively).

CONCLUSION

EPI reduced the scan duration, improved visual image quality, and was associated with more accurate net flow volume than CS. However, the flow velocity, WSS, and EL values obtained with EPI and other sequences may not be directly comparable.

Diagnostic value of synthetic diffusion-weighted imaging on breast magnetic resonance imaging assessment: comparison with conventional diffusion-weighted imaging

PURPOSE

To compare images generated by synthetic diffusion-weighted imaging (sDWI) with those from conventional DWI in terms of their diagnostic performance in detecting breast lesions when performing breast magnetic resonance imaging (MRI).

METHODS

A total of 128 consecutive patients with 135 enhanced lesions who underwent dynamic MRI between 2018 and 2021 were included. The sDWI and DWI signals were compared by three radiologists with at least 10 years of experience in breast radiology.

RESULTS

Of the 82 malignant lesions, 91.5% were hyperintense on sDWI and 73.2% were hyperintense on DWI. Of the 53 benign lesions, 71.7% were isointense on sDWI and 37.7% were isointense on DWI. sDWI provides accurate signal intensity data with statistical significance compared with DWI (P < 0.05). The diagnostic performance of DWI and sDWI to differentiate malignant breast masses from benign masses was as follows: sensitivity 73.1% [95% confidence interval (CI): 62-82], specificity 37.7% (95% CI: 24-52); sensitivity 91.5% (95% CI: 83-96), specificity 71.7% (95% CI: 57-83), respectively. The diagnostic accuracy of DWI and sDWI was 59.2% and 83.7%, respectively. However, when the DWI images were evaluated with apparent diffusion coefficient mapping and compared with the sDWI images, the sensitivity was 92.68% (95% CI: 84-97) and the specificity was 79.25% (95% CI: 65-89) with no statistically significant difference. The inter-reader agreement was almost perfect (P < 0.001).

CONCLUSION

Synthetic DWI is superior to DWI for lesion visibility with no additional acquisition time and should be taken into consideration when conducting breast MRI scans. The evaluation of sDWI in routine MRI reporting will increase diagnostic accuracy.

Measurement of changes in cerebrospinal fluid pulsation after traumatic brain injury using echo-planar imaging-based functional MRI

Traumatic brain injury (TBI) is a major public health concern worldwide, with a high incidence and a significant impact on morbidity and mortality. The alteration of cerebrospinal fluid (CSF) dynamics after TBI is a well-known phenomenon; however, the underlying mechanisms and their implications for cognitive function are not fully understood. In this study, we propose a new approach to studying the alteration of CSF dynamics in TBI patients. Our approach involves using conventional echo-planar imaging-based functional MRI with no additional scan, allowing for simultaneous assessment of functional CSF dynamics and blood oxygen level-dependent-based functional brain activities. We utilized two previously suggested indices of (i) CSFpulse, and (ii) correlation between global brain activity and CSF inflow. Using CSFpulse, we demonstrated a significant decrease in CSF pulsation following TBI (p < 0.05), which was consistent with previous studies. Furthermore, we confirmed that the decrease in CSF pulsation was most prominent in the early months after TBI, which could be explained by ependymal ciliary loss, intracranial pressure increment, or aquaporin-4 dysregulation. We also observed a decreasing trend in the correlation between global brain activity and CSF inflow in TBI patients (p < 0.05). Our findings suggest that the decreased CSF pulsation after TBI could lead to the accumulation of toxic substances in the brain and an adverse effect on brain function. Further longitudinal studies with larger sample sizes, TBI biomarker data, and various demographic information are needed to investigate the association between cognitive decline and CSF dynamics after TBI. Overall, this study sheds light on the potential role of altered CSF dynamics in TBI-induced neurologic symptoms and may contribute to the development of novel therapeutic interventions.

Evaluating Reproducibility of the ADC and Distortion in Diffusion-weighted Imaging (DWI) with Reverse Encoding Distortion Correction (RDC)

PURPOSE

To compare image distortion and reproducibility of quantitative values between reverse encoding distortion correction (RDC) diffusion-weighted imaging (DWI) and conventional DWI techniques in a phantom study and in healthy volunteers.

METHODS

This prospective study was conducted with the approval of our institutional review board. Written informed consent was obtained from each participant. RDC-DWIs were created from images obtained at 3T in three orthogonal directions in a phantom and in 10 participants (mean age, 70.9 years; age range, 63-83 years). Images without distortion correction (noDC-DWI) and those corrected with B0 (B0c-DWI) were also created. To evaluate distortion, coefficients of variation were calculated for each voxel and ROIs were placed at four levels of the brain. To evaluate the reproducibility of apparent diffusion coefficient (ADC) measurements, intra- and inter-scan variability (%CVADC) were calculated from repeated scans of the phantom. Analysis was performed using Wilcoxon signed-rank test with Bonferroni correction, and P < 0.05 was considered statistically significant.

RESULTS

In the phantom, distortion was less in RDC-DWI than in B0c-DWI (P < 0.006), and was less in B0c-DWI than in noDC-DWI (P < 0.006). Intra-scan %CVADC was within 1.30%, and inter-scan %CVADC was within 2.99%. In the volunteers, distortion was less in RDC-DWI than in B0c-DWI in three of four locations (P < 0.006), and was less in B0c-DWI than in noDC-DWI (P < 0.006). At the middle cerebellar peduncle, distortion was less in RDC-DWI than in noDC-DWI (P < 0.006), and was less in noDC-DWI than in B0c-DWI (P < 0.0177).

CONCLUSION

In both the phantom and in volunteers, distortion was the least in RDC-DWI than in B0c-DWI and noDC-DWI.

Estimation of undistorted images in brain echo-planar images with distortions using the conjugate gradient method with anatomical regularization

PURPOSE

Although echo-planar imaging (EPI) is widely used for diffusion magnetic resonance (MR) imaging, EPI images suffer from susceptibility-induced geometric distortions. We herein propose a new estimation method for undistorted EPI images using anatomical T₁ -weighted images (T₁ WIs) based on the physics of MR imaging.

METHODS

Our proposed method estimates the undistorted EPI image in the image domain while estimating the magnetic field inhomogeneity map using the conjugate gradient method with anatomical regularization. Our method synthesizes the distorted image to match the measured EPI image containing geometric distortions by alternately updating the undistorted EPI image and the magnetic field inhomogeneity map. We evaluated our proposed method and compared it with a nonrigid registration-based distortion correction method using simulated data and using real data. In the evaluation of the estimation of the magnetic field inhomogeneity map, we used the normalized root-mean-squared error (NRMSE) between the estimated results and the ground truth. In the evaluation of the estimation of undistorted images, we used mutual information (MI) between the undistorted EPI image and the anatomical T₁ WI.

RESULTS

Using the simulated data, the means and standard deviations of the NRMSE values in the nonrigid registration-based method and proposed method were 1.29 ± 0.63 and 0.64 ± 0.30, respectively. The MI values in the proposed method were larger than those in the nonrigid registration-based method in all evaluated slices. For the real data, the proposed method improved the distortion, and the MI values in the proposed method were larger than those in the nonrigid registration-based method. In the estimation of the magnetic field inhomogeneity map, the NRMSE values in our method were smaller than those in the nonrigid registration-based method.

CONCLUSIONS

We demonstrated that our proposed method can estimate the regions with compressed distortions that are not well represented by the nonrigid registration-based methods. The results suggest that the proposed method could be useful in analyses combining EPI images with T₁ WIs.

Acceleration of Magnetic Resonance Fingerprinting Reconstruction Using Denoising and Self-Attention Pyramidal Convolutional Neural Network

Magnetic resonance fingerprinting (MRF) based on echo-planar imaging (EPI) enables whole-brain imaging to rapidly obtain T1 and T2* relaxation time maps. Reconstructing parametric maps from the MRF scanned baselines by the inner-product method is computationally expensive. We aimed to accelerate the reconstruction of parametric maps for MRF-EPI by using a deep learning model. The proposed approach uses a two-stage model that first eliminates noise and then regresses the parametric maps. Parametric maps obtained by dictionary matching were used as a reference and compared with the prediction results of the two-stage model. MRF-EPI scans were collected from 32 subjects. The signal-to-noise ratio increased significantly after the noise removal by the denoising model. For prediction with scans in the testing dataset, the mean absolute percentage errors between the standard and the final two-stage model were 3.1%, 3.2%, and 1.9% for T1, and 2.6%, 2.3%, and 2.8% for T2* in gray matter, white matter, and lesion locations, respectively. Our proposed two-stage deep learning model can effectively remove noise and accurately reconstruct MRF-EPI parametric maps, increasing the speed of reconstruction and reducing the storage space required by dictionaries.

Recommendations of Choice of Head Coil and Prescan Normalize Filter Depend on Region of Interest and Task

The performance of MRI head coils together with the influence of the prescan normalize filter in different brain regions was evaluated. Functional and structural data were recorded from 26 participants performing motor, auditory, and visual tasks in different conditions: with the 20- and 64-channel Siemens head/neck coil and the prescan normalize filter turned ON or OFF. Data were analyzed with the MRIQC tool to evaluate data quality differences. The functional data were statistically evaluated by comparison of the β estimates and the time-course signal-to-noise ratio (tSNR) in four regions of interest, i.e., the auditory, visual, and motor cortices and the thalamus. The MRIQC tool indicated a better data quality for both functional and structural data with the prescan normalize filter, with an advantage for the 20-channel head coil in functional data and an advantage for the 64-channel head coil in structural measurements. Nevertheless, recommendations for the functional data regarding choice of head coils and prescan normalize filter depend on the brain regions of interest. Higher β estimates and tSNR values occurred in the auditory cortex and thalamus with the prescan normalize filter, whereas the contrary was true for the visual and motor cortices. Due to higher β estimates in the visual cortex in the 64-channel head coil, this head coil is recommended for studies investigating the visual cortex. For most of the research questions, the 20-channel head coil is better suited for functional experiments, with the prescan normalize filter, especially when investigating deep brain areas. For anatomical studies, the 64-channel head coil seemed to be the better choice.

Diffusion-weighted Imaging of the Abdomen during a Single Breath-hold Using Simultaneous-multislice Echo-planar Imaging

PURPOSE

This multi-scanner study aimed to investigate the validity of single breath-hold (BH) diffusion-weighted imaging (DWI) using simultaneous-multislice (SMS) echo-planar imaging in multiple abdominal organs to enable faster acquisition and reliable quantification of apparent diffusion coefficient (ADC).

METHODS

SNR, geometric distortion (GD), and ADC in a phantom; the ADC in the liver, renal cortex, paraspinal muscle, spleen, and pancreas; and the signal intensity ratio of the portal vein-to-muscle (SIR_PV-M) in healthy volunteers were compared between BH- and respiratory-triggered (RT) DWI with b-values of 0 and 800 s/mm² in two different MRI scanners.

RESULTS

The phantom study showed that the SNR of BH-DWI was significantly lower than that of the RT-DWI (P < 0.05 for both scanners), whereas the GD and ADC of BH-DWI did not differ significantly from those of the RT-DWI (P = 0.09-0.60). In the volunteer study, the scan times were 23 seconds for BH-DWI and 184±33 seconds for RT-DWI, respectively. The ADC of the liver in BH-DWI was significantly lower than that in RT-DWI (P < 0.05 for both scanners), whereas there were no significant differences in the ADCs of the renal cortex, paraspinal muscle, spleen, or pancreas between BH-DWI and RT-DWI (P = 0.07-0.86). The SIR_PV-M in BH-DWI was significantly smaller than in RT-DWI (P < 0.05 for both scanners).

CONCLUSION

The proposed method enables the acquisition of abdominal diffusion-weighted images in a single BH.

Comparison between turbo spin-echo and echo planar diffusion-weighted imaging of the female pelvis with 3T MRI

Background

Echo-planar imaging (EPI)-diffusion-weighted imaging (DWI) may take unclear image affected by susceptibility, geometric distortions and chemical shift artifacts.

Purpose

To compare the image quality and usefulness of EPI-DWI and turbo spin echo (TSE)-DWI in female patients who required imaging of the pelvis.

Material and Methods

All 57 patients were examined with a 3.0-T MR scanner. Both TSE- and EPI-DWI were performed with b values of 0 and 1000 s/mm². We compared geometric distortion, the contrast ratio (CR) of the myometrium to the muscle and the apparent diffusion coefficient (ADC) values for the myometrium and lesion. Two radiologists scored the TSE- and EPI-DWI of each patient for qualitative evaluation.

Results

The mean percent distortion was significantly smaller with TSE- than EPI-DWI (p = 0.00). The CR was significantly higher with TSE- than EPI-DWI (p = 0.003). There was a significant difference in the ADC value for the uterus and lesions between the EPI- and TSE-DWI (p < 0.05). Finally, the ADC values of cancer were significantly different from those for the uterus and benign with both the two sequences (p < 0.05). The scores for ghosting artifacts were higher with TSE- than EPI-DWI (p = 0.019). But there were no significant differences between TSE- and EPI-DWI with regard to image contrast and overall image quality.

Conclusion

TSE-DWI on the female pelvis by 3T MRI produces less distortion and higher CR than EPI-DWI, but there is no difference in contrast and image quality.

Robust autocalibrated structured low-rank EPI ghost correction

PURPOSE

We propose and evaluate a new structured low-rank method for echo-planar imaging (EPI) ghost correction called Robust Autocalibrated LORAKS (RAC-LORAKS). The method can be used to suppress EPI ghosts arising from the differences between different readout gradient polarities and/or the differences between different shots. It does not require conventional EPI navigator signals, and is robust to imperfect autocalibration data.

METHODS

Autocalibrated LORAKS is a previous structured low-rank method for EPI ghost correction that uses GRAPPA-type autocalibration data to enable high-quality ghost correction. This method works well when the autocalibration data are pristine, but performance degrades substantially when the autocalibration information is imperfect. RAC-LORAKS generalizes Autocalibrated LORAKS in two ways. First, it does not completely trust the information from autocalibration data, and instead considers the autocalibration and EPI data simultaneously when estimating low-rank matrix structure. Second, it uses complementary information from the autocalibration data to improve EPI reconstruction in a multi-contrast joint reconstruction framework. RAC-LORAKS is evaluated using simulations and in vivo data, including comparisons to state-of-the-art methods.

RESULTS

RAC-LORAKS is demonstrated to have good ghost elimination performance compared to state-of-the-art methods in several complicated EPI acquisition scenarios (including gradient-echo brain imaging, diffusion-encoded brain imaging, and cardiac imaging).

CONCLUSIONS

RAC-LORAKS provides effective suppression of EPI ghosts and is robust to imperfect autocalibration data.

On the signal-to-noise ratio benefit of spiral acquisition in diffusion MRI

PURPOSE

Spiral readouts combine several favorable properties that promise superior net sensitivity for diffusion imaging. The purpose of this study is to verify the signal-to-noise ratio (SNR) benefit of spiral acquisition in comparison with current echo-planar imaging (EPI) schemes.

METHODS

Diffusion-weighted in vivo brain data from three subjects were acquired with a single-shot spiral sequence and several variants of single-shot EPI, including full-Fourier and partial-Fourier readouts as well as different diffusion-encoding schemes. Image reconstruction was based on an expanded signal model including field dynamics obtained by concurrent field monitoring. The effective resolution of each sequence was matched to that of full-Fourier EPI with 1 mm nominal resolution. SNR maps were generated by determining the noise statistics of the raw data and analyzing the propagation of equivalent synthetic noise through image reconstruction. Using the same approach, maps of noise amplification due to parallel imaging (g-factor) were calculated for different acceleration factors.

RESULTS

Relative to full-Fourier EPI at b = 0 s/mm² , spiral acquisition yielded SNR gains of 42-88% and 40-89% in white and gray matter, respectively, depending on the diffusion-encoding scheme. Relative to partial-Fourier EPI, the gains were 36-44% and 34-42%. Spiral g-factor maps exhibited less spatial variation and lower maxima than their EPI counterparts.

CONCLUSION

Spiral readouts achieve significant SNR gains in the order of 40-80% over EPI in diffusion imaging at 3T. Combining systematic effects of shorter echo time, readout efficiency, and favorable g-factor behavior, similar benefits are expected across clinical and neurosciences uses of diffusion imaging.

Deconvolution-based distortion correction of EPI using analytic single-voxel point-spread functions

PURPOSE

To develop a postprocessing algorithm that corrects geometric distortions due to spatial variations of the static magnetic field amplitude, B₀ , and effects from relaxation during signal acquisition in EPI.

THEORY AND METHODS

An analytic, complex point-spread function is deduced for k-space trajectories of EPI variants and applied to corresponding acquisitions in a resolution phantom and in human volunteers at 3 T. With the analytic point-spread function and experimental maps of B₀ (and, optionally, the effective transverse relaxation time, T 2 * ) as input, a point-spread function matrix operator is devised for distortion correction by a Thikonov-regularized deconvolution in image space. The point-spread function operator provides additional information for an appropriate correction of the signal intensity distribution. A previous image combination algorithm for acquisitions with opposite phase blip polarities is adapted to the proposed method to recover destructively interfering signal contributions.

RESULTS

Applications of the proposed deconvolution-based distortion correction ("DecoDisCo") algorithm demonstrate excellent distortion corrections and superior performance regarding the recovery of an undistorted intensity distribution in comparison to a multifrequency reconstruction. Examples include full and partial Fourier standard EPI scans as well as double-shot center-out trajectories. Compared with other distortion-correction approaches, DecoDisCo permits additional deblurring to obtain sharper images in cases of significant T 2 * effects.

CONCLUSION

Robust distortion corrections in EPI acquisitions are feasible with high quality by regularized deconvolution with an analytic point-spread function. The general algorithm, which is publicly released on GitHub, can be straightforwardly adapted for specific EPI variants or other acquisition schemes.

Simultaneous Transcranial Magnetic Stimulation and Functional Magnetic Resonance Imaging: Aspects of Technical Implementation

The simultaneous transcranial magnetic stimulation (TMS) and functional magnetic resonance imaging (fMRI) offers a unique opportunity to non-invasively stimulate brain circuits while simultaneously monitoring changes in brain activity. However, to take advantage of this multimodal technique, some technical issues need to be considered/addressed. In this work, we evaluated technical issues associated with the setup and utilization of this multimodal tool, such as the use of a large single-channel radio frequency (rf) coil, and the artifacts induced by TMS when interleaved with the echo-planar imaging (EPI) sequence. We demonstrated that good image quality can be achieved with this rf coil and that the adoption of axial imaging orientation in conjunction with a safe interval of 100 ms, between the TMS pulse and imaging acquisition, is a suitable combination to eliminate potential image artifacts when using the combined TMS-fMRI technique in 3-T MRI scanners.

An equal-TE ultrafast 3D gradient-echo imaging method with high tolerance to magnetic susceptibility artifacts: Application to BOLD functional MRI

PURPOSE

To develop an ultrafast 3D gradient echo-based MRI method with constant TE and high tolerance to B₀ inhomogeneity, dubbed ERASE (equal-TE rapid acquisition with sequential excitation), and to introduce its use in BOLD functional MRI (fMRI).

THEORY AND METHODS

Essential features of ERASE, including spin behavior, were characterized, and a comparison study was conducted with conventional EPI. To demonstrate high tolerance to B₀ inhomogeneity, in vivo imaging of the mouse brain with a fiber-optic implant was performed at 9.4 T, and human brain imaging (including the orbitofrontal cortex) was performed at 3 T and 7 T. To evaluate the performance of ERASE in BOLD-fMRI, the characteristics of SNR and temporal SNR were analyzed for in vivo rat brains at 9.4 T in comparison with multislice gradient-echo EPI. Percent signal changes and t-scores are also presented.

RESULTS

For both mouse brain and human brain imaging, ERASE exhibited a high tolerance to magnetic susceptibility artifacts, showing much lower distortion and signal dropout, especially in the regions involving large magnetic susceptibility effects. For BOLD-fMRI, ERASE provided higher temporal SNR and t-scores than EPI, but exhibited similar percent signal changes in in vivo rat brains at 9.4 T.

CONCLUSION

When compared with conventional EPI, ERASE is much less sensitive, not only to EPI-related artifacts such as Nyquist ghosting, but also to B₀ inhomogeneity including magnetic susceptibility effects. It is promising for use in BOLD-fMRI, providing higher temporal SNR and t-scores with constant TE when compared with EPI, although further optimization is needed for human fMRI.

Prospective motion correction for diffusion weighted EPI of the brain using an optical markerless tracker

PURPOSE

To enable motion-robust diffusion weighted imaging of the brain using well-established imaging techniques.

METHODS

An optical markerless tracking system was used to estimate and correct for rigid body motion of the head in real time during scanning. The imaging coordinate system was updated before each excitation pulse in a single-shot EPI sequence accelerated by GRAPPA with motion-robust calibration. Full Fourier imaging was used to reduce effects of motion during diffusion encoding. Subjects were imaged while performing prescribed motion patterns, each repeated with prospective motion correction on and off.

RESULTS

Prospective motion correction with dynamic ghost correction enabled high quality DWI in the presence of fast and continuous motion within a 10° range. Images acquired without motion were not degraded by the prospective correction. Calculated diffusion tensors tolerated the motion well, but ADC values were slightly increased.

CONCLUSIONS

Prospective correction by markerless optical tracking minimizes patient interaction and appears to be well suited for EPI-based DWI of patient groups unable to remain still including those who are not compliant with markers.

Imaging quality of PROPELLER diffusion-weighted MR imaging and its diagnostic performance in distinguishing pleomorphic adenomas from Warthin tumors of the parotid gland

The aim of this study was to evaluate the imaging quality and diagnostic performance of fast spin echo diffusion-weighted imaging with periodically rotated overlapping parallel lines with enhanced reconstruction (FSE-PROP-DWI) in distinguishing parotid pleomorphic adenoma (PMA) from Warthin tumor (WT). This retrospective study enrolled 44 parotid gland tumors from 34 patients, including 15 PMAs and 29 WTs with waived written informed consent. All participants underwent 1.5 T diffusion-weighted imaging including FSE-PROP-DWI and single-shot echo-planar diffusion-weighted imaging (SS-EP-DWI). After imaging resizing and registration among T2WI, FSE-PROP-DWI and SS-EP-DWI, imaging distortion was quantitatively analyzed by using the Dice coefficient. Signal-to-noise ratio and contrast-to-noise ratio were qualitatively evaluated. The mean apparent diffusion coefficient (ADC) of parotid gland tumors was calculated. Wilcoxon signed-rank test was used for paired comparison between FSE-PROP-DWI versus SS-EP-DWI. Mann-Whitney U test was used for independent group comparison between PMAs versus WTs. Diagnostic performance was evaluated by receiver operating characteristics curve analysis. P < 0.05 was considered statistically significant. The Dice coefficient was statistically significantly higher on FSE-PROP-DWI than SS-EP-DWI for both tumors (P < 0.005). Mean ADC was statistically significantly higher in PMAs than WTs on both FSE-PROP-DWI and SS-EP-DWI (P < 0.005). FSE-PROP-DWI and SS-EP-DWI successfully distinguished PMAs from WTs with an AUC of 0.880 and 0.945, respectively (P < 0.05). Sensitivity, specificity, positive predictive value, negative predictive value and accuracy in diagnosing PMAs were 100%, 69.0%, 62.5%, 100% and 79.5% for FSE-PROP-DWI, and 100%, 82.8%, 75%, 100% and 88.6% for SS-EP-DWI, respectively. FSE-PROP-DWI is useful to distinguish parotid PMAs from WTs with less distortion of tumors but lower AUC than SS-EP-DWI.

Elimination of residual aliasing artifact that resembles brain lesion on multi-oblique diffusion-weighted echo-planar imaging with parallel imaging using virtual coil acquisition

BACKGROUND

Single-shot diffusion-weighted echo-planar imaging (ssDW-EPI) acquired with parallel imaging and a multi-oblique scan plane may suffer from residual aliasing artifacts, resembling lesions on the calculated apparent diffusion coefficient (ADC) map.

PURPOSE

To combine ssDW-EPI and virtual coil acquisition and develop a self-reference reconstruction method to eliminate the residual aliasing artifact on multi-oblique ssDW-EPI sequence with parallel imaging and multiple signal averaging.

STUDY TYPE

Prospective.

SUBJECTS

Three healthy subjects and 50 stroke patients.

FIELD STRENGTH/SEQUENCE

Conventional ssDW-EPI with parallel imaging, and proposed ssDW-EPI with virtual coil acquisition at 1.5T.

ASSESSMENT

The efficacy of the proposed method was evaluated in 50 stroke patients by comparing the ssDW-EPI with conventional parallel imaging reconstructions. The extent of residual aliasing artifacts were rated on a 5-point Likert scale by three independent raters. Only the data without residual aliasing artifacts on conventional ssDW-EPI were included for the assessment of signal-to-noise ratio (SNR), ghost-to-signal ratio (GSR), and ADC.

STATISTICAL TESTS

The interobserver agreements for examining residual aliasing artifacts were measured by the intraclass correlation coefficient (ICC). A two-sample t-test was performed for comparing SNR, GSR, and ADC.

RESULTS

There was a perfect agreement (ICC = 1.00) in the examination of residual aliasing artifacts on images obtained using the proposed method. The incidence rates of the residual aliasing artifact on the ADC maps obtained from the scanner console and proposed method were 60% (ie, 30 out of 50) and 0%, respectively. The proposed method offers significantly lower GSR than conventional parallel imaging reconstruction (P < 0.001). There was no significant difference in SNR (P = 0.20-0.51) and ADC values (P = 0.20-0.94) between conventional parallel imaging reconstructions and the proposed method.

DATA CONCLUSION

It appears that our method could effectively eliminate artifacts and significantly improve the GSR of b = 0 T₂ WI and b > 0 DWI, as well as permit ADC measurement consistent with conventional techniques. Our method may be beneficial to clinical assessment of the brain that utilizes multi-oblique ssDW-EPI.

LEVEL OF EVIDENCE

1 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2020;51:1442-1453.

A Comparison of Techniques for Correcting Eddy-current and Motion-induced Distortions in Diffusion-weighted Echo-planar Images

The purpose of this study was to show the efficacy of dynamic field correction (DFC), a technique provided by the scanner software, in comparison to the FMRIB Software Library (FSL) post-processing “eddy” tool. DFC requires minimal additional scan time for the correction of eddy-current and motion-induced geometrical image distortions in diffusion-weighted echo-planar images. The fractional anisotropy derived from images corrected with DFC were comparable to images corrected with the “eddy” tool and significantly higher than images without correction, which demonstrates the utility of DFC.

Data Quality and Optimal Background Correction Order of Respiratory-Gated k-Space Segmented Spoiled Gradient Echo (SGRE) and Echo Planar Imaging (EPI)-Based 4D Flow MRI

Background

A reduction in scan time of 4D Flow MRI would facilitate clinical application. A recent study indicates that echo‐planar imaging (EPI) 4D Flow MRI allows for a reduction in scan time and better data quality than the recommended k‐space segmented spoiled gradient echo (SGRE) sequence. It was argued that the poor data quality of SGRE was related to the nonrecommended absence of respiratory motion compensation. However, data quality can also be affected by the background offset compensation.

Purpose

To compare the data quality of respiratory motion‐compensated SGRE and EPI 4D Flow MRI and their dependence on background correction (BC) order.

Study Type

Retrospective.

Subjects

Eighteen healthy subjects (eight female, mean age 32 ± 5 years).

Field Strength and Sequence

1.5 T. [Correction added on July 26, 2019, after first online publication: The preceding field strength was corrected.] SGRE and EPI‐based 4D Flow MRI.

Assessment

Data quality was investigated visually and by comparing flows through the cardiac valves and aorta. Measurements were obtained from transvalvular flow and pathline analysis.

Statistical Tests

Linear regression and Bland–Altman analysis were used. Wilcoxon test was used for comparison of visual scoring. Student's t‐test was used for comparison of flow volumes.

Results

No significant difference was found by visual inspection (P = 0.08). Left ventricular (LV) flows were strongly and very strongly associated with SGRE and EPI, respectively (R² = 0.86–0.94 SGRE; 0.71–0.79 EPI, BC0–4). LV and right ventricular (RV) outflows and LV pathline flows were very strongly associated (R² = 0.93–0.95 SGRE; 0.88–0.91 EPI, R² = 0.91–0.95 SGRE; 0.91–0.93 EPI, BC1–4). EPI LV outflow was lower than the short‐axis‐based stroke volume. EPI RV outflow and proximal descending aortic flow were lower than SGREs.

Data Conclusion

Both sequences yielded good internal data consistency when an adequate background correction was applied. Second and first BC order were considered sufficient for transvalvular flow analysis in SGRE and EPI, respectively. Higher BC orders were preferred for particle tracing.

Level of Evidence 4

Technical Efficacy Stage 1

J. Magn. Reson. Imaging 2020;51:885–896.

Reduced field of view echo-planar imaging diffusion tensor MRI for pediatric spinal tumors

OBJECTIVE

Spine MRI is a diagnostic modality for evaluating pediatric CNS tumors. Applying diffusion-weighted MRI (DWI) or diffusion tensor imaging (DTI) to the spine poses challenges due to intrinsic spinal anatomy that exacerbates various image-related artifacts, such as signal dropouts or pileups, geometrical distortions, and incomplete fat suppression. The zonal oblique multislice (ZOOM)-echo-planar imaging (EPI) technique reduces geometric distortion and image blurring by reducing the field of view (FOV) without signal aliasing into the FOV. The authors hypothesized that the ZOOM-EPI method for spine DTI in concert with conventional spinal MRI is an efficient method for augmenting the evaluation of pediatric spinal tumors.

METHODS

Thirty-eight consecutive patients (mean age 8 years) who underwent ZOOM-EPI spine DTI for CNS tumor workup were retrospectively identified. Patients underwent conventional spine MRI and ZOOM-EPI DTI spine MRI. Two blinded radiologists independently reviewed two sets of randomized images: conventional spine MRI without ZOOM-EPI DTI, and conventional spine MRI with ZOOM-EPI DTI. For both image sets, the reviewers scored the findings based on lesion conspicuity and diagnostic confidence using a 5-point Likert scale. The reviewers also recorded presence of tumors. Quantitative apparent diffusion coefficient (ADC) measurements of various spinal tumors were extracted. Tractography was performed in a subset of patients undergoing presurgical evaluation.

RESULTS

Sixteen patients demonstrated spinal tumor lesions. The readers were in moderate agreement (kappa = 0.61, 95% CI 0.30-0.91). The mean scores for conventional MRI and combined conventional MRI and DTI were as follows, respectively: 3.0 and 4.0 for lesion conspicuity (p = 0.0039), and 2.8 and 3.9 for diagnostic confidence (p < 0.001). ZOOM-EPI DTI identified new lesions in 3 patients. In 3 patients, tractography used for neurosurgical planning showed characteristic fiber tract projections. The mean weighted ADCs of low- and high-grade tumors were 1201 × 10-6 and 865 × 10-6 mm2/sec (p = 0.002), respectively; the mean minimum weighted ADCs were 823 × 10-6 and 474 × 10-6 mm2/sec (p = 0.0003), respectively.

CONCLUSIONS

Diffusion MRI with ZOOM-EPI can improve the detection of spinal lesions while providing quantitative diffusion information that helps distinguish low- from high-grade tumors. By adding a 2-minute DTI scan, quantitative diffusion information and tract profiles can reliably be obtained and serve as a useful adjunct to presurgical planning for pediatric spinal tumors.

Value of High-Resolution DWI in Combination With Texture Analysis for the Evaluation of Tumor Response After Preoperative Chemoradiotherapy for Locally Advanced Rectal Cancer

OBJECTIVE. The purpose of this study is to determine the performance of the apparent diffusion coefficient (ADC) value calculated from high-resolution DWI using readout-segmented echo-planar imaging (rs-EPI) and to assess the texture parameters of T2-weighted MR images in identifying pathologic complete response (pCR) after patients with locally advanced rectal cancer (LARC) undergo preoperative chemoradiotherapy (CRT). MATERIALS AND METHODS. A total of 76 patients with LARC who underwent preoperative CRT and subsequent surgery were enrolled in the study retrospectively. All patients underwent post-CRT MRI, which included acquisition of a DWI sequence with use of the rs-EPI technique. The histopathologic tumor regression grade was the reference standard. Patients were subdivided into pCR and non-pCR groups. Two radiologists independently drew whole-tumor ROIs on DW images and T2-weighted MR images to calculate the mean ADC value and first-order texture parameters. RESULTS. Interobserver agreement was good to excellent (intraclass correlation coefficient [ICC], 0.79-0.993) for imaging analysis. Calculated from high-resolution DWI, the mean post-CRT ADC value was significantly higher in the pCR group (p < 0.001). The pCR group also showed lower uniformity (p < 0.001) of the T2-weighted image. The mean ADC value and uniformity were significantly correlated with the tumor regression grade. The mean ADC value was a good indicator for differentiating pCR from absence of pCR (ROC AUC value, 0.912). Uniformity (ROC AUC value, 0.776) showed a moderate ability to identify pCR. Combining the mean ADC value and uniformity yielded an ROC AUC value comparable to that of the mean ADC value (p = 0.125). CONCLUSION. Mean post-CRT ADC values calculated from high-resolution DWI using rs-EPI could effectively select for patients with LARC who have a pCR after preoperative CRT. First-order texture parameters of T2-weighted MR images could also identify patients with pCR by reflecting tumor heterogeneity, even though they could not significantly improve the diagnostic performance.

Comparison between 8- and 32-channel phased-array receive coils for in vivo hyperpolarized 13 C imaging of the human brain

PURPOSE

To compare the performance of an 8-channel surface coil/clamshell transmitter and 32-channel head array coil/birdcage transmitter for hyperpolarized ¹³ C brain metabolic imaging.

METHODS

To determine the field homogeneity of the radiofrequency transmitters, B₁ + mapping was performed on an ethylene glycol head phantom and evaluated by means of the double angle method. Using a 3D echo-planar imaging sequence, coil sensitivity and noise-only phantom data were acquired with the 8- and 32-channel receiver arrays, and compared against data from the birdcage in transceiver mode. Multislice frequency-specific ¹³ C dynamic echo-planar imaging was performed on a patient with a brain tumor for each hardware configuration following injection of hyperpolarized [1-¹³ C]pyruvate. Signal-to-noise ratio (SNR) was evaluated from pre-whitened phantom and temporally summed patient data after coil combination based on optimal weights.

RESULTS

The birdcage transmitter produced more uniform B₁ + compared with the clamshell: 0.07 versus 0.12 (fractional error). Phantom experiments conducted with matched lateral housing separation demonstrated 8- versus 32-channel mean transceiver-normalized SNR performance: 0.91 versus 0.97 at the head center; 6.67 versus 2.08 on the sides; 0.66 versus 2.73 at the anterior; and 0.67 versus 3.17 on the posterior aspect. While the 8-channel receiver array showed SNR benefits along lateral aspects, the 32-channel array exhibited greater coverage and a more uniform coil-combined profile. Temporally summed, parameter-normalized patient data showed SNR_mean,slice ratios (8-channel/32-channel) ranging 0.5-2.00 from apical to central brain. White matter lactate-to-pyruvate ratios were conserved across hardware: 0.45 ± 0.12 (8-channel) versus 0.43 ± 0.14 (32-channel).

CONCLUSION

The 8- and 32-channel hardware configurations each have advantages in particular brain anatomy.

Reducing distortions in echo-planar breast imaging at ultrahigh field with high-resolution off-resonance maps

PURPOSE

DWI is a promising modality in breast MRI, but its clinical acceptance is slow. Analysis of DWI is hampered by geometric distortion artifacts, which are caused by off-resonant spins in combination with the low phase-encoding bandwidth of the EPI sequence used. Existing correction methods assume smooth off-resonance fields, which we show to be invalid in the human breast, where high discontinuities arise at tissue interfaces.

METHODS

We developed a distortion correction method that incorporates high-resolution off-resonance maps to better solve for severe distortions at tissue interfaces. The method was evaluated quantitatively both ex vivo in a porcine tissue phantom and in vivo in 5 healthy volunteers. The added value of high-resolution off-resonance maps was tested using a Wilcoxon signed rank test comparing the quantitative results obtained with a low-resolution off-resonance map with those obtained with a high-resolution map.

RESULTS

Distortion correction using low-resolution off-resonance maps corrected most of the distortions, as expected. Still, all quantitative comparison metrics showed increased conformity between the corrected EPI images and a high-bandwidth reference scan for both the ex vivo and in vivo experiments. All metrics showed a significant improvement when a high-resolution off-resonance map was used (P < 0.05), in particular at tissue boundaries.

CONCLUSION

The use of off-resonance maps of a resolution higher than EPI scans significantly improves upon existing distortion correction techniques, specifically by superior correction at glandular tissue boundaries.

Diffusion-weighted MRI of the lung at 3T evaluated using echo-planar-based and single-shot turbo spin-echo-based acquisition techniques for radiotherapy applications

Purpose

To compare single‐shot echo‐planar (SS‐EPI)‐based and turbo spin‐echo (SS‐TSE)‐based diffusion‐weighted imaging (DWI) in Non‐Small Cell Lung Cancer (NSCLC) patients and to characterize the distributions of apparent diffusion coefficient (ADC) values generated by the two techniques.

Methods

Ten NSCLC patients were enrolled in a prospective IRB‐approved study to compare and optimize DWI using EPI and TSE‐based techniques for radiotherapy planning. The imaging protocol included axial T2w, EPI‐based DWI and TSE‐based DWI on a 3 T Philips scanner. Both EPI‐based and TSE‐based DWI sequences used three b values (0, 400, and 800 s/mm²). The acquisition times for EPI‐based and TSE‐based DWI were 5 and 8 min, respectively. DW‐MR images were manually coregistered with axial T2w images, and tumor volume contoured on T2w images were mapped onto the DWI scans. A pixel‐by‐pixel fit of tumor ADC was calculated based on monoexponential signal behavior. Tumor ADC mean, standard deviation, kurtosis, and skewness were calculated and compared between EPI and TSE‐based DWI. Image distortion and ADC values between the two techniques were also quantified using fieldmap analysis and a NIST traceable ice‐water diffusion phantom, respectively.

Results

The mean ADC for EPI and TSE‐based DWI were 1.282 ± 0.42 × 10⁻³ and 1.211 ± 0.31 × 10⁻³ mm²/s. The average skewness and kurtosis were 0.14 ± 0.4 and 2.43 ± 0.40 for DWI‐EPI and −0.06 ± 0.69 and 2.89 ± 0.62 for DWI‐TSE. Fieldmap analysis showed a mean distortion of 13.72 ± 8.12 mm for GTV for DWI‐EPI and 0.61 ± 0.4 mm for DWI‐TSE. ADC values obtained using the diffusion phantom for the two techniques were within 0.03 × 10⁻³ mm²/s with respect to each other as well as the established values.

Conclusions

Diffusion‐weighted turbo spin‐echo shows better geometrical accuracy compared to DWI‐EPI. Mean ADC values were similar with both acquisitions but the shape of the histograms was different based on the skewness and kurtosis values. The impact of differences in respiratory technique on ADC values requires further investigation.

Navigator-Free EPI Ghost Correction With Structured Low-Rank Matrix Models: New Theory and Methods

Structured low-rank matrix models have previously been introduced to enable calibrationless MR image reconstruction from sub-Nyquist data, and such ideas have recently been extended to enable navigator-free echo-planar imaging (EPI) ghost correction. This paper presents a novel theoretical analysis which shows that, because of uniform subsampling, the structured low-rank matrix optimization problems for EPI data will always have either undesirable or non-unique solutions in the absence of additional constraints. This theory leads us to recommend and investigate problem formulations for navigator-free EPI that incorporate side information from either image-domain or k-space domain parallel imaging methods. The importance of using nonconvex low-rank matrix regularization is also identified. We demonstrate using phantom and in vivo data that the proposed methods are able to eliminate ghost artifacts for several navigator-free EPI acquisition schemes, obtaining better performance in comparison with the state-of-the-art methods across a range of different scenarios. Results are shown for both single-channel acquisition and highly accelerated multi-channel acquisition.

Tilted-CAIPI for highly accelerated distortion-free EPI with point spread function (PSF) encoding

PURPOSE

To develop a method for fast distortion- and blurring-free imaging.

THEORY

EPI with point-spread-function (PSF) mapping can achieve distortion- and blurring-free imaging at a cost of long acquisition time. In this study, an acquisition/reconstruction technique, termed "tilted-CAIPI," is proposed to achieve >20× acceleration for PSF-EPI. The proposed method systematically optimized the k-space sampling trajectory with B₀ -inhomogeneity-informed reconstruction, to exploit the inherent signal correlation in PSF-EPI and take full advantage of coil sensitivity. Susceptibility-induced phase accumulation is regarded as an additional encoding that is estimated by calibration data and integrated into reconstruction. Self-navigated phase correction was developed to correct shot-to-shot phase variation in diffusion imaging.

METHODS

Tilted-CAIPI was implemented at 3T, with incorporation of partial Fourier and simultaneous multislice to achieve further accelerations. T₂ -weighted, T₂ ^* -weighted, and diffusion-weighted imaging experiments were conducted to evaluate the proposed method.

RESULTS

The ability of tilted-CAIPI to provide highly accelerated imaging without distortion and blurring was demonstrated through in vivo brain experiments, where only 8 shots per simultaneous slice group were required to provide high-quality, high-SNR imaging at 0.8-1 mm resolution.

CONCLUSION

Tilted-CAIPI achieved fast distortion- and blurring-free imaging with high SNR. Whole-brain T₂ -weighted, T₂ ^* -weighted, and diffusion imaging can be obtained in just 15-60 s.

ROBUST AUTOCALIBRATED LORAKS FOR EPI GHOST CORRECTION

Nyquist ghosts are a longstanding problem in a variety of fast MRI experiments that use echo-planar imaging (EPI). Recently, several structured low-rank matrix modeling approaches have been proposed that achieve state-of-the-art ghost-elimination, although the performance of these approaches is still inadequate in several important scenarios. We present a new structured low-rank matrix recovery ghost correction method that we call Robust Autocalibrated LORAKS (RAC-LORAKS). RAC-LORAKS incorporates constraints from autocalibration data to avoid ill-posedness, but allows adaptation of these constraints to gain robustness against possible autocalibration imperfections. RAC-LORAKS is tested in two challenging scenarios: highly-undersampled multi-channel EPI of the brain, and cardiac EPI with a double-oblique slice orientation. Results show that RAC-LORAKS can provide substantial improvements over existing ghost correction methods, and potentially enables new imaging applications that were previously confounded by ghost artifacts.

Volumetric MRI thermometry using a three-dimensional stack-of-stars echo-planar imaging pulse sequence

PURPOSE

To measure temperature over a large brain volume with fine spatiotemporal resolution.

METHODS

A three-dimensional stack-of-stars echo-planar imaging sequence combining echo-planar imaging and radial sampling with golden angle spacing was implemented at 3T for proton resonance frequency-shift temperature imaging. The sequence acquires a 188x188x43 image matrix with 1.5x1.5x2.75 mm3 spatial resolution. Temperature maps were reconstructed using sensitivity encoding (SENSE) image reconstruction followed by the image domain hybrid method, and using the k-space hybrid method. In vivo temperature maps were acquired without heating to measure temperature precision in the brain, and in a phantom during high-intensity focused ultrasound sonication.

RESULTS

In vivo temperature standard deviation was less than 1°C at dynamic scan times down to 0.75 s. For a given frame rate, scanning at a minimum repetition time (TR) with minimum acceleration yielded the lowest standard deviation. With frame rates around 3 s, the scan was tolerant to a small number of receive coils, and temperature standard deviation was 48% higher than a standard two-dimensional Fourier transform temperature mapping scan, but provided whole-brain coverage. Phantom temperature maps with no visible aliasing were produced for dynamic scan times as short as 0.38 s. k-Space hybrid reconstructions were more tolerant to acceleration.

CONCLUSION

Three-dimensional stack-of-stars echo-planar imaging temperature mapping provides volumetric brain coverage and fine spatiotemporal resolution. Magn Reson Med 79:2003-2013, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Comparison of DWI Methods in the Pediatric Brain: PROPELLER Turbo Spin-Echo Imaging Versus Readout-Segmented Echo-Planar Imaging Versus Single-Shot Echo-Planar Imaging

OBJECTIVE

The purpose of this study is to compare DWI for pediatric brain evaluation using single-shot echo-planar imaging (EPI), periodically rotated overlapping parallel lines with enhanced reconstruction (Blade), and readout-segmented EPI (Resolve).

MATERIALS AND METHODS

Blade, Resolve, and single-shot EPI were performed for 27 pediatric patients (median age, 9 years), and three datasets were independently reviewed by two radiologists. Qualitative analyses were performed for perceptive coarseness, image distortion, susceptibility-related changes, motion artifacts, and lesion conspicuity using a 5-point Likert scale. Quantitative analyses were conducted for spatial distortion and signal uniformity of each sequence.

RESULTS

Mean scores were 2.13, 3.17, and 3.76 for perceptive coarseness; 4.85, 3.96, and 2.19 for image distortion; 4.76, 3.96, and 2.30 for susceptibility-related change; 4.96, 3.83, and 4.69 for motion artifacts; and 2.71, 3.75, and 1.92 for lesion conspicuity, for Blade, Resolve, and single-shot EPI, respectively. Blade and Resolve showed better quality than did single-shot EPI for image distortion, susceptibility-related changes, and lesion conspicuity. Blade showed less image distortion, fewer susceptibility-related changes, and fewer motion artifacts than did Resolve, whereas lesion conspicuity was better with Resolve. Blade showed increased signal variation compared with Resolve and single-shot EPI (coefficients of variation were 0.10, 0.08, and 0.05 for lateral ventricle; 0.13, 0.09, and 0.05 for centrum semiovale; and 0.16, 0.09, and 0.06 for pons in Blade, Resolve, and single-shot EPI, respectively).

CONCLUSION

DWI with Resolve or Blade yields better quality regarding distortion, susceptibility-related changes, and lesion conspicuity, compared with single-shot EPI. Blade is less susceptible to motion artifacts than is Resolve, whereas Resolve yields less noise and better lesion conspicuity than does Blade.

Proton change of parotid glands after gustatory stimulation examined by magnetic resonance imaging

The aim of this study was to investigate proton changes of the parotid gland after gustatory stimulation by semi-quantitative parameters and an empirical mathematical model (EMM) using high-temporal-resolution, double-echo, echo-planar imaging (EPI). Approved by a local institutional review board, this study examined 20 parotid glands from 10 healthy volunteers (male:female = 6: 4; age ± standard deviation =35.1 ± 14.1 years) with written informed consent obtained. All participants underwent 1.5-T, double-echo EPI with gustatory stimulation. Semi-quantitative parameters, including maximal drop ratio (MDR), time to peak (TTP), drop slope (DS), recovery slope (RS) and recovery ratio (RR), were calculated. The effect of temporal resolution on parotid functional parameters was evaluated. An EMM comprising an output function ( Sot=Aoe-kot+B) and an input function ( Sint=Ain1-e-kint) was also applied to fit all dynamic curves. Kruskal-Wallis test, Wilcoxon test, linear regression analysis and goodness of fit were used for statistical analysis. p < 0.05 was considered to be statistically significant. The signal intensity dropped significantly after gustatory stimulation on the proton density (PD) image (p < 0.01). MDR was 8.26% in the PD image. MDR and RR were negatively associated with time interval, whereas DS and TTP were significantly positively associated with time interval (all p < 0.05). EMM parametric values derived from PD-time curves of parotid glands were 12.04 ± 6.81%, 6.43 ± 4.23 min^-1 , 88.73 ± 6.18%, 8.41 ± 4.86 min^-1 and 1.09 ± 1.35 for A_o , k_o , B, A_in and k_in , respectively. Semi-quantitative functional parameters and EMM parameters using high-temporal-resolution, double-echo EPI allow the quantification of parotid proton changes after gustatory stimulation.

Connectome analysis for pre-operative brain mapping in neurosurgery

Object: Brain mapping has entered a new era focusing on complex network connectivity. Central to this is the search for the connectome or the brains ‘wiring diagram’. Graph theory analysis of the connectome allows understanding of the importance of regions to network function, and the consequences of their impairment or excision. Our goal was to apply connectome analysis in patients with brain tumours to characterise overall network topology and individual patterns of connectivity alterations.

Methods: Resting-state functional MRI data were acquired using multi-echo, echo planar imaging pre-operatively from five participants each with a right temporal–parietal–occipital glioblastoma. Complex networks analysis was initiated by parcellating the brain into anatomically regions amongst which connections were identified by retaining the most significant correlations between the respective wavelet decomposed time-series.

Results: Key characteristics of complex networks described in healthy controls were preserved in these patients, including ubiquitous small world organization. An exponentially truncated power law fit to the degree distribution predicted findings of general network robustness to injury but with a core of hubs exhibiting disproportionate vulnerability. Tumours produced a consistent reduction in local and long-range connectivity with distinct patterns of connection loss depending on lesion location.

Conclusions: Connectome analysis is a feasible and novel approach to brain mapping in individual patients with brain tumours. Applications to pre-surgical planning include identifying regions critical to network function that should be preserved and visualising connections at risk from tumour resection. In the future one could use such data to model functional plasticity and recovery of cognitive deficits.

The pulsatility volume index: an indicator of cerebrovascular compliance based on fast magnetic resonance imaging of cardiac and respiratory pulsatility

The influence of cardiac activity on the viscoelastic properties of intracranial tissue is one of the mechanisms through which brain-heart interactions take place, and is implicated in cerebrovascular disease. Cerebrovascular disease risk is not fully explained by current risk factors, including arterial compliance. Cerebrovascular compliance is currently estimated indirectly through Doppler sonography and magnetic resonance imaging (MRI) measures of blood velocity changes. In order to meet the need for novel cerebrovascular disease risk factors, we aimed to design and validate an MRI indicator of cerebrovascular compliance based on direct endogenous measures of blood volume changes. We implemented a fast non-gated two-dimensional MRI pulse sequence based on echo-planar imaging (EPI) with ultra-short repetition time (approx. 30-50 ms), which stepped through slices every approximately 20 s. We constrained the solution of the Bloch equations for spins moving faster than a critical speed to produce an endogenous contrast primarily dependent on spin volume changes, and an approximately sixfold signal gain compared with Ernst angle acquisitions achieved by the use of a 90° flip angle. Using cardiac and respiratory peaks detected on physiological recordings, average cardiac and respiratory MRI pulse waveforms in several brain compartments were obtained at 7 Tesla, and used to derive a compliance indicator, the pulsatility volume index (pVI). The pVI, evaluated in larger cerebral arteries, displayed significant variation within and across vessels. Multi-echo EPI showed the presence of significant pulsatility effects in both S0 and [Formula: see text] signals, compatible with blood volume changes. Lastly, the pVI dynamically varied during breath-holding compared with normal breathing, as expected for a compliance indicator. In summary, we characterized and performed an initial validation of a novel MRI indicator of cerebrovascular compliance, which might prove useful to investigate brain-heart interactions in cerebrovascular disease and other disorders.

Single-scan MRI with exceptional resilience to field heterogeneities

PURPOSE

Single-scan two-dimensional MRI has been generally constrained to acquisitions in high quality magnets. This study introduces a methodology, cross-term spatiotemporal encoding (xSPEN), that delivers such images under much poorer external field conditions.

METHODS

xSPEN departs from conventional k-space scanning, by relying on spatiotemporally encoding the image being sought. Unlike hitherto proposed SPEN methods, however, xSPEN's image readout does not take place using a field gradient along the direction being probed, but rather with the aid of an ancillary source of inhomogeneous frequency broadening. This ancillary dimension was here imposed by an orthogonal field gradient; for example, images along the "y" axis were read out by application of a "z" gradient. The principles and characteristics of this new approach, compatible with existing scanners and free from the need to collect auxiliary information such as field maps, are presented and discussed.

RESULTS

Single- and multi-slice in vitro, ex vivo, and in vivo MRI experiments, confirmed the unusual resilience of this new single-shot MRI method to multiple chemical sites on phantoms, animals and humans.

CONCLUSION

xSPEN can deliver single-scan MRI with good sensitivity and exceptional resilience to field inhomogeneities. This could enable investigations that have hitherto escaped from MRI's scope. Magn Reson Med 77:623-634, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Reference-free single-pass EPI Nyquist ghost correction using annihilating filter-based low rank Hankel matrix (ALOHA)

PURPOSE

MR measurements from an echo-planar imaging (EPI) sequence produce Nyquist ghost artifacts that originate from inconsistencies between odd and even echoes. Several reconstruction algorithms have been proposed to reduce such artifacts, but most of these methods require either additional reference scans or multipass EPI acquisition. This article proposes a novel and accurate single-pass EPI ghost artifact correction method that does not require any additional reference data.

THEORY AND METHODS

After converting a ghost correction problem into separate k-space data interpolation problems for even and odd phase encoding, our algorithm exploits an observation that the differential k-space data between the even and odd echoes is a Fourier transform of an underlying sparse image. Accordingly, we can construct a rank-deficient Hankel structured matrix, whose missing data can be recovered using an annihilating filter-based low rank Hankel structured matrix completion approach.

RESULTS

The proposed method was applied to EPI data for both single and multicoil acquisitions. Experimental results using in vivo data confirmed that the proposed method can completely remove ghost artifacts successfully without prescan echoes.

CONCLUSION

Owing to the discovery of the annihilating filter relationship from the intrinsic EPI image property, the proposed method successfully suppresses ghost artifacts without a prescan step. Magn Reson Med 76:1775-1789, 2016. © 2016 International Society for Magnetic Resonance in Medicine.

Accelerated 3D echo-planar imaging with compressed sensing for time-resolved hyperpolarized 13 C studies

PURPOSE

To enable large field-of-view, time-resolved volumetric coverage in hyperpolarized ¹³ C metabolic imaging by implementing a novel data acquisition and image reconstruction method based on the compressed sensing framework.

METHODS

A spectral-spatial pulse for single-resonance excitation followed by a symmetric echo-planar imaging (EPI) readout was implemented for encoding a 72 × 18 cm² field of view at 5 × 5 mm² resolution. Random undersampling was achieved with blipped z-gradients during the ramp portion of the echo-planar imaging readout. The sequence and reconstruction were tested with phantom studies and consecutive in vivo hyperpolarized ¹³ C scans in rats. Retrospectively and prospectively undersampled data were compared on the basis of structural similarity in the reconstructed images and the quantification of the lactate-to-pyruvate ratio in rat kidneys.

RESULTS

No artifacts or loss of resolution are evident in the compressed sensing reconstructed images acquired with the proposed sequence. Structural similarity analysis indicate that compressed sensing reconstructions can accurately recover spatial features in the metabolic images evaluated.

CONCLUSION

A novel z-blip acquisition sequence for compressed sensing accelerated hyperpolarized ¹³ C 3D echo-planar imaging was developed and demonstrated. The close agreement in lactate-to-pyruvate ratios from both retrospectively and prospectively undersampled data from rats shows that metabolic information is preserved with acceleration factors up to 3-fold with the developed method. Magn Reson Med 77:538-546, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Motion immune diffusion imaging using augmented MUSE for high-resolution multi-shot EPI

PURPOSE

To develop new techniques for reducing the effects of microscopic and macroscopic patient motion in diffusion imaging acquired with high-resolution multishot echo-planar imaging.

THEORY

The previously reported multiplexed sensitivity encoding (MUSE) algorithm is extended to account for macroscopic pixel misregistrations, as well as motion-induced phase errors in a technique called augmented MUSE (AMUSE). Furthermore, to obtain more accurate quantitative diffusion-tensor imaging measures in the presence of subject motion, we also account for the altered diffusion encoding among shots arising from macroscopic motion.

METHODS

MUSE and AMUSE were evaluated on simulated and in vivo motion-corrupted multishot diffusion data. Evaluations were made both on the resulting imaging quality and estimated diffusion tensor metrics.

RESULTS

AMUSE was found to reduce image blurring resulting from macroscopic subject motion compared to MUSE but yielded inaccurate tensor estimations when neglecting the altered diffusion encoding. Including the altered diffusion encoding in AMUSE produced better estimations of diffusion tensors.

CONCLUSION

The use of AMUSE allows for improved image quality and diffusion tensor accuracy in the presence of macroscopic subject motion during multishot diffusion imaging. These techniques should facilitate future high-resolution diffusion imaging.

Dual optimization method of radiofrequency and quasistatic field simulations for reduction of eddy currents generated on 7T radiofrequency coil shielding

PURPOSE

To optimize the design of radiofrequency (RF) shielding of transmit coils at 7T and reduce eddy currents generated on the RF shielding when imaging with rapid gradient waveforms.

METHODS

One set of a four-element, 2 × 2 Tic-Tac-Toe head coil structure was selected and constructed to study eddy currents on the RF coil shielding. The generated eddy currents were quantitatively studied in the time and frequency domains. The RF characteristics were studied using the finite difference time domain method. Five different kinds of RF shielding were tested on a 7T MRI scanner with phantoms and in vivo human subjects.

RESULTS

The eddy current simulation method was verified by the measurement results. Eddy currents induced by solid/intact and simple-structured slotted RF shielding significantly distorted the gradient fields. Echo-planar images, B1+ maps, and S matrix measurements verified that the proposed slot pattern suppressed the eddy currents while maintaining the RF characteristics of the transmit coil.

CONCLUSION

The presented dual-optimization method could be used to design RF shielding and reduce the gradient field-induced eddy currents while maintaining the RF characteristics of the transmit coil.

Interleaved diffusion-weighted improved by adaptive partial-Fourier and multiband multiplexed sensitivity-encoding reconstruction

PURPOSE

We report a series of techniques to reliably eliminate artifacts in interleaved echo-planar imaging (EPI) based diffusion-weighted imaging (DWI).

METHODS

First, we integrate the previously reported multiplexed sensitivity encoding (MUSE) algorithm with a new adaptive Homodyne partial-Fourier reconstruction algorithm, so that images reconstructed from interleaved partial-Fourier DWI data are free from artifacts even in the presence of either (a) motion-induced k-space energy peak displacement, or (b) susceptibility field gradient induced fast phase changes. Second, we generalize the previously reported single-band MUSE framework to multiband MUSE, so that both through-plane and in-plane aliasing artifacts in multiband multishot interleaved DWI data can be effectively eliminated.

RESULTS

The new adaptive Homodyne-MUSE reconstruction algorithm reliably produces high-quality and high-resolution DWI, eliminating residual artifacts in images reconstructed with previously reported methods. Furthermore, the generalized MUSE algorithm is compatible with multiband and high-throughput DWI.

CONCLUSION

The integration of the multiband and adaptive Homodyne-MUSE algorithms significantly improves the spatial-resolution, image quality, and scan throughput of interleaved DWI. We expect that the reported reconstruction framework will play an important role in enabling high-resolution DWI for both neuroscience research and clinical uses.

Readout-segmented echo-planar imaging for diffusion-weighted imaging in the pelvis at 3T-A feasibility study

RATIONALE AND OBJECTIVES

Diffusion-weighted imaging (DWI) of the pelvis at 3T is prone to artifacts that diminish the image quality. Readout-segmented echo-planar imaging (RS-EPI) is a new DWI technique that can reduce the artifacts associated with standard single-shot echo-planar imaging (SS-EPI) DWI. The purpose of this study was to evaluate the feasibility and image quality of RS-EPI in pelvic DWI compared to SS-EPI on a 3T imaging system.

MATERIALS AND METHODS

Thirty patients underwent pelvic DWI on a 3T scanner with SS-EPI and RS-EPI techniques. Two blinded readers independently assessed each set of images for geometric distortion, image blurring, ghosting artifacts, lesion conspicuity, and overall image quality on a 7-point scale. Qualitative image scores were compared using paired Wilcoxon signed rank test. Interreader correlation was assessed by Spearman rank correlation.

RESULTS

Geometric distortion, imaging blurring, ghosting artifacts, lesion conspicuity, and overall image quality were rated significantly better by both readers for RS-EPI technique (P < .01 for all parameters). There was moderate-high correlation between the readers (r = 0.649-0.752) for all parameters apart from lesion conspicuity (r = 0.351). Both readers preferred the RS-EPI set of DWI images in most of the cases (reader 1: 0.87, 95% CI 0.74-0.99; reader 2: 0.77, 95% CI 0.61-0.93). Mean difference and limits of agreement between apparent diffusion coefficient (ADC) values obtained from the two methods were 0.01 (-0.08, 0.10) × 10(-3) mm(2)/s.

CONCLUSIONS

RS-EPI DWI images showed improved image quality compared to SS-EPI technique at 3T. RS-EPI is a feasible technique in the pelvis for producing high-resolution DWI.

Frequency correction method for improved spatial correlation of hyperpolarized 13C metabolites and anatomy

Blip-reversed echo-planar imaging (EPI) is investigated as a method for measuring and correcting the spatial shifts that occur due to bulk frequency offsets in (13)C metabolic imaging in vivo. By reversing the k-space trajectory for every other time point, the direction of the spatial shift for a given frequency is reversed. Here, mutual information is used to find the 'best' alignment between images and thereby measure the frequency offset. Time-resolved 3D images of pyruvate/lactate/urea were acquired with 5 s temporal resolution over a 1 min duration in rats (N = 6). For each rat, a second injection was performed with the demodulation frequency purposely mis-set by +35 Hz, to test the correction for erroneous shifts in the images. Overall, the shift induced by the 35 Hz frequency offset was 5.9 ± 0.6 mm (mean ± standard deviation). This agrees well with the expected 5.7 mm shift based on the 2.02 ms delay between k-space lines (giving 30.9 Hz per pixel). The 0.6 mm standard deviation in the correction corresponds to a frequency-detection accuracy of 4 Hz. A method was presented for ensuring the spatial registration between (13)C metabolic images and conventional anatomical images when long echo-planar readouts are used. The frequency correction method was shown to have an accuracy of 4 Hz. Summing the spatially corrected frames gave a signal-to-noise ratio (SNR) improvement factor of 2 or greater, compared with the highest single frame.

Strategies for assessing diffusion anisotropy on the basis of magnetic resonance images: comparison of systematic errors

Diffusion weighted imaging uses the signal loss associated with the random thermal motion of water molecules in the presence of magnetic field gradients to derive a number of parameters that reflect the translational mobility of the water molecules in tissues. With a suitable experimental set-up, it is possible to calculate all the elements of the local diffusion tensor (DT) and derived parameters describing the behavior of the water molecules in each voxel. One of the emerging applications of the information obtained is an interpretation of the diffusion anisotropy in terms of the architecture of the underlying tissue. These interpretations can only be made provided the experimental data which are sufficiently accurate. However, the DT results are susceptible to two systematic error sources: On one hand, the presence of signal noise can lead to artificial divergence of the diffusivities. In contrast, the use of a simplified model for the interaction of the protons with the diffusion weighting and imaging field gradients (b matrix calculation), common in the clinical setting, also leads to deviation in the derived diffusion characteristics. In this paper, we study the importance of these two sources of error on the basis of experimental data obtained on a clinical magnetic resonance imaging system for an isotropic phantom using a state of the art single-shot echo planar imaging sequence. Our results show that optimal diffusion imaging require combining a correct calculation of the b-matrix and a sufficiently large signal to noise ratio.

Improved B0 -distortion correction in diffusion MRI using interlaced q-space sampling and constrained reconstruction

PURPOSE

To enable high-quality correction of susceptibility-induced geometric distortion artifacts in diffusion magnetic resonance imaging (MRI) images without increasing scan time.

THEORY AND METHODS

A new method for distortion correction is proposed based on subsampling a generalized version of the state-of-the-art reversed-gradient distortion correction method. Rather than acquire each q-space sample multiple times with different distortions (as in the conventional reversed-gradient method), we sample each q-space point once with an interlaced sampling scheme that measures different distortions at different q-space locations. Distortion correction is achieved using a novel constrained reconstruction formulation that leverages the smoothness of diffusion data in q-space.

RESULTS

The effectiveness of the proposed method is demonstrated with simulated and in vivo diffusion MRI data. The proposed method is substantially faster than the reversed-gradient method, and can also provide smaller intensity errors in the corrected images and smaller errors in derived quantitative diffusion parameters.

CONCLUSION

The proposed method enables state-of-the-art distortion correction performance without increasing data acquisition time.

Is there a cerebellar compensatory effort in first-episode, treatment-naive major depressive disorder at rest?

BACKGROUND

This study was undertaken to explore whether there is a cerebellar compensatory response in patients with first-episode, treatment-naive major depressive disorder (MDD). The cerebellar compensatory response is defined as a cerebellar hyperactivity which would be inversely correlated with both the activation of the functionally connected cerebral regions and the depression severity.

METHODS

Resting-state functional magnetic resonance imaging (fMRI) data of 24 patients with MDD and 24 healthy subjects were analyzed with the fractional amplitude of low-frequency fluctuations (fALFF) and functional connectivity (FC) methods. The structural images were processed with the voxel-based morphometry (VBM) method.

RESULTS

Compared to healthy controls, depressed patients had significantly increased fALFF in the left Crus I and the left cerebellar lobule VI. FC analysis of these two seeded regions found that depressed patients had increased FC between the left Crus I and the right hippocampus, but had decreased FC between the left Crus I and the left inferior parietal lobule (IPL), and between the left cerebellar lobule VI and bilateral inferior temporal gyrus. No correlation was observed between the abnormal fALFF of the seeds and their connected regions and the depression severity or the executive function. The VBM results did not show significant reduction in gray or white matter volume in any above-mentioned region.

CONCLUSIONS

Our findings suggest that increased cerebellar activity at resting state may be a disease state phenomenon but not a compensatory response to the dysfunction of the default mode network (DMN) in MDD.

Direct design of 2D RF pulses using matrix inversion

Multi-dimensional pulses are frequently used in MRI for applications such as targeted excitation, fat-water separation or metabolic imaging with hyperpolarised (13)C compounds. For the design, the problem is typically separated into the different dimensions. In this work, a method to directly design two-dimensional pulses using the small-tip angle approximation is introduced based on a direct matrix representation. The numerical problem is solved in a single step directly in two dimensions by matrix inversion. Exemplary spectral-spatial excitation and spatio-temporal encoding (SPEN) pulses are designed and validated. The main benefits of the direct design approach include a reduction of artefacts in case of spectral-spatial pulses, a simple and straightforward computer implementation and high flexibility in the pulse design.

Quasi-periodic patterns (QPP): large-scale dynamics in resting state fMRI that correlate with local infraslow electrical activity

Functional connectivity measurements from resting state blood-oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) are proving a powerful tool to probe both normal brain function and neuropsychiatric disorders. However, the neural mechanisms that coordinate these large networks are poorly understood, particularly in the context of the growing interest in network dynamics. Recent work in anesthetized rats has shown that the spontaneous BOLD fluctuations are tightly linked to infraslow local field potentials (LFPs) that are seldom recorded but comparable in frequency to the slow BOLD fluctuations. These findings support the hypothesis that long-range coordination involves low frequency neural oscillations and establishes infraslow LFPs as an excellent candidate for probing the neural underpinnings of the BOLD spatiotemporal patterns observed in both rats and humans. To further examine the link between large-scale network dynamics and infraslow LFPs, simultaneous fMRI and microelectrode recording were performed in anesthetized rats. Using an optimized filter to isolate shared components of the signals, we found that time-lagged correlation between infraslow LFPs and BOLD is comparable in spatial extent and timing to a quasi-periodic pattern (QPP) found from BOLD alone, suggesting that fMRI-measured QPPs and the infraslow LFPs share a common mechanism. As fMRI allows spatial resolution and whole brain coverage not available with electroencephalography, QPPs can be used to better understand the role of infraslow oscillations in normal brain function and neurological or psychiatric disorders.

Investigating univariate temporal patterns for intrinsic connectivity networks based on complexity and low-frequency oscillation: a test-retest reliability study

Intrinsic connectivity networks (ICNs) are composed of spatial components and time courses. The spatial components of ICNs were discovered with moderate-to-high reliability. So far as we know, few studies focused on the reliability of the temporal patterns for ICNs based their individual time courses. The goals of this study were twofold: to investigate the test-retest reliability of temporal patterns for ICNs, and to analyze these informative univariate metrics. Additionally, a correlation analysis was performed to enhance interpretability. Our study included three datasets: (a) short- and long-term scans, (b) multi-band echo-planar imaging (mEPI), and (c) eyes open or closed. Using dual regression, we obtained the time courses of ICNs for each subject. To produce temporal patterns for ICNs, we applied two categories of univariate metrics: network-wise complexity and network-wise low-frequency oscillation. Furthermore, we validated the test-retest reliability for each metric. The network-wise temporal patterns for most ICNs (especially for default mode network, DMN) exhibited moderate-to-high reliability and reproducibility under different scan conditions. Network-wise complexity for DMN exhibited fair reliability (ICC<0.5) based on eyes-closed sessions. Specially, our results supported that mEPI could be a useful method with high reliability and reproducibility. In addition, these temporal patterns were with physiological meanings, and certain temporal patterns were correlated to the node strength of the corresponding ICN. Overall, network-wise temporal patterns of ICNs were reliable and informative and could be complementary to spatial patterns of ICNs for further study.

The Importance of k-Space Trajectory on Off-Resonance Artifact in Segmented Echo-Planar Imaging

Segmented interleaved echo planar imaging (EPI) is a highly efficient data acquisition technique; however, EPI is sensitive to artifacts from off-resonance spins. The choice of k-space trajectories is important in determining how off-resonance spins contribute to image artifacts. Top-down and center-out trajectories are theoretically analyzed, simulated, implemented, and tested in phantom and volunteer experiments. Theoretical results show off-resonance artifact manifests as a simple positional shift for the top-down trajectory, while for the center-out trajectory off-resonance artifact manifests as a splitting of the object, which entails both shift and blurring. These results were validated using simulation and phantom scan data where a frequency-offset was introduced ranging from -300 Hz to +300 Hz. As predicted by the theoretical results, inferior image quality was observed for the center-out trajectory in a single volunteer. Off-resonance produces more severe and complex artifacts with the center-out trajectory than the top-down trajectory.

Robust GRAPPA-accelerated diffusion-weighted readout-segmented (RS)-EPI

Readout segmentation (RS-EPI) has been suggested as a promising variant to echo-planar imaging (EPI) for high-resolution imaging, particularly when combined with parallel imaging. This work details some of the technical aspects of diffusion-weighted (DW)-RS-EPI, outlining a set of reconstruction methods and imaging parameters that can both minimize the scan time and afford high-resolution diffusion imaging with reduced distortions. These methods include an efficient generalized autocalibrating partially parallel acquisition (GRAPPA) calibration for DW-RS-EPI data without scan time penalty, together with a variant for the phase correction of partial Fourier RS-EPI data. In addition, the role of pulsatile and rigid-body brain motion in DW-RS-EPI was assessed. Corrupt DW-RS-EPI data arising from pulsatile nonlinear brain motion had a prevalence of approximately 7% and were robustly identified via k-space entropy metrics. For DW-RS-EPI data corrupted by rigid-body motion, we showed that no blind overlap was required. The robustness of RS-EPI toward phase errors and motion, together with its minimized distortions compared with EPI, enables the acquisition of exquisite 3 T DW images with matrix sizes close to 512(2).