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.

7 Tesla MRI Followed by Histological 3D Reconstructions in Whole-Brain Specimens.

Post mortem magnetic resonance imaging (MRI) studies on the human brain are of great interest for the validation of in vivo MRI. It facilitates a link between functional and anatomical information available from MRI in vivo and neuroanatomical knowledge available from histology/immunocytochemistry. However, linking in vivo and post mortem MRI to microscopy techniques poses substantial challenges. Fixation artifacts and tissue deformation of extracted brains, as well as co registration of 2D histology to 3D MRI volumes complicate direct comparison between modalities. Moreover, post mortem brain tissue does not have the same physical properties as in vivo tissue, and therefore MRI approaches need to be adjusted accordingly. Here, we present a pipeline in which whole-brain human post mortem in situ MRI is combined with subsequent tissue processing of the whole human brain, providing a 3-dimensional reconstruction via blockface imaging. To this end, we adapted tissue processing procedures to allow both post mortem MRI and subsequent histological and immunocytochemical processing. For MRI, tissue was packed in a susceptibility matched solution, tailored to fit the dimensions of the MRI coil. Additionally, MRI sequence parameters were adjusted to accommodate T1 and T2∗ shortening, and scan time was extended, thereby benefiting the signal-to-noise-ratio that can be achieved using extensive averaging without motion artifacts. After MRI, the brain was extracted from the skull and subsequently cut while performing optimized blockface imaging, thereby allowing three-dimensional reconstructions. Tissues were processed for Nissl and silver staining, and co-registered with the blockface images. The combination of these techniques allows direct comparisons across modalities.

Superficial white matter imaging: Contrast mechanisms and whole-brain in vivo mapping.

Superficial white matter (SWM) contains the most cortico-cortical white matter connections in the human brain encompassing the short U-shaped association fibers. Despite its importance for brain connectivity, very little is known about SWM in humans, mainly due to the lack of noninvasive imaging methods. Here, we lay the groundwork for systematic in vivo SWM mapping using ultrahigh resolution 7 T magnetic resonance imaging. Using biophysical modeling informed by quantitative ion beam microscopy on postmortem brain tissue, we demonstrate that MR contrast in SWM is driven by iron and can be linked to the microscopic iron distribution. Higher SWM iron concentrations were observed in U-fiber-rich frontal, temporal, and parietal areas, potentially reflecting high fiber density or late myelination in these areas. Our SWM mapping approach provides the foundation for systematic studies of interindividual differences, plasticity, and pathologies of this crucial structure for cortico-cortical connectivity in humans.