Parallel imaging performance as a function of field strength--an experimental investigation using electrodynamic scaling

Lipid suppression via double inversion recovery with symmetric frequency sweep for robust 2D-GRAPPA-accelerated MRSI of the brain at 7 T

This work presents a new approach for high-resolution MRSI of the brain at 7 T in clinically feasible measurement times. Two major problems of MRSI are the long scan times for large matrix sizes and the possible spectral contamination by the transcranial lipid signal. We propose a combination of free induction decay (FID)-MRSI with a short acquisition delay and acceleration via in-plane two-dimensional generalised autocalibrating partially parallel acquisition (2D-GRAPPA) with adiabatic double inversion recovery (IR)-based lipid suppression to allow robust high-resolution MRSI. We performed Bloch simulations to evaluate the magnetisation pathways of lipids and metabolites, and compared the results with phantom measurements. Acceleration factors in the range 2-25 were tested in a phantom. Five volunteers were scanned to verify the value of our MRSI method in vivo. GRAPPA artefacts that cause fold-in of transcranial lipids were suppressed via double IR, with a non-selective symmetric frequency sweep. The use of long, low-power inversion pulses (100 ms) reduced specific absorption rate requirements. The symmetric frequency sweep over both pulses provided good lipid suppression (>90%), in addition to a reduced loss in metabolite signal-to-noise ratio (SNR), compared with conventional IR suppression (52-70%). The metabolic mapping over the whole brain slice was not limited to a rectangular region of interest. 2D-GRAPPA provided acceleration up to a factor of nine for in vivo FID-MRSI without a substantial increase in g-factors (<1.1). A 64 × 64 matrix can be acquired with a common repetition time of ~1.3 s in only 8 min without lipid artefacts caused by acceleration. Overall, we present a fast and robust MRSI method, using combined double IR fat suppression and 2D-GRAPPA acceleration, which may be used in (pre)clinical studies of the brain at 7 T.

Clinically feasible NODDI characterization of glioma using multiband EPI at 7 T

Recent technological progress in the multiband echo planer imaging (MB EPI) technique enables accelerated MR diffusion weighted imaging (DWI) and allows whole brain, multi-b-value diffusion imaging to be acquired within a clinically feasible time. However, its applications at 7 T have been limited due to B1 field inhomogeneity and increased susceptibility artifact. It is an ongoing debate whether DWI at 7 T can be performed properly in patients, and a systematic SNR comparison for multiband spin-echo EPI between 3 T and 7 T has not been methodically studied. The goal of this study was to use MB EPI at 7 T in order to obtain 90-directional multi-shell DWI within a clinically feasible acquisition time for patients with glioma. This study included an SNR comparison between 3 T and 7 T, and the application of B1 mapping and distortion correction procedures for reducing the impact of variations in B0 and B1. The optimized multiband sequence was applied in 20 patients with glioma to generate both DTI and NODDI maps for comparison of values in tumor and normal appearing white matter (NAWM). Our SNR analysis showed that MB EPI at 7 T was comparable to that at 3 T, and the data quality acquired in patients was clinically acceptable. NODDI maps provided unique contrast within the T2 lesion that was not seen in anatomical images or DTI maps. Such contrast may reflect the complexity of tissue compositions associated with disease progression and treatment effects. The ability to consistently obtain high quality diffusion data at 7 T will contribute towards the implementation of a comprehensive brain MRI examination at ultra-high field.

•

NODDI characterization of glioma at 7 T within 6 min with multiband EPI

•

An SNR comparison was performed between 7 T and 3 T SE-EPI.

•

SNR was comparable between 7 T and 3 T multiband SE-EPI.

•

NODDI maps provided unique contrast within the glioma T2 lesion.

Rapid multi-orientation quantitative susceptibility mapping

Three-dimensional gradient echo (GRE) is the main workhorse sequence used for susceptibility weighted imaging (SWI), quantitative susceptibility mapping (QSM), and susceptibility tensor imaging (STI). Achieving optimal phase signal-to-noise ratio requires late echo times, thus necessitating a long repetition time (TR). Combined with the large encoding burden of whole-brain coverage with high resolution, this leads to increased scan time. Further, the dipole kernel relating the tissue phase to the underlying susceptibility distribution undersamples the frequency content of the susceptibility map. Scans at multiple head orientations along with calculation of susceptibility through multi-orientation sampling (COSMOS) are one way to effectively mitigate this issue. Additionally, STI requires a minimum of 6 head orientations to solve for the independent tensor elements in each voxel. The requirements of high-resolution imaging with long TR at multiple orientations substantially lengthen the acquisition of COSMOS and STI. The goal of this work is to dramatically speed up susceptibility mapping at multiple head orientations. We demonstrate highly efficient acquisition using 3D-GRE with Wave-CAIPI and dramatically reduce the acquisition time of these protocols. Using R=15-fold acceleration with Wave-CAIPI permits acquisition per head orientation in 90s at 1.1mm isotropic resolution, and 5:35min at 0.5mm isotropic resolution. Since Wave-CAIPI fully harnesses the 3D spatial encoding capability of receive arrays, the maximum g-factor noise amplification remains below 1.30 at 3T and 1.12 at 7T. This allows a 30-min exam for STI with 12 orientations, thus paving the way to its clinical application.

Real diffusion-weighted MRI enabling true signal averaging and increased diffusion contrast

This project aims to characterize the impact of underlying noise distributions on diffusion-weighted imaging. The noise floor is a well-known problem for traditional magnitude-based diffusion-weighted MRI (dMRI) data, leading to biased diffusion model fits and inaccurate signal averaging. Here, we introduce a total-variation-based algorithm to eliminate shot-to-shot phase variations of complex-valued diffusion data with the intention to extract real-valued dMRI datasets. The obtained real-valued diffusion data are no longer superimposed by a noise floor but instead by a zero-mean Gaussian noise distribution, yielding dMRI data without signal bias. We acquired high-resolution dMRI data with strong diffusion weighting and, thus, low signal-to-noise ratio. Both the extracted real-valued and traditional magnitude data were compared regarding signal averaging, diffusion model fitting and accuracy in resolving crossing fibers. Our results clearly indicate that real-valued diffusion data enables idealized conditions for signal averaging. Furthermore, the proposed method enables unbiased use of widely employed linear least squares estimators for model fitting and demonstrates an increased sensitivity to detect secondary fiber directions with reduced angular error. The use of phase-corrected, real-valued data for dMRI will therefore help to clear the way for more detailed and accurate studies of white matter microstructure and structural connectivity on a fine scale.

Signal-to-noise ratio and MR tissue parameters in human brain imaging at 3, 7, and 9.4 tesla using current receive coil arrays

PURPOSE

Relaxation times, transmit homogeneity, signal-to-noise ratio (SNR) and parallel imaging g-factor were determined in the human brain at 3T, 7T, and 9.4T, using standard, tight-fitting coil arrays.

METHODS

The same human subjects were scanned at all three field strengths, using identical sequence parameters and similar 31- or 32-channel receive coil arrays. The SNR of three-dimensional (3D) gradient echo images was determined using a multiple replica approach and corrected with measured flip angle and T2 (*) distributions and the T1 of white matter to obtain the intrinsic SNR. The g-factor maps were derived from 3D gradient echo images with several GRAPPA accelerations.

RESULTS

As expected, T1 values increased, T2 (*) decreased and the B1 -homogeneity deteriorated with increasing field. The SNR showed a distinctly supralinear increase with field strength by a factor of 3.10 ± 0.20 from 3T to 7T, and 1.76 ± 0.13 from 7T to 9.4T over the entire cerebrum. The g-factors did not show the expected decrease, indicating a dominating role of coil design.

CONCLUSION

In standard experimental conditions, SNR increased supralinearly with field strength (SNR ∼ B0 (1.65) ). To take full advantage of this gain, the deteriorating B1 -homogeneity and the decreasing T2 (*) have to be overcome.

Evaluation of highly accelerated simultaneous multi-slice EPI for fMRI

Echo planar imaging (EPI) is the MRI technique that is most widely used for blood oxygen level-dependent (BOLD) functional MRI (fMRI). Recent advances in EPI speed have been made possible with simultaneous multi-slice (SMS) methods which combine acceleration factors M from multiband (MB) radiofrequency pulses and S from simultaneous image refocusing (SIR) to acquire a total of N=S×M images in one echo train, providing up to N times speed-up in total acquisition time over conventional EPI. We evaluated accelerations as high as N=48 using different combinations of S and M which allow for whole brain imaging in as little as 100ms at 3T with a 32 channel head coil. The various combinations of acceleration parameters were evaluated by tSNR as well as BOLD contrast-to-noise ratio (CNR) and information content from checkerboard and movie clips in fMRI experiments. We found that at low acceleration factors (N≤6), setting S=1 and varying M alone yielded the best results in all evaluation metrics, while at acceleration N=8 the results were mixed using both S=1 and S=2 sequences. At higher acceleration factors (N>8), using S=2 yielded maximal BOLD CNR and information content as measured by classification of movie clip frames. Importantly, we found significantly greater BOLD information content using relatively fast TRs in the range of 300ms-600ms compared to a TR of 2s, suggesting that faster TRs capture more information per unit time in task based fMRI.

A 16-channel receive, forced current excitation dual-transmit coil for breast imaging at 7T

Purpose

To enable high spatial and temporal breast imaging resolution via combined use of high field MRI, array coils, and forced current excitation (FCE) multi channel transmit.

Materials and Methods

A unilateral 16-channel receive array insert was designed for use in a transmit volume coil optimized for quadrature operation with dual-transmit RF shimming at 7T. Signal-to-noise ratio (SNR) maps, g-factor maps, and high spatial and temporal resolution in vivo images were acquired to demonstrate the utility of the coil architecture.

Results

The dual-transmit FCE coil provided homogeneous excitation and the array provided an increase in average SNR of 3.3 times (max 10.8, min 1.5) compared to the volume coil in transmit/receive mode. High resolution accelerated in vivo breast imaging demonstrated the ability to achieve isotropic spatial resolution of 0.5 mm within clinically relevant 90 s scan times, as well as the ability to perform 1.0 mm isotropic resolution imaging, 7 s per dynamics, with the use of bidirectional SENSE acceleration of up to R = 9.

Conclusion

The FCE design of the transmit coil easily accommodates the addition of a sixteen channel array coil. The improved spatial and temporal resolution provided by the high-field array coil with FCE dual-channel transmit will ultimately be beneficial in lesion detection and characterization.

Design of parallel transmission radiofrequency pulses robust against respiration in cardiac MRI at 7 Tesla

PURPOSE

Two-spoke parallel transmission (pTX) radiofrequency (RF) pulses have been demonstrated in cardiac MRI at 7T. However, current pulse designs rely on a single set of B1(+)/B0 maps that may not be valid for subsequent scans acquired at another phase of the respiration cycle because of organ displacement. Such mismatches may yield severe excitation profile degradation.

METHODS

B1(+)/B0 maps were obtained, using 16 transmit channels at 7T, at three breath-hold positions: exhale, half-inhale, and inhale. Standard and robust RF pulses were designed using maps obtained at exhale only, and at multiple respiratory positions, respectively. Excitation patterns were analyzed for all positions using Bloch simulations. Flip-angle homogeneity was compared in vivo in cardiac CINE acquisitions.

RESULTS

Standard one- and two-spoke pTX RF pulses are sensitive to breath-hold position, primarily due to B1(+) alterations, with high dependency on excitation trajectory for two spokes. In vivo excitation inhomogeneity varied from nRMSE = 8.2% (exhale) up to 32.5% (inhale) with the standard design; much more stable results were obtained with the robust design with nRMSE = 9.1% (exhale) and 10.6% (inhale).

CONCLUSION

A new pTX RF pulse design robust against respiration induced variations of B1(+)/B0 maps is demonstrated and is expected to have a positive impact on cardiac MRI in breath-hold, free-breathing, and real-time acquisitions.

MS lesions are better detected with 3D T1 gradient-echo than with 2D T1 spin-echo gadolinium-enhanced imaging at 3T

BACKGROUND AND PURPOSE

In multiple sclerosis, gadolinium enhancement is used to classify lesions as active. Regarding the need for a standardized and accurate method for detection of multiple sclerosis activity, we compared 2D-spin-echo with 3D-gradient-echo T1WI for the detection of gadolinium-enhancing MS lesions.

MATERIALS AND METHODS

Fifty-eight patients with MS were prospectively imaged at 3T by using both 2D-spin-echo and 3D-gradient recalled-echo T1WI in random order after the injection of gadolinium. Blinded and independent evaluation was performed by a junior and a senior reader to count gadolinium-enhancing lesions and to characterize their location, size, pattern of enhancement, and the relative contrast between enhancing lesions and the adjacent white matter. Finally, the SNR and relative contrast of gadolinium-enhancing lesions were computed for both sequences by using simulations.

RESULTS

Significantly more gadolinium-enhancing lesions were reported on 3D-gradient recalled-echo than on 2D-spin-echo (n = 59 versus n = 30 for the junior reader, P = .021; n = 77 versus n = 61 for the senior reader, P = .017). The difference between the 2 readers was significant on 2D-spin-echo (P = .044), for which images were less reproducible (κ = 0.51) than for 3D-gradient recalled-echo (κ = 0.65). Further comparisons showed that there were statistically more small lesions (<5 mm) on 3D-gradient recalled-echo than on 2D-spin-echo (P = .04), while other features were similar. Theoretic results from simulations predicted SNR and lesion contrast for 3D-gradient recalled-echo to be better than for 2D-spin-echo for visualization of small enhancing lesions and were, therefore, consistent with clinical observations.

CONCLUSIONS

At 3T, 3D-gradient recalled-echo provides a higher detection rate of gadolinium-enhancing lesions, especially those with smaller size, with a better reproducibility; this finding suggests using 3D-gradient recalled-echo to detect MS activity, with potential impact in initiation, monitoring, and optimization of therapy.

Metamaterial magnetoinductive lens performance as a function of field strength

Metamaterials are artificial composites that exhibit exotic electromagnetic properties, as the ability of metamaterial slabs to behave like lenses with sub-wavelength resolution for the electric or the magnetic field. In previous works, the authors investigated magnetic resonance imaging (MRI) applications of metamaterial slabs that behave like lenses for the radiofrequency magnetic field. In particular, the authors investigated the ability of MRI metamaterial lenses to increase the signal-to-noise ratio (SNR) of surface coils, and to localize the field of view (FOV) of the coils, which is of interest for parallel MRI (pMRI) applications. A metamaterial lens placed between a surface coil and the tissue enhances the sensitivity of the coil. Although the metamaterial lens introduces losses which add to the losses of the tissue, the enhancement of the sensitivity can compensate these additional losses and the SNR of the coil is increased. In a previous work, an optimization procedure was followed to find a metamaterial structure with minimum losses that will maximize the SNR. This structure was termed magnetoinductive (MI) lens by the authors. The properties of surface coils in the presence of MI lenses were investigated in previous works at the proton frequency of 1.5 T systems. The different frequency dependence of the losses in both the MI lenses and the tissue encouraged us to investigate the performance of MI lenses at different frequencies. Thus, in the present work, the SNR and the pMRI ability of MI lenses are investigated as a function of field strength. A numerical analysis is carried out with an algorithm developed by the authors to predict the SNR behavior of a surface coil loaded with a MI lens at the proton frequencies of 0.5 T, 1.5 T and 3 T systems. The results show that, at 0.5 T, there is a gain in the SNR for short distances, but the SNR is highly degraded at deeper distances. However, at 1.5 T and 3T, the MI lenses provide a gain in the SNR up to a certain penetration depth, which is deeper at 3T, and do not degrade the SNR at deeper distances. These numerical results are checked by means of an experiment. Moreover, a second experiment developed with two-channel arrays of surface coils loaded with MI lenses shows that the pMRI ability of the lenses also improves from 1.5 T to 3 T. This improvement was quantified by means of the calculation of the GRAPPA g-factor.

Magnetic resonance imaging at ultrahigh fields

Since the introduction of 4 T human systems in three academic laboratories circa 1990, rapid progress in imaging and spectroscopy studies in humans at 4 T and animal model systems at 9.4 T have led to the introduction of 7 T and higher magnetic fields for human investigation at about the turn of the century. Work conducted on these platforms has demonstrated the existence of significant advantages in SNR and biological information content at these ultrahigh fields, as well as the presence of numerous challenges. Primary difference from lower fields is the deviation from the near field regime; at the frequencies corresponding to hydrogen resonance conditions at ultrahigh fields, the RF is characterized by attenuated traveling waves in the human body, which leads to image nonuniformities for a given sample-coil configuration because of interferences. These nonuniformities were considered detrimental to the progress of imaging at high field strengths. However, they are advantageous for parallel imaging for signal reception and parallel transmission, two critical technologies that account, to a large extend, for the success of ultrahigh fields. With these technologies, and improvements in instrumentation and imaging methods, ultrahigh fields have provided unprecedented gains in imaging of brain function and anatomy, and started to make inroads into investigation of the human torso and extremities. As extensive as they are, these gains still constitute a prelude to what is to come given the increasingly larger effort committed to ultrahigh field research and development of ever better instrumentation and techniques.

Coil combination of multichannel MRSI data at 7 T: MUSICAL

The goal of this study was to evaluate a new method of combining multi-channel (1)H MRSI data by direct use of a matching imaging scan as a reference, rather than computing sensitivity maps. Seven healthy volunteers were measured on a 7-T MR scanner using a head coil with a 32-channel array coil for receive-only and a volume coil for receive/transmit. The accuracy of prediction of the phase of the (1)H MRSI data with a fast imaging pre-scan was investigated with the volume coil. The array coil (1)H MRSI data were combined using matching imaging data as coil combination weights. The signal-to-noise ratio (SNR), spectral quality, metabolic map quality and Cramér-Rao lower bounds were then compared with the data obtained by two standard methods, i.e. using sensitivity maps and the first free induction decay (FID) data point. Additional noise decorrelation was performed to further optimize the SNR gain. The new combination method improved significantly the SNR (+29%), overall spectral quality and visual appearance of metabolic maps, and lowered the Cramér-Rao lower bounds (-34%), compared with the combination method based on the first FID data point. The results were similar to those obtained by the combination method using sensitivity maps, but the new method increased the SNR slightly (+1.7%), decreased the algorithm complexity, required no reference coil and pre-phased all spectra correctly prior to spectral processing. Noise decorrelation further increased the SNR by 13%. The proposed method is a fast, robust and simple way to improve the coil combination in (1)H MRSI of the human brain at 7 T, and could be extended to other (1)H MRSI techniques.

Pushing spatial and temporal resolution for functional and diffusion MRI in the Human Connectome Project

The Human Connectome Project (HCP) relies primarily on three complementary magnetic resonance (MR) methods. These are: 1) resting state functional MR imaging (rfMRI) which uses correlations in the temporal fluctuations in an fMRI time series to deduce 'functional connectivity'; 2) diffusion imaging (dMRI), which provides the input for tractography algorithms used for the reconstruction of the complex axonal fiber architecture; and 3) task based fMRI (tfMRI), which is employed to identify functional parcellation in the human brain in order to assist analyses of data obtained with the first two methods. We describe technical improvements and optimization of these methods as well as instrumental choices that impact speed of acquisition of fMRI and dMRI images at 3T, leading to whole brain coverage with 2 mm isotropic resolution in 0.7 s for fMRI, and 1.25 mm isotropic resolution dMRI data for tractography analysis with three-fold reduction in total dMRI data acquisition time. Ongoing technical developments and optimization for acquisition of similar data at 7 T magnetic field are also presented, targeting higher spatial resolution, enhanced specificity of functional imaging signals, mitigation of the inhomogeneous radio frequency (RF) fields, and reduced power deposition. Results demonstrate that overall, these approaches represent a significant advance in MR imaging of the human brain to investigate brain function and structure.

Visualization of intra-thalamic nuclei with optimized white-matter-nulled MPRAGE at 7T

Novel MR image acquisition strategies have been investigated to elicit contrast within the thalamus, but direct visualization of individual thalamic nuclei remains a challenge because of their small size and the low intrinsic contrast between adjacent nuclei. We present a step-by-step specific optimization of the 3D MPRAGE pulse sequence at 7T to visualize the intra-thalamic nuclei. We first measured T1 values within different sub-regions of the thalamus at 7T in 5 individuals. We used these to perform simulations and sequential experimental measurements (n=17) to tune the parameters of the MPRAGE sequence. The optimal set of parameters was used to collect high-quality data in 6 additional volunteers. Delineation of thalamic nuclei was performed twice by one rater and MR-defined nuclei were compared to the classic Morel histological atlas. T1 values within the thalamus ranged from 1400ms to 1800ms for adjacent nuclei. Using these values for theoretical evaluations combined with in vivo measurements, we showed that a short inversion time (TI) close to the white matter null regime (TI=670ms) enhanced the contrast between the thalamus and the surrounding tissues, and best revealed intra-thalamic contrast. At this particular nulling regime, lengthening the time between successive inversion pulses (TS=6000ms) increased the thalamic signal and contrast and lengthening the α pulse train time (N*TR) further increased the thalamic signal. Finally, a low flip angle during the gradient echo acquisition (α=4°) was observed to mitigate the blur induced by the evolution of the magnetization along the α pulse train. This optimized set of parameters enabled the 3D delineation of 15 substructures in all 6 individuals; these substructures corresponded well with the known anatomical structures of the thalamus based on the classic Morel atlas. The mean Euclidean distance between the centers of mass of MR- and Morel atlas-defined nuclei was 2.67mm (±1.02mm). The reproducibility of the MR-defined nuclei was excellent with intraclass correlation coefficient measured at 0.997 and a mean Euclidean distance between corresponding centers of mass found at first versus second readings of 0.69mm (±0.38mm). This 7T strategy paves the way to better identification of thalamic nuclei for neurosurgical planning and investigation of regional changes in neurological disorders.

Transmit and receive RF fields determination from a single low-tip-angle gradient-echo scan by scaling of SVD data

PURPOSE

A new method, called Transmit and Receive Patterns from Low-Tip-angle gradient-Echo Images (TRIPLET), is described which simultaneously maps the B1+ and B1- fields of a transmit/receive radiofrequency coil array. The input data are low-tip-angle gradient-echo images, which can be acquired in a relatively short scanning time.

THEORY AND METHODS

For each voxel in the field of view, a matrix can be assembled with the low-tip-angle gradient-echo image values of the radiofrequency coil array. Applying the singular value decomposition to those matrices, datasets are obtained which show a high resemblance with the true B1+ and B1- fields. These datasets are a voxel-wise scaled version of the true radiofrequency maps. The channel independent scaling parameters can be found by implicitly forcing the reconstructed fields to be solutions of the Maxwell equations. This is achieved by introducing a multipole expansion consisting of Bessel/Fourier functions.

RESULTS

Two FDTD simulated radiofrequency fields for two coil array combinations at 7 T and a measured, in vivo dataset at 7 T are investigated to illustrate the singular value decomposition analysis of the low-tip-angle gradient-echo images and to show how the B1+ and B1- fields can be reconstructed by Transmit and Receive Patterns from Low-Tip-angle gradient-Echo Images.

CONCLUSION

The Transmit and Receive Patterns from Low-Tip-angle gradient-Echo Images algorithm can convert the datasets from singular value decomposition analysis of low-tip-angle gradient-echo images to true B1+ and B1- fields.

Nineteen-channel receive array and four-channel transmit array coil for cervical spinal cord imaging at 7T

PURPOSE

To design and validate a radiofrequency (RF) array coil for cervical spinal cord imaging at 7T.

METHODS

A 19-channel receive array with a four-channel transmit array was developed on a close-fitting coil former at 7T. Transmit efficiency and specific absorption rate were evaluated in a B1 (+) mapping study and an electromagnetic model. Receive signal-to-noise ratio (SNR) and noise amplification for parallel imaging were evaluated and compared with a commercial 3T 19-channel head-neck array and a 7T four-channel spine array. The performance of the array was qualitatively demonstrated in human volunteers using high-resolution imaging (down to 300 μm in-plane).

RESULTS

The transmit and receive arrays showed good bench performance. The SNR was approximately 4.2-fold higher in the 7T receive array at the location of the cord with respect to the 3T coil. The g-factor results showed an additional acceleration was possible with the 7T array. In vivo imaging was feasible and showed high SNR and tissue contrast.

CONCLUSION

The highly parallel transmit and receive arrays were demonstrated to be fit for spinal cord imaging at 7T. The high sensitivity of the receive coil combined with ultra-high field will likely improve investigations of microstructure and tissue segmentation in the healthy and pathological spinal cord.

Evaluation of slice accelerations using multiband echo planar imaging at 3 T

We evaluate residual aliasing among simultaneously excited and acquired slices in slice accelerated multiband (MB) echo planar imaging (EPI). No in-plane accelerations were used in order to maximize and evaluate achievable slice acceleration factors at 3 T. We propose a novel leakage (L-) factor to quantify the effects of signal leakage between simultaneously acquired slices. With a standard 32-channel receiver coil at 3 T, we demonstrate that slice acceleration factors of up to eight (MB=8) with blipped controlled aliasing in parallel imaging (CAIPI), in the absence of in-plane accelerations, can be used routinely with acceptable image quality and integrity for whole brain imaging. Spectral analyses of single-shot fMRI time series demonstrate that temporal fluctuations due to both neuronal and physiological sources were distinguishable and comparable up to slice-acceleration factors of nine (MB=9). The increased temporal efficiency could be employed to achieve, within a given acquisition period, higher spatial resolution, increased fMRI statistical power, multiple TEs, faster sampling of temporal events in a resting state fMRI time series, increased sampling of q-space in diffusion imaging, or more quiet time during a scan.

Ultra high spatial and temporal resolution breast imaging at 7T

There is a need to obtain higher specificity in the detection of breast lesions using MRI. To address this need, Dynamic Contrast-Enhanced (DCE) MRI has been combined with other structural and functional MRI techniques. Unfortunately, owing to time constraints structural images at ultra-high spatial resolution can generally not be obtained during contrast uptake, whereas the relatively low spatial resolution of functional imaging (e.g. diffusion and perfusion) limits the detection of small lesions. To be able to increase spatial as well as temporal resolution simultaneously, the sensitivity of MR detection needs to increase as well as the ability to effectively accelerate the acquisition. The required gain in signal-to-noise ratio (SNR) can be obtained at 7T, whereas acceleration can be obtained with high-density receiver coil arrays. In this case, morphological imaging can be merged with DCE-MRI, and other functional techniques can be obtained at higher spatial resolution, and with less distortion [e.g. Diffusion Weighted Imaging (DWI)]. To test the feasibility of this concept, we developed a unilateral breast coil for 7T. It comprises a volume optimized dual-channel transmit coil combined with a 30-channel receive array coil. The high density of small coil elements enabled efficient acceleration in any direction to acquire ultra high spatial resolution MRI of close to 0.6 mm isotropic detail within a temporal resolution of 69 s, high spatial resolution MRI of 1.5 mm isotropic within an ultra high temporal resolution of 6.7 s and low distortion DWI at 7T, all validated in phantoms, healthy volunteers and a patient with a lesion in the right breast classified as Breast Imaging Reporting and Data System (BI-RADS) IV.

Noise amplification in parallel whole-head ultra-low-field magnetic resonance imaging using 306 detectors

In ultra-low-field magnetic resonance imaging, arrays of up to hundreds of highly sensitive superconducting quantum interference devices (SQUIDs) can be used to detect the weak magnetic fields emitted by the precessing magnetization. Here, we investigate the noise amplification in sensitivity-encoded ultra-low-field MRI at various acceleration rates using a SQUID array consisting of 102 magnetometers, 102 gradiometers, or 306 magnetometers and gradiometers, to cover the whole head. Our results suggest that SQUID arrays consisting of 102 magnetometers and 102 gradiometers are similar in g-factor distribution. A SQUID array of 306 sensors (102 magnetometers and 204 gradiometers) only marginally improves the g-factor. Corroborating with previous studies, the g-factor in 2D sensitivity-encoded ultra-low-field MRI with 9 to 16-fold 2D accelerations using the SQUID array studied here may be acceptable.

Readout-segmented echo-planar imaging improves the diagnostic performance of diffusion-weighted MR breast examinations at 3.0 T

PURPOSE

To qualitatively and quantitatively compare the diagnostic value of diffusion-weighted (DW) magnetic resonance (MR) imaging based on standard single-shot echo-planar imaging and readout-segmented echo-planar imaging in patients with breast cancer at 3.0 T.

MATERIALS AND METHODS

Institutional review board approval and written informed consent were obtained. Forty-seven patients with 49 histopathologically verified lesions were included in this study. In all patients, DW imaging, with single-shot echo-planar imaging and readout-segmented echo-planar imaging with comparable imaging parameters, was performed with a 3.0-T MR imager. Two independent readers visually assessed image quality and lesion conspicuity, and image properties (ie, signal-to-noise ratio, contrast, geometric distortions) were quantified. Regions of interest were drawn in all lesions (28 malignant, 21 benign) and in the normal breast parenchyma to investigate differences in apparent diffusion coefficient (ADC). Diagnostic accuracy was calculated on the basis of an ADC threshold of 1.25 × 10(-3) mm(2)/sec.

RESULTS

Each reader found a higher diagnostic accuracy for readout-segmented (96%) than for single-shot (90%) echo-planar imaging. The area under the curve for readout-segmented echo-planar imaging (0.981) was significantly larger than for single-shot echo-planar imaging (0.867) (P = .026). There was no significant difference in the ADC obtained by using either DW imaging method. Lesion conspicuity and image quality of readout-segmented echo-planar imaging were rated superior to those of single-shot echo-planar imaging (P < .001). Readout-segmented echo-planar imaging reduced geometric distortions by a factor of three.

CONCLUSION

DW imaging based on readout-segmented echo-planar imaging provided significantly higher image quality and lesion conspicuity than single-shot echo-planar imaging by reducing geometric distortions, image blurring, and artifact level with a clinical high-field-strength MR imager. Thereby, readout-segmented echo-planar imaging reached a higher diagnostic accuracy for the differentiation of benign and malignant breast lesions.

Implementation of vascular-space-occupancy MRI at 7T

Vascular-space-occupancy (VASO) MRI exploits the difference between blood and tissue T1 to null blood signal and measure cerebral blood volume changes using the residual tissue signal. VASO imaging is more difficult at higher field because of sensitivity loss due to the convergence of tissue and blood T1 values and increased contamination from blood-oxygenation-level-dependent (BOLD) effects. In addition, compared to 3T, 7T MRI suffers from increased geometrical distortions, e.g., when using echo-planar-imaging, and from increased power deposition, the latter especially problematic for the spin-echo-train sequences commonly used for VASO MRI. Third, non-steady-state blood spin effects become substantial at 7T when only a head coil is available for radiofrequency transmit. In this study, the magnetization-transfer-enhanced-VASO approach was applied to maximize tissue-blood signal difference, which boosted signal-to-noise ratio by 149% ± 13% (n = 7) compared to VASO. Second, a 3D fast gradient-echo sequence with low flip-angle (7°) and short echo-time (1.8 ms) was used to minimize the BOLD effect and to reduce image distortion and power deposition. Finally, a magnetization-reset technique was combined with a motion-sensitized-driven-equilibrium approach to suppress three types of non-steady-state spins. Our initial functional MRI results in normal human brains at 7T with this optimized VASO sequence showed better signal-to-noise ratio than at 3T.

The road to functional imaging and ultrahigh fields

The Center for Magnetic Resonance (CMRR) at the University of Minnesota was one of the laboratories where the work that simultaneously and independently introduced functional magnetic resonance imaging (fMRI) of human brain activity was carried out. However, unlike other laboratories pursuing fMRI at the time, our work was performed at 4T magnetic field and coincided with the effort to push human magnetic resonance imaging to field strength significantly beyond 1.5T which was the high-end standard of the time. The human fMRI experiments performed in CMRR were planned between two colleagues who had known each other and had worked together previously in Bell Laboratories, namely Seiji Ogawa and myself, immediately after the Blood Oxygenation Level Dependent (BOLD) contrast was developed by Seiji. We were waiting for our first human system, a 4T system, to arrive in order to attempt at imaging brain activity in the human brain and these were the first experiments we performed on the 4T instrument in CMRR when it became marginally operational. This was a prelude to a subsequent systematic push we initiated for exploiting higher magnetic fields to improve the accuracy and sensitivity of fMRI maps, first going to 9.4T for animal model studies and subsequently developing a 7T human system for the first time. Steady improvements in high field instrumentation and ever expanding armamentarium of image acquisition and engineering solutions to challenges posed by ultrahigh fields have brought fMRI to submillimeter resolution in the whole brain at 7T, the scale necessary to reach cortical columns and laminar differentiation in the whole brain. The solutions that emerged in response to technological challenges posed by 7T also propagated and continues to propagate to lower field clinical systems, a major advantage of the ultrahigh fields effort that is underappreciated. Further improvements at 7T are inevitable. Further translation of these improvements to lower field clinical systems to achieve new capabilities and to magnetic fields significantly higher than 7T to enable human imaging is inescapable.

k-space and q-space: combining ultra-high spatial and angular resolution in diffusion imaging using ZOOPPA at 7 T

There is ongoing debate whether using a higher spatial resolution (sampling k-space) or a higher angular resolution (sampling q-space angles) is the better way to improve diffusion MRI (dMRI) based tractography results in living humans. In both cases, the limiting factor is the signal-to-noise ratio (SNR), due to the restricted acquisition time. One possible way to increase the spatial resolution without sacrificing either SNR or angular resolution is to move to a higher magnetic field strength. Nevertheless, dMRI has not been the preferred application for ultra-high field strength (7 T). This is because single-shot echo-planar imaging (EPI) has been the method of choice for human in vivo dMRI. EPI faces several challenges related to the use of a high resolution at high field strength, for example, distortions and image blurring. These problems can easily compromise the expected SNR gain with field strength. In the current study, we introduce an adapted EPI sequence in conjunction with a combination of ZOOmed imaging and Partially Parallel Acquisition (ZOOPPA). We demonstrate that the method can produce high quality diffusion-weighted images with high spatial and angular resolution at 7 T. We provide examples of in vivo human dMRI with isotropic resolutions of 1 mm and 800 μm. These data sets are particularly suitable for resolving complex and subtle fiber architectures, including fiber crossings in the white matter, anisotropy in the cortex and fibers entering the cortex.

Comparison between eight- and sixteen-channel TEM transceive arrays for body imaging at 7 T

Eight- and sixteen-channel transceive stripline/TEM body arrays were compared at 7 T (297 MHz) both in simulation and experiment. Despite previous demonstrations of similar arrays for use in body applications, a quantitative comparison of the two configurations has not been undertaken to date. Results were obtained on a male pelvis for assessing transmit, signal to noise ratio, and parallel imaging performance and to evaluate local power deposition versus transmit B(1) (B(1) (+) ). All measurements and simulations were conducted after performing local B(1) (+) phase shimming in the region of the prostate. Despite the additional challenges of decoupling immediately adjacent coils, the sixteen-channel array demonstrated improved or nearly equivalent performance to the eight-channel array based on the evaluation criteria. Experimentally, transmit performance and signal to noise ratio were 22% higher for the sixteen-channel array while significantly increased reduction factors were achievable in the left-right direction for parallel imaging. Finite difference time domain simulations demonstrated similar results with respect to transmit and parallel imaging performance, however, a higher transmit efficiency advantage of 33% was predicted. Simulations at both 3 and 7 T verified the expected parallel imaging improvements with increasing field strength and showed that, for a specific B(1) (+) shimming strategy used, the sixteen-channel array exhibited lower local and global specific absorption rate for a given B(1) (+) .

High-field FMRI for human applications: an overview of spatial resolution and signal specificity

In the last decade, dozens of 7 Tesla scanners have been purchased or installed around the world, while 3 Tesla systems have become a standard. This increased interest in higher field strengths is driven by a demonstrated advantage of high fields for available signal-to-noise ratio (SNR) in the magnetic resonance signal. Functional imaging studies have additional advantages of increases in both the contrast and the spatial specificity of the susceptibility based BOLD signal. One use of this resultant increase in the contrast to noise ratio (CNR) for functional MRI studies at high field is increased image resolution. However, there are many factors to consider in predicting exactly what kind of resolution gains might be made at high fields, and what the opportunity costs might be. The first part of this article discusses both hardware and image quality considerations for higher resolution functional imaging. The second part draws distinctions between image resolution, spatial specificity, and functional specificity of the fMRI signals that can be acquired at high fields, suggesting practical limitations for attainable resolutions of fMRI experiments at a given field, given the current state of the art in imaging techniques. Finally, practical resolution limitations and pulse sequence options for studies in human subjects are considered.

Recent progress in high-resolution functional MRI

PURPOSE OF REVIEW

For functional MRI (fMRI), as for any imaging technique, the higher the spatial resolution, the more the details it can reveal. This review will discuss the factors restricting the spatial resolution of fMRI, describe high-resolution fMRI (HR-fMRI) applications in neuroscience and outline a few research areas for future HR-fMRI studies.

RECENT FINDINGS

HR-fMRI has been successfully used to map fine cortical architectures and reveal cortical laminar structures and subcortical structures. HR-fMRI has also played important roles in resolving controversies regarding modular representations in the ventral visual pathway and interpretations of multivariate pattern analysis results.

SUMMARY

Real-time HR-fMRI as well as high-resolution anatomical MRI may emerge as indispensable tools for surgical planning, diagnosis of neurological diseases and targeting of deep brain stimulation.

Diffusion MR imaging: basic principles

From their origin as simple techniques primarily used for detecting acute cerebral ischemia, diffusion MR imaging techniques have rapidly evolved into a versatile set of tools that provide the only noninvasive means of characterizing brain microstructure and connectivity, becoming a mainstay of both clinical and investigational brain MR imaging. In this article, the basic principles required for understanding diffusion MR imaging techniques are reviewed with clinical neuroradiologists in mind.

Parallel traveling-wave MRI: a feasibility study

Traveling-wave magnetic resonance imaging utilizes far fields of a single-piece patch antenna in the magnet bore to generate radio frequency fields for imaging large-size samples, such as the human body. In this work, the feasibility of applying the "traveling-wave" technique to parallel imaging is studied using microstrip patch antenna arrays with both the numerical analysis and experimental tests. A specific patch array model is built and each array element is a microstrip patch antenna. Bench tests show that decoupling between two adjacent elements is better than -26-dB while matching of each element reaches -36-dB, demonstrating excellent isolation performance and impedance match capability. The sensitivity patterns are simulated and g-factors are calculated for both unloaded and loaded cases. The results on B 1- sensitivity patterns and g-factors demonstrate the feasibility of the traveling-wave parallel imaging. Simulations also suggest that different array configuration such as patch shape, position and orientation leads to different sensitivity patterns and g-factor maps, which provides a way to manipulate B(1) fields and improve the parallel imaging performance. The proposed method is also validated by using 7T MR imaging experiments.

Clinical application of 7.0 T magnetic resonance images in Gamma Knife radiosurgery for a patient with brain metastases

In the study we assessed the distortion of 7.0 T magnetic resonance (MR) images in reference to 1.5 T MR images in the radiosurgery of metastatic brain tumors. Radiosurgery with Gamma Knife Perfexion was performed for the treatment of a 54-yr-old female patient with multiple brain metastases by the co-registered images of the 7.0 T and 1.5 T magnetic resonance images (MRI). There was no significant discrepancy in the positions of anterior and posterior commissures as well as the locations of four metastatic brain tumors in the co-registered images between 7.0 T and 1.5 T MRI with better visualization of the anatomical details in 7.0 T MR images. This study demonstrates for the first time that 7.0 T MR images can be safely utilized in Perfexion Gamma Knife radiosurgery for the treatment of metastatic brain tumors. Furthermore 7.0 T MR images provide better visualization of brain tumors without image distortion in comparison to 1.5 T MR images.

Calculation of radiofrequency electromagnetic fields and their effects in MRI of human subjects

Radiofrequency magnetic fields are critical to nuclear excitation and signal reception in magnetic resonance imaging. The interactions between these fields and human tissues in anatomical geometries results in a variety of effects regarding image integrity and safety of the human subject. In recent decades, numerical methods of calculation have been used increasingly to understand the effects of these interactions and aid in engineering better, faster, and safer equipment and methods. As magnetic resonance imaging techniques and technology have evolved through the years, so to have the requirements for meaningful interpretation of calculation results. Here, we review the basic physics of radiofrequency electromagnetics in magnetic resonance imaging and discuss a variety of ways radiofrequency field calculations are used in magnetic resonance imaging in engineering and safety assurance from simple systems and sequences through advanced methods of development for the future.

Traveling-wave RF shimming and parallel MRI

At sufficiently high Larmor frequencies, traveling electromagnetic waves along a magnet bore can be used for remote magnetic resonance excitation and detection, effectively using the bore as a waveguide. So far, this approach has relied only on the lowest waveguide modes and thus has not supported multiple-channel operation for radiofrequency shimming and parallel imaging. In this work, this limitation is addressed by establishing a larger number of propagating modes and tapping their spatial field diversity with multiple waveguide ports. The number of available modes is increased by loading with dielectric inserts; the ports are implemented by stub and loop couplers at the end of a waveguide extension. The resulting traveling-wave array, operated at 298 MHz in a 7T whole-body magnet, is shown to enable radiofrequency shimming as well as parallel imaging with commonly used acceleration factors. The last part of the study concerns the amount of dielectric loading that is required. For the given Larmor frequency and bore dimensions, it is found that rather few water-filled inserts, occupying ∼5% of the bore cross-section, are sufficient for effective parallel imaging.

Multiplexed echo planar imaging for sub-second whole brain FMRI and fast diffusion imaging

Echo planar imaging (EPI) is an MRI technique of particular value to neuroscience, with its use for virtually all functional MRI (fMRI) and diffusion imaging of fiber connections in the human brain. EPI generates a single 2D image in a fraction of a second; however, it requires 2–3 seconds to acquire multi-slice whole brain coverage for fMRI and even longer for diffusion imaging. Here we report on a large reduction in EPI whole brain scan time at 3 and 7 Tesla, without significantly sacrificing spatial resolution, and while gaining functional sensitivity. The multiplexed-EPI (M-EPI) pulse sequence combines two forms of multiplexing: temporal multiplexing (m) utilizing simultaneous echo refocused (SIR) EPI and spatial multiplexing (n) with multibanded RF pulses (MB) to achieve m×n images in an EPI echo train instead of the normal single image. This resulted in an unprecedented reduction in EPI scan time for whole brain fMRI performed at 3 Tesla, permitting TRs of 400 ms and 800 ms compared to a more conventional 2.5 sec TR, and 2–4 times reductions in scan time for HARDI imaging of neuronal fibertracks. The simultaneous SE refocusing of SIR imaging at 7 Tesla advantageously reduced SAR by using fewer RF refocusing pulses and by shifting fat signal out of the image plane so that fat suppression pulses were not required. In preliminary studies of resting state functional networks identified through independent component analysis, the 6-fold higher sampling rate increased the peak functional sensitivity by 60%. The novel M-EPI pulse sequence resulted in a significantly increased temporal resolution for whole brain fMRI, and as such, this new methodology can be used for studying non-stationarity in networks and generally for expanding and enriching the functional information.

High-resolution magnetization-prepared 3D-FLAIR imaging at 7.0 Tesla

The aim of the present study is to develop a submillimeter volumetric (three-dimensional) fluid-attenuated inversion recovery sequence at 7T. Implementation of the fluid-attenuated inversion recovery sequence is difficult as increased T(1) weighting from prolonged T(1) constants at 7T dominate the desired T(2) contrast and yield suboptimal signal-to-noise ratio. Magnetization preparation was used to reduce T(1) weighting and improve the T(2) weighting. Also, practical challenges limit the implementation. Long refocusing trains with low flip angles were used to mitigate the specific absorption rate constraints. This resulted in a three-dimensional magnetization preparation fluid-attenuated inversion recovery sequence with 0.8 x 0.8 x 0.8 = 0.5 mm(3) resolution in a clinically acceptable scan time. The contrast-to-noise ratio between gray matter and white matter (contrast-to-noise ratio = signal-to-noise ratio [gray matter] - signal-to-noise ratio [white matter]) increased from 12 +/- 9 without magnetization preparation to 28 +/- 8 with magnetization preparation (n = 12). The signal-to-noise ratio increased for white matter by 13 +/- 6% and for gray matter by 48 +/- 15%. In conclusion, three-dimensional fluid-attenuated inversion recovery with high resolution and full brain coverage is feasible at 7T. Magnetization preparation reduces the T(1) weighting, thereby improving the T(2) weighted contrast and signal-to-noise ratio.

A 32-channel lattice transmission line array for parallel transmit and receive MRI at 7 tesla

Transmit and receive RF coil arrays have proven to be particularly beneficial for ultra-high-field MR. Transmit coil arrays enable such techniques as B(1) (+) shimming to substantially improve transmit B(1) homogeneity compared to conventional volume coil designs, and receive coil arrays offer enhanced parallel imaging performance and SNR. Concentric coil arrangements hold promise for developing transceiver arrays incorporating large numbers of coil elements. At magnetic field strengths of 7 tesla and higher where the Larmor frequencies of interest can exceed 300 MHz, the coil array design must also overcome the problem of the coil conductor length approaching the RF wavelength. In this study, a novel concentric arrangement of resonance elements built from capacitively-shortened half-wavelength transmission lines is presented. This approach was utilized to construct an array with whole-brain coverage using 16 transceiver elements and 16 receive-only elements, resulting in a coil with a total of 16 transmit and 32 receive channels.

Performance of external and internal coil configurations for prostate investigations at 7 T

Three different coil configurations were evaluated through simulation and experimentally to determine safe operating limits and evaluate subject size-dependent performance for prostate imaging at 7 T. The coils included a transceiver endorectal coil (trERC), a 16-channel transceiver external surface array (trESA) and a trESA combined with a receive-only ERC (trESA+roERC). Although the transmit B(1) (B 1+) homogeneity was far superior for the trESA, the maximum achievable B 1+ is subject size dependent and limited by transmit chain losses and amplifier performance. For the trERC, limitations in transmit homogeneity greatly compromised image quality and limited coverage of the prostate. Despite these challenges, the high peak B 1+ close to the trERC and subject size-independent performance provides potential advantages especially for spectroscopic localization where high-bandwidth radiofrequency pulses are required. On the receive side, the combined trESA+roERC provided the highest signal-to-noise ratio and improved homogeneity over the trERC resulting in better visualization of the prostate and surrounding anatomy. In addition, the parallel imaging performance of the trESA+roERC holds strong promise for diffusion-weighted imaging and dynamic contrast-enhanced MRI.

Multiband multislice GE-EPI at 7 tesla, with 16-fold acceleration using partial parallel imaging with application to high spatial and temporal whole-brain fMRI

Parallel imaging in the form of multiband radiofrequency excitation, together with reduced k-space coverage in the phase-encode direction, was applied to human gradient echo functional MRI at 7 T for increased volumetric coverage and concurrent high spatial and temporal resolution. Echo planar imaging with simultaneous acquisition of four coronal slices separated by 44mm and simultaneous 4-fold phase-encoding undersampling, resulting in 16-fold acceleration and up to 16-fold maximal aliasing, was investigated. Task/stimulus-induced signal changes and temporal signal behavior under basal conditions were comparable for multiband and standard single-band excitation and longer pulse repetition times. Robust, whole-brain functional mapping at 7 T, with 2 x 2 x 2mm(3) (pulse repetition time 1.25 sec) and 1 x 1 x 2mm(3) (pulse repetition time 1.5 sec) resolutions, covering fields of view of 256 x 256 x 176 mm(3) and 192 x 172 x 176 mm(3), respectively, was demonstrated with current gradient performance.

Shielded microstrip array for 7T human MR imaging

The high-frequency transceiver array based on the microstrip transmission line design is a promising technique for ultrahigh field magnetic resonance imaging (MRI) signal excitation and reception. However, with the increase of radio-frequency (RF) channels, the size of the ground plane in each microstrip coil element is usually not sufficient to provide a perfect ground. Consequently, the transceiver array may suffer from cable resonance, lower Q-factors, and imaging quality degradations. In this paper, we present an approach to improving the performance of microstrip transceiver arrays by introducing RF shielding outside the microstrip array and the feeding coaxial cables. This improvement reduced interactions among cables, increased resonance stability, and Q-factors, and thus improved imaging quality. An experimental method was also introduced and utilized for quantitative measurement and evaluation of RF coil resonance stability or "cable resonance" behavior.

GPU accelerated FDTD solver and its application in MRI

The finite difference time domain (FDTD) method is a popular technique for computational electromagnetics (CEM). The large computational power often required, however, has been a limiting factor for its applications. In this paper, we will present a graphics processing unit (GPU)-based parallel FDTD solver and its successful application to the investigation of a novel B1 shimming scheme for high-field magnetic resonance imaging (MRI). The optimized shimming scheme exhibits considerably improved transmit B(1) profiles. The GPU implementation dramatically shortened the runtime of FDTD simulation of electromagnetic field compared with its CPU counterpart. The acceleration in runtime has made such investigation possible, and will pave the way for other studies of large-scale computational electromagnetic problems in modern MRI which were previously impractical.