The BigBrainWarp toolbox for integration of BigBrain 3D histology with multimodal neuroimaging

A sensorimotor-association axis of thalamocortical connection development

Human cortical development follows a sensorimotor-to-association sequence during childhood and adolescence1-6. The brain's capacity to enact this sequence over decades indicates that it relies on intrinsic mechanisms to regulate inter-regional differences in the timing of cortical maturation, yet regulators of human developmental chronology are not well understood. Given evidence from animal models that thalamic axons modulate windows of cortical plasticity7-12, here we evaluate the overarching hypothesis that structural connections between the thalamus and cortex help to coordinate cortical maturational heterochronicity during youth. We first introduce, cortically annotate, and anatomically validate a new atlas of human thalamocortical connections using diffusion tractography. By applying this atlas to three independent youth datasets (ages 8-23 years; total N = 2,676), we reproducibly demonstrate that thalamocortical connections develop along a maturational gradient that aligns with the cortex's sensorimotor-association axis. Associative cortical regions with thalamic connections that take longest to mature exhibit protracted expression of neurochemical, structural, and functional markers indicative of higher circuit plasticity as well as heightened environmental sensitivity. This work highlights a central role for the thalamus in the orchestration of hierarchically organized and environmentally sensitive windows of cortical developmental malleability.

Functional connectivity development along the sensorimotor-association axis enhances the cortical hierarchy

Human cortical maturation has been posited to be organized along the sensorimotor-association axis, a hierarchical axis of brain organization that spans from unimodal sensorimotor cortices to transmodal association cortices. Here, we investigate the hypothesis that the development of functional connectivity during childhood through adolescence conforms to the cortical hierarchy defined by the sensorimotor-association axis. We tested this pre-registered hypothesis in four large-scale, independent datasets (total n = 3355; ages 5-23 years): the Philadelphia Neurodevelopmental Cohort (n = 1207), Nathan Kline Institute-Rockland Sample (n = 397), Human Connectome Project: Development (n = 625), and Healthy Brain Network (n = 1126). Across datasets, the development of functional connectivity systematically varied along the sensorimotor-association axis. Connectivity in sensorimotor regions increased, whereas connectivity in association cortices declined, refining and reinforcing the cortical hierarchy. These consistent and generalizable results establish that the sensorimotor-association axis of cortical organization encodes the dominant pattern of functional connectivity development.

Molecular signatures of cortical expansion in the human fetal brain

The third trimester of human gestation is characterised by rapid increases in brain volume and cortical surface area. A growing catalogue of cells in the prenatal brain has revealed remarkable molecular diversity across cortical areas.1,2 Despite this, little is known about how this translates into the patterns of differential cortical expansion observed in humans during the latter stages of gestation. Here we present a new resource, μBrain, to facilitate knowledge translation between molecular and anatomical descriptions of the prenatal developing brain. Built using generative artificial intelligence, μBrain is a three-dimensional cellular-resolution digital atlas combining publicly-available serial sections of the postmortem human brain at 21 weeks gestation3 with bulk tissue microarray data, sampled across 29 cortical regions and 5 transient tissue zones.4 Using μBrain, we evaluate the molecular signatures of preferentially-expanded cortical regions during human gestation, quantified in utero using magnetic resonance imaging (MRI). We find that differences in the rates of expansion across cortical areas during gestation respect anatomical and evolutionary boundaries between cortical types5 and are founded upon extended periods of upper-layer cortical neuron migration that continue beyond mid-gestation. We identify a set of genes that are upregulated from mid-gestation and highly expressed in rapidly expanding neocortex, which are implicated in genetic disorders with cognitive sequelae. Our findings demonstrate a spatial coupling between areal differences in the timing of neurogenesis and rates of expansion across the neocortical sheet during the prenatal epoch. The μBrain atlas is available from: https://garedaba.github.io/micro-brain/ and provides a new tool to comprehensively map early brain development across domains, model systems and resolution scales.

How brain structure-function decoupling supports individual cognition and its molecular mechanism

Functional signals emerge from the structural network, supporting multiple cognitive processes through underlying molecular mechanism. The link between human brain structure and function is region-specific and hierarchical across the neocortex. However, the relationship between hierarchical structure-function decoupling and the manifestation of individual behavior and cognition, along with the significance of the functional systems involved, and the specific molecular mechanism underlying structure-function decoupling remain incompletely characterized. Here, we used the structural-decoupling index (SDI) to quantify the dependency of functional signals on the structural connectome using a significantly larger cohort of healthy subjects. Canonical correlation analysis (CCA) was utilized to assess the general multivariate correlation pattern between region-specific SDIs across the whole brain and multiple cognitive traits. Then, we predicted five composite cognitive scores resulting from multivariate analysis using SDIs in primary networks, association networks, and all networks, respectively. Finally, we explored the molecular mechanism related to SDI by investigating its genetic factors and relationship with neurotransmitter receptors/transporters. We demonstrated that structure-function decoupling is hierarchical across the neocortex, spanning from primary networks to association networks. We revealed better performance in cognition prediction is achieved by using high-level hierarchical SDIs, with varying significance of different brain regions in predicting cognitive processes. We found that the SDIs were associated with the gene expression level of several receptor-related terms, and we also found the spatial distributions of four receptors/transporters significantly correlated with SDIs, namely D2, NET, MOR, and mGluR5, which play an important role in the flexibility of neuronal function. Collectively, our findings corroborate the association between hierarchical macroscale structure-function decoupling and individual cognition and provide implications for comprehending the molecular mechanism of structure-function decoupling. PRACTITIONER POINTS: Structure-function decoupling is hierarchical across the neocortex, spanning from primary networks to association networks. High-level hierarchical structure-function decoupling contributes much more than low-level decoupling to individual cognition. Structure-function decoupling could be regulated by genes associated with pivotal receptors that are crucial for neuronal function flexibility.

An energy costly architecture of neuromodulators for human brain evolution and cognition

In comparison to other species, the human brain exhibits one of the highest energy demands relative to body metabolism. It remains unclear whether this heightened energy demand uniformly supports an enlarged brain or if specific signaling mechanisms necessitate greater energy. We hypothesized that the regional distribution of energy demands will reveal signaling strategies that have contributed to human cognitive development. We measured the energy distribution within the brain functional connectome using multimodal brain imaging and found that signaling pathways in evolutionarily expanded regions have up to 67% higher energetic costs than those in sensory-motor regions. Additionally, histology, transcriptomic data, and molecular imaging independently reveal an up-regulation of signaling at G-protein-coupled receptors in energy-demanding regions. Our findings indicate that neuromodulator activity is predominantly involved in cognitive functions, such as reading or memory processing. This study suggests that an up-regulation of neuromodulator activity, alongside increased brain size, is a crucial aspect of human brain evolution.

Whole-brain Optical Imaging: A Powerful Tool for Precise Brain Mapping at the Mesoscopic Level

The mammalian brain is a highly complex network that consists of millions to billions of densely-interconnected neurons. Precise dissection of neural circuits at the mesoscopic level can provide important structural information for understanding the brain. Optical approaches can achieve submicron lateral resolution and achieve "optical sectioning" by a variety of means, which has the natural advantage of allowing the observation of neural circuits at the mesoscopic level. Automated whole-brain optical imaging methods based on tissue clearing or histological sectioning surpass the limitation of optical imaging depth in biological tissues and can provide delicate structural information in a large volume of tissues. Combined with various fluorescent labeling techniques, whole-brain optical imaging methods have shown great potential in the brain-wide quantitative profiling of cells, circuits, and blood vessels. In this review, we summarize the principles and implementations of various whole-brain optical imaging methods and provide some concepts regarding their future development.

The regional variation of laminar thickness in the human isocortex is related to cortical hierarchy and interregional connectivity

The human isocortex consists of tangentially organized layers with unique cytoarchitectural properties. These layers show spatial variations in thickness and cytoarchitecture across the neocortex, which is thought to support function through enabling targeted corticocortical connections. Here, leveraging maps of the 6 cortical layers based on 3D human brain histology, we aimed to quantitatively characterize the systematic covariation of laminar structure in the cortex and its functional consequences. After correcting for the effect of cortical curvature, we identified a spatial pattern of changes in laminar thickness covariance from lateral frontal to posterior occipital regions, which differentiated the dominance of infra- versus supragranular layer thickness. Corresponding to the laminar regularities of cortical connections along cortical hierarchy, the infragranular-dominant pattern of laminar thickness was associated with higher hierarchical positions of regions, mapped based on resting-state effective connectivity in humans and tract-tracing of structural connections in macaques. Moreover, we show that regions with similar laminar thickness patterns have a higher likelihood of structural connections and strength of functional connections. In sum, here we characterize the organization of laminar thickness in the human isocortex and its association with cortico-cortical connectivity, illustrating how laminar organization may provide a foundational principle of cortical function.

The biological role of local and global fMRI BOLD signal variability in human brain organization

Variability drives the organization and behavior of complex systems, including the human brain. Understanding the variability of brain signals is thus necessary to broaden our window into brain function and behavior. Few empirical investigations of macroscale brain signal variability have yet been undertaken, given the difficulty in separating biological sources of variance from artefactual noise. Here, we characterize the temporal variability of the most predominant macroscale brain signal, the fMRI BOLD signal, and systematically investigate its statistical, topographical and neurobiological properties. We contrast fMRI acquisition protocols, and integrate across histology, microstructure, transcriptomics, neurotransmitter receptor and metabolic data, fMRI static connectivity, and empirical and simulated magnetoencephalography data. We show that BOLD signal variability represents a spatially heterogeneous, central property of multi-scale multi-modal brain organization, distinct from noise. Our work establishes the biological relevance of BOLD signal variability and provides a lens on brain stochasticity across spatial and temporal scales.

Neurophysiological signatures of cortical micro-architecture

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6800 time-series features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is co-localized with multiple micro-architectural features, including gene expression gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

Cortical microstructural gradients capture memory network reorganization in temporal lobe epilepsy

Temporal lobe epilepsy (TLE), one of the most common pharmaco-resistant epilepsies, is associated with pathology of paralimbic brain regions, particularly in the mesiotemporal lobe. Cognitive dysfunction in TLE is frequent, and particularly affects episodic memory. Crucially, these difficulties challenge the quality of life of patients, sometimes more than seizures, underscoring the need to assess neural processes of cognitive dysfunction in TLE to improve patient management. Our work harnessed a novel conceptual and analytical approach to assess spatial gradients of microstructural differentiation between cortical areas based on high-resolution MRI analysis. Gradients track region-to-region variations in intracortical lamination and myeloarchitecture, serving as a system-level measure of structural and functional reorganization. Comparing cortex-wide microstructural gradients between 21 patients and 35 healthy controls, we observed a reorganization of this gradient in TLE driven by reduced microstructural differentiation between paralimbic cortices and the remaining cortex with marked abnormalities in ipsilateral temporopolar and dorsolateral prefrontal regions. Findings were replicated in an independent cohort. Using an independent post-mortem dataset, we observed that in vivo findings reflected topographical variations in cortical cytoarchitecture. We indeed found that macroscale changes in microstructural differentiation in TLE reflected increased similarity of paralimbic and primary sensory/motor regions. Disease-related transcriptomics could furthermore show specificity of our findings to TLE over other common epilepsy syndromes. Finally, microstructural dedifferentiation was associated with cognitive network reorganization seen during an episodic memory functional MRI paradigm and correlated with interindividual differences in task accuracy. Collectively, our findings showing a pattern of reduced microarchitectural differentiation between paralimbic regions and the remaining cortex provide a structurally-grounded explanation for large-scale functional network reorganization and cognitive dysfunction characteristic of TLE.

Integrating multimodal and multiscale connectivity blueprints of the human cerebral cortex in health and disease

The brain is composed of disparate neural populations that communicate and interact with one another. Although fiber bundles, similarities in molecular architecture, and synchronized neural activity all reflect how brain regions potentially interact with one another, a comprehensive study of how all these interregional relationships jointly reflect brain structure and function remains missing. Here, we systematically integrate 7 multimodal, multiscale types of interregional similarity ("connectivity modes") derived from gene expression, neurotransmitter receptor density, cellular morphology, glucose metabolism, haemodynamic activity, and electrophysiology in humans. We first show that for all connectivity modes, feature similarity decreases with distance and increases when regions are structurally connected. Next, we show that connectivity modes exhibit unique and diverse connection patterns, hub profiles, spatial gradients, and modular organization. Throughout, we observe a consistent primacy of molecular connectivity modes-namely correlated gene expression and receptor similarity-that map onto multiple phenomena, including the rich club and patterns of abnormal cortical thickness across 13 neurological, psychiatric, and neurodevelopmental disorders. Finally, to construct a single multimodal wiring map of the human cortex, we fuse all 7 connectivity modes and show that the fused network maps onto major organizational features of the cortex including structural connectivity, intrinsic functional networks, and cytoarchitectonic classes. Altogether, this work contributes to the integrative study of interregional relationships in the human cerebral cortex.

Spatially heterogeneous structure-function coupling in haemodynamic and electromagnetic brain networks

The relationship between structural and functional connectivity in the brain is a key question in connectomics. Here we quantify patterns of structure-function coupling across the neocortex, by comparing structural connectivity estimated using diffusion MRI with functional connectivity estimated using both neurophysiological (MEG-based) and haemodynamic (fMRI-based) recordings. We find that structure-function coupling is heterogeneous across brain regions and frequency bands. The link between structural and functional connectivity is generally stronger in multiple MEG frequency bands compared to resting state fMRI. Structure-function coupling is greater in slower and intermediate frequency bands compared to faster frequency bands. We also find that structure-function coupling systematically follows the archetypal sensorimotor-association hierarchy, as well as patterns of laminar differentiation, peaking in granular layer IV. Finally, structure-function coupling is better explained using structure-informed inter-regional communication metrics than using structural connectivity alone. Collectively, these results place neurophysiological and haemodynamic structure-function relationships in a common frame of reference and provide a starting point for a multi-modal understanding of structure-function coupling in the brain.

Multimodal mapping of regional brain vulnerability to focal cortical dysplasia

Focal cortical dysplasia (FCD) type II is a highly epileptogenic developmental malformation and a common cause of surgically treated drug-resistant epilepsy. While clinical observations suggest frequent occurrence in the frontal lobe, mechanisms for such propensity remain unexplored. Here, we hypothesized that cortex-wide spatial associations of FCD distribution with cortical cytoarchitecture, gene expression and organizational axes may offer complementary insights into processes that predispose given cortical regions to harbour FCD. We mapped the cortex-wide MRI distribution of FCDs in 337 patients collected from 13 sites worldwide. We then determined its associations with (i) cytoarchitectural features using histological atlases by Von Economo and Koskinas and BigBrain; (ii) whole-brain gene expression and spatiotemporal dynamics from prenatal to adulthood stages using the Allen Human Brain Atlas and PsychENCODE BrainSpan; and (iii) macroscale developmental axes of cortical organization. FCD lesions were preferentially located in the prefrontal and fronto-limbic cortices typified by low neuron density, large soma and thick grey matter. Transcriptomic associations with FCD distribution uncovered a prenatal component related to neuroglial proliferation and differentiation, likely accounting for the dysplastic makeup, and a postnatal component related to synaptogenesis and circuit organization, possibly contributing to circuit-level hyperexcitability. FCD distribution showed a strong association with the anterior region of the antero-posterior axis derived from heritability analysis of interregional structural covariance of cortical thickness, but not with structural and functional hierarchical axes. Reliability of all results was confirmed through resampling techniques. Multimodal associations with cytoarchitecture, gene expression and axes of cortical organization indicate that prenatal neurogenesis and postnatal synaptogenesis may be key points of developmental vulnerability of the frontal lobe to FCD. Concordant with a causal role of atypical neuroglial proliferation and growth, our results indicate that FCD-vulnerable cortices display properties indicative of earlier termination of neurogenesis and initiation of cell growth. They also suggest a potential contribution of aberrant postnatal synaptogenesis and circuit development to FCD epileptogenicity.

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.

Intrinsic activity development unfolds along a sensorimotor-association cortical axis in youth

Animal studies of neurodevelopment have shown that recordings of intrinsic cortical activity evolve from synchronized and high amplitude to sparse and low amplitude as plasticity declines and the cortex matures. Leveraging resting-state functional MRI (fMRI) data from 1,033 youths (ages 8-23 years), we find that this stereotyped refinement of intrinsic activity occurs during human development and provides evidence for a cortical gradient of neurodevelopmental change. Declines in the amplitude of intrinsic fMRI activity were initiated heterochronously across regions and were coupled to the maturation of intracortical myelin, a developmental plasticity regulator. Spatiotemporal variability in regional developmental trajectories was organized along a hierarchical, sensorimotor-association cortical axis from ages 8 to 18. The sensorimotor-association axis furthermore captured variation in associations between youths' neighborhood environments and intrinsic fMRI activity; associations suggest that the effects of environmental disadvantage on the maturing brain diverge most across this axis during midadolescence. These results uncover a hierarchical neurodevelopmental axis and offer insight into the progression of cortical plasticity in humans.

Mapping Cytoarchitectonics and Receptor Architectonics to Understand Brain Function and Connectivity

This review focuses on cytoarchitectonics and receptor architectonics as biological correlates of function and connectivity. It introduces the 3-dimensional cytoarchitectonic probabilistic maps of cortical areas and nuclei of the Julich-Brain Atlas, available at EBRAINS, to study structure-function relationships. The maps are linked to the BigBrain as microanatomical reference model and template space. The siibra software tool suite enables programmatic access to the maps and to receptor architectonic data that are anchored to brain areas. Such cellular and molecular data are tools for studying magnetic resonance connectivity including modeling and simulation. At the end, we highlight perspectives of the Julich-Brain as well as methodological considerations. Thus, microstructural maps as part of a multimodal atlas help elucidate the biological correlates of large-scale networks and brain function with a high level of anatomical detail, which provides a basis to study brains of patients with psychiatric disorders.

BigBrain-MR: a new digital phantom with anatomically-realistic magnetic resonance properties at 100-µm resolution for magnetic resonance methods development

The benefits, opportunities and growing availability of ultra-high field magnetic resonance imaging (MRI) for humans have prompted an expansion in research and development efforts towards increasingly more advanced high-resolution imaging techniques. To maximize their effectiveness, these efforts need to be supported by powerful computational simulation platforms that can adequately reproduce the biophysical characteristics of MRI, with high spatial resolution. In this work, we have sought to address this need by developing a novel digital phantom with realistic anatomical detail up to 100-µm resolution, including multiple MRI properties that affect image generation. This phantom, termed BigBrain-MR, was generated from the publicly available BigBrain histological dataset and lower-resolution in-vivo 7T-MRI data, using a newly-developed image processing framework that allows mapping the general properties of the latter into the fine anatomical scale of the former. Overall, the mapping framework was found to be effective and robust, yielding a diverse range of realistic "in-vivo-like" MRI contrasts and maps at 100-µm resolution. BigBrain-MR was then tested in three imaging applications (motion effects and interpolation, super-resolution imaging, and parallel imaging reconstruction) to investigate its properties, value and validity as a simulation platform. The results consistently showed that BigBrain-MR can closely approximate the behavior of real in-vivo data, more realistically and with more extensive features than a more classic option such as the Shepp-Logan phantom. Its flexibility in simulating different contrast mechanisms and artifacts may also prove valuable for educational applications. BigBrain-MR is therefore deemed a favorable choice to support methodological development and demonstration in brain MRI, and has been made freely available to the community.

Neurophysiological signatures of cortical micro-architecture

Systematic spatial variation in micro-architecture is observed across the cortex. These micro-architectural gradients are reflected in neural activity, which can be captured by neurophysiological time-series. How spontaneous neurophysiological dynamics are organized across the cortex and how they arise from heterogeneous cortical micro-architecture remains unknown. Here we extensively profile regional neurophysiological dynamics across the human brain by estimating over 6 800 timeseries features from the resting state magnetoencephalography (MEG) signal. We then map regional time-series profiles to a comprehensive multi-modal, multi-scale atlas of cortical micro-architecture, including microstructure, metabolism, neurotransmitter receptors, cell types and laminar differentiation. We find that the dominant axis of neurophysiological dynamics reflects characteristics of power spectrum density and linear correlation structure of the signal, emphasizing the importance of conventional features of electromagnetic dynamics while identifying additional informative features that have traditionally received less attention. Moreover, spatial variation in neurophysiological dynamics is colocalized with multiple micro-architectural features, including genomic gradients, intracortical myelin, neurotransmitter receptors and transporters, and oxygen and glucose metabolism. Collectively, this work opens new avenues for studying the anatomical basis of neural activity.

Quantitative optical coherence microscopy of neuron morphology in human entorhinal cortex

Introduction

The size and shape of neurons are important features indicating aging and the pathology of neurodegenerative diseases. Despite the significant advances of optical microscopy, quantitative analysis of the neuronal features in the human brain remains largely incomplete. Traditional histology on thin slices bears tremendous distortions in three-dimensional reconstruction, the magnitude of which are often greater than the structure of interest. Recently development of tissue clearing techniques enable the whole brain to be analyzed in small animals; however, the application in the human remains challenging.

Methods

In this study, we present a label-free quantitative optical coherence microscopy (OCM) technique to obtain the morphological parameters of neurons in human entorhinal cortex (EC). OCM uses the intrinsic back-scattering property of tissue to identify individual neurons in 3D. The area, length, width, and orientation of individual neurons are quantified and compared between layer II and III in EC.

Results

The high-resolution mapping of neuron size, shape, and orientation shows significant differences between layer II and III neurons in EC. The results are validated by standard Nissl staining of the same samples.

Discussion

The quantitative OCM technique in our study offers a new solution to analyze variety of neurons and their organizations in the human brain, which opens new insights in advancing our understanding of neurodegenerative diseases.

Asymmetric signaling across the hierarchy of cytoarchitecture within the human connectome

Cortical variations in cytoarchitecture form a sensory-fugal axis that shapes regional profiles of extrinsic connectivity and is thought to guide signal propagation and integration across the cortical hierarchy. While neuroimaging work has shown that this axis constrains local properties of the human connectome, it remains unclear whether it also shapes the asymmetric signaling that arises from higher-order topology. Here, we used network control theory to examine the amount of energy required to propagate dynamics across the sensory-fugal axis. Our results revealed an asymmetry in this energy, indicating that bottom-up transitions were easier to complete compared to top-down. Supporting analyses demonstrated that asymmetries were underpinned by a connectome topology that is wired to support efficient bottom-up signaling. Lastly, we found that asymmetries correlated with differences in communicability and intrinsic neuronal time scales and lessened throughout youth. Our results show that cortical variation in cytoarchitecture may guide the formation of macroscopic connectome topology.

Micapipe: A pipeline for multimodal neuroimaging and connectome analysis

Multimodal magnetic resonance imaging (MRI) has accelerated human neuroscience by fostering the analysis of brain microstructure, geometry, function, and connectivity across multiple scales and in living brains. The richness and complexity of multimodal neuroimaging, however, demands processing methods to integrate information across modalities and to consolidate findings across different spatial scales. Here, we present micapipe, an open processing pipeline for multimodal MRI datasets. Based on BIDS-conform input data, micapipe can generate i) structural connectomes derived from diffusion tractography, ii) functional connectomes derived from resting-state signal correlations, iii) geodesic distance matrices that quantify cortico-cortical proximity, and iv) microstructural profile covariance matrices that assess inter-regional similarity in cortical myelin proxies. The above matrices can be automatically generated across established 18 cortical parcellations (100-1000 parcels), in addition to subcortical and cerebellar parcellations, allowing researchers to replicate findings easily across different spatial scales. Results are represented on three different surface spaces (native, conte69, fsaverage5), and outputs are BIDS-conform. Processed outputs can be quality controlled at the individual and group level. micapipe was tested on several datasets and is available at https://github.com/MICA-MNI/micapipe, documented at https://micapipe.readthedocs.io/, and containerized as a BIDS App http://bids-apps.neuroimaging.io/apps/. We hope that micapipe will foster robust and integrative studies of human brain microstructure, morphology, function, cand connectivity.

Multiscale neural gradients reflect transdiagnostic effects of major psychiatric conditions on cortical morphology

It is increasingly recognized that multiple psychiatric conditions are underpinned by shared neural pathways, affecting similar brain systems. Here, we carried out a multiscale neural contextualization of shared alterations of cortical morphology across six major psychiatric conditions (autism spectrum disorder, attention deficit/hyperactivity disorder, major depression disorder, obsessive-compulsive disorder, bipolar disorder, and schizophrenia). Our framework cross-referenced shared morphological anomalies with respect to cortical myeloarchitecture and cytoarchitecture, as well as connectome and neurotransmitter organization. Pooling disease-related effects on MRI-based cortical thickness measures across six ENIGMA working groups, including a total of 28,546 participants (12,876 patients and 15,670 controls), we identified a cortex-wide dimension of morphological changes that described a sensory-fugal pattern, with paralimbic regions showing the most consistent alterations across conditions. The shared disease dimension was closely related to cortical gradients of microstructure as well as neurotransmitter axes, specifically cortex-wide variations in serotonin and dopamine. Multiple sensitivity analyses confirmed robustness with respect to slight variations in analytical choices. Our findings embed shared effects of common psychiatric conditions on brain structure in multiple scales of brain organization, and may provide insights into neural mechanisms of transdiagnostic vulnerability.

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic-haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

Closing the mechanistic gap: the value of microarchitecture in understanding cognitive networks

Cognitive neuroscience aims to provide biologically relevant accounts of cognition. Contemporary research linking spatial patterns of neural activity to psychological constructs describes 'where' hypothesised functions occur, but not 'how' these regions contribute to cognition. Technological, empirical, and conceptual advances allow this mechanistic gap to be closed by embedding patterns of functional activity in macro- and microscale descriptions of brain organisation. Recent work on the default mode network (DMN) and the multiple demand network (MDN), for example, highlights a microarchitectural landscape that may explain how activity in these networks integrates varied information, thus providing an anatomical foundation that will help to explain how these networks contribute to many different cognitive states. This perspective highlights emerging insights into how microarchitecture can constrain network accounts of human cognition.

Brain and cognitive ageing: The present, and some predictions (…about the future)

Experiencing decline in one's cognitive abilities is among the most feared aspects of growing old [53]. Age-related cognitive decline carries a huge personal, societal, and financial cost both in pathological ageing (such as dementias) and also within the non-clinical majority of the population. A projected 152 million people worldwide will suffer from dementia by 2050 [3]. The early stages of cognitive decline are much more prevalent than dementia, and can still impose serious limitations of performance on everyday activities, independence, and quality of life in older age [5], [60], [80]. Cognitive decline also predicts poorer health, adherence to medical regimens, and financial decision-making, and can herald dementia, illness, and death [6], [40]. Of course, when seeking to understand why some people experience more severe cognitive ageing than others, researchers have turned to the organ of thinking for clues about the nature, possible mechanisms, and determinants that might underpin more and less successful cognitive agers. However, that organ is relatively inaccessible, a limitation partly alleviated by advances in neuroimaging. Here we discuss lessons for cognitive and brain ageing that have come from neuroimaging research (especially structural brain imaging), what neuroimaging still has left to teach us, and our views on possible ways forward in this multidisciplinary field.