Motion-Correction Enabled Ultra-High Resolution In-Vivo 7T-MRI of the Brain

Leveraging ultra-high field (7T) MRI in psychiatric research

Non-invasive brain imaging has played a critical role in establishing our understanding of the neural properties that contribute to the emergence of psychiatric disorders. However, characterizing core neurobiological mechanisms of psychiatric symptomatology requires greater structural, functional, and neurochemical specificity than is typically obtainable with standard field strength MRI acquisitions (e.g., 3T). Ultra-high field (UHF) imaging at 7 Tesla (7T) provides the opportunity to identify neurobiological systems that confer risk, determine etiology, and characterize disease progression and treatment outcomes of major mental illnesses. Increases in scanner availability, regulatory approval, and sequence availability have made the application of UHF to clinical cohorts more feasible than ever before, yet the application of UHF approaches to the study of mental health remains nascent. In this technical review, we describe core neuroimaging methodologies which benefit from UHF acquisition, including high resolution structural and functional imaging, single (1H) and multi-nuclear (e.g., 31P) MR spectroscopy, and quantitative MR techniques for assessing brain tissue iron and myelin. We discuss advantages provided by 7T MRI, including higher signal- and contrast-to-noise ratio, enhanced spatial resolution, increased test-retest reliability, and molecular and neurochemical specificity, and how these have begun to uncover mechanisms of psychiatric disorders. Finally, we consider current limitations of UHF in its application to clinical cohorts, and point to ongoing work that aims to overcome technical hurdles through the continued development of UHF hardware, software, and protocols.

B0-insensitive image navigators for prospective motion-corrected MRS with localized second-order shimming

PURPOSE

Localized shimming in single-voxel MRS often results in large B0 inhomogeneity outside the volume-of-interest. This causes unacceptable degradation in motion navigator images. Switching back and forth between whole-brain shim and localized shim is possible for linear shims, but not for higher-order shims. Here we propose motion navigators largely insensitive to B0 inhomogeneity for prospective motion-corrected MRS with localized higher-order shimming.

METHODS

A recent fast high-resolution motion navigator based on spiral-in/out k-space trajectories and multislice-to-volume registration was modified by splitting the readout into multiple shot interleaves which shortened the echo time and reduced the effect of B0 inhomogeneity. The performance of motion correction was assessed in healthy subjects in the prefrontal cortex using a sLASER sequence at 3T (N = 5) and 7T (N = 5).

RESULTS

With multiple spatial interleaves, excellent quality navigator images were acquired in the whole brain in spite of large B0 inhomogeneity outside the MRS voxel. The total duration of the navigator in sLASER remained relatively short even with multiple shots (3T: 10 spatial interleaves 94 ms per slice; 7T: 15 spatial interleaves 103 ms per slice). Prospective motion correction using the multi-shot navigators yielded comparable spectral quality (water linewidth and metabolite SNR) with and without subject motion.

CONCLUSION

B0-insensitive motion navigators enable prospective motion correction for MRS with all first- and second-order shims adjusted in the MRS voxel, providing optimal spectral linewidth.

BigBrain-MR: a new digital phantom with anatomically-realistic magnetic resonance properties at 100-µm resolution for magnetic resonance methods development

The benefits, opportunities and growing availability of ultra-high field magnetic resonance imaging (MRI) for humans have prompted an expansion in research and development efforts towards increasingly more advanced high-resolution imaging techniques. To maximize their effectiveness, these efforts need to be supported by powerful computational simulation platforms that can adequately reproduce the biophysical characteristics of MRI, with high spatial resolution. In this work, we have sought to address this need by developing a novel digital phantom with realistic anatomical detail up to 100-µm resolution, including multiple MRI properties that affect image generation. This phantom, termed BigBrain-MR, was generated from the publicly available BigBrain histological dataset and lower-resolution in-vivo 7T-MRI data, using a newly-developed image processing framework that allows mapping the general properties of the latter into the fine anatomical scale of the former. Overall, the mapping framework was found to be effective and robust, yielding a diverse range of realistic "in-vivo-like" MRI contrasts and maps at 100-µm resolution. BigBrain-MR was then tested in three imaging applications (motion effects and interpolation, super-resolution imaging, and parallel imaging reconstruction) to investigate its properties, value and validity as a simulation platform. The results consistently showed that BigBrain-MR can closely approximate the behavior of real in-vivo data, more realistically and with more extensive features than a more classic option such as the Shepp-Logan phantom. Its flexibility in simulating different contrast mechanisms and artifacts may also prove valuable for educational applications. BigBrain-MR is therefore deemed a favorable choice to support methodological development and demonstration in brain MRI, and has been made freely available to the community.

Magnetic Resonance Imaging Studies of Neurodegenerative Disease: From Methods to Translational Research

Neurodegenerative diseases (NDs) have become a significant threat to an aging human society. Numerous studies have been conducted in the past decades to clarify their pathologic mechanisms and search for reliable biomarkers. Magnetic resonance imaging (MRI) is a powerful tool for investigating structural and functional brain alterations in NDs. With the advantages of being non-invasive and non-radioactive, it has been frequently used in both animal research and large-scale clinical investigations. MRI may serve as a bridge connecting micro- and macro-level analysis and promoting bench-to-bed translational research. Nevertheless, due to the abundance and complexity of MRI techniques, exploiting their potential is not always straightforward. This review aims to briefly introduce research progress in clinical imaging studies and discuss possible strategies for applying MRI in translational ND research.

Mesoscopic in vivo human T2* dataset acquired using quantitative MRI at 7 Tesla

Mesoscopic (0.1-0.5 mm) interrogation of the living human brain is critical for advancing neuroscience and bridging the resolution gap with animal models. Despite the variety of MRI contrasts measured in recent years at the mesoscopic scale, in vivo quantitative imaging of T₂^* has not been performed. Here we provide a dataset containing empirical T₂^* measurements acquired at 0.35 × 0.35 × 0.35 mm³ voxel resolution using 7 Tesla MRI. To demonstrate unique features and high quality of this dataset, we generate flat map visualizations that reveal fine-scale cortical substructures such as layers and vessels, and we report quantitative depth-dependent T₂^* (as well as R₂^*) values in primary visual cortex and auditory cortex that are highly consistent across subjects. This dataset is freely available at https://doi.org/10.17605/OSF.IO/N5BJ7, and may prove useful for anatomical investigations of the human brain, as well as for improving our understanding of the basis of the T₂^*-weighted (f)MRI signal.

Toward evaluation of multiresolution cortical thickness estimation with FreeSurfer, MaCRUISE, and BrainSuite

Advances in Magnetic Resonance Imaging hardware and methodologies allow for promoting the cortical morphometry with submillimeter spatial resolution. In this paper, we generated 3D self-enhanced high-resolution (HR) MRI imaging, by adapting 1 deep learning architecture, and 3 standard pipelines, FreeSurfer, MaCRUISE, and BrainSuite, have been collectively employed to evaluate the cortical thickness. We systematically investigated the differences in cortical thickness estimation for MRI sequences at multiresolution homologously originated from the native image. It has been revealed that there systematically exhibited the preferences in determining both inner and outer cortical surfaces at higher resolution, yielding most deeper cortical surface placements toward GM/WM or GM/CSF boundaries, which directs a consistent reduction tendency of mean cortical thickness estimation; on the contrary, the lower resolution data will most probably provide a more coarse and rough evaluation in cortical surface reconstruction, resulting in a relatively thicker estimation. Although the differences of cortical thickness estimation at the diverse spatial resolution varied with one another, almost all led to roughly one-sixth to one-fifth significant reduction across the entire brain at the HR, independent to the pipelines we applied, which emphasizes on generally coherent improved accuracy in a data-independent manner and endeavors to cost-efficiency with quantitative opportunities.

Neuroimaging at 7 Tesla: a pictorial narrative review

Neuroimaging using the 7-Tesla (7T) human magnetic resonance (MR) system is rapidly gaining popularity after being approved for clinical use in the European Union and the USA. This trend is the same for functional MR imaging (MRI). The primary advantages of 7T over lower magnetic fields are its higher signal-to-noise and contrast-to-noise ratios, which provide high-resolution acquisitions and better contrast, making it easier to detect lesions and structural changes in brain disorders. Another advantage is the capability to measure a greater number of neurochemicals by virtue of the increased spectral resolution. Many structural and functional studies using 7T have been conducted to visualize details in the white matter and layers of the cortex and hippocampus, the subnucleus or regions of the putamen, the globus pallidus, thalamus and substantia nigra, and in small structures, such as the subthalamic nucleus, habenula, perforating arteries, and the perivascular space, that are difficult to observe at lower magnetic field strengths. The target disorders for 7T neuroimaging range from tumoral diseases to vascular, neurodegenerative, and psychiatric disorders, including Alzheimer's disease, Parkinson's disease, multiple sclerosis, epilepsy, major depressive disorder, and schizophrenia. MR spectroscopy has also been used for research because of its increased chemical shift that separates overlapping peaks and resolves neurochemicals more effectively at 7T than a lower magnetic field. This paper presents a narrative review of these topics and an illustrative presentation of images obtained at 7T. We expect 7T neuroimaging to provide a new imaging biomarker of various brain disorders.

Tracking of rigid head motion during MRI using an EEG system

Purpose

To demonstrate a novel method for tracking of head movements during MRI using electroencephalography (EEG) hardware for recording signals induced by native imaging gradients.

Theory and Methods

Gradient switching during simultaneous EEG–fMRI induces distortions in EEG signals, which depend on subject head position and orientation. When EEG electrodes are interconnected with high‐impedance carbon wire loops, the induced voltages are linear combinations of the temporal gradient waveform derivatives. We introduce head tracking based on these signals (CapTrack) involving 3 steps: (1) phantom scanning is used to characterize the target sequence and a fast calibration sequence; (2) a linear relation between changes of induced signals and head pose is established using the calibration sequence; and (3) induced signals recorded during target sequence scanning are used for tracking and retrospective correction of head movement without prolonging the scan time of the target sequence. Performance of CapTrack is compared directly to interleaved navigators.

Results

Head‐pose tracking at 27.5 Hz during echo planar imaging (EPI) was demonstrated with close resemblance to rigid body alignment (mean absolute difference: [0.14 0.38 0.15]‐mm translation, [0.30 0.27 0.22]‐degree rotation). Retrospective correction of 3D gradient‐echo imaging shows an increase of average edge strength of 12%/−0.39% for instructed/uninstructed motion with CapTrack pose estimates, with a tracking interval of 1561 ms and high similarity to interleaved navigator estimates (mean absolute difference: [0.13 0.33 0.12] mm, [0.28 0.15 0.22] degrees).

Conclusion

Motion can be estimated from recordings of gradient switching with little or no sequence modification, optionally in real time at low computational burden and synchronized to image acquisition, using EEG equipment already found at many research institutions.

Automatic segmentation of Deep Grey Nuclei using a high-resolution 7T MRI Atlas - quantification of T1 values in healthy volunteers

We present a new consensus atlas of deep grey nuclei obtained by shape-based averaging of manual segmentation of two experienced neuroradiologists and optimized from 7T MP2RAGE images acquired at (0.6mm)³ in 60 healthy subjects. A group-wise normalization method was used to build a high-contrast and high-resolution T₁ -weighted brain template (0.5mm)³ using data from 30 out of the 60 controls. Delineation of 24 deep grey nuclei per hemisphere, including the claustrum and twelve thalamic nuclei, was then performed by two expert neuroradiologists and reviewed by a third neuroradiologist according to tissue contrast and external references based on the Morel atlas. Corresponding deep grey matter structures were also extracted from the Morel and CIT168 atlases. The data-derived, Morel and CIT168 atlases were all applied at the individual level using non-linear registration to fit the subject reference and to extract absolute mean quantitative T₁ values derived from the 3D-MP2RAGE volumes, after correction for residual B₁ ⁺ biases. Three metrics (The Dice and the volumetric similarity coefficients, and a novel Hausdorff distance) were used to estimate the inter-rater agreement of manual MRI segmentation and inter-atlas variability, and these metrics were measured to quantify biases due to image registration and their impact on the measurements of the quantitative T₁ values was highlighted. This represents a fully-automated segmentation process permitting the extraction of unbiased normative T₁ values in a population of young healthy controls as a reference for characterizing subtle structural alterations of deep grey nuclei relevant to a range of neurological diseases.

7 Tesla and Beyond: Advanced Methods and Clinical Applications in Magnetic Resonance Imaging

ABSTRACT

Ultrahigh magnetic fields offer significantly higher signal-to-noise ratio, and several magnetic resonance applications additionally benefit from a higher contrast-to-noise ratio, with static magnetic field strengths of B0 ≥ 7 T currently being referred to as ultrahigh fields (UHFs). The advantages of UHF can be used to resolve structures more precisely or to visualize physiological/pathophysiological effects that would be difficult or even impossible to detect at lower field strengths. However, with these advantages also come challenges, such as inhomogeneities applying standard radiofrequency excitation techniques, higher energy deposition in the human body, and enhanced B0 field inhomogeneities. The advantages but also the challenges of UHF as well as promising advanced methodological developments and clinical applications that particularly benefit from UHF are discussed in this review article.

High-resolution Structural Magnetic Resonance Imaging and Quantitative Susceptibility Mapping

High-resolution 7-T imaging and quantitative susceptibility mapping produce greater anatomic detail compared with conventional strengths because of improvements in signal/noise ratio and contrast. The exquisite anatomic details of deep structures, including delineation of microscopic architecture using advanced techniques such as quantitative susceptibility mapping, allows improved detection of abnormal findings thought to be imperceptible on clinical strengths. This article reviews caveats and techniques for translating sequences commonly used on 1.5 or 3 T to high-resolution 7-T imaging. It discusses for several broad disease categories how high-resolution 7-T imaging can advance the understanding of various diseases, improve diagnosis, and guide management.

Can 7T MPRAGE match MP2RAGE for gray-white matter contrast?

Ultra-High Field (UHF) MRI provides a significant increase in Signal-to-Noise Ratio (SNR) and gains in contrast weighting in several functional and structural acquisitions. Unfortunately, an increase in field strength also induces non-uniformities in the transmit field (B1+) that can compromise image contrast non-uniformly. The MPRAGE is one of the most common T1 weighted (T1w) image acquisitions for structural imaging. It provides excellent contrast between gray and white matter and is widely used for brain segmentation. At 7T, the signal non-uniformities tend to complicate this and therefore, the self-bias-field corrected MP2RAGE is often used there. In both MPRAGE and MP2RAGE, more homogeneous image contrast can be achieved with adiabatic pulses, like the TR-FOCI inversion pulse, or special pulse design on parallel transmission systems, like Universal Pulses (UP). In the present study, we investigate different strategies to improve the bias-field for MPRAGE at 7T, comparing the contrast and GM/WM segmentability against MP2RAGE. The higher temporal efficiency of MPRAGE combined with the potential of the user-friendly UPs was the primary motivation for this MPRAGE-MP2RAGE comparison. We acquired MPRAGE data in six volunteers, adding a k-space shutter to reduce scan time, a kt-point UP approach for homogeneous signal excitation, and a TR-FOCI pulse for homogeneous inversion. Our results show remarkable signal contrast improvement throughout the brain, including regions of low B1⁺ such as the cerebellum. The improvements in the MPRAGE were largest following the introduction of the UPs. In addition to the CNR, both SNR and GM/WM segmentability were also assessed. Among the MPRAGEs, the combined strategy (UP + TR-FOCI) yielded highest SNR and showed highest spatial similarity between GM segments to the MP2RAGE. Interestingly, the distance between gray and white matter peaks in the intensity histograms did not increase, as better pulses and higher SNR especially benefitted the (cerebellar) gray matter. Overall, the gray-white matter contrast from MP2RAGE is higher, with higher CNR and higher intensity peak distances, even when scaled to scan time. Hence, the extra acquisition time for MP2RAGE is justified by the improved segmentability.

Navigator-Guided Motion and B0 Correction of T2*-Weighted Magnetic Resonance Imaging Improves Multiple Sclerosis Cortical Lesion Detection

BACKGROUND

Cortical lesions are common in multiple sclerosis (MS). T2*-weighted (T2*w) imaging at 7 T is relatively sensitive for cortical lesions, but quality is often compromised by motion and main magnetic field (B0) fluctuations.

PURPOSE

The aim of this study was to determine whether motion and B0 correction with a navigator-guided gradient-recalled echo sequence can improve cortical lesion detection in T2*w magnetic resonance imaging.

MATERIALS AND METHODS

In this prospective study, a gradient-recalled echo sequence incorporating a navigator allowing for motion and B0 field correction was applied to collect T2*w images at 7 T from adults with MS between August 2019 and March 2020. T2*-weighted images were acquired in 1 to 3 partially overlapping scans per individual and were reconstructed using global average B0 correction ("uncorrected") or motion correction and spatially linear B0 correction ("corrected"). Image quality rating and manual segmentation of cortical lesions were performed on uncorrected and corrected images. Lesions seen on a single scan were retrospectively evaluated on the complementary scan. The association of cortical lesions with clinical disability was assessed. Mixed models were used to determine the effect of correction on lesion detection as well as on the relationship between disability and lesion count.

RESULTS

A total of 22 T2*w scans were performed on 11 adults with MS (mean [SD] age, 49 [11] years; 8 women). Quality improved for 20 of 22 scans (91%) after correction. A total of 69 cortical lesions were identified on uncorrected images (median per scan, 2; range, 0-11) versus 148 on corrected images (median per scan, 4.5; range, 0-25; rate ratio [RR], 2.1; P < 0.0001). For low-quality uncorrected scans with moderate to severe motion artifact (18/22, 82%), there was an improvement in cortical lesion detection with correction (RR, 2.5; P < 0.0001), whereas there was no significant change in cortical lesion detection for high-quality scans (RR, 1.3; P = 0.43).

CONCLUSIONS

Navigator-guided motion and B0 correction substantially improves the overall image quality of T2*w magnetic resonance imaging at 7 T and increases its sensitivity for cortical lesions.

Quantitative T1 mapping using multi-slice multi-shot inversion recovery EPI

An efficient multi-slice inversion–recovery EPI (MS-IR-EPI) sequence for fast, high spatial resolution, quantitative T₁ mapping is presented, using a segmented simultaneous multi-slice acquisition, combined with slice order shifting across multiple acquisitions. The segmented acquisition minimises the effective TE and readout duration compared to a single-shot EPI scheme, reducing geometric distortions to provide high quality T₁ maps with a narrow point-spread function. The precision and repeatability of MS-IR-EPI T₁ measurements are assessed using both T₁-calibrated and T₂-calibrated ISMRM/NIST phantom spheres at 3 and 7 T and compared with single slice IR and MP2RAGE methods. Magnetization transfer (MT) effects of the spectrally-selective fat-suppression (FS) pulses required for in vivo imaging are shown to shorten the measured in-vivo T₁ values. We model the effect of these fat suppression pulses on T₁ measurements and show that the model can remove their MT contribution from the measured T₁, thus providing accurate T₁ quantification. High spatial resolution T₁ maps of the human brain generated with MS-IR-EPI at 7 T are compared with those generated with the widely implemented MP2RAGE sequence. Our MS-IR-EPI sequence provides high SNR per unit time and sharper T₁ maps than MP2RAGE, demonstrating the potential for ultra-high resolution T₁ mapping and the improved discrimination of functionally relevant cortical areas in the human brain.

Quantity and quality: Normative open-access neuroimaging databases

The focus of this article is to compare twenty normative and open-access neuroimaging databases based on quantitative measures of image quality, namely, signal-to-noise (SNR) and contrast-to-noise ratios (CNR). We further the analysis through discussing to what extent these databases can be used for the visualization of deeper regions of the brain, such as the subcortex, as well as provide an overview of the types of inferences that can be drawn. A quantitative comparison of contrasts including T1-weighted (T1w) and T2-weighted (T2w) images are summarized, providing evidence for the benefit of ultra-high field MRI. Our analysis suggests a decline in SNR in the caudate nuclei with increasing age, in T1w, T2w, qT1 and qT2* contrasts, potentially indicative of complex structural age-dependent changes. A similar decline was found in the corpus callosum of the T1w, qT1 and qT2* contrasts, though this relationship is not as extensive as within the caudate nuclei. These declines were accompanied by a declining CNR over age in all image contrasts. A positive correlation was found between scan time and the estimated SNR as well as a negative correlation between scan time and spatial resolution. Image quality as well as the number and types of contrasts acquired by these databases are important factors to take into account when selecting structural data for reuse. This article highlights the opportunities and pitfalls associated with sampling existing databases, and provides a quantitative backing for their usage.

Perivascular Space Imaging at Ultrahigh Field MR Imaging

The recent Food and Drug Administration approval of 7 T MR imaging scanners for clinical use has introduced the possibility to study the brain not only in physiologic but also in pathologic conditions at ultrahigh field (UHF). Because UHF MR imaging offers higher signal-to-noise ratio and spatial resolution compared with lower field clinical scanners, the benefits of UHF MR imaging are particularly evident for imaging small anatomic structures, such as the cerebral perivascular spaces (PVS). In this article, the authors describe the application of UHF MR imaging for the investigation of PVS.

Sharpness in motion corrected quantitative imaging at 7T

Sub-millimeter imaging at 7T has opened new possibilities for qualitatively and quantitatively studying brain structure as it evolves throughout the life span. However, subject motion introduces image blurring on the order of magnitude of the spatial resolution and is thus detrimental to image quality. Such motion can be corrected for, but widespread application has not yet been achieved and quantitative evaluation is lacking. This raises a need to quantitatively measure image sharpness throughout the brain. We propose a method to quantify sharpness of brain structures at sub-voxel resolution, and use it to assess to what extent limited motion is related to image sharpness. The method was evaluated in a cohort of 24 healthy volunteers with a wide and uniform age range, aiming to arrive at results that largely generalize to larger populations. Using 3D fat-excited motion navigators, quantitative R₁, R₂^* and Quantitative Susceptibility Maps and T₁-weighted images were retrospectively corrected for motion. Sharpness was quantified in all modalities for selected regions of interest (ROI) by fitting the sigmoidally shaped error function to data within locally homogeneous clusters. A strong, almost linear correlation between motion and sharpness improvement was observed, and motion correction significantly improved sharpness. Overall, the Full Width at Half Maximum reduced from 0.88 mm to 0.70 mm after motion correction, equivalent to a 2.0 times smaller voxel volume. Motion and sharpness were not found to correlate with the age of study participants. We conclude that in our data, motion correction using fat navigators is overall able to restore the measured sharpness to the imaging resolution, irrespective of the amount of motion observed during scanning.

Spatial Resolution and the Magnitude of Infarct Volume Measurement Error in DWI in Acute Ischemic Stroke

BACKGROUND AND PURPOSE

Infarct volume in acute ischemic stroke is an important prognostic marker and determines endovascular treatment decisions. This study evaluates the magnitude and potential clinical impact of the error related to partial volume effects in infarct volume measurement on diffusion-weighted MR imaging in acute stroke and explores how increasing spatial resolution could reduce this error.

MATERIALS AND METHODS

Diffusion-weighted imaging of 393 patients with acute stroke, of whom 56 had anterior circulation large-vessel occlusion, was coregistered to standard space. Lesion boundaries were manually segmented. A 3D lesion-volume model was resampled for voxel sizes from 4 × 4 × 8 to 1 × 1 × 2 mm, and the surface-volume, corresponding to the partial volume error, was calculated. The number of cases with anterior circulation large-vessel occlusion, in which the endovascular therapy core threshold of 70 mL was contained within the margin of error, was calculated as a function of imaging resolution.

RESULTS

The mean infarct core volume was 27.2 ± 49.9 mL. The mean surface volume was 14.7 ± 20.8 mL for 2 × 2 × 4 mm resolution and 7.4 ± 10.7 mL for 1 × 1 × 2 mm resolution. With a resolution of 2 × 2 × 4 mm, 70 mL was contained within the margin of error in 7/56 cases (12.5%) with large-vessel occlusion, while with a 1 × 1 × 2 mm voxel size, the margin of error was 3/56 (5%). The lesion-volume range of potentially misclassified lesions dropped from 46.5-94.1 mL for a 2 × 2 × 4 mm resolution to 64.4-80.1 mL for a 1 × 1 × 2 mm resolution.

CONCLUSIONS

Partial volume effect is an important source of error in infarct volume measurement in acute stroke. Increasing spatial resolution substantially decreases the mean error. Standard use of high-resolution DWI should be considered to increase the reliability of infarct volume measurements.

A within-coil optical prospective motion-correction system for brain imaging at 7T

PURPOSE

Motion artifact limits the clinical translation of high-field MR. We present an optical prospective motion correction system for 7 Tesla MRI using a custom-built, within-coil camera to track an optical marker mounted on a subject.

METHODS

The camera was constructed to fit between the transmit-receive coils with direct line of sight to a forehead-mounted marker, improving upon prior mouthpiece work at 7 Tesla MRI. We validated the system by acquiring a 3D-IR-FSPGR on a phantom with deliberate motion applied. The same 3D-IR-FSPGR and a 2D gradient echo were then acquired on 7 volunteers, with/without deliberate motion and with/without motion correction. Three neuroradiologists blindly assessed image quality. In 1 subject, an ultrahigh-resolution 2D gradient echo with 4 averages was acquired with motion correction. Four single-average acquisitions were then acquired serially, with the subject allowed to move between acquisitions. A fifth single-average 2D gradient echo was acquired following subject removal and reentry.

RESULTS

In both the phantom and human subjects, deliberate and involuntary motion were well corrected. Despite marked levels of motion, high-quality images were produced without spurious artifacts. The quantitative ratings confirmed significant improvements in image quality in the absence and presence of deliberate motion across both acquisitions (P < .001). The system enabled ultrahigh-resolution visualization of the hippocampus during a long scan and robust alignment of serially acquired scans with interspersed movement.

CONCLUSION

We demonstrate the use of a within-coil camera to perform optical prospective motion correction and ultrahigh-resolution imaging at 7 Tesla MRI. The setup does not require a mouthpiece, which could improve accessibility of motion correction during 7 Tesla MRI exams.

MR spectroscopy using static higher order shimming with dynamic linear terms (HOS-DLT) for improved water suppression, interleaved MRS-fMRI, and navigator-based motion correction at 7T

PURPOSE

To interleave global and local higher order shimming for single voxel MRS. Single voxel MR spectroscopy requires optimization of the B₀ field homogeneity in the region of the voxel to obtain a narrow linewidth and provide high data quality. However, the optimization of local higher order fields on a localized MRS voxel typically leads to large field offsets outside that volume. This compromises interleaved MR sequence elements that benefit from global field homogeneity such as water suppression, interleaved MRS-fMRI, and MR motion correction.

METHODS

A shimming algorithm was developed to optimize the MRS voxel homogeneity and the whole brain homogeneity for interleaved sequence elements, using static higher order shims and dynamic linear terms (HOS-DLT). Shimming performance was evaluated using 6 brain regions and 10 subjects. Furthermore, the benefits of HOS-DLT was demonstrated for water suppression, MRS-fMRI, and motion corrected MRS using fat-navigators.

RESULTS

The HOS-DLT algorithm was shown to improve the whole brain homogeneity compared to an MRS voxel-based shim, without compromising the MRS voxel homogeneity. Improved water suppression over the brain, reduced image distortions in MRS-fMRI, and improved quality of motion navigators were demonstrated using the HOS-DLT method.

CONCLUSION

HOS-DLT shimming allowed for both local and global field homogeneity, providing excellent MR spectroscopy data quality, as well as good field homogeneity for interleaved sequence elements, even without the need for dynamic higher order shimming capabilities.

Synchrotron Radiation-Based Three-Dimensional Visualization of Angioarchitectural Remodeling in Hippocampus of Epileptic Rats

Characterizing the three-dimensional (3D) morphological alterations of microvessels under both normal and seizure conditions is crucial for a better understanding of epilepsy. However, conventional imaging techniques cannot detect microvessels on micron/sub-micron scales without angiography. In this study, synchrotron radiation (SR)-based X-ray in-line phase-contrast imaging (ILPCI) and quantitative 3D characterization were used to acquire high-resolution, high-contrast images of rat brain tissue under both normal and seizure conditions. The number of blood microvessels was markedly increased on days 1 and 14, but decreased on day 60 after seizures. The surface area, diameter distribution, mean tortuosity, and number of bifurcations and network segments also showed similar trends. These pathological changes were confirmed by histological tests. Thus, SR-based ILPCI provides systematic and detailed views of cerebrovascular anatomy at the micron level without using contrast-enhancing agents. This holds considerable promise for better diagnosis and understanding of the pathogenesis and development of epilepsy.

Denoising High-Field Multi-Dimensional MRI With Local Complex PCA

Modern high field and ultra high field magnetic resonance imaging (MRI) experiments routinely collect multi-dimensional data with high spatial resolution, whether multi-parametric structural, diffusion or functional MRI. While diffusion and functional imaging have benefited from recent advances in multi-dimensional signal analysis and denoising, structural MRI has remained untouched. In this work, we propose a denoising technique for multi-parametric quantitative MRI, combining a highly popular denoising method from diffusion imaging, over-complete local PCA, with a reconstruction of the complex-valued MR signal in order to define stable estimates of the noise in the decomposition. With this approach, we show signal to noise ratio (SNR) improvements in high resolution MRI without compromising the spatial accuracy or generating spurious perceptual boundaries.

Low-field MRI: An MR physics perspective

Historically, clinical MRI started with main magnetic field strengths in the ∼0.05-0.35T range. In the past 40 years there have been considerable developments in MRI hardware, with one of the primary ones being the trend to higher magnetic fields. While resulting in large improvements in data quality and diagnostic value, such developments have meant that conventional systems at 1.5 and 3T remain relatively expensive pieces of medical imaging equipment, and are out of the financial reach for much of the world. In this review we describe the current state-of-the-art of low-field systems (defined as 0.25-1T), both with respect to its low cost, low foot-print, and subject accessibility. Furthermore, we discuss how low field could potentially benefit from many of the developments that have occurred in higher-field MRI. In the first section, the signal-to-noise ratio (SNR) dependence on the static magnetic field and its impact on the achievable contrast, resolution, and acquisition times are discussed from a theoretical perspective. In the second section, developments in hardware (eg, magnet, gradient, and RF coils) used both in experimental low-field scanners and also those that are currently in the market are reviewed. In the final section the potential roles of new acquisition readouts, motion tracking, and image reconstruction strategies, currently being developed primarily at higher fields, are presented. Level of Evidence: 5 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2019.

Real-Time Filtering with Sparse Variations for Head Motion in Magnetic Resonance Imaging

Estimating a time-varying signal, such as head motion from magnetic resonance imaging data, becomes particularly challenging in the face of other temporal dynamics such as functional activation. This paper describes a new Kalman filter-like framework that includes a sparse residual term in the measurement model. This additional term allows the extended Kalman filter to generate real-time motion estimates suitable for prospective motion correction when such dynamics occur. An iterative augmented Lagrangian algorithm similar to the alterating direction method of multipliers implements the update step for this Kalman filter. This paper evaluates the accuracy and convergence rate of this iterative method for small and large motion in terms of its sensitivity to parameter selection. The included experiment on a simulated functional magnetic resonance imaging acquisition demonstrates that the resulting method improves the maximum Youden's J index of the time series analysis by 2-3% versus retrospective motion correction, while the sensitivity index increases from 4.3 to 5.4 when combining prospective and retrospective correction.

Anatomical and Technical Reappraisal of the Pallidothalamic Tractotomy With the Incisionless Transcranial MR-Guided Focused Ultrasound. A Technical Note

Background: MR-guided focused ultrasound (MRgFUS) offers new perspectives for safe and efficient lesioning inside the brain. The issue of target coverage remains primordial and sub-optimally addressed or solved in the field of functional neurosurgery. Objective: To provide an optimized planning and operative strategy to perform a pallidothalamic tractotomy (PTT) in chronic therapy-resistant Parkinson's disease (PD) with the technology of MRgFUS. Methods and results: Histological sections and maps from 6 human brain hemispheres were analyzed and outlines of the pallidothalamic tract on Myelin-stained sections were drawn and superimposed. We determined a standardized PTT target coverage characterized by 5 to 7 preplanned target lesion sub-units of 1.5 × 1.5 × 3.0 mm, which were placed using focal point displacements and shortest possible times, under thermal dose control. Conclusion: We hereby present our current approach to the MRgFUS PTT on the basis of a histological reappraisal and optimized heat application to the pallidothalamic tract in the H1 field of Forel.

Pros and cons of ultra-high-field MRI/MRS for human application

Magnetic resonance imaging and spectroscopic techniques are widely used in humans both for clinical diagnostic applications and in basic research areas such as cognitive neuroimaging. In recent years, new human MR systems have become available operating at static magnetic fields of 7 T or higher (≥300 MHz proton frequency). Imaging human-sized objects at such high frequencies presents several challenges including non-uniform radiofrequency fields, enhanced susceptibility artifacts, and higher radiofrequency energy deposition in the tissue. On the other side of the scale are gains in signal-to-noise or contrast-to-noise ratio that allow finer structures to be visualized and smaller physiological effects to be detected. This review presents an overview of some of the latest methodological developments in human ultra-high field MRI/MRS as well as associated clinical and scientific applications. Emphasis is given to techniques that particularly benefit from the changing physical characteristics at high magnetic fields, including susceptibility-weighted imaging and phase-contrast techniques, imaging with X-nuclei, MR spectroscopy, CEST imaging, as well as functional MRI. In addition, more general methodological developments such as parallel transmission and motion correction will be discussed that are required to leverage the full potential of higher magnetic fields, and an overview of relevant physiological considerations of human high magnetic field exposure is provided.

Prospective motion correction improves high-resolution quantitative susceptibility mapping at 7T

PURPOSE

Recent literature has shown the potential of high-resolution quantitative susceptibility mapping (QSM) with ultra-high field MRI for imaging the anatomy, the vasculature, and investigating their magnetostatic properties. Higher spatial resolutions, however, translate to longer scans resulting, therefore, in higher vulnerability to, and likelihood of, subject movement. We propose a gradient-recalled echo sequence with prospective motion correction (PMC) to address such limitation.

METHODS

Data from 4 subjects were acquired at 7T. The effect of small and large motion on QSM with and without PMC was assessed qualitatively and quantitatively. Full brain QSM and QSM-based venograms with up to 0.33 mm isotropic voxel size were reconstructed.

RESULTS

With PMC, motion artifacts in QSM and QSM-based venograms were largely eliminated, enabling-in both large- and small-amplitude motion regimes-accurate depiction of the cortex, vasculature, and other small anatomical structures that are often blurred as a result of head movement or indiscernible at lower image resolutions. Quantitative analyses demonstrated that uncorrected motion could bias regional susceptibility distributions, a trend that was greatly reduced with PMC.

CONCLUSION

Qualitatively, PMC prevented image degradation because of motion artifacts, providing highly detailed QSM images and venograms. Quantitatively, PMC increased the reproducibility of susceptibility measures.

Age differences in head motion and estimates of cortical morphology

Cortical morphology is known to differ with age, as measured by cortical thickness, fractal dimensionality, and gyrification. However, head motion during MRI scanning has been shown to influence estimates of cortical thickness as well as increase with age. Studies have also found task-related differences in head motion and relationships between body–mass index (BMI) and head motion. Here I replicated these prior findings, as well as several others, within a large, open-access dataset (Centre for Ageing and Neuroscience, CamCAN). This is a larger dataset than these results have been demonstrated previously, within a sample size of more than 600 adults across the adult lifespan. While replicating prior findings is important, demonstrating these key findings concurrently also provides an opportunity for additional related analyses: critically, I test for the influence of head motion on cortical fractal dimensionality and gyrification; effects were statistically significant in some cases, but small in magnitude.

Visualizing the Human Subcortex Using Ultra-high Field Magnetic Resonance Imaging

With the recent increased availability of ultra-high field (UHF) magnetic resonance imaging (MRI), substantial progress has been made in visualizing the human brain, which can now be done in extraordinary detail. This review provides an extensive overview of the use of UHF MRI in visualizing the human subcortex for both healthy and patient populations. The high inter-subject variability in size and location of subcortical structures limits the usability of atlases in the midbrain. Fortunately, the combined results of this review indicate that a large number of subcortical areas can be visualized in individual space using UHF MRI. Current limitations and potential solutions of UHF MRI for visualizing the subcortex are also discussed.

Diagnosis of multiple sclerosis through the lens of ultra-high-field MRI

The long-standing relationship between ultra-high-field (7 T) MRI and multiple sclerosis (MS) has brought new insights to our understanding of lesion evolution and its associated pathology. With the recent FDA approval of a commercially available scanner, 7 T MRI is finally entering the clinic with great expectations about its potential added value. By looking through the prism of MS diagnosis, this perspective article discusses current limitations and prospects of 7 T MRI techniques relevant to helping clinicians diagnose patients encountered in daily practice.

Analysis strategies for high-resolution UHF-fMRI data

Functional MRI (fMRI) benefits from both increased sensitivity and specificity with increasing magnetic field strength, making it a key application for Ultra-High Field (UHF) MRI scanners. Most UHF-fMRI studies utilize the dramatic increases in sensitivity and specificity to acquire high-resolution data reaching sub-millimeter scales, which enable new classes of experiments to probe the functional organization of the human brain. This review article surveys advanced data analysis strategies developed for high-resolution fMRI at UHF. These include strategies designed to mitigate distortion and artifacts associated with higher fields in ways that attempt to preserve spatial resolution of the fMRI data, as well as recently introduced analysis techniques that are enabled by these extremely high-resolution data. Particular focus is placed on anatomically-informed analyses, including cortical surface-based analysis, which are powerful techniques that can guide each step of the analysis from preprocessing to statistical analysis to interpretation and visualization. New intracortical analysis techniques for laminar and columnar fMRI are also reviewed and discussed. Prospects for single-subject individualized analyses are also presented and discussed. Altogether, there are both specific challenges and opportunities presented by UHF-fMRI, and the use of proper analysis strategies can help these valuable data reach their full potential.

Transcranial MRI-guided high-intensity focused ultrasound for treatment of essential tremor: A pilot study on the correlation between lesion size, lesion location, thermal dose, and clinical outcome

BACKGROUND

Transcranial MR-guided high-intensity focused ultrasound (tcMRgFUS) is a promising noninvasive method to treat medication-refractory essential tremor.

PURPOSE/HYPOTHESIS

To define the correlation between lesion size after ablation, thermal dose, and clinical outcome in tcMRgFUS treatment of essential tremor.

STUDY TYPE

Retrospective.

POPULATION/SUBJECTS/PHANTOM/SPECIMEN/ANIMAL MODEL

Eight patients with medication-refractory essential tremor were treated using a tcMRgFUS system at 3T.

FIELD STRENGTH/SEQUENCE

T₂ -weighted images were acquired immediately and at 1 year posttreatment at 3T.

ASSESSMENT

An atlas of the thalamic nuclei and dose maps were warped to the posttreatment images. The thermal dose, the immediate posttreatment lesion volume and 1-year final lesion volume, and the volumes confined inside the ventral division of the ventral lateral posterior thalamic nucleus (VLpv) were correlated to clinical outcome at 1 month and 1 year using Pearson's coefficient. The spatial region of treatment correlating with maximal clinical outcome was derived in a normalized space from average maps of clinical tremor score improvement at 1 year.

STATISTICAL TESTS

Statistical significance was assessed using the Wilcoxon two-tailed rank test.

RESULTS

The correlations between thermal dose, lesion volume posttreatment and at 1 year, and outcome at 1 year were good (r = 0.73, 0.65, 0.73, respectively), and were slightly better than at 1 month (r = 0.57, 0.49, 0.65). Reducing the measurement to include only the portion within the VLpv did not significantly modify the correlations (P = 0.09). The center of the spatial region of treatment was found in the anterior commissure - posterior commissure plane, 14.3 mm lateral from the midline, and 8.3 mm rostral to the posterior commissure.

DATA CONCLUSION

In this pilot study a good correlation was found between the size of the lesion, the thermal dose, and the clinical outcome in patients treated for essential tremor with ablation of the VLpv with tcMRgFUS.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 4 J. Magn. Reson. Imaging 2017.

Advantages of cortical surface reconstruction using submillimeter 7 T MEMPRAGE

Recent advances in MR technology have enabled increased spatial resolution for routine functional and anatomical imaging, which has created demand for software tools that are able to process these data. The availability of high-resolution data also raises the question of whether higher resolution leads to substantial gains in accuracy of quantitative morphometric neuroimaging procedures, in particular the cortical surface reconstruction and cortical thickness estimation. In this study we adapted the FreeSurfer cortical surface reconstruction pipeline to process structural data at native submillimeter resolution. We then quantified the differences in surface placement between meshes generated from (0.75 mm)³ isotropic resolution data acquired in 39 volunteers and the same data downsampled to the conventional 1 mm³ voxel size. We find that when processed at native resolution, cortex is estimated to be thinner in most areas, but thicker around the Cingulate and the Calcarine sulci as well as in the posterior bank of the Central sulcus. Thickness differences are driven by two kinds of effects. First, the gray-white surface is found closer to the white matter, especially in cortical areas with high myelin content, and thus low contrast, such as the Calcarine and the Central sulci, causing local increases in thickness estimates. Second, the gray-CSF surface is placed more interiorly, especially in the deep sulci, contributing to local decreases in thickness estimates. We suggest that both effects are due to reduced partial volume effects at higher spatial resolution. Submillimeter voxel sizes can therefore provide improved accuracy for measuring cortical thickness.

The impact of vessel size, orientation and intravascular contribution on the neurovascular fingerprint of BOLD bSSFP fMRI

Monte Carlo simulations have been used to analyze oxygenation-related signal changes in pass-band balanced steady state free precession (bSSFP) as well as in gradient echo (GE) and spin echo (SE) sequences. Signal changes were calculated for artificial cylinders and neurovascular networks acquired from the mouse parietal cortex by two-photon laser scanning microscopy at 1 μm isotropic resolution. Signal changes as a function of vessel size, blood volume, vessel orientation to the main magnetic field B₀ as well as relations of intra- and extravascular and of micro- and macrovascular contributions have been analyzed. The results show that bSSFP is highly sensitive to extravascular and microvascular components. Furthermore, GE and bSSFP, and to a lesser extent SE, exhibit a strong dependence of their signal change on the orientation of the vessel network to B₀.

Effects of aging on T₁, T₂*, and QSM MRI values in the subcortex

The aging brain undergoes several anatomical changes that can be measured with Magnetic Resonance Imaging (MRI). Early studies using lower field strengths have assessed changes in tissue properties mainly qualitatively, using [Formula: see text]- or [Formula: see text]- weighted images to provide image contrast. With the development of higher field strengths (7 T and above) and more advanced MRI contrasts, quantitative measures can be acquired even of small subcortical structures. This study investigates volumetric, spatial, and quantitative MRI parameter changes associated with healthy aging in a range of subcortical nuclei, including the basal ganglia, red nucleus, and the periaqueductal grey. The results show that aging has a heterogenous effects across regions. Across the subcortical areas an increase of [Formula: see text] values is observed, most likely indicating a loss of myelin. Only for a number of areas, a decrease of [Formula: see text] and increase of QSM is found, indicating an increase of iron. Aging also results in a location shift for a number of structures indicating the need for visualization of the anatomy of individual brains.

Test-retest reliability of brain morphology estimates

Metrics of brain morphology are increasingly being used to examine inter-individual differences, making it important to evaluate the reliability of these structural measures. Here we used two open-access datasets to assess the intersession reliability of three cortical measures (thickness, gyrification, and fractal dimensionality) and two subcortical measures (volume and fractal dimensionality). Reliability was generally good, particularly with the gyrification and fractal dimensionality measures. One dataset used a sequence previously optimized for brain morphology analyses and had particularly high reliability. Examining the reliability of morphological measures is critical before the measures can be validly used to investigate inter-individual differences.

T1-weighted in vivo human whole brain MRI dataset with an ultrahigh isotropic resolution of 250 μm

We present an ultrahigh resolution in vivo human brain magnetic resonance imaging (MRI) dataset. It consists of T₁-weighted whole brain anatomical data acquired at 7 Tesla with a nominal isotropic resolution of 250 μm of a single young healthy Caucasian subject and was recorded using prospective motion correction. The raw data amounts to approximately 1.2 TB and was acquired in eight hours total scan time. The resolution of this dataset is far beyond any previously published in vivo structural whole brain dataset. Its potential use is to build an in vivo MR brain atlas. Methods for image reconstruction and image restoration can be improved as the raw data is made available. Pre-processing and segmentation procedures can possibly be enhanced for high magnetic field strength and ultrahigh resolution data. Furthermore, potential resolution induced changes in quantitative data analysis can be assessed, e.g., cortical thickness or volumetric measures, as high quality images with an isotropic resolution of 1 and 0.5 mm of the same subject are included in the repository as well.

Insights from Ultrahigh Field Imaging in Multiple Sclerosis

Ultrahigh-field (≥7 T) magnetic resonance (MR) imaging is being used at many leading academic medical centers to study neurologic disorders. The improved spatial resolution and anatomic detail are due to the increase in signal-to-noise and contrast-to-noise ratio at higher magnetic field strengths. Ultrahigh-field MR imaging improves multiple sclerosis (MS) lesion detection, with particular sensitivity to detect cortical lesions. The increase in magnetic susceptibility effects inherent to ultrahigh field can be used to detect pathologic features of MS lesions, including a central vein, potentially useful for diagnostic considerations, and heterogeneity among MS lesions, potentially useful in determining lesion outcomes.