A multi-scale cortical wiring space links cellular architecture and functional dynamics in the human brain

MULTIMODAL NEURAL CORRELATES OF CHILDHOOD PSYCHOPATHOLOGY

Complex structural and functional changes occurring in typical and atypical development necessitate multidimensional approaches to better understand the risk of developing psychopathology. Here, we simultaneously examined structural and functional brain network patterns in relation to dimensions of psychopathology in the Adolescent Brain Cognitive Development dataset. Several components were identified, recapitulating the psychopathology hierarchy, with the general psychopathology (p) factor explaining most covariance with multimodal imaging features, while the internalizing, externalizing, and neurodevelopmental dimensions were each associated with distinct morphological and functional connectivity signatures. Connectivity signatures associated with the p factor and neurodevelopmental dimensions followed the sensory-to-transmodal axis of cortical organization, which is related to the emergence of complex cognition and risk for psychopathology. Results were consistent in two separate data subsamples, supporting generalizability, and robust to variations in analytical parameters. Our findings help in better understanding biological mechanisms underpinning dimensions of psychopathology, and could provide brain-based vulnerability markers.

Transcriptional expression patterns of the cortical morphometric similarity network in progressive supranuclear palsy

BACKGROUND

It has been demonstrated that progressive supranuclear palsy (PSP) correlates with structural abnormalities in several distinct regions of the brain. However, whether there are changes in the morphological similarity network (MSN) and the relationship between changes in brain structure and gene expression remain largely unknown.

METHODS

We used two independent cohorts (discovery dataset: PSP: 51, healthy controls (HC): 82; replication dataset: PSP: 53, HC: 55) for MSN analysis and comparing the longitudinal changes in the MSN of PSP. Then, we applied partial least squares regression to determine the relationships between changes in MSN and spatial transcriptional features and identified specific genes associated with MSN differences in PSP. We further investigated the biological processes enriched in PSP-associated genes and the cellular characteristics of these genes, and finally, we performed an exploratory analysis of the relationship between MSN changes and neurotransmitter receptors.

RESULTS

We found that the MSN in PSP patients was mainly decreased in the frontal and temporal cortex but increased in the occipital cortical region. This difference is replicable. In longitudinal studies, MSN differences are mainly manifested in the frontal and parietal regions. Furthermore, the expression pattern associated with MSN changes in PSP involves genes implicated in astrocytes and excitatory and inhibitory neurons and is functionally enriched in neuron-specific biological processes related to synaptic signaling. Finally, we found that the changes in MSN were mainly negatively correlated with the levels of serotonin, norepinephrine, and opioid receptors.

CONCLUSIONS

These results have enhanced our understanding of the microscale genetic and cellular mechanisms responsible for large-scale morphological abnormalities in PSP patients, suggesting potential targets for future therapeutic trials.

Volume electron microscopy analysis of synapses in primary regions of the human cerebral cortex

Functional and structural studies investigating macroscopic connectivity in the human cerebral cortex suggest that high-order associative regions exhibit greater connectivity compared to primary ones. However, the synaptic organization of these brain regions remains unexplored. In the present work, we conducted volume electron microscopy to investigate the synaptic organization of the human brain obtained at autopsy. Specifically, we examined layer III of Brodmann areas 17, 3b, and 4, as representative areas of primary visual, somatosensorial, and motor cortex. Additionally, we conducted comparative analyses with our previous datasets of layer III from temporopolar and anterior cingulate associative cortical regions (Brodmann areas 24, 38, and 21). 9,690 synaptic junctions were 3D reconstructed, showing that certain synaptic characteristics are specific to particular regions. The number of synapses per volume, the proportion of the postsynaptic targets, and the synaptic size may distinguish one region from another, regardless of whether they are associative or primary cortex. By contrast, other synaptic characteristics were common to all analyzed regions, such as the proportion of excitatory and inhibitory synapses, their shapes, their spatial distribution, and a higher proportion of synapses located on dendritic spines. The present results provide further insights into the synaptic organization of the human cerebral cortex.

Cortical gene expression architecture links healthy neurodevelopment to the imaging, transcriptomics and genetics of autism and schizophrenia

Human brain organization involves the coordinated expression of thousands of genes. For example, the first principal component (C1) of cortical transcription identifies a hierarchy from sensorimotor to association regions. In this study, optimized processing of the Allen Human Brain Atlas revealed two new components of cortical gene expression architecture, C2 and C3, which are distinctively enriched for neuronal, metabolic and immune processes, specific cell types and cytoarchitectonics, and genetic variants associated with intelligence. Using additional datasets (PsychENCODE, Allen Cell Atlas and BrainSpan), we found that C1-C3 represent generalizable transcriptional programs that are coordinated within cells and differentially phased during fetal and postnatal development. Autism spectrum disorder and schizophrenia were specifically associated with C1/C2 and C3, respectively, across neuroimaging, differential expression and genome-wide association studies. Evidence converged especially in support of C3 as a normative transcriptional program for adolescent brain development, which can lead to atypical supragranular cortical connectivity in people at high genetic risk for schizophrenia.

Neural correlates of obesity across the lifespan

Associations between brain and obesity are bidirectional: changes in brain structure and function underpin over-eating, while chronic adiposity leads to brain atrophy. Investigating brain-obesity interactions across the lifespan can help better understand these relationships. This study explores the interaction between obesity and cortical morphometry in children, young adults, adults, and older adults. We also investigate the genetic, neurochemical, and cognitive correlates of the brain-obesity associations. Our findings reveal a pattern of lower cortical thickness in fronto-temporal brain regions associated with obesity across all age cohorts and varying age-dependent patterns in the remaining brain regions. In adults and older adults, obesity correlates with neurochemical changes and expression of inflammatory and mitochondrial genes. In children and older adults, adiposity is associated with modifications in brain regions involved in emotional and attentional processes. Thus, obesity might originate from cognitive changes during early adolescence, leading to neurodegeneration in later life through mitochondrial and inflammatory mechanisms.

Body mass index-dependent shifts along large-scale gradients in human cortical organization explain dietary regulatory success

Making healthy dietary choices is essential for keeping weight within a normal range. Yet many people struggle with dietary self-control despite good intentions. What distinguishes neural processing in those who succeed or fail to implement healthy eating goals? Does this vary by weight status? To examine these questions, we utilized an analytical framework of gradients that characterize systematic spatial patterns of large-scale neural activity, which have the advantage of considering the entire suite of processes subserving self-control and potential regulatory tactics at the whole-brain level. Using an established laboratory food task capturing brain responses in natural and regulatory conditions (N = 123), we demonstrate that regulatory changes of dietary brain states in the gradient space predict individual differences in dietary success. Better regulators required smaller shifts in brain states to achieve larger goal-consistent changes in dietary behaviors, pointing toward efficient network organization. This pattern was most pronounced in individuals with lower weight status (low-BMI, body mass index) but absent in high-BMI individuals. Consistent with prior work, regulatory goals increased activity in frontoparietal brain circuits. However, this shift in brain states alone did not predict variance in dietary success. Instead, regulatory success emerged from combined changes along multiple gradients, showcasing the interplay of different large-scale brain networks subserving dietary control and possible regulatory strategies. Our results provide insights into how the brain might solve the problem of dietary control: Dietary success may be easier for people who adopt modes of large-scale brain activation that do not require significant reconfigurations across contexts and goals.

Whole-brain structural connectome asymmetry in autism

Autism spectrum disorder is a common neurodevelopmental condition that manifests as a disruption in sensory and social skills. Although it has been shown that the brain morphology of individuals with autism is asymmetric, how this differentially affects the structural connectome organization of each hemisphere remains under-investigated. We studied whole-brain structural connectivity-based brain asymmetry in individuals with autism using diffusion magnetic resonance imaging obtained from the Autism Brain Imaging Data Exchange initiative. By leveraging dimensionality reduction techniques, we constructed low-dimensional representations of structural connectivity and calculated their asymmetry index. Comparing the asymmetry index between individuals with autism and neurotypical controls, we found atypical structural connectome asymmetry in the sensory and default-mode regions, particularly showing weaker asymmetry towards the right hemisphere in autism. Network communication provided topological underpinnings by demonstrating that the inferior temporal cortex and limbic and frontoparietal regions showed reduced global network communication efficiency and decreased send-receive network navigation in the inferior temporal and lateral visual cortices in individuals with autism. Finally, supervised machine learning revealed that structural connectome asymmetry could be used as a measure for predicting communication-related autistic symptoms and nonverbal intelligence. Our findings provide insights into macroscale structural connectome alterations in autism and their topological underpinnings.

Connectome-wide structure-function coupling models implicate polysynaptic alterations in autism

Autism spectrum disorder (ASD) is one of the most common neurodevelopmental diagnoses. Although incompletely understood, structural and functional network alterations are increasingly recognized to be at the core of the condition. We utilized multimodal imaging and connectivity modeling to study structure-function coupling in ASD and probed mono- and polysynaptic mechanisms on structurally-governed network function. We examined multimodal magnetic resonance imaging data in 80 ASD and 61 neurotypical controls from the Autism Brain Imaging Data Exchange (ABIDE) II initiative. We predicted intrinsic functional connectivity from structural connectivity data in each participant using a Riemannian optimization procedure that varies the times that simulated signals can unfold along tractography-derived personalized connectomes. In both ASD and neurotypical controls, we observed improved structure-function prediction at longer diffusion time scales, indicating better modeling of brain function when polysynaptic mechanisms are accounted for. Prediction accuracy differences (∆prediction accuracy) were marked in transmodal association systems, such as the default mode network, in both neurotypical controls and ASD. Differences were, however, lower in ASD in a polysynaptic regime at higher simulated diffusion times. We compared regional differences in ∆prediction accuracy between both groups to assess the impact of polysynaptic communication on structure-function coupling. This analysis revealed that between-group differences in ∆prediction accuracy followed a sensory-to-transmodal cortical hierarchy, with an increased gap between controls and ASD in transmodal compared to sensory/motor systems. Multivariate associative techniques revealed that structure-function differences reflected inter-individual differences in autistic symptoms and verbal as well as non-verbal intelligence. Our network modeling approach sheds light on atypical structure-function coupling in autism, and suggests that polysynaptic network mechanisms are implicated in the condition and that these can help explain its wide range of associated symptoms.

Pattern of ventral temporal lobe interconnections in rhesus macaques

The entorhinal cortex (EC, A28) is linked through reciprocal pathways with nearby perirhinal and visual, auditory, and multimodal association cortices in the temporal lobe, in pathways associated with the flow of information for memory processing. The density and laminar organization of these pathways is not well understood in primates. We studied interconnections within the ventral temporal lobe in young adult rhesus monkeys of both sexes with the aid of neural tracers injected in temporal areas (Ts1, Ts2, TE1, area 36, temporal polar area TPro, and area 28) to determine the density and laminar distribution of projection neurons within the temporal lobe. These temporal areas can be categorized into three different cortical types based on their laminar architecture: the sensory association areas Ts1, Ts2, and TE1 have six layers (eulaminate); the perirhinal limbic areas TPro and area 36 have an incipient layer IV (dysgranular); and area 28 lacks layer IV (agranular). We found that (1) temporal areas that are similar in laminar architecture by cortical type are strongly interconnected, and (2) the laminar pattern of connections is dependent on the difference in cortical laminar structure between linked areas. Thus, agranular A28 is more strongly connected with other agranular/dysgranular areas than with eulaminate cortices. Further, A28 predominantly projected via feedback-like pathways that originated in the deep layers, and received feedforward-like projections from areas of greater laminar differentiation, which emanated from the upper layers. Our results are consistent with the Structural Model, which relates the density and laminar distribution of connections to the relationship of the laminar structure between the linked areas. These connections were viewed in the context of the inhibitory microenvironment of A28, which is the key recipient of pathways from the cortex and of the output of hippocampus. Our findings revealed a higher population of calretinin (CR)-expressing neurons in EC, with a significantly higher density in its lateral division. Medial EC had a higher density of CR neurons in the deep layers, particularly in layer Va. In contrast, parvalbumin (PV) neurons were more densely distributed in the deep layers of the lateral subdivisions of rostral EC, especially in layer Va, whereas the densities of calbindin (CB) neurons in the medial and lateral EC were comparable in all layers, except for layer IIIa, in which medial EC had a higher CB population than the lateral. The pattern of connections in the inhibitory microenvironment of EC, which sends and receives input from the hippocampus, may shed light on signal propagation in this network associated with diverse aspects of memory, and disruptions in neurologic and psychiatric diseases that affect this region.

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.

Schizophrenia Polygenic Risk During Typical Development Reflects Multiscale Cortical Organization

Background

Schizophrenia is widely recognized as a neurodevelopmental disorder. Abnormal cortical development in otherwise typically developing children and adolescents may be revealed using polygenic risk scores for schizophrenia (PRS-SCZ).

Methods

We assessed PRS-SCZ and cortical morphometry in typically developing children and adolescents (3-21 years, 46.8% female) using whole-genome genotyping and T1-weighted magnetic resonance imaging (n = 390) from the PING (Pediatric Imaging, Neurocognition, and Genetics) cohort. We contextualized the findings using 1) age-matched transcriptomics, 2) histologically defined cytoarchitectural types and functionally defined networks, and 3) case-control differences of schizophrenia and other major psychiatric disorders derived from meta-analytic data of 6 ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) working groups, including a total of 12,876 patients and 15,670 control participants.

Results

Higher PRS-SCZ was associated with greater cortical thickness, which was most prominent in areas with heightened gene expression of dendrites and synapses. PRS-SCZ-related increases in vertexwise cortical thickness were mainly distributed in association cortical areas, particularly the ventral attention network, while relatively sparing koniocortical type cortex (i.e., primary sensory areas). The large-scale pattern of cortical thickness increases related to PRS-SCZ mirrored the pattern of cortical thinning in schizophrenia and mood-related psychiatric disorders derived from the ENIGMA consortium. Age group models illustrate a possible trajectory from PRS-SCZ-associated cortical thickness increases in early childhood toward thinning in late adolescence, with the latter resembling the adult brain phenotype of schizophrenia.

Conclusions

Collectively, combining imaging genetics with multiscale mapping, our work provides novel insight into how genetic risk for schizophrenia affects the cortex early in life.

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.

3D synaptic organization of layer III of the human anterior cingulate and temporopolar cortex

The human anterior cingulate and temporopolar cortices have been proposed as highly connected nodes involved in high-order cognitive functions, but their synaptic organization is still basically unknown due to the difficulties involved in studying the human brain. Using Focused Ion Beam/Scanning Electron Microscopy (FIB/SEM) to study the synaptic organization of the human brain obtained with a short post-mortem delay allows excellent results to be obtained. We have used this technology to analyze layer III of the anterior cingulate cortex (Brodmann area 24) and the temporopolar cortex, including the temporal pole (Brodmann area 38 ventral and dorsal) and anterior middle temporal gyrus (Brodmann area 21). Our results, based on 6695 synaptic junctions fully reconstructed in 3D, revealed that Brodmann areas 24, 21 and ventral area 38 showed similar synaptic density and synaptic size, whereas dorsal area 38 displayed the highest synaptic density and the smallest synaptic size. However, the proportion of the different types of synapses (excitatory and inhibitory), the postsynaptic targets, and the shapes of excitatory and inhibitory synapses were similar, regardless of the region examined. These observations indicate that certain aspects of the synaptic organization are rather homogeneous, whereas others show specific variations across cortical regions.

The Cellular Underpinnings of the Human Cortical Connectome

The functional properties of the human brain arise, in part, from the vast assortment of cell types that pattern the cortex. The cortical sheet can be broadly divided into distinct networks, which are further embedded into processing streams, or gradients, that extend from unimodal systems through higher-order association territories. Here, using transcriptional data from the Allen Human Brain Atlas, we demonstrate that imputed cell type distributions are spatially coupled to the functional organization of cortex, as estimated through fMRI. Cortical cellular profiles follow the macro-scale organization of the functional gradients as well as the associated large-scale networks. Distinct cellular fingerprints were evident across networks, and a classifier trained on post-mortem cell-type distributions was able to predict the functional network allegiance of cortical tissue samples. These data indicate that the in vivo organization of the cortical sheet is reflected in the spatial variability of its cellular composition.

In vivo mapping of pharmacologically induced functional reorganization onto the human brain's neurotransmitter landscape

To understand how pharmacological interventions can exert their powerful effects on brain function, we need to understand how they engage the brain's rich neurotransmitter landscape. Here, we bridge microscale molecular chemoarchitecture and pharmacologically induced macroscale functional reorganization, by relating the regional distribution of 19 neurotransmitter receptors and transporters obtained from positron emission tomography, and the regional changes in functional magnetic resonance imaging connectivity induced by 10 different mind-altering drugs: propofol, sevoflurane, ketamine, lysergic acid diethylamide (LSD), psilocybin, N,N-Dimethyltryptamine (DMT), ayahuasca, 3,4-methylenedioxymethamphetamine (MDMA), modafinil, and methylphenidate. Our results reveal a many-to-many mapping between psychoactive drugs' effects on brain function and multiple neurotransmitter systems. The effects of both anesthetics and psychedelics on brain function are organized along hierarchical gradients of brain structure and function. Last, we show that regional co-susceptibility to pharmacological interventions recapitulates co-susceptibility to disorder-induced structural alterations. Collectively, these results highlight rich statistical patterns relating molecular chemoarchitecture and drug-induced reorganization of the brain's functional architecture.

Language network lateralization is reflected throughout the macroscale functional organization of cortex

Hemispheric specialization is a fundamental feature of human brain organization. However, it is not yet clear to what extent the lateralization of specific cognitive processes may be evident throughout the broad functional architecture of cortex. While the majority of people exhibit left-hemispheric language dominance, a substantial minority of the population shows reverse lateralization. Using twin and family data from the Human Connectome Project, we provide evidence that atypical language dominance is associated with global shifts in cortical organization. Individuals with atypical language organization exhibit corresponding hemispheric differences in the macroscale functional gradients that situate discrete large-scale networks along a continuous spectrum, extending from unimodal through association territories. Analyses reveal that both language lateralization and gradient asymmetries are, in part, driven by genetic factors. These findings pave the way for a deeper understanding of the origins and relationships linking population-level variability in hemispheric specialization and global properties of cortical organization.

Atypical connectome topography and signal flow in temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is one of the most common pharmaco-resistant epilepsies in adults. While hippocampal pathology is the hallmark of this condition, emerging evidence indicates that brain alterations extend beyond the mesiotemporal epicenter and affect macroscale brain function and cognition. We studied macroscale functional reorganization in TLE, explored structural substrates, and examined cognitive associations. We investigated a multisite cohort of 95 patients with pharmaco-resistant TLE and 95 healthy controls using state-of-the-art multimodal 3T magnetic resonance imaging (MRI). We quantified macroscale functional topographic organization using connectome dimensionality reduction techniques and estimated directional functional flow using generative models of effective connectivity. We observed atypical functional topographies in patients with TLE relative to controls, manifesting as reduced functional differentiation between sensory/motor networks and transmodal systems such as the default mode network, with peak alterations in bilateral temporal and ventromedial prefrontal cortices. TLE-related topographic changes were consistent in all three included sites and reflected reductions in hierarchical flow patterns between cortical systems. Integration of parallel multimodal MRI data indicated that these findings were independent of TLE-related cortical grey matter atrophy, but mediated by microstructural alterations in the superficial white matter immediately beneath the cortex. The magnitude of functional perturbations was robustly associated with behavioral markers of memory function. Overall, this work provides converging evidence for macroscale functional imbalances, contributing microstructural alterations, and their associations with cognitive dysfunction in TLE.

Connectional Hierarchy in Human Brain Revealed by Individual Variability of Functional Network Edges

The human cerebral cortex is connected by intricate inter-areal wiring at the macroscale. The cortical hierarchy from primary sensorimotor to higher-order association areas is a unifying organizational principle across various neurobiological properties; however, previous studies have not clarified whether the connections between cortical regions exhibit a similar hierarchical pattern. Here, we identify a connectional hierarchy indexed by inter-individual variability of functional connectivity edges, which continuously progresses along a hierarchical gradient from within-network connections to between-network edges connecting sensorimotor and association networks. We found that this connectional hierarchy of variability aligns with both hemodynamic and electromagnetic connectivity strength and is constrained by structural connectivity strength. Moreover, the patterning of connectional hierarchy is related to inter-regional similarity in transcriptional and neurotransmitter receptor profiles. Using the Neurosynth cognitive atlas and cortical vulnerability maps in 13 brain disorders, we found that the connectional hierarchy of variability is associated with similarity networks of cognitive relevance and that of disorder vulnerability. Finally, we found that the prominence of this hierarchical gradient of connectivity variability declines during youth. Together, our results reveal a novel hierarchal organizational principle at the connectional level that links multimodal and multiscale human connectomes to individual variability in functional connectivity.

Integration of white matter architecture to stereo-EEG better describes epileptic spike propagation

OBJECTIVE

Stereo-electroencephalography (SEEG)-derived epilepsy networks are used to better understand a patient's epilepsy; however, a unimodal approach provides an incomplete picture. We combine tractography and SEEG to determine the relationship between spike propagation and the white matter architecture and to improve our understanding of spike propagation mechanisms.

METHODS

Probablistic tractography from diffusion imaging (dMRI) of matched subjects from the Human Connectome Project (HCP) was combined with patient-specific SEEG-derived spike propagation networks. Two regions-of-interest (ROIs) with a significant spike propagation relationship constituted a Propagation Pair.

RESULTS

In 56 of 59 patients, Propagation Pairs were more often tract-connected as compared to all ROI pairs (p < 0.01; d = -1.91). The degree of spike propagation between tract-connected ROIs was greater (39 ± 21%) compared to tract-unconnected ROIs (31 ± 18%; p < 0.0001). Within the same network, ROIs receiving propagation earlier were more often tract-connected to the source (59.7%) as compared to late receivers (25.4%; p < 0.0001).

CONCLUSIONS

Brain regions involved in spike propagation are more likely to be connected by white matter tracts. Between nodes, presence of tracts suggests a direct course of propagation, whereas the absence of tracts suggests an indirect course of propagation.

SIGNIFICANCE

We demonstrate a logical and consistent relationship between spike propagation and the white matter architecture.

Transcriptional signatures of the cortical morphometric similarity network gradient in first-episode, treatment-naive major depressive disorder

Recent studies have shown that major depressive disorder (MDD) is accompanied by alterations in functional and structural network gradients. However, whether changes are present in the cortical morphometric similarity (MS) network gradient, and the relationship between alterations of the gradient and gene expression remains largely unknown. In this study, the MS network was constructed, and its gradient was calculated in 71 patients with first-episode, treatment-naive MDD, and 69 demographically matched healthy controls. Between-group comparisons were performed to investigate abnormalities in the MS network gradient, and partial least squares regression analysis was conducted to explore the association between gene expression profiles and MS network gradient-based alternations in MDD. We found that the gradient was primarily significantly decreased in sensorimotor regions in patients with MDD compared with healthy controls, and increased in visual-related regions. In addition, the altered principal MS network gradient in the left postcentral cortex and right lingual cortex exhibited significant correlations with symptom severity. The abnormal gradient pattern was spatially correlated with the brain-wide expression of genes enriched for neurobiologically relevant pathways, downregulated in the MDD postmortem brain, and preferentially expressed in different cell types and cortical layers. These results demonstrated alterations of the principal MS network gradient in MDD and suggested the molecular mechanisms for structural alternations underlying MDD.

Distributed harmonic patterns of structure-function dependence orchestrate human consciousness

A central question in neuroscience is how consciousness arises from the dynamic interplay of brain structure and function. Here we decompose functional MRI signals from pathological and pharmacologically-induced perturbations of consciousness into distributed patterns of structure-function dependence across scales: the harmonic modes of the human structural connectome. We show that structure-function coupling is a generalisable indicator of consciousness that is under bi-directional neuromodulatory control. We find increased structure-function coupling across scales during loss of consciousness, whether due to anaesthesia or brain injury, capable of discriminating between behaviourally indistinguishable sub-categories of brain-injured patients, tracking the presence of covert consciousness. The opposite harmonic signature characterises the altered state induced by LSD or ketamine, reflecting psychedelic-induced decoupling of brain function from structure and correlating with physiological and subjective scores. Overall, connectome harmonic decomposition reveals how neuromodulation and the network architecture of the human connectome jointly shape consciousness and distributed functional activation across scales.

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.

Coordinated cortical thickness alterations across six neurodevelopmental and psychiatric disorders

Neuropsychiatric disorders are increasingly conceptualized as overlapping spectra sharing multi-level neurobiological alterations. However, whether transdiagnostic cortical alterations covary in a biologically meaningful way is currently unknown. Here, we studied co-alteration networks across six neurodevelopmental and psychiatric disorders, reflecting pathological structural covariance. In 12,024 patients and 18,969 controls from the ENIGMA consortium, we observed that co-alteration patterns followed normative connectome organization and were anchored to prefrontal and temporal disease epicenters. Manifold learning revealed frontal-to-temporal and sensory/limbic-to-occipitoparietal transdiagnostic gradients, differentiating shared illness effects on cortical thickness along these axes. The principal gradient aligned with a normative cortical thickness covariance gradient and established a transcriptomic link to cortico-cerebello-thalamic circuits. Moreover, transdiagnostic gradients segregated functional networks involved in basic sensory, attentional/perceptual, and domain-general cognitive processes, and distinguished between regional cytoarchitectonic profiles. Together, our findings indicate that shared illness effects occur in a synchronized fashion and along multiple levels of hierarchical cortical organization.

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.

Cortical gradient of a human functional similarity network captured by the geometry of cytoarchitectonic organization

Mapping the functional topology from a multifaceted perspective and relating it to underlying cross-scale structural principles is crucial for understanding the structural-functional relationships of the cerebral cortex. Previous works have described a sensory-association gradient axis in terms of coupling relationships between structure and function, but largely based on single specific feature, and the mesoscopic underpinnings are rarely determined. Here we show a gradient pattern encoded in a functional similarity network based on data from Human Connectome Project and further link it to cytoarchitectonic organizing principles. The spatial distribution of the primary gradient follows an inferior-anterior to superior-posterior axis. The primary gradient demonstrates converging relationships with layer-specific microscopic gene expression and mesoscopic cortical layer thickness, and is captured by the geometric representation of a myelo- and cyto-architecture based laminar differentiation theorem, involving a dual origin theory. Together, these findings provide a gradient, which describes the functional topology, and more importantly, linking the macroscale functional landscape with mesoscale laminar differentiation principles.

Analysis of functional MRI data from the Human Connectome Project reveals the presence of a cortical gradient, which reflects functional topology and links the macroscale functional landscape with mesoscale laminar differentiation principles.

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.

An Open MRI Dataset For Multiscale Neuroscience

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

The cortical spectrum: A robust structural continuum in primate cerebral cortex revealed by histological staining and magnetic resonance imaging

High-level characterizations of the primate cerebral cortex sit between two extremes: on one end the cortical mantle is seen as a mosaic of structurally and functionally unique areas, and on the other it is seen as a uniform six-layered structure in which functional differences are defined solely by extrinsic connections. Neither of these extremes captures the crucial neuroanatomical finding: that the cortex exhibits systematic gradations in architectonic structure. These gradations have been shown to predict cortico-cortical connectivity, which in turn suggests powerful ways to ground connectomics in anatomical structure, and by extension cortical function. A challenge to widespread use of this concept is the labor-intensive and invasive nature of histological staining, which is the primary means of recognizing anatomical gradations. Here we show that a novel computational analysis technique can provide a coarse-grained picture of cortical variation. For each of 78 cortical areas spanning the entire cortical mantle of the rhesus macaque, we created a high dimensional set of anatomical features derived from captured images of cortical tissue stained for myelin and SMI-32. The method involved semi-automated de-noising of images, and enabled comparison of brain areas without hand-labeling of features such as layer boundaries. We applied multidimensional scaling (MDS) to the dataset to visualize similarity among cortical areas. This analysis shows a systematic variation between weakly laminated (limbic) cortices and sharply laminated (eulaminate) cortices. We call this smooth continuum the “cortical spectrum”. We also show that this spectrum is visible within subsystems of the cortex: the occipital, parietal, temporal, motor, prefrontal, and insular cortices. We compared the MDS-derived spectrum with a spectrum produced using T1- and T2-weighted magnetic resonance imaging (MRI) data derived from macaque, and found close agreement of the two coarse-graining methods. This suggests that T1w/T2w data, routinely obtained in human MRI studies, can serve as an effective proxy for data derived from high-resolution histological methods. More generally, this approach shows that the cortical spectrum is robust to the specific method used to compare cortical areas, and is therefore a powerful tool to understand the principles of organization of the primate cortex.

Empirical transmit field bias correction of T1w/T2w myelin maps

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

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.

Graded Variation in T1w/T2w Ratio during Adolescence: Measurement, Caveats, and Implications for Development of Cortical Myelin

Adolescence is characterized by the maturation of cortical microstructure and connectivity supporting complex cognition and behavior. Axonal myelination influences brain connectivity during development by enhancing neural signaling speed and inhibiting plasticity. However, the maturational timing of cortical myelination during human adolescence remains poorly understood. Here, we take advantage of recent advances in high-resolution cortical T1w/T2w mapping methods, including principled correction of B1+ transmit field effects, using data from the Human Connectome Project in Development (HCP-D; N = 628, ages 8-21). We characterize microstructural changes relevant to myelination by estimating age-related differences in T1w/T2w throughout the cerebral neocortex from childhood to early adulthood. We apply Bayesian spline models and clustering analysis to demonstrate graded variation in age-dependent cortical T1w/T2w differences that are correlated with the sensorimotor-association (S-A) axis of cortical organization reported by others. In sensorimotor areas, T1w/T2w ratio measures start at high levels at early ages, increase at a fast pace, and decelerate at later ages (18-21). In intermediate multimodal areas along the S-A axis, T1w/T2w starts at intermediate levels and increases linearly at an intermediate pace. In transmodal/paralimbic association areas, T1w/T2w starts at low levels and increases linearly at the slowest pace. These data provide evidence for graded variation of the T1w/T2w ratio along the S-A axis that may reflect cortical myelination changes during adolescence underlying the development of complex information processing and psychological functioning. We discuss the implications of these results as well as caveats in interpreting magnetic resonance imaging (MRI)-based estimates of myelination.SIGNIFICANCE STATEMENT Myelin is a lipid membrane that is essential to healthy brain function. Myelin wraps axons to increase neural signaling speed, enabling complex neuronal functioning underlying learning and cognition. Here, we characterize the developmental timing of myelination across the cerebral cortex during adolescence using a noninvasive proxy measure, T1w/T2w mapping. Our results provide new evidence demonstrating graded variation across the cortex in the timing of T1w/T2w changes during adolescence, with rapid T1w/T2w increases in lower-order sensory areas and gradual T1w/T2w increases in higher-order association areas. This spatial pattern of microstructural brain development closely parallels the sensorimotor-to-association axis of cortical organization and plasticity during ontogeny.

Adolescent development of multiscale structural wiring and functional interactions in the human connectome

Adolescence is a time of profound changes in the physical wiring and function of the brain. Here, we analyzed structural and functional brain network development in an accelerated longitudinal cohort spanning 14 to 25 y (n = 199). Core to our work was an advanced in vivo model of cortical wiring incorporating MRI features of corticocortical proximity, microstructural similarity, and white matter tractography. Longitudinal analyses assessing age-related changes in cortical wiring identified a continued differentiation of multiple corticocortical structural networks in youth. We then assessed structure-function coupling using resting-state functional MRI measures in the same participants both via cross-sectional analysis at baseline and by studying longitudinal change between baseline and follow-up scans. At baseline, regions with more similar structural wiring were more likely to be functionally coupled. Moreover, correlating longitudinal structural wiring changes with longitudinal functional connectivity reconfigurations, we found that increased structural differentiation, particularly between sensory/unimodal and default mode networks, was reflected by reduced functional interactions. These findings provide insights into adolescent development of human brain structure and function, illustrating how structural wiring interacts with the maturation of macroscale functional hierarchies.

Structural insight into the individual variability architecture of the functional brain connectome

Human cognition and behaviors depend upon the brain's functional connectomes, which vary remarkably across individuals. However, whether and how the functional connectome individual variability architecture is structurally constrained remains largely unknown. Using tractography- and morphometry-based network models, we observed the spatial convergence of structural and functional connectome individual variability, with higher variability in heteromodal association regions and lower variability in primary regions. We demonstrated that functional variability is significantly predicted by a unifying structural variability pattern and that this prediction follows a primary-to-heteromodal hierarchical axis, with higher accuracy in primary regions and lower accuracy in heteromodal regions. We further decomposed group-level connectome variability patterns into individual unique contributions and uncovered the structural-functional correspondence that is associated with individual cognitive traits. These results advance our understanding of the structural basis of individual functional variability and suggest the importance of integrating multimodal connectome signatures for individual differences in cognition and behaviors.

Modeling spatial, developmental, physiological, and topological constraints on human brain connectivity

The complex connectivity of nervous systems is thought to have been shaped by competitive selection pressures to minimize wiring costs and support adaptive function. Accordingly, recent modeling work indicates that stochastic processes, shaped by putative trade-offs between the cost and value of each connection, can successfully reproduce many topological properties of macroscale human connectomes measured with diffusion magnetic resonance imaging. Here, we derive a new formalism that more accurately captures the competing pressures of wiring cost minimization and topological complexity. We further show that model performance can be improved by accounting for developmental changes in brain geometry and associated wiring costs, and by using interregional transcriptional or microstructural similarity rather than topological wiring rules. However, all models struggled to capture topographical (i.e., spatial) network properties. Our findings highlight an important role for genetics in shaping macroscale brain connectivity and indicate that stochastic models offer an incomplete account of connectome organization.

A synergistic core for human brain evolution and cognition

How does the organization of neural information processing enable humans' sophisticated cognition? Here we decompose functional interactions between brain regions into synergistic and redundant components, revealing their distinct information-processing roles. Combining functional and structural neuroimaging with meta-analytic results, we demonstrate that redundant interactions are predominantly associated with structurally coupled, modular sensorimotor processing. Synergistic interactions instead support integrative processes and complex cognition across higher-order brain networks. The human brain leverages synergistic information to a greater extent than nonhuman primates, with high-synergy association cortices exhibiting the highest degree of evolutionary cortical expansion. Synaptic density mapping from positron emission tomography and convergent molecular and metabolic evidence demonstrate that synergistic interactions are supported by receptor diversity and human-accelerated genes underpinning synaptic function. This information-resolved approach provides analytic tools to disentangle information integration from coupling, enabling richer, more accurate interpretations of functional connectivity, and illuminating how the human neurocognitive architecture navigates the trade-off between robustness and integration.

Cytoarchitectonic Maps of the Human Metathalamus in 3D Space

The human metathalamus plays an important role in processing visual and auditory information. Understanding its layers and subdivisions is important to gain insights in its function as a subcortical relay station and involvement in various pathologies. Yet, detailed histological references of the microanatomy in 3D space are still missing. We therefore aim at providing cytoarchitectonic maps of the medial geniculate body (MGB) and its subdivisions in the BigBrain – a high-resolution 3D-reconstructed histological model of the human brain, as well as probabilistic cytoarchitectonic maps of the MGB and lateral geniculate body (LGB). Therefore, histological sections of ten postmortem brains were studied. Three MGB subdivisions (MGBv, MGBd, MGBm) were identified on every 5th BigBrain section, and a deep-learning based tool was applied to map them on every remaining section. The maps were 3D-reconstructed to show the shape and extent of the MGB and its subdivisions with cellular precision. The LGB and MGB were additionally identified in nine other postmortem brains. Probabilistic cytoarchitectonic maps in the MNI “Colin27” and MNI ICBM152 reference spaces were computed which reveal an overall low interindividual variability in topography and extent. The probabilistic maps were included into the Julich-Brain atlas, and are freely available. They can be linked to other 3D data of human brain organization and serve as an anatomical reference for diagnostic, prognostic and therapeutic neuroimaging studies of healthy brains and patients. Furthermore, the high-resolution MGB BigBrain maps provide a basis for data integration, brain modeling and simulation to bridge the larger scale involvement of thalamocortical and local subcortical circuits.

The highways and byways of the brain

Activity patterns and complex functions of the brain rely on the characteristic communication network formed by axonal fiber networks, but how many axons actually connect different brain regions? This Primer explores a study in PLOS Biology which finds that most areas of the human cerebral cortex are linked by an astoundingly small number of fibers.

Edges in brain networks: Contributions to models of structure and function

Network models describe the brain as sets of nodes and edges that represent its distributed organization. So far, most discoveries in network neuroscience have prioritized insights that highlight distinct groupings and specialized functional contributions of network nodes. Importantly, these functional contributions are determined and expressed by the web of their interrelationships, formed by network edges. Here, we underscore the important contributions made by brain network edges for understanding distributed brain organization. Different types of edges represent different types of relationships, including connectivity and similarity among nodes. Adopting a specific definition of edges can fundamentally alter how we analyze and interpret a brain network. Furthermore, edges can associate into collectives and higher order arrangements, describe time series, and form edge communities that provide insights into brain network topology complementary to the traditional node-centric perspective. Focusing on the edges, and the higher order or dynamic information they can provide, discloses previously underappreciated aspects of structural and functional network organization.

Impact of In Utero Exposure to Antiepileptic Drugs on Neonatal Brain Function

In utero brain development underpins brain health across the lifespan but is vulnerable to physiological and pharmacological perturbation. Here, we show that antiepileptic medication during pregnancy impacts on cortical activity during neonatal sleep, a potent indicator of newborn brain health. These effects are evident in frequency-specific functional brain networks and carry prognostic information for later neurodevelopment. Notably, such effects differ between different antiepileptic drugs that suggest neurodevelopmental adversity from exposure to antiepileptic drugs and not maternal epilepsy per se. This work provides translatable bedside metrics of brain health that are sensitive to the effects of antiepileptic drugs on postnatal neurodevelopment and carry direct prognostic value.

Connectivity alterations in autism reflect functional idiosyncrasy

Autism spectrum disorder (ASD) is commonly understood as an alteration of brain networks, yet case-control analyses against typically-developing controls (TD) have yielded inconsistent results. Here, we devised a novel approach to profile the inter-individual variability in functional network organization and tested whether such idiosyncrasy contributes to connectivity alterations in ASD. Studying a multi-centric dataset with 157 ASD and 172 TD, we obtained robust evidence for increased idiosyncrasy in ASD relative to TD in default mode, somatomotor and attention networks, but also reduced idiosyncrasy in lateral temporal cortices. Idiosyncrasy increased with age and significantly correlated with symptom severity in ASD. Furthermore, while patterns of functional idiosyncrasy were not correlated with ASD-related cortical thickness alterations, they co-localized with the expression patterns of ASD risk genes. Notably, we could demonstrate that patterns of atypical idiosyncrasy in ASD closely overlapped with connectivity alterations that are measurable with conventional case-control designs and may, thus, be a principal driver of inconsistency in the autism connectomics literature. These findings support important interactions between inter-individual heterogeneity in autism and functional signatures. Our findings provide novel biomarkers to study atypical brain development and may consolidate prior research findings on the variable nature of connectome level anomalies in autism.

Neurodevelopment of the association cortices: Patterns, mechanisms, and implications for psychopathology

The human brain undergoes a prolonged period of cortical development that spans multiple decades. During childhood and adolescence, cortical development progresses from lower-order, primary and unimodal cortices with sensory and motor functions to higher-order, transmodal association cortices subserving executive, socioemotional, and mentalizing functions. The spatiotemporal patterning of cortical maturation thus proceeds in a hierarchical manner, conforming to an evolutionarily rooted, sensorimotor-to-association axis of cortical organization. This developmental program has been characterized by data derived from multimodal human neuroimaging and is linked to the hierarchical unfolding of plasticity-related neurobiological events. Critically, this developmental program serves to enhance feature variation between lower-order and higher-order regions, thus endowing the brain's association cortices with unique functional properties. However, accumulating evidence suggests that protracted plasticity within late-maturing association cortices, which represents a defining feature of the human developmental program, also confers risk for diverse developmental psychopathologies.

Cortical patterning of morphometric similarity gradient reveals diverged hierarchical organization in sensory-motor cortices

The topological organization of the cerebral cortex provides hierarchical axes, namely gradients, which reveal systematic variations of brain structure and function. However, the hierarchical organization of macroscopic brain morphology and how it constrains cortical function along the organizing axes remains unclear. We map the gradient of cortical morphometric similarity (MS) connectome, combining multiple features conceptualized as a "fingerprint" of an individual's brain. The principal MS gradient is anchored by motor and sensory cortices at two extreme ends, which are reliable and reproducible. Notably, divergences between motor and sensory hierarchies are consistent with the laminar histological thickness gradient but contrary to the canonical functional connectivity (FC) gradient. Moreover, the MS dissociates with FC gradients in the higher-order association cortices. The MS gradient recapitulates fundamental properties of cortical organization, from gene expression and cyto- and myeloarchitecture to evolutionary expansion. Collectively, our findings provide a heuristic hierarchical organization of cortical morphological neuromarkers.

Comparing spatial null models for brain maps

Technological and data sharing advances have led to a proliferation of high-resolution structural and functional maps of the brain. Modern neuroimaging research increasingly depends on identifying correspondences between the topographies of these maps; however, most standard methods for statistical inference fail to account for their spatial properties. Recently, multiple methods have been developed to generate null distributions that preserve the spatial autocorrelation of brain maps and yield more accurate statistical estimates. Here, we comprehensively assess the performance of ten published null frameworks in statistical analyses of neuroimaging data. To test the efficacy of these frameworks in situations with a known ground truth, we first apply them to a series of controlled simulations and examine the impact of data resolution and spatial autocorrelation on their family-wise error rates. Next, we use each framework with two empirical neuroimaging datasets, investigating their performance when testing (1) the correspondence between brain maps (e.g., correlating two activation maps) and (2) the spatial distribution of a feature within a partition (e.g., quantifying the specificity of an activation map within an intrinsic functional network). Finally, we investigate how differences in the implementation of these null models may impact their performance. In agreement with previous reports, we find that naive null models that do not preserve spatial autocorrelation consistently yield elevated false positive rates and unrealistically liberal statistical estimates. While spatially-constrained null models yielded more realistic, conservative estimates, even these frameworks suffer from inflated false positive rates and variable performance across analyses. Throughout our results, we observe minimal impact of parcellation and resolution on null model performance. Altogether, our findings highlight the need for continued development of statistically-rigorous methods for comparing brain maps. The present report provides a harmonised framework for benchmarking and comparing future advancements.

Bridging the gap between single receptor type activity and whole-brain dynamics

What is the effect of activating a single modulatory neuronal receptor type on entire brain network dynamics? Can such effect be isolated at all? These are important questions because characterizing elementary neuronal processes that influence network activity across the given anatomical backbone is fundamental to guide theories of brain function. Here, we introduce the concept of the cortical 'receptome' taking into account the distribution and densities of expression of different modulatory receptor types across the brain's anatomical connectivity matrix. By modelling whole-brain dynamics in silico, we suggest a bidirectional coupling between modulatory neurotransmission and neuronal connectivity hardware exemplified by the impact of single serotonergic (5-HT) receptor types on cortical dynamics. As experimental support of this concept, we show how optogenetic tools enable specific activation of a single 5-HT receptor type across the cortex as well as in vivo measurement of its distinct effects on cortical processing. Altogether, we demonstrate how the structural neuronal connectivity backbone and its modulation by a single neurotransmitter system allow access to a rich repertoire of different brain states that are fundamental for flexible behaviour. We further propose that irregular receptor expression patterns-genetically predisposed or acquired during a lifetime-may predispose for neuropsychiatric disorders like addiction, depression and anxiety along with distinct changes in brain state. Our long-term vision is that such diseases could be treated through rationally targeted therapeutic interventions of high specificity to eventually recover natural transitions of brain states.

Topographic gradients of intrinsic dynamics across neocortex

The intrinsic dynamics of neuronal populations are shaped by both microscale attributes and macroscale connectome architecture. Here we comprehensively characterize the rich temporal patterns of neural activity throughout the human brain. Applying massive temporal feature extraction to regional haemodynamic activity, we systematically estimate over 6000 statistical properties of individual brain regions' time-series across the neocortex. We identify two robust spatial gradients of intrinsic dynamics, one spanning a ventromedial-dorsolateral axis and dominated by measures of signal autocorrelation, and the other spanning a unimodal-transmodal axis and dominated by measures of dynamic range. These gradients reflect spatial patterns of gene expression, intracortical myelin and cortical thickness, as well as structural and functional network embedding. Importantly, these gradients are correlated with patterns of meta-analytic functional activation, differentiating cognitive versus affective processing and sensory versus higher-order cognitive processing. Altogether, these findings demonstrate a link between microscale and macroscale architecture, intrinsic dynamics, and cognition.