Highfield imaging of the subgenual anterior cingulate cortex in uni- and bipolar depression.

Background: The subgenual Anterior Cingulate Cortex (sgACC), as a part of the Anterior Cingulate Cortex and the limbic system plays a crucial role in mood regulation. Previous structural and functional brain imaging studies of the sgACC have revealed alterations of Gray Matter (GM) volumes and Blood Oxygenation Level Dependent signals (BOLD) in patients with Major Depressive Disorder (MDD) and Bipolar Disorder (BD), suggesting potential biomarker traits for affective disorders. Method: In this study we investigated the gray matter volume of the sgACC in 3 different patient groups: 40 MDD patients, of which 20 were medicated (MDDm) and 20 were unmedicated (MDDu), and 21 medicated BD patients, and compared them with 23 healthy volunteers. We examined GM volume alteration using high-resolution 7T Magnetic Resonance Imaging (MRI) which produced quantitative maps of the spin-lattice relaxation time (T1). T1 maps provide high contrast between gray and white matter, and at 7 Tesla voxels with submillimeter resolution can be acquired in a reasonable scan time. We developed a semi-automatic segmentation protocol based on refined landmarks derived from previous volumetric studies using quantitative T1 maps as raw input data for automatic tissue segmentation of GM, WM and CSF (cerebrospinal fluid) tissue. The sgACC ROI was then superimposed on these tissue probability maps and traced manually by two independent raters (F.B., M.M.) following our semi-automatic segmentation protocol. Interrater reliability was calculated on a subset of 10 brain scans for each hemisphere, showing an Intra-Class Correlation coefficient (ICC) r = 0.96 for left sgACC and r = 0.84 for right sgACC respectively. In summary, we have developed a reproducible and reliable semi-automatic segmentation protocol to measure gray matter volume in the sgACC. Based on previous findings from meta-analyses on morphometric studies of the sgACC, we hypothesized that patients with MDD would have lower gray matter sgACC volumes compared to healthy subjects. Results: Post-hoc analysis revealed smaller subgenual volumes for the left hemisphere in both the medicated (MDDm) and non-medicated (MDDu) group versus healthy controls (p = .001, p = .008) respectively. For the right hemisphere, the (MDDu) and BD group exhibited significantly lower subgenual volumes than healthy controls (p < .001, p = .004) respectively. Conclusion: To our knowledge, this is the first morphometric MRI study using T1 maps obtained in high-resolution 7 Tesla MRI to compare MDD and BD patients with healthy controls.

A high-resolution 7 Tesla resting-state fMRI dataset optimized for studying the subcortex.

To achieve a comprehensive understanding of spontaneous brain dynamics in humans, in vivo acquisition of intrinsic activity across both cortical and subcortical regions is necessary. Here we present advanced whole-brain, resting-state functional magnetic resonance imaging (rs-fMRI) data acquired at 7 Tesla with 1.5 mm isotropic voxel resolution. Functional images were obtained from 56 healthy adults (33 females, ages 19-39 years) in two runs of 15 min eyes-open wakeful rest. The high spatial resolution and short echo times of the multiband echo-planar imaging (EPI) protocol optimizes blood oxygen level-dependent (BOLD)-sensitivity for the subcortex while concurrent respiratory and cardiac measures enable retrospective correction of physiological noise, resulting in data that is highly suitable for researchers interested in subcortical BOLD signal. Functional timeseries were coregistered to high-resolution T1-weighted structural data (0.75 mm isotropic voxels) acquired during the same scanning session. To accommodate data reutilization, functional and structural images were formatted to the Brain Imaging Data Structure (BIDS) and preprocessed with fMRIPrep.

Decreased long-range temporal correlations in the resting-state functional magnetic resonance imaging blood-oxygen-level-dependent signal reflect motor sequence learning up to 2 weeks following training.

Decreased long-range temporal correlations (LRTC) in brain signals can be used to measure cognitive effort during task execution. Here, we examined how learning a motor sequence affects long-range temporal memory within resting-state functional magnetic resonance imaging signal. Using the Hurst exponent (HE), we estimated voxel-wise LRTC and assessed changes over 5 consecutive days of training, followed by a retention scan 12 days later. The experimental group learned a complex visuomotor sequence while a complementary control group performed tightly matched movements. An interaction analysis revealed that HE decreases were specific to the complex sequence and occurred in well-known motor sequence learning associated regions including left supplementary motor area, left premotor cortex, left M1, left pars opercularis, bilateral thalamus, and right striatum. Five regions exhibited moderate to strong negative correlations with overall behavioral performance improvements. Following learning, HE values returned to pretraining levels in some regions, whereas in others, they remained decreased even 2 weeks after training. Our study presents new evidence of HE's possible relevance for functional plasticity during the resting-state and suggests that a cortical subset of sequence-specific regions may continue to represent a functional signature of learning reflected in decreased long-range temporal dependence after a period of inactivity.

Echoes from Intrinsic Connectivity Networks in the Subcortex.

Decades of research have greatly improved our understanding of intrinsic human brain organization in terms of functional networks and the transmodal hubs within the cortex at which they converge. However, substrates of multinetwork integration in the human subcortex are relatively uncharted. Here, we leveraged recent advances in subcortical atlasing and ultra-high field (7 T) imaging optimized for the subcortex to investigate the functional architecture of 14 individual structures in healthy adult males and females with a fully data-driven approach. We revealed that spontaneous neural activity in subcortical regions can be decomposed into multiple independent subsignals that correlate with, or "echo," the activity in functional networks across the cortex. Distinct subregions of the thalamus, striatum, claustrum, and hippocampus showed a varied pattern of echoes from attention, control, visual, somatomotor, and default mode networks, demonstrating evidence for a heterogeneous organization supportive of functional integration. Multiple network activity furthermore converged within the globus pallidus externa, substantia nigra, and ventral tegmental area but was specific to one subregion, while the amygdala and pedunculopontine nucleus preferentially affiliated with a single network, showing a more homogeneous topography. Subregional connectivity of the globus pallidus interna, subthalamic nucleus, red nucleus, periaqueductal gray, and locus coeruleus did not resemble patterns of cortical network activity. Together, these finding describe potential mechanisms through which the subcortex participates in integrated and segregated information processing and shapes the spontaneous cognitive dynamics during rest.SIGNIFICANCE STATEMENT Despite the impact of subcortical dysfunction on brain health and cognition, large-scale functional mapping of subcortical structures severely lags behind that of the cortex. Recent developments in subcortical atlasing and imaging at ultra-high field provide new avenues for studying the intricate functional architecture of the human subcortex. With a fully data-driven analysis, we reveal subregional connectivity profiles of a large set of noncortical structures, including those rarely studied in fMRI research. The results have implications for understanding how the functional organization of the subcortex facilitates integrative processing through cross-network information convergence, paving the way for future work aimed at improving our knowledge of subcortical contributions to intrinsic brain dynamics and spontaneous cognition.

Modeling venous bias in resting state functional MRI metrics.

Resting-state (rs) functional magnetic resonance imaging (fMRI) is used to detect low-frequency fluctuations in the blood oxygen-level dependent (BOLD) signal across brain regions. Correlations between temporal BOLD signal fluctuations are commonly used to infer functional connectivity. However, because BOLD is based on the dilution of deoxyhemoglobin, it is sensitive to veins of all sizes, and its amplitude is biased by draining veins. These biases affect local BOLD signal location and amplitude, and may also influence BOLD-derived connectivity measures, but the magnitude of this venous bias and its relation to vein size and proximity is unknown. Here, veins were identified using high-resolution quantitative susceptibility maps and utilized in a biophysical model to investigate systematic venous biases on common local rsfMRI-derived measures. Specifically, we studied the impact of vein diameter and distance to veins on the amplitude of low-frequency fluctuations (ALFF), fractional ALFF (fALFF), Hurst exponent (HE), regional homogeneity (ReHo), and eigenvector centrality values in the grey matter. Values were higher across all distances in smaller veins, and decreased with increasing vein diameter. Additionally, rsfMRI values associated with larger veins decrease with increasing distance from the veins. ALFF and ReHo were the most biased by veins, while HE and fALFF exhibited the smallest bias. Across all metrics, the amplitude of the bias was limited in voxel-wise data, confirming that venous structure is not the dominant source of contrast in these rsfMRI metrics. Finally, the models presented can be used to correct this venous bias in rsfMRI metrics.

Mixed methodology in human brain research: integrating MRI and histology.

Postmortem magnetic resonance imaging (MRI) can provide a bridge between histological observations and the in vivo anatomy of the human brain. Approaches aimed at the co-registration of data derived from the two techniques are gaining interest. Optimal integration of the two research fields requires detailed knowledge of the tissue property requirements for individual research techniques, as well as a detailed understanding of the consequences of tissue fixation steps on the imaging quality outcomes for both MRI and histology. Here, we provide an overview of existing studies that bridge between state-of-the-art imaging modalities, and discuss the background knowledge incorporated into the design, execution and interpretation of postmortem studies. A subset of the discussed challenges transfer to animal studies as well. This insight can contribute to furthering our understanding of the normal and diseased human brain, and to facilitate discussions between researchers from the individual disciplines.

A unified model for reconstruction and R2* mapping of accelerated 7T data using the quantitative recurrent inference machine.

Quantitative MRI (qMRI) acquired at the ultra-high field of 7 Tesla has been used in visualizing and analyzing subcortical structures. qMRI relies on the acquisition of multiple images with different scan settings, leading to extended scanning times. Data redundancy and prior information from the relaxometry model can be exploited by deep learning to accelerate the imaging process. We propose the quantitative Recurrent Inference Machine (qRIM), with a unified forward model for joint reconstruction and R2*-mapping from sparse data, embedded in a Recurrent Inference Machine (RIM), an iterative inverse problem-solving network. To study the dependency of the proposed extension of the unified forward model to network architecture, we implemented and compared a quantitative End-to-End Variational Network (qE2EVN). Experiments were performed with high-resolution multi-echo gradient echo data of the brain at 7T of a cohort study covering the entire adult life span. The error in reconstructed R2* from undersampled data relative to reference data significantly decreased for the unified model compared to sequential image reconstruction and parameter fitting using the RIM. With increasing acceleration factor, an increasing reduction in the reconstruction error was observed, pointing to a larger benefit for sparser data. Qualitatively, this was following an observed reduction of image blurriness in R2*-maps. In contrast, when using the U-Net as network architecture, a negative bias in R2* in selected regions of interest was observed. Compressed Sensing rendered accurate, but less precise estimates of R2*. The qE2EVN showed slightly inferior reconstruction quality compared to the qRIM but better quality than the U-Net and Compressed Sensing. Subcortical maturation over age measured by a linearly increasing interquartile range of R2* in the striatum was preserved up to an acceleration factor of 9. With the integrated prior of the unified forward model, the proposed qRIM can exploit the redundancy among repeated measurements and shared information between tasks, facilitating relaxometry in accelerated MRI.

The influence of white matter lesions on the electric field in transcranial electric stimulation.

BACKGROUND: Transcranial direct current stimulation (tDCS) is a promising tool to enhance therapeutic efforts, for instance, after a stroke. The achieved stimulation effects exhibit high inter-subject variability, primarily driven by perturbations of the induced electric field (EF). Differences are further elevated in the aging brain due to anatomical changes such as atrophy or lesions. Informing tDCS protocols by computer-based, individualized EF simulations is a suggested measure to mitigate this variability. OBJECTIVE: While brain anatomy in general and specifically atrophy as well as stroke lesions are deemed influential on the EF in simulation studies, the influence of the uncertainty in the change of the electrical properties of the white matter due to white matter lesions (WMLs) has not been quantified yet. METHODS: A group simulation study with 88 subjects assigned into four groups of increasing lesion load was conducted. Due to the lack of information about the electrical conductivity of WMLs, an uncertainty analysis was employed to quantify the variability in the simulation when choosing an arbitrary conductivity value for the lesioned tissue. RESULTS: The contribution of WMLs to the EF variance was on average only one tenth to one thousandth of the contribution of the other modeled tissues. While the contribution of the WMLs significantly increased (p≪.01) in subjects exhibiting a high lesion load compared to low lesion load subjects, typically by a factor of 10 and above, the total variance of the EF didnot change with the lesion load. CONCLUSION: Our results suggest that WMLs do not perturb the EF globally and can thus be omitted when modeling subjects with low to medium lesion load. However, for high lesion load subjects, the omission of WMLs may yield less robust local EF estimations in the vicinity of the lesioned tissue. Our results contribute to the efforts of accurate modeling of tDCS for treatment planning.

A high-resolution multi-shell 3T diffusion magnetic resonance imaging dataset as part of the Amsterdam Ultra-high field adult lifespan database (AHEAD).

In order to further our understanding of brain function and the underlying networks, more advanced diffusion weighted magnetic resonance imaging (DWI MRI) data are essential. Here we present freely available high-resolution multi-shell multi-directional 3 Tesla (T) DWI MRI data as part of the 'Amsterdam Ultra-high field adult lifespan database' (AHEAD). The 3T DWI AHEAD dataset include 1.28mm isotropic whole brain DWI data of 49 healthy adult participants between 18 and 90 years old. The acquired data include DWIs at three non-zero b-values (48 directions, b-value 700 s/mm2; 56 directions, b-value 1000 s/mm2; 64 directions, b-value 1600 s/mm2) including a total of twelve volumes with a b-value of 0 s/mm2 (b0 volumes). In addition, eight b0 volumes with a reversed phase encoding direction were acquired to correct for distortions. To facilitate future use, the DWI data have been denoised, corrected for eddy currents, susceptibility-induced off-resonance field distortions, bias fields, and are skull stripped.

A unified 3D map of microscopic architecture and MRI of the human brain.

We present the first three-dimensional (3D) concordance maps of cyto- and fiber architecture of the human brain, combining histology, immunohistochemistry, and 7-T quantitative magnetic resonance imaging (MRI), in two individual specimens. These 3D maps each integrate data from approximately 800 microscopy sections per brain, showing neuronal and glial cell bodies, nerve fibers, and interneuronal populations, as well as ultrahigh-field quantitative MRI, all coaligned at the 200-µm scale to the stacked blockface images obtained during sectioning. These unprecedented 3D multimodal datasets are shared without any restrictions and provide a unique resource for the joint study of cell and fiber architecture of the brain, detailed anatomical atlasing, or modeling of the microscopic underpinnings of MRI contrasts.

Predictive coding during action observation - A depth-resolved intersubject functional correlation study at 7T.

While the brain regions involved in action observation are relatively well documented in humans and primates, how these regions communicate to help understand and predict actions remains poorly understood. Traditional views emphasized a feed-forward architecture in which visual features are organized into increasingly complex representations that feed onto motor programs in parietal and then premotor cortices where the matching of observed actions upon the observer's own motor programs contributes to action understanding. Predictive coding models place less emphasis on feed-forward connections and propose that feed-back connections from premotor regions back to parietal and visual neurons represent predictions about upcoming actions that can supersede visual inputs when actions become predictable, with visual input then merely representing prediction errors. Here we leverage the notion that feed-back connections target specific cortical layers to help adjudicate across these views. Specifically, we test whether observing sequences of hand actions in their natural order, which permits participants to predict upcoming actions, triggers more feed-back input to parietal regions than seeing the same actions in a scrambled sequence that hinders making predictions. Using submillimeter fMRI acquisition at 7T, we find that watching predictable sequences triggers more action-related activity (as measured using intersubject functional correlation) in the parietal cortical area PFt at depths receiving feed-back connections (layers III and V/VI) than watching the exact same actions in scrambled and hence unpredictable sequence. In addition, functional connectivity analysis performed using intersubject functional connectivity confirms that these increased action-related signals in PFt could originate from ventral premotor region BA44. This data showcases the utility of intersubject functional correlation in combination with 7T MRI to explore the architecture of social cognition under more naturalistic conditions, and provides evidence for models that emphasize the importance of feed-back connections in action prediction.

Charting human subcortical maturation across the adult lifespan with in vivo 7 T MRI.

The human subcortex comprises hundreds of unique structures. Subcortical functioning is crucial for behavior, and disrupted function is observed in common neurodegenerative diseases. Despite their importance, human subcortical structures continue to be difficult to study in vivo. Here we provide a detailed account of 17 prominent subcortical structures and ventricles, describing their approximate iron and myelin contents, morphometry, and their age-related changes across the normal adult lifespan. The results provide compelling insights into the heterogeneity and intricate age-related alterations of these structures. They also show that the locations of many structures shift across the lifespan, which is of direct relevance for the use of standard magnetic resonance imaging atlases. The results further our understanding of subcortical morphometry and neuroimaging properties, and of normal aging processes which ultimately can improve our understanding of neurodegeneration.

Motor sequences; separating the sequence from the motor. A longitudinal rsfMRI study.

In motor learning, sequence specificity, i.e. the learning of specific sequential associations, has predominantly been studied using task-based fMRI paradigms. However, offline changes in resting state functional connectivity after sequence-specific motor learning are less well understood. Previous research has established that plastic changes following motor learning can be divided into stages including fast learning, slow learning and retention. A description of how resting state functional connectivity after sequence-specific motor sequence learning (MSL) develops across these stages is missing. This study aimed to identify plastic alterations in whole-brain functional connectivity after learning a complex motor sequence by contrasting an active group who learned a complex sequence with a control group who performed a control task matched for motor execution. Resting state fMRI and behavioural performance were collected in both groups over the course of 5 consecutive training days and at follow-up after 12 days to encompass fast learning, slow learning, overall learning and retention. Between-group interaction analyses showed sequence-specific decreases in functional connectivity during overall learning in the right supplementary motor area (SMA). We found that connectivity changes in a key region of the motor network, the superior parietal cortex (SPC) were not a result of sequence-specific learning but were instead linked to motor execution. Our study confirms the sequence-specific role of SMA that has previously been identified in online task-based learning studies, and extends it to resting state network changes after sequence-specific MSL.

Mapping the human lateral geniculate nucleus and its cytoarchitectonic subdivisions using quantitative MRI.

The human lateral geniculate nucleus (LGN) of the visual thalamus is a key subcortical processing site for visual information analysis. Due to its small size and deep location within the brain, a non-invasive characterization of the LGN and its microstructurally distinct magnocellular (M) and parvocellular (P) subdivisions in humans is challenging. Here, we investigated whether structural quantitative MRI (qMRI) methods that are sensitive to underlying microstructural tissue features enable MR-based mapping of human LGN M and P subdivisions. We employed high-resolution 7 Tesla in-vivo qMRI in N = 27 participants and ultra-high resolution 7 Tesla qMRI of a post-mortem human LGN specimen. We found that a quantitative assessment of the LGN and its subdivisions is possible based on microstructure-informed qMRI contrast alone. In both the in-vivo and post-mortem qMRI data, we identified two components of shorter and longer longitudinal relaxation time (T1) within the LGN that coincided with the known anatomical locations of a dorsal P and a ventral M subdivision, respectively. Through ground-truth histological validation, we further showed that the microstructural MRI contrast within the LGN pertains to cyto- and myeloarchitectonic tissue differences between its subdivisions. These differences were based on cell and myelin density, but not on iron content. Our qMRI-based mapping strategy paves the way for an in-depth understanding of LGN function and microstructure in humans. It further enables investigations into the selective contributions of LGN subdivisions to human behavior in health and disease.

White matter microstructural changes in short-term learning of a continuous visuomotor sequence.

Efficient neural transmission is crucial for optimal brain function, yet the plastic potential of white matter (WM) has long been overlooked. Growing evidence now shows that modifications to axons and myelin occur not only as a result of long-term learning, but also after short training periods. Motor sequence learning (MSL), a common paradigm used to study neuroplasticity, occurs in overlapping learning stages and different neural circuits are involved in each stage. However, most studies investigating short-term WM plasticity have used a pre-post design, in which the temporal dynamics of changes across learning stages cannot be assessed. In this study, we used multiple magnetic resonance imaging (MRI) scans at 7 T to investigate changes in WM in a group learning a complex visuomotor sequence (LRN) and in a control group (SMP) performing a simple sequence, for five consecutive days. Consistent with behavioral results, where most improvements occurred between the two first days, structural changes in WM were observed only in the early phase of learning (d1-d2), and in overall learning (d1-d5). In LRNs, WM microstructure was altered in the tracts underlying the primary motor and sensorimotor cortices. Moreover, our structural findings in WM were related to changes in functional connectivity, assessed with resting-state functional MRI data in the same cohort, through analyses in regions of interest (ROIs). Significant changes in WM microstructure were found in a ROI underlying the right supplementary motor area. Together, our findings provide evidence for highly dynamic WM plasticity in the sensorimotor network during short-term MSL.

Quantity and quality: Normative open-access neuroimaging databases.

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

A probabilistic atlas of the human ventral tegmental area (VTA) based on 7 Tesla MRI data.

Functional magnetic resonance imaging (fMRI) BOLD signal is commonly localized by using neuroanatomical atlases, which can also serve for region of interest analyses. Yet, the available MRI atlases have serious limitations when it comes to imaging subcortical structures: only 7% of the 455 subcortical nuclei are captured by current atlases. This highlights the general difficulty in mapping smaller nuclei deep in the brain, which can be addressed using ultra-high field 7 Tesla (T) MRI. The ventral tegmental area (VTA) is a subcortical structure that plays a pivotal role in reward processing, learning and memory. Despite the significant interest in this nucleus in cognitive neuroscience, there are currently no available, anatomically precise VTA atlases derived from 7 T MRI data that cover the full region of the VTA. Here, we first provide a protocol for multimodal VTA imaging and delineation. We then provide a data description of a probabilistic VTA atlas based on in vivo 7 T MRI data.

Multi-contrast anatomical subcortical structures parcellation.

The human subcortex is comprised of more than 450 individual nuclei which lie deep in the brain. Due to their small size and close proximity, up until now only 7% have been depicted in standard MRI atlases. Thus, the human subcortex can largely be considered as terra incognita. Here, we present a new open-source parcellation algorithm to automatically map the subcortex. The new algorithm has been tested on 17 prominent subcortical structures based on a large quantitative MRI dataset at 7 Tesla. It has been carefully validated against expert human raters and previous methods, and can easily be extended to other subcortical structures and applied to any quantitative MRI dataset. In sum, we hope this novel parcellation algorithm will facilitate functional and structural neuroimaging research into small subcortical nuclei and help to chart terra incognita.

3 versus 7 Tesla magnetic resonance imaging for parcellations of subcortical brain structures in clinical settings.

7 Tesla (7T) magnetic resonance imaging holds great promise for improved visualization of the human brain for clinical purposes. To assess whether 7T is superior regarding localization procedures of small brain structures, we compared manual parcellations of the red nucleus, subthalamic nucleus, substantia nigra, globus pallidus interna and externa. These parcellations were created on a commonly used clinical anisotropic clinical 3T with an optimized isotropic (o)3T and standard 7T scan. The clinical 3T MRI scans did not allow delineation of an anatomically plausible structure due to its limited spatial resolution. o3T and 7T parcellations were directly compared. We found that 7T outperformed the o3T MRI as reflected by higher Dice scores, which were used as a measurement of interrater agreement for manual parcellations on quantitative susceptibility maps. This increase in agreement was associated with higher contrast to noise ratios for smaller structures, but not for the larger globus pallidus segments. Additionally, control-analyses were performed to account for potential biases in manual parcellations by assessing semi-automatic parcellations. These results showed a higher consistency for structure volumes for 7T compared to optimized 3T which illustrates the importance of the use of isotropic voxels for 3D visualization of the surgical target area. Together these results indicate that 7T outperforms c3T as well as o3T given the constraints of a clinical setting.