The myeloarchitectonic studies on the human cerebral cortex of the Vogt-Vogt school, and their significance for the interpretation of functional neuroimaging data

Quantifying myelin density in the feline auditory cortex

The cerebral cortex comprises many distinct regions that differ in structure, function, and patterns of connectivity. Current approaches to parcellating these regions often take advantage of functional neuroimaging approaches that can identify regions involved in a particular process with reasonable spatial resolution. However, neuroanatomical biomarkers are also very useful in identifying distinct cortical regions either in addition to, or in place of functional measures. For example, differences in myelin density are thought to relate to functional differences between regions, are sensitive to individual patterns of experience, and have been shown to vary across functional hierarchies in a predictable manner. Accordingly, the current study provides quantitative stereological estimates of myelin density for each of the 13 regions that make up the feline auditory cortex. We demonstrate that significant differences can be observed between auditory cortical regions, with the highest myelin density observed in the regions that comprise the auditory core (i.e., the primary auditory cortex and anterior auditory field). Moreover, our myeloarchitectonic map suggests that myelin density varies in a hierarchical fashion that conforms to the traditional model of spatial organization in auditory cortex. Taken together, these results establish myelin as a useful biomarker for parcellating auditory cortical regions, and provide detailed estimates against which other, less invasive methods of quantifying cortical myelination may be compared.

Quantitative MRI at 7-Tesla reveals novel frontocortical myeloarchitecture anomalies in major depressive disorder

Whereas meta-analytical data highlight abnormal frontocortical macrostructure (thickness/surface area/volume) in Major Depressive Disorder (MDD), the underlying microstructural processes remain uncharted, due to the use of conventional MRI scanners and acquisition techniques. We uniquely combined Ultra-High Field MRI at 7.0 Tesla with Quantitative Imaging to map intracortical myelin (proxied by longitudinal relaxation time T1) and iron concentration (proxied by transverse relaxation time T2*), microstructural processes deemed particularly germane to cortical macrostructure. Informed by meta-analytical evidence, we focused specifically on orbitofrontal and rostral anterior cingulate cortices among adult MDD patients (N = 48) and matched healthy controls (HC; N = 10). Analyses probed the association of MDD diagnosis and clinical profile (severity, medication use, comorbid anxiety disorders, childhood trauma) with aforementioned microstructural properties. MDD diagnosis (p's < 0.05, Cohen's D = 0.55-0.66) and symptom severity (p's < 0.01, r = 0.271-0.267) both related to decreased intracortical myelination (higher T1 values) within the lateral orbitofrontal cortex, a region tightly coupled to processing negative affect and feelings of sadness in MDD. No relations were found with local iron concentrations. These findings allow uniquely fine-grained insights on frontocortical microstructure in MDD, and cautiously point to intracortical demyelination as a possible driver of macroscale cortical disintegrity in MDD.

The cytoarchitectonic landscape revealed by deep learning method facilitated precise positioning in mouse neocortex

Neocortex is a complex structure with different cortical sublayers and regions. However, the precise positioning of cortical regions can be challenging due to the absence of distinct landmarks without special preparation. To address this challenge, we developed a cytoarchitectonic landmark identification pipeline. The fluorescence micro-optical sectioning tomography method was employed to image the whole mouse brain stained by general fluorescent nucleotide dye. A fast 3D convolution network was subsequently utilized to segment neuronal somas in entire neocortex. By approach, the cortical cytoarchitectonic profile and the neuronal morphology were analyzed in 3D, eliminating the influence of section angle. And the distribution maps were generated that visualized the number of neurons across diverse morphological types, revealing the cytoarchitectonic landscape which characterizes the landmarks of cortical regions, especially the typical signal pattern of barrel cortex. Furthermore, the cortical regions of various ages were aligned using the generated cytoarchitectonic landmarks suggesting the structural changes of barrel cortex during the aging process. Moreover, we observed the spatiotemporally gradient distributions of spindly neurons, concentrated in the deep layer of primary visual area, with their proportion decreased over time. These findings could improve structural understanding of neocortex, paving the way for further exploration with this method.

Transcriptomic contributions to a modern cytoarchitectonic parcellation of the human cerebral cortex

Transcriptomic contributions to the anatomical, functional, and network layout of the human cerebral cortex (HCC) have become a major interest in cognitive and systems neuroscience. Here, we tested if transcriptomic differences support a modern, algorithmic cytoarchitectonic parcellation of HCC. Using a data-driven approach, we identified a sparse subset of genes that differentially contributed to the cytoarchitectonic parcellation of HCC. A combined metric of cortical thickness and myelination (CT/M ratio), as well as cell density, correlated with gene expression. Enrichment analyses showed that genes specific to the cytoarchitectonic parcellation of the HCC were related to molecular functions such as transmembrane transport and ion channel activity. Together, the relationship between transcriptomics and cytoarchitecture bridges the gap among (i) gradients at the macro-scale (including thickness and myelination), (ii) areas at the meso-scale, and (iii) cell density at the microscale, as well as supports the recently proposed cortical spectrum theory and structural model.

7T MP2RAGE for cortical myelin segmentation: Impact of aging

BACKGROUND

Myelin and iron are major contributors to the cortical MR signal. The aim of this study was to investigate 1. Can MP2RAGE-derived contrasts at 7T in combination with k-means clustering be used to distinguish between heavily and sparsely myelinated layers in cortical gray matter (GM)? 2. Does this approach provide meaningful biological information?

METHODS

The following contrasts were generated from the 7T MP2RAGE images from 45 healthy controls (age: 19-75, f/m = 23/22) from the ATAG data repository: 1. T1 weighted image (UNI). 2. T1 relaxation image (T1map). 3. INVC/T1map ratio (RATIO). K-means clustering identified 6 clusters/tissue maps (csf, csf/gm-transition, wm, wm/gm transition, heavily myelinated cortical GM (dGM), sparsely myelinated cortical GM (sGM)). These tissue maps were then processed with SPM/DARTEL (volume-based analyses) and Freesurfer (surface-based analyses) and dGM and sGM volume/thickness of young adults (n = 27, 19-27 years) compared to those of older adults (n = 18, 42-75 years) at p<0.001 uncorrected.

RESULTS

The resulting maps showed good agreement with histological maps in the literature. Volume- and surface analyses found age-related dGM loss/thinning in the mid-posterior cingulate and parahippocampal/entorhinal gyrus and age-related sGM losses in lateral, mesial and orbitofrontal frontal, insular cortex and superior temporal gyrus.

CONCLUSION

The MP2RAGE derived UNI, T1map and RATIO contrasts can be used to identify dGM and sGM. Considering the close relationship between cortical myelo- and cytoarchitecture, the findings reported here indicate that this new technique might provide new insights into the nature of cortical GM loss in physiological and pathological conditions.

Local gradient analysis of human brain function using the Vogt-Bailey Index

In this work, we take a closer look at the Vogt-Bailey (VB) index, proposed in Bajada et al. (NeuroImage 221:117140, 2020) as a tool for studying local functional homogeneity in the human cortex. We interpret the VB index in terms of the minimum ratio cut, a scaled cut-set weight that indicates whether a network can easily be disconnected into two parts having a comparable number of nodes. In our case, the nodes of the network consist of a brain vertex/voxel and its neighbours, and a given edge is weighted according to the affinity of the nodes it connects (as reflected by the modified Pearson correlation between their fMRI time series). Consequently, the minimum ratio cut quantifies the degree of small-scale similarity in brain activity: the greater the similarity, the 'heavier' the edges and the more difficult it is to disconnect the network, hence the higher the value of the minimum ratio cut. We compare the performance of the VB index with that of the Regional Homogeneity (ReHo) algorithm, commonly used to assess whether voxels in close proximity have synchronised fMRI signals, and find that the VB index is uniquely placed to detect sharp changes in the (local) functional organization of the human cortex.

Cortical myelin and thickness mapping provide insights into whole-brain tumor burden in diffuse midline glioma

Systemic infiltration is a hallmark of diffuse midline glioma pathogenesis, which can trigger distant disturbances in cortical structure. However, the existence and effects of these changes have been underexamined. This study aimed to investigate whole-brain cortical myelin and thickness alternations induced by diffuse midline glioma. High-resolution T1- and T2-weighted images were acquired from 90 patients with diffuse midline glioma with H3 K27-altered and 64 patients with wild-type and 86 healthy controls. Cortical thickness and myelin content was calculated using Human Connectome Project pipeline. Significant differences in cortical thickness and myelin content were detected among groups. Short-term survival prediction model was constructed using automated machine learning. Compared with healthy controls, diffuse midline glioma with H3 K27-altered patients showed significantly reduced cortical myelin in bilateral precentral gyrus, postcentral gyrus, insular, parahippocampal gyrus, fusiform gyrus, and cingulate gyrus, whereas diffuse midline glioma with H3 K27 wild-type patients exhibited well-preserved myelin content. Furtherly, when comparing diffuse midline glioma with H3 K27-altered and diffuse midline glioma with H3 K27 wild-type, the decreased cortical thickness in parietal and occipital regions along with demyelination in medial orbitofrontal cortex was observed in diffuse midline glioma with H3 K27-altered. Notably, a combination of cortical features and tumor radiomics allowed short-term survival prediction with accuracy 0.80 and AUC 0.84. These findings may aid clinicians in tailoring therapeutic approaches based on cortical characteristics, potentially enhancing the efficacy of current and future treatment modalities.

Analyzing the cortical fine structure as revealed by ex-vivo anatomical MRI

The human cortex has a rich fiber structure as revealed by myelin-staining of histological slices. Myelin also contributes to the image contrast in Magnetic Resonance Imaging (MRI). Recent advances in Magnetic Resonance (MR) scanner and imaging technology allowed the acquisition of an ex-vivo data set at an isotropic resolution of 100 µm. This study focused on a computational analysis of this data set with the aim of bridging between histological knowledge and MRI-based results. This work highlights: (1) the design and implementation of a processing chain that extracts intracortical features from a high-resolution MR image; (2) a demonstration of the correspondence between MRI-based cortical intensity profiles and the myelo-architectonic layering of the cortex; (3) the characterization and classification of four basic myelo-architectonic profile types; (4) the distinction of cortical regions based on myelo-architectonic features; and (5) the segmentation of cortical modules in the entorhinal cortex.

High-Resolution Tractography Protocol to Investigate the Pathways between Human Mediodorsal Thalamic Nucleus and Prefrontal Cortex

Animal studies have established that the mediodorsal nucleus (MD) of the thalamus is heavily and reciprocally connected with all areas of the prefrontal cortex (PFC). In humans, however, these connections are difficult to investigate. High-resolution imaging protocols capable of reliably tracing the axonal tracts linking the human MD with each of the PFC areas may thus be key to advance our understanding of the variation, development, and plastic changes of these important circuits, in health and disease. Here, we tested in adult female and male humans the reliability of a new reconstruction protocol based on in vivo diffusion MRI to trace, measure, and characterize the fiber tracts interconnecting the MD with 39 human PFC areas per hemisphere. Our protocol comprised the following three components: (1) defining regions of interest; (2) preprocessing diffusion data; and, (3) modeling white matter tracts and tractometry. This analysis revealed largely separate PFC territories of reciprocal MD-PFC tracts bearing striking resemblance with the topographic layout observed in macaque connection-tracing studies. We then examined whether our protocol could reliably reconstruct each of these MD-PFC tracts and their profiles across test and retest sessions. Results revealed that this protocol was able to trace and measure, in both left and right hemispheres, the trajectories of these 39 area-specific axon bundles with good-to-excellent test-retest reproducibility. This protocol, which has been made publicly available, may be relevant for cognitive neuroscience and clinical studies of normal and abnormal PFC function, development, and plasticity.SIGNIFICANCE STATEMENT Reciprocal MD-PFC interactions are critical for complex human cognition and learning. Reliably tracing, measuring and characterizing MD-PFC white matter tracts using high-resolution noninvasive methods is key to assess individual variation of these systems in humans. Here, we propose a high-resolution tractography protocol that reliably reconstructs 39 area-specific MD-PFC white matter tracts per hemisphere and quantifies structural information from diffusion MRI data. This protocol revealed a detailed mapping of thalamocortical and corticothalamic MD-PFC tracts in four different PFC territories (dorsal, medial, orbital/frontal pole, inferior frontal) showing structural connections resembling those observed in tracing studies with macaques. Furthermore, our automated protocol revealed high test-retest reproducibility and is made publicly available, constituting a step forward in mapping human MD-PFC circuits in clinical and academic research.

Transcriptional regulatory network for neuron-glia interactions and its implication for DYT6 dystonia

Advances in sequencing technologies have identified novel genes associated with inherited forms of dystonia, providing valuable insights into its genetic basis and revealing diverse genetic pathways and mechanisms involved in its pathophysiology. Since identifying genetic variation in the transcription factor coding THAP1 gene linked to isolated dystonia, numerous investigations have employed transcriptomic studies in DYT-THAP1 models to uncover pathogenic molecular mechanisms underlying dystonia. This review examines key findings from transcriptomic studies conducted on in vivo and in vitro DYT-THAP1 models, which demonstrate that the THAP1-regulated transcriptome is diverse and cell-specific, yet it is bound and co-regulated by a common set of proteins. Prominent among its functions, THAP1 and its co-regulatory network target molecular pathways critical for generating myelinating oligodendrocytes that ensheath axons and generate white matter in the central nervous system. Several lines of investigation have demonstrated the importance of myelination and oligodendrogenesis in motor function during development and in adults, emphasizing the non-cell autonomous contributions of glial cells to neural circuits involved in motor function. Further research on the role of myelin abnormalities in motor deficits in DYT6 models will enhance our understanding of axon-glia interactions in dystonia pathophysiology and provide potential therapeutic interventions targeting these pathways.

Controversies and progress on standardization of large-scale brain network nomenclature

Progress in scientific disciplines is accompanied by standardization of terminology. Network neuroscience, at the level of macroscale organization of the brain, is beginning to confront the challenges associated with developing a taxonomy of its fundamental explanatory constructs. The Workgroup for HArmonized Taxonomy of NETworks (WHATNET) was formed in 2020 as an Organization for Human Brain Mapping (OHBM)-endorsed best practices committee to provide recommendations on points of consensus, identify open questions, and highlight areas of ongoing debate in the service of moving the field toward standardized reporting of network neuroscience results. The committee conducted a survey to catalog current practices in large-scale brain network nomenclature. A few well-known network names (e.g., default mode network) dominated responses to the survey, and a number of illuminating points of disagreement emerged. We summarize survey results and provide initial considerations and recommendations from the workgroup. This perspective piece includes a selective review of challenges to this enterprise, including (1) network scale, resolution, and hierarchies; (2) interindividual variability of networks; (3) dynamics and nonstationarity of networks; (4) consideration of network affiliations of subcortical structures; and (5) consideration of multimodal information. We close with minimal reporting guidelines for the cognitive and network neuroscience communities to adopt.

Human microsaccade cueing modulation in visual- and memory-delay saccade tasks after theta burst transcranial magnetic stimulation over the frontal eye field

Microsaccades that occur during periods of fixation are modulated by various cognitive processes and have an impact on visual processing. A network of brain areas is involved in microsaccade generation, including the superior colliculus and frontal eye field (FEF) which are involved in modulating microsaccade rate and direction after the appearance of a visual cue (referred to as microsaccade cueing modulation). Although the neural mechanisms underlying microsaccade cueing modulations have been demonstrated in monkeys, limited research has investigated a causal role of these areas in humans. By applying continuous theta-burst transcranial magnetic stimulation (cTBS) over the right FEF and vertex, we investigated the role of human FEF in modulating microsaccade responses after the appearance of a visual target in a visual- and memory-delay saccade task. After target appearance, microsaccade rate was initially suppressed but then increased in both cTBS conditions. More importantly, in the visual-delay task, microsaccades after target appearance were directed to the ipsilateral side more often with FEF, compared to vertex stimulation. Moreover, microsaccades were directed towards the target location, then to the opposite location of the target in both tasks, with larger effects in the visual-, compared to, memory-delay task. This microsaccade direction modulation was delayed after FEF stimulation in the memory-delay task. Overall, some microsaccade cueing modulations were moderately disrupted after FEF cTBS, suggesting a causal role for involvement of the human FEF in microsaccade generation after presentation of salient stimuli.

A new 3D myeloarchitectonic map of the human neocortex based on data from the Vogt-Vogt school

During the period extending from 1900 to 1970, Oskar and Cécile Vogt and their numerous collaborators ('the Vogt-Vogt school') published a large number of studies on the myeloarchitecture of the human cerebral cortex. During the last decade, we have concerned ourselves with a detailed meta-analysis of these now almost totally forgotten studies, with the aim to bringing them into the modern era of science. This scrutiny yielded inter alia a myeloarchitectonic map of the human neocortex, showing a parcellation into 182 areas (Nieuwenhuys et al. in Brain Struct Funct 220:2551-2573, 2015; Erratum in Brain Struct Funct 220: 3753-3755, 2015). This map, termed 2D'15, which is based on data derived from all of the 20 publications constituting the myeloarchitectonic legacy of the Vogt-Vogt school, has the limitation that it is two-dimensional i.e. it shows only the parts of the cortex exposed at the free surface of the cerebral hemispheres and not the extensive stretches of cortex hidden in the cortical sulci. However, a limited set of data, derived from four of the 20 publications available, has enabled us to create a 3D map, showing the myeloarchitectonic parcellation of the entire human neocortex. This map, designated as 3D'23, contains 182 areas: 64 frontal, 30 parietal, 6 insular, 19 occipital and 63 temporal. We have also prepared a 2D version (2D'23), of this 3D'23 map to serve as a link between the latter and our original 2D'15 map. Detailed comparison of the parcellations visualized in our three maps (2D'15, 2D'23 and 3D'23) warrants the conclusion that our new 3D'23 map may be considered as representative for the entire myeloarchitectural legacy of the Vogt-Vogt School. Hence it is now possible to compare the rich amount of myeloarchitectonic data assembled by that school directly with the results of current 3D analyses of the architecture of the human cortex, such as the meticulous quantitative cyto- and receptor architectonic studies of Zilles, Amunts and their numerous associates (Amunts et al. in Science 369:988-992, 2020), and the multimodal parcellation of the human cortex based on magnetic resonance images from the Human Connectome Project, performed by Glasser et al. in Nature 536:171-178, 2016).

The intracortical myelin content of impulsive choices: results from T1- and T2-weighted MRI myelin mapping

Delay discounting (DD) refers to a phenomenon that humans tend to choose small-sooner over large-later rewards during intertemporal choices. Steep discounting of delayed outcome is related to a variety of maladaptive behaviors and is considered as a transdiagnostic process across psychiatric disorders. Previous studies have investigated the association between brain structure (e.g. gray matter volume) and DD; however, it is unclear whether the intracortical myelin (ICM) influences DD. Here, based on a sample of 951 healthy young adults drawn from the Human Connectome Project, we examined the relationship between ICM, which was measured by the contrast of T1w and T2w images, and DD and further tested whether the identified associations were mediated by the regional homogeneity (ReHo) of brain spontaneous activity. Vertex-wise regression analyses revealed that steeper DD was significantly associated with lower ICM in the left temporoparietal junction (TPJ) and right middle-posterior cingulate cortex. Region-of-interest analysis revealed that the ReHo values in the left TPJ partially mediated the association of its myelin content with DD. Our findings provide the first evidence that cortical myelination is linked with individual differences in decision impulsivity and suggest that the myelin content affects cognitive performances partially through altered local brain synchrony.

Myeloarchitectonic cortical parcellation data for contemporary neuroimaging - the Vogt-Vogt legacy in the 21st century

Obtaining precise and detailed parcellations of the human brain has been a major focus of neuroscience research. Here, we present a multimodal dataset, MYATLAS, based on histology-derived myeloarchitectonic parcellations for use with contemporary neuroimaging analyses software. The core of MYATLAS is a novel 3D neocortical, surface-based atlas derived from legacy myeloarchitectonic histology studies. Additionally, we provide digitized quantitative laminar profiles of intracortical myelin content derived from postmortem photometric data, cross-correlated with in vivo myeloarchitectonic features obtained by quantitative MRI mapping. Moreover, congregated, digitized and quality-improved Vogt-Vogt legacy histology data is made available. Finally, to allow for cross-modality correlations, maps of quantitative myelin estimates and corresponding von Economo-Koskinas' cytoarchitectonic features are also included. We share all necessary surface and volume-based registration files as well as shell scripts to facilitate applications of MYATLAS to future in vivo MRI studies.

High-resolution magnetization-transfer imaging of post-mortem marmoset brain: Comparisons with relaxometry and histology

Cell membranes and macromolecules or paramagnetic compounds interact with water proton spins, which modulates magnetic resonance imaging (MRI) contrast providing information on tissue composition. For a further investigation, quantitative magnetization transfer (qMT) parameters (at 3T), including the ratio of the macromolecular and water proton pools, F, and the exchange-rate constant as well as the (observed) longitudinal and the effective transverse relaxation rates (at 3T and 7T), R₁^obs and R₂^*, respectively, were measured at high spatial resolution (200 µm) in a slice of fixed marmoset brain and compared to histology results obtained with Gallyas' myelin stain and Perls' iron stain. R₁^obs and R₂^* were linearly correlated with the iron content for the entire slice, whereas distinct differences were obtained between gray and white matter for correlations of relaxometry and qMT parameters with myelin content. The combined results suggest that the macromolecular pool interacting with water consists of myelin and (less efficient) non-myelin contributions. Despite strong correlation of F and R₁^obs, none of these parameters was uniquely specific to myelination. Due to additional sensitivity to iron stores, R₁^obs and R₂^* were more sensitive for depicting microstructural differences between cortical layers than F.

Towards automated in vivo parcellation of the human cerebral cortex using supervised classification of magnetic resonance fingerprinting residuals

BACKGROUND

Accurate parcellation of the cerebral cortex in an individual is a guide to its underlying organization. The most promising in vivo quantitative magnetic resonance (MR)-based microstructural cortical mapping methods are yet to achieve a level of parcellation accuracy comparable to quantitative histology.

METHODS

We scanned 6 participants using a 3D echo-planar imaging MR fingerprinting (EPI-MRF) sequence on a 7T Siemens scanner. After projecting MRF signals to the individual-specific inflated model of the cortical surface, normalized autocorrelations of MRF residuals of vertices of 8 microstructurally distinct areas (BA1, BA2, BA4a, BA6, BA44, BA45, BA17, and BA18) from 3 cortical regions were used as feature vector inputs into linear support vector machine (SVM), radial basis function SVM (RBF-SVM), random forest, and k-nearest neighbors supervised classification algorithms. The algorithms' prediction performance was compared using: (i) features from each vertex or (ii) features from neighboring vertices.

RESULTS

The neighborhood-based RBF-SVM classifier achieved the highest prediction score of 0.85 for classification of MRF residuals in the central region from a held-out participant.

CONCLUSIONS

We developed an automated method of cortical parcellation using a combination of MR fingerprinting residual analysis and machine learning classification. Our findings provide the basis for employing unsupervised learning algorithms for whole-cortex structural parcellation in individuals.

Mapping myelin in white matter with T1-weighted/T2-weighted maps: discrepancy with histology and other myelin MRI measures

The ratio of T1-weighted/T2-weighted magnetic resonance images (T1w/T2w MRI) has been successfully applied at the cortical level since 2011 and is now one of the most used myelin mapping methods. However, no reports have explored the histological validity of T1w/T2w myelin mapping in white matter. Here we compare T1w/T2w with ex vivo postmortem histology and in vivo MRI methods, namely quantitative susceptibility mapping (QSM) and multi-echo T2 myelin water fraction (MWF) mapping techniques. We report a discrepancy between T1w/T2w myelin maps of the human corpus callosum and the histology and analyse the putative causes behind such discrepancy. T1w/T2w does not positively correlate with Luxol Fast Blue (LFB)-Optical Density but shows a weak to moderate, yet significant, negative correlation. On the contrary, MWF is strongly and positively correlated with LFB, whereas T1w/T2w and MWF maps are weakly negatively correlated. The discrepancy between T1w/T2w MRI maps, MWF and histological myelin maps suggests caution in using T1w/T2w as a white matter mapping method at the callosal level. While T1w/T2w imaging may correlate with myelin content at the cortical level, it is not a specific method to map myelin density in white matter.

Myeloarchitectonic maps of the human cerebral cortex registered to surface and sections of a standard atlas brain

C. and O. Vogt had set up a research program with the aim of establishing a detailed cartography of the medullary fiber distribution of the human brain. As part of this program, around 200 cortical fields were differentiated based on their myeloarchitectural characteristics and mapped with regard to their exact location in the isocortex. The typical features were graphically documented and classified by a sophisticated linguistic coding. Their results have only recently received adequate attention and applications. The reasons for the revival of this spectrum of their research include interest in the myeloarchitecture of the cortex as a differentiating feature of the cortex architecture and function, as well as the importance for advanced imaging methodologies, particularly tractography and molecular imaging. Here, we describe our approach to exploit the original work of the Vogts and their co-workers to construct a myeloarchitectonic map that is referenced to the Atlas of the Human Brain (AHB) in standard space. We developed a semi-automatic pipeline for processing and integrating the various original maps into a single coherent map. To optimize the precision of the registration between the published maps and the AHB, we augmented the maps with topographic landmarks of the brains that were originally analyzed. Registration of all maps into the AHB opened several possibilities. First, for the majority of the fields, multiple maps from different authors are available, which allows for sophisticated statistical integration, for example, unification with a label-fusion technique. Second, each field in the myeloarchitectonic surface map can be visualized on the myelin-stained cross-section of the AHB at the best possible correspondence. The features of each field can be correlated with the fiber-stained cross-sections in the AHB and with the extensive published materials from the Vogt school and, if necessary, corrected. Third, mapping to the AHB allows the relationship between fiber characteristics of the cortex and the subcortex to be examined. Fourth, the cytoarchitectonic maps from Brodmann and von Economo and Koskinas, which are also registered to the AHB, can be compared. This option allows the study of the correspondence between cyto- and myeloarchitecture in each field. Finally, by using our "stripe" technology - where any other feature registered to the same space can be directly compared owing to the linear and parallel representation of the correlated cortex segments - this map becomes part of a multidimensional co-registration platform.

Effects of non-modifiable risk factors of Alzheimer's disease on intracortical myelin content

BACKGROUND

Non-modifiable risk factors of Alzheimer's disease (AD) have lifelong effects on cortical integrity that could be mitigated if identified at early stages. However, it remains unknown whether cortical microstructure is affected in older individuals with non-modifiable AD risk factors and whether altered cortical tissue integrity produces abnormalities in brain functional networks in this AD-risk population.

METHODS

Using relative T1w/T2w (rT1w/T2w) ratio maps, we have compared tissue integrity of normal-appearing cortical GM between controls and cognitively normal older adults with either APOE4 (N = 50), with a first-degree family history (FH) of AD (N = 52), or with the co-occurrence of both AD risk factors (APOE4+FH) (N = 35). Additionally, individuals with only one risk factor (APOE4 or FH) were combined into one group (N = 102) and compared with controls. The same number of controls matched in age, sex, and years of education was employed for each of these comparisons. Group differences in resting state functional connectivity (rs-FC) patterns were also investigated, using as FC seeds those cortical regions showing significant changes in rT1w/T2w ratios.

RESULTS

Overall, individuals with non-modifiable AD risk factors exhibited significant variations in rT1w/T2w ratios compared to controls, being APOE4 and APOE4+FH at opposite ends of a continuum. The co-occurrence of APOE4 and FH was further accompanied by altered patterns of rs-FC.

CONCLUSIONS

These findings may have practical implications for early detection of cortical abnormalities in older populations with APOE4 and/or FH of AD and open new avenues to monitor changes in cortical tissue integrity associated with non-modifiable AD risk factors.

Stress and the brain: Emotional support mediates the association between myelination in the right supramarginal gyrus and perceived chronic stress

Perceived stress, which refers to people's evaluation of a stressful event and their ability to cope with it, has emerged as a stable predictor for physical and mental health outcomes. Increasing evidence has suggested the buffering effect of social support on perceived stress. Although previous studies have investigated the brain structural features (e.g., gray matter volume) associated with perceived stress, less is known about the association between perceived chronic stress and intra-cortical myelin (ICM), which is an important microstructure of brain and is essential for healthy brain functions, and the role of social support in this association. Using a sample of 1076 healthy young adults drawn from the Human Connectome Project, we quantified the ICMby the contrast of T1w and T2w images and examined its association with perceived chronic stress during the last month and social support. Behavioral results showed that perceived chronic stress was negatively associated with both emotional support and instrumental support. Vertex-wise multiple regression analyses revealed that higher level of perceived chronic stress was significantly associated with lower ICM content of a cluster in the right supramarginal gyrus (rSMG). Interestingly, the emotional support, but not the instrumental support, significantly mediated the association of perceived chronic stress with ICM in the rSMG. Overall, the present study provides novel evidence for the cortical myelination of perceived chronic stress in humans and highlights the essential role of the rSMG in perceived chronic stress and emotional support.

COnstrained Reference frame diffusion TEnsor Correlation Spectroscopic (CORTECS) MRI: A practical framework for high-resolution diffusion tensor distribution imaging

High-resolution imaging studies have consistently shown that in cortical tissue water diffuses preferentially along radial and tangential orientations with respect to the cortical surface, in agreement with histology. These dominant orientations do not change significantly even if the relative contributions from microscopic water pools to the net voxel signal vary across experiments that use different diffusion times, b-values, TEs, and TRs. With this in mind, we propose a practical new framework for imaging non-parametric diffusion tensor distributions (DTDs) by constraining the microscopic diffusion tensors of the DTD to be diagonalized using the same orthonormal reference frame of the mesoscopic voxel. In each voxel, the constrained DTD (cDTD) is completely determined by the correlation spectrum of the microscopic principal diffusivities associated with the axes of the voxel reference frame. Consequently, all cDTDs are inherently limited to the domain of positive definite tensors and can be reconstructed efficiently using Inverse Laplace Transform methods. Moreover, the cDTD reconstruction can be performed using only data acquired efficiently with single diffusion encoding, although it also supports datasets with multiple diffusion encoding. In tissues with a well-defined architecture, such as the cortex, we can further constrain the cDTD to contain only cylindrically symmetric diffusion tensors and measure the 2D correlation spectra of principal diffusivities along the radial and tangential orientation with respect to the cortical surface. To demonstrate this framework, we perform numerical simulations and analyze high-resolution dMRI data from a fixed macaque monkey brain. We estimate 2D cDTDs in the cortex and derive, in each voxel, the marginal distributions of the microscopic principal diffusivities, the corresponding distributions of the microscopic fractional anisotropies and mean diffusivities along with their 2D correlation spectra to quantify the cDTD shape-size characteristics. Signal components corresponding to specific bands in these cDTD-derived spectra show high specificity to cortical laminar structures observed with histology. Our framework drastically simplifies the measurement of non-parametric DTDs in high-resolution datasets with mesoscopic voxel sizes much smaller than the radius of curvature of the underlying anatomy, e.g., cortical surface, and can be applied retrospectively to analyze existing diffusion MRI data from fixed cortical tissues.

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

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

Cortical myelin profile variations in healthy aging brain: A T1w/T2w ratio study

Demyelination is observed in both healthy aging and age-related neurodegenerative disorders. While the significance of myelin within the cortex is well acknowledged, studies focused on intracortical demyelination and depth-specific structural alterations in normal aging are lacking. Using the recently available Human Connectome Project Aging dataset, we investigated intracortical myelin in a normal aging population using the T1w/T2w ratio. To capture the fine changes across cortical depths, we employed a surface-based approach by constructing cortical profiles traveling perpendicularly through the cortical ribbon and sampling T1w/T2w values. The curvatures of T1w/T2w cortical profiles may be influenced by differences in local myeloarchitecture and other tissue properties, which are known to vary across cortical regions. To quantify the shape of these profiles, we parametrized the level of curvature using a nonlinearity index (NLI) that measures the deviation of the profile from a straight line. We showed that NLI exhibited a steep decline in aging that was independent of local cortical thinning. Further examination of the profiles revealed that lower T1w/T2w near the gray-white matter boundary and superficial cortical depths were major contributors to the apparent NLI variations with age. These findings suggest that demyelination and changes in other T1w/T2w related tissue properties in normal aging may be depth-specific and highlight the potential of NLI as a unique marker of microstructural alterations within the cerebral cortex.

Cortical adaptation of the night monkey to a nocturnal niche environment: a comparative non-invasive T1w/T2w myelin study

Night monkeys (Aotus) are the only genus of monkeys within the Simian lineage that successfully occupy a nocturnal environmental niche. Their behavior is supported by their sensory organs' distinctive morphological features; however, little is known about their evolutionary adaptations in sensory regions of the cerebral cortex. Here, we investigate this question by exploring the cortical organization of night monkeys using high-resolution in-vivo brain MRI and comparative cortical-surface T1w/T2w myeloarchitectonic mapping. Our results show that the night monkey cerebral cortex has a qualitatively similar but quantitatively different pattern of cortical myelin compared to the diurnal macaque and marmoset monkeys. T1w/T2w myelin and its gradient allowed us to parcellate high myelin areas, including the middle temporal complex (MT +) and auditory cortex, and a low-myelin area, Brodmann area 7 (BA7) in the three species, despite species differences in cortical convolutions. Relative to the total cortical-surface area, those of MT + and the auditory cortex are significantly larger in night monkeys than diurnal monkeys, whereas area BA7 occupies a similar fraction of the cortical sheet in all three species. We propose that the selective expansion of sensory areas dedicated to visual motion and auditory processing in night monkeys may reflect cortical adaptations to a nocturnal environment.

HOA2.0-ComPaRe: A next generation Harvard-Oxford Atlas comparative parcellation reasoning method for human and macaque individual brain parcellation and atlases of the cerebral cortex

Comparative structural neuroanatomy is a cornerstone for understanding human brain structure and function. A parcellation framework that relates systematically to fundamental principles of histological organization is an essential step in generating structural comparisons between species. In the present investigation, we developed a comparative parcellation reasoning system (ComPaRe), which is a formal ontological system in human and non-human primate brains based on the cortical cytoarchitectonic mapping used for both species as detailed by Brodmann. ComPaRe provides a theoretical foundation for mapping neural systems in humans and other species using neuroimaging. Based on this approach, we revised the methodology of the original Harvard-Oxford Atlas (HOA) system of brain parcellation to produce a comparative framework for the human (hHOA) and the rhesus monkey (mHOA) brains, which we refer to as HOA2.0-ComPaRe. In addition, we used dedicated segmentation software in the publicly available 3D Slicer platform to parcellate an individual human and rhesus monkey brain. This method produces quantitative morphometric parcellations in the individual brains. Based on these parcellations we created a representative template and 3D brain atlas for the two species, each based on a single subject. Thus, HOA2.0-ComPaRe provides a theoretical foundation for mapping neural systems in humans and other species using neuroimaging, while also representing a significant revision of the original human and macaque monkey HOA parcellation schemas. The methodology and atlases presented here can be used in basic and clinical neuroimaging for morphometric (volumetric) analysis, further generation of atlases, as well as localization of function and structural lesions.

The large numbers of minicolumns in the primary visual cortex of humans, chimpanzees and gorillas are related to high visual acuity

Minicolumns are thought to be a fundamental neural unit in the neocortex and their replication may have formed the basis of the rapid cortical expansion that occurred during primate evolution. We sought evidence of minicolumns in the primary visual cortex (V-1) of three great apes, three rodents and representatives from three other mammalian orders: Eulipotyphla (European hedgehog), Artiodactyla (domestic pig) and Carnivora (ferret). Minicolumns, identified by the presence of a long bundle of radial, myelinated fibers stretching from layer III to the white matter of silver-stained sections, were found in the human, chimpanzee, gorilla and guinea pig V-1. Shorter bundles confined to one or two layers were found in the other species but represent modules rather than minicolumns. The inter-bundle distance, and hence density of minicolumns, varied systematically both within a local area that might represent a hypercolumn but also across the whole visual field. The distance between all bundles had a similar range for human, chimpanzee, gorilla, ferret and guinea pig: most bundles were 20-45 μm apart. By contrast, the space between bundles was greater for the hedgehog and pig (20-140 μm). The mean density of minicolumns was greater in tangential sections of the gorilla and chimpanzee (1,243-1,287 bundles/mm2) than in human (314-422 bundles/mm2) or guinea pig (643 bundles/mm2). The minicolumnar bundles did not form a hexagonal lattice but were arranged in thin curving and branched bands separated by thicker bands of neuropil/somata. Estimates of the total number of modules/minicolumns within V-1 were strongly correlated with visual acuity.

Quantitative MRI reveals differences in striatal myelin in children with DLD

Developmental language disorder (DLD) is a common neurodevelopmental disorder characterised by receptive or expressive language difficulties or both. While theoretical frameworks and empirical studies support the idea that there may be neural correlates of DLD in frontostriatal loops, findings are inconsistent across studies. Here, we use a novel semiquantitative imaging protocol - multi-parameter mapping (MPM) - to investigate microstructural neural differences in children with DLD. The MPM protocol allows us to reproducibly map specific indices of tissue microstructure. In 56 typically developing children and 33 children with DLD, we derived maps of (1) longitudinal relaxation rate R1 (1/T1), (2) transverse relaxation rate R2* (1/T2*), and (3) Magnetization Transfer saturation (MTsat). R1 and MTsat predominantly index myelin, while R2* is sensitive to iron content. Children with DLD showed reductions in MTsat values in the caudate nucleus bilaterally, as well as in the left ventral sensorimotor cortex and Heschl's gyrus. They also had globally lower R1 values. No group differences were noted in R2* maps. Differences in MTsat and R1 were coincident in the caudate nucleus bilaterally. These findings support our hypothesis of corticostriatal abnormalities in DLD and indicate abnormal levels of myelin in the dorsal striatum in children with DLD.

An Open MRI Dataset For Multiscale Neuroscience

Multimodal neuroimaging grants a powerful window into the structure and function of the human brain at multiple scales. Recent methodological and conceptual advances have enabled investigations of the interplay between large-scale spatial trends (also referred to as gradients) in brain microstructure and connectivity, offering an integrative framework to study multiscale brain organization. Here, we share a multimodal MRI dataset for Microstructure-Informed Connectomics (MICA-MICs) acquired in 50 healthy adults (23 women; 29.54 ± 5.62 years) who underwent high-resolution T1-weighted MRI, myelin-sensitive quantitative T1 relaxometry, diffusion-weighted MRI, and resting-state functional MRI at 3 Tesla. In addition to raw anonymized MRI data, this release includes brain-wide connectomes derived from (i) resting-state functional imaging, (ii) diffusion tractography, (iii) microstructure covariance analysis, and (iv) geodesic cortical distance, gathered across multiple parcellation scales. Alongside, we share large-scale gradients estimated from each modality and parcellation scale. Our dataset will facilitate future research examining the coupling between brain microstructure, connectivity, and function. MICA-MICs is available on the Canadian Open Neuroscience Platform data portal ( https://portal.conp.ca ) and the Open Science Framework ( https://osf.io/j532r/ ).

A whole-brain 3D myeloarchitectonic atlas: Mapping the Vogt-Vogt legacy to the cortical surface

Building precise and detailed parcellations of anatomically and functionally distinct brain areas has been a major focus in Neuroscience. Pioneer anatomists parcellated the cortical manifold based on extensive histological studies of post-mortem brain, harnessing local variations in cortical cyto- and myeloarchitecture to define areal boundaries. Compared to the cytoarchitectonic field, where multiple neuroimaging studies have recently translated this old legacy data into useful analytical resources, myeloarchitectonics, which parcellate the cortex based on the organization of myelinated fibers, has received less attention. Here, we present the neocortical surface-based myeloarchitectonic atlas based on the histology-derived maps of the Vogt-Vogt school and its 2D translation by Nieuwenhuys. In addition to a myeloarchitectonic parcellation, our package includes intracortical laminar profiles of myelin content based on Vogt-Vogt-Hopf original publications. Histology-derived myelin density mapped on our atlas demonstrated a close overlap with in vivo quantitative MRI markers for myelin and relates to cytoarchitectural features. Complementing the existing battery of approaches for digital cartography, the whole-brain myeloarchitectonic atlas offers an opportunity to validate imaging surrogate markers of myelin in both health and disease.

Relationship of neurite architecture to brain activity during task-based fMRI

Functional MRI (fMRI) has been widely used to examine changes in neuronal activity during cognitive tasks. Commonly used measures of gray matter macrostructure (e.g., cortical thickness, surface area, volume) do not consistently appear to serve as structural correlates of brain function. In contrast, gray matter microstructure, measured using neurite orientation dispersion and density imaging (NODDI), enables the estimation of indices of neurite density (neurite density index; NDI) and organization (orientation dispersion index; ODI) in gray matter. Our study explored the relationship among neurite architecture, BOLD (blood-oxygen-level-dependent) fMRI, and cognition, using a large sample (n = 750) of young adults of the human connectome project (HCP) and two tasks that index more cortical (working memory) and more subcortical (emotion processing) targeting of brain functions. Using NODDI, fMRI, structural MRI and task performance data, hierarchical regression analyses revealed that higher working memory- and emotion processing-evoked BOLD activity was related to lower ODI in the right DLPFC, and lower ODI and NDI values in the right and left amygdala, respectively. Common measures of brain macrostructure (i.e., DLPFC thickness/surface area and amygdala volume) did not explain any additional variance (beyond neurite architecture) in BOLD activity. A moderating effect of neurite architecture on the relationship between emotion processing task-evoked BOLD response and performance was observed. Our findings provide evidence that neuro-/social-affective cognition-related BOLD activity is partially driven by the local neurite organization and density with direct impact on emotion processing. In vivo gray matter microstructure represents a new target of investigation providing strong potential for clinical translation.

Parental socioeconomic status is linked to cortical microstructure and language abilities in children and adolescents

Gradients in parental socioeconomic status (SES) are closely linked to important life outcomes in children and adolescents, such as cognitive abilities, school achievement, and mental health. Parental SES may also influence brain development, with several magnetic resonance imaging (MRI) studies reporting associations with youth brain morphometry. However, MRI signal intensity metrics have not been assessed, but could offer a microstructural correlate, thereby increasing our understanding of SES influences on neurobiology. We computed a parental SES score from family income, parental education and parental occupation, and assessed relations with cortical microstructure as measured by T1w/T2w ratio (n = 504, age = 3-21 years). We found negative age-stabile relations between parental SES and T1w/T2w ratio, indicating that youths from lower SES families have higher ratio in widespread frontal, temporal, medial parietal and occipital regions, possibly indicating a more developed cortex. Effect sizes were small, but larger than for conventional morphometric properties i.e. cortical surface area and thickness, which were not significantly associated with parental SES. Youths from lower SES families had poorer language related abilities, but microstructural differences did not mediate these relations. T1w/T2w ratio appears to be a sensitive imaging marker for further exploring the association between parental SES and child brain development.

Decoding transcriptional signatures of the association between free water and macroscale organizations in healthy adolescents

We leveraged a novel index of diffusion MRI to investigate the relationships among cortical free water, macro-organizations and gene expression in healthy adults. Few research has been conducted to investigate the role of free water in the healthy adults due to it can easily be affected also by aging diseases. High quality data of 350 subjects from Human Connectome Project were used in our study. Cortical free water was estimated by using a bi-tensor model. The free water was high in the limbic, insular and somatosensory cortex, while being lower in motor and association cortex. The negative correlation between the free water and cortical thickness has been consistently identified in almost all the cortical regions. Negative correlation between the cortical free water and structural covariance (rho=-0.38, p_spin=0.005) revealed the free water was sensitive to cortical heterogeneity. Using human gene expression dataset, we found the gene expression pattern of the relationship between the free water and cortical thickness spatially coupled with primary gradient of structural covariance network (rho=0.40, p_spin=0.004). Our findings indicated the free water was sensitive to the cortical cellular status. The relationship between free water and macroscale organization also reflected hierarchal structures of cerebral cortex.

Apolipoprotein E ε4 Mediates Myelin Breakdown by Targeting Oligodendrocytes in Sporadic Alzheimer Disease

White matter degradation in the frontal lobe is one of the earliest detectable changes in aging and Alzheimer disease. The ε4 allele of apolipoprotein E (APOE4) is strongly associated with such myelin pathology but the underlying cellular mechanisms remain obscure. We hypothesized that, as a lipid transporter, APOE4 directly triggers pathology in the cholesterol-rich myelin sheath independent of AD pathology. To test this, we performed immunohistochemistry on brain tissues from healthy controls, sporadic, and familial Alzheimer disease subjects. While myelin basic protein expression was largely unchanged, in frontal cortex the number of oligodendrocytes (OLs) was significantly reduced in APOE4 brains independent of their Braak stage or NIA-RI criteria. This high vulnerability of OLs was confirmed in humanized APOE3 or APOE4 transgenic mice. A gradual decline of OL numbers was found in the aging brain without associated neuronal loss. Importantly, the application of lipidated human APOE4, but not APOE3, proteins significantly reduced the formation of myelinating OL in primary cell culture derived from Apoe-knockout mice, especially in cholesterol-depleted conditions. Our findings suggest that the disruption of myelination in APOE4 carriers may represent a direct OL pathology, rather than an indirect consequence of amyloid plaque formation or neuronal loss.

Linking Plasma Amyloid Beta and Neurofilament Light Chain to Intracortical Myelin Content in Cognitively Normal Older Adults

Evidence suggests that lightly myelinated cortical regions are vulnerable to aging and Alzheimer’s disease (AD). However, it remains unknown whether plasma markers of amyloid and neurodegeneration are related to deficits in intracortical myelin content, and whether this relationship, in turn, is associated with altered patterns of resting-state functional connectivity (rs-FC). To shed light into these questions, plasma levels of amyloid-β fragment 1–42 (Aβ_1–42) and neurofilament light chain (NfL) were measured using ultra-sensitive single-molecule array (Simoa) assays, and the intracortical myelin content was estimated with the ratio T1-weigthed/T2-weighted (T1w/T2w) in 133 cognitively normal older adults. We assessed: (i) whether plasma Aβ_1–42 and/or NfL levels were associated with intracortical myelin content at different cortical depths and (ii) whether cortical regions showing myelin reductions also exhibited altered rs-FC patterns. Surface-based multiple regression analyses revealed that lower plasma Aβ_1–42 and higher plasma NfL were associated with lower myelin content in temporo-parietal-occipital regions and the insular cortex, respectively. Whereas the association with Aβ_1–42 decreased with depth, the NfL-myelin relationship was most evident in the innermost layer. Older individuals with higher plasma NfL levels also exhibited altered rs-FC between the insula and medial orbitofrontal cortex. Together, these findings establish a link between plasma markers of amyloid/neurodegeneration and intracortical myelin content in cognitively normal older adults, and support the role of plasma NfL in boosting aberrant FC patterns of the insular cortex, a central brain hub highly vulnerable to aging and neurodegeneration.

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.

Anatomy of the temporal lobe: From macro to micro

The temporal cortex encompasses a large number of different areas ranging from the six-layered isocortex to the allocortex. The areas support auditory, visual, and language processing, as well as emotions and memory. The primary auditory cortex is found at the Heschl gyri, which develop early in ontogeny with the Sylvian fissure, a deep and characteristic fissure that separates the temporal lobe from the parietal and frontal lobes. Gyri and sulci as well as brain areas vary between brains and between hemispheres, partly linked to the functional organization of language and lateralization. Interindividual variability in anatomy makes a direct comparison between different brains in structure-functional analysis often challenging, but can be addressed by applying cytoarchitectonic probability maps of the Julich-Brain atlas. We review the macroanatomy of the temporal lobe, its variability and asymmetry at the macro- and the microlevel, discuss the relationship to brain areas and their microstructure, and emphasize the advantage of a multimodal approach to address temporal lobe organization. We review recent data on combined cytoarchitectonic and molecular architectonic studies of temporal areas, and provide links to their function.

Prefrontal connectomics: from anatomy to human imaging

The fundamental importance of prefrontal cortical connectivity to information processing and, therefore, disorders of cognition, emotion, and behavior has been recognized for decades. Anatomic tracing studies in animals have formed the basis for delineating the direct monosynaptic connectivity, from cells of origin, through axon trajectories, to synaptic terminals. Advances in neuroimaging combined with network science have taken the lead in developing complex wiring diagrams or connectomes of the human brain. A key question is how well these magnetic resonance imaging (MRI)-derived networks and hubs reflect the anatomic "hard wiring" first proposed to underlie the distribution of information for large-scale network interactions. In this review, we address this challenge by focusing on what is known about monosynaptic prefrontal cortical connections in non-human primates and how this compares to MRI-derived measurements of network organization in humans. First, we outline the anatomic cortical connections and pathways for each prefrontal cortex (PFC) region. We then review the available MRI-based techniques for indirectly measuring structural and functional connectivity, and introduce graph theoretical methods for analysis of hubs, modules, and topologically integrative features of the connectome. Finally, we bring these two approaches together, using specific examples, to demonstrate how monosynaptic connections, demonstrated by tract-tracing studies, can directly inform understanding of the composition of PFC nodes and hubs, and the edges or pathways that connect PFC to cortical and subcortical areas.

Using high spatial resolution fMRI to understand representation in the auditory network

Following rapid methodological advances, ultra-high field (UHF) functional and anatomical magnetic resonance imaging (MRI) has been repeatedly and successfully used for the investigation of the human auditory system in recent years. Here, we review this work and argue that UHF MRI is uniquely suited to shed light on how sounds are represented throughout the network of auditory brain regions. That is, the provided gain in spatial resolution at UHF can be used to study the functional role of the small subcortical auditory processing stages and details of cortical processing. Further, by combining high spatial resolution with the versatility of MRI contrasts, UHF MRI has the potential to localize the primary auditory cortex in individual hemispheres. This is a prerequisite to study how sound representation in higher-level auditory cortex evolves from that in early (primary) auditory cortex. Finally, the access to independent signals across auditory cortical depths, as afforded by UHF, may reveal the computations that underlie the emergence of an abstract, categorical sound representation based on low-level acoustic feature processing. Efforts on these research topics are underway. Here we discuss promises as well as challenges that come with studying these research questions using UHF MRI, and provide a future outlook.

Recent developments in representations of the connectome

Research into the human connectome (i.e., all connections in the human brain) with the use of resting state functional MRI has rapidly increased in popularity in recent years, especially with the growing availability of large-scale neuroimaging datasets. The goal of this review article is to describe innovations in functional connectome representations that have come about in the past 8 years, since the 2013 NeuroImage special issue on 'Mapping the Connectome'. In the period, research has shifted from group-level brain parcellations towards the characterization of the individualized connectome and of relationships between individual connectomic differences and behavioral/clinical variation. Achieving subject-specific accuracy in parcel boundaries while retaining cross-subject correspondence is challenging, and a variety of different approaches are being developed to meet this challenge, including improved alignment, improved noise reduction, and robust group-to-subject mapping approaches. Beyond the interest in the individualized connectome, new representations of the data are being studied to complement the traditional parcellated connectome representation (i.e., pairwise connections between distinct brain regions), such as methods that capture overlapping and smoothly varying patterns of connectivity ('gradients'). These different connectome representations offer complimentary insights into the inherent functional organization of the brain, but challenges for functional connectome research remain. Interpretability will be improved by future research towards gaining insights into the neural mechanisms underlying connectome observations obtained from functional MRI. Validation studies comparing different connectome representations are also needed to build consensus and confidence to proceed with clinical trials that may produce meaningful clinical translation of connectome insights.

Magnetic resonance fingerprinting residual signals can disassociate human grey matter regions

The importance of accurate structural discrimination of the human grey matter regions has motivated the development of observer-independent reproducible methods that account for inter-individual architectonic variations. We introduce a non-invasive statistical residual analysis framework, employing unique tissue-specific magnetic resonance fingerprinting (MRF) signals after adjusting for the effect of T₁ and T₂^* MR relaxometry parameters (here termed MRF residuals). A 7 T Siemens MR scanner was used to acquire MRF signals, quantitative transmit magnetic field (B₁⁺) maps and T₁-weighted anatomical images of eleven cortical areas (5L, 5M, 5Ci, 7A, 7P, 7PC, hIP3, BA2, BA4a, BA4p and BA6) from six female participants. MRF residual signal for each voxel was calculated as the difference between the actual and best matching MRF signal evolutions from a precomputed MRF dictionary covering a range of T₁, T₂^* and B₁⁺ values. To compare MRF residuals between regions of interest, normalised autocorrelation was used as a shape-based statistical signal characterisation method and the Euclidean distance between autocorrelation profiles of residuals was used to measure the interareal dissimilarity. In the eleven cortical areas in both cerebral hemispheres of six participants, the proposed MRF residual analysis consistently showed interareal dissimilarity profiles that concorded with histological studies, indicating that MRF residuals potentially contain tissue microstructural information. MRF residual signals provide additional area-specific information that is complementary to the MR relaxometry-based (T₁, T₂^*) information used previously for distinguishing microstructural differences between human cerebral cortex regions in vivo. The proposed approach led to more accurate identification of structural variations across cortical areas of interest.

Quantitative Macromolecular Proton Fraction Mapping Reveals Altered Cortical Myelin Profile in Schizophrenia Spectrum Disorders

Myelin abnormalities have been reported in schizophrenia spectrum disorders (SSD) in white matter. However, in vivo examinations of cortical myeloarchitecture in SSD, especially those using quantitative measures, are limited. Here, we employed macromolecular proton fraction (MPF) obtained from quantitative magnetization transfer imaging to characterize intracortical myelin organization in 30 SSD patients versus 34 healthy control (HC) participants. We constructed cortical myelin profiles by extracting MPF values at various cortical depths and quantified their shape using a nonlinearity index (NLI). To delineate the association of illness duration with myelin changes, SSD patients were further divided into 3 duration groups. Between-group comparisons revealed reduced NLI in the SSD group with the longest illness duration (>5.5 years) compared with HC predominantly in bilateral prefrontal areas. Within the SSD group, cortical NLI decreased with disease duration and was positively associated with a measure of spatial working memory capacity as well as with cortical thickness (CT). Layer-specific analyses suggested that NLI decreases in the long-duration SSD group may arise in part from significantly increased MPF values in the midcortical layers. The current study reveals cortical myelin profile changes in SSD with illness progression, which may reflect an abnormal compensatory mechanism of the disorder.

New insights into the dynamic development of the cerebral cortex in childhood and adolescence: Integrating macro- and microstructural MRI findings

Through dynamic transactional processes between genetic and environmental factors, childhood and adolescence involve reorganization and optimization of the cerebral cortex. The cortex and its development plays a crucial role for prototypical human cognitive abilities. At the same time, many common mental disorders appear during these critical phases of neurodevelopment. Magnetic resonance imaging (MRI) can indirectly capture several multifaceted changes of cortical macro- and microstructure, of high relevance to further our understanding of the neural foundation of cognition and mental health. Great progress has been made recently in mapping the typical development of cortical morphology. Moreover, newer less explored MRI signal intensity and specialized quantitative T2 measures have been applied to assess microstructural cortical development. We review recent findings of typical postnatal macro- and microstructural development of the cerebral cortex from early childhood to young adulthood. We cover studies of cortical volume, thickness, area, gyrification, T1-weighted (T1w) tissue contrasts such a grey/white matter contrast, T1w/T2w ratio, magnetization transfer and myelin water fraction. Finally, we integrate imaging studies with cortical gene expression findings to further our understanding of the underlying neurobiology of the developmental changes, bridging the gap between ex vivo histological- and in vivo MRI studies.

Revisiting the effects of transcranial direct current stimulation on pattern-reversal visual evoked potentials

Despite increasing growth of interest in transcranial direct current stimulation (tDCS), its underlying mechanisms are still unclear. With many claims based on the anodal-excitation and cathodal-inhibition dichotomy originally observed in the motor cortex, surprisingly few studies have examined these fundamental polarity-specific effects beyond the motor cortex. The after-effects of tDCS on the visual cortex are of particular interest because of their potential application to vision restoration and migraine treatment. Yet the limited studies revealed conflicting results. Here we investigated whether polarity-specific tDCS effects exist in the visual cortex. In a counterbalanced within-subject crossover design, 20 healthy subjects each completed three sessions of anodal, cathodal and sham tDCS (2 mA for 20 min) applied over the visual cortex. Pattern-reversal visual evoked potentials (VEP) and their habituation slopes were measured at five time-points immediately before, after and every 15 min following the end of tDCS. Compared to sham, we found no significant tDCS induced after-effects on VEP amplitudes or habituation slopes, supported by strong evidence from Bayesian statistics. Neither were there any after-effects of tDCS on EEG power of the frequency of stimulus presentation, theta or alpha band. In conclusion, our results challenge previous findings of robust polarity-dependent after-effects of tDCS over the visual cortex.

Dissociable Auditory Cortico-Cerebellar Pathways in the Human Brain Estimated by Intrinsic Functional Connectivity

The cerebellum, a structure historically associated with motor control, has more recently been implicated in several higher-order auditory-cognitive functions. However, the exact functional pathways that mediate cerebellar influences on auditory cortex (AC) remain unclear. Here, we sought to identify auditory cortico-cerebellar pathways based on intrinsic functional connectivity magnetic resonance imaging. In contrast to previous connectivity studies that principally consider the AC as a single functionally homogenous unit, we mapped the cerebellar connectivity across different parts of the AC. Our results reveal that auditory subareas demonstrating different levels of interindividual functional variability are functionally coupled with distinct cerebellar regions. Moreover, auditory and sensorimotor areas show divergent cortico-cerebellar connectivity patterns, although sensorimotor areas proximal to the AC are often functionally grouped with the AC in previous connectivity-based network analyses. Lastly, we found that the AC can be functionally segmented into highly similar subareas based on either cortico-cerebellar or cortico-cortical functional connectivity, suggesting the existence of multiple parallel auditory cortico-cerebellar circuits that involve different subareas of the AC. Overall, the present study revealed multiple auditory cortico-cerebellar pathways and provided a fine-grained map of AC subareas, indicative of the critical role of the cerebellum in auditory processing and multisensory integration.

The developing relations between networks of cortical myelin and neurophysiological connectivity

Recent work identified that patterns of distributed brain regions sharing similar myeloarchitecture are related to underlying functional connectivity, demonstrating cortical myelin's plasticity to changes in functional demand. However, the changing relations between functional and structural architecture throughout child and adulthood are poorly understood. We show that structural covariance connectivity (T1-weighted/T2-weighted ratio) and functional connectivity (magnetoencephalography) exhibit nonlinear developmental changes. We then show significant relations between structural and functional connectivity, which have shared and distinct characteristics dependent on the neural oscillatory frequency. Increases in structure-function coupling are visible during the protracted myelination observed throughout childhood and adolescence and are followed by decreases near the onset of adulthood. Our work lays the foundation for understanding the mechanisms by which myeloarchitecture supports brain function, enabling future investigations into how clinical populations may deviate from normative patterns.

Layer-specific parameters of intracortical microstimulation of the somatosensory cortex

Objective. Intracortical microstimulation of the primary somatosensory cortex (S1) has shown great progress in restoring touch sensations to patients with paralysis. Stimulation parameters such as amplitude, phase duration, and frequency can influence the quality of the evoked percept as well as the amount of charge necessary to elicit a response. Previous studies in V1 and auditory cortices have shown that the behavioral responses to stimulation amplitude and phase duration change across cortical depth. However, this depth-dependent response has yet to be investigated in S1. Similarly, to our knowledge, the response to microstimulation frequency across cortical depth remains unexplored.Approach. To assess these questions, we implanted rats in S1 with a microelectrode with electrode-sites spanning all layers of the cortex. A conditioned avoidance behavioral paradigm was used to measure detection thresholds and responses to phase duration and frequency across cortical depth.Main results. Analogous to other cortical areas, the sensitivity to charge and strength-duration chronaxies in S1 varied across cortical layers. Likewise, the sensitivity to microstimulation frequency was layer dependent.Significance. These findings suggest that cortical depth can play an important role in the fine-tuning of stimulation parameters and in the design of intracortical neuroprostheses for clinical applications.

The Myelin Content of the Human Precentral Hand Knob Reflects Interindividual Differences in Manual Motor Control at the Physiological and Behavioral Level

The primary motor cortex hand area (M1_HAND) and adjacent dorsal premotor cortex (PMd) form the so-called motor hand knob in the precentral gyrus. M1_HAND and PMd are critical for dexterous hand use and are densely interconnected via corticocortical axons, lacking a sharp demarcating border. In 24 young right-handed volunteers, we performed multimodal mapping to delineate the relationship between structure and function in the right motor hand knob. Quantitative structural magnetic resonance imaging (MRI) at 3 tesla yielded regional R1 maps as a proxy of cortical myelin content. Participants also underwent functional MRI (fMRI). We mapped task-related activation and temporal precision, while they performed a visuomotor synchronization task requiring visually cued abduction movements with the left index or little finger. We also performed sulcus-aligned transcranial magnetic stimulation of the motor hand knob to localize the optimal site (hotspot) for evoking a motor evoked potential (MEP) in two intrinsic hand muscles. Individual motor hotspot locations varied along the rostrocaudal axis. The more rostral the motor hotspot location in the precentral crown, the longer were corticomotor MEP latencies. “Hotspot rostrality” was associated with the regional myelin content in the precentral hand knob. Cortical myelin content also correlated positively with task-related activation of the precentral crown and temporal precision during the visuomotor synchronization task. Together, our results suggest a link among cortical myelination, the spatial cortical representation, and temporal precision of finger movements. We hypothesize that the myelination of cortical axons facilitates neuronal integration in PMd and M1_HAND and, hereby, promotes the precise timing of movements.

SIGNIFICANCE STATEMENT Here we used magnetic resonance imaging and transcranial magnetic stimulation of the precentral motor hand knob to test for a link among cortical myelin content, functional corticomotor representations, and manual motor control. A higher myelin content of the precentral motor hand knob was associated with more rostral corticomotor presentations, with stronger task-related activation and a higher precision of movement timing during a visuomotor synchronization task. We propose that a high precentral myelin content enables fast and precise neuronal integration in M1 (primary motor cortex) and dorsal premotor cortex, resulting in higher temporal precision during dexterous hand use. Our results identify the degree of myelination as an important structural feature of the neocortex that is tightly linked to the function and behavior supported by the cortical area.