Functional magnetic resonance imaging (fMRI) has become an indispensable tool for investigating the human brain. However, the inherently poor signal-to-noise-ratio (SNR) of the fMRI measurement represents a major barrier to expanding its spatiotemporal scale as well as its utility and ultimate impact. Here we introduce a denoising technique that selectively suppresses the thermal noise contribution to the fMRI experiment. Using 7-Tesla, high-resolution human brain data, we demonstrate improvements in key metrics of functional mapping (temporal-SNR, the detection and reproducibility of stimulus-induced signal changes, and accuracy of functional maps) while leaving the amplitude of the stimulus-induced signal changes, spatial precision, and functional point-spread-function unaltered. We demonstrate that the method enables the acquisition of ultrahigh resolution (0.5 mm isotropic) functional maps but is also equally beneficial for a large variety of fMRI applications, including supra-millimeter resolution 3- and 7-Tesla data obtained over different cortical regions with different stimulation/task paradigms and acquisition strategies.

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.

Determining the acquisition parameters in diffusion magnetic resonance imaging (dMRI) is governed by a series of trade-offs. Images of lower resolution have less spatial specificity but higher signal to noise ratio (SNR). At the same time higher angular contrast, important for resolving complex fibre patterns, also yields lower SNR. Considering these trade-offs, the Human Connectome Project (HCP) acquires high quality dMRI data for the same subjects at different field strengths (3T and 7T), which are publically released. Due to differences in the signal behavior and in the underlying scanner hardware, the HCP 3T and 7T data have complementary features in k- and q-space. The 3T dMRI has higher angular contrast and resolution, while the 7T dMRI has higher spatial resolution. Given the availability of these datasets, we explore the idea of fusing them together with the aim of combining their benefits. We extend a previously proposed data-fusion framework and apply it to integrate both datasets from the same subject into a single joint analysis. We use a generative model for performing parametric spherical deconvolution and estimate fibre orientations by simultaneously using data acquired under different protocols. We illustrate unique features from each dataset and how they are retained after fusion. We further show that this allows us to complement benefits and improve brain connectivity analysis compared to analyzing each of the datasets individually.

To increase granularity in human neuroimaging science, we designed and built a next-generation 7 Tesla magnetic resonance imaging scanner to reach ultra-high resolution by implementing several advances in hardware. To improve spatial encoding and increase the image signal-to-noise ratio, we developed a head-only asymmetric gradient coil (200 mT m-1, 900 T m-1s-1) with an additional third layer of windings. We integrated a 128-channel receiver system with 64- and 96-channel receiver coil arrays to boost signal in the cerebral cortex while reducing g-factor noise to enable higher accelerations. A 16-channel transmit system reduced power deposition and improved image uniformity. The scanner routinely performs functional imaging studies at 0.35-0.45 mm isotropic spatial resolution to reveal cortical layer functional activity, achieves high angular resolution in diffusion imaging and reduces acquisition time for both functional and structural imaging.

Functional magnetic resonance imaging (fMRI) has emerged as an essential tool for exploring human brain function. Submillimeter fMRI, in particular, has emerged as a tool to study mesoscopic computations. The inherently low signal-to-noise ratio (SNR) at submillimeter resolutions warrants the use of denoising approaches tailored at reducing thermal noise - the dominant contributing noise component in high resolution fMRI. NORDIC PCA is one of such approaches, and has been benchmarked against other approaches in several applications. Here, we investigate the effects that two versions of NORDIC denoising have on auditory submillimeter data. As investigating auditory functional responses poses unique challenges, we anticipated that the benefit of this technique would be especially pronounced. Our results show that NORDIC denoising improves the detection sensitivity and the reliability of estimates in submillimeter auditory fMRI data. These effects can be explained by the reduction of the noise-induced signal variability. However, we also observed a reduction in the average response amplitude (percent signal), which may suggest that a small amount of signal was also removed. We conclude that, while evaluating the effects of the signal reduction induced by NORDIC may be necessary for each application, using NORDIC in high resolution auditory fMRI studies may be advantageous because of the large reduction in variability of the estimated responses.

Functional magnetic resonance imaging (fMRI) allows studying human brain function non-invasively up to the spatial resolution of cortical columns and layers. Most fMRI acquisitions rely on the blood oxygenation level dependent (BOLD) contrast employing T(*) 2 weighted 2D multi-slice echo-planar imaging (EPI). At ultra-high magnetic field (i.e., 7 T and above), it has been shown experimentally and by simulation, that T2 weighted acquisitions yield a signal that is spatially more specific to the site of neuronal activity at the cost of functional sensitivity. This study compared two T2 weighted imaging sequences, inner-volume 3D Gradient-and-Spin-Echo (3D-GRASE) and 2D Spin-Echo EPI (SE-EPI), with evaluation of their imaging point-spread function (PSF), functional specificity, and functional sensitivity at sub-millimeter resolution. Simulations and measurements of the imaging PSF revealed that the strongest anisotropic blurring in 3D-GRASE (along the second phase-encoding direction) was about 60% higher than the strongest anisotropic blurring in 2D SE-EPI (along the phase-encoding direction). In a visual paradigm, the BOLD sensitivity of 3D-GRASE was found to be superior due to its higher temporal signal-to-noise ratio (tSNR). High resolution cortical depth profiles suggested that the contrast mechanisms are similar between the two sequences, however, 2D SE-EPI had a higher surface bias owing to the higher T(*) 2 contribution of the longer in-plane EPI echo-train for full field of view compared to the reduced field of view of zoomed 3D-GRASE.

The quest to map brain connectivity is being pursued worldwide using diffusion imaging, among other techniques. Even so, we know little about how brain connectivity measures depend on the magnetic field strength of the scanner. To investigate this, we scanned 10 healthy subjects at 7 and 3 tesla-using 128-gradient high-angular resolution diffusion imaging. For each subject and scan, whole-brain tractography was used to estimate connectivity between 113 cortical and subcortical regions. We examined how scanner field strength affects (i) the signal-to-noise ratio (SNR) of the non-diffusion-sensitized reference images (b(0)); (ii) diffusion tensor imaging (DTI)-derived fractional anisotropy (FA), mean, radial, and axial diffusivity (MD/RD/AD), in atlas-defined regions; (iii) whole-brain tractography; (iv) the 113 × 113 brain connectivity maps; and (v) five commonly used network topology measures. We also assessed effects of the multi-channel reconstruction methods (sum-of-squares, SOS, at 7T; adaptive recombine, AC, at 3T). At 7T with SOS, the b0 images had 18.3% higher SNR than with 3T-AC. FA was similar for most regions of interest (ROIs) derived from an online DTI atlas (ICBM81), but higher at 7T in the cerebral peduncle and internal capsule. MD, AD, and RD were lower at 7T for most ROIs. The apparent fiber density between some subcortical regions was greater at 7T-SOS than 3T-AC, with a consistent connection pattern overall. Suggesting the need for caution, the recovered brain network was apparently more efficient at 7T, which cannot be biologically true as the same subjects were assessed. Care is needed when comparing network measures across studies, and when interpreting apparently discrepant findings.

Simultaneous multi-slice (SMS) or multiband (MB) imaging has recently been attempted for arterial spin labeled (ASL) perfusion MRI in conjunction with echo-planar imaging (EPI) readout. It was found that SMS-EPI can reduce the T1 relaxation effect of the label and improve image coverage and resolution with little penalty in signal-to-noise ratio (SNR). However, EPI still suffers from geometric distortion and signal dropout from field inhomogeneity effects especially at high and ultrahigh magnetic fields. Here we present a novel scheme for achieving high fidelity distortion-free quantitative perfusion imaging by combining pseudo-continuous ASL (pCASL) with SMS Turbo-FLASH (TFL) readout at both 3 and 7 T. Bloch equation simulation was performed to characterize and optimize the TFL-based pCASL perfusion signal. Two MB factors (3 and 5) were implemented in SMS-TFL pCASL and compared with standard 2D TFL and EPI pCASL sequences. The temporal SNR of SMS-TFL pCASL relative to that of standard TFL pCASL was 0.76 ± 0.10 and 0.74 ± 0.11 at 7 T and 0.70 ± 0.05 and 0.65 ± 0.05 at 3T for MB factor of 3 and 5, respectively. By implementing background suppression in conjunction with SMS-TFL at 3T, the relative temporal SNR improved to 0.84 ± 0.09 and 0.79 ± 0.10 for MB factor of 3 and 5, respectively. Compared to EPI pCASL, significantly increased temporal SNR (p<0.001) and improved visualization of orbitofrontal cortex were achieved using SMS-TFL pCASL. By combining SMS acceleration with TFL pCASL, we demonstrated the feasibility for whole brain distortion-free quantitative mapping of cerebral blood flow at high and ultrahigh magnetic fields.

Basal ganglia circuits are affected in neurological disorders such as Parkinson's disease (PD), essential tremor, dystonia and Tourette syndrome. Understanding the structural and functional connectivity of these circuits is critical for elucidating the mechanisms of the movement and neuropsychiatric disorders, and is vital for developing new therapeutic strategies such as deep brain stimulation (DBS). Knowledge about the connectivity of the human basal ganglia and thalamus has rapidly evolved over recent years through non-invasive imaging techniques, but has remained incomplete because of insufficient resolution and sensitivity of these techniques. Here, we present an imaging and computational protocol designed to generate a comprehensive in vivo and subject-specific, three-dimensional model of the structure and connections of the human basal ganglia. High-resolution structural and functional magnetic resonance images were acquired with a 7-Tesla magnet. Capitalizing on the enhanced signal-to-noise ratio (SNR) and enriched contrast obtained at high-field MRI, detailed structural and connectivity representations of the human basal ganglia and thalamus were achieved. This unique combination of multiple imaging modalities enabled the in-vivo visualization of the individual human basal ganglia and thalamic nuclei, the reconstruction of seven white-matter pathways and their connectivity probability that, to date, have only been reported in animal studies, histologically, or group-averaged MRI population studies. Also described are subject-specific parcellations of the basal ganglia and thalamus into sub-territories based on their distinct connectivity patterns. These anatomical connectivity findings are supported by functional connectivity data derived from resting-state functional MRI (R-fMRI). This work demonstrates new capabilities for studying basal ganglia circuitry, and opens new avenues of investigation into the movement and neuropsychiatric disorders, in individual human subjects.

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.