Overcoming false-positive gene-category enrichment in the analysis of spatially resolved transcriptomic brain atlas data

Transcriptomic decoding of regional cortical vulnerability to major depressive disorder

Previous studies in small samples have identified inconsistent cortical abnormalities in major depressive disorder (MDD). Despite genetic influences on MDD and the brain, it is unclear how genetic risk for MDD is translated into spatially patterned cortical vulnerability. Here, we initially examined voxel-wise differences in cortical function and structure using the largest multi-modal MRI data from 1660 MDD patients and 1341 controls. Combined with the Allen Human Brain Atlas, we then adopted transcription-neuroimaging spatial correlation and the newly developed ensemble-based gene category enrichment analysis to identify gene categories with expression related to cortical changes in MDD. Results showed that patients had relatively circumscribed impairments in local functional properties and broadly distributed disruptions in global functional connectivity, consistently characterized by hyper-function in associative areas and hypo-function in primary regions. Moreover, the local functional alterations were correlated with genes enriched for biological functions related to MDD in general (e.g., endoplasmic reticulum stress, mitogen-activated protein kinase, histone acetylation, and DNA methylation); and the global functional connectivity changes were associated with not only MDD-general, but also brain-relevant genes (e.g., neuron, synapse, axon, glial cell, and neurotransmitters). Our findings may provide important insights into the transcriptomic signatures of regional cortical vulnerability to MDD.

Quantitative susceptibility mapping of multiple system atrophy and Parkinson's disease correlates with neurotransmitter reference maps

BACKGROUND

Multiple system atrophy (MSA) and Parkinson's disease (PD) are neurodegenerative disorders characterized by α-synuclein pathology, disrupted iron homeostasis and impaired neurochemical transmission. Considering the critical role of iron in neurotransmitter synthesis and transport, our study aims to identify distinct patterns of whole-brain iron accumulation in MSA and PD, and to elucidate the corresponding neurochemical substrates.

METHODS

A total of 122 PD patients, 58 MSA patients and 78 age-, sex-matched health controls underwent multi-echo gradient echo sequences and neurological evaluations. We conducted voxel-wise and regional analyses using quantitative susceptibility mapping to explore MSA or PD-specific alterations in cortical and subcortical iron concentrations. Spatial correlation approaches were employed to examine the topographical alignment of cortical iron accumulation patterns with normative atlases of neurotransmitter receptor and transporter densities. Furthermore, we assessed the associations between the colocalization strength of neurochemical systems and disease severity.

RESULTS

MSA patients exhibited increased susceptibility in the striatal, midbrain, cerebellar nuclei, as well as the frontal, temporal, occipital lobes, and anterior cingulate gyrus. In contrast, PD patients displayed elevated iron levels in the left inferior occipital gyrus, precentral gyrus, and substantia nigra. The excessive iron accumulation in MSA or PD correlated with the spatial distribution of cholinergic, noradrenaline, glutamate, serotonin, cannabinoids, and opioid neurotransmitters, and the degree of this alignment was related to motor deficits.

CONCLUSIONS

Our findings provide evidence of the interaction between iron accumulation and non-dopamine neurotransmitters in the pathogenesis of MSA and PD, which inspires research on potential targets for pharmacotherapy.

A shifting role of thalamocortical connectivity in the emergence of cortical functional organization

The cortical patterning principle has been a long-standing question in neuroscience, yet how this translates to macroscale functional specialization in the human brain remains largely unknown. Here we examine age-dependent differences in resting-state thalamocortical connectivity to investigate its role in the emergence of large-scale functional networks during early life, using a primarily cross-sectional but also longitudinal approach. We show that thalamocortical connectivity during infancy reflects an early differentiation of sensorimotor networks and genetically influenced axonal projection. This pattern changes in childhood, when connectivity is established with the salience network, while decoupling externally and internally oriented functional systems. A developmental simulation using generative network models corroborated these findings, demonstrating that thalamic connectivity contributes to developing key features of the mature brain, such as functional segregation and the sensory-association axis, especially across 12-18 years of age. Our study suggests that the thalamus plays an important role in functional specialization during development, with potential implications for studying conditions with compromised internal and external processing.

Functional magnetic resonance imaging alternations in suicide attempts individuals and their association with gene expression

BACKGROUND

Functional Magnetic Resonance Imaging (fMRI) has shown brain activity alterations in individuals with a history of attempted suicide (SA) who are diagnosed with depression disorder (DD) or bipolar disorder (BD). However, patterns of spontaneous brain activity and their genetic correlations need further investigation.

METHODS

A voxel-based meta-analysis of 19 studies including 26 datasets, involving 742 patients with a history of SA and 978 controls (both nonsuicidal patients and healthy controls) was conducted. We examined fMRI changes in SA patients and analyzed the association between these changes and gene expression profiles using data from the Allen Human Brain Atlas by partial least squares regression analysis.

RESULTS

SA patients demonstrated increased spontaneous brain activity in several brain regions including the bilateral inferior temporal gyrus, hippocampus, fusiform gyrus, and right insula, and decreased activity in areas like the bilateral paracentral lobule and inferior frontal gyrus. Additionally, 5,077 genes were identified, exhibiting expression patterns associated with SA-related fMRI alterations. Functional enrichment analyses demonstrated that these SA-related genes were enriched for biological functions including glutamatergic synapse and mitochondrial structure. Concurrently, specific expression analyses showed that these genes were specifically expressed in the brain tissue, in neurons cells, and during early developmental periods.

CONCLUSION

Our findings suggest a neurobiological basis for fMRI abnormalities in SA patients with DD or BD, potentially guiding future genetic and therapeutic research.

Synaptic signaling modeled by functional connectivity predicts metabolic demands of the human brain

PURPOSE

The human brain is characterized by interacting large-scale functional networks fueled by glucose metabolism. Since former studies could not sufficiently clarify how these functional connections shape glucose metabolism, we aimed to provide a neurophysiologically-based approach.

METHODS

51 healthy volunteers underwent simultaneous PET/MRI to obtain BOLD functional connectivity and [18F]FDG glucose metabolism. These multimodal imaging proxies of fMRI and PET were combined in a whole-brain extension of metabolic connectivity mapping. Specifically, functional connectivity of all brain regions were used as input to explain glucose metabolism of a given target region. This enabled the modeling of postsynaptic energy demands by incoming signals from distinct brain regions.

RESULTS

Functional connectivity input explained a substantial part of metabolic demands but with pronounced regional variations (34 - 76%). During cognitive task performance this multimodal association revealed a shift to higher network integration compared to resting state. In healthy aging, a dedifferentiation (decreased segregated/modular structure of the brain) of brain networks during rest was observed. Furthermore, by including data from mRNA maps, [11C]UCB-J synaptic density and aerobic glycolysis (oxygen-to-glucose index from PET data), we show that whole-brain functional input reflects non-oxidative, on-demand metabolism of synaptic signaling. The metabolically-derived directionality of functional inputs further marked them as top-down predictions. In addition, the approach uncovered formerly hidden networks with superior efficiency through metabolically informed network partitioning.

CONCLUSIONS

Applying multimodal imaging, we decipher a crucial part of the metabolic and neurophysiological basis of functional connections in the brain as interregional on-demand synaptic signaling fueled by anaerobic metabolism. The observed task- and age-related effects indicate promising future applications to characterize human brain function and clinical alterations.

Morphological Brain Networks of White Matter: Mapping, Evaluation, Characterization, and Application

Although white matter (WM) accounts for nearly half of adult brain, its wiring diagram is largely unknown. Here, an approach is developed to construct WM networks by estimating interregional morphological similarity based on structural magnetic resonance imaging. It is found that morphological WM networks showed nontrivial topology, presented good-to-excellent test-retest reliability, accounted for phenotypic interindividual differences in cognition, and are under genetic control. Through integration with multimodal and multiscale data, it is further showed that morphological WM networks are able to predict the patterns of hamodynamic coherence, metabolic synchronization, gene co-expression, and chemoarchitectonic covariance, and associated with structural connectivity. Moreover, the prediction followed WM functional connectomic hierarchy for the hamodynamic coherence, is related to genes enriched in the forebrain neuron development and differentiation for the gene co-expression, and is associated with serotonergic system-related receptors and transporters for the chemoarchitectonic covariance. Finally, applying this approach to multiple sclerosis and neuromyelitis optica spectrum disorders, it is found that both diseases exhibited morphological dysconnectivity, which are correlated with clinical variables of patients and are able to diagnose and differentiate the diseases. Altogether, these findings indicate that morphological WM networks provide a reliable and biologically meaningful means to explore WM architecture in health and disease.

Transcriptomics : Approaches to Quantifying Gene Expression and Their Application to Studying the Human Brain

To date, the field of transcriptomics has been characterized by rapid methods development and technological advancement, with new technologies continuously rendering older ones obsolete.This chapter traces the evolution of approaches to quantifying gene expression and provides an overall view of the current state of the field of transcriptomics, its applications to the study of the human brain, and its place in the broader emerging multiomics landscape.

Genes associated with cortical thickness alterations in behavioral addiction

Behavioral addiction (BA) is a conceptually new addictive phenotype characterized by compulsive reward-seeking behaviors despite adverse consequences. Currently, its underlying neurogenetic mechanism remains unclear. Here, this study aimed to investigate the association between cortical thickness (CTh) and genetic phenotypes in BA. We conducted a systematic search in five databases and extracted gene expression data from the Allen Human Brain Atlas. Meta-analysis of 10 studies (343 addicted individuals and 355 controls) revealed that the BA group showed thinner CTh in the precuneus, postcentral gyrus, orbital-frontal cortex, and dorsolateral prefrontal cortex (P < 0.005). Meta-regression showed that the CTh in the precuneus and postcentral gyrus were negatively associated with the addiction severity (P < 0.0005). More importantly, the CTh phenotype of BA was spatially correlated with the expression of 12 genes (false discovery rate [FDR] < 0.05), and the dopamine D2 receptor had the highest correlation (rho = 0.55). Gene enrichment analysis further revealed that the 12 genes were involved in the biological processes of behavior regulation and response to stimulus (FDR < 0.05). In conclusion, our findings demonstrated the thinner CTh in cognitive control-related brain areas in BA, which could be associated with the expression of genes involving dopamine metabolism and behavior regulation.

Association between hearing ability and cortical morphology in the elderly: multiparametric mapping, cognitive relevance, and neurobiological underpinnings

BACKGROUND

Hearing impairment is a common condition in the elderly. However, a comprehensive understanding of its neural correlates is still lacking.

METHODS

We recruited 284 elderly adults who underwent structural MRI, magnetic resonance spectroscopy, audiometry, and cognitive assessments. Individual hearing abilities indexed by pure tone average (PTA) were correlated with multiple structural MRI-derived cortical morphological indices. For regions showing significant correlations, mediation analyses were performed to examine their role in the relationship between hearing ability and cognitive function. Finally, the correlation maps between hearing ability and cortical morphology were linked with publicly available connectomic gradient, transcriptomic, and neurotransmitter maps.

FINDINGS

Poorer hearing was related to cortical thickness (CT) reductions in widespread regions and gyrification index (GI) reductions in the right Area 52 and Insular Granular Complex. The GI in the right Area 52 mediated the relationship between hearing ability and executive function. This mediating effect was further modulated by glutamate and N-acetylaspartate levels in the right auditory region. The PTA-CT correlation map followed microstructural connectomic hierarchy, were related to genes involved in certain biological processes (e.g., glutamate metabolic process), cell types (e.g., excitatory neurons and astrocytes), and developmental stages (i.e., childhood to young adulthood), and covaried with dopamine receptor 1, dopamine transporter, and fluorodopa. The PTA-GI correlation map was related to 5-hydroxytryptamine receptor 2a.

INTERPRETATION

Poorer hearing is associated with cortical thinning and folding reductions, which may be engaged in the relationship between hearing impairment and cognitive decline in the elderly and have different neurobiological substrates.

FUNDING

See the Acknowledgements section.

Unraveling the molecular relevance of brain phenotypes: A comparative analysis of null models and test statistics

Correlating transcriptional profiles with imaging-derived phenotypes has the potential to reveal possible molecular architectures associated with cognitive functions, brain development and disorders. Competitive null models built by resampling genes and self-contained null models built by spinning brain regions, along with varying test statistics, have been used to determine the significance of transcriptional associations. However, there has been no systematic evaluation of their performance in imaging transcriptomics analyses. Here, we evaluated the performance of eight different test statistics (mean, mean absolute value, mean squared value, max mean, median, Kolmogorov-Smirnov (KS), Weighted KS and the number of significant correlations) in both competitive null models and self-contained null models. Simulated brain maps (n = 1,000) and gene sets (n = 500) were used to calculate the probability of significance (Psig) for each statistical test. Our results suggested that competitive null models may result in false positive results driven by co-expression within gene sets. Furthermore, we demonstrated that the self-contained null models may fail to account for distribution characteristics (e.g., bimodality) of correlations between all available genes and brain phenotypes, leading to false positives. These two confounding factors interacted differently with test statistics, resulting in varying outcomes. Specifically, the sign-sensitive test statistics (i.e., mean, median, KS, Weighted KS) were influenced by co-expression bias in the competitive null models, while median and sign-insensitive test statistics were sensitive to the bimodality bias in the self-contained null models. Additionally, KS-based statistics produced conservative results in the self-contained null models, which increased the risk of false negatives. Comprehensive supplementary analyses with various configurations, including realistic scenarios, supported the results. These findings suggest utilizing sign-insensitive test statistics such as mean absolute value, max mean in the competitive null models and the mean as the test statistic for the self-contained null models. Additionally, adopting the confounder-matched (e.g., coexpression-matched) null models as an alternative to standard null models can be a viable strategy. Overall, the present study offers insights into the selection of statistical tests for imaging transcriptomics studies, highlighting areas for further investigation and refinement in the evaluation of novel and commonly used tests.

Spatial transcriptomic brain imaging reveals the effects of immunomodulation therapy on specific regional brain cells in a mouse dementia model

Increasing evidence of brain-immune crosstalk raises expectations for the efficacy of novel immunotherapies in Alzheimer's disease (AD), but the lack of methods to examine brain tissues makes it difficult to evaluate therapeutics. Here, we investigated the changes in spatial transcriptomic signatures and brain cell types using the 10x Genomics Visium platform in immune-modulated AD models after various treatments. To proceed with an analysis suitable for barcode-based spatial transcriptomics, we first organized a workflow for segmentation of neuroanatomical regions, establishment of appropriate gene combinations, and comprehensive review of altered brain cell signatures. Ultimately, we investigated spatial transcriptomic changes following administration of immunomodulators, NK cell supplements and an anti-CD4 antibody, which ameliorated behavior impairment, and designated brain cells and regions showing probable associations with behavior changes. We provided the customized analytic pipeline into an application named STquantool. Thus, we anticipate that our approach can help researchers interpret the real action of drug candidates by simultaneously investigating the dynamics of all transcripts for the development of novel AD therapeutics.

Temporal dissociation between local and global functional adaptations of the maternal brain to childbirth: a longitudinal assessment

The maternal brain undergoes significant reorganization during birth and the postpartum period. However, the temporal dynamics of these changes remain unclear. Using resting-state functional magnetic resonance imaging, we report on local and global brain function alterations in 75 mothers in their first postpartum week, compared to 23 nulliparous women. In a subsample followed longitudinally for the next six months, we observed a temporal and spatial dissociation between changes observed at baseline (cluster mass permutation: pFWE < 0.05). Local activity and connectivity changes in widespread neocortical regions persisted throughout the studied time period (ANCOVAs vs. controls: pFDR < 0.05), with preliminary evidence linking these alterations to behavioral and psychological adaptations (interaction effect with postpartum time: uncorrected p < 0.05). In contrast, the initially reduced whole-brain connectivity of putamen-centered subcortical areas returned to control levels within six to nine weeks postpartum (linear and quadratic mixed linear models: pFDR < 0.05). The whole-brain spatial colocalization with hormone receptor distributions (Spearman correlations: pFDR < 0.05) and preliminary blood hormone associations (interaction effect with postpartum time: uncorrected p < 0.05) suggested that the postpartum restoration of progesterone levels may underlie this rapid normalization. These observations enhance our understanding of healthy maternal brain function, contributing to the identification of potential markers for pathological postpartum adaptation processes, which in turn could underlie postpartum psychiatric disorders.

Disproportionate neuroanatomical effects of DCC haploinsufficiency in adolescence compared with adulthood: links to dopamine, connectivity, covariance, and gene expression brain maps in mice

BACKGROUND

Critical adolescent neural refinement is controlled by the DCC (deleted in colorectal cancer) protein, a receptor for the netrin-1 guidance cue. We sought to describe the effects of reduced DCC on neuroanatomy in the adolescent and adult mouse brain.

METHODS

We examined neuronal connectivity, structural covariance, and molecular processes in a DCC-haploinsufficient mouse model, compared with wild-type mice, using new, custom analytical tools designed to leverage publicly available databases from the Allen Institute.

RESULTS

We included 11 DCC-haploinsufficient mice and 16 wild-type littermates. Neuroanatomical effects of DCC haploinsufficiency were more severe in adolescence than adulthood and were largely restricted to the mesocorticolimbic dopamine system. The latter finding was consistent whether we identified the regions of the mesocorticolimbic dopamine system a priori or used connectivity data from the Allen Brain Atlas to determine de novo where these dopamine axons terminated. Covariance analyses found that DCC haploinsufficiency disrupted the coordinated development of the brain regions that make up the mesocorticolimbic dopamine system. Gene expression maps pointed to molecular processes involving the expression of DCC, UNC5C (encoding DCC's co-receptor), and NTN1 (encoding its ligand, netrin-1) as underlying our structural findings.

LIMITATIONS

Our study involved a single sex (males) at only 2 ages.

CONCLUSION

The neuroanatomical phenotype of DCC haploinsufficiency described in mice parallels that observed in DCC-haploinsufficient humans. It is critical to understand the DCC-haploinsufficient mouse as a clinically relevant model system.

Childhood maltreatment influences adult brain structure through its effects on immune, metabolic, and psychosocial factors

Childhood maltreatment (CM) leads to a lifelong susceptibility to mental ill-health which might be reflected by its effects on adult brain structure, perhaps indirectly mediated by its effects on adult metabolic, immune, and psychosocial systems. Indexing these systemic factors via body mass index (BMI), C-reactive protein (CRP), and rates of adult trauma (AT), respectively, we tested three hypotheses: (H1) CM has direct or indirect effects on adult trauma, BMI, and CRP; (H2) adult trauma, BMI, and CRP are all independently related to adult brain structure; and (H3) childhood maltreatment has indirect effects on adult brain structure mediated in parallel by BMI, CRP, and AT. Using path analysis and data from N = 116,887 participants in UK Biobank, we find that CM is related to greater BMI and AT levels, and that these two variables mediate CM's effects on CRP [H1]. Regression analyses on the UKB MRI subsample (N = 21,738) revealed that greater CRP and BMI were both independently related to a spatially convergent pattern of cortical effects (Spearman's ρ = 0.87) characterized by fronto-occipital increases and temporo-parietal reductions in thickness. Subcortically, BMI was associated with greater volume, AT with lower volume and CPR with effects in both directions [H2]. Finally, path models indicated that CM has indirect effects in a subset of brain regions mediated through its direct effects on BMI and AT and indirect effects on CRP [H3]. Results provide evidence that childhood maltreatment can influence brain structure decades after exposure by increasing individual risk toward adult trauma, obesity, and inflammation.

Homotopic functional connectivity disruptions in schizophrenia and their associated gene expression

It has been revealed that abnormal voxel-mirrored homotopic connectivity (VMHC) is present in patients with schizophrenia, yet there are inconsistencies in the relevant findings. Moreover, little is known about their association with brain gene expression profiles. In this study, transcription-neuroimaging association analyses using gene expression data from Allen Human Brain Atlas and case-control VMHC differences from both the discovery (meta-analysis, including 9 studies with a total of 386 patients and 357 controls) and replication (separate group-level comparisons within two datasets, including a total of 258 patients and 287 controls) phases were performed to identify genes associated with VMHC alterations. Enrichment analyses were conducted to characterize the biological functions and specific expression of identified genes, and Neurosynth decoding analysis was performed to examine the correlation between cognitive-related processes and VMHC alterations in schizophrenia. In the discovery and replication phases, patients with schizophrenia exhibited consistent VMHC changes compared to controls, which were correlated with a series of cognitive-related processes; meta-regression analysis revealed that illness duration was negatively correlated with VMHC abnormalities in the cerebellum and postcentral/precentral gyrus. The abnormal VMHC patterns were stably correlated with 1287 genes enriched for fundamental biological processes like regulation of cell communication, nervous system development, and cell communication. In addition, these genes were overexpressed in astrocytes and immune cells, enriched in extensive cortical regions and wide developmental time windows. The present findings may contribute to a more comprehensive understanding of the molecular mechanisms underlying VMHC alterations in patients with schizophrenia.

A latent clinical-anatomical dimension relating metabolic syndrome to brain structure and cognition

The link between metabolic syndrome (MetS) and neurodegenerative as well as cerebrovascular conditions holds substantial implications for brain health in at-risk populations. This study elucidates the complex relationship between MetS and brain health by conducting a comprehensive examination of cardiometabolic risk factors, brain morphology, and cognitive function in 40,087 individuals. Multivariate, data-driven statistics identified a latent dimension linking more severe MetS to widespread brain morphological abnormalities, accounting for up to 71% of shared variance in the data. This dimension was replicable across sub-samples. In a mediation analysis, we could demonstrate that MetS-related brain morphological abnormalities mediated the link between MetS severity and cognitive performance in multiple domains. Employing imaging transcriptomics and connectomics, our results also suggest that MetS-related morphological abnormalities are linked to the regional cellular composition and macroscopic brain network organization. By leveraging extensive, multi-domain data combined with a dimensional stratification approach, our analysis provides profound insights into the association of MetS and brain health. These findings can inform effective therapeutic and risk mitigation strategies aimed at maintaining brain integrity.

Global-brain functional connectivity related with trait anxiety and its association with neurotransmitters and gene expression profiles

BACKGROUND

Numerous studies have explored the neural correlates of trait anxiety, a predisposing factor for several stress-related disorders. However, the findings from previous studies are inconsistent, which might be due to the limited regions of interest (ROI). A recent approach, named global-brain functional connectivity (GBC), has been demonstrated to address the shortcomings of ROI-based analysis. Furthermore, research on the transcriptome-connectome association has provided an approach to link the microlevel transcriptome profile with the macroscale brain network. In this paper, we aim to explore the neurobiology of trait anxiety with an imaging transcriptomic approach using GBC, biological neurotransmitters, and transcriptome profiles.

METHODS

Using a sample of resting-state fMRI data, we investigated trait anxiety-related alteration in GBC. We further used behavioral analysis, spatial correlation analysis, and postmortem gene expression to separately assess the cognitive functions, neurotransmitters, and transcriptional profiles related to alteration in GBC in individuals with trait anxiety.

RESULTS

GBC values in the ventromedial prefrontal cortex and the precuneus were negatively correlated with levels of trait anxiety. This alteration was correlated with behavioral terms including social cognition, emotion, and memory. A strong association was revealed between trait anxiety-related alteration in GBC and neurotransmitters, including dopaminergic, serotonergic, GABAergic, and glutamatergic systems in the ventromedial prefrontal cortex and the precuneus. The transcriptional profiles explained the functional connectivity, with correlated genes enriched in transmembrane signaling.

LIMITATIONS

Several limitations should be taken into account in this research. For example, future research should consider using some different approaches based on dynamic or task-based functional connectivity analysis, include more neurotransmitter receptors, additional gene expression data from different samples or more genes related to other stress-related disorders. Meanwhile, it is of great significance to include a larger sample size of individuals with a diagnosis of major depression disorder or other disorders for analysis and comparison and apply stricter multiple-comparison correction and threshold settings in future research.

CONCLUSIONS

Our research employed multimodal data to investigate GBC in the context of trait anxiety and to establish its associations with neurotransmitters and transcriptome profiles. This approach may improve understanding of the neural mechanism, together with the biological and molecular genetic foundations of GBC in trait anxiety.

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.

Molecular signatures of attention networks

Attention network theory proposes three distinct types of attention-alerting, orienting, and control-that are supported by separate brain networks and modulated by different neurotransmitters, that is, norepinephrine, acetylcholine, and dopamine. Here, we explore the extent of cortical, genetic, and molecular dissociation of these three attention systems using multimodal neuroimaging. We evaluated the spatial overlap between fMRI activation maps from the attention network test (ANT) and cortex-wide gene expression data from the Allen Human Brain Atlas. The goal was to identify genes associated with each of the attention networks in order to determine whether specific groups of genes were co-expressed with the corresponding attention networks. Furthermore, we analyzed publicly available PET-maps of neurotransmitter receptors and transporters to investigate their spatial overlap with the attention networks. Our analyses revealed a substantial number of genes (3871 for alerting, 6905 for orienting, 2556 for control) whose cortex-wide expression co-varied with the activation maps, prioritizing several molecular functions such as the regulation of protein biosynthesis, phosphorylation, and receptor binding. Contrary to the hypothesized associations, the ANT activation maps neither aligned with the distribution of norepinephrine, acetylcholine, and dopamine receptor and transporter molecules, nor with transcriptomic profiles that would suggest clearly separable networks. Independence of the attention networks appeared additionally constrained by a high level of spatial dependency between the network maps. Future work may need to reconceptualize the attention networks in terms of their segregation and reevaluate the presumed independence at the neural and neurochemical level.

Neural Correlates of Early-Life Urbanization and Their Spatial Relationships with Gene Expression, Neurotransmitter, and Behavioral Domain Atlases

Previous neuroimaging research has established associations between urban exposure during early life and alterations in brain function and structure. However, the molecular mechanisms and behavioral relevance of these associations remain largely unknown. Here, we aimed to address this question using a combined analysis of multimodal data. Initially, we calculated amplitude of low-frequency fluctuations (ALFF) and gray matter volume (GMV) using resting-state functional and structural MRI to investigate their associations with early-life urbanization in a large sample of 511 healthy young adults. Then, we examined the spatial relationships of the identified neural correlates of early-life urbanization with gene expression, neurotransmitter, and behavioral domain atlases. Results showed that higher early-life urbanization scores were correlated with increased ALFF of the right fusiform gyrus and decreased GMV of the left dorsal medial prefrontal cortex and left precuneus. Remarkably, the identified neural correlates of early-life urbanization were spatially correlated with expression of gene categories primarily involving immune system process, signal transduction, and cellular metabolic process. Concurrently, there were significant associations between the neural correlates and specific neurotransmitter systems including dopamine, acetylcholine, and serotonin. Finally, we found that the ALFF correlates were associated with behavioral terms including "perception," "sensory," "cognitive control," and "reasoning." Apart from expanding existing knowledge of early-life urban environmental risk for mental disorders and health in general, our findings may contribute to an emerging framework for integrating social science, neuroscience, genetics, and public policy to respond to the major health challenge of world urbanization.

Macroscale intrinsic dynamics are associated with microcircuit function in focal and generalized epilepsies

Epilepsies are a group of neurological disorders characterized by abnormal spontaneous brain activity, involving multiscale changes in brain functional organizations. However, it is not clear to what extent the epilepsy-related perturbations of spontaneous brain activity affect macroscale intrinsic dynamics and microcircuit organizations, that supports their pathological relevance. We collect a sample of patients with temporal lobe epilepsy (TLE) and genetic generalized epilepsy with tonic-clonic seizure (GTCS), as well as healthy controls. We extract massive temporal features of fMRI BOLD time-series to characterize macroscale intrinsic dynamics, and simulate microcircuit neuronal dynamics used a large-scale biological model. Here we show whether macroscale intrinsic dynamics and microcircuit dysfunction are differed in epilepsies, and how these changes are linked. Differences in macroscale gradient of time-series features are prominent in the primary network and default mode network in TLE and GTCS. Biophysical simulations indicate reduced recurrent connection within somatomotor microcircuits in both subtypes, and even more reduced in GTCS. We further demonstrate strong spatial correlations between differences in the gradient of macroscale intrinsic dynamics and microcircuit dysfunction in epilepsies. These results emphasize the impact of abnormal neuronal activity on primary network and high-order networks, suggesting a systematic abnormality of brain hierarchical organization.

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.

The Variegation of Human Brain Vulnerability to Rare Genetic Disorders and Convergence With Behaviorally Defined Disorders

BACKGROUND

Diverse gene dosage disorders (GDDs) increase risk for psychiatric impairment, but characterization of GDD effects on the human brain has so far been piecemeal, with few simultaneous analyses of multiple brain features across different GDDs.

METHODS

Here, through multimodal neuroimaging of 3 aneuploidy syndromes (XXY [total n = 191, 92 control participants], XYY [total n = 81, 47 control participants], and trisomy 21 [total n = 69, 41 control participants]), we systematically mapped the effects of supernumerary X, Y, and chromosome 21 dosage across a breadth of 15 different macrostructural, microstructural, and functional imaging-derived phenotypes (IDPs).

RESULTS

The results revealed considerable diversity in cortical changes across GDDs and IDPs. This variegation of IDP change underlines the limitations of studying GDD effects unimodally. Integration across all IDP change maps revealed highly distinct architectures of cortical change in each GDD along with partial coalescence onto a common spatial axis of cortical vulnerability that is evident in all 3 GDDs. This common axis shows strong alignment with shared cortical changes in behaviorally defined psychiatric disorders and is enriched for specific molecular and cellular signatures.

CONCLUSIONS

Use of multimodal neuroimaging data in 3 aneuploidies indicates that different GDDs impose unique fingerprints of change in the human brain that differ widely depending on the imaging modality that is being considered. Embedded in this variegation is a spatial axis of shared multimodal change that aligns with shared brain changes across psychiatric disorders and therefore represents a major high-priority target for future translational research in neuroscience.

Model-based hypothesis tests for the causal mediation of semi-competing risks

Analyzing the causal mediation of semi-competing risks has become important in medical research. Semi-competing risks refers to a scenario wherein an intermediate event may be censored by a primary event but not vice versa. Causal mediation analyses decompose the effect of an exposure on the primary outcome into an indirect (mediation) effect: an effect mediated through a mediator, and a direct effect: an effect not through the mediator. Here we proposed a model-based testing procedure to examine the indirect effect of the exposure on the primary event through the intermediate event. Under the counterfactual outcome framework, we defined a causal mediation effect using counting process. To assess statistical evidence for the mediation effect, we proposed two tests: an intersection-union test (IUT) and a weighted log-rank test (WLR). The test statistic was developed from a semi-parametric estimator of the mediation effect using a Cox proportional hazards model for the primary event and a series of logistic regression models for the intermediate event. We built a connection between the IUT and WLR. Asymptotic properties of the two tests were derived, and the IUT was determined to be a size [Formula: see text] test and statistically more powerful than the WLR. In numerical simulations, both the model-based IUT and WLR can properly adjust for confounding covariates, and the Type I error rates of the proposed methods are well protected, with the IUT being more powerful than the WLR. Our methods demonstrate the strongly significant effects of hepatitis B or C on the risk of liver cancer mediated through liver cirrhosis incidence in a prospective cohort study. The proposed method is also applicable to surrogate endpoint analyses in clinical trials.

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.

A latent clinical-anatomical dimension relating metabolic syndrome to brain structure and cognition

The link between metabolic syndrome (MetS) and neurodegenerative as well cerebrovascular conditions holds substantial implications for brain health in at-risk populations. This study elucidates the complex relationship between MetS and brain health by conducting a comprehensive examination of cardiometabolic risk factors, cortical morphology, and cognitive function in 40,087 individuals. Multivariate, data-driven statistics identified a latent dimension linking more severe MetS to widespread brain morphological abnormalities, accounting for up to 71% of shared variance in the data. This dimension was replicable across sub-samples. In a mediation analysis we could demonstrate that MetS-related brain morphological abnormalities mediated the link between MetS severity and cognitive performance in multiple domains. Employing imaging transcriptomics and connectomics, our results also suggest that MetS-related morphological abnormalities are linked to the regional cellular composition and macroscopic brain network organization. By leveraging extensive, multi-domain data combined with a dimensional stratification approach, our analysis provides profound insights into the association of MetS and brain health. These findings can inform effective therapeutic and risk mitigation strategies aimed at maintaining brain integrity.

A bipolar taxonomy of adult human brain sulcal morphology related to timing of fetal sulcation and trans-sulcal gene expression gradients

We developed a computational pipeline (now provided as a resource) for measuring morphological similarity between cortical surface sulci to construct a sulcal phenotype network (SPN) from each magnetic resonance imaging (MRI) scan in an adult cohort (N=34,725; 45-82 years). Networks estimated from pairwise similarities of 40 sulci on 5 morphological metrics comprised two clusters of sulci, represented also by the bipolar distribution of sulci on a linear-to-complex dimension. Linear sulci were more heritable and typically located in unimodal cortex; complex sulci were less heritable and typically located in heteromodal cortex. Aligning these results with an independent fetal brain MRI cohort (N=228; 21-36 gestational weeks), we found that linear sulci formed earlier, and the earliest and latest-forming sulci had the least between-adult variation. Using high-resolution maps of cortical gene expression, we found that linear sulcation is mechanistically underpinned by trans-sulcal gene expression gradients enriched for developmental processes.

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.

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.

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.

Genes associated with spontaneous brain activity changes in clinically different patients with major depressive disorder: A transcription-neuroimaging association study

AIMS

The amplitude of low-frequency fluctuations (ALFF) of resting-state functional MRI signals is a reliable neuroimaging measure of spontaneous brain activity. Inconsistent ALFF alterations have been reported in major depressive disorder (MDD) possibly due to clinical heterogeneity. This study was designed to investigate clinically sensitive and insensitive genes associated with ALFF alterations in MDD and the potential mechanisms.

METHODS

Transcription-neuroimaging association analyses of case-control ALFF differences from two independent neuroimaging datasets with gene expression data from Allen Human Brain Atlas were performed to identify the two gene sets. Various enrichment analyses were conducted to characterize their preference in biological functions, cell types, temporal stages, and shared effects with other psychiatric disorders.

RESULTS

Compared with controls, first-episode and drug-naïve patients showed more extensive ALFF alterations than patients with varied clinical features. We identified 903 clinically sensitive genes and 633 clinically insensitive genes, and the former was enriched for genes with down-regulated expression in the cerebral cortex of MDD patients. Despite shared functions of cell communication, signaling, and transport, clinically sensitive genes were enriched for cell differentiation and development whereas clinically insensitive genes were for ion transport and synaptic signaling. Clinically sensitive genes showed enrichment for microglia and macrophage from childhood to young adulthood in contrast to clinically insensitive genes for neurons before early infancy. Clinically sensitive genes (15.2%) were less likely correlated with ALFF alterations in schizophrenia than clinically insensitive genes (66.8%), and both were not relevant to bipolar disorder and adult attention deficit and hyperactivity disorder based on a third independent neuroimaging dataset.

CONCLUSIONS

Present results provide novel insights into the molecular mechanisms of spontaneous brain activity changes in clinically different patients with MDD.

Molecular mechanisms underlying resting-state brain functional correlates of behavioral inhibition

Previous literature has established the presence of sex differences in behavioral inhibition as well as its neural substrates and related disease risk. However, there is limited evidence that speaks directly to the question of whether or not there are sex-dependent associations between behavioral inhibition and resting-state brain function and, if so, how they are modulated by the underlying molecular mechanisms. We computed functional connectivity density (FCD) using resting-state functional MRI data to examine their associations with behavioral inhibition ability measured using a Go/No-Go task across a large cohort of 510 healthy young adults. Then, we examined the spatial relationships of the FCD correlates of behavioral inhibition with gene expression and neurotransmitter atlases to explore their potential genetic architecture and neurochemical basis. A significant negative correlation between behavioral inhibition and FCD in the left superior parietal lobule was found in females but not males. Further spatial correlation analyses demonstrated that the identified neural correlates of behavioral inhibition were associated with expression of gene categories predominantly implicating essential components of the cerebral cortex (glial cell, neuron, axon, dendrite, and synapse) and ion channel activity, as well as were linked to the serotonergic system. Our findings may not only yield important insights into the molecular mechanisms underlying the female-specific neural substrates of behavioral inhibition, but also provide a critical context for understanding how biological sex might contribute to variation in behavioral inhibition and its related disease risk.

The Neurogenetics of Functional Connectivity Alterations in Autism: Insights From Subtyping in 657 Individuals

BACKGROUND

There is little consensus and controversial evidence on anatomical alterations in the brains of people with autism spectrum disorder (ASD), due in part to the large heterogeneity present in ASD, which in turn is a major drawback for developing therapies. One strategy to characterize this heterogeneity in ASD is to cluster large-scale functional brain connectivity profiles.

METHODS

A subtyping approach based on consensus clustering of functional brain connectivity patterns was applied to a population of 657 autistic individuals with quality-assured neuroimaging data. We then used high-resolution gene transcriptomic data to characterize the molecular mechanism behind each subtype by performing enrichment analysis of the set of genes showing a high spatial similarity with the profiles of functional connectivity alterations between each subtype and a group of typically developing control participants.

RESULTS

Two major stable subtypes were found: subtype 1 exhibited hypoconnectivity (less average connectivity than typically developing control participants) and subtype 2, hyperconnectivity. The 2 subtypes did not differ in structural imaging metrics in any of the analyzed regions (68 cortical and 14 subcortical) or in any of the behavioral scores (including IQ, Autism Diagnostic Interview, and Autism Diagnostic Observation Schedule). Finally, only subtype 2, comprising about 43% of ASD participants, led to significant enrichments after multiple testing corrections. Notably, the dominant enrichment corresponded to excitation/inhibition imbalance, a leading well-known primary mechanism in the pathophysiology of ASD.

CONCLUSIONS

Our results support a link between excitation/inhibition imbalance and functional connectivity alterations, but only in one ASD subtype, overall characterized by brain hyperconnectivity and major alterations in somatomotor and default mode networks.

Indirect evidence for altered dopaminergic neurotransmission in very premature-born adults

While animal models indicate altered brain dopaminergic neurotransmission after premature birth, corresponding evidence in humans is scarce due to missing molecular imaging studies. To overcome this limitation, we studied dopaminergic neurotransmission changes in human prematurity indirectly by evaluating the spatial co-localization of regional alterations in blood oxygenation fluctuations with the distribution of adult dopaminergic neurotransmission. The study cohort comprised 99 very premature-born (<32 weeks of gestation and/or birth weight below 1500 g) and 107 full-term born young adults, being assessed by resting-state functional MRI (rs-fMRI) and IQ testing. Normative molecular imaging dopamine neurotransmission maps were derived from independent healthy control groups. We computed the co-localization of local (rs-fMRI) activity alterations in premature-born adults with respect to term-born individuals to different measures of dopaminergic neurotransmission. We performed selectivity analyses regarding other neuromodulatory systems and MRI measures. In addition, we tested if the strength of the co-localization is related to perinatal measures and IQ. We found selectively altered co-localization of rs-fMRI activity in the premature-born cohort with dopamine-2/3-receptor availability in premature-born adults. Alterations were specific for the dopaminergic system but not for the used MRI measure. The strength of the co-localization was negatively correlated with IQ. In line with animal studies, our findings support the notion of altered dopaminergic neurotransmission in prematurity which is associated with cognitive performance.

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.

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.

Progressive trajectories of schizophrenia across symptoms, genes, and the brain

BACKGROUND

Schizophrenia is characterized by complex psychiatric symptoms and unclear pathological mechanisms. Most previous studies have focused on the morphological changes that occur over the development of the disease; however, the corresponding functional trajectories remain unclear. In the present study, we aimed to explore the progressive trajectories of patterns of dysfunction after diagnosis.

METHODS

Eighty-six patients with schizophrenia and 120 healthy controls were recruited as the discovery dataset. Based on multiple functional indicators of resting-state brain functional magnetic resonance imaging, we conducted a duration-sliding dynamic analysis framework to investigate trajectories in association with disease progression. Neuroimaging findings were associated with clinical symptoms and gene expression data from the Allen Human Brain Atlas database. A replication cohort of patients with schizophrenia from the University of California, Los Angeles, was used as the replication dataset for the validation analysis.

RESULTS

Five stage-specific phenotypes were identified. A symptom trajectory was characterized by positive-dominated, negative ascendant, negative-dominated, positive ascendant, and negative surpassed stages. Dysfunctional trajectories from primary and subcortical regions to higher-order cortices were recognized; these are associated with abnormal external sensory gating and a disrupted internal excitation-inhibition equilibrium. From stage 1 to stage 5, the importance of neuroimaging features associated with behaviors gradually shifted from primary to higher-order cortices and subcortical regions. Genetic enrichment analysis identified that neurodevelopmental and neurodegenerative factors may be relevant as schizophrenia progresses and highlighted multiple synaptic systems.

CONCLUSIONS

Our convergent results indicate that progressive symptoms and functional neuroimaging phenotypes are associated with genetic factors in schizophrenia. Furthermore, the identification of functional trajectories complements previous findings of structural abnormalities and provides potential targets for drug and non-drug interventions in different stages of schizophrenia.

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.

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.

Molecular and network-level mechanisms explaining individual differences in autism spectrum disorder

The mechanisms underlying phenotypic heterogeneity in autism spectrum disorder (ASD) are not well understood. Using a large neuroimaging dataset, we identified three latent dimensions of functional brain network connectivity that predicted individual differences in ASD behaviors and were stable in cross-validation. Clustering along these three dimensions revealed four reproducible ASD subgroups with distinct functional connectivity alterations in ASD-related networks and clinical symptom profiles that were reproducible in an independent sample. By integrating neuroimaging data with normative gene expression data from two independent transcriptomic atlases, we found that within each subgroup, ASD-related functional connectivity was explained by regional differences in the expression of distinct ASD-related gene sets. These gene sets were differentially associated with distinct molecular signaling pathways involving immune and synapse function, G-protein-coupled receptor signaling, protein synthesis and other processes. Collectively, our findings delineate atypical connectivity patterns underlying different forms of ASD that implicate distinct molecular signaling mechanisms.

Adolescent neurodevelopment and psychopathology: The interplay between adversity exposure and genetic risk for accelerated brain ageing

In adulthood, stress exposure and genetic risk heighten psychological vulnerability by accelerating neurobiological senescence. To investigate whether molecular and brain network maturation processes play a similar role in adolescence, we analysed genetic, as well as longitudinal task neuroimaging (inhibitory control, incentive processing) and early life adversity (i.e., material deprivation, violence) data from the Adolescent Brain and Cognitive Development study (N = 980, age range: 9-13 years). Genetic risk was estimated separately for Major Depressive Disorder (MDD) and Alzheimer's Disease (AD), two pathologies linked to stress exposure and allegedly sharing a causal connection (MDD-to-AD). Adversity and genetic risk for MDD/AD jointly predicted functional network segregation patterns suggestive of accelerated (GABA-linked) visual/attentional, but delayed (dopamine [D2]/glutamate [GLU5R]-linked) somatomotor/association system development. A positive relationship between brain maturation and psychopathology emerged only among the less vulnerable adolescents, thereby implying that normatively maladaptive neurodevelopmental alterations could foster adjustment among the more exposed and genetically more stress susceptible youths. Transcriptomic analyses suggested that sensitivity to stress may underpin the joint neurodevelopmental effect of adversity and genetic risk for MDD/AD, in line with the proposed role of negative emotionality as a precursor to AD, likely to account for the alleged causal impact of MDD on dementia onset.

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.

Revealing the neurobiology underlying interpersonal neural synchronization with multimodal data fusion

Humans synchronize with one another to foster successful interactions. Here, we use a multimodal data fusion approach with the aim of elucidating the neurobiological mechanisms by which interpersonal neural synchronization (INS) occurs. Our meta-analysis of 22 functional magnetic resonance imaging and 69 near-infrared spectroscopy hyperscanning experiments (740 and 3721 subjects) revealed robust brain regional correlates of INS in the right temporoparietal junction and left ventral prefrontal cortex. Integrating this meta-analytic information with public databases, biobehavioral and brain-functional association analyses suggested that INS involves sensory-integrative hubs with functional connections to mentalizing and attention networks. On the molecular and genetic levels, we found INS to be associated with GABAergic neurotransmission and layer IV/V neuronal circuits, protracted developmental gene expression patterns, and disorders of neurodevelopment. Although limited by the indirect nature of phenotypic-molecular association analyses, our findings generate new testable hypotheses on the neurobiological basis of INS.

Bioinformatic prediction of the molecular links between Alzheimer's disease and diabetes mellitus

Background

Alzheimer's disease (AD) and type 2 diabetes mellitus (DM2) are chronic degenerative diseases with complex molecular processes that are potentially interconnected. The aim of this work was to predict the potential molecular links between AD and DM2 from different sources of biological information.

Materials and Methods

In this work, data mining of nine databases (DisGeNET, Ensembl, OMIM, Protein Data Bank, The Human Protein Atlas, UniProt, Gene Expression Omnibus, Human Cell Atlas, and PubMed) was performed to identify gene and protein information that was shared in AD and DM2. Next, the information was mapped to human protein-protein interaction (PPI) networks based on experimental data using the STRING web platform. Then, gene ontology biological process (GOBP) and pathway analyses with EnrichR showed its specific and shared biological process and pathway deregulations. Finally, potential biomarkers and drug targets were predicted with the Metascape platform.

Results

A total of 1,551 genes shared in AD and DM2 were identified. The highest average degree of nodes within the PPI was for DM2 (average = 2.97), followed by AD (average degree = 2.35). GOBP for AD was related to specific transcriptional and translation genetic terms occurring in neurons cells. The GOBP and pathway information for the association AD-DM2 were linked mainly to bioenergetics and cytokine signaling. Within the AD-DM2 association, 10 hub proteins were identified, seven of which were predicted to be present in plasma and exhibit pharmacological interaction with monoclonal antibodies in use, anticancer drugs, and flavonoid derivatives.

Conclusion

Our data mining and analysis strategy showed that there are a plenty of biological information based on experiments that links AD and DM2, which could provide a rational guide to design further diagnosis and treatment for AD and DM2.

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

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

Cortical profiles of numerous psychiatric disorders and normal development share a common pattern

The neurobiological bases of the association between development and psychopathology remain poorly understood. Here, we identify a shared spatial pattern of cortical thickness (CT) in normative development and several psychiatric and neurological disorders. Principal component analysis (PCA) was applied to CT of 68 regions in the Desikan-Killiany atlas derived from three large-scale datasets comprising a total of 41,075 neurotypical participants. PCA produced a spatially broad first principal component (PC1) that was reproducible across datasets. Then PC1 derived from healthy adult participants was compared to the pattern of CT differences associated with psychiatric and neurological disorders comprising a total of 14,886 cases and 20,962 controls from seven ENIGMA disease-related working groups, normative maturation and aging comprising a total of 17,697 scans from the ABCD Study® and the IMAGEN developmental study, and 17,075 participants from the ENIGMA Lifespan working group, as well as gene expression maps from the Allen Human Brain Atlas. Results revealed substantial spatial correspondences between PC1 and widespread lower CT observed in numerous psychiatric disorders. Moreover, the PC1 pattern was also correlated with the spatial pattern of normative maturation and aging. The transcriptional analysis identified a set of genes including KCNA2, KCNS1 and KCNS2 with expression patterns closely related to the spatial pattern of PC1. The gene category enrichment analysis indicated that the transcriptional correlations of PC1 were enriched to multiple gene ontology categories and were specifically over-represented starting at late childhood, coinciding with the onset of significant cortical maturation and emergence of psychopathology during the prepubertal-to-pubertal transition. Collectively, the present study reports a reproducible latent pattern of CT that captures interregional profiles of cortical changes in both normative brain maturation and a spectrum of psychiatric disorders. The pubertal timing of the expression of PC1-related genes implicates disrupted neurodevelopment in the pathogenesis of the spectrum of psychiatric diseases emerging during adolescence.

Radiomic features of gray matter in never-treated first-episode schizophrenia

Alterations of radiomic features (RFs) in gray matter are observed in schizophrenia, of which the results may be limited by small study samples and confounding effects of drug therapies. We tested for RFs alterations of gray matter in never-treated first-episode schizophrenia (NT-FES) patients and examined their associations with known gene expression profiles. RFs were examined in the first sample with 197 NT-FES and 178 healthy controls (HCs) and validated in the second independent sample (90 NT-FES and 74 HCs). One-year follow-up data were available from 87 patients to determine whether RFs were associated with treatment outcomes. Associations between identified RFs in NT-FES and gene expression profiles were evaluated. NT-FES exhibited alterations of 30 RFs, with the greatest involvement of microstructural heterogeneity followed by measures of brain region shape. The identified RFs were mainly located in the central executive network, frontal-temporal network, and limbic system. Two baseline RFs with the involvement of microstructural heterogeneity predicted treatment response with moderate accuracy (78% for the first sample, 70% for the second sample). Exploratory analyses indicated that RF alterations were spatially related to the expression of schizophrenia risk genes. In summary, the present findings link brain abnormalities in schizophrenia with molecular features and treatment response.

Regional gene expression patterns are associated with task-specific brain activation during reward and emotion processing measured with functional MRI

The exploration of the spatial relationship between gene expression profiles and task-evoked response patterns known to be altered in neuropsychiatric disorders, for example depression, can guide the development of more targeted therapies. Here, we estimated the correlation between human transcriptome data and two different brain activation maps measured with functional magnetic resonance imaging (fMRI) in healthy subjects. Whole-brain activation patterns evoked during an emotional face recognition task were associated with topological mRNA expression of genes involved in cellular transport. In contrast, fMRI activation patterns related to the acceptance of monetary rewards were associated with genes implicated in cellular localization processes, metabolism, translation, and synapse regulation. An overlap of these genes with risk genes from major depressive disorder genome-wide association studies revealed the involvement of the master regulators TCF4 and PAX6 in emotion and reward processing. Overall, the identification of stable relationships between spatial gene expression profiles and fMRI data may reshape the prospects for imaging transcriptomics studies.