Tradeoffs in pushing the spatial resolution of fMRI for the 7T Human Connectome Project

7 T and beyond: toward a synergy between fMRI-based presurgical mapping at ultrahigh magnetic fields, AI, and robotic neurosurgery

Presurgical evaluation with functional magnetic resonance imaging (fMRI) can reduce postsurgical morbidity. Here, we discuss presurgical fMRI mapping at ultra-high magnetic fields (UHF), i.e., ≥ 7 T, in the light of the current growing interest in artificial intelligence (AI) and robot-assisted neurosurgery. The potential of submillimetre fMRI mapping can help better appreciate uncertainty on resection margins, though geometric distortions at UHF might lessen the accuracy of fMRI maps. A useful trade-off for UHF fMRI is to collect data with 1-mm isotropic resolution to ensure high sensitivity and subsequently a low risk of false negatives. Scanning at UHF might yield a revival interest in slow event-related fMRI, thereby offering a richer depiction of the dynamics of fMRI responses. The potential applications of AI concern denoising and artefact removal, generation of super-resolution fMRI maps, and accurate fusion or coregistration between anatomical and fMRI maps. The latter can benefit from the use of T1-weighted echo-planar imaging for better visualization of brain activations. Such AI-augmented fMRI maps would provide high-quality input data to robotic surgery systems, thereby improving the accuracy and reliability of robot-assisted neurosurgery. Ultimately, the advancement in fMRI at UHF would promote clinically useful synergies between fMRI, AI, and robotic neurosurgery.Relevance statement This review highlights the potential synergies between fMRI at UHF, AI, and robotic neurosurgery in improving the accuracy and reliability of fMRI-based presurgical mapping.Key points• Presurgical fMRI mapping at UHF improves spatial resolution and sensitivity.• Slow event-related designs offer a richer depiction of fMRI responses dynamics.• AI can support denoising, artefact removal, and generation of super-resolution fMRI maps.• AI-augmented fMRI maps can provide high-quality input data to robotic surgery systems.

RF coil design strategies for improving SNR at the ultrahigh magnetic field of 10.5 Tesla

Purpose

To develop multichannel transmit and receive arrays towards capturing the ultimate-intrinsic-SNR (uiSNR) at 10.5 Tesla (T) and to demonstrate the feasibility and potential of whole-brain, high-resolution human brain imaging at this high field strength.

Methods

A dual row 16-channel self-decoupled transmit (Tx) array was converted to a 16Tx/Rx transceiver using custom transmit/receive switches. A 64-channel receive-only (64Rx) array was built to fit into the 16Tx/Rx array. Electromagnetic modeling and experiments were employed to define safe operation limits of the resulting 16Tx/80Rx array and obtain FDA approval for human use.

Results

The 64Rx array alone captured approximately 50% of the central uiSNR at 10.5T while the identical 7T 64Rx array captured ∼76% of uiSNR at this lower field strength. The 16Tx/80Rx configuration brought the fraction of uiSNR captured at 10.5T to levels comparable to the performance of the 64Rx array at 7T. SNR data obtained at the two field strengths with these arrays displayed dependent increases over a large central region. Whole-brain high resolution T 2 * and T 1 weighted anatomical and gradient-recalled echo EPI BOLD fMRI images were obtained at 10.5T for the first time with such an advanced array, illustrating the promise of >10T fields in studying the human brain.

Conclusion

We demonstrated the ability to approach the uiSNR at 10.5T over the human brain with a novel, high channel count array, achieving large SNR gains over 7T, currently the most commonly employed ultrahigh field platform, and demonstrate high resolution and high contrast anatomical and functional imaging at 10.5T.

Functional Localization of the Human Auditory and Visual Thalamus Using a Thalamic Localizer Functional Magnetic Resonance Imaging Task

Functional magnetic resonance imaging (fMRI) of the auditory and visual sensory systems of the human brain is an active area of investigation in the study of human health and disease. The medial geniculate nucleus (MGN) and lateral geniculate nucleus (LGN) are key thalamic nuclei involved in the processing and relay of auditory and visual information, respectively, and are the subject of blood-oxygen-level-dependent (BOLD) fMRI studies of neural activation and functional connectivity in human participants. However, localization of BOLD fMRI signal originating from neural activity in MGN and LGN remains a technical challenge, due in part to the poor definition of boundaries of these thalamic nuclei in standard T1-weighted and T2-weighted magnetic resonance imaging sequences. Here, we report the development and evaluation of an auditory and visual sensory thalamic localizer (TL) fMRI task that produces participant-specific functionally-defined regions of interest (fROIs) of both MGN and LGN, using 3 Tesla multiband fMRI and a clustered-sparse temporal acquisition sequence, in less than 16 minutes of scan time. We demonstrate the use of MGN and LGN fROIs obtained from the TL fMRI task in standard resting-state functional connectivity (RSFC) fMRI analyses in the same participants. In RSFC analyses, we validated the specificity of MGN and LGN fROIs for signals obtained from primary auditory and visual cortex, respectively, and benchmark their performance against alternative atlas- and segmentation-based localization methods. The TL fMRI task and analysis code (written in Presentation and MATLAB, respectively) have been made freely available to the wider research community.

Disentangling neuroplasticity mechanisms in post-stroke language recovery

A major objective in post-stroke aphasia research is to gain a deeper understanding of neuroplastic mechanisms that drive language recovery, with the ultimate goal of enhancing treatment outcomes. Subsequent to recent advances in neuroimaging techniques, we now have the ability to examine more closely how neural activity patterns change after a stroke. However, the way these neural activity changes relate to language impairments and language recovery is still debated. The aim of this review is to provide a theoretical framework to better investigate and interpret neuroplasticity mechanisms underlying language recovery in post-stroke aphasia. We detail two sets of neuroplasticity mechanisms observed at the synaptic level that may explain functional neuroimaging findings in post-stroke aphasia recovery at the network level: feedback-based homeostatic plasticity and associative Hebbian plasticity. In conjunction with these plasticity mechanisms, higher-order cognitive control processes dynamically modulate neural activity in other regions to meet communication demands, despite reduced neural resources. This work provides a network-level neurobiological framework for understanding neural changes observed in post-stroke aphasia and can be used to define guidelines for personalized treatment development.

Frequency-specific segregation and integration of human cerebral cortex: An intrinsic functional atlas

The cognitive and behavioral functions of the human brain are supported by its frequency multiplexing mechanism. However, there is limited understanding of the dynamics of the functional network topology. This study aims to investigate the frequency-specific topology of the functional human brain using 7T rs-fMRI data. Frequency-specific parcellations were first performed, revealing frequency-dependent dynamics within the frontoparietal control, parietal memory, and visual networks. An intrinsic functional atlas containing 456 parcels was proposed and validated using stereo-EEG. Graph theory analysis suggested that, in addition to the task-positive vs. task-negative organization observed in static networks, there was a cognitive control system additionally from a frequency perspective. The reproducibility and plausibility of the identified hub sets were confirmed through 3T fMRI analysis, and their artificial removal had distinct effects on network topology. These results indicate a more intricate and subtle dynamics of the functional human brain and emphasize the significance of accurate topography.

Early excitatory-inhibitory cortical modifications following skill learning are associated with motor memory consolidation and plasticity overnight

Consolidation of motor memories is vital to offline enhancement of new motor skills and involves short and longer-term offline processes following learning. While emerging evidence link glutamate and GABA dynamics in the primary motor cortex (M1) to online motor skill practice, its relationship with offline consolidation processes in humans is unclear. Using two-day repeated measures of behavioral and multimodal neuroimaging data before and following motor sequence learning, we show that short-term glutamatergic and GABAergic responses in M1 within minutes after learning were associated with longer-term learning-induced functional, structural, and behavioral modifications overnight. Furthermore, Glutamatergic and GABAergic modifications were differentially associated with different facets of motor memory consolidation. Our results point to unique and distinct roles of Glutamate and GABA in motor memory consolidation processes in the human brain across timescales and mechanistic levels, tying short-term changes on the neurochemical level to overnight changes in macroscale structure, function, and behavior.

The current state of spinal cord functional magnetic resonance imaging and its application in clinical research

Since its development, spinal cord functional magnetic resonance imaging (fMRI) has utilized various methodologies and stimulation protocols to develop a deeper understanding of a healthy human spinal cord that lays a foundation for its use in clinical research and practice. In this review, we conducted a comprehensive literature search on spinal cord fMRI studies and summarized the recent advancements and resulting scientific achievements of spinal cord fMRI in the following three aspects: the current state of spinal cord fMRI methodologies and stimulation protocols, knowledge about the healthy spinal cord's functions obtained via spinal cord fMRI, and fMRI's exemplary usage in spinal cord diseases and injuries. We conclude with a discussion that, while technical challenges exist, novel fMRI technologies for and new knowledge about the healthy human spinal cord have been established. Empowered by these developments, investigations of pathological and injury states within the spinal cord have become the next important direction of spinal cord fMRI. Recent clinical investigations into spinal cord pathologies, for example, fibromyalgia, multiple sclerosis, spinal cord injury, and cervical spondylotic myelopathy, have already provided deep insights into spinal cord impairments and the time course of impairment-caused changes. We expect that future spinal cord fMRI advancement and research development will further enhance our understanding of various spinal cord diseases and provide the foundation for evaluating existing and developing new treatment plans.

Evaluating functional brain organization in individuals and identifying contributions to network overlap

Individual differences in the spatial organization of resting state networks have received increased attention in recent years. Measures of individual-specific spatial organization of brain networks and overlapping network organization have been linked to important behavioral and clinical traits and are therefore potential biomarker targets for personalized psychiatry approaches. To better understand individual-specific spatial brain organization, this paper addressed three key goals. First, we determined whether it is possible to reliably estimate weighted (non-binarized) resting state network maps using data from only a single individual, while also maintaining maximum spatial correspondence across individuals. Second, we determined the degree of spatial overlap between distinct networks, using test-retest and twin data. Third, we systematically tested multiple hypotheses (spatial mixing, temporal switching, and coupling) as candidate explanations for why networks overlap spatially. To estimate weighted network organization, we adopt the Probabilistic Functional Modes (PROFUMO) algorithm, which implements a Bayesian framework with hemodynamic and connectivity priors to supplement optimization for spatial sparsity/independence. Our findings showed that replicable individual-specific estimates of weighted resting state networks can be derived using high quality fMRI data within individual subjects. Network organization estimates using only data from each individual subject closely resembled group-informed network estimates (which was not explicitly modeled in our individual-specific analyses), suggesting that cross-subject correspondence was largely maintained. Furthermore, our results confirmed the presence of spatial overlap in network organization, which was replicable across sessions within individuals and in monozygotic twin pairs. Intriguingly, our findings provide evidence that network overlap is indicative of linear additive coupling. These results suggest that regions of network overlap concurrently process information from all contributing networks, potentially pointing to the role of overlapping network organization in the integration of information across multiple brain systems.

Individual Variability of Human Cortical Spontaneous Activity by 3T/7T fMRI

Mapping variability in cortical spontaneous activity (CSA) is an essential goal of understanding various sources of dark brain energy in human neuroscience. CSA was traditionally characterized using resting-state functional MRI (rfMRI) at 1.5T or 3T magnets while recently with 7T-rfMRI. However, the utility and interpretability of 7T-rfMRI must first be established for its variability. By leveraging rfMRI data from the Human Connectome Project (HCP), we derived CSA metrics with 3T-rfMRI and 7T-rfMRI for the same 84 healthy participants (52 females). The 7T-rfMRI produces different CSA metrics at multiple spatial-scales and their variability from the 3T-rfMRI. These differences were spatially dependent and varied according to specific cortical organization. For the amplitude metric, 7T-rfMRI enhanced its spatial contrasts in the anterior cortex but weakened it in the posterior cortex. An opposite pattern was observed for the connectivity metrics. The reliability changes of these metrics were scale dependent, indicating enhanced reliability for connectivity but weakened reliability for amplitude by 7T-rfMRI. These effects were primarily located in the high-order associate cortex, parsing the corresponding changes in individual differences with respect to 7T-rfMRI: (1) higher connectivity variability between participants and the lower connectivity variability within individual participants, and (2) lower amplitude variability between participants and higher amplitude variability within participants. Our work, for the first time, demonstrated the variability of the human CSA across space, rfMRI settings/platforms, and individuals. We discussed the statistical implications of our findings on CSA-based experimental designs and reproducible neuroscience as well as their translational value for personalized applications.

Evaluation of consciousness rehabilitation via neuroimaging methods

Accurate evaluation of patients with disorders of consciousness (DoC) is crucial for personalized treatment. However, misdiagnosis remains a serious issue. Neuroimaging methods could observe the conscious activity in patients who have no evidence of consciousness in behavior, and provide objective and quantitative indexes to assist doctors in their diagnosis. In the review, we discussed the current research based on the evaluation of consciousness rehabilitation after DoC using EEG, fMRI, PET, and fNIRS, as well as the advantages and limitations of each method. Nowadays single-modal neuroimaging can no longer meet the researchers` demand. Considering both spatial and temporal resolution, recent studies have attempted to focus on the multi-modal method which can enhance the capability of neuroimaging methods in the evaluation of DoC. As neuroimaging devices become wireless, integrated, and portable, multi-modal neuroimaging methods will drive new advancements in brain science research.

Resting-state functional connectivity-based parcellation of the human dentate nucleus: new findings and clinical relevance

For years, the cerebellum was left out of functional magnetic resonance imaging (fMRI) studies due to technological limitations. The advent of novel data acquisition and reconstruction strategies (e.g., whole-brain simultaneous multi-slice imaging) employing multi-channel array coils has overcome such limitations, ushering unprecedented improvements in temporal signal-to-noise ratio and spatiotemporal resolution. Here, we aim to provide a brief report on the deep cerebellar nuclei, specifically focusing on the dentate nuclei, the primary output nuclei, situated within both cognitive and motor cerebello-cerebral circuits. We highlight the importance of functional parcellation in refining our understanding of broad resting-state functional connectivity (RSFC) in both health and disease. First, we review work relevant to the functional topography of the dentate nuclei, including recent advances in functional parcellation. Next, we review RSFC studies using the dentate nuclei as seed regions of interest in neurological and psychiatric populations and discuss the potential benefits of applying functionally defined subdivisions. Finally, we discuss recent technological advances and underscore ultrahigh-field neuroimaging as a tool to potentiate functionally parcellated RSFC analyses in clinical populations.

Age-related differences in the functional topography of the locus coeruleus and their implications for cognitive and affective functions

The locus coeruleus (LC) is an important noradrenergic nucleus that has recently attracted a lot of attention because of its emerging role in cognitive and psychiatric disorders. Although previous histological studies have shown that the LC has heterogeneous connections and cellular features, no studies have yet assessed its functional topography in vivo, how this heterogeneity changes over aging, and whether it is associated with cognition and mood. Here, we employ a gradient-based approach to characterize the functional heterogeneity in the organization of the LC over aging using 3T resting-state fMRI in a population-based cohort aged from 18 to 88 years of age (Cambridge Centre for Ageing and Neuroscience cohort, n=618). We show that the LC exhibits a rostro-caudal functional gradient along its longitudinal axis, which was replicated in an independent dataset (Human Connectome Project [HCP] 7T dataset, n=184). Although the main rostro-caudal direction of this gradient was consistent across age groups, its spatial features varied with increasing age, emotional memory, and emotion regulation. More specifically, a loss of rostral-like connectivity, more clustered functional topography, and greater asymmetry between right and left LC gradients was associated with higher age and worse behavioral performance. Furthermore, participants with higher-than-normal Hospital Anxiety and Depression Scale (HADS) ratings exhibited alterations in the gradient as well, which manifested in greater asymmetry. These results provide an in vivo account of how the functional topography of the LC changes over aging, and imply that spatial features of this organization are relevant markers of LC-related behavioral measures and psychopathology.

Multiple routes of communication within the amygdala-mPFC network: A comparative approach in humans and macaques

The network formed by the amygdala (AMG) and the medial Prefrontal Cortex (mPFC), at the interface between our internal and external environment, has been shown to support some important aspects of behavioral adaptation. Whether and how the anatomo-functional organization of this network evolved across primates remains unclear. Here, we compared AMG nuclei morphological characteristics and their functional connectivity with the mPFC in humans and macaques to identify potential homologies and differences between these species. Based on selected studies, we highlight two subsystems within the AMG-mPFC circuits, likely involved in distinct temporal dynamics of integration during behavioral adaptation. We also show that whereas the mPFC displays a large expansion but a preserved intrinsic anatomo-functional organization, the AMG displays a volume reduction and morphological changes related to specific nuclei. We discuss potential commonalities and differences in the dialogue between AMG nuclei and mPFC in humans and macaques based on available data.

Age-related differences in the functional topography of the locus coeruleus: implications for cognitive and affective functions

The locus coeruleus (LC) is an important noradrenergic nucleus that has recently attracted a lot of attention because of its emerging role in cognitive and psychiatric disorders. Although previous histological studies have shown that the LC has heterogeneous connections and cellular features, no studies have yet assessed its functional topography in vivo, how this heterogeneity changes over aging and whether it is associated with cognition and mood. Here we employ a gradient-based approach to characterize the functional heterogeneity in the organization of the LC over aging using 3T resting-state fMRI in a population-based cohort aged from 18 to 88 years old (Cambridge Centre for Ageing and Neuroscience cohort, n=618). We show that the LC exhibits a rostro-caudal functional gradient along its longitudinal axis, which was replicated in an independent dataset (Human Connectome Project 7T dataset, n=184). Although the main rostro-caudal direction of this gradient was consistent across age groups, its spatial features varied with increasing age, emotional memory and emotion regulation. More specifically, a loss of rostral-like connectivity, more clustered functional topography and greater asymmetry between right and left LC gradients was associated with higher age and worse behavioral performance. Furthermore, participants with higher-than-normal Hospital Anxiety and Depression Scale ratings exhibited alterations in the gradient as well, which manifested in greater asymmetry. These results provide an in vivo account of how the functional topography of the LC changes over aging, and imply that spatial features of this organization are relevant markers of LC-related behavioral measures and psychopathology.

Improved Seizure Onset-Zone Lateralization in Temporal Lobe Epilepsy using 7T Resting-State fMRI: A Direct Comparison with 3T

Objective

Resting-state functional magnetic resonance imaging (rs-fMRI) at ultra high-field strengths (≥7T) is known to provide superior signal-to-noise and statistical power than comparable acquisitions at lower field strengths. In this study, we aim to provide a direct comparison of the seizure onset-zone (SOZ) lateralizing ability of 7T rs-fMRI and 3T rs-fMRI.

Methods

We investigated a cohort of 70 temporal lobe epilepsy (TLE) patients. A paired cohort of 19 patients had 3T and 7T rs-fMRI acquisitions for direct comparison between the two field strengths. Forty-three patients had only 3T, and 8 patients had only 7T rs-fMRI acquisitions. We quantified the functional connectivity between the hippocampus and other nodes within the default mode network (DMN) using seed-to-voxel connectivity, and measured how hippocampo-DMN connectivity could inform SOZ lateralization at 7T and 3T field strengths.

Results

Differences between hippocampo-DMN connectivity ipsilateral and contralateral to the SOZ were significantly higher at 7T (pFDR=0.008) than at 3T (pFDR=0.80) when measured in the same subjects. We found that our ability to lateralize the SOZ, by distinguishing subjects with left TLE from subjects with right TLE, was superior at 7T (AUC = 0.97) than 3T (AUC = 0.68). Our findings were reproduced in extended cohorts of subjects scanned at either 3T or 7T. Our rs-fMRI findings at 7T, but not 3T, are consistent and highly correlated (Spearman Rho=0.65) with clinical FDG-PET lateralizing hypometabolism.

Significance

We show superior SOZ lateralization in TLE patients when using 7T relative to 3T rs-fMRI, supporting the adoption of high-field strength functional imaging in the epilepsy presurgical evaluation.

White matter BOLD signals at 7 Tesla reveal visual field maps in optic radiation and vertical occipital fasciculus

There is growing evidence that blood-oxygen-level-dependent (BOLD) activity in the white matter (WM) can be detected by functional magnetic resonance imaging (fMRI). However, the functional relevance and significance of WM BOLD signals remain controversial. Here we investigated whether 7T BOLD fMRI can reveal fine-scale functional organizations of a WM bundle. Population receptive field (pRF) analyses of the 7T retinotopy dataset from the Human Connectome Project revealed clear contralateral retinotopic organizations of two visual WM bundles: the optic radiation (OR) and the vertical occipital fasciculus (VOF). The retinotopic maps of OR are highly consistent with post-mortem dissections and diffusion tractographies, while the VOF maps are compatible with the dorsal and ventral visual areas connected by the WM. Similar to the grey matter (GM) visual areas, both WM bundles show over-representations of the central visual field and increasing pRF size with eccentricity. Hemodynamic response functions of visual WM were slower and wider compared with those of GM areas. These findings clearly demonstrate that WM BOLD at 7 Tesla is closely coupled with neural activity related to axons, encoding highly specific information that can be used to characterize fine-scale functional organizations of a WM bundle.

Phenotyping Superagers Using Resting-State fMRI

BACKGROUND AND PURPOSE

Superagers are defined as older adults with episodic memory performance similar or superior to that in middle-aged adults. This study aimed to investigate the key differences in discriminative networks and their main nodes between superagers and cognitively average elderly controls. In addition, we sought to explore differences in sensitivity in detecting these functional activities across the networks at 3T and 7T MR imaging fields.

MATERIALS AND METHODS

Fifty-five subjects 80 years of age or older were screened using a detailed neuropsychological protocol, and 31 participants, comprising 14 superagers and 17 cognitively average elderly controls, were included for analysis. Participants underwent resting-state-fMRI at 3T and 7T MR imaging. A prediction classification algorithm using a penalized regression model on the measurements of the network was used to calculate the probabilities of a healthy older adult being a superager. Additionally, ORs quantified the influence of each node across preselected networks.

RESULTS

The key networks that differentiated superagers and elderly controls were the default mode, salience, and language networks. The most discriminative nodes (ORs > 1) in superagers encompassed areas in the precuneus posterior cingulate cortex, prefrontal cortex, temporoparietal junction, temporal pole, extrastriate superior cortex, and insula. The prediction classification model for being a superager showed better performance using the 7T compared with 3T resting-state-fMRI data set.

CONCLUSIONS

Our findings suggest that the functional connectivity in the default mode, salience, and language networks can provide potential imaging biomarkers for predicting superagers. The 7T field holds promise for the most appropriate study setting to accurately detect the functional connectivity patterns in superagers.

The Psychosis Human Connectome Project: Design and rationale for studies of visual neurophysiology

Visual perception is abnormal in psychotic disorders such as schizophrenia. In addition to hallucinations, laboratory tests show differences in fundamental visual processes including contrast sensitivity, center-surround interactions, and perceptual organization. A number of hypotheses have been proposed to explain visual dysfunction in psychotic disorders, including an imbalance between excitation and inhibition. However, the precise neural basis of abnormal visual perception in people with psychotic psychopathology (PwPP) remains unknown. Here, we describe the behavioral and 7 tesla MRI methods we used to interrogate visual neurophysiology in PwPP as part of the Psychosis Human Connectome Project (HCP). In addition to PwPP (n = 66) and healthy controls (n = 43), we also recruited first-degree biological relatives (n = 44) in order to examine the role of genetic liability for psychosis in visual perception. Our visual tasks were designed to assess fundamental visual processes in PwPP, whereas MR spectroscopy enabled us to examine neurochemistry, including excitatory and inhibitory markers. We show that it is feasible to collect high-quality data across multiple psychophysical, functional MRI, and MR spectroscopy experiments with a sizable number of participants at a single research site. These data, in addition to those from our previously described 3 tesla experiments, will be made publicly available in order to facilitate further investigations by other research groups. By combining visual neuroscience techniques and HCP brain imaging methods, our experiments offer new opportunities to investigate the neural basis of abnormal visual perception in PwPP.

Impact of autocalibration method on accelerated EPI of the cervical spinal cord at 7 T

PURPOSE

The spinal cord contains sensorimotor neural circuits of scientific and clinical interest. However, spinal cord functional MRI (fMRI) is significantly more technically demanding than brain fMRI, due primarily to its proximity to the lungs. Accelerated echo-planar imaging (EPI) at 7 T is particularly vulnerable to k-space phase inconsistencies induced by motion or B0 fluctuation, during either autocalibration signal (ACS) or time-series acquisition. For 7 T brain fMRI, sensitivity to motion and B0 fluctuation can be reduced using a re-ordered segmented EPI ACS based on the fast low-angle excitation echo-planar technique (FLEET). However, respiration-induced B0 fluctuations (exceeding 100 Hz at C7) are greater, and fewer k-space lines per slice are required for cervical spinal cord fMRI at 7 T, necessitating a separate evaluation of ACS methods.

METHODS

We compared 24-line single-shot EPI with 48-line two-shot segmented EPI, two-shot FLEET, and gradient echo (GRE)-based ACS acquisition methods, performed under various physiological conditions, in terms of temporal signal-to-noise ratio and prevalence of artifacts in generalized autocalibrating partially parallel acquisition (GRAPPA)-accelerated EPI of the cervical spinal cord at 7 T.

RESULTS

Segmented EPI and FLEET ACS produce images with nearly identical patterns of severe image artifacts. GRE and single-shot EPI ACS consistently produce images free from significant artifacts, and temporal signal-to-noise ratio is significantly greater for GRE ACS, particularly in lower slices where through-slice dephasing is most severe.

CONCLUSIONS

GRE and single-shot EPI-ACS acquisition methods, which are robust to respiration-induced phase errors between k-space segments, produce images with fewer and less severe artifacts than either FLEET or conventionally segmented EPI for accelerated EPI of the cervical spinal cord at 7 T.

Variability of the Surface Area of the V1, V2, and V3 Maps in a Large Sample of Human Observers

How variable is the functionally defined structure of early visual areas in human cortex and how much variability is shared between twins? Here we quantify individual differences in the best understood functionally defined regions of cortex: V1, V2, V3. The Human Connectome Project 7T Retinotopy Dataset includes retinotopic measurements from 181 subjects (109 female, 72 male), including many twins. We trained four "anatomists" to manually define V1-V3 using retinotopic features. These definitions were more accurate than automated anatomical templates and showed that surface areas for these maps varied more than threefold across individuals. This threefold variation was little changed when normalizing visual area size by the surface area of the entire cerebral cortex. In addition to varying in size, we find that visual areas vary in how they sample the visual field. Specifically, the cortical magnification function differed substantially among individuals, with the relative amount of cortex devoted to central vision varying by more than a factor of 2. To complement the variability analysis, we examined the similarity of visual area size and structure across twins. Whereas the twin sample sizes are too small to make precise heritability estimates (50 monozygotic pairs, 34 dizygotic pairs), they nonetheless reveal high correlations, consistent with strong effects of the combination of shared genes and environment on visual area size. Collectively, these results provide the most comprehensive account of individual variability in visual area structure to date, and provide a robust population benchmark against which new individuals and developmental and clinical populations can be compared.SIGNIFICANCE STATEMENT Areas V1, V2, and V3 are among the best studied functionally defined regions in human cortex. Using the largest retinotopy dataset to date, we characterized the variability of these regions across individuals and the similarity between twin pairs. We find that the size of visual areas varies dramatically (up to 3.5×) across healthy young adults, far more than the variability of the cerebral cortex size as a whole. Much of this variability appears to arise from inherited factors, as we find very high correlations in visual area size between monozygotic twin pairs, and lower but still substantial correlations between dizygotic twin pairs. These results provide the most comprehensive assessment of how functionally defined visual cortex varies across the population to date.

Denoise Functional Magnetic Resonance Imaging With Random Matrix Theory Based Principal Component Analysis

High-resolution functional MRI (fMRI) is largely hindered by random thermal noise. Random matrix theory (RMT)-based principal component analysis (PCA) is promising to reduce such noise in fMRI data. However, there is no consensus about the optimal strategy and practice in implementation. In this work, we propose a comprehensive RMT-based denoising method that consists of 1) rank and noise estimation based on a set of newly derived multiple criteria, and 2) optimal singular value shrinkage, with each module explained and implemented based on the RMT. By incorporating the variance stabilizing approach, the denoising method can deal with low signal-to-noise ratio (SNR) (such as <5) magnitude fMRI data with favorable performance compared to other state-of-the-art methods. Results from both simulation and in-vivo high-resolution fMRI data show that the proposed denoising method dramatically improves image restoration quality, promoting functional sensitivity at the same level of functional mapping blurring compared to existing denoising methods. Moreover, the denoising method can serve as a drop-in step in data preprocessing pipelines along with other procedures aimed at removal of structured physiological noises. We expect that the proposed denoising method will play an important role in leveraging high-quality, high-resolution task fMRI, which is desirable in many neuroscience and clinical applications.

Structural Covariance and Heritability of the Optic Tract and Primary Visual Cortex in Living Human Brains

Individual differences among human brains exist at many scales, spanning gene expression, white matter tissue properties, and the size and shape of cortical areas. One notable example is an approximately 3-fold range in the size of human primary visual cortex (V1), a much larger range than is found in overall brain size. A previous study (Andrews et al., 1997) reported a correlation between optic tract (OT) cross-section area and V1 size in postmortem human brains, suggesting that there may be a common developmental mechanism for multiple components of the visual pathways. We evaluated the relationship between properties of the OT and V1 in a much larger sample of living human brains by analyzing the Human Connectome Project (HCP) 7 Tesla Retinotopy Dataset (including 107 females and 71 males). This dataset includes retinotopic maps measured with functional MRI (fMRI) and fiber tract data measured with diffusion MRI (dMRI). We found a negative correlation between OT fractional anisotropy (FA) and V1 surface area (r = -0.19). This correlation, although small, was consistent across multiple dMRI datasets differing in acquisition parameters. Further, we found that both V1 surface area and OT properties were correlated among twins, with higher correlations for monozygotic (MZ) than dizygotic (DZ) twins, indicating a high degree of heritability for both properties. Together, these results demonstrate covariation across individuals in properties of the retina (OT) and cortex (V1) and show that each is influenced by genetic factors.SIGNIFICANCE STATEMENT The size of human primary visual cortex (V1) has large interindividual differences. These differences do not scale with overall brain size. A previous postmortem study reported a correlation between the size of the human optic tract (OT) and V1. In this study, we evaluated the relationship between the OT and V1 in living humans by analyzing a neuroimaging dataset that included functional MRI (fMRI) and diffusion MRI (dMRI) data. We found a small, but robust correlation between OT tissue properties and V1 size, supporting the existence of structural covariance between the OT and V1 in living humans. The results suggest that characteristics of retinal ganglion cells (RGCs), reflected in OT measurements, are correlated with individual differences in human V1.

Spatiotemporal patterns of spontaneous brain activity: a mini-review

The brain exists in a state of constant activity in the absence of any external sensory input. The spatiotemporal patterns of this spontaneous brain activity have been studied using various recording and imaging techniques. This has enabled considerable progress to be made in elucidating the cellular and network mechanisms that are involved in the observed spatiotemporal dynamics. This mini-review outlines different spatiotemporal dynamic patterns that have been identified in four commonly used modalities: electrophysiological recordings, optical imaging, functional magnetic resonance imaging, and electroencephalography. Signal sources for each modality, possible sources of the observed dynamics, and future directions are also discussed.

Age-Related Differences in Amygdala Activation Associated With Face Trustworthiness but No Evidence of Oxytocin Modulation

The amygdala has been shown to be responsive to face trustworthiness. While older adults typically give higher face trustworthiness ratings than young adults, a direct link between amygdala response and age-related differences in face trustworthiness evaluation has not yet been confirmed. Additionally, there is a possible modulatory role of the neuropeptide oxytocin in face trustworthiness evaluation, but the results are mixed and effects unexplored in aging. To address these research gaps, young, and older adults were randomly assigned to oxytocin or placebo self-administration via a nasal spray before rating faces on trustworthiness while undergoing functional magnetic resonance imaging. There was no overall age-group difference in face trustworthiness ratings, but older compared to young participants gave higher trustworthiness ratings to ambivalently untrustworthy-looking faces. In both age groups, lower face trustworthiness ratings were associated with higher left amygdala activity. A comparable negative linear association was observed in right amygdala but only among young participants. Also, in the right amygdala, lower and higher, compared to moderate, face trustworthiness ratings were associated with greater right amygdala activity (i.e., positive quadratic (U-shaped) association) for both age groups. Neither the behavioral nor the brain effects were modulated by a single dose of intranasal oxytocin administration, however. These results suggest dampened response to faces with lower trustworthiness among older compared to young adults, supporting the notion of reduced sensitivity to cues of untrustworthiness in aging. The findings also extend evidence of an age-related positivity effect to the evaluation of face trustworthiness.

Altered hippocampus and amygdala subregion connectome hierarchy in major depressive disorder

The hippocampus and amygdala limbic structures are critical to the etiology of major depressive disorder (MDD). However, there are no high-resolution characterizations of the role of their subregions in the whole brain network (connectome). Connectomic examination of these subregions can uncover disorder-related patterns that are otherwise missed when treated as single structures. 38 MDD patients and 40 healthy controls (HC) underwent anatomical and diffusion imaging using 7-Tesla MRI. Whole-brain segmentation was performed along with hippocampus and amygdala subregion segmentation, each representing a node in the connectome. Graph theory analysis was applied to examine the importance of the limbic subregions within the brain network using centrality features measured by node strength (sum of weights of the node's connections), Betweenness (number of shortest paths that traverse the node), and clustering coefficient (how connected the node's neighbors are to one another and forming a cluster). Compared to HC, MDD patients showed decreased node strength of the right hippocampus cornu ammonis (CA) 3/4, indicating decreased connectivity to the rest of the brain, and decreased clustering coefficient of the right dentate gyrus, implying it is less embedded in a cluster. Additionally, within the MDD group, the greater the embedding of the right amygdala central nucleus (CeA) in a cluster, the greater the severity of depressive symptoms. The altered role of these limbic subregions in the whole-brain connectome is related to diagnosis and depression severity, contributing to our understanding of the limbic system involvement in MDD and may elucidate the underlying mechanisms of depression.

Mitigating transmit‐B 1 artifacts by predicting parallel transmission images with deep learning: A feasibility study using high‐resolution whole‐brain diffusion at 7 Tesla

Purpose

To propose a novel deep learning (DL) approach to transmit‐B₁ (B₁⁺)‐artifact mitigation without direct use of parallel transmission (pTx), by predicting pTx images from single‐channel transmission (sTx) images.

Methods

A deep encoder–decoder convolutional neural network was constructed and trained to learn the mapping from sTx to pTx images. The feasibility was demonstrated using 7 T Human‐Connectome Project (HCP)‐style diffusion MRI. The training dataset comprised images acquired on 5 healthy subjects using commercial Nova RF coils. Relevant hyperparameters were tuned with a nested cross‐validation, and the generalization performance evaluated using a regular cross‐validation.

Results

Our DL method effectively improved the image quality for sTx images by restoring the signal dropout, with quality measures (including normalized root‐mean‐square error, peak SNR, and structural similarity index measure) improved in most brain regions. The improved image quality was translated into improved performances for diffusion tensor imaging analysis; our method improved accuracy for fractional anisotropy and mean diffusivity estimations, reduced the angular errors of principal eigenvectors, and improved the fiber orientation delineation relative to sTx images. Moreover, the final DL model trained on data of all 5 subjects was successfully used to predict pTx images for unseen new subjects (randomly selected from the 7 T HCP database), effectively recovering the signal dropout and improving color‐coded fractional anisotropy maps with largely reduced noise levels.

Conclusion

The proposed DL method has potential to provide images with reduced B1⁺ artifacts in healthy subjects even when pTx resources are inaccessible on the user side.

Discovering hidden brain network responses to naturalistic stimuli via tensor component analysis of multi-subject fMRI data

The study of brain network interactions during naturalistic stimuli facilitates a deeper understanding of human brain function. To estimate large-scale brain networks evoked with naturalistic stimuli, a tensor component analysis (TCA) based framework was used to characterize shared spatio-temporal patterns across subjects in a purely data-driven manner. In this framework, a third-order tensor is constructed from the timeseries extracted from all brain regions from a given parcellation, for all participants, with modes of the tensor corresponding to spatial distribution, time series and participants. TCA then reveals spatially and temporally shared components, i.e., evoked networks with the naturalistic stimuli, their time courses of activity and subject loadings of each component. To enhance the reproducibility of the estimation with the adaptive TCA algorithm, a novel spectral clustering method, tensor spectral clustering, was proposed and applied to evaluate the stability of the TCA algorithm. We demonstrated the effectiveness of the proposed framework via simulations and real fMRI data collected during a motor task with a traditional fMRI study design. We also applied the proposed framework to fMRI data collected during passive movie watching to illustrate how reproducible brain networks are evoked by naturalistic movie viewing.

Mapping of whole-cerebrum resting-state networks using ultra-high resolution acquisition protocols

Resting‐state functional magnetic resonance imaging (fMRI) has been used in numerous studies to map networks in the brain that employ spatially disparate regions. However, attempts to map networks with high spatial resolution have been hampered by conflicting technical demands and associated problems. Results from recent fMRI studies have shown that spatial resolution remains around 0.7 × 0.7 × 0.7 mm³, with only partial brain coverage. Therefore, this work aims to present a novel fMRI technique that was developed based on echo‐planar‐imaging with keyhole (EPIK) combined with repetition‐time‐external (TR‐external) EPI phase correction. Each technique has been previously shown to be effective in enhancing the spatial resolution of fMRI, and in this work, the combination of the two techniques into TR‐external EPIK provided a nominal spatial resolution of 0.51 × 0.51 × 1.00 mm³ (0.26 mm³ voxel) with whole‐cerebrum coverage. Here, the feasibility of using half‐millimetre in‐plane TR‐external EPIK for resting‐state fMRI was validated using 13 healthy subjects and the corresponding reproducible mapping of resting‐state networks was demonstrated. Furthermore, TR‐external EPIK enabled the identification of various resting‐state networks distributed throughout the brain from a single fMRI session, with mapping fidelity onto the grey matter at 7T. The high‐resolution functional image further revealed mesoscale anatomical structures, such as small cerebral vessels and the internal granular layer of the cortex within the postcentral gyrus.

From temporal to spatial topography: hierarchy of neural dynamics in higher- and lower-order networks shapes their complexity

The brain shows a topographical hierarchy along the lines of lower- and higher-order networks. The exact temporal dynamics characterization of this lower-higher-order topography at rest and its impact on task states remains unclear, though. Using 2 functional magnetic resonance imaging data sets, we investigate lower- and higher-order networks in terms of the signal compressibility, operationalized by Lempel-Ziv complexity (LZC). As we assume that this degree of complexity is related to the slow-fast frequency balance, we also compute the median frequency (MF), an estimation of frequency distribution. We demonstrate (i) topographical differences at rest between higher- and lower-order networks, showing lower LZC and MF in the former; (ii) task-related and task-specific changes in LZC and MF in both lower- and higher-order networks; (iii) hierarchical relationship between LZC and MF, as MF at rest correlates with LZC rest-task change along the lines of lower- and higher-order networks; and (iv) causal and nonlinear relation between LZC at rest and LZC during task, with MF at rest acting as mediator. Together, results show that the topographical hierarchy of lower- and higher-order networks converges with their temporal hierarchy, with these neural dynamics at rest shaping their range of complexity during task states in a nonlinear way.

Medial temporal lobe contributions to resting-state networks

The medial temporal lobe (MTL) is a set of interconnected brain regions that have been shown to play a central role in behavior as well as in neurological disease. Recent studies using resting-state functional magnetic resonance imaging (rsfMRI) have attempted to understand the MTL in terms of its functional connectivity with the rest of the brain. However, the exact characterization of the whole-brain networks that co-activate with the MTL as well as how the various sub-regions of the MTL are associated with these networks remains poorly understood. Here, we attempted to advance these issues by exploiting the high spatial resolution 7T rsfMRI dataset from the Human Connectome Project with a data-driven analysis approach that relied on independent component analysis (ICA) restricted to the MTL. We found that four different well-known resting-state networks co-activated with a unique configuration of MTL subcomponents. Specifically, we found that different sections of the parahippocampal cortex were involved in the default mode, visual and dorsal attention networks; sections of the hippocampus in the somatomotor and default mode networks; and the lateral entorhinal cortex in the dorsal attention network. We replicated this set of results in a validation sample. These results provide new insight into how the MTL and its subcomponents contribute to known resting-state networks. The participation of the MTL in an expanded range of resting-state networks is in line with recent proposals on MTL function.

Ultra-high field neuroimaging in psychosis: A narrative review

Schizophrenia and related psychoses are complex neuropsychiatric diseases representing dysconnectivity across multiple scales, through the micro (cellular), meso (brain network), manifest (behavioral), and social (interpersonal) levels. In vivo human neuroimaging, particularly at ultra-high field (UHF), offers unprecedented opportunity to examine multiscale dysconnectivity in psychosis. In this review, we provide an overview of the literature to date on UHF in psychosis, focusing on microscale findings from magnetic resonance spectroscopy (MRS), mesoscale studies on structural and functional magnetic resonance imaging (fMRI), and multiscale studies assessing multiple neuroimaging modalities and relating UHF findings to behavior. We highlight key insights and considerations from multiscale and longitudinal studies and provide recommendations for future research on UHF neuroimaging in psychosis.

Test–Retest Reliability of Synchrony and Metastability in Resting State fMRI

In recent years, interest has been growing in dynamic characteristic of brain signals from resting-state functional magnetic resonance imaging (rs-fMRI). Synchrony and metastability, as neurodynamic indexes, are considered as one of methods for analyzing dynamic characteristics. Although much research has studied the analysis of neurodynamic indices, few have investigated its reliability. In this paper, the datasets from the Human Connectome Project have been used to explore the test–retest reliabilities of synchrony and metastability from multiple angles through intra-class correlation (ICC). The results showed that both of these indexes had fair test–retest reliability, but they are strongly affected by the field strength, the spatial resolution, and scanning interval, less affected by the temporal resolution. Denoising processing can help improve their ICC values. In addition, the reliability of neurodynamic indexes was affected by the node definition strategy, but these effects were not apparent. In particular, by comparing the test–retest reliability of different resting-state networks, we found that synchrony of different networks was basically stable, but the metastability varied considerably. Among these, DMN and LIM had a relatively higher test–retest reliability of metastability than other networks. This paper provides a methodological reference for exploring the brain dynamic neural activity by using synchrony and metastability in fMRI signals.

fMRI-Informed EEG for brain mapping of imagined lower limb movement: Feasibility of a brain computer interface

BACKGROUND

EEG and fMRI have contributed greatly to our understanding of brain activity and its link to behaviors by helping to identify both when and where the activity occurs. This is particularly important in the development of brain-computer interfaces (BCIs), where feed forward systems gather data from imagined brain activity and then send that information to an effector. The purpose of this study was to develop and evaluate a computational approach that enables an accurate mapping of spatial brain activity (fMRI) in relation to the temporal receptors (EEG electrodes) associated with imagined lower limb movement.

NEW METHOD

EEG and fMRI data from 16 healthy, male participants while imagining lower limb movement were used for this purpose. A combined analysis of fMRI data and EEG electrode locations was developed to identify EEG electrodes with a high likelihood of capturing imagined lower limb movement originating from various clusters of brain activity. This novel feature selection tool was used to develop an artificial neural network model to classify right and left lower limb movement.

RESULTS

Results showed that left versus right lower limb imagined movement could be classified with 66.5% accuracy using this approach. Comparison with existing methods: Adopting a purely data-driven approach for feature selection to use in the right/left classification task resulted in the same accuracy (66.6%) but with reduced interpretability.

CONCLUSIONS

The developed fMRI-informed EEG approach could pave the way towards improved brain computer interfaces for lower limb movement while also being applicable to other systems where fMRI could be helpful to inform EEG acquisition and processing.

Cognitive control affects motor learning through local variations in GABA within the primary motor cortex

The primary motor cortex (M1) is crucial for motor learning; however, its interaction with other brain areas during motor learning remains unclear. We hypothesized that the fronto-parietal execution network (FPN) provides learning-related information critical for the flexible cognitive control that is required for practice. We assessed network-level changes during sequential finger tapping learning under speed pressure by combining magnetic resonance spectroscopy and task and resting-state functional magnetic resonance imaging. There was a motor learning-related increase in preparatory activity in the fronto-parietal regions, including the right M1, overlapping the FPN and sensorimotor network (SMN). Learning-related increases in M1-seeded functional connectivity with the FPN, but not the SMN, were associated with decreased GABA/glutamate ratio in the M1, which were more prominent in the parietal than the frontal region. A decrease in the GABA/glutamate ratio in the right M1 was positively correlated with improvements in task performance (p = 0.042). Our findings indicate that motor learning driven by cognitive control is associated with local variations in the GABA/glutamate ratio in the M1 that reflects remote connectivity with the FPN, representing network-level motor sequence learning formations.

Neuroimaging at 3T vs 7T: Is It Really Worth It?

Food and Drug Administration approval of 7T MR imaging allows ultrahigh-field neuroimaging to extend from the research realm into the clinical realm. Increased signal is clinically advantageous for smaller voxels and thereby high spatial resolution imaging, with additional advantages of increased tissue contrast. Susceptibility, time-of-flight signal, and blood oxygen level-dependent signal also have favorable clinical benefit from 7T. This article provides a survey of clinical cases showcasing some advantages of 7T.

Statistical power or more precise insights into neuro-temporal dynamics? Assessing the benefits of rapid temporal sampling in fMRI

Functional magnetic resonance imaging (fMRI), a non-invasive and widely used human neuroimaging method, is most known for its spatial precision. However, there is a growing interest in its temporal sensitivity. This is despite the temporal blurring of neuronal events by the blood oxygen level dependent (BOLD) signal, the peak of which lags neuronal firing by 4-6 seconds. Given this, the goal of this review is to answer a seemingly simple question - "What are the benefits of increased temporal sampling for fMRI?". To answer this, we have combined fMRI data collected at multiple temporal scales, from 323 to 1000 milliseconds, with a review of both historical and contemporary temporal literature. After a brief discussion of technological developments that have rekindled interest in temporal research, we next consider the potential statistical and methodological benefits. Most importantly, we explore how fast fMRI can uncover previously unobserved neuro-temporal dynamics - effects that are entirely missed when sampling at conventional 1 to 2 second rates. With the intrinsic link between space and time in fMRI, this temporal renaissance also delivers improvements in spatial precision. Far from producing only statistical gains, the array of benefits suggest that the continued temporal work is worth the effort.

Learning brain dynamics for decoding and predicting individual differences

Insights from functional Magnetic Resonance Imaging (fMRI), as well as recordings of large numbers of neurons, reveal that many cognitive, emotional, and motor functions depend on the multivariate interactions of brain signals. To decode brain dynamics, we propose an architecture based on recurrent neural networks to uncover distributed spatiotemporal signatures. We demonstrate the potential of the approach using human fMRI data during movie-watching data and a continuous experimental paradigm. The model was able to learn spatiotemporal patterns that supported 15-way movie-clip classification (∼90%) at the level of brain regions, and binary classification of experimental conditions (∼60%) at the level of voxels. The model was also able to learn individual differences in measures of fluid intelligence and verbal IQ at levels comparable to that of existing techniques. We propose a dimensionality reduction approach that uncovers low-dimensional trajectories and captures essential informational (i.e., classification related) properties of brain dynamics. Finally, saliency maps and lesion analysis were employed to characterize brain-region/voxel importance, and uncovered how dynamic but consistent changes in fMRI activation influenced decoding performance. When applied at the level of voxels, our framework implements a dynamic version of multivariate pattern analysis. Our approach provides a framework for visualizing, analyzing, and discovering dynamic spatially distributed brain representations during naturalistic conditions.

Minimal specifications for non-human primate MRI: Challenges in standardizing and harmonizing data collection

Recent methodological advances in MRI have enabled substantial growth in neuroimaging studies of non-human primates (NHPs), while open data-sharing through the PRIME-DE initiative has increased the availability of NHP MRI data and the need for robust multi-subject multi-center analyses. Streamlined acquisition and analysis protocols would accelerate and improve these efforts. However, consensus on minimal standards for data acquisition protocols and analysis pipelines for NHP imaging remains to be established, particularly for multi-center studies. Here, we draw parallels between NHP and human neuroimaging and provide minimal guidelines for harmonizing and standardizing data acquisition. We advocate robust translation of widely used open-access toolkits that are well established for analyzing human data. We also encourage the use of validated, automated pre-processing tools for analyzing NHP data sets. These guidelines aim to refine methodological and analytical strategies for small and large-scale NHP neuroimaging data. This will improve reproducibility of results, and accelerate the convergence between NHP and human neuroimaging strategies which will ultimately benefit fundamental and translational brain science.

Reducing SAR in 7T brain fMRI by circumventing fat suppression while removing the lipid signal through a parallel acquisition approach

Ultra-high-field functional magnetic resonance imaging (fMRI) offers a way to new insights while increasing the spatial and temporal resolution. However, a crucial concern in 7T human MRI is the increase in power deposition, supervised through the specific absorption rate (SAR). The SAR limitation can restrict the brain coverage or the minimal repetition time of fMRI experiments. In the majority of today’s studies fMRI relies on the well-known gradient-echo echo-planar imaging (GRE-EPI) sequence, which offers ultrafast acquisition. Commonly, the GRE-EPI sequence comprises two pulses: fat suppression and excitation. This work provides the means for a significant reduction in the SAR by circumventing the fat-suppression pulse. Without this fat-suppression, however, lipid signal can result in artifacts due to the chemical shift between the lipid and water signals. Our approach exploits a reconstruction similar to the simultaneous-multi-slice method to separate the lipid and water images, thus avoiding undesired lipid artifacts in brain images. The lipid-water separation is based on the known spatial shift of the lipid signal, which can be detected by the multi-channel coils sensitivity profiles. Our study shows robust human imaging, offering greater flexibility to reduce the SAR, shorten the repetition time or increase the volume coverage with substantial benefit for brain functional studies.

Structural and resting state functional connectivity beyond the cortex

Mapping the structural and functional connectivity of the central nervous system has become a key area within neuroimaging research. While detailed network structures across the entire brain have been probed using animal models, non-invasive neuroimaging in humans has thus far been dominated by cortical investigations. Beyond the cortex, subcortical nuclei have traditionally been less accessible due to their smaller size and greater distance from radio frequency coils. However, major neuroimaging developments now provide improved signal and the resolution required to study these structures. Here, we present an overview of the connectivity between the amygdala, brainstem, cerebellum, spinal cord and the rest of the brain. While limitations to their imaging and analyses remain, we also provide some recommendations and considerations for mapping brain connectivity beyond the cortex.

Assessing methods for geometric distortion compensation in 7 T gradient echo functional MRI data

Echo planar imaging (EPI) is widely used in functional and diffusion‐weighted MRI, but suffers from significant geometric distortions in the phase encoding direction caused by inhomogeneities in the static magnetic field (B₀). This is a particular challenge for EPI at very high field (≥7 T), as distortion increases with higher field strength. A number of techniques for distortion correction exist, including those based on B₀ field mapping and acquiring EPI scans with opposite phase encoding directions. However, few quantitative comparisons of distortion compensation methods have been performed using human EPI data, especially at very high field. Here, we compared distortion compensation using B₀ field maps and opposite phase encoding scans in two different software packages (FSL and AFNI) applied to 7 T gradient echo (GE) EPI data from 31 human participants. We assessed distortion compensation quality by quantifying alignment to anatomical reference scans using Dice coefficients and mutual information. Performance between FSL and AFNI was equivalent. In our whole‐brain analyses, we found superior distortion compensation using GE scans with opposite phase encoding directions, versus B₀ field maps or spin echo (SE) opposite phase encoding scans. However, SE performed better when analyses were limited to ventromedial prefrontal cortex, a region with substantial dropout. Matching the type of opposite phase encoding scans to the EPI data being corrected (e.g., SE‐to‐SE) also yielded better distortion correction. While the ideal distortion compensation approach likely varies depending on methodological differences across experiments, this study provides a framework for quantitative comparison of different distortion compensation methods.

ULTRAHIGH FIELD and ULTRAHIGH RESOLUTION fMRI

Functional magnetic resonance imaging (fMRI) has become one of the most powerful tools for investigating the human brain. Ultrahigh magnetic field (UHF) of 7 Tesla has played a critical role in enabling higher resolution and more accurate (relative to the neuronal activity) functional maps. However, even with these gains, the fMRI approach is challenged relative to the spatial scale over which brain function is organized. Therefore, going forward, significant advances in fMRI are still needed. Such advances will predominantly come from magnetic fields significantly higher than 7 Tesla, which is the most commonly used UHF platform today, and additional technologies that will include developments in pulse sequences, image reconstruction, noise suppression, and image analysis in order to further enhance and augment the gains than can be realized by going to higher magnetic fields.

Cortical response to naturalistic stimuli is largely predictable with deep neural networks

Naturalistic stimuli, such as movies, activate a substantial portion of the human brain, invoking a response shared across individuals. Encoding models that predict neural responses to arbitrary stimuli can be very useful for studying brain function. However, existing models focus on limited aspects of naturalistic stimuli, ignoring the dynamic interactions of modalities in this inherently context-rich paradigm. Using movie-watching data from the Human Connectome Project, we build group-level models of neural activity that incorporate several inductive biases about neural information processing, including hierarchical processing, temporal assimilation, and auditory-visual interactions. We demonstrate how incorporating these biases leads to remarkable prediction performance across large areas of the cortex, beyond the sensory-specific cortices into multisensory sites and frontal cortex. Furthermore, we illustrate that encoding models learn high-level concepts that generalize to task-bound paradigms. Together, our findings underscore the potential of encoding models as powerful tools for studying brain function in ecologically valid conditions.

Identifying Biological Signatures of N-Acetylcysteine for Non-Suicidal Self-Injury in Adolescents and Young Adults

The prevalence of non-suicidal self-injury (NSSI) is high in adolescents and young adults. However, there is a paucity of evidence-based treatments to address this clinical problem. An open-label, pilot study in the target population showed that treatment with oral N-acetylcysteine (NAC), a widely available dietary supplement, was associated with reduction in NSSI frequency. In preparation for a biologically informed design of an efficacy trial, a critical preliminary step is to clarify NAC's biological signatures, or measures of the mechanisms underlying its clinical effects. Toward that end, we propose a 2-stage project to investigate NAC's biological signatures (changes in glutathione (GSH) and/or glutamate (Glu)) in women with NSSI. The first stage; a double-blind randomized placebo-controlled study will focus on identifying the optimal dose to achieve meaningful change in GSH and Glu during short-term (4 weeks) NAC treatment in 36 women aged 16-24 years with NSSI. Go/No-go criteria to determine if the study will progress to the second stage include pre-specified changes in brain and blood measures of GSH. Changes in the brain GSH are measured through magnetic resonance spectroscopy (MRS). The dose for the stage 2 will be selected based on the biological changes and the tolerability observed in the stage 1. The stage 2 will seek to replicate the biological signature findings in an 8-week trial in a new patient cohort, and examine the relationships among biological signatures, NAC pharmacokinetics and clinical response. This 2-stage project is unique as it unifies clinical psychiatric measurements, quantitative MRS and pharmacological approaches in the first placebo-controlled clinical trial of NAC in young women with NSSI.

TRIAL REGISTRATION

The stage 1 trial protocol has been registered on https://clinicaltrials.gov/ with ClinicalTrials.gov ID "NCT04005053" (Registered on 02 July 2019. Available from: https://clinicaltrials.gov/ct2/show/NCT04005053).

Clarifying the role of higher-level cortices in resolving perceptual ambiguity using ultra high field fMRI

The brain is organized into distinct, flexible networks. Within these networks, cognitive variables such as attention can modulate sensory representations in accordance with moment-to-moment behavioral requirements. These modulations can be studied by varying task demands; however, the tasks employed are often incongruent with the postulated functions of a sensory system, limiting the characterization of the system in relation to natural behaviors. Here we combine domain-specific task manipulations and ultra-high field fMRI to study the nature of top-down modulations. We exploited faces, a visual category underpinned by a complex cortical network, and instructed participants to perform either a stimulus-relevant/domain-specific or a stimulus-irrelevant task in the scanner. We found that 1. perceptual ambiguity (i.e. difficulty of achieving a stable percept) is encoded in top-down modulations from higher-level cortices; 2. the right inferior-temporal lobe is active under challenging conditions and uniquely encodes trial-by-trial variability in face perception.

Ultra-high field (10.5 T) resting state fMRI in the macaque

Resting state functional connectivity refers to the temporal correlations between spontaneous hemodynamic signals obtained using functional magnetic resonance imaging. This technique has demonstrated that the structure and dynamics of identifiable networks are altered in psychiatric and neurological disease states. Thus, resting state network organizations can be used as a diagnostic, or prognostic recovery indicator. However, much about the physiological basis of this technique is unknown. Thus, providing a translational bridge to an optimal animal model, the macaque, in which invasive circuit manipulations are possible, is of utmost importance. Current approaches to resting state measurements in macaques face unique challenges associated with signal-to-noise, the need for contrast agents limiting translatability, and within-subject designs. These limitations can, in principle, be overcome through ultra-high magnetic fields. However, imaging at magnetic fields above 7T has yet to be adapted for fMRI in macaques. Here, we demonstrate that the combination of high channel count transmitter and receiver arrays, optimized pulse sequences, and careful anesthesia regimens, allows for detailed single-subject resting state analysis at high resolutions using a 10.5 Tesla scanner. In this study, we uncover thirty spatially detailed resting state components that are highly robust across individual macaques and closely resemble the quality and findings of connectomes from large human datasets. This detailed map of the rsfMRI 'macaque connectome' will be the basis for future neurobiological circuit manipulation work, providing valuable biological insights into human connectomics.

Higher temporal resolution multiband fMRI provides improved presurgical language maps

PURPOSE

We investigated the hypothesis that increasing fMRI temporal resolution using a multiband (MB) gradient echo-echo planar imaging (GRE-EPI) pulse sequence provides fMRI language maps of higher statistical quality than those acquired with a traditional GRE-EPI sequence.

METHODS

This prospective study enrolled 29 consecutive patients receiving language fMRI prior to a potential brain resection for tumor, AVM, or epilepsy. A 4-min rhyming task was performed at 3.0 Tesla with a traditional GRE-EPI pulse sequence (TR = 2000, TE = 30, matrix = 64/100%, slice = 4/0, FOV = 24, slices = 30, time points = 120) and an additional MB GRE-EPI pulse sequence with an acceleration factor of 6 (TR = 333, TE = 30, matrix 64/100%, slice = 4/0, FOV = 24, time points = 720). Spatially filtered t statistical maps were generated. Volumes of interest (VOIs) were drawn around activations at Broca's, dorsolateral prefrontal cortex, Wernicke's, and the visual word form areas. The t value maxima were measured for the overall brain and each of the VOIs. A paired t test was performed for the corresponding traditional and MB GRE-EPI measurements.

RESULTS

The mean age of subjects was 42.6 years old (18-75). Sixty-two percent were male. The average overall brain t statistic maxima for the MB pulse sequence (t = 15.4) was higher than for the traditional pulse sequence (t = 9.3, p = < .0001). This also held true for Broca's area (p < 0.0001), Wernicke's area (p < .0001), dorsolateral prefrontal cortex (p < .0001), and the visual word form area (p < .0001).

CONCLUSION

A MB GRE-EPI fMRI pulse sequence employing high temporal resolution provides clinical fMRI language maps of greater statistical significance than those obtained with a traditional GRE-EPI sequence.

Topographic organization of the human subcortex unveiled with functional connectivity gradients

Brain atlases are fundamental to understanding the topographic organization of the human brain, yet many contemporary human atlases cover only the cerebral cortex, leaving the subcortex a terra incognita. We use functional MRI (fMRI) to map the complex topographic organization of the human subcortex, revealing large-scale connectivity gradients and new areal boundaries. We unveil four scales of subcortical organization that recapitulate well-known anatomical nuclei at the coarsest scale and delineate 27 new bilateral regions at the finest. Ultrahigh field strength fMRI corroborates and extends this organizational structure, enabling the delineation of finer subdivisions of the hippocampus and the amygdala, while task-evoked fMRI reveals a subtle subcortical reorganization in response to changing cognitive demands. A new subcortical atlas is delineated, personalized to represent individual differences and used to uncover reproducible brain-behavior relationships. Linking cortical networks to subcortical regions recapitulates a task-positive to task-negative axis. This new atlas enables holistic connectome mapping and characterization of cortico-subcortical connectivity.

New Cognitive Neurotechnology Facilitates Studies of Cortical-Subcortical Interactions

Most of the studies employing neuroimaging have focused on cortical and subcortical signals individually to obtain neurophysiological signatures of cognitive functions. However, understanding the dynamic communication between the cortex and subcortical structures is essential for unraveling the neural correlates of cognition. In this quest, magnetoencephalography (MEG) and electroencephalography (EEG) are the methods of choice because they are noninvasive electrophysiological recording techniques with high temporal resolution. Sophisticated MEG/EEG source estimation techniques and network analysis methods, developed recently, can provide a more comprehensive understanding of the neurophysiological mechanisms of fundamental cognitive processes. Used together with noninvasive modulation of cortical-subcortical communication, these approaches may open up new possibilities for expanding the repertoire of noninvasive cognitive neurotechnology.