High efficiency, low distortion 3D diffusion tensor imaging with variable density spiral fast spin echoes (3D DW VDS RARE)

3D diffusion MRI with twin navigator-based GRASE and comparison with 2D EPI for tractography in the human brain

PURPOSE

3D pulse sequences enable high-resolution acquisition with a high SNR and ideal slice profiles, which, however, is particularly difficult for diffusion MRI (dMRI) due to the additional phase errors from diffusion encoding.

METHODS

We proposed a twin navigator-based 3D diffusion-weighted gradient spin-echo (GRASE) sequence to correct the phase errors between shots and between odd and even spin echoes for human whole-brain acquisition. We then compared the SNR of 3D GRASE and 2D simultaneous multi-slice EPI within the same acquisition time. We further tested the performance of 2D versus 3D acquisition at equivalent SNR on fiber tracking and microstructural mapping, using the diffusion tensor and high-order fiber orientation density-based metrics.

RESULTS

The proposed twin navigator approach removed multi-shot phase errors to some extent in the whole brain dMRI, and the 2D navigator performed better than the 1D navigator. Comparisons of SNR between the 2D simultaneous multi-slice EPI and 3D GRASE sequences demonstrated that the SNR of the GRASE sequence was 1.4-1.5-fold higher than the EPI sequence at an equivalent scan time. More importantly, we found a significantly higher fiber cross-section in the cerebrospinal tract, as well as richer subcortical fibers (U-fibers) using the 3D GRASE sequence compared to 2D EPI.

CONCLUSION

The twin navigator-based 3D diffusion-weighted-GRASE sequence minimized the multishot phase error and effectively improved the SNR for whole-brain dMRI acquisition. We found differences in fiber tracking and microstructural mapping between 2D and 3D acquisitions, possibly due to the different slice profiles.

Sampling strategies and integrated reconstruction for reducing distortion and boundary slice aliasing in high-resolution 3D diffusion MRI

PURPOSE

To develop a new method for high-fidelity, high-resolution 3D multi-slab diffusion MRI with minimal distortion and boundary slice aliasing.

METHODS

Our method modifies 3D multi-slab imaging to integrate blip-reversed acquisitions for distortion correction and oversampling in the slice direction (kz ) for reducing boundary slice aliasing. Our aim is to achieve robust acceleration to keep the scan time the same as conventional 3D multi-slab acquisitions, in which data are acquired with a single direction of blip traversal and without kz -oversampling. We employ a two-stage reconstruction. In the first stage, the blip-up/down images are respectively reconstructed and analyzed to produce a field map for each diffusion direction. In the second stage, the blip-reversed data and the field map are incorporated into a joint reconstruction to produce images that are corrected for distortion and boundary slice aliasing.

RESULTS

We conducted experiments at 7T in six healthy subjects. Stage 1 reconstruction produces images from highly under-sampled data (R = 7.2) with sufficient quality to provide accurate field map estimation. Stage 2 joint reconstruction substantially reduces distortion artifacts with comparable quality to fully-sampled blip-reversed results (2.4× scan time). Whole-brain in-vivo results acquired at 1.22 mm and 1.05 mm isotropic resolutions demonstrate improved anatomical fidelity compared to conventional 3D multi-slab imaging. Data demonstrate good reliability and reproducibility of the proposed method over multiple subjects.

CONCLUSION

The proposed acquisition and reconstruction framework provide major reductions in distortion and boundary slice aliasing for 3D multi-slab diffusion MRI without increasing the scan time, which can potentially produce high-quality, high-resolution diffusion MRI.

Mono-planar T-Hex: Speed and flexibility for high-resolution 3D imaging

Purpose

The aim of this work is the reconciliation of high spatial and temporal resolution for MRI. For this purpose, a novel sampling strategy for 3D encoding is proposed, which provides flexible k‐space segmentation along with uniform sampling density and benign filtering effects related to signal decay.

Methods

For time‐critical MRI applications such as functional MRI (fMRI), 3D k‐space is usually sampled by stacking together 2D trajectories such as echo planar imaging (EPI) or spiral readouts, where each shot covers one k‐space plane. For very high temporal and medium to low spatial resolution, tilted hexagonal sampling (T‐Hex) was recently proposed, which allows the acquisition of a larger k‐space volume per excitation than can be covered with a planar readout. Here, T‐Hex is described in a modified version where it instead acquires a smaller k‐space volume per shot for use with medium temporal and high spatial resolution.

Results

Mono‐planar T‐Hex sampling provides flexibility in the choice of speed, signal‐to‐noise ratio (SNR), and contrast for rapid MRI acquisitions. For use with a conventional gradient system, it offers the greatest benefit in a regime of high in‐plane resolution <1 mm. The sampling scheme is combined with spirals for high sampling speed as well as with more conventional EPI trajectories.

Conclusion

Mono‐planar T‐Hex sampling combines fast 3D encoding with SNR efficiency and favorable depiction characteristics regarding noise amplification and filtering effects from T2∗ decay, thereby providing flexibility in the choice of imaging parameters. It is attractive both for high‐resolution time series such as fMRI and for applications that require rapid anatomical imaging.

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.

Diffusion Tensor Imaging-Based Studies at the Group-Level Applied to Animal Models of Neurodegenerative Diseases

The understanding of human and non-human microstructural brain alterations in the course of neurodegenerative diseases has substantially improved by the non-invasive magnetic resonance imaging (MRI) technique of diffusion tensor imaging (DTI). Animal models (including disease or knockout models) allow for a variety of experimental manipulations, which are not applicable to humans. Thus, the DTI approach provides a promising tool for cross-species cross-sectional and longitudinal investigations of the neurobiological targets and mechanisms of neurodegeneration. This overview with a systematic review focuses on the principles of DTI analysis as used in studies at the group level in living preclinical models of neurodegeneration. The translational aspect from in-vivo animal models toward (clinical) applications in humans is covered as well as the DTI-based research of the non-human brains' microstructure, the methodological aspects in data processing and analysis, and data interpretation at different abstraction levels. The aim of integrating DTI in multiparametric or multimodal imaging protocols will allow the interrogation of DTI data in terms of directional flow of information and may identify the microstructural underpinnings of neurodegeneration-related patterns.

The quest for high spatial resolution diffusion-weighted imaging of the human brain in vivo

Diffusion-weighted imaging, a contrast unique to MRI, is used for assessment of tissue microstructure in vivo. However, this exquisite sensitivity to finer scales far above imaging resolution comes at the cost of vulnerability to errors caused by sources of motion other than diffusion motion. Addressing the issue of motion has traditionally limited diffusion-weighted imaging to a few acquisition techniques and, as a consequence, to poorer spatial resolution than other MRI applications. Advances in MRI imaging methodology have allowed diffusion-weighted MRI to push to ever higher spatial resolution. In this review we focus on the pulse sequences and associated techniques under development that have pushed the limits of image quality and spatial resolution in diffusion-weighted MRI.

Simultaneous and inherent correction of B0 and eddy-current induced distortions in high-resolution diffusion MRI using reversed polarity gradients and multiplexed sensitivity encoding (RPG-MUSE)

In diffusion MRI (dMRI), static magnetic field (B0) inhomogeneity and time varying gradient eddy currents induce spatial distortions in reconstructed images. These distortions are exacerbated when high spatial resolutions are used, and many field-mapping based correction techniques often only acquire maps of static B0 distortion, which are not adequate for correcting eddy current induced image distortions. This report presents a novel technique, termed RPG-MUSE, for achieving distortion-free high-resolution diffusion MRI by integrating reversed polarity gradients (RPG) into the multi-shot echo planar imaging acquisition scheme used in multiplexed sensitivity encoding (MUSE). By alternating the phase encoding direction between shots in both baseline and diffusion-weighted acquisitions, maps of both static B0 and eddy current induced field inhomogeneities can be inherently derived, without the need for additional data acquisition. Through both 2D and 3D encoded dMRI acquisitions, it is shown that an RPG-MUSE reconstruction can simultaneously achieve high spatial resolution, high spatial fidelity, and subsequently, high accuracy in diffusion metrics.

Exploration of male urethral sphincter complex using diffusion tensor imaging (DTI)-based fiber-tracking

BACKGROUND

Urinary incontinence is a major clinical problem arising primarily from age-related degenerative changes to the sphincter muscles. However, the precise anatomy of the normal male sphincter muscles has yet to be established. Diffusion tensor imaging (DTI) may offer a unique insight into muscle microstructure and fiber architecture.

PURPOSE

To explore the anatomy of the urethral sphincter muscles pertinent to urinary continence function using DT-MRI.

STUDY TYPE

Prospective cohort study.

SUBJECTS

Eleven normal male subjects (mean age: 25.4 years); two subjects were scanned in three separate sessions to assess reproducibility.

FIELD STRENGTH/SEQUENCE

3T; using a diffusion-weighted spin echo planar sequence.

ASSESSMENT

DT parameters including fractional anisotropy (FA), primary (λ1 ), secondary (λ2 ), and tertiary (λ3 ) eigenvalues, Apparent diffusion coefficient and radial diffusivity were analyzed statistically, while tracked muscle fibers were assessed visually.

STATISTICAL TESTS

Regional differences (sphincters and longitudinal muscle of the urethra) in the DTI indices were assessed by one-way analysis of variance. A Tukey post-hoc test was used to identify significant differences between muscle regions.

RESULTS

Two sphincter muscles, one proximal near the base of the bladder, corresponding to the lisso-sphincter, and the other distal to the end of the prostate corresponding to the rhabdo-sphincter, surrounding a central urethral muscle fiber bundle, were clearly identified. FA was higher and λ3 lower in the proximal sphincter muscle compared to the central urethral muscle and the distal sphincter (P < 0.05). The average coefficient of variation ranged from 5-12% for the DTI indices.

DATA CONCLUSION

Since DTI values are known to reflect underlying tissue microarchitecture, significant differences in DTI indices identified here between the muscles of the urethral complex may potentially arise from differences in tissue microarchitecture that may in turn be related to the specific function of the sphincter and other muscles.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018;48:1002-1011.

A diffusion-matched principal component analysis (DM-PCA) based two-channel denoising procedure for high-resolution diffusion-weighted MRI

Over the past several years, significant efforts have been made to improve the spatial resolution of diffusion-weighted imaging (DWI), aiming at better detecting subtle lesions and more reliably resolving white-matter fiber tracts. A major concern with high-resolution DWI is the limited signal-to-noise ratio (SNR), which may significantly offset the advantages of high spatial resolution. Although the SNR of DWI data can be improved by denoising in post-processing, existing denoising procedures may potentially reduce the anatomic resolvability of high-resolution imaging data. Additionally, non-Gaussian noise induced signal bias in low-SNR DWI data may not always be corrected with existing denoising approaches. Here we report an improved denoising procedure, termed diffusion-matched principal component analysis (DM-PCA), which comprises 1) identifying a group of (not necessarily neighboring) voxels that demonstrate very similar magnitude signal variation patterns along the diffusion dimension, 2) correcting low-frequency phase variations in complex-valued DWI data, 3) performing PCA along the diffusion dimension for real- and imaginary-components (in two separate channels) of phase-corrected DWI voxels with matched diffusion properties, 4) suppressing the noisy PCA components in real- and imaginary-components, separately, of phase-corrected DWI data, and 5) combining real- and imaginary-components of denoised DWI data. Our data show that the new two-channel (i.e., for real- and imaginary-components) DM-PCA denoising procedure performs reliably without noticeably compromising anatomic resolvability. Non-Gaussian noise induced signal bias could also be reduced with the new denoising method. The DM-PCA based denoising procedure should prove highly valuable for high-resolution DWI studies in research and clinical uses.

A 2D spiral turbo-spin-echo technique

PURPOSE

2D turbo-spin-echo (TSE) is widely used in the clinic for neuroimaging. However, the long refocusing radiofrequency pulse train leads to high specific absorption rate (SAR) and alters the contrast compared to conventional spin-echo. The purpose of this work is to develop a robust 2D spiral TSE technique for fast T₂ -weighted imaging with low SAR and improved contrast.

METHODS

A spiral-in/out readout is incorporated into 2D TSE to fully take advantage of the acquisition efficiency of spiral sampling while avoiding potential off-resonance-related artifacts compared to a typical spiral-out readout. A double encoding strategy and a signal demodulation method are proposed to mitigate the artifacts because of the T₂ -decay-induced signal variation. An adapted prescan phase correction as well as a concomitant phase compensation technique are implemented to minimize the phase errors.

RESULTS

Phantom data demonstrate the efficacy of the proposed double encoding/signal demodulation, as well as the prescan phase correction and concomitant phase compensation. Volunteer data show that the proposed 2D spiral TSE achieves fast scan speed with high SNR, low SAR, and improved contrast compared to conventional Cartesian TSE.

CONCLUSION

A robust 2D spiral TSE technique is feasible and provides a potential alternative to conventional 2D Cartesian TSE for T₂ -weighted neuroimaging.

3D-MB-MUSE: A robust 3D multi-slab, multi-band and multi-shot reconstruction approach for ultrahigh resolution diffusion MRI

Recent advances in achieving ultrahigh spatial resolution (e.g. sub-millimeter) diffusion MRI (dMRI) data have proven highly beneficial in characterizing tissue microstructures in organs such as the brain. However, the routine acquisition of in-vivo dMRI data at such high spatial resolutions has been largely prohibited by factors that include prolonged acquisition times, motion induced artifacts, and low SNR. To overcome these limitations, we present here a framework for acquiring and reconstructing 3D multi-slab, multi-band and interleaved multi-shot EPI data, termed 3D-MB-MUSE. Through multi-band excitations, the simultaneous acquisition of multiple 3D slabs enables whole brain dMRI volumes to be acquired in-vivo on a 3 T clinical MRI scanner at high spatial resolution within a reasonably short amount of time. Representing a true 3D model, 3D-MB-MUSE reconstructs an entire 3D multi-band, multi-shot dMRI slab at once while simultaneously accounting for coil sensitivity variations across the slab as well as motion induced artifacts commonly associated with both 3D and multi-shot diffusion imaging. Such a reconstruction fully preserves the SNR advantages of both 3D and multi-shot acquisitions in high resolution dMRI images by removing both motion and aliasing artifacts across multiple dimensions. By enabling ultrahigh resolution dMRI for routine use, the 3D-MB-MUSE framework presented here may prove highly valuable in both clinical and research applications.

Image formation in diffusion MRI: A review of recent technical developments

Diffusion magnetic resonance imaging (MRI) is a standard imaging tool in clinical neurology, and is becoming increasingly important for neuroscience studies due to its ability to depict complex neuroanatomy (eg, white matter connectivity). Single-shot echo-planar imaging is currently the predominant formation method for diffusion MRI, but suffers from blurring, distortion, and low spatial resolution. A number of methods have been proposed to address these limitations and improve diffusion MRI acquisition. Here, the recent technical developments for image formation in diffusion MRI are reviewed. We discuss three areas of advance in diffusion MRI: improving image fidelity, accelerating acquisition, and increasing the signal-to-noise ratio.

LEVEL OF EVIDENCE

5 Technical Efficacy: Stage 1 J. MAGN. RESON. IMAGING 2017;46:646-662.

Reducing slab boundary artifacts in three-dimensional multislab diffusion MRI using nonlinear inversion for slab profile encoding (NPEN)

PURPOSE

To propose a method to reduce the slab boundary artifacts in three-dimensional multislab diffusion MRI.

METHODS

Bloch simulation is used to investigate the effects of multiple factors on slab boundary artifacts, including characterization of residual errors on diffusion quantification. A nonlinear inversion method is proposed to simultaneously estimate the slab profile and the underlying (corrected) image.

RESULTS

Correction results of numerical phantom and in vivo data demonstrate that the method can effectively remove slab boundary artifacts for diffusion data. Notably, the nonlinear inversion is also successful at short TR, a regimen where previously proposed methods (slab profile encoding and weighted average) retain residual artifacts in both diffusion-weighted images and diffusion metrics (mean diffusion coefficient and fractional anisotropy).

CONCLUSION

The nonlinear inversion for removing slab boundary artifacts provides improvements over existing methods, particularly at the short TRs required to maximize SNR efficiency. Magn Reson Med 76:1183-1195, 2016. © 2015 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.

High-resolution diffusion MRI at 7T using a three-dimensional multi-slab acquisition

High-resolution diffusion MRI can provide the ability to resolve small brain structures, enabling investigations of detailed white matter architecture. A major challenge for in vivo high-resolution diffusion MRI is the low signal-to-noise ratio. In this work, we combine two highly compatible methods, ultra-high field and three-dimensional multi-slab acquisition to improve the SNR of high-resolution diffusion MRI. As each k_z plane is encoded using a single-shot echo planar readout, scan speeds of the proposed technique are similar to the commonly used two-dimensional diffusion MRI. In-plane parallel acceleration is applied to reduce image distortions. To reduce the sensitivity of auto-calibration signal data to subject motion and respiration, several new adaptions of the fast low angle excitation echo-planar technique (FLEET) that are suitable for 3D multi-slab echo planar imaging are proposed and evaluated. A modified reconstruction scheme is proposed for auto-calibration with the most robust method, Slice-FLEET acquisition, to make it compatible with navigator correction of motion induced phase errors. Slab boundary artefacts are corrected using the nonlinear slab profile encoding method recently proposed by our group. In vivo results demonstrate that using 7T and three-dimensional multi-slab acquisition with improved auto-calibration signal acquisition and nonlinear slab boundary artefacts correction, high-quality diffusion MRI data with ~1mm isotropic resolution can be achieved.

Joint correction of Nyquist artifact and minuscule motion-induced aliasing artifact in interleaved diffusion weighted EPI data using a composite two-dimensional phase correction procedure

Diffusion-weighted imaging (DWI) obtained with interleaved echo-planar imaging (EPI) pulse sequence has great potential of characterizing brain tissue properties at high spatial-resolution. However, interleaved EPI based DWI data may be corrupted by various types of aliasing artifacts. First, inconsistencies in k-space data obtained with opposite readout gradient polarities result in Nyquist artifact, which is usually reduced with 1D phase correction in post-processing. When there exist eddy current cross terms (e.g., in oblique-plane EPI), 2D phase correction is needed to effectively reduce Nyquist artifact. Second, minuscule motion induced phase inconsistencies in interleaved DWI scans result in image-domain aliasing artifact, which can be removed with reconstruction procedures that take shot-to-shot phase variations into consideration. In existing interleaved DWI reconstruction procedures, Nyquist artifact and minuscule motion-induced aliasing artifact are typically removed subsequently in two stages. Although the two-stage phase correction generally performs well for non-oblique plane EPI data obtained from well-calibrated system, the residual artifacts may still be pronounced in oblique-plane EPI data or when there exist eddy current cross terms. To address this challenge, here we report a new composite 2D phase correction procedure, which effective removes Nyquist artifact and minuscule motion induced aliasing artifact jointly in a single step. Our experimental results demonstrate that the new 2D phase correction method can much more effectively reduce artifacts in interleaved EPI based DWI data as compared with the existing two-stage artifact correction procedures. The new method robustly enables high-resolution DWI, and should prove highly valuable for clinical uses and research studies of DWI.

High-resolution multishot spiral diffusion tensor imaging with inherent correction of motion-induced phase errors

PURPOSE

To develop and compare three novel reconstruction methods designed to inherently correct for motion-induced phase errors in multishot spiral diffusion tensor imaging without requiring a variable-density spiral trajectory or a navigator echo.

THEORY AND METHODS

The first method simply averages magnitude images reconstructed with sensitivity encoding from each shot, whereas the second and third methods rely on sensitivity encoding to estimate the motion-induced phase error for each shot and subsequently use either a direct phase subtraction or an iterative conjugate gradient algorithm, respectively, to correct for the resulting artifacts. Numerical simulations and in vivo experiments on healthy volunteers were performed to assess the performance of these methods.

RESULTS

The first two methods suffer from a low signal-to-noise ratio or from residual artifacts in the reconstructed diffusion-weighted images and fractional anisotropy maps. In contrast, the third method provides high-quality, high-resolution diffusion tensor imaging results, revealing fine anatomical details such as a radial diffusion anisotropy in cortical gray matter.

CONCLUSION

The proposed sensitivity encoding + conjugate gradient method can inherently and effectively correct for phase errors, signal loss, and aliasing artifacts caused by both rigid and nonrigid motion in multishot spiral diffusion tensor imaging, without increasing the scan time or reducing the signal-to-noise ratio.

PSF mapping-based correction of eddy-current-induced distortions in diffusion-weighted echo-planar imaging

PURPOSE

To accurately correct diffusion-encoding direction-dependent eddy-current-induced geometric distortions in diffusion-weighted echo-planar imaging (DW-EPI) and to minimize the calibration time at 7 Tesla (T).

METHODS

A point spread function (PSF) mapping based eddy-current calibration method is newly presented to determine eddy-current-induced geometric distortions even including nonlinear eddy-current effects within the readout acquisition window. To evaluate the temporal stability of eddy-current maps, calibration was performed four times within 3 months. Furthermore, spatial variations of measured eddy-current maps versus their linear superposition were investigated to enable correction in DW-EPIs with arbitrary diffusion directions without direct calibration. For comparison, an image-based eddy-current correction method was additionally applied. Finally, this method was combined with a PSF-based susceptibility-induced distortion correction approach proposed previously to correct both susceptibility and eddy-current-induced distortions in DW-EPIs.

RESULTS

Very fast eddy-current calibration in a three-dimensional volume is possible with the proposed method. The measured eddy-current maps are very stable over time and very similar maps can be obtained by linear superposition of principal-axes eddy-current maps. High resolution in vivo brain results demonstrate that the proposed method allows more efficient eddy-current correction than the image-based method.

CONCLUSION

The combination of both PSF-based approaches allows distortion-free images, which permit reliable analysis in diffusion tensor imaging applications at 7T.

Human brain diffusion tensor imaging at submillimeter isotropic resolution on a 3Tesla clinical MRI scanner

The advantages of high-resolution diffusion tensor imaging (DTI) have been demonstrated in a recent post-mortem human brain study (Miller et al., NeuroImage 2011;57(1):167-181), showing that white matter fiber tracts can be much more accurately detected in data at a submillimeter isotropic resolution. To our knowledge, in vivo human brain DTI at a submillimeter isotropic resolution has not been routinely achieved yet because of the difficulty in simultaneously achieving high resolution and high signal-to-noise ratio (SNR) in DTI scans. Here we report a 3D multi-slab interleaved EPI acquisition integrated with multiplexed sensitivity encoded (MUSE) reconstruction, to achieve high-quality, high-SNR and submillimeter isotropic resolution (0.85×0.85×0.85mm(3)) in vivo human brain DTI on a 3Tesla clinical MRI scanner. In agreement with the previously reported post-mortem human brain DTI study, our in vivo data show that the structural connectivity networks of human brains can be mapped more accurately and completely with high-resolution DTI as compared with conventional DTI (e.g., 2×2×2mm(3)).

Sliding-slab three-dimensional TSE imaging with a spiral-In/Out readout

PURPOSE

T2 -weighted imaging is of great diagnostic value in neuroimaging. Three-dimensional (3D) Cartesian turbo spin echo (TSE) scans provide high signal-to-noise ratio (SNR) and contiguous slice coverage. The purpose of this preliminary work is to implement a novel 3D spiral TSE technique with image quality comparable to 2D/3D Cartesian TSE.

METHODS

The proposed technique uses multislab 3D TSE imaging. To mitigate the slice boundary artifacts, a sliding-slab method is extended to spiral imaging. A spiral-in/out readout is adopted to minimize the artifacts that may be present with the conventional spiral-out readout. Phase errors induced by B0 eddy currents are measured and compensated to allow for the combination of the spiral-in and spiral-out images. A nonuniform slice encoding scheme is used to reduce the truncation artifacts while preserving the SNR performance.

RESULTS

Preliminary results show that each of the individual measures contributes to the overall performance, and the image quality of the results obtained with the proposed technique is, in general, comparable to that of 2D or 3D Cartesian TSE.

CONCLUSION

3D sliding-slab TSE with a spiral-in/out readout provides good-quality T2 -weighted images, and, therefore, may become a promising alternative to Cartesian TSE.

POCS-based reconstruction of multiplexed sensitivity encoded MRI (POCSMUSE): A general algorithm for reducing motion-related artifacts

PURPOSE

A projection onto convex sets reconstruction of multiplexed sensitivity encoded MRI (POCSMUSE) is developed to reduce motion-related artifacts, including respiration artifacts in abdominal imaging and aliasing artifacts in interleaved diffusion-weighted imaging.

THEORY

Images with reduced artifacts are reconstructed with an iterative projection onto convex sets (POCS) procedure that uses the coil sensitivity profile as a constraint. This method can be applied to data obtained with different pulse sequences and k-space trajectories. In addition, various constraints can be incorporated to stabilize the reconstruction of ill-conditioned matrices.

METHODS

The POCSMUSE technique was applied to abdominal fast spin-echo imaging data, and its effectiveness in respiratory-triggered scans was evaluated. The POCSMUSE method was also applied to reduce aliasing artifacts due to shot-to-shot phase variations in interleaved diffusion-weighted imaging data corresponding to different k-space trajectories and matrix condition numbers.

RESULTS

Experimental results show that the POCSMUSE technique can effectively reduce motion-related artifacts in data obtained with different pulse sequences, k-space trajectories and contrasts.

CONCLUSION

POCSMUSE is a general post-processing algorithm for reduction of motion-related artifacts. It is compatible with different pulse sequences, and can also be used to further reduce residual artifacts in data produced by existing motion artifact reduction methods.

Improved diffusion tensor imaging of the optic nerve using multishot two-dimensional navigated acquisitions

PURPOSE

A diffusion-weighted multishot echo-planar imaging approach combined with SENSE and a two-dimensional (2D) navigated motion correction was investigated as an alternative to conventional single-shot counterpart to obtain optic nerve images at higher spatial resolution with reduced artifacts.

METHODS

Fifteen healthy subjects were enrolled in the study. Six of these subjects underwent a repeated acquisition at least 2 weeks after the initial scan session to address reproducibility. Both single-shot and multishot diffusion tensor imaging studies of the human optic nerve were performed with matched scan time. Effect of subject motions were corrected using 2D phase navigator during multishot image reconstruction. Tensor-derived indices from proposed multishot were compared against conventional single-shot approach. Image resolution difference, right-left optic nerve asymmetry, and test-retest reproducibility were also assessed.

RESULTS

In vivo results of acquired multishot images and quantitative maps of diffusion properties of the optic nerve showed significantly reduced image artifacts (e.g., distortions and blurring), and the derived diffusion indices were comparable to those from other studies. Single-shot scans presented larger variability between right and left optic nerves than multishot scans. Multishot scans also presented smaller variations across scans at different time points when compared with single-shot counterparts.

CONCLUSION

The multishot technique has considerable potential for providing improved information on optic nerve pathology and may also be translated to higher fields.

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.

Cortical depth dependence of the diffusion anisotropy in the human cortical gray matter in vivo

Diffusion tensor imaging (DTI) is typically used to study white matter fiber pathways, but may also be valuable to assess the microstructure of cortical gray matter. Although cortical diffusion anisotropy has previously been observed in vivo, its cortical depth dependence has mostly been examined in high-resolution ex vivo studies. This study thus aims to investigate the cortical depth dependence of the diffusion anisotropy in the human cortex in vivo on a clinical 3 T scanner. Specifically, a novel multishot constant-density spiral DTI technique with inherent correction of motion-induced phase errors was used to achieve a high spatial resolution (0.625×0.625×3 mm) and high spatial fidelity with no scan time penalty. The results show: (i) a diffusion anisotropy in the cortical gray matter, with a primarily radial diffusion orientation, as observed in previous ex vivo and in vivo studies, and (ii) a cortical depth dependence of the fractional anisotropy, with consistently higher values in the middle cortical lamina than in the deep and superficial cortical laminae, as observed in previous ex vivo studies. These results, which are consistent across subjects, demonstrate the feasibility of this technique for investigating the cortical depth dependence of the diffusion anisotropy in the human cortex in vivo.

High-resolution reduced field of view diffusion tensor imaging using spatially selective RF pulses

PURPOSE

Diffusion tensor imaging (DTI) plays a vital role in identifying white matter fiber bundles. Achievable imaging resolution and imaging time demands remain the major challenges in detecting small fiber bundles with current clinical DTI sequences.

METHODS

A novel reduced field of view ultra-high-resolution DTI technique named eZOOM (elliptically refocused zonally oblique multislice) was developed. A small circular disk was imaged using spatially selective radiofrequency (RF) pulses, reducing the imaging matrix size. The frequency profile of the spectral-spatial refocusing RF pulse provided intrinsic fat suppression, eliminating the need for fat saturation pulses.

RESULTS

Multislice DTI at a resolution of 0.35 × 0.35 mm in a celery fiber phantom was successfully performed by scanning an 8-cm field of view at 3T. An adequate diffusion-to-noise ratio (DNR >20) was achieved for a 25-min acquisition using a direct-sampling RF receiver. Human subjects (n = 7) were scanned at resolutions of 0.47 × 0.47 mm having a DNR <20 within a 75-min scanning time, requiring further enhancements to increase the signal-to-noise ratio.

CONCLUSIONS

The new eZOOM-DTI method offers multislice DTI at ultra-high imaging resolutions substantially exceeding those available with current echo-planar DTI techniques. Parallel and fast spin echo methods can be combined with eZOOM to improve SNR and DNR in humans.

Accelerated multi-shot diffusion imaging

PURPOSE

To reduce image distortion in MR diffusion imaging using an accelerated multi-shot method.

METHODS

The proposed method exploits the fact that diffusion-encoded data tend to be sparse when represented in the kb-kd space, where kb and kd are the Fourier transform duals of b and d, the b-factor and the diffusion direction, respectively. Aliasing artifacts are displaced toward under-used regions of the kb-kd plane, allowing nonaliased signals to be recovered. A main characteristic of the proposed approach is how thoroughly the navigator information gets used during reconstruction: The phase of navigator images is used for motion correction, while the magnitude of the navigator signal in kb-kd space is used for regularization purposes. As opposed to most acceleration methods based on compressed sensing, the proposed method reduces the number of ky lines needed for each diffusion-encoded image, but not the total number of images required. Consequently, it tends to be most effective at reducing image distortion rather than reducing total scan time.

RESULTS

Results are presented for three volunteers with acceleration factors ranging from 4 to 8, with and without the inclusion of parallel imaging.

CONCLUSION

An accelerated motion-corrected diffusion imaging method was introduced that achieves good image quality at relatively high acceleration factors.

Myofiber prestretch magnitude determines regional systolic function during ectopic activation in the tachycardia-induced failing canine heart

Electrical dyssynchrony leads to prestretch in late-activated regions and alters the sequence of mechanical contraction, although prestretch and its mechanisms are not well defined in the failing heart. We hypothesized that in heart failure, fiber prestretch magnitude increases with the amount of early-activated tissue and results in increased end-systolic strains, possibly due to length-dependent muscle properties. In five failing dog hearts with scars, three-dimensional strains were measured at the anterolateral left ventricle (LV). Prestretch magnitude was varied via ventricular pacing at increasing distances from the measurement site and was found to increase with activation time at various wall depths. At the subepicardium, prestretch magnitude positively correlated with the amount of early-activated tissue. At the subendocardium, local end-systolic strains (fiber shortening, radial wall thickening) increased proportionally to prestretch magnitude, resulting in greater mean strain values in late-activated compared with early-activated tissue. Increased fiber strains at end systole were accompanied by increases in preejection fiber strain, shortening duration, and the onset of fiber relengthening, which were all positively correlated with local activation time. In a dog-specific computational failing heart model, removal of length and velocity dependence on active fiber stress generation, both separately and together, alter the correlations between local electrical activation time and timing of fiber strains but do not primarily account for these relationships.

Patient-Specific Models of Cardiac Biomechanics

Patient-specific models of cardiac function have the potential to improve diagnosis and management of heart disease by integrating medical images with heterogeneous clinical measurements subject to constraints imposed by physical first principles and prior experimental knowledge. We describe new methods for creating three-dimensional patient-specific models of ventricular biomechanics in the failing heart. Three-dimensional bi-ventricular geometry is segmented from cardiac CT images at end-diastole from patients with heart failure. Human myofiber and sheet architecture is modeled using eigenvectors computed from diffusion tensor MR images from an isolated, fixed human organ-donor heart and transformed to the patient-specific geometric model using large deformation diffeomorphic mapping. Semi-automated methods were developed for optimizing the passive material properties while simultaneously computing the unloaded reference geometry of the ventricles for stress analysis. Material properties of active cardiac muscle contraction were optimized to match ventricular pressures measured by cardiac catheterization, and parameters of a lumped-parameter closed-loop model of the circulation were estimated with a circulatory adaptation algorithm making use of information derived from echocardiography. These components were then integrated to create a multi-scale model of the patient-specific heart. These methods were tested in five heart failure patients from the San Diego Veteran's Affairs Medical Center who gave informed consent. The simulation results showed good agreement with measured echocardiographic and global functional parameters such as ejection fraction and peak cavity pressures.

A robust multi-shot scan strategy for high-resolution diffusion weighted MRI enabled by multiplexed sensitivity-encoding (MUSE)

Diffusion weighted magnetic resonance imaging (DWI) data have been mostly acquired with single-shot echo-planar imaging (EPI) to minimize motion induced artifacts. The spatial resolution, however, is inherently limited in single-shot EPI, even when the parallel imaging (usually at an acceleration factor of 2) is incorporated. Multi-shot acquisition strategies could potentially achieve higher spatial resolution and fidelity, but they are generally susceptible to motion-induced phase errors among excitations that are exacerbated by diffusion sensitizing gradients, rendering the reconstructed images unusable. It has been shown that shot-to-shot phase variations may be corrected using navigator echoes, but at the cost of imaging throughput. To address these challenges, a novel and robust multi-shot DWI technique, termed multiplexed sensitivity-encoding (MUSE), is developed here to reliably and inherently correct nonlinear shot-to-shot phase variations without the use of navigator echoes. The performance of the MUSE technique is confirmed experimentally in healthy adult volunteers on 3Tesla MRI systems. This newly developed technique should prove highly valuable for mapping brain structures and connectivities at high spatial resolution for neuroscience studies.

An atlas-based geometry pipeline for cardiac Hermite model construction and diffusion tensor reorientation

Here we present a novel atlas-based geometry pipeline for constructing three-dimensional cubic Hermite finite element meshes of the whole human heart from tomographic patient image data. To build the cardiac atlas, two superior atria, two inferior ventricles as well as the aorta and the pulmonary trunk are first segmented, and epicardial and endocardial boundary surfaces are extracted and smoothed. Critical points and skeletons (or central-line paths) are identified, following the cardiac topology. The surface model and the path tree are used to construct a hexahedral control mesh via a skeleton-based sweeping method. Derivative parameters are computed from the control mesh, defining cubic Hermite finite elements. The thickness of the atria and the ventricles is obtained using segmented epicardial boundaries or via offsetting from the endocardial surfaces in regions where the image resolution is insufficient. We also develop a robust optical flow approach to deform the constructed atlas and align it with the image from a second patient. This registration method is fully-automatic, and avoids manual operations required by segmentation and path extraction. Moreover, we demonstrate that this method can also be used to deformably map diffusion tensor MRI data with patient geometries to include fiber and sheet orientations in the finite element model.

Inherent correction of motion-induced phase errors in multishot spiral diffusion-weighted imaging

Multishot spiral imaging is a promising alternative to echo-planar imaging for high-resolution diffusion-weighted imaging and diffusion tensor imaging. However, subject motion in the presence of diffusion-weighting gradients causes phase inconsistencies among different shots, resulting in signal loss and aliasing artifacts in the reconstructed images. Such artifacts can be reduced using a variable-density spiral trajectory or a navigator echo, however at the cost of a longer scan time. Here, a novel iterative phase correction method is proposed to inherently correct for the motion-induced phase errors without requiring any additional scan time. In this initial study, numerical simulations and in vivo experiments are performed to demonstrate that the proposed method can effectively and efficiently correct for spatially linear phase errors caused by rigid-body motion in multishot spiral diffusion-weighted imaging of the human brain.

Motion-induced phase error estimation and correction in 3D diffusion tensor imaging

A multishot data acquisition strategy is one way to mitigate B0 distortion and T2∗ blurring for high-resolution diffusion-weighted magnetic resonance imaging experiments. However, different object motions that take place during different shots cause phase inconsistencies in the data, leading to significant image artifacts. This work proposes a maximum likelihood estimation and k-space correction of motion-induced phase errors in 3D multishot diffusion tensor imaging. The proposed error estimation is robust, unbiased, and approaches the Cramer-Rao lower bound. For rigid body motion, the proposed correction effectively removes motion-induced phase errors regardless of the k-space trajectory used and gives comparable performance to the more computationally expensive 3D iterative nonlinear phase error correction method. The method has been extended to handle multichannel data collected using phased-array coils. Simulation and in vivo data are shown to demonstrate the performance of the method.

Dynamic correction of artifacts due to susceptibility effects and time-varying eddy currents in diffusion tensor imaging

In diffusion tensor imaging (DTI), spatial and temporal variations of the static magnetic field (B(0)) caused by susceptibility effects and time-varying eddy currents result in severe distortions, blurring, and misregistration artifacts, which in turn lead to errors in DTI metrics and in fiber tractography. Various correction methods have been proposed, but typically assume that the eddy current-induced magnetic field can be modeled as a constant or a single exponential decay within the DTI readout window. Here, we show that its temporal dependence is more complex because of the interaction of multiple eddy currents with different time constants, but that it remains very consistent over time. As such, we propose a novel dynamic B(0) mapping and off-resonance correction method that measures the exact spatial, temporal, and diffusion-weighting direction dependence of the susceptibility- and eddy current-induced magnetic fields to effectively and efficiently correct for artifacts caused by both susceptibility effects and time-varying eddy currents, thereby resulting in a high spatial fidelity and accuracy.

Head injury or head motion? Assessment and quantification of motion artifacts in diffusion tensor imaging studies

The relationship between head motion and diffusion values such as fractional anisotropy (FA) and mean diffusivity (MD) is currently not well understood. Simulation studies suggest that head motion may introduce either a positive or negative bias, but this has not been quantified in clinical studies. Moreover, alternative measures for removing bias as result of head motion, such as the removal of problematic gradients, has been suggested but not carefully evaluated. The current study examined the impact of head motion on FA and MD across three common pipelines (tract-based spatial statistics, voxelwise, and region of interest analyses) and determined the impact of removing diffusion weighted images. Our findings from a large cohort of healthy controls indicate that while head motion was associated with a positive bias for both FA and MD, the effect was greater for MD. The positive bias was observed across all three analysis pipelines and was present following established protocols for data processing, suggesting that current techniques (i.e., correction of both image and gradient table) for removing motion bias are likely insufficient. However, the removal of images with gross artifacts did not fundamentally change the relationship between motion and DTI scalar values. In addition, Monte Carlo simulations suggested that the random removal of images increases the bias and reduces the precision of both FA and MD. Finally, we provide an example of how head motion can be quantified across different neuropsychiatric populations, which should be implemented as part of any diffusion tensor imaging quality assurance protocol.

Use of High Resolution 3D Diffusion Tensor Imaging to Study Brain White Matter Development in Live Neonatal Rats

High resolution diffusion tensor imaging (DTI) can provide important information on brain development, yet it is challenging in live neonatal rats due to the small size of neonatal brain and motion-sensitive nature of DTI. Imaging in live neonatal rats has clear advantages over fixed brain scans, as longitudinal and functional studies would be feasible to understand neuro-developmental abnormalities. In this study, we developed imaging strategies that can be used to obtain high resolution 3D DTI images in live neonatal rats at postnatal day 5 (PND5) and PND14, using only 3 h of imaging acquisition time. An optimized 3D DTI pulse sequence and appropriate animal setup to minimize physiological motion artifacts are the keys to successful high resolution 3D DTI imaging. Thus, a 3D rapid acquisition relaxation enhancement DTI sequence with twin navigator echoes was implemented to accelerate imaging acquisition time and minimize motion artifacts. It has been suggested that neonatal mammals possess a unique ability to tolerate mild-to-moderate hypothermia and hypoxia without long term impact. Thus, we additionally utilized this ability to minimize motion artifacts in magnetic resonance images by carefully suppressing the respiratory rate to around 15/min for PND5 and 30/min for PND14 using mild-to-moderate hypothermia. These imaging strategies have been successfully implemented to study how the effect of cocaine exposure in dams might affect brain development in their rat pups. Image quality resulting from this in vivo DTI study was comparable to ex vivo scans. fractional anisotropy values were also similar between the live and fixed brain scans. The capability of acquiring high quality in vivo DTI imaging offers a valuable opportunity to study many neurological disorders in brain development in an authentic living environment.