A transcriptional signature of hub connectivity in the mouse connectome

Integration of estimated regional gene expression with neuroimaging and clinical phenotypes at biobank scale

An understanding of human brain individuality requires the integration of data on brain organization across people and brain regions, molecular and systems scales, as well as healthy and clinical states. Here, we help advance this understanding by leveraging methods from computational genomics to integrate large-scale genomic, transcriptomic, neuroimaging, and electronic-health record data sets. We estimated genetically regulated gene expression (gr-expression) of 18,647 genes, across 10 cortical and subcortical regions of 45,549 people from the UK Biobank. First, we showed that patterns of estimated gr-expression reflect known genetic-ancestry relationships, regional identities, as well as inter-regional correlation structure of directly assayed gene expression. Second, we performed transcriptome-wide association studies (TWAS) to discover 1,065 associations between individual variation in gr-expression and gray-matter volumes across people and brain regions. We benchmarked these associations against results from genome-wide association studies (GWAS) of the same sample and found hundreds of novel associations relative to these GWAS. Third, we integrated our results with clinical associations of gr-expression from the Vanderbilt Biobank. This integration allowed us to link genes, via gr-expression, to neuroimaging and clinical phenotypes. Fourth, we identified associations of polygenic gr-expression with structural and functional MRI phenotypes in the Human Connectome Project (HCP), a small neuroimaging-genomic data set with high-quality functional imaging data. Finally, we showed that estimates of gr-expression and magnitudes of TWAS were generally replicable and that the p-values of TWAS were replicable in large samples. Collectively, our results provide a powerful new resource for integrating gr-expression with population genetics of brain organization and disease.

The neocortical infrastructure for language involves region-specific patterns of laminar gene expression

The language network of the human brain has core components in the inferior frontal cortex and superior/middle temporal cortex, with left-hemisphere dominance in most people. Functional specialization and interconnectivity of these neocortical regions is likely to be reflected in their molecular and cellular profiles. Excitatory connections between cortical regions arise and innervate according to layer-specific patterns. Here, we generated a gene expression dataset from human postmortem cortical tissue samples from core language network regions, using spatial transcriptomics to discriminate gene expression across cortical layers. Integration of these data with existing single-cell expression data identified 56 genes that showed differences in laminar expression profiles between the frontal and temporal language cortex together with upregulation in layer II/III and/or layer V/VI excitatory neurons. Based on data from large-scale genome-wide screening in the population, DNA variants within these 56 genes showed set-level associations with interindividual variation in structural connectivity between the left-hemisphere frontal and temporal language cortex, and with the brain-related disorders dyslexia and schizophrenia which often involve affected language. These findings identify region-specific patterns of laminar gene expression as a feature of the brain's language network.

A network control theory pipeline for studying the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains the dynamics of a system. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter the dynamics of a system in a desired way. An interesting development for NCT in the neuroscience field is its application to study behavior and mental health symptoms. To date, NCT has been validated to study different aspects of the human structural connectome. NCT outputs can be monitored throughout developmental stages to study the effects of connectome topology on neural dynamics and, separately, to test the coherence of empirical datasets with brain function and stimulation. Here, we provide a comprehensive pipeline for applying NCT to structural connectomes by following two procedures. The main procedure focuses on computing the control energy associated with the transitions between specific neural activity states. The second procedure focuses on computing average controllability, which indexes nodes' general capacity to control the dynamics of the system. We provide recommendations for comparing NCT outputs against null network models, and we further support this approach with a Python-based software package called 'network control theory for python'. The procedures in this protocol are appropriate for users with a background in network neuroscience and experience in dynamical systems theory.

BrainTACO: an explorable multi-scale multi-modal brain transcriptomic and connectivity data resource

Exploring the relationships between genes and brain circuitry can be accelerated by joint analysis of heterogeneous datasets from 3D imaging data, anatomical data, as well as brain networks at varying scales, resolutions, and modalities. Generating an integrated view, beyond the individual resources' original purpose, requires the fusion of these data to a common space, and a visualization that bridges the gap across scales. However, despite ever expanding datasets, few platforms for integration and exploration of this heterogeneous data exist. To this end, we present the BrainTACO (Brain Transcriptomic And Connectivity Data) resource, a selection of heterogeneous, and multi-scale neurobiological data spatially mapped onto a common, hierarchical reference space, combined via a holistic data integration scheme. To access BrainTACO, we extended BrainTrawler, a web-based visual analytics framework for spatial neurobiological data, with comparative visualizations of multiple resources. This enables gene expression dissection of brain networks with, to the best of our knowledge, an unprecedented coverage and allows for the identification of potential genetic drivers of connectivity in both mice and humans that may contribute to the discovery of dysconnectivity phenotypes. Hence, BrainTACO reduces the need for time-consuming manual data aggregation often required for computational analyses in script-based toolboxes, and supports neuroscientists by directly leveraging the data instead of preparing it.

Disrupted cerebellar structural connectome in spinocerebellar ataxia type 3 and its association with transcriptional profiles

Spinocerebellar ataxia type 3 (SCA3) is primarily characterized by progressive cerebellar degeneration, including gray matter atrophy and disrupted anatomical and functional connectivity. The alterations of cerebellar white matter structural network in SCA3 and the underlying neurobiological mechanism remain unknown. Using a cohort of 20 patients with SCA3 and 20 healthy controls, we constructed cerebellar structural networks from diffusion MRI and investigated alterations of topological organization. Then, we mapped the alterations with transcriptome data from the Allen Human Brain Atlas to identify possible biological mechanisms for regional selective vulnerability to white matter damage. Compared with healthy controls, SCA3 patients exhibited reduced global and nodal efficiency, along with a widespread decrease in edge strength, particularly affecting edges connected to hub regions. The strength of inter-module connections was lower in SCA3 group and negatively correlated with the Scale for the Assessment and Rating of Ataxia score, International Cooperative Ataxia Rating Scale score, and cytosine-adenine-guanine repeat number. Moreover, the transcriptome-connectome association study identified the expression of genes involved in synapse-related and metabolic biological processes. These findings suggest a mechanism of white matter vulnerability and a potential image biomarker for the disease severity, providing insights into neurodegeneration and pathogenesis in this disease.

The mouse motor system contains multiple premotor areas and partially follows human organizational principles

While humans are known to have several premotor cortical areas, secondary motor cortex (M2) is often considered to be the only higher-order motor area of the mouse brain and is thought to combine properties of various human premotor cortices. Here, we show that axonal tracer, functional connectivity, myelin mapping, gene expression, and optogenetics data contradict this notion. Our analyses reveal three premotor areas in the mouse, anterior-lateral motor cortex (ALM), anterior-lateral M2 (aM2), and posterior-medial M2 (pM2), with distinct structural, functional, and behavioral properties. By using the same techniques across mice and humans, we show that ALM has strikingly similar functional and microstructural properties to human anterior ventral premotor areas and that aM2 and pM2 amalgamate properties of human pre-SMA and cingulate cortex. These results provide evidence for the existence of multiple premotor areas in the mouse and chart a comparative map between the motor systems of humans and mice.

Transcriptomic Similarity Informs Neuromorphic Deviations in Depression Biotypes

BACKGROUND

Major depressive disorder (MDD) is complicated by population heterogeneity, motivating the investigation of biotypes through imaging-derived phenotypes. However, neuromorphic heterogeneity in MDD remains unclear, and how the correlated gene expression (CGE) connectome constrains these neuromorphic anomalies in MDD biotypes has not yet been studied.

METHODS

Here, we related cortical thickness deviations in MDD biotypes to a pattern of CGE connectome. Cortical thickness was estimated from 3-dimensional T1-weighted magnetic resonance images in 2 independent cohorts (discovery cohort: N = 425; replication cohort: N = 217). The transcriptional activity was measured according to Allen Human Brain Atlas. A density peak-based clustering algorithm was used to identify MDD biotypes.

RESULTS

We found that patients with MDD were clustered into 2 replicated biotypes based on single-patient regional deviations from healthy control participants across 2 datasets. Biotype 1 mainly exhibited cortical thinning across the brain, whereas biotype 2 mainly showed cortical thickening in the brain. Using brainwide gene expression data, we found that deviations of transcriptionally connected neighbors predicted regional deviation for both biotypes. Furthermore, putative CGE-informed epicenters of biotype 1 were concentrated on the cognitive control circuit, whereas biotype 2 epicenters were located in the social perception circuit. The patterns of epicenter likelihood were separately associated with depression- and anxiety-response maps, suggesting that epicenters of MDD biotypes may be associated with clinical efficacies.

CONCLUSIONS

Our findings linked the CGE connectome and neuromorphic deviations to identify distinct epicenters in MDD biotypes, providing insight into how microscale gene expressions informed MDD biotypes.

Antagonists of the stress and opioid systems restore the functional connectivity of the prefrontal cortex during alcohol withdrawal through divergent mechanisms

Chronic alcohol consumption leads to dependence and withdrawal symptoms upon cessation, contributing to persistent use. However, the brain network mechanisms by which the brain orchestrates alcohol withdrawal and how these networks are affected by pharmacological treatments remain elusive. Recent work revealed that alcohol withdrawal produces a widespread increase in coordinated brain activity and a decrease in modularity of the whole-brain functional network using single-cell whole-brain imaging of immediate early genes. This decreased modularity and functional hyperconnectivity are hypothesized to be novel biomarkers of alcohol withdrawal in alcohol dependence, which could potentially be used to evaluate the efficacy of new medications for alcohol use disorder. However, there is no evidence that current FDA-approved medications or experimental treatments known to reduce alcohol drinking in animal models can normalize the changes in whole-brain functional connectivity. In this report, we tested the effect of R121919, a CRF1 antagonist, and naltrexone, an FDA-approved treatment for alcohol use disorder, on whole-brain functional connectivity using the cellular marker FOS combined with graph theory and advanced network analyses. Results show that both R121919 and naltrexone restored the functional connectivity of the prefrontal cortex during alcohol withdrawal, but through divergent mechanisms. Specifically, R121919 increased FOS activation in the prefrontal cortex, partially restored modularity, and normalized connectivity, particularly in CRF1-rich regions, including the prefrontal, pallidum, and extended amygdala circuits. On the other hand, naltrexone decreased FOS activation throughout the brain, decreased modularity, and increased connectivity overall except for the Mu opioid receptor-rich regions, including the thalamus. These results identify the brain networks underlying the pharmacological effects of R121919 and naltrexone and demonstrate that these drugs restored different aspects of functional connectivity of the prefrontal cortex, pallidum, amygdala, and thalamus during alcohol withdrawal. Notably, these effects were particularly prominent in CRF1- and Mu opioid receptors-rich regions highlighting the potential of whole-brain functional connectivity using FOS as a tool for identifying neuronal network mechanisms underlying the pharmacological effects of existing and new medications for alcohol use disorder.

Homeostasis to Allostasis: Prefrontal Astrocyte Roles in Cognitive Flexibility and Stress Biology

In the intricate landscape of neurophysiology, astrocytes have been traditionally cast as homeostatic cells; however, their mechanistic involvement in allostasis-particularly how they modulate the adaptive response to stress and its accumulative impact that disrupts cognitive functions and precipitates psychiatric disorders-is now starting to be unraveled. Here, we address the gap by positing astrocytes as crucial allostatic players whose molecular adaptations underlie cognitive flexibility in stress-related neuropsychiatric conditions. We review how astrocytes, responding to stress mediators such as glucocorticoid and epinephrine/norepinephrine, undergo morphological and functional transformations that parallel the maladaptive changes. Our synthesis of recent findings reveals that these glial changes, especially in the metabolically demanding prefrontal cortex, may underlie some of the neuropsychiatric mechanisms characterized by the disruption of energy metabolism and astrocytic networks, compromised glutamate clearance, and diminished synaptic support. We argue that astrocytes extend beyond their homeostatic role, actively participating in the brain's allostatic response, especially by modulating energy substrates critical for cognitive functions.

Can hubs of the human connectome be identified consistently with diffusion MRI?

Recent years have seen a surge in the use of diffusion MRI to map connectomes in humans, paralleled by a similar increase in processing and analysis choices. Yet these different steps and their effects are rarely compared systematically. Here, in a healthy young adult population (n = 294), we characterized the impact of a range of analysis pipelines on one widely studied property of the human connectome: its degree distribution. We evaluated the effects of 40 pipelines (comparing common choices of parcellation, streamline seeding, tractography algorithm, and streamline propagation constraint) and 44 group-representative connectome reconstruction schemes on highly connected hub regions. We found that hub location is highly variable between pipelines. The choice of parcellation has a major influence on hub architecture, and hub connectivity is highly correlated with regional surface area in most of the assessed pipelines (ρ > 0.70 in 69% of the pipelines), particularly when using weighted networks. Overall, our results demonstrate the need for prudent decision-making when processing diffusion MRI data, and for carefully considering how different processing choices can influence connectome organization.

Towards a biologically annotated brain connectome

The brain is a network of interleaved neural circuits. In modern connectomics, brain connectivity is typically encoded as a network of nodes and edges, abstracting away the rich biological detail of local neuronal populations. Yet biological annotations for network nodes - such as gene expression, cytoarchitecture, neurotransmitter receptors or intrinsic dynamics - can be readily measured and overlaid on network models. Here we review how connectomes can be represented and analysed as annotated networks. Annotated connectomes allow us to reconceptualize architectural features of networks and to relate the connection patterns of brain regions to their underlying biology. Emerging work demonstrates that annotated connectomes help to make more veridical models of brain network formation, neural dynamics and disease propagation. Finally, annotations can be used to infer entirely new inter-regional relationships and to construct new types of network that complement existing connectome representations. In summary, biologically annotated connectomes offer a compelling way to study neural wiring in concert with local biological features.

Functional neuroimaging as a catalyst for integrated neuroscience

Functional magnetic resonance imaging (fMRI) enables non-invasive access to the awake, behaving human brain. By tracking whole-brain signals across a diverse range of cognitive and behavioural states or mapping differences associated with specific traits or clinical conditions, fMRI has advanced our understanding of brain function and its links to both normal and atypical behaviour. Despite this headway, progress in human cognitive neuroscience that uses fMRI has been relatively isolated from rapid advances in other subdomains of neuroscience, which themselves are also somewhat siloed from one another. In this Perspective, we argue that fMRI is well-placed to integrate the diverse subfields of systems, cognitive, computational and clinical neuroscience. We first summarize the strengths and weaknesses of fMRI as an imaging tool, then highlight examples of studies that have successfully used fMRI in each subdomain of neuroscience. We then provide a roadmap for the future advances that will be needed to realize this integrative vision. In this way, we hope to demonstrate how fMRI can help usher in a new era of interdisciplinary coherence in neuroscience.

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.

Spatial cell-type enrichment predicts mouse brain connectivity

A fundamental neuroscience topic is the link between the brain's molecular, cellular, and cytoarchitectonic properties and structural connectivity. Recent studies relate inter-regional connectivity to gene expression, but the relationship to regional cell-type distributions remains understudied. Here, we utilize whole-brain mapping of neuronal and non-neuronal subtypes via the matrix inversion and subset selection algorithm to model inter-regional connectivity as a function of regional cell-type composition with machine learning. We deployed random forest algorithms for predicting connectivity from cell-type densities, demonstrating surprisingly strong prediction accuracy of cell types in general, and particular non-neuronal cells such as oligodendrocytes. We found evidence of a strong distance dependency in the cell connectivity relationship, with layer-specific excitatory neurons contributing the most for long-range connectivity, while vascular and astroglia were salient for short-range connections. Our results demonstrate a link between cell types and connectivity, providing a roadmap for examining this relationship in other species, including humans.

Network connectivity and local transcriptomic vulnerability underpin cortical atrophy progression in Parkinson's disease

Parkinson's disease pathology is hypothesized to spread through the brain via axonal connections between regions and is further modulated by local vulnerabilities within those regions. The resulting changes to brain morphology have previously been demonstrated in both prodromal and de novo Parkinson's disease patients. However, it remains unclear whether the pattern of atrophy progression in Parkinson's disease over time is similarly explained by network-based spreading and local vulnerability. We address this gap by mapping the trajectory of cortical atrophy rates in a large, multi-centre cohort of Parkinson's disease patients and relate this atrophy progression pattern to network architecture and gene expression profiles. Across 4-year follow-up visits, increased atrophy rates were observed in posterior, temporal, and superior frontal cortices. We demonstrated that this progression pattern was shaped by network connectivity. Regional atrophy rates were strongly related to atrophy rates across structurally and functionally connected regions. We also found that atrophy progression was associated with specific gene expression profiles. The genes whose spatial distribution in the brain was most related to atrophy rate were those enriched for mitochondrial and metabolic function. Taken together, our findings demonstrate that both global and local brain features influence vulnerability to neurodegeneration in Parkinson's disease.

Predicting longitudinal brain atrophy in Parkinson's disease using a Susceptible-Infected-Removed agent-based model

Parkinson's disease is a progressive neurodegenerative disorder characterized by accumulation of abnormal isoforms of alpha-synuclein. Alpha-synuclein is proposed to act as a prion in Parkinson's disease: In its misfolded pathologic state, it favors the misfolding of normal alpha-synuclein molecules, spreads trans-neuronally, and causes neuronal damage as it accumulates. This theory remains controversial. We have previously developed a Susceptible-Infected-Removed (SIR) computational model that simulates the templating, propagation, and toxicity of alpha-synuclein molecules in the brain. In this study, we test this model with longitudinal MRI collected over 4 years from the Parkinson's Progression Markers Initiative (1,068 T1 MRI scans, 790 Parkinson's disease scans, and 278 matched control scans). We find that brain deformation progresses in subcortical and cortical regions. The SIR model recapitulates the spatiotemporal distribution of brain atrophy observed in Parkinson's disease. We show that connectome topology and geometry significantly contribute to model fit. We also show that the spatial expression of two genes implicated in alpha-synuclein synthesis and clearance, SNCA and GBA, also influences the atrophy pattern. We conclude that the progression of atrophy in Parkinson's disease is consistent with the prion-like hypothesis and that the SIR model is a promising tool to investigate multifactorial neurodegenerative diseases over time.

A phylogenetically-conserved axis of thalamocortical connectivity in the human brain

The thalamus enables key sensory, motor, emotive, and cognitive processes via connections to the cortex. These projection patterns are traditionally considered to originate from discrete thalamic nuclei, however recent work showing gradients of molecular and connectivity features in the thalamus suggests the organisation of thalamocortical connections occurs along a continuous dimension. By performing a joint decomposition of densely sampled gene expression and non-invasive diffusion tractography in the adult human thalamus, we define a principal axis of genetic and connectomic variation along a medial-lateral thalamic gradient. Projections along this axis correspond to an anterior-posterior cortical pattern and are aligned with electrophysiological properties of the cortex. The medial-lateral axis demonstrates phylogenetic conservation, reflects transitions in neuronal subtypes, and shows associations with neurodevelopment and common brain disorders. This study provides evidence for a supra-nuclear axis of thalamocortical organisation characterised by a graded transition in molecular properties and anatomical connectivity.

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.

Using network control theory to study the dynamics of the structural connectome

Network control theory (NCT) is a simple and powerful tool for studying how network topology informs and constrains dynamics. Compared to other structure-function coupling approaches, the strength of NCT lies in its capacity to predict the patterns of external control signals that may alter dynamics in a desired way. We have extensively developed and validated the application of NCT to the human structural connectome. Through these efforts, we have studied (i) how different aspects of connectome topology affect neural dynamics, (ii) whether NCT outputs cohere with empirical data on brain function and stimulation, and (iii) how NCT outputs vary across development and correlate with behavior and mental health symptoms. In this protocol, we introduce a framework for applying NCT to structural connectomes following two main pathways. Our primary pathway focuses on computing the control energy associated with transitioning between specific neural activity states. Our second pathway focuses on computing average controllability, which indexes nodes' general capacity to control dynamics. We also provide recommendations for comparing NCT outputs against null network models. Finally, we support this protocol with a Python-based software package called network control theory for python (nctpy).

Multimodal and multiscale evidence for network-based cortical thinning in major depressive disorder

BACKGROUND

Major depressive disorder (MDD) is associated with widespread, irregular cortical thickness (CT) reductions across the brain. However, little is known regarding mechanisms that govern spatial distribution of the reductions.

METHODS

We combined multimodal MRI and genetic, cytoarchitectonic and chemoarchitectonic data to examine structural covariance, functional synchronization, gene co-expression, cytoarchitectonic similarity and chemoarchitectonic covariance between regions atrophied in MDD.

RESULTS

Regions atrophied in MDD were associated with significantly higher structural covariance, functional synchronization, gene co-expression and chemoarchitectonic covariance. These results were robust against methodological variations in brain parcellation and null model, reproducible in patients and controls, and independent of age at onset of MDD. Despite no significant differences in the cytoarchitectonic similarity, MDD-related CT reductions were susceptible to specific cytoarchitectonic class of association cortex. Further, we found that nodal shortest path lengths to disease epicenters derived from structural (right supramarginal gyrus) and chemoarchitectonic covariance (right sulcus intermedius primus) networks of healthy brains were correlated with the extent to which a region was atrophied in MDD, supporting the transneuronal spread hypothesis that regions closer to the epicenters are more susceptible to MDD. Finally, we showed that structural covariance and functional synchronization among regions atrophied in MDD were mainly related to genes enriched in metabolic and membrane-related processes, driven by genes in excitatory neurons, and associated with specific neurotransmitter transporters and receptors.

CONCLUSIONS

Altogether, our findings provide empirical evidence for and genetic and molecular insights into connectivity-constrained CT thinning in MDD.

Gradients of neurotransmitter receptor expression in the macaque cortex

Dynamics and functions of neural circuits depend on interactions mediated by receptors. Therefore, a comprehensive map of receptor organization across cortical regions is needed. In this study, we used in vitro receptor autoradiography to measure the density of 14 neurotransmitter receptor types in 109 areas of macaque cortex. We integrated the receptor data with anatomical, genetic and functional connectivity data into a common cortical space. We uncovered a principal gradient of receptor expression per neuron. This aligns with the cortical hierarchy from sensory cortex to higher cognitive areas. A second gradient, driven by serotonin 5-HT1A receptors, peaks in the anterior cingulate, default mode and salience networks. We found a similar pattern of 5-HT1A expression in the human brain. Thus, the macaque may be a promising translational model of serotonergic processing and disorders. The receptor gradients may enable rapid, reliable information processing in sensory cortical areas and slow, flexible integration in higher cognitive areas.

Hyperconnectivity of Two Separate Long-Range Cholinergic Systems Contributes to the Reorganization of the Brain Functional Connectivity during Nicotine Withdrawal in Male Mice

Chronic nicotine results in dependence with withdrawal symptoms on discontinuation of use, through desensitization of nicotinic acetylcholine receptors and altered cholinergic neurotransmission. Nicotine withdrawal is associated with increased whole-brain functional connectivity and decreased network modularity; however, the role of cholinergic neurons in those changes is unknown. To identify the contribution of nicotinic receptors and cholinergic regions to changes in the functional network, we analyzed the contribution of the main cholinergic regions to brain-wide activation of the immediate early-gene Fos during withdrawal in male mice and correlated these changes with the expression of nicotinic receptor mRNA throughout the brain. We show that the main functional connectivity modules included the main long-range cholinergic regions, which were highly synchronized with the rest of the brain. However, despite this hyperconnectivity, they were organized into two anticorrelated networks that were separated into basal forebrain-projecting and brainstem-thalamic-projecting cholinergic regions, validating a long-standing hypothesis of the organization of the brain cholinergic systems. Moreover, baseline (without nicotine) expression of Chrna2, Chrna3, Chrna10, and Chrnd mRNA of each brain region correlated with withdrawal-induced changes in Fos expression. Finally, by mining the Allen Brain mRNA expression database, we were able to identify 1755 gene candidates and three pathways (Sox2-Oct4-Nanog, JAK-STAT, and MeCP2-GABA) that may contribute to nicotine withdrawal-induced Fos expression. These results identify the dual contribution of the basal forebrain and brainstem-thalamic cholinergic systems to whole-brain functional connectivity during withdrawal; and identify nicotinic receptors and novel cellular pathways that may be critical for the transition to nicotine dependence.

Assortative mixing in micro-architecturally annotated brain connectomes

The wiring of the brain connects micro-architecturally diverse neuronal populations, but the conventional graph model, which encodes macroscale brain connectivity as a network of nodes and edges, abstracts away the rich biological detail of each regional node. Here, we annotate connectomes with multiple biological attributes and formally study assortative mixing in annotated connectomes. Namely, we quantify the tendency for regions to be connected based on the similarity of their micro-architectural attributes. We perform all experiments using four cortico-cortical connectome datasets from three different species, and consider a range of molecular, cellular, and laminar annotations. We show that mixing between micro-architecturally diverse neuronal populations is supported by long-distance connections and find that the arrangement of connections with respect to biological annotations is associated to patterns of regional functional specialization. By bridging scales of cortical organization, from microscale attributes to macroscale connectivity, this work lays the foundation for next-generation annotated connectomics.

A whole new world: embracing the systems-level to understand the indirect impact of pathology in neurodegenerative disorders

The direct link between neuropathology and the symptoms that emerge from damage to the brain is often difficult to discern. In this perspective, we argue that a satisfying account of neurodegenerative symptoms most naturally emerges from the consideration of the brain from the systems-level. Specifically, we will highlight the role of the neuromodulatory arousal system, which is uniquely positioned to coordinate the brain's ability to flexibly integrate the otherwise segregated structures required to support higher cognitive functions. Importantly, the neuromodulatory arousal system is highly heterogeneous, encompassing structures that are common sites of neurodegeneration across Alzheimer's and Parkinson's disease. We will review studies that implicate the dysfunctional interactions amongst distributed brain regions as a side-effect of pathological involvement of the neuromodulatory arousal system in these neurodegenerative disorders. From this perspective, we will argue that future work in clinical neuroscience should attempt to consider the inherent complexity in the brain and employ analytic techniques that do not solely focus on regional functional impairments, but rather captures the brain as an inherently dynamic, distributed, multi-scale system. Through this lens, we hope that we will devise new and improved diagnostic markers and interventional approaches to aid in the treatment of neurodegenerative disorders.

Hyperconnectivity of two separate long-range cholinergic systems contributes to the reorganization of the brain functional connectivity during nicotine withdrawal in male mice

Chronic nicotine results in dependence with withdrawal symptoms upon discontinuation of use, through desensitization of nicotinic acetylcholine receptors and altered cholinergic neurotransmission. Nicotine withdrawal is associated with increased whole-brain functional connectivity and decreased network modularity, however, the role of cholinergic neurons in those changes is unknown. To identify the contribution of nicotinic receptors and cholinergic regions to changes in the functional network, we analyzed the contribution of the main cholinergic regions to brain-wide activation of the immediate early-gene FOS during withdrawal in male mice and correlated these changes with the expression of nicotinic receptor mRNA throughout the brain. We show that the main functional connectivity modules included the main long-range cholinergic regions, which were highly synchronized with the rest of the brain. However, despite this hyperconnectivity they were organized into two anticorrelated networks that were separated into basal forebrain projecting and brainstem-thalamic projecting cholinergic regions, validating a long-standing hypothesis of the organization of the brain cholinergic systems. Moreover, baseline (without nicotine) expression of Chrna2 , Chrna3 , Chrna10 , and Chrnd mRNA of each brain region correlated with withdrawal-induced changes in FOS expression. Finally, by mining the Allen Brain mRNA expression database, we were able to identify 1755 gene candidates and three pathways (Sox2-Oct4-Nanog, JAK-STAT, and MeCP2-GABA) that may contribute to nicotine withdrawal-induced FOS expression. These results identify the dual contribution of the basal forebrain and brainstem-thalamic cholinergic systems to whole-brain functional connectivity during withdrawal; and identify nicotinic receptors and novel cellular pathways that may be critical for the transition to nicotine dependence.

Significance Statement

Discontinuation of nicotine use in dependent users is associated with increased whole-brain activation and functional connectivity and leads to withdrawal symptoms. Here we investigated the contribution of the nicotinic cholinergic receptors and main cholinergic projecting brain areas in the whole-brain changes associated with withdrawal. This not only allowed us to visualize and confirm the previously described duality of the cholinergic brain system using this novel methodology, but also identify nicotinic receptors together with 1751 other genes that contribute, and could thus be targets for treatments against, nicotine withdrawal and dependence.

Altered global signal topography in Alzheimer's disease

BACKGROUND

Alzheimer's disease (AD) is a neurodegenerative disease associated with widespread disruptions in intrinsic local specialization and global integration in the functional system of the brain. These changes in integration may further disrupt the global signal (GS) distribution, which might represent the local relative contribution to global activity in functional magnetic resonance imaging (fMRI).

METHODS

fMRI scans from a discovery dataset (n = 809) and a validated dataset (n = 542) were used in the analysis. We investigated the alteration of GS topography using the GS correlation (GSCORR) in patients with mild cognitive impairment (MCI) and AD. The association between GS alterations and functional network properties was also investigated based on network theory. The underlying mechanism of GSCORR alterations was elucidated using imaging-transcriptomics.

FINDINGS

Significantly increased GS topography in the frontal lobe and decreased GS topography in the hippocampus, cingulate gyrus, caudate, and middle temporal gyrus were observed in patients with AD (P_adj < 0.05). Notably, topographical GS changes in these regions correlated with cognitive ability (P < 0.05). The changes in GS topography also correlated with the changes in functional network segregation (ρ = 0.5). Moreover, the genes identified based on GS topographical changes were enriched in pathways associated with AD and neurodegenerative diseases.

INTERPRETATION

Our findings revealed significant changes in GS topography and its molecular basis, confirming the informative role of GS in AD and further contributing to the understanding of the relationship between global and local neuronal activities in patients with AD.

FUNDING

Beijing Natural Science Funds for Distinguished Young Scholars, China; Fundamental Research Funds for the Central Universities, China; National Natural Science Foundation, China.

Toward Best Practices for Imaging Transcriptomics of the Human Brain

Modern brainwide transcriptional atlases provide unprecedented opportunities for investigating the molecular correlates of brain organization, as quantified using noninvasive neuroimaging. However, integrating neuroimaging data with transcriptomic measures is not straightforward, and careful consideration is required to make valid inferences. In this article, we review recent work exploring how various methodological choices affect 3 main phases of imaging transcriptomic analyses, including 1) processing of transcriptional atlas data; 2) relating transcriptional measures to independently derived neuroimaging phenotypes; and 3) evaluating the functional implications of identified associations through gene enrichment analyses. Our aim is to facilitate the development of standardized and reproducible approaches for this rapidly growing field. We identify sources of methodological variability, key choices that can affect findings, and considerations for mitigating false positive and/or spurious results. Finally, we provide an overview of freely available open-source toolboxes implementing current best-practice procedures across all 3 analysis phases.

Long-range functional connections mirror and link microarchitectural and cognitive hierarchies in the human brain

BACKGROUND

Higher-order cognition is hypothesized to be implemented via distributed cortical networks that are linked via long-range connections. However, it is unknown how computational advantages of long-range connections reflect cortical microstructure and microcircuitry.

METHODS

We investigated this question by (i) profiling long-range cortical connectivity using resting-state functional magnetic resonance imaging (MRI) and cortico-cortical geodesic distance mapping, (ii) assessing how long-range connections reflect local brain microarchitecture, and (iii) examining the microarchitectural similarity of regions connected through long-range connections.

RESULTS

Analysis of 2 independent datasets indicated that sensory/motor areas had more clustered short-range connections, while transmodal association systems hosted distributed, long-range connections. Meta-analytical decoding suggested that this topographical difference mirrored shifts in cognitive function, from perception/action towards emotional/social processing. Analysis of myelin-sensitive in vivo MRI as well as postmortem histology and transcriptomics datasets established that gradients in functional connectivity distance are paralleled by those present in cortical microarchitecture. Notably, long-range connections were found to link spatially remote regions of association cortex with an unexpectedly similar microarchitecture.

CONCLUSIONS

By mapping covarying topographies of long-range functional connections and cortical microcircuits, the current work provides insights into structure-function relations in human neocortex.

Shared and Distinct Brain Regions Targeted for Immediate Early Gene Expression by Ketamine and Psilocybin

Psilocybin is a psychedelic with therapeutic potential. While there is growing evidence that psilocybin exerts its beneficial effects through enhancing neural plasticity, the exact brain regions involved are not completely understood. Determining the impact of psilocybin on plasticity-related gene expression throughout the brain can broaden our understanding of the neural circuits involved in psychedelic-evoked neural plasticity. In this study, whole-brain serial two-photon microscopy and light sheet microscopy were employed to map the expression of the immediate early gene, c-Fos, in male and female mice. The drug-induced c-Fos expression following psilocybin administration was compared to that of subanesthetic ketamine and saline control. Psilocybin and ketamine produced acutely comparable elevations in c-Fos expression in numerous brain regions, including anterior cingulate cortex, locus coeruleus, primary visual cortex, central and basolateral amygdala, medial and lateral habenula, and claustrum. Select regions exhibited drug-preferential differences, such as dorsal raphe and insular cortex for psilocybin and the CA1 subfield of hippocampus for ketamine. To gain insights into the contributions of receptors and cell types, the c-Fos expression maps were related to brain-wide in situ hybridization data. The transcript analyses showed that the endogenous levels of Grin2a and Grin2b predict whether a cortical region is sensitive to drug-evoked neural plasticity for both ketamine and psilocybin. Collectively, the systematic mapping approach produced an unbiased list of brain regions impacted by psilocybin and ketamine. The data are a resource that highlights previously underappreciated regions for future investigations. Furthermore, the robust relationships between drug-evoked c-Fos expression and endogenous transcript distributions suggest glutamatergic receptors as a potential convergent target for how psilocybin and ketamine produce their rapid-acting and long-lasting therapeutic effects.

Gene Expressions Preferentially Influence Cortical Thickness of Human Connectome Project Atlas Parcellated Regions in First-Episode Antipsychotic-Naïve Psychoses

Altered gene expressions may mechanistically link genetic factors with brain morphometric alterations. Existing gene expression studies have examined selected morphometric features using low-resolution atlases in medicated schizophrenia. We examined the relationship of gene expression with cortical thickness (CT), surface area (SA), and gray matter volume (GMV) of first-episode antipsychotic-naïve psychosis patients (FEAP = 85) and 81 controls, hypothesizing that gene expressions often associated with psychosis will differentially associate with different morphometric features. We explored such associations among schizophrenia and non-schizophrenia subgroups within FEAP group compared to controls. We mapped 360 Human Connectome Project atlas-based parcellations on brain MRI on to the publicly available brain gene expression data from the Allen Brain Institute collection. Significantly correlated genes were investigated using ingenuity pathway analysis to elucidate molecular pathways. CT but not SA or GMV correlated with expression of 1137 out of 15 633 genes examined controlling for age, sex, and average CT. Among these ≈19%, ≈39%, and 8% of genes were unique to FEAP, schizophrenia, and non-schizophrenia, respectively. Variants of 10 among these 1137 correlated genes previously showed genome-wide-association with schizophrenia. Molecular pathways associated with CT were axonal guidance and sphingosine pathways (common to FEAP and controls), selected inflammation pathways (unique to FEAP), synaptic modulation (unique to schizophrenia), and telomere extension (common to NSZ and healthy controls). We demonstrate that different sets of genes and molecular pathways may preferentially influence CT in different diagnostic groups. Genes with altered expressions correlating with CT and associated pathways may be targets for pathophysiological investigations and novel treatment designs.

Correspondence between gene expression and neurotransmitter receptor and transporter density in the human brain

Neurotransmitter receptors modulate signaling between neurons. Thus, neurotransmitter receptors and transporters play a key role in shaping brain function. Due to the lack of comprehensive neurotransmitter receptor/transporter density datasets, microarray gene expression measuring mRNA transcripts is often used as a proxy for receptor densities. In the present report, we comprehensively test the spatial correlation between gene expression and protein density for a total of 27 neurotransmitter receptors, receptor binding-sites, and transporters across 9 different neurotransmitter systems, using both PET and autoradiography radioligand-based imaging modalities. We find poor spatial correspondences between gene expression and density for all neurotransmitter receptors and transporters except four single-protein metabotropic receptors (5-HT_1A, CB₁, D₂, and MOR). These expression-density associations are related to gene differential stability and can vary between cortical and subcortical structures. Altogether, we recommend using direct measures of receptor and transporter density when relating neurotransmitter systems to brain structure and function.

Alteration of a brain network with stable and strong functional connections in subjects with schizophrenia

It is widely accepted that there are some common network patterns in the human brain. However, the existence of stable and strong functional connections in the human brain and whether they change in schizophrenia is still a question. By setting 1% connections with the smallest coefficient of variation, we found a widespread brain functional network (frame network) in healthy people(n = 380, two datasets from public databases). We then explored the alterations in a medicated group (60 subjects with schizophrenia vs 71 matched controls) and a drug-naive first-episode group (68 subjects with schizophrenia vs 45 matched controls). A linear support vector classifier (SVC) was constructed to distinguish patients and controls using the medicated patients' frame network. We found most frame connections of healthy people had high strength, which were symmetrical and connected the left and right hemispheres. Conversely, significant differences in frame connections were observed in both patient groups, which were positively correlated with negative symptoms (mainly language dysfunction). Additionally, patients' frame network were more left-lateralized, concentrating on the left frontal lobe, and was quite accurate at distinguishing medicated patients from controls (classifier accuracy was 78.63%, sensitivity was 86.67%, specificity was 76.06%, and the area under the curve (AUC) was 0.83). Furthermore, the results were repeated in the drug-naive set (accuracy was 84.96%, sensitivity was 85.29%, specificity was 88.89%, and AUC was 0.93). These findings indicate that the abnormal pattern of frame network in subjects with schizophrenia might provide new insights into the dysconnectivity in schizophrenia.

Neural Mechanisms and Psychology of Psychedelic Ego Dissolution

Neuroimaging studies of psychedelics have advanced our understanding of hierarchical brain organization and the mechanisms underlying their subjective and therapeutic effects. The primary mechanism of action of classic psychedelics is binding to serotonergic 5-HT2A receptors. Agonist activity at these receptors leads to neuromodulatory changes in synaptic efficacy that can have a profound effect on hierarchical message-passing in the brain. Here, we review the cognitive and neuroimaging evidence for the effects of psychedelics: in particular, their influence on selfhood and subject-object boundaries-known as ego dissolution-surmised to underwrite their subjective and therapeutic effects. Agonism of 5-HT2A receptors, located at the apex of the cortical hierarchy, may have a particularly powerful effect on sentience and consciousness. These effects can endure well after the pharmacological half-life, suggesting that psychedelics may have effects on neural plasticity that may play a role in their therapeutic efficacy. Psychologically, this may be accompanied by a disarming of ego resistance that increases the repertoire of perceptual hypotheses and affords alternate pathways for thought and behavior, including those that undergird selfhood. We consider the interaction between serotonergic neuromodulation and sentience through the lens of hierarchical predictive coding, which speaks to the value of psychedelics in understanding how we make sense of the world and specific predictions about effective connectivity in cortical hierarchies that can be tested using functional neuroimaging. SIGNIFICANCE STATEMENT: Classic psychedelics bind to serotonergic 5-HT2A receptors. Their agonist activity at these receptors leads to neuromodulatory changes in synaptic efficacy, resulting in a profound effect on information processing in the brain. Here, we synthesize an abundance of brain imaging research with pharmacological and psychological interpretations informed by the framework of predictive coding. Moreover, predictive coding is suggested to offer more sophisticated interpretations of neuroimaging findings by bridging the role between the 5-HT2A receptors and large-scale brain networks.

Local molecular and global connectomic contributions to cross-disorder cortical abnormalities

Numerous brain disorders demonstrate structural brain abnormalities, which are thought to arise from molecular perturbations or connectome miswiring. The unique and shared contributions of these molecular and connectomic vulnerabilities to brain disorders remain unknown, and has yet to be studied in a single multi-disorder framework. Using MRI morphometry from the ENIGMA consortium, we construct maps of cortical abnormalities for thirteen neurodevelopmental, neurological, and psychiatric disorders from N = 21,000 participants and N = 26,000 controls, collected using a harmonised processing protocol. We systematically compare cortical maps to multiple micro-architectural measures, including gene expression, neurotransmitter density, metabolism, and myelination (molecular vulnerability), as well as global connectomic measures including number of connections, centrality, and connection diversity (connectomic vulnerability). We find a relationship between molecular vulnerability and white-matter architecture that drives cortical disorder profiles. Local attributes, particularly neurotransmitter receptor profiles, constitute the best predictors of both disorder-specific cortical morphology and cross-disorder similarity. Finally, we find that cross-disorder abnormalities are consistently subtended by a small subset of network epicentres in bilateral sensory-motor, inferior temporal lobe, precuneus, and superior parietal cortex. Collectively, our results highlight how local molecular attributes and global connectivity jointly shape cross-disorder cortical abnormalities.

Bidirectionally connected cores in a mouse connectome: towards extracting the brain subnetworks essential for consciousness

Where in the brain consciousness resides remains unclear. It has been suggested that the subnetworks supporting consciousness should be bidirectionally (recurrently) connected because both feed-forward and feedback processing are necessary for conscious experience. Accordingly, evaluating which subnetworks are bidirectionally connected and the strength of these connections would likely aid the identification of regions essential to consciousness. Here, we propose a method for hierarchically decomposing a network into cores with different strengths of bidirectional connection, as a means of revealing the structure of the complex brain network. We applied the method to a whole-brain mouse connectome. We found that cores with strong bidirectional connections consisted of regions presumably essential to consciousness (e.g. the isocortical and thalamic regions, and claustrum) and did not include regions presumably irrelevant to consciousness (e.g. cerebellum). Contrarily, we could not find such correspondence between cores and consciousness when we applied other simple methods that ignored bidirectionality. These findings suggest that our method provides a novel insight into the relation between bidirectional brain network structures and consciousness.

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.

Null models in network neuroscience

Recent advances in imaging and tracing technology provide increasingly detailed reconstructions of brain connectomes. Concomitant analytic advances enable rigorous identification and quantification of functionally important features of brain network architecture. Null models are a flexible tool to statistically benchmark the presence or magnitude of features of interest, by selectively preserving specific architectural properties of brain networks while systematically randomizing others. Here we describe the logic, implementation and interpretation of null models of connectomes. We introduce randomization and generative approaches to constructing null networks, and outline a taxonomy of network methods for statistical inference. We highlight the spectrum of null models - from liberal models that control few network properties, to conservative models that recapitulate multiple properties of empirical networks - that allow us to operationalize and test detailed hypotheses about the structure and function of brain networks. We review emerging scenarios for the application of null models in network neuroscience, including for spatially embedded networks, annotated networks and correlation-derived networks. Finally, we consider the limits of null models, as well as outstanding questions for the field.

Extracting Dynamical Understanding From Neural-Mass Models of Mouse Cortex

New brain atlases with high spatial resolution and whole-brain coverage have rapidly advanced our knowledge of the brain's neural architecture, including the systematic variation of excitatory and inhibitory cell densities across the mammalian cortex. But understanding how the brain's microscale physiology shapes brain dynamics at the macroscale has remained a challenge. While physiologically based mathematical models of brain dynamics are well placed to bridge this explanatory gap, their complexity can form a barrier to providing clear mechanistic interpretation of the dynamics they generate. In this work, we develop a neural-mass model of the mouse cortex and show how bifurcation diagrams, which capture local dynamical responses to inputs and their variation across brain regions, can be used to understand the resulting whole-brain dynamics. We show that strong fits to resting-state functional magnetic resonance imaging (fMRI) data can be found in surprisingly simple dynamical regimes-including where all brain regions are confined to a stable fixed point-in which regions are able to respond strongly to variations in their inputs, consistent with direct structural connections providing a strong constraint on functional connectivity in the anesthetized mouse. We also use bifurcation diagrams to show how perturbations to local excitatory and inhibitory coupling strengths across the cortex, constrained by cell-density data, provide spatially dependent constraints on resulting cortical activity, and support a greater diversity of coincident dynamical regimes. Our work illustrates methods for visualizing and interpreting model performance in terms of underlying dynamical mechanisms, an approach that is crucial for building explanatory and physiologically grounded models of the dynamical principles that underpin large-scale brain activity.

Caught in vicious circles: a perspective on dynamic feed-forward loops driving oxidative stress in schizophrenia

A growing body of evidence has emerged demonstrating a pathological link between oxidative stress and schizophrenia. This evidence identifies oxidative stress as a convergence point or "central hub" for schizophrenia genetic and environmental risk factors. Here we review the existing experimental and translational research pinpointing the complex dynamics of oxidative stress mechanisms and their modulation in relation to schizophrenia pathophysiology. We focus on evidence supporting the crucial role of either redox dysregulation, N-methyl-D-aspartate receptor hypofunction, neuroinflammation or mitochondria bioenergetics dysfunction, initiating "vicious circles" centered on oxidative stress during neurodevelopment. These processes would amplify one another in positive feed-forward loops, leading to persistent impairments of the maturation and function of local parvalbumin-GABAergic neurons microcircuits and myelinated fibers of long-range macrocircuitry. This is at the basis of neural circuit synchronization impairments and cognitive, emotional, social and sensory deficits characteristic of schizophrenia. Potential therapeutic approaches that aim at breaking these different vicious circles represent promising strategies for timely and safe interventions. In order to improve early detection and increase the signal-to-noise ratio for adjunctive trials of antioxidant, anti-inflammatory and NMDAR modulator drugs, a reverse translation of validated circuitry approach is needed. The above presented processes allow to identify mechanism based biomarkers guiding stratification of homogenous patients groups and target engagement required for successful clinical trials, paving the way towards precision medicine in psychiatry.

Network structure and transcriptomic vulnerability shape atrophy in frontotemporal dementia

Connections among brain regions allow pathological perturbations to spread from a single source region to multiple regions. Patterns of neurodegeneration in multiple diseases, including behavioural variant of frontotemporal dementia (bvFTD), resemble the large-scale functional systems, but how bvFTD-related atrophy patterns relate to structural network organization remains unknown. Here we investigate whether neurodegeneration patterns in sporadic and genetic bvFTD are conditioned by connectome architecture. Regional atrophy patterns were estimated in both genetic bvFTD (75 patients, 247 controls) and sporadic bvFTD (70 patients, 123 controls). First, we identified distributed atrophy patterns in bvFTD, mainly targeting areas associated with the limbic intrinsic network and insular cytoarchitectonic class. Regional atrophy was significantly correlated with atrophy of structurally- and functionally-connected neighbours, demonstrating that network structure shapes atrophy patterns. The anterior insula was identified as the predominant group epicentre of brain atrophy using data-driven and simulation-based methods, with some secondary regions in frontal ventromedial and antero-medial temporal areas. We found that FTD-related genes, namely C9orf72 and TARDBP, confer local transcriptomic vulnerability to the disease, modulating the propagation of pathology through the connectome. Collectively, our results demonstrate that atrophy patterns in sporadic and genetic bvFTD are jointly shaped by global connectome architecture and local transcriptomic vulnerability, providing an explanation as to how heterogenous pathological entities can lead to the same clinical syndrome.

Expression of CRY2 Gene in the Brain Is Related to Human Navigation

Navigation is a complex cognitive process. CRY2 gene has been proposed to play an important role in navigation behaviors in various non-human animal species. Utilizing a recently developed neuroimaging-transcriptomics approach, the present study reported a tentative link between the CRY2 gene and human navigation. Specifically, we showed a significant pattern similarity between CRY2 gene expression in the human brain and navigation-related neural activation in functional magnetic resonance imaging. To further illuminate the functionality of CRY2 in human navigation, we examined the correlation between CRY2 expression and various cognitive processes underlying navigation, and found high correlation of CRY2 expression with neural activity of multiple cognitive domains, particularly object and shape perception and spatial memory. Further analyses on the relation between the neural activity of human navigation and the expression maps of genes of two CRY2-related pathways, i.e., the magnetoreceptive and circadian-related functions, found a trend of correlation for the CLOCK gene, a core circadian regulator gene, suggesting that CRY2 may modulate human navigation through its role in circadian rhythm. This observation was further confirmed by a behavioral study where individuals with better circadian regularity in daily life showed better sense of direction. Taken together, our study presents the first neural evidence that links CRY2 with human navigation, possibly through the modulation of circadian rhythm.

Where the genome meets the connectome: Understanding how genes shape human brain connectivity

The integration of modern neuroimaging methods with genetically informative designs and data can shed light on the molecular mechanisms underlying the structural and functional organization of the human connectome. Here, we review studies that have investigated the genetic basis of human brain network structure and function through three complementary frameworks: (1) the quantification of phenotypic heritability through classical twin designs; (2) the identification of specific DNA variants linked to phenotypic variation through association and related studies; and (3) the analysis of correlations between spatial variations in imaging phenotypes and gene expression profiles through the integration of neuroimaging and transcriptional atlas data. We consider the basic foundations, strengths, limitations, and discoveries associated with each approach. We present converging evidence to indicate that anatomical connectivity is under stronger genetic influence than functional connectivity and that genetic influences are not uniformly distributed throughout the brain, with phenotypic variation in certain regions and connections being under stronger genetic control than others. We also consider how the combination of imaging and genetics can be used to understand the ways in which genes may drive brain dysfunction in different clinical disorders.

Characteristics of Neural Network Changes in Normal Aging and Early Dementia

To understand the mechanisms underlying preserved and impaired cognitive function in healthy aging and dementia, respectively, the spatial relationships of brain networks and mechanisms of their resilience should be understood. The hub regions of the brain, such as the multisensory integration and default mode networks, are critical for within- and between-network communication, remain well-preserved during aging, and play an essential role in compensatory processes. On the other hand, these brain hubs are the preferred sites for lesions in neurodegenerative dementias, such as Alzheimer's disease. Disrupted primary information processing networks, such as the auditory, visual, and sensorimotor networks, may lead to overactivity of the multisensory integration networks and accumulation of pathological proteins that cause dementia. At the cellular level, the brain hub regions contain many synapses and require a large amount of energy. These regions are rich in ATP-related gene expression and had high glucose metabolism as demonstrated on positron emission tomography (PET). Importantly, the number and function of mitochondria, which are the center of ATP production, decline by about 8% every 10 years. Dementia patients often have dysfunction of the ubiquitin-proteasome and autophagy-lysosome systems, which require large amounts of ATP. If there is low energy supply but the demand is high, the risk of disease can be high. Imbalance between energy supply and demand may cause accumulation of pathological proteins and play an important role in the development of dementia. This energy imbalance may explain why brain hub regions are vulnerable to damage in different dementias. Here, we review (1) the characteristics of gray matter network, white matter network, and resting state functional network changes related to resilience in healthy aging, (2) the mode of resting state functional network disruption in neurodegenerative dementia, and (3) the cellular mechanisms associated with the disruption.

NuMorph: Tools for cortical cellular phenotyping in tissue-cleared whole-brain images

Tissue-clearing methods allow every cell in the mouse brain to be imaged without physical sectioning. However, the computational tools currently available for cell quantification in cleared tissue images have been limited to counting sparse cell populations in stereotypical mice. Here, we introduce NuMorph, a group of analysis tools to quantify all nuclei and nuclear markers within the mouse cortex after clearing and imaging by light-sheet microscopy. We apply NuMorph to investigate two distinct mouse models: a Topoisomerase 1 (Top1) model with severe neurodegenerative deficits and a Neurofibromin 1 (Nf1) model with a more subtle brain overgrowth phenotype. In each case, we identify differential effects of gene deletion on individual cell-type counts and distribution across cortical regions that manifest as alterations of gross brain morphology. These results underline the value of whole-brain imaging approaches, and the tools are widely applicable for studying brain structure phenotypes at cellular resolution.

Imaging Transcriptomics of Brain Disorders

Noninvasive neuroimaging is a powerful tool for quantifying diverse aspects of brain structure and function in vivo, and it has been used extensively to map the neural changes associated with various brain disorders. However, most neuroimaging techniques offer only indirect measures of underlying pathological mechanisms. The recent development of anatomically comprehensive gene expression atlases has opened new opportunities for studying the transcriptional correlates of noninvasively measured neural phenotypes, offering a rich framework for evaluating pathophysiological hypotheses and putative mechanisms. Here, we provide an overview of some fundamental methods in imaging transcriptomics and outline their application to understanding brain disorders of neurodevelopment, adulthood, and neurodegeneration. Converging evidence indicates that spatial variations in gene expression are linked to normative changes in brain structure during age-related maturation and neurodegeneration that are in part associated with cell-specific gene expression markers of gene expression. Transcriptional correlates of disorder-related neuroimaging phenotypes are also linked to transcriptionally dysregulated genes identified in ex vivo analyses of patient brains. Modeling studies demonstrate that spatial patterns of gene expression are involved in regional vulnerability to neurodegeneration and the spread of disease across the brain. This growing body of work supports the utility of transcriptional atlases in testing hypotheses about the molecular mechanism driving disease-related changes in macroscopic neuroimaging phenotypes.

Mapping gene transcription and neurocognition across human neocortex

Regulation of gene expression drives protein interactions that govern synaptic wiring and neuronal activity. The resulting coordinated activity among neuronal populations supports complex psychological processes, yet how gene expression shapes cognition and emotion remains unknown. Here, we directly bridge the microscale and macroscale by mapping gene expression patterns to functional activation patterns across the cortical sheet. Applying unsupervised learning to the Allen Human Brain Atlas and Neurosynth databases, we identify a ventromedial-dorsolateral gradient of gene assemblies that separate affective and perceptual domains. This topographic molecular-psychological signature reflects the hierarchical organization of the neocortex, including systematic variations in cell type, myeloarchitecture, laminar differentiation and intrinsic network affiliation. In addition, this molecular-psychological signature strengthens over neurodevelopment and can be replicated in two independent repositories. Collectively, our results reveal spatially covarying transcriptomic and cognitive architectures, highlighting the influence that molecular mechanisms exert on psychological processes.

The human connectome in Alzheimer disease - relationship to biomarkers and genetics

The pathology of Alzheimer disease (AD) damages structural and functional brain networks, resulting in cognitive impairment. The results of recent connectomics studies have now linked changes in structural and functional network organization in AD to the patterns of amyloid-β and tau accumulation and spread, providing insights into the neurobiological mechanisms of the disease. In addition, the detection of gene-related connectome changes might aid in the early diagnosis of AD and facilitate the development of personalized therapeutic strategies that are effective at earlier stages of the disease spectrum. In this article, we review studies of the associations between connectome changes and amyloid-β and tau pathologies as well as molecular genetics in different subtypes and stages of AD. We also highlight the utility of connectome-derived computational models for replicating empirical findings and for tracking and predicting the progression of biomarker-indicated AD pathophysiology.