BRAIN NETWORKS. Correlated gene expression supports synchronous activity in brain networks

Global spatiotemporal synchronizing structures of spontaneous neural activities in different cell types

Increasing evidence has revealed the large-scale nonstationary synchronizations as traveling waves in spontaneous neural activity. However, the interplay of various cell types in fine-tuning these spatiotemporal patters remains unclear. Here, we performed comprehensive exploration of spatiotemporal synchronizing structures across different cell types, states (awake, anesthesia, motion) and developmental axis in male mice. We found traveling waves in glutamatergic neurons exhibited greater variety than those in GABAergic neurons. Moreover, the synchronizing structures of GABAergic neurons converged toward those of glutamatergic neurons during development, but the evolution of waves exhibited varying timelines for different sub-type interneurons. Functional connectivity arises from both standing and traveling waves, and negative connections can be elucidated by the spatial propagation of waves. In addition, some traveling waves were correlated with the spatial distribution of gene expression. Our findings offer further insights into the neural underpinnings of traveling waves, functional connectivity, and resting-state networks, with cell-type specificity and developmental perspectives.

The resting-state brain activity signatures for addictive disorders

BACKGROUND

Addiction is a chronic and relapsing brain disorder. Despite numerous neuroimaging and neurophysiological studies on individuals with substance use disorder (SUD) or behavioral addiction (BEA), currently a clear neural activity signature for the addicted brain is lacking.

METHODS

We first performed systemic coordinate-based meta-analysis and partial least-squares regression to identify shared or distinct brain regions across multiple addictive disorders, with abnormal resting-state activity in SUD and BEA based on 46 studies (55 contrasts), including regional homogeneity (ReHo) and low-frequency fluctuation amplitude (ALFF) or fractional ALFF. We then combined Neurosynth, postmortem gene expression, and receptor/transporter distribution data to uncover the potential molecular mechanisms underlying these neural activity signatures.

FINDINGS

The overall comparison between addiction cohorts and healthy subjects indicated significantly increased ReHo and ALFF in the right striatum (putamen) and bilateral supplementary motor area, as well as decreased ReHo and ALFF in the bilateral anterior cingulate cortex and ventral medial prefrontal cortex, in the addiction group. On the other hand, neural activity in cingulate cortex, ventral medial prefrontal cortex, and orbitofrontal cortex differed between SUD and BEA subjects. Using molecular analyses, the altered resting activity recapitulated the spatial distribution of dopaminergic, GABAergic, and acetylcholine system in SUD, while this also includes the serotonergic system in BEA.

CONCLUSIONS

These results indicate both common and distinctive neural substrates underlying SUD and BEA, which validates and supports targeted neuromodulation against addiction.

FUNDING

This work was supported by the National Natural Science Foundation of China and Intramural Research Program of the National Institute on Drug Abuse, National Institutes of Health.

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.

Spatio-molecular profiles shape the human cerebellar hierarchy along the sensorimotor-association axis

Cerebellar involvement in both motor and non-motor functions manifests in specific regions of the human cerebellum, revealing the functional heterogeneity within it. One compelling theory places the heterogeneity within the cerebellar functional hierarchy along the sensorimotor-association (SA) axis. Despite extensive neuroimaging studies, evidence for the cerebellar SA axis from different modalities and scales was lacking. Thus, we establish a significant link between the cerebellar SA axis and spatio-molecular profiles. Utilizing the gene set variation analysis, we find the intermediate biological principles the significant genes leveraged to scaffold the cerebellar SA axis. Interestingly, we find these spatio-molecular profiles notably associated with neuropsychiatric dysfunction and recent evolution. Furthermore, cerebello-cerebral interactions at genetic and functional connectivity levels mirror the cerebral cortex and cerebellum's SA axis. These findings can provide a deeper understanding of how the human cerebellar SA axis is shaped and its role in transitioning from sensorimotor to association functions.

Spatial transcriptomic patterns underlying amyloid-β and tau pathology are associated with cognitive dysfunction in Alzheimer's disease

Amyloid-β (Aβ) and tau proteins accumulate within distinct neuronal systems in Alzheimer's disease (AD). Although it is not clear why certain brain regions are more vulnerable to Aβ and tau pathologies than others, gene expression may play a role. We study the association between brain-wide gene expression profiles and regional vulnerability to Aβ (gene-to-Aβ associations) and tau (gene-to-tau associations) pathologies by leveraging two large independent AD cohorts. We identify AD susceptibility genes and gene modules in a gene co-expression network with expression profiles specifically related to regional vulnerability to Aβ and tau pathologies in AD. In addition, we identify distinct biochemical pathways associated with the gene-to-Aβ and the gene-to-tau associations. These findings may explain the discordance between regional Aβ and tau pathologies. Finally, we propose an analytic framework, linking the identified gene-to-pathology associations to cognitive dysfunction in AD at the individual level, suggesting potential clinical implication of the gene-to-pathology associations.

Parkinson's disease and schizophrenia interactomes contain temporally distinct gene clusters underlying comorbid mechanisms and unique disease processes

Genome-wide association studies suggest significant overlaps in Parkinson's disease (PD) and schizophrenia (SZ) risks, but the underlying mechanisms remain elusive. The protein-protein interaction network ('interactome') plays a crucial role in PD and SZ and can incorporate their spatiotemporal specificities. Therefore, to study the linked biology of PD and SZ, we compiled PD- and SZ-associated genes from the DisGeNET database, and constructed their interactomes using BioGRID and HPRD. We examined the interactomes using clustering and enrichment analyses, in conjunction with the transcriptomic data of 26 brain regions spanning foetal stages to adulthood available in the BrainSpan Atlas. PD and SZ interactomes formed four gene clusters with distinct temporal identities (Disease Gene Networks or 'DGNs'1-4). DGN1 had unique SZ interactome genes highly expressed across developmental stages, corresponding to a neurodevelopmental SZ subtype. DGN2, containing unique SZ interactome genes expressed from early infancy to adulthood, correlated with an inflammation-driven SZ subtype and adult SZ risk. DGN3 contained unique PD interactome genes expressed in late infancy, early and late childhood, and adulthood, and involved in mitochondrial pathways. DGN4, containing prenatally-expressed genes common to both the interactomes, involved in stem cell pluripotency and overlapping with the interactome of 22q11 deletion syndrome (comorbid psychosis and Parkinsonism), potentially regulates neurodevelopmental mechanisms in PD-SZ comorbidity. Our findings suggest that disrupted neurodevelopment (regulated by DGN4) could expose risk windows in PD and SZ, later elevating disease risk through inflammation (DGN2). Alternatively, variant clustering in DGNs may produce disease subtypes, e.g., PD-SZ comorbidity with DGN4, and early/late-onset SZ with DGN1/DGN2.

Identification of gene biomarkers for brain diseases via multi-network topological semantics extraction and graph convolutional network

BACKGROUND

Brain diseases pose a significant threat to human health, and various network-based methods have been proposed for identifying gene biomarkers associated with these diseases. However, the brain is a complex system, and extracting topological semantics from different brain networks is necessary yet challenging to identify pathogenic genes for brain diseases.

RESULTS

In this study, we present a multi-network representation learning framework called M-GBBD for the identification of gene biomarker in brain diseases. Specifically, we collected multi-omics data to construct eleven networks from different perspectives. M-GBBD extracts the spatial distributions of features from these networks and iteratively optimizes them using Kullback-Leibler divergence to fuse the networks into a common semantic space that represents the gene network for the brain. Subsequently, a graph consisting of both gene and large-scale disease proximity networks learns representations through graph convolution techniques and predicts whether a gene is associated which brain diseases while providing associated scores. Experimental results demonstrate that M-GBBD outperforms several baseline methods. Furthermore, our analysis supported by bioinformatics revealed CAMP as a significantly associated gene with Alzheimer's disease identified by M-GBBD.

CONCLUSION

Collectively, M-GBBD provides valuable insights into identifying gene biomarkers for brain diseases and serves as a promising framework for brain networks representation learning.

Molecular signatures of cortical expansion in the human fetal brain

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

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

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

Structural connectome architecture shapes the maturation of cortical morphology from childhood to adolescence

Cortical thinning is an important hallmark of the maturation of brain morphology during childhood and adolescence. However, the connectome-based wiring mechanism that underlies cortical maturation remains unclear. Here, we show cortical thinning patterns primarily located in the lateral frontal and parietal heteromodal nodes during childhood and adolescence, which are structurally constrained by white matter network architecture and are particularly represented using a network-based diffusion model. Furthermore, connectome-based constraints are regionally heterogeneous, with the largest constraints residing in frontoparietal nodes, and are associated with gene expression signatures of microstructural neurodevelopmental events. These results are highly reproducible in another independent dataset. These findings advance our understanding of network-level mechanisms and the associated genetic basis that underlies the maturational process of cortical morphology during childhood and adolescence.

DSM: Deep sequential model for complete neuronal morphology representation and feature extraction

The full morphology of single neurons is indispensable for understanding cell types, the basic building blocks in brains. Projecting trajectories are critical to extracting biologically relevant information from neuron morphologies, as they provide valuable information for both connectivity and cell identity. We developed an artificial intelligence method, deep sequential model (DSM), to extract concise, cell-type-defining features from projections across brain regions. DSM achieves more than 90% accuracy in classifying 12 major neuron projection types without compromising performance when spatial noise is present. Such remarkable robustness enabled us to efficiently manage and analyze several major full-morphology data sources, showcasing how characteristic long projections can define cell identities. We also succeeded in applying our model to both discovering previously unknown neuron subtypes and analyzing exceptional co-expressed genes involved in neuron projection circuits.

Engineered serum markers for non-invasive monitoring of gene expression in the brain

Measurement of gene expression in the brain requires invasive analysis of brain tissue or non-invasive methods that are limited by low sensitivity. Here we introduce a method for non-invasive, multiplexed, site-specific monitoring of endogenous gene or transgene expression in the brain through engineered reporters called released markers of activity (RMAs). RMAs consist of an easily detectable reporter and a receptor-binding domain that enables transcytosis across the brain endothelium. RMAs are expressed in the brain but exit into the blood, where they can be easily measured. We show that expressing RMAs at a single mouse brain site representing approximately 1% of the brain volume provides up to a 100,000-fold signal increase over the baseline. Expression of RMAs in tens to hundreds of neurons is sufficient for their reliable detection. We demonstrate that chemogenetic activation of cells expressing Fos-responsive RMA increases serum RMA levels >6-fold compared to non-activated controls. RMAs provide a non-invasive method for repeatable, multiplexed monitoring of gene expression in the intact animal brain.

Gray matter volume abnormalities in obsessive-compulsive disorder correlate with molecular and transcriptional profiles

Neuroimaging studies have consistently established altered brain structure in obsessive-compulsive disorder (OCD). However, the molecular and genetic mechanisms underlying structural brain abnormalities remain unclear. In this study, we aimed to investigate altered gray matter volume and its underlying molecular and genetic mechanisms in patients with OCD. Gray matter morphological abnormalities measured with voxel based morphometry analysis were identified in patients with OCD in comparison to sex- and age-matched healthy controls (HCs). Spatial correlations between gray matter morphological abnormalities and neurotransmitter maps were calculated to identify neurotransmitters relating to structural abnormalities. Structural abnormalities related genes were identified by conducting transcriptome-neuroimaging spatial correlations. Compared with HCs, patients with OCD demonstrated significant morphological abnormalities in distributed brain areas, including gray matter atrophy in the anterior cingulate and increased gray matter volume in the thalamus, caudate and precentral and postcentral gyrus. The morphological abnormalities were significantly associated with dopamine synthesis capacity and expression profiles of 1110 genes enriched for trans-synaptic signaling, regulation of membrane potential, modulation of chemical synaptic transmission, brain development, synapse organization and regulation of neurotransmitter levels. These results elucidate the molecular and transcriptional basis of altered gray matter morphology and build linking between molecular, transcriptional and neuroimaging information facilitating an integrative understanding of OCD.

Organization of thalamocortical structural covariance and a corresponding 3D atlas of the mouse thalamus

For information from sensory organs to be processed by the brain, it is usually passed to appropriate areas of the cerebral cortex. Almost all of this information passes through the thalamus, a relay structure that reciprocally connects to the vast majority of the cortex. The thalamus facilitates this information transfer through a set of thalamocortical connections that vary in cellular structure, molecular profiles, innervation patterns, and firing rates. Additionally, corticothalamic connections allow for intracortical information transfer through the thalamus. These efferent and afferent connections between the thalamus and cortex have been the focus of many studies, and the importance of cortical connectivity in defining thalamus anatomy is demonstrated by multiple studies that parcellate the thalamus based on cortical connectivity profiles. Here, we examine correlated morphological variation between the thalamus and cortex, or thalamocortical structural covariance. For each voxel in the thalamus as a seed, we construct a cortical structural covariance map that represents correlated cortical volume variation, and examine whether high structural covariance is observed in cortical areas that are functionally relevant to the seed. Then, using these cortical structural covariance maps as features, we subdivide the thalamus into six non-overlapping regions (clusters of voxels), and assess whether cortical structural covariance is associated with cortical connectivity that specifically originates from these regions. We show that cortical structural covariance is high in areas of the cortex that are functionally related to the seed voxel, cortical structural covariance varies along cortical depth, and sharp transitions in cortical structural covariance profiles are observed when varying seed locations in the thalamus. Subdividing the thalamus based on structural covariance, we additionally demonstrate that the six thalamic clusters of voxels stratify cortical structural covariance along the dorsal-ventral, medial-lateral, and anterior-posterior axes. These cluster-associated structural covariance patterns are prominently detected in cortical regions innervated by fibers projecting out of their related thalamic subdivisions. Together, these results advance our understanding of how the thalamus and the cortex couple in their volumes. Our results indicate that these volume correlations reflect functional organization and structural connectivity, and further provides a novel segmentation of the mouse thalamus that can be used to examine thalamic structural variation and thalamocortical structural covariation in disease models.

Functional connectome hierarchy dysfunction in Alzheimer's disease and its relationship with cognition and gene expression profiling

Numerous researches have shown that the human brain organizes as a continuum axis crossing from sensory motor to transmodal cortex. Functional network alterations were commonly found in Alzheimer's disease (AD). Whether the hierarchy of AD brain networks has changed and how these changes related to gene expression profiling and cognition is unclear. Using resting-state functional magnetic resonance imaging data from 233 subjects (185 AD patients and 48 healthy controls), we studied the changes in the functional network gradients in AD. Moreover, we investigated the relationships between gradient alterations and cognition, and gene expression profiling, respectively. We found that the second gradient organizes as a continuum axis crossing from the sensory motor to the transmodal cortex. Compared to the healthy controls, the secondary gradient scores of the visual and somatomotor network (SOM) increased significantly in AD, and the secondary gradient scores of default mode and frontoparietal network decreased significantly in AD. The secondary gradient scores of SOM and salience network (SAL) significantly positively correlated with memory function in AD. The secondary gradient in SAL also significantly positively correlated with language function. The AD-related second gradient alterations were spatially associated with the gene expression and the relevant genes enriched in neurobiology-related pathways, specially expressed in various tissues, cell types, and developmental stages. These findings suggested the changes in the functional network gradients in AD and deepened our understanding of the correlation between macroscopic gradient structure and microscopic gene expression profiling in AD.

Protosequences in human cortical organoids model intrinsic states in the developing cortex

Neuronal firing sequences are thought to be the basic building blocks of neural coding and information broadcasting within the brain. However, when sequences emerge during neurodevelopment remains unknown. We demonstrate that structured firing sequences are present in spontaneous activity of human brain organoids and ex vivo neonatal brain slices from the murine somatosensory cortex. We observed a balance between temporally rigid and flexible firing patterns that are emergent phenomena in human brain organoids and early postnatal murine somatosensory cortex, but not in primary dissociated cortical cultures. Our findings suggest that temporal sequences do not arise in an experience-dependent manner, but are rather constrained by an innate preconfigured architecture established during neurogenesis. These findings highlight the potential for brain organoids to further explore how exogenous inputs can be used to refine neuronal circuits and enable new studies into the genetic mechanisms that govern assembly of functional circuitry during early human brain development.

Functional Connectome Hierarchy in Schizotypy and Its Associations With Expression of Schizophrenia-Related Genes

BACKGROUND AND HYPOTHESIS

Schizotypy has been conceptualized as a continuum of symptoms with marked genetic, neurobiological, and sensory-cognitive overlaps to schizophrenia. Hierarchical organization represents a general organizing principle for both the cortical connectome supporting sensation-to-cognition continuum and gene expression variability across the cortex. However, a mapping of connectome hierarchy to schizotypy remains to be established. Importantly, the underlying changes of the cortical connectome hierarchy that mechanistically link gene expressions to schizotypy are unclear.

STUDY DESIGN

The present study applied novel connectome gradient on resting-state fMRI data from 1013 healthy young adults to investigate schizotypy-associated sensorimotor-to-transmodal connectome hierarchy and assessed its similarity with the connectome hierarchy of schizophrenia. Furthermore, normative and differential postmortem gene expression data were utilized to examine transcriptional profiles linked to schizotypy-associated connectome hierarchy.

STUDY RESULTS

We found that schizotypy was associated with a compressed functional connectome hierarchy. Moreover, the pattern of schizotypy-related hierarchy exhibited a positive correlation with the connectome hierarchy observed in schizophrenia. This pattern was closely colocated with the expression of schizophrenia-related genes, with the correlated genes being enriched in transsynaptic, receptor signaling and calcium ion binding.

CONCLUSIONS

The compressed connectome hierarchy suggests diminished functional system differentiation, providing a novel and holistic system-level basis for various sensory-cognition deficits in schizotypy. Importantly, its linkage with schizophrenia-altered hierarchy and schizophrenia-related gene expression yields new insights into the neurobiological continuum of psychosis. It also provides mechanistic insight into how gene variation may drive alterations in functional hierarchy, mediating biological vulnerability of schizotypy to schizophrenia.

Transcriptional signatures of the whole-brain voxel-wise resting-state functional network centrality alterations in schizophrenia

Neuroimaging studies have revealed that patients with schizophrenia exhibit disrupted resting-state functional connectivity. However, the inconsistent findings across these studies have hindered our comprehensive understanding of the functional connectivity changes associated with schizophrenia, and the molecular mechanisms associated with these alterations remain largely unclear. A quantitative meta-analysis was first conducted on 21 datasets, involving 1057 patients and 1186 healthy controls, to examine disrupted resting-state functional connectivity in schizophrenia, as measured by whole-brain voxel-wise functional network centrality (FNC). Subsequently, partial least squares regression analysis was employed to investigate the relationship between FNC changes and gene expression profiles obtained from the Allen Human Brain Atlas database. Finally, gene enrichment analysis was performed to unveil the biological significance of the altered FNC-related genes. Compared with healthy controls, patients with schizophrenia show consistently increased FNC in the right inferior parietal cortex extending to the supramarginal gyrus, angular gyrus, bilateral medial prefrontal cortex, and right dorsolateral prefrontal cortex, while decreased FNC in the bilateral insula, bilateral postcentral gyrus, and right inferior temporal gyrus. Meta-regression analysis revealed that increased FNC in the right inferior parietal cortex was positively correlated with clinical score. In addition, these observed functional connectivity changes were found to be spatially associated with the brain-wide expression of specific genes, which were enriched in diverse biological pathways and cell types. These findings highlight the aberrant functional connectivity observed in schizophrenia and its potential molecular underpinnings, providing valuable insights into the neuropathology of dysconnectivity associated with this disorder.

Transcriptional correlates of frequency-dependent brain functional activity associated with symptom severity in degenerative cervical myelopathy

BACKGROUND

Neuroimaging techniques provide insights into the brain abnormalities secondary to degenerative cervical myelopathy (DCM) and their association with neurological deficits. However, the neural correlates underlying the discrepancy between symptom severity and the degree of spinal cord compression, as well as the transcriptional correlates of these cortical abnormalities, remain unknown in DCM patients.

METHODS

In this cross-sectional study, which collected resting-state functional MRI (rs-fMRI) images and the Japanese Orthopedic Association (JOA) score, enrolled 104 participants (54 patients and 50 healthy controls). The frequency-dependent amplitude of low-frequency fluctuation (ALFF) was obtained for all participants. We investigated the ALFF differences between mild-symptom DCM patients and severe-symptom DCM patients while carefully matching the degree of compression between these two groups via both univariate comparison and searchlight classification for three frequency bands (e.g., Slow-4, Slow-5, and Full-band). Additionally, we identified genes associated with symptom severity in DCM patients by linking the spatial patterns of gene expression of Allen Human Brain Atlas and brain functional differences between mild symptom and severe symptom groups.

RESULTS

(1) We found that the frequency-specific brain activities within the sensorimotor network (SMN), visual network (VN), and default mode network (DMN) were associated with the varying degrees of functional impairment in DCM patients; (2) the frequency-specific brain activity within the SMN correlated with the functional recovery in patients with DCM; (3) a spatial correlation between the brain-wide expression of genes involved in neuronal migration and the brain functional activities associated with symptom severity was identified in DCM patients.

CONCLUSION

In conclusion, our study bridges gaps between genes, cell classes, biological processes, and brain functional correlates of DCM. While our findings are correlational in nature, they suggest that the neural activities of sensorimotor cortices in DCM are associated with the severity of symptoms and might be associated with neuronal migration within the brain.

Spatiotemporal expression patterns of anxiety disorder-associated genes

Anxiety disorders (ADs) are the most common form of mental disorder that affects millions of individuals worldwide. Although physiological studies have revealed the neural circuits related to AD symptoms, how AD-associated genes are spatiotemporally expressed in the human brain still remains unclear. In this study, we integrated genome-wide association studies of four human AD subtypes-generalized anxiety disorder, social anxiety disorder, panic disorder, and obsessive-compulsive disorder-with spatial gene expression patterns. Our investigation uncovered a novel division among AD-associated genes, marked by significant and distinct expression enrichments in the cerebral nuclei, limbic, and midbrain regions. Each gene cluster was associated with specific anxiety-related behaviors, signaling pathways, region-specific gene networks, and cell types. Notably, we observed a significant negative correlation in the temporal expression patterns of these gene clusters during various developmental stages. Moreover, the specific brain regions enriched in each gene group aligned with neural circuits previously associated with negative decision-making and anxious temperament. These results suggest that the two distinct gene clusters may underlie separate neural systems involved in anxiety. As a result, our findings bridge the gap between genes and neural circuitry, shedding light on the mechanisms underlying AD-associated behaviors.

Longitudinal development of the human white matter structural connectome and its association with brain transcriptomic and cellular architecture

From childhood to adolescence, the spatiotemporal development pattern of the human brain white matter connectome and its underlying transcriptomic and cellular mechanisms remain largely unknown. With a longitudinal diffusion MRI cohort of 604 participants, we map the developmental trajectory of the white matter connectome from global to regional levels and identify that most brain network properties followed a linear developmental trajectory. Importantly, connectome-transcriptomic analysis reveals that the spatial development pattern of white matter connectome is potentially regulated by the transcriptomic architecture, with positively correlated genes involve in ion transport- and development-related pathways expressed in excitatory and inhibitory neurons, and negatively correlated genes enriches in synapse- and development-related pathways expressed in astrocytes, inhibitory neurons and microglia. Additionally, the macroscale developmental pattern is also associated with myelin content and thicknesses of specific laminas. These findings offer insights into the underlying genetics and neural mechanisms of macroscale white matter connectome development from childhood to adolescence.

Culture shapes spontaneous brain dynamics - Shared versus idiosyncratic neural features among Chinese versus Canadian subjects

Environmental factors, such as culture, are known to shape individual variation in brain activity including spontaneous activity, but less is known about their population-level effects. Eastern and Western cultures differ strongly in their cultural norms about relationships between individuals. For example, the collectivism, interdependence and tightness of Eastern cultures relative to the individualism, independence and looseness of Western cultures, promote interpersonal connectedness and coordination. Do such cultural contexts therefore influence the group-level variability of their cultural members' spontaneous brain activity? Using novel methods adapted from studies of inter-subject neural synchrony, we compare the group-level variability of resting state EEG dynamics in Chinese and Canadian samples. We observe that Chinese subjects show significantly higher inter-subject correlation and lower inter-subject distance in their EEG power spectra than Canadian subjects, as well as lower variability in theta power and alpha peak frequency. We demonstrate, for the first time, different relationships among subjects' resting state brain dynamics in Chinese and Canadian samples. These results point to more idiosyncratic neural dynamics among Canadian participants, compared with more shared neural features in Chinese participants.

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.

Single-subject cortical morphological brain networks: Phenotypic associations and neurobiological substrates

Although single-subject morphological brain networks provide an important way for human connectome studies, their roles and origins are poorly understood. Combining cross-sectional and repeated structural magnetic resonance imaging scans from adults, children and twins with behavioral and cognitive measures and brain-wide transcriptomic, cytoarchitectonic and chemoarchitectonic data, this study examined phenotypic associations and neurobiological substrates of single-subject morphological brain networks. We found that single-subject morphological brain networks explained inter-individual variance and predicted individual outcomes in Motor and Cognition domains, and distinguished individuals from each other. The performance can be further improved by integrating different morphological indices for network construction. Low-moderate heritability was observed for single-subject morphological brain networks with the highest heritability for sulcal depth-derived networks and higher heritability for inter-module connections. Furthermore, differential roles of genetic, cytoarchitectonic and chemoarchitectonic factors were observed for single-subject morphological brain networks. Cortical thickness-derived networks were related to the three factors with contributions from genes enriched in membrane and transport related functions, genes preferentially located in supragranular and granular layers, overall thickness in the molecular layer and thickness of wall in the infragranular layers, and metabotropic glutamate receptor 5 and dopamine transporter; fractal dimension-, gyrification index- and sulcal depth-derived networks were only associated with the chemoarchitectonic factor with contributions from different sets of neurotransmitter receptors. Most results were reproducible across different parcellation schemes and datasets. Altogether, this study demonstrates phenotypic associations and neurobiological substrates of single-subject morphological brain networks, which provide intermediate endophenotypes to link molecular and cellular architecture and behavior and cognition.

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.

Functional genomics and systems biology in human neuroscience

Neuroscience research has entered a phase of key discoveries in the realm of neurogenomics owing to strong financial and intellectual support for resource building and tool development. The previous challenge of tissue heterogeneity has been met with the application of techniques that can profile individual cells at scale. Moreover, the ability to perturb genes, gene regulatory elements and neuronal activity in a cell-type-specific manner has been integrated with gene expression studies to uncover the functional underpinnings of the genome at a systems level. Although these insights have necessarily been grounded in model systems, we now have the opportunity to apply these approaches in humans and in human tissue, thanks to advances in human genetics, brain imaging and tissue collection. We acknowledge that there will probably always be limits to the extent to which we can apply the genomic tools developed in model systems to human neuroscience; however, as we describe in this Perspective, the neuroscience field is now primed with an optimal foundation for tackling this ambitious challenge. The application of systems-level network analyses to these datasets will facilitate a deeper appreciation of human neurogenomics that cannot otherwise be achieved from directly observable phenomena.

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.

Deciphering the neural mechanisms of miR-134 in major depressive disorder with population-based and person-specific imaging transcriptomic techniques

MiR-134 has emerged as a potential molecular biomarker for the detection and management of major depressive disorder (MDD). Nevertheless, the specific effects of miR-134 as a regulatory element on brain function and its implications for the clinical presentation of MDD are not yet fully understood. In order to investigate the potential neural mechanisms that contribute to the relationship between miR-134 and MDD, we employed a parallel two-stage cross-scale multi-omics approach. This involved utilizing the anterior cingulate cortex (ACC) functional connectivity as a means to connect microscopic molecular structures with macroscopic brain function in two separate cohorts: the MDD-I dataset (56 MDD patients and 51 healthy controls) and the MDD-II dataset (57 MDD patients and 52 healthy controls). We found a stable ACC functional dysconnectivity pattern of MDD and established the hierarchical cross-scale association from molecular organizations of miR-134 target genes to macroscopic brain functional dysconnectivity and associated behavior, as revealed by population-based analysis. Additionally, our person-specific imaging transcriptomic study revealed that individual exosomal miR-134 expression levels impact on individual clinical symptoms of MDD by modulating ACC-related functional dysconnectivity. Together, our findings provide compelling evidence of the correlation between miR-134 and depression across multi scales within the gene-brain-behavior context.

Interindividual variation in human cortical cell type abundance and expression

Single-cell transcriptomic studies have identified a conserved set of neocortical cell types from small postmortem cohorts. We extended these efforts by assessing cell type variation across 75 adult individuals undergoing epilepsy and tumor surgeries. Nearly all nuclei map to one of 125 robust cell types identified in the middle temporal gyrus. However, we found interindividual variance in abundances and gene expression signatures, particularly in deep-layer glutamatergic neurons and microglia. A minority of donor variance is explainable by age, sex, ancestry, disease state, and cell state. Genomic variation was associated with expression of 150 to 250 genes for most cell types. This characterization of cellular variation provides a baseline for cell typing in health and 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.

Connectome-based modelling of neurodegenerative diseases: towards precision medicine and mechanistic insight

Neurodegenerative diseases are the most common cause of dementia. Although their underlying molecular pathologies have been identified, there is substantial heterogeneity in the patterns of progressive brain alterations across and within these diseases. Recent advances in neuroimaging methods have revealed that pathological proteins accumulate along specific macroscale brain networks, implicating the network architecture of the brain in the system-level pathophysiology of neurodegenerative diseases. However, the extent to which 'network-based neurodegeneration' applies across the wide range of neurodegenerative disorders remains unclear. Here, we discuss the state-of-the-art of neuroimaging-based connectomics for the mapping and prediction of neurodegenerative processes. We review findings supporting brain networks as passive conduits through which pathological proteins spread. As an alternative view, we also discuss complementary work suggesting that network alterations actively modulate the spreading of pathological proteins between connected brain regions. We conclude this Perspective by proposing an integrative framework in which connectome-based models can be advanced along three dimensions of innovation: incorporating parameters that modulate propagation behaviour on the basis of measurable biological features; building patient-tailored models that use individual-level information and allowing model parameters to interact dynamically over time. We discuss promises and pitfalls of these strategies for improving disease insights and moving towards precision medicine.

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.

Network models to enhance the translational impact of cross-species studies

Neuroscience studies are often carried out in animal models for the purpose of understanding specific aspects of the human condition. However, the translation of findings across species remains a substantial challenge. Network science approaches can enhance the translational impact of cross-species studies by providing a means of mapping small-scale cellular processes identified in animal model studies to larger-scale inter-regional circuits observed in humans. In this Review, we highlight the contributions of network science approaches to the development of cross-species translational research in neuroscience. We lay the foundation for our discussion by exploring the objectives of cross-species translational models. We then discuss how the development of new tools that enable the acquisition of whole-brain data in animal models with cellular resolution provides unprecedented opportunity for cross-species applications of network science approaches for understanding large-scale brain networks. We describe how these tools may support the translation of findings across species and imaging modalities and highlight future opportunities. Our overarching goal is to illustrate how the application of network science tools across human and animal model studies could deepen insight into the neurobiology that underlies phenomena observed with non-invasive neuroimaging methods and could simultaneously further our ability to translate findings across species.

Neuroimaging based biotypes for precision diagnosis and prognosis in cancer-related cognitive impairment

Cancer related cognitive impairment (CRCI) is commonly associated with cancer and its treatments, yet the present binary diagnostic approach fails to capture the full spectrum of this syndrome. Cognitive function is highly complex and exists on a continuum that is poorly characterized by dichotomous categories. Advanced statistical methodologies applied to symptom assessments have demonstrated that there are multiple subclasses of CRCI. However, studies suggest that relying on symptom assessments alone may fail to account for significant differences in the neural mechanisms that underlie a specific cognitive phenotype. Treatment plans that address the specific physiologic mechanisms involved in an individual patient's condition is the heart of precision medicine. In this narrative review, we discuss how biotyping, a precision medicine framework being utilized in other mental disorders, could be applied to CRCI. Specifically, we discuss how neuroimaging can be used to determine biotypes of CRCI, which allow for increased precision in prediction and diagnosis of CRCI via biologic mechanistic data. Biotypes may also provide more precise clinical endpoints for intervention trials. Biotyping could be made more feasible with proxy imaging technologies or liquid biomarkers. Large cross-sectional phenotyping studies are needed in addition to evaluation of longitudinal trajectories, and data sharing/pooling is highly feasible with currently available digital infrastructures.

Spatial transcriptomic patterns underlying regional vulnerability to amyloid-β and tau pathologies and their relationships to cognitive dysfunction in Alzheimer's disease

Amyloid-β (Aβ) and tau proteins accumulate within distinct neuronal systems in Alzheimer's disease (AD). Although it is not clear why certain brain regions are more vulnerable to Aβ and tau pathologies than others, gene expression may play a role. We studied the association between brain-wide gene expression profiles and regional vulnerability to Aβ (gene-to-Aβ associations) and tau (gene-to-tau associations) pathologies leveraging two large independent cohorts (n = 715) of participants along the AD continuum. We identified several AD susceptibility genes and gene modules in a gene co-expression network with expression profiles related to regional vulnerability to Aβ and tau pathologies in AD. In particular, we found that the positive APOE -to-tau association was only seen in the AD cohort, whereas patients with AD and frontotemporal dementia shared similar positive MAPT -to-tau association. Some AD candidate genes showed sex-dependent negative gene-to-Aβ and gene-to-tau associations. In addition, we identified distinct biochemical pathways associated with the gene-to-Aβ and the gene-to-tau associations. Finally, we proposed a novel analytic framework, linking the identified gene-to-pathology associations to cognitive dysfunction in AD at the individual level, suggesting potential clinical implication of the gene-to-pathology associations. Taken together, our study identified distinct gene expression profiles and biochemical pathways that may explain the discordance between regional Aβ and tau pathologies, and filled the gap between gene-to-pathology associations and cognitive dysfunction in individual AD patients that may ultimately help identify novel personalized pathogenetic biomarkers and therapeutic targets.

One Sentence Summary

We identified replicable cognition-related associations between regional gene expression profiles and selectively regional vulnerability to amyloid-β and tau pathologies in AD.

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.

Brain mitochondrial diversity and network organization predict anxiety-like behavior in male mice

The brain and behavior are under energetic constraints, limited by mitochondrial energy transformation capacity. However, the mitochondria-behavior relationship has not been systematically studied at a brain-wide scale. Here we examined the association between multiple features of mitochondrial respiratory chain capacity and stress-related behaviors in male mice with diverse behavioral phenotypes. Miniaturized assays of mitochondrial respiratory chain enzyme activities and mitochondrial DNA (mtDNA) content were deployed on 571 samples across 17 brain areas, defining specific patterns of mito-behavior associations. By applying multi-slice network analysis to our brain-wide mitochondrial dataset, we identified three large-scale networks of brain areas with shared mitochondrial signatures. A major network composed of cortico-striatal areas exhibited the strongest mitochondria-behavior correlations, accounting for up to 50% of animal-to-animal behavioral differences, suggesting that this mito-based network is functionally significant. The mito-based brain networks also overlapped with regional gene expression and structural connectivity, and exhibited distinct molecular mitochondrial phenotype signatures. This work provides convergent multimodal evidence anchored in enzyme activities, gene expression, and animal behavior that distinct, behaviorally-relevant mitochondrial phenotypes exist across the male mouse brain.

Robust estimation of cortical similarity networks from brain MRI

Structural similarity is a growing focus for magnetic resonance imaging (MRI) of connectomes. Here we propose Morphometric INverse Divergence (MIND), a new method to estimate within-subject similarity between cortical areas based on the divergence between their multivariate distributions of multiple MRI features. Compared to the prior approach of morphometric similarity networks (MSNs) on n > 11,000 scans spanning three human datasets and one macaque dataset, MIND networks were more reliable, more consistent with cortical cytoarchitectonics and symmetry and more correlated with tract-tracing measures of axonal connectivity. MIND networks derived from human T1-weighted MRI were more sensitive to age-related changes than MSNs or networks derived by tractography of diffusion-weighted MRI. Gene co-expression between cortical areas was more strongly coupled to MIND networks than to MSNs or tractography. MIND network phenotypes were also more heritable, especially edges between structurally differentiated areas. MIND network analysis provides a biologically validated lens for cortical connectomics using readily available MRI data.

A Molecular Basis of Human Brain Connectivity

Neuroimaging is commonly used to infer human brain connectivity, but those measurements are far-removed from the molecular underpinnings at synapses. To uncover the molecular basis of human brain connectivity, we analyzed a unique cohort of 98 individuals who provided neuroimaging and genetic data contemporaneous with dendritic spine morphometric, proteomic, and gene expression data from the superior frontal and inferior temporal gyri. Through cellular contextualization of the molecular data with dendritic spine morphology, we identified hundreds of proteins related to synapses, energy metabolism, and RNA processing that explain between-individual differences in functional connectivity and structural covariation. By integrating data at the genetic, molecular, subcellular, and tissue levels, we bridged the divergent fields of molecular biology and neuroimaging to identify a molecular basis of brain connectivity.

One-Sentence Summary

Dendritic spine morphometry and synaptic proteins unite the divergent fields of molecular biology and neuroimaging.

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.

Noncoding transcripts are linked to brain resting-state activity in non-human primates

Brain-derived transcriptomes are known to correlate with resting-state brain activity in humans. Whether this association holds in nonhuman primates remains uncertain. Here, we search for such molecular correlates by integrating 757 transcriptomes derived from 100 macaque cortical regions with resting-state activity in separate conspecifics. We observe that 150 noncoding genes explain variations in resting-state activity at a comparable level with protein-coding genes. In-depth analysis of these noncoding genes reveals that they are connected to the function of nonneuronal cells such as oligodendrocytes. Co-expression network analysis finds that the modules of noncoding genes are linked to both autism and schizophrenia risk genes. Moreover, genes associated with resting-state noncoding genes are highly enriched in human resting-state functional genes and memory-effect genes, and their links with resting-state functional magnetic resonance imaging (fMRI) signals are altered in the brains of patients with autism. Our results highlight the potential for noncoding RNAs to explain resting-state activity in the nonhuman primate brain.

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.

Transcriptomic and Macroscopic Architectures of Multimodal Covariance Network Reveal Molecular-Structural-Functional Co-alterations

Human cognition is usually underpinned by intrinsic structure and functional neural co-activation in spatially distributed brain regions. Owing to lacking an effective approach to quantifying the covarying of structure and functional responses, how the structural-functional circuits interact and how genes encode the relationships, to deepen our knowledge of human cognition and disease, are still unclear. Here, we propose a multimodal covariance network (MCN) construction approach to capture interregional covarying of the structural skeleton and transient functional activities for a single individual. We further explored the potential association between brain-wide gene expression patterns and structural-functional covarying in individuals involved in a gambling task and individuals with major depression disorder (MDD), adopting multimodal data from a publicly available human brain transcriptomic atlas and 2 independent cohorts. MCN analysis showed a replicable cortical structural-functional fine map in healthy individuals, and the expression of cognition- and disease phenotype-related genes was found to be spatially correlated with the corresponding MCN differences. Further analysis of cell type-specific signature genes suggests that the excitatory and inhibitory neuron transcriptomic changes could account for most of the observed correlation with task-evoked MCN differences. In contrast, changes in MCN of MDD patients were enriched for biological processes related to synapse function and neuroinflammation in astrocytes, microglia, and neurons, suggesting its promising application in developing targeted therapies for MDD patients. Collectively, these findings confirmed the correlations of MCN-related differences with brain-wide gene expression patterns, which captured genetically validated structural-functional differences at the cellular level in specific cognitive processes and psychiatric patients.

Structural covariances of prefrontal subregions selectively associate with dopamine-related gene coexpression and schizophrenia

Evidence highlights that dopamine (DA) system dysregulation and prefrontal cortex (PFC) dysfunction may underlie the pathophysiology of schizophrenia. However, the associations among DA genes, PFC morphometry, and schizophrenia have not yet been fully clarified. Based on the brain gene expression dataset from Allen Human Brain Atlas and structural magnetic resonance imaging data (NDIS = 1727, NREP = 408), we first identified 10 out of 22 PFC subregions whose gray matter volume (GMV) covariance profiles were reliably associated with their DA genes coexpression profiles, then four out of the identified 10 PFC subregions demonstrated abnormally increased GMV covariance with the hippocampus, insula, and medial frontal areas in schizophrenia patients (NCASE = 100; NCONTROL = 102). Moreover, based on a schizophrenia postmortem expression dataset, we found that the DA genes coexpression of schizophrenia was significantly reduced between the middle frontal gyrus and hippocampus, in which 21 DA genes showed significantly unsynchronized expression changes, and the 21 genes' brain expression were enriched in brain activity invoked by working memory, reward, speech production, and episodic memory. Our findings indicate the DA genes selectively regulate the structural covariance of PFC subregions by their coexpression profiles, which may underlie the disrupted GMV covariance and impaired cognitive functions in schizophrenia.

Acoustically-Targeted Measurement of Transgene Expression in the Brain

Gene expression is a critical component of brain physiology and activity, but monitoring this expression in the living brain represents a significant challenge. Here, we introduce a new paradigm called Recovery of Markers through InSonation (REMIS) for noninvasive measurement of gene expression in the brain with cell-type, spatial, and temporal specificity. Our approach relies on engineered protein markers that are designed to be expressed in neurons and exit into the interstitium. By applying ultrasound to targeted brain regions, these markers are released into the bloodstream, where they can be readily detected using biochemical techniques. REMIS can noninvasively confirm gene delivery and measure endogenous signaling in specific brain sites through a simple insonation and a subsequent blood test. Using REMIS, we successfully measured chemogenetic induction of neuronal activity in ultrasound-tar-geted brain regions. REMIS recovery of markers is reliable and demonstrated improved recovery of markers from the brain into the blood in every tested animal. Overall, our work establishes a noninvasive, spatially-specific means of monitoring gene delivery outcomes and endogenous signaling in mammalian brains, opening up possibilities for brain research and noninvasive monitoring of gene therapies in the brain.

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.

Cross-sectional and longitudinal neuroanatomical profiles of distinct clinical (adaptive) outcomes in autism

Individuals with autism spectrum disorder (henceforth referred to as autism) display significant variation in clinical outcome. For instance, across age, some individuals' adaptive skills naturally improve or remain stable, while others' decrease. To pave the way for 'precision-medicine' approaches, it is crucial to identify the cross-sectional and, given the developmental nature of autism, longitudinal neurobiological (including neuroanatomical and linked genetic) correlates of this variation. We conducted a longitudinal follow-up study of 333 individuals (161 autistic and 172 neurotypical individuals, aged 6-30 years), with two assessment time points separated by ~12-24 months. We collected behavioural (Vineland Adaptive Behaviour Scale-II, VABS-II) and neuroanatomical (structural magnetic resonance imaging) data. Autistic participants were grouped into clinically meaningful "Increasers", "No-changers", and "Decreasers" in adaptive behaviour (based on VABS-II scores). We compared each clinical subgroup's neuroanatomy (surface area and cortical thickness at T1, ∆T (intra-individual change) and T2) to that of the neurotypicals. Next, we explored the neuroanatomical differences' potential genomic associates using the Allen Human Brain Atlas. Clinical subgroups had distinct neuroanatomical profiles in surface area and cortical thickness at baseline, neuroanatomical development, and follow-up. These profiles were enriched for genes previously associated with autism and for genes previously linked to neurobiological pathways implicated in autism (e.g. excitation-inhibition systems). Our findings suggest that distinct clinical outcomes (i.e. intra-individual change in clinical profiles) linked to autism core symptoms are associated with atypical cross-sectional and longitudinal, i.e. developmental, neurobiological profiles. If validated, our findings may advance the development of interventions, e.g. targeting mechanisms linked to relatively poorer outcomes.

Cortical gradients during naturalistic processing are hierarchical and modality-specific

Understanding cortical topographic organization and how it supports complex perceptual and cognitive processes is a fundamental question in neuroscience. Previous work has characterized functional gradients that demonstrate large-scale principles of cortical organization. How these gradients are modulated by rich ecological stimuli remains unknown. Here, we utilize naturalistic stimuli via movie-fMRI to assess macroscale functional organization. We identify principal movie gradients that delineate separate hierarchies anchored in sensorimotor, visual, and auditory/language areas. At the opposite/heteromodal end of these perception-to-cognition axes, we find a more central role for the frontoparietal network along with the default network. Even across different movie stimuli, movie gradients demonstrated good reliability, suggesting that these hierarchies reflect a brain state common across different naturalistic conditions. The relative position of brain areas within movie gradients showed stronger and more numerous correlations with cognitive behavioral scores compared to resting state gradients. Together, these findings provide an ecologically valid representation of the principles underlying cortical organization while the brain is active and engaged in multimodal, dynamic perceptual and cognitive processing.

From Neurotransmitters to Networks: Transcending Organisational Hierarchies with Molecular-informed Functional Imaging

The human brain exhibits complex interactions across micro, meso-, and macro-scale organisational principles. Recent synergistic multi-modal approaches have begun to link micro-scale information to systems level dynamics, transcending organisational hierarchies and offering novel perspectives into the brain's function and dysfunction. Specifically, the distribution of micro-scale properties (such as receptor density or gene expression) can be mapped onto macro-scale measures from functional MRI to provide novel neurobiological insights. Methodological approaches to enrich functional imaging analyses with molecular information are rapidly evolving, with several streams of research having developed relatively independently, each offering unique potential to explore the trans-hierarchical functioning of the brain. Here, we address the three principal streams of research - spatial correlation, molecular-enriched network, and in-silico whole brain modelling analyses - to provide a critical overview of the different sources of molecular information, how this information can be utilised within analyses of fMRI data, the merits and pitfalls of each methodology, and, through the use of key examples, highlight their promise to shed new light on key domains of neuroscientific inquiry.