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

Normative structural connectome constrains spreading transient brain activity in generalized epilepsy

BACKGROUND

Genetic generalized epilepsy is characterized by transient episodes of spontaneous abnormal neural activity in anatomically distributed brain regions that ultimately propagate to wider areas. However, the connectome-based mechanisms shaping these abnormalities remain largely unknown. We aimed to investigate how the normative structural connectome constrains abnormal brain activity spread in genetic generalized epilepsy with generalized tonic-clonic seizure (GGE-GTCS).

METHODS

Abnormal transient activity patterns between individuals with GGE-GTCS (n = 97) and healthy controls (n = 141) were estimated from the amplitude of low-frequency fluctuations measured by resting-state functional MRI. The normative structural connectome was derived from diffusion-weighted images acquired in an independent cohort of healthy adults (n = 326). Structural neighborhood analysis was applied to assess the degree of constraints between activity vulnerability and structural connectome. Dominance analysis was used to determine the potential molecular underpinnings of these constraints. Furthermore, a network-based diffusion model was utilized to simulate the spread of pathology and identify potential disease epicenters.

RESULTS

Brain activity abnormalities among patients with GGE-GTCS were primarily located in the temporal, cingulate, prefrontal, and parietal cortices. The collective abnormality of structurally connected neighbors significantly predicted regional activity abnormality, indicating that white matter network architecture constrains aberrant activity patterns. Molecular fingerprints, particularly laminar differentiation and neurotransmitter receptor profiles, constituted key predictors of these connectome-constrained activity abnormalities. Network-based diffusion modeling effectively replicated transient pathological activity spreading patterns, identifying the limbic-temporal, dorsolateral prefrontal, and occipital cortices as putative disease epicenters. These results were robust across different clinical factors and individual patients.

CONCLUSIONS

Our findings suggest that the structural connectome shapes the spatial patterning of brain activity abnormalities, advancing our understanding of the network-level mechanisms underlying vulnerability to abnormal brain activity onset and propagation in GGE-GTCS.

Two Axes of White Matter Development

Despite decades of neuroimaging research, how white matter develops along the length of major tracts in humans remains unknown. Here, we identify fundamental patterns of white matter maturation by examining developmental variation along major, long-range cortico-cortical tracts in youth ages 5-23 years using diffusion MRI from three large-scale, cross-sectional datasets (total N = 2,710). Across datasets, we delineate two replicable axes of human white matter development. First, we find a deep-to-superficial axis, in which superficial tract regions near the cortical surface exhibit greater age-related change than deep tract regions. Second, we demonstrate that the development of superficial tract regions aligns with the cortical hierarchy defined by the sensorimotor-association axis, with tract ends adjacent to sensorimotor cortices maturing earlier than those adjacent to association cortices. These results reveal developmental variation along tracts that conventional tract-average analyses have previously obscured, challenging the implicit assumption that white matter tracts mature uniformly along their length. Such developmental variation along tracts may have functional implications, including mitigating ephaptic coupling in densely packed deep tract regions and tuning neural synchrony through hierarchical development in superficial tract regions - ultimately refining neural transmission in youth.

Temporal Lobe Epilepsy Perturbs the Brain-Wide Excitation-Inhibition Balance: Associations with Microcircuit Organization, Clinical Parameters, and Cognitive Dysfunction

Excitation-inhibition (E/I) imbalance is theorized as a key mechanism in the pathophysiology of epilepsy, with ample research focusing on elucidating its cellular manifestations. However, few studies investigate E/I imbalance at the macroscale, whole-brain level, and its microcircuit-level mechanisms and clinical significance remain incompletely understood. Here, the Hurst exponent, an index of the E/I ratio, is computed from resting-state fMRI time series, and microcircuit parameters are simulated using biophysical models. A broad decrease in the Hurst exponent is observed in pharmaco-resistant temporal lobe epilepsy (TLE), suggesting more excitable network dynamics. Connectome decoders point to temporolimbic and frontocentral cortices as plausible network epicenters of E/I imbalance. Furthermore, computational simulations reveal that enhancing cortical excitability in TLE reflects atypical increases in recurrent connection strength of local neuronal ensembles. Mixed cross-sectional and longitudinal analyses show stronger E/I ratio elevation in patients with longer disease duration, more frequent electroclinical seizures as well as interictal epileptic spikes, and worse cognitive functioning. Hurst exponent-informed classifiers discriminate patients from healthy controls with high accuracy (72.4% [57.5%-82.5%]). Replicated in an independent dataset, this work provides in vivo evidence of a macroscale shift in E/I balance in TLE patients and points to progressive functional imbalances that relate to cognitive decline.

The architecture of the human default mode network explored through cytoarchitecture, wiring and signal flow

The default mode network (DMN) is implicated in many aspects of complex thought and behavior. Here, we leverage postmortem histology and in vivo neuroimaging to characterize the anatomy of the DMN to better understand its role in information processing and cortical communication. Our results show that the DMN is cytoarchitecturally heterogenous, containing cytoarchitectural types that are variably specialized for unimodal, heteromodal and memory-related processing. Studying diffusion-based structural connectivity in combination with cytoarchitecture, we found the DMN contains regions receptive to input from sensory cortex and a core that is relatively insulated from environmental input. Finally, analysis of signal flow with effective connectivity models showed that the DMN is unique amongst cortical networks in balancing its output across the levels of sensory hierarchies. Together, our study establishes an anatomical foundation from which accounts of the broad role the DMN plays in human brain function and cognition can be developed.

From histology to macroscale function in the human amygdala

The amygdala is a subcortical region in the mesiotemporal lobe that plays a key role in emotional and sensory functions. Conventional neuroimaging experiments treat this structure as a single, uniform entity, but there is ample histological evidence for subregional heterogeneity in microstructure and function. The current study characterized subregional structure-function coupling in the human amygdala, integrating post-mortem histology and in vivo MRI at ultra-high fields. Core to our work was a novel neuroinformatics approach that leveraged multiscale texture analysis as well as non-linear dimensionality reduction techniques to identify salient dimensions of microstructural variation in a 3D post-mortem histological reconstruction of the human amygdala. We observed two axes of subregional variation in this region, describing inferior-superior as well as mediolateral trends in microstructural differentiation that in part recapitulated established atlases of amygdala subnuclei. Translating our approach to in vivo MRI data acquired at 7 Tesla, we could demonstrate the generalizability of these spatial trends across 10 healthy adults. We then cross-referenced microstructural axes with functional blood-oxygen-level dependent (BOLD) signal analysis obtained during task-free conditions, and revealed a close association of structural axes with macroscale functional network embedding, notably the temporo-limbic, default mode, and sensory-motor networks. Our novel multiscale approach consolidates descriptions of amygdala anatomy and function obtained from histological and in vivo imaging techniques.

Geometric influences on the regional organization of the mammalian brain

The mammalian brain is comprised of anatomically and functionally distinct regions. Substantial work over the past century has pursued the generation of ever-more accurate maps of regional boundaries, using either expert judgement or data-driven clustering of functional, connectional, and/or architectonic properties. However, these approaches are often purely descriptive, have limited generalizability, and do not elucidate the underlying generative mechanisms that shape the regional organization of the brain. Here, we develop a novel approach that leverages a simple, hierarchical principle for generating a multiscale parcellation of any brain structure in any mammalian species using only its geometry. We show that this approach yields regions at any resolution scale that are more homogeneous than those defined in nearly all existing benchmark brain parcellations in use today across hundreds of anatomical, functional, cellular, and molecular brain properties measured in humans, macaques, marmosets, and mice. We additionally show how our method can be generalized to previously unstudied mammalian species for which no parcellations exist. Finally, we demonstrate how our approach captures the essence of a simple, hierarchical reaction-diffusion mechanism, in which the geometry of a brain structure shapes the spatial expression of putative patterning molecules linked to the formation of distinct regions through development. Our findings point to a highly conserved and universal influence of geometry on the regional organization of the mammalian brain.

Heterochronous laminar maturation in the human prefrontal cortex

The human prefrontal cortex (PFC) exhibits markedly protracted developmental plasticity, yet whether reductions in plasticity occur synchronously across prefrontal cortical layers is unclear. Animal studies have shown that intracortical myelin consolidates neural circuits to close periods of plasticity. Here, we use quantitative myelin imaging collected from youth (ages 10-32 years) at ultra-high field (7T) to investigate whether deep and superficial PFC layers exhibit different timeframes of plasticity. We find that myelin matures along a deep-to-superficial axis in the PFC; this axis of maturational timing is expressed to a different extent in cytoarchitecturally distinct regions along the frontal cortical hierarchy. By integrating myelin mapping with electroencephalogram and cognitive phenotyping, we provide evidence that deep and superficial prefrontal myelin dissociably impact timescales of neural activity, task learning rates, and cognitive processing speed. Heterochronous maturation across deep and superficial layers is an underrecognized mechanism through which association cortex balances cognitively-relevant increases in circuit stability and efficiency with extended neuroplasticity.

Brain maps of general cognitive function and spatial correlations with neurobiological cortical profiles

In this paper, we attempt to answer two questions: 1) which regions of the human brain, in terms of morphometry, are most strongly related to individual differences in domain-general cognitive functioning (g)? and 2) what are the underlying neurobiological properties of those regions? We meta-analyse vertex-wise g-cortical morphometry (volume, surface area, thickness, curvature and sulcal depth) associations using data from 3 cohorts: the UK Biobank (UKB), Generation Scotland (GenScot), and the Lothian Birth Cohort 1936 (LBC1936), with the meta-analytic N = 38,379 (age range = 44 to 84 years old). These g-morphometry associations vary in magnitude and direction across the cortex (|β| range = -0.12 to 0.17 across morphometry measures) and show good cross-cohort agreement (mean spatial correlation r = 0.57, SD = 0.18). Then, to address (2), we bring together existing - and derive new - cortical maps of 33 neurobiological characteristics from multiple modalities (including neurotransmitter receptor densities, gene expression, functional connectivity, metabolism, and cytoarchitectural similarity). We discover that these 33 profiles spatially covary along four major dimensions of cortical organisation (accounting for 65.9% of the variance) and denote aspects of neurobiological scaffolding that underpin the spatial patterning of MRI-cognitive associations we observe (significant |r| range = 0.21 to 0.56). Alongside the cortical maps from these analyses, which we make openly accessible, we provide a compendium of cortex-wide and within-region spatial correlations among general and specific facets of brain cortical organisation and higher order cognitive functioning, which we hope will serve as a framework for analysing other aspects of behaviour-brain MRI associations.

Microstructural asymmetry in the human cortex

The human cerebral cortex shows hemispheric asymmetry, yet the microstructural basis of this asymmetry remains incompletely understood. Here, we probe layer-specific microstructural asymmetry using one post-mortem male brain. Overall, anterior and posterior regions show leftward and rightward asymmetry respectively, but this pattern varies across cortical layers. A similar anterior-posterior pattern is observed using in vivo Human Connectome Project (N = 1101) T1w/T2w microstructural data, with average cortical asymmetry showing the strongest similarity with post-mortem-based asymmetry of layer III. Moreover, microstructural asymmetry is found to be heritable, varies as a function of age and sex, and corresponds to intrinsic functional asymmetry. We also observe a differential association of language and markers of mental health with microstructural asymmetry patterns at the individual level, illustrating a functional divergence between inferior-superior and anterior-posterior microstructural axes, possibly anchored in development. Last, we could show concordant evidence with alternative in vivo microstructural measures: magnetization transfer (N = 286) and quantitative T1 (N = 50). Together, our study highlights microstructural asymmetry in the human cortex and its functional and behavioral relevance.

Molecular signatures of cortical expansion in the human foetal brain

The third trimester of human gestation is characterised by rapid increases in brain volume and cortical surface area. Recent studies have revealed a remarkable molecular diversity across the prenatal cortex but little is known about how this diversity translates into the differential rates of cortical expansion observed during gestation. We present a digital resource, μBrain, to facilitate knowledge translation between molecular and anatomical descriptions of the prenatal brain. Using μBrain, we evaluate the molecular signatures of preferentially-expanded cortical regions, quantified in utero using magnetic resonance imaging. Our findings demonstrate a spatial coupling between areal differences in the timing of neurogenesis and rates of neocortical expansion during gestation. We identify genes, upregulated from mid-gestation, that are highly expressed in rapidly expanding neocortex and implicated in genetic disorders with cognitive sequelae. The μBrain atlas provides a tool to comprehensively map early brain development across domains, model systems and resolution scales.

Direct segmentation of cortical cytoarchitectonic domains using ultra-high-resolution whole-brain diffusion MRI

We assess the potential of detecting cortical laminar patterns and areal borders by directly clustering voxel values of microstructural parameters derived from high-resolution mean apparent propagator (MAP) magnetic resonance imaging (MRI), as an alternative to conventional template-warping-based cortical parcellation methods. We acquired MAP-MRI data with 200μm resolution in a fixed macaque monkey brain. To improve the sensitivity to cortical layers, we processed the data with a local anisotropic Gaussian filter determined voxel-wise by the plane tangent to the cortical surface. We directly clustered all cortical voxels using only the MAP-derived microstructural imaging biomarkers, with no information regarding their relative spatial location or dominant diffusion orientations. MAP-based 3D cytoarchitectonic segmentation revealed laminar patterns similar to those observed in the corresponding histological images. Moreover, transition regions between these laminar patterns agreed more accurately with histology than the borders between cortical areas estimated using conventional atlas/template-warping cortical parcellation. By cross-tabulating all cortical labels in the atlas- and MAP-based segmentations, we automatically matched the corresponding MAP-derived clusters (i.e., cytoarchitectonic domains) across the left and right hemispheres. Our results demonstrate that high-resolution MAP-MRI biomarkers can effectively delineate three-dimensional cortical cytoarchitectonic domains in single individuals. Their intrinsic tissue microstructural contrasts enable the construction of whole-brain mesoscopic cortical atlases.

Gradients of Brain Organization: Smooth Sailing from Methods Development to User Community

Multimodal neuroimaging grants a powerful in vivo window into the structure and function of the human brain. Recent methodological and conceptual advances have enabled investigations of the interplay between large-scale spatial trends - or gradients - in brain structure and function, offering a framework to unify principles of brain organization across multiple scales. Strong community enthusiasm for these techniques has been instrumental in their widespread adoption and implementation to answer key questions in neuroscience. Following a brief review of current literature on this framework, this perspective paper will highlight how pragmatic steps aiming to make gradient methods more accessible to the community propelled these techniques to the forefront of neuroscientific inquiry. More specifically, we will emphasize how interest for gradient methods was catalyzed by data sharing, open-source software development, as well as the organization of dedicated workshops led by a diverse team of early career researchers. To this end, we argue that the growing excitement for brain gradients is the result of coordinated and consistent efforts to build an inclusive community and can serve as a case in point for future innovations and conceptual advances in neuroinformatics. We close this perspective paper by discussing challenges for the continuous refinement of neuroscientific theory, methodological innovation, and real-world translation to maintain our collective progress towards integrated models of brain organization.

Abnormal structural‒functional coupling patterning in progressive supranuclear palsy is associated with diverse gradients and histological features

The anatomy of the brain supports inherent processes, fostering mental abilities and eventually facilitating adaptive behavior. Recent studies have shown that progressive supranuclear palsy (PSP) is accompanied by alterations in functional and structural networks. However, how the structure and function of PSP coordinates change is not clear, and the relationships between structural‒functional coupling (SFC) and the gradient of hierarchical structure and cellular histology remain largely unknown. Here, we use neuroimaging data from two independent cohorts and a public histological dataset to investigate the relationships among the cellular histology, hierarchical structure, and SFC of PSP patients. We find that the SFC of the entire cortex in PSP is severely disrupted, with higher coupling in the visual network (VN). Moreover, coupling differences in PSP follow a macroscopic organizational principle from unimodal to transmodal gradients. Finally, we elucidate greater laminar differentiation in VN regions sensitive to SFC changes in PSP, which is related mainly to the higher cellular density and smaller size of the internal-granular layer. In conclusion, our findings provide an interpretable framework for understanding SFC changes in PSP and provide new insights into the consistency of structural and functional changes in PSP regarding hierarchical structure and cellular histology.

Contracted functional connectivity profiles in autism

OBJECTIVE

Autism spectrum disorder (ASD) is a neurodevelopmental condition that is associated with atypical brain network organization, with prior work suggesting differential connectivity alterations with respect to functional connection length. Here, we tested whether functional connectopathy in ASD specifically relates to disruptions in long- relative to short-range functional connections. Our approach combined functional connectomics with geodesic distance mapping, and we studied associations to macroscale networks, microarchitectural patterns, as well as socio-demographic and clinical phenotypes.

METHODS

We studied 211 males from three sites of the ABIDE-I dataset comprising 103 participants with an ASD diagnosis (mean ± SD age = 20.8 ± 8.1 years) and 108 neurotypical controls (NT, 19.2 ± 7.2 years). For each participant, we computed cortex-wide connectivity distance (CD) measures by combining geodesic distance mapping with resting-state functional connectivity profiling. We compared CD between ASD and NT participants using surface-based linear models, and studied associations with age, symptom severity, and intelligence scores. We contextualized CD alterations relative to canonical networks and explored spatial associations with functional and microstructural cortical gradients as well as cytoarchitectonic cortical types.

RESULTS

Compared to NT, ASD participants presented with widespread reductions in CD, generally indicating shorter average connection length and thus suggesting reduced long-range connectivity but increased short-range connections. Peak reductions were localized in transmodal systems (i.e., heteromodal and paralimbic regions in the prefrontal, temporal, and parietal and temporo-parieto-occipital cortex), and effect sizes correlated with the sensory-transmodal gradient of brain function. ASD-related CD reductions appeared consistent across inter-individual differences in age and symptom severity, and we observed a positive correlation of CD to IQ scores.

LIMITATIONS

Despite rigorous harmonization across the three different acquisition sites, heterogeneity in autism poses a potential limitation to the generalizability of our results. Additionally, we focussed male participants, warranting future studies in more balanced cohorts.

CONCLUSIONS

Our study showed reductions in CD as a relatively stable imaging phenotype of ASD that preferentially impacted paralimbic and heteromodal association systems. CD reductions in ASD corroborate previous reports of ASD-related imbalance between short-range overconnectivity and long-range underconnectivity.

A next-generation, histological atlas of the human brain and its application to automated brain MRI segmentation

Magnetic resonance imaging (MRI) is the standard tool to image the human brain in vivo. In this domain, digital brain atlases are essential for subject-specific segmentation of anatomical regions of interest (ROIs) and spatial comparison of neuroanatomy from different subjects in a common coordinate frame. High-resolution, digital atlases derived from histology (e.g., Allen atlas [7], BigBrain [13], Julich [15]), are currently the state of the art and provide exquisite 3D cytoarchitectural maps, but lack probabilistic labels throughout the whole brain. Here we present NextBrain, a next-generation probabilistic atlas of human brain anatomy built from serial 3D histology and corresponding highly granular delineations of five whole brain hemispheres. We developed AI techniques to align and reconstruct ~10,000 histological sections into coherent 3D volumes with joint geometric constraints (no overlap or gaps between sections), as well as to semi-automatically trace the boundaries of 333 distinct anatomical ROIs on all these sections. Comprehensive delineation on multiple cases enabled us to build the first probabilistic histological atlas of the whole human brain. Further, we created a companion Bayesian tool for automated segmentation of the 333 ROIs in any in vivo or ex vivo brain MRI scan using the NextBrain atlas. We showcase two applications of the atlas: automated segmentation of ultra-high-resolution ex vivo MRI and volumetric analysis of Alzheimer's disease and healthy brain ageing based on ~4,000 publicly available in vivo MRI scans. We publicly release: the raw and aligned data (including an online visualisation tool); the probabilistic atlas; the segmentation tool; and ground truth delineations for a 100 μm isotropic ex vivo hemisphere (that we use for quantitative evaluation of our segmentation method in this paper). By enabling researchers worldwide to analyse brain MRI scans at a superior level of granularity without manual effort or highly specific neuroanatomical knowledge, NextBrain holds promise to increase the specificity of MRI findings and ultimately accelerate our quest to understand the human brain in health and disease.

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.