Correlation between cortical gene expression and resting-state functional network centrality in healthy young adults

Dysfunction in the hierarchy of morphometric similarity network in Alzheimer's disease and its correlation with cognitive performance and gene expression profiles

BACKGROUND

Previous research has shown abnormal functional network gradients in Alzheimer's disease (AD). Structural network gradient is capable of capturing continuous changes in brain morphology and has the ability to elucidate the underlying processes of neurodevelopment. However, it remains unclear whether structural network gradients are altered in AD and what associations exist between these changes and cognitive function, and gene expression profiles.

METHODS

By constructing an individualized structural network gradient decomposition framework, we calculated the morphological similarity network (MSN) gradients for 404 subjects (186 AD patients and 218 normal controls). We investigated AD-related alterations in MSN gradients, along with the associations between MSN gradients and cognitive function, MSN topological properties, and gene expression profiles.

RESULTS

Our findings indicated that the principal MSN gradient alterations in AD were primarily characterized by an increase in the primary and secondary sensory cortices and a decrease in the association cortex 1. The primary and higher-order cortices exhibited opposite associations with cognition, including executive function, language skills, and memory processes. Moreover, the principal MSN gradients were found to significantly predict cognitive function in AD. The altered gradient pattern was 14.8% attributable to gene expression profiles, and the genes demonstrating the highest correlation are involved in metabolic activity and synaptic signaling.

CONCLUSIONS

Our results offered novel insights into the underlying mechanisms of structural brain network impairment in AD patients, enhancing our understanding of the neurobiological processes responsible for impaired cognition in patients with AD, and offering a new dimensional structural biomarker for AD.

Disparate and shared transcriptomic signatures associated with cortical atrophy in genetic behavioral variant frontotemporal degeneration

BACKGROUND

Cortical atrophy is a common manifestation in behavioral variant frontotemporal degeneration (bvFTD), exhibiting spatial heterogeneity across various genetic subgroups, which may be driven by distinct biological mechanisms.

METHODS

We employed an integrative imaging transcriptomics approach to identify both disparate and shared transcriptomic signatures associated with cortical thickness in bvFTD with C9orf72 repeat expansions or pathogenic variants in GRN or MAPT. Functional enrichment analyses were conducted on each gene list significantly associated with cortical thickness. Additionally, we mapped neurotransmitter receptor/transporter density maps to the cortical thickness maps, to uncover different correlation patterns for each genetic form. Furthermore, we examined whether the identified genes were enriched for pathology-related genes by using previously identified genes linked to TDP-43 positive neurons and genes associated with tau pathology.

RESULTS

For each genetic form of bvFTD, we identified cortical thickness signatures and gene sets associated with them. The cortical thickness associated genes for GRN-bvFTD were significantly involved in neurotransmitter system and circadian entrainment. The different patterns of spatial correlations between synaptic density and cortical thinning, further confirmed the critical role of neurotransmission and synaptic signaling in shaping brain structure, especially in the GRN-bvFTD group. Furthermore, we observed significant overlap between genes linked to TDP-43 pathology and the gene sets associated with cortical thickness in C9orf72-bvFTD and GRN-bvFTD but not the MAPT-bvFTD group providing specificity for our associations. C9orf72-bvFTD and GRN-bvFTD also shared genes displaying consistent directionality, with those exhibiting either positive or negative correlations with cortical thickness in C9orf72-bvFTD showing the same direction (positive or negative) in GRN-bvFTD. MAPT-bvFTD displayed more pronounced differences in transcriptomic signatures compared to the other two genetic forms. The genes that exhibited significantly positive or negative correlations with cortical thickness in MAPT-bvFTD showed opposing directionality in C9orf72-bvFTD and GRN-bvFTD.

CONCLUSIONS

Overall, this integrative transcriptomic approach identified several new shared and disparate genes associated with regional vulnerability with increased biological interpretation including overlap with synaptic density maps and pathologically-specific gene expression. These findings illuminated the intricate molecular underpinnings contributing to the heterogeneous nature of disease distribution in bvFTD with distinct genetic backgrounds.

Atypical dynamic neural configuration in autism spectrum disorder and its relationship to gene expression profiles

Although it is well recognized that autism spectrum disorder (ASD) is associated with atypical dynamic functional connectivity patterns, the dynamic changes in brain intrinsic activity over each time point and the potential molecular mechanisms associated with atypical dynamic temporal characteristics in ASD remain unclear. Here, we employed the Hidden Markov Model (HMM) to explore the atypical neural configuration at every scanning time point in ASD, based on resting-state functional magnetic resonance imaging (rs-fMRI) data from the Autism Brain Imaging Data Exchange. Subsequently, partial least squares regression and pathway enrichment analysis were employed to explore the potential molecular mechanism associated with atypical neural dynamics in ASD. 8 HMM states were inferred from rs-fMRI data. Compared to typically developing, individuals on the autism spectrum showed atypical state-specific temporal characteristics, including number of states and occurrences, mean life time and transition probability between states. Moreover, these atypical temporal characteristics could predict communication difficulties of ASD, and states assoicated with negative activation in default mode network and frontoparietal network, and positive activation in somatomotor network, ventral attention network, and limbic network, had higher predictive contribution. Furthermore, a total of 321 genes was revealed to be significantly associated with atypical dynamic brain states of ASD, and these genes are mainly enriched in neurodevelopmental pathways. Our study provides new insights into characterizing the atypical neural dynamics from a moment-to-moment perspective, and indicates a linkage between atypical neural configuration and gene expression in ASD.

Mechanisms underlying the spontaneous reorganization of depression network after stroke

Exploring the causal relationship between focal brain lesions and post-stroke depression (PSD) can provide therapeutic insights. However, a gap exists between causal and therapeutic information. Exploring post-stroke brain repair processes post-stroke could bridge this gap. We defined a depression network using the normative connectome and investigated the predictive capacity of lesion-induced network damage on depressive symptoms in discovery cohort of 96 patients, at baseline and six months post-stroke. Stepwise functional connectivity (SFC) was used to examine topological changes in the depression network over time to identify patterns of network reorganization. The predictive value of reorganization information was evaluated for follow-up symptoms in discovery and validation cohort 1 (22 worsening PSD patients) as well as for treatment responsiveness in validation cohort 2 (23 antidepressant-treated patients). We evaluated the consistency of significant reorganization areas with neuromodulation targets. Spatial correlations of network reorganization patterns with gene expression and neurotransmitter maps were analyzed. The predictive power of network damage for symptoms diminished at follow-up compared to baseline (Δadjusted R2 = -0.070, p < 0.001). Reorganization information effectively predicted symptoms at follow-up in the discovery cohort (adjust R2 = 0.217, 95 %CI: 0.010 to 0.431), as well as symptom exacerbation (r = 0.421, p = 0.033) and treatment responsiveness (r = 0.587, p = 0.012) in the validation cohorts. Regions undergoing significant reorganization overlapped with neuromodulatory targets known to be effective in treating depression. The reorganization of the depression network was associated with immune-inflammatory responses gene expressions and gamma-aminobutyric acid. Our findings may yield important insights into the repair mechanisms of PSD and provide a critical context for developing post-stroke treatment strategies.

Spatial transcriptomics reveals modulation of transcriptional networks across brain regions after auditory threat conditioning

Prior research has demonstrated genome-wide transcriptional changes related to fear and anxiety across species, often focusing on individual brain regions or cell types. However, the extent of gene expression differences across brain regions and how these changes interact at the level of transcriptional connectivity remains unclear. To address this, we performed spatial transcriptomics RNAseq analyses in an auditory threat conditioning paradigm in mice. We generated a spatial transcriptomic atlas of a coronal mouse brain section covering cortical and subcortical regions, corresponding to histologically defined regions. Our finding revealed widespread transcriptional responses across all brain regions examined, particularly in the medial and lateral habenula, and the choroid plexus. Network analyses highlighted altered transcriptional connectivity between cortical and subcortical regions, emphasizing the role of steroidogenic factor 1. These results provide new insights into the transcriptional networks involved in auditory threat conditioning, enhancing our understanding of molecular and neural mechanisms underlying fear and anxiety disorders.

Age-dependent functional development pattern in neonatal brain: An fMRI-based brain entropy study

The relationship between brain entropy (BEN) and early brain development has been established through animal studies. However, it remains unclear whether the BEN can be used to identify age-dependent functional changes in human neonatal brains and the genetic underpinning of the new neuroimaging marker remains to be elucidated. In this study, we analyzed resting-state fMRI data from the Developing Human Connectome Project, including 280 infants who were scanned at 37.5-43.5 weeks postmenstrual age. The BEN maps were calculated for each subject, and a voxel-wise analysis was conducted using a general linear model to examine the effects of age, sex, and preterm birth on BEN. Additionally, we evaluated the correlation between regional BEN and gene expression levels. Our results demonstrated that the BEN in the sensorimotor-auditory and association cortices, along the 'S-A' axis, was significantly positively correlated with postnatal age (PNA), and negatively correlated with gestational age (GA), respectively. Meanwhile, the BEN in the right rolandic operculum correlated significantly with both GA and PNA. Preterm-born infants exhibited increased BEN values in widespread cortical areas, particularly in the visual-motor cortex, when compared to term-born infants. Moreover, we identified five BEN-related genes (DNAJC12, FIG4, STX12, CETN2, and IRF2BP2), which were involved in protein folding, synaptic vesicle transportation and cell division. These findings suggest that the fMRI-based BEN can serve as an indicator of age-dependent brain functional development in human neonates, which may be influenced by specific genes.

Dynamic functional connectivity and gene expression correlates in temporal lobe epilepsy: insights from hidden markov models

BACKGROUD

Temporal lobe epilepsy (TLE) is associated with abnormal dynamic functional connectivity patterns, but the dynamic changes in brain activity at each time point remain unclear, as does the potential molecular mechanisms associated with the dynamic temporal characteristics of TLE.

METHODS

Resting-state functional magnetic resonance imaging (rs-fMRI) was acquired for 84 TLE patients and 35 healthy controls (HCs). The data was then used to conduct HMM analysis on rs-fMRI data from TLE patients and an HC group in order to explore the intricate temporal dynamics of brain activity in TLE patients with cognitive impairment (TLE-CI). Additionally, we aim to examine the gene expression profiles associated with the dynamic modular characteristics in TLE patients using the Allen Human Brain Atlas (AHBA) database.

RESULTS

Five HMM states were identified in this study. Compared with HCs, TLE and TLE-CI patients exhibited distinct changes in dynamics, including fractional occupancy, lifetimes, mean dwell time and switch rate. Furthermore, transition probability across HMM states were significantly different between TLE and TLE-CI patients (p < 0.05). The temporal reconfiguration of states in TLE and TLE-CI patients was associated with several brain networks (including the high-order default mode network (DMN), subcortical network (SCN), and cerebellum network (CN). Furthermore, a total of 1580 genes were revealed to be significantly associated with dynamic brain states of TLE, mainly enriched in neuronal signaling and synaptic function.

CONCLUSIONS

This study provides new insights into characterizing dynamic neural activity in TLE. The brain network dynamics defined by HMM analysis may deepen our understanding of the neurobiological underpinnings of TLE and TLE-CI, indicating a linkage between neural configuration and gene expression in TLE.

Gene expression and brain imaging association study reveals gene signatures in major depressive disorder

Major depressive disorder is often characterized by changes in the structure and function of the brain, which are influenced by modifications in gene expression profiles. How the depression-related genes work together within the scope of time and space to cause pathological changes remains unclear. By integrating the brain-wide gene expression data and imaging data in major depressive disorder, we identified gene signatures of major depressive disorder and explored their temporal-spatial expression specificity, network properties, function annotations and sex differences systematically. Based on correlation analysis with permutation testing, we found 345 depression-related genes significantly correlated with functional and structural alteration of brain images in major depressive disorder and separated them by directional effects. The genes with negative effect for grey matter density and positive effect for functional indices are enriched in downregulated genes in the post-mortem brain samples of patients with depression and risk genes identified by genome-wide association studies than genes with positive effect for grey matter density and negative effect for functional indices and control genes, confirming their potential association with major depressive disorder. By introducing a parameter of dispersion measure on the gene expression data of developing human brains, we revealed higher spatial specificity and lower temporal specificity of depression-related genes than control genes. Meanwhile, we found depression-related genes tend to be more highly expressed in females than males, which may contribute to the difference in incidence rate between male and female patients. In general, we found the genes with negative effect have lower network degree, more specialized function, higher spatial specificity, lower temporal specificity and more sex differences than genes with positive effect, indicating they may play different roles in the occurrence and development of major depressive disorder. These findings can enhance the understanding of molecular mechanisms underlying major depressive disorder and help develop tailored diagnostic and treatment strategies for patients of depression of different sex.

Disparate and shared transcriptomic signatures associated with cortical atrophy in genetic bvFTD

Cortical atrophy in behavioral variant frontotemporal degeneration (bvFTD) exhibits spatial heterogeneity across genetic subgroups, potentially driven by distinct biological mechanisms. Using an integrative imaging-transcriptomics approach, we identified disparate and shared transcriptomic signatures associated with cortical thickness in C9orf72 , GRN or MAPT -related bvFTD. Genes associated with cortical thinning in GRN -bvFTD were implicated in neurotransmission, further supported by mapping synaptic density maps to cortical thickness maps. Previously identified genes linked to TDP-43 positive neurons were significantly overlapped with genes associated with C9orf72 -bvFTD and GRN -bvFTD, but not MAPT -bvFTD providing specificity for our associations. C9orf72 -bvFTD and GRN -bvFTD shared genes displaying consistent directionality of correlations with cortical thickness, while MAPT -bvFTD displayed more pronounced differences in transcriptomic signatures with opposing directionality. Overall, we identified disparate and shared genes tied to regional vulnerability with increased biological interpretation including overlap with synaptic density maps and pathologically-specific gene expression, illuminating intricate molecular underpinnings contributing to heterogeneities in bvFTD.

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.

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.

Molecular basis underlying default mode network functional abnormalities in postpartum depression with and without anxiety

Although Postpartum depression (PPD) and PPD with anxiety (PPD-A) have been well characterized as functional disruptions within or between multiple brain systems, however, how to quantitatively delineate brain functional system irregularity and the molecular basis of functional abnormalities in PPD and PPD-A remains unclear. Here, brain sample entropy (SampEn), resting-state functional connectivity (RSFC), transcriptomic and neurotransmitter density data were used to investigate brain functional system irregularity, functional connectivity abnormalities and associated molecular basis for PPD and PPD-A. PPD-A exhibited higher SampEn in medial prefrontal cortex (MPFC) and posterior cingulate cortex (PPC) than healthy postnatal women (HPW) and PPD while PPD showed lower SampEn in PPC compared to HPW and PPD-A. The functional connectivity analysis with MPFC and PPC as seed areas revealed decreased functional couplings between PCC and paracentral lobule and between MPFC and angular gyrus in PPD compared to both PPD-A and HPW. Moreover, abnormal SampEn and functional connectivity were associated with estrogenic level and clinical symptoms load. Importantly, spatial association analyses between functional changes and transcriptome and neurotransmitter density maps revealed that these functional changes were primarily associated with synaptic signaling, neuron projection, neurotransmitter level regulation, amino acid metabolism, cyclic adenosine monophosphate (cAMP) signaling pathways, and neurotransmitters of 5-hydroxytryptamine (5-HT), norepinephrine, glutamate, dopamine and so on. These results reveal abnormal brain entropy and functional connectivities primarily in default mode network (DMN) and link these changes to transcriptome and neurotransmitters to establish the molecular basis for PPD and PPD-A for the first time. Our findings highlight the important role of DMN in neuropathology of PPD and PPD-A.

Molecular basis underlying changes of brain entropy and functional connectivity in major depressive disorders after electroconvulsive therapy

INTRODUCTION

Electroconvulsive therapy (ECT) is widely used for treatment-resistant depression. However, it is unclear whether/how ECT can be targeted to affect brain regions and circuits in the brain to dynamically regulate mood and cognition.

METHODS

This study used brain entropy (BEN) to measure the irregular levels of brain systems in 46 major depressive disorder (MDD) patients before and after ECT treatment. Functional connectivity (FC) was further adopted to reveal changes of functional couplings. Moreover, transcriptomic and neurotransmitter receptor data were used to reveal genetic and molecular basis of the changes of BEN and functional connectivities.

RESULTS

Compared to pretreatment, the BEN in the posterior cerebellar lobe (PCL) significantly decreased and FC between the PCL and the right temporal pole (TP) significantly increased in MDD patients after treatment. Moreover, we found that these changes of BEN and FC were closely associated with genes' expression profiles involved in MAPK signaling pathway, GABAergic synapse, and dopaminergic synapse and were significantly correlated with the receptor/transporter density of 5-HT, norepinephrine, glutamate, etc. CONCLUSION: These findings suggest that loops in the cerebellum and TP are crucial for ECT regulation of mood and cognition, which provides new evidence for the antidepressant effects of ECT and the potential molecular mechanism leading to cognitive impairment.

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.

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.

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.

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.

Trade-offs among cost, integration, and segregation in the human connectome

The human brain structural network is thought to be shaped by the optimal trade-off between cost and efficiency. However, most studies on this problem have focused on only the trade-off between cost and global efficiency (i.e., integration) and have overlooked the efficiency of segregated processing (i.e., segregation), which is essential for specialized information processing. Direct evidence on how trade-offs among cost, integration, and segregation shape the human brain network remains lacking. Here, adopting local efficiency and modularity as segregation factors, we used a multiobjective evolutionary algorithm to investigate this problem. We defined three trade-off models, which represented trade-offs between cost and integration (Dual-factor model), and trade-offs among cost, integration, and segregation (local efficiency or modularity; Tri-factor model), respectively. Among these, synthetic networks with optimal trade-off among cost, integration, and modularity (Tri-factor model [Q]) showed the best performance. They had a high recovery rate of structural connections and optimal performance in most network features, especially in segregated processing capacity and network robustness. Morphospace of this trade-off model could further capture the variation of individual behavioral/demographic characteristics in a domain-specific manner. Overall, our results highlight the importance of modularity in the formation of the human brain structural network and provide new insights into the original cost-efficiency trade-off hypothesis.

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.

Neutrophil-to-Lymphocyte Ratio Amplifies the Effects of Aging on Decrements in Grip Strength and Its Functional Neural Underpinnings

The neutrophil-to-lymphocyte ratio (NLR) is a trans-prognostic biomarker of physiologic stress and inflammation linked to muscle weakness in older adults. Generation of grip force coincides with sustained activity in the primary sensorimotor cortex (SM1). The current study investigates whether whole-brain functional connectivity, that is, degree centrality (CD) of SM1 relates to grip strength and whether both functional measures are predicted by advancing age as a function of the NLR. A structural regression model investigated the main and interactive effects of age and NLR on grip strength and CD of SM1 in 589 adults aged 21-85 years (M = 45.87, SD = 18.06). The model including the entire sample had a good fit (χ 2(4) = 1.63, p = .804). In individuals aged 50 years and older, age predicted lower grip strength and SM1 CD as a function of increasing NLR. In a model stratified by sex, the effect of age, NLR, and their interaction on grip strength are significant for older men but not older women. Analyses support CD of SM1 at rest as a neural biomarker of grip strength. Grip and its neural underpinnings decrease with advancing age and increasing NLR in mid to late life. Age-related decrements in grip strength and functional connectivity of brain regions involved in the generation of dynamic grip appear to be accelerated as a function of systemic physiological stress and inflammation, particularly in older men.

Changes in Functional Connectivity of Hippocampal Subregions in Patients with Obstructive Sleep Apnea after Six Months of Continuous Positive Airway Pressure Treatment

Previous studies have shown that the structural and functional impairments of hippocampal subregions in patients with obstructive sleep apnea (OSA) are related to cognitive impairment. Continuous positive airway pressure (CPAP) treatment can improve the clinical symptoms of OSA. Therefore, this study aimed to investigate functional connectivity (FC) changes in hippocampal subregions of patients with OSA after six months of CPAP treatment (post-CPAP) and its relationship with neurocognitive function. We collected and analyzed baseline (pre-CPAP) and post-CPAP data from 20 patients with OSA, including sleep monitoring, clinical evaluation, and resting-state functional magnetic resonance imaging. The results showed that compared with pre-CPAP OSA patients, the FC between the right anterior hippocampal gyrus and multiple brain regions, and between the left anterior hippocampal gyrus and posterior central gyrus were reduced in post-CPAP OSA patients. By contrast, the FC between the left middle hippocampus and the left precentral gyrus was increased. The changes in FC in these brain regions were closely related to cognitive dysfunction. Therefore, our findings suggest that CPAP treatment can effectively change the FC patterns of hippocampal subregions in patients with OSA, facilitating a better understanding of the neural mechanisms of cognitive function improvement, and emphasizing the importance of early diagnosis and timely treatment of OSA.

Neurofilament light-associated connectivity in young-adult Huntington's disease is related to neuronal genes

Upregulation of functional network connectivity in the presence of structural degeneration is seen in the premanifest stages of Huntington's disease (preHD) 10-15 years from clinical diagnosis. However, whether widespread network connectivity changes are seen in gene carriers much further from onset has yet to be explored. We characterized functional network connectivity throughout the brain and related it to a measure of disease pathology burden (CSF neurofilament light, NfL) and measures of structural connectivity in asymptomatic gene carriers, on average 24 years from onset. We related these measurements to estimates of cortical and subcortical gene expression. We found no overall differences in functional (or structural) connectivity anywhere in the brain comparing control and preHD participants. However, increased functional connectivity, particularly between posterior cortical areas, correlated with increasing CSF NfL level in preHD participants. Using the Allen Human Brain Atlas and expression-weighted cell-type enrichment analysis, we demonstrated that this functional connectivity upregulation occurred in cortical regions associated with regional expression of genes specific to neuronal cells. This relationship was validated using single-nucleus RNAseq data from post-mortem Huntington's disease and control brains showing enrichment of neuronal-specific genes that are differentially expressed in Huntington's disease. Functional brain networks in asymptomatic preHD gene carriers very far from disease onset show evidence of upregulated connectivity correlating with increased disease burden. These changes occur among brain areas that show regional expression of genes specific to neuronal GABAergic and glutamatergic cells.

Neural signatures of default mode network in major depression disorder after electroconvulsive therapy

Functional abnormalities of default mode network (DMN) have been well documented in major depressive disorder (MDD). However, the association of DMN functional reorganization with antidepressant treatment and gene expression is unclear. Moreover, whether the functional interactions of DMN could predict treatment efficacy is also unknown. Here, we investigated the link of treatment response with functional alterations of DMN and gene expression with a comparably large sample including 46 individuals with MDD before and after electroconvulsive therapy (ECT) and 46 age- and sex-matched healthy controls. Static and dynamic functional connectivity (dFC) analyses showed increased intrinsic/static but decreased dynamic functional couplings of inter- and intra-subsystems and between nodes of DMN. The changes of static functional connections of DMN were spatially correlated with brain gene expression profiles. Moreover, static and dFC of the DMN before treatment as features could predict depressive symptom improvement following ECT. Taken together, these results shed light on the underlying neural and genetic basis of antidepressant effect of ECT and the intrinsic functional connectivity of DMN have the potential to serve as prognostic biomarkers to guide accurate personalized treatment.

Transcriptional substrates underlying functional connectivity profiles of subregions within the human sensorimotor cortex

The human sensorimotor cortex has multiple subregions showing functional commonalities and differences, likely attributable to their connectivity profiles. However, the molecular substrates underlying such connectivity profiles are unclear. Here, transcriptome-neuroimaging spatial correlation analyses were performed between transcriptomic data from the Allen human brain atlas and resting-state functional connectivity (rsFC) of 24 fine-grained sensorimotor subregions from 793 healthy subjects. Results showed that rsFC of six sensorimotor subregions were associated with expression measures of six gene sets that were specifically expressed in brain tissue. These sensorimotor subregions could be classified into the polygenic- and oligogenic-modulated subregions, whose rsFC were related to gene sets diverging on their numbers (hundreds vs. dozens) and functional characteristics. First, the former were specifically expressed in multiple types of neurons and immune cells, yet the latter were not specifically expressed in any cortical cell types. Second, the former were preferentially expressed during the middle and late stages of cortical development, while the latter showed no preferential expression during any stages. Third, the former were prone to be enriched for general biological functions and pathways, but the latter for specialized biological functions and pathways. Fourth, the former were enriched for neuropsychiatric disorders, whereas this enrichment was absent for the latter. Finally, although the identified genes were commonly associated with sensorimotor behavioral processes, the polygenic-modulated subregions associated genes were additionally related to vision and dementia. These findings may advance our understanding of the functional homogeneity and heterogeneity of the human sensorimotor cortex from the perspective of underlying genetic architecture.

Genetic mechanisms underlying brain functional homotopy: a combined transcriptome and resting-state functional MRI study

Functional homotopy, the high degree of spontaneous activity synchrony and functional coactivation between geometrically corresponding interhemispheric regions, is a fundamental characteristic of the intrinsic functional architecture of the brain. However, little is known about the genetic mechanisms underlying functional homotopy. Resting-state functional magnetic resonance imaging data from a discovery dataset (656 healthy subjects) and 2 independent cross-race, cross-scanner validation datasets (103 and 329 healthy subjects) were used to calculate voxel-mirrored homotopic connectivity (VMHC) indexing brain functional homotopy. In combination with the Allen Human Brain Atlas, transcriptome-neuroimaging spatial correlation analysis was conducted to identify genes linked to VMHC. We found 1,001 genes whose expression measures were spatially associated with VMHC. Functional enrichment analyses demonstrated that these VMHC-related genes were enriched for biological functions including protein kinase activity, ion channel regulation, and synaptic function as well as many neuropsychiatric disorders. Concurrently, specific expression analyses showed that these genes were specifically expressed in the brain tissue, in neurons and immune cells, and during nearly all developmental periods. In addition, the VMHC-associated genes were linked to multiple behavioral domains, including vision, execution, and attention. Our findings suggest that interhemispheric communication and coordination involve a complex interaction of polygenes with a rich range of functional features.

Alzheimer’s Disease-Related Genes Identified by Linking Spatial Patterns of Pathology and Gene Expression

Background

Alzheimer’s Disease (AD) is an age-related neurodegenerative disease with a poorly understood etiology, shown to be partly genetic. Glucose hypometabolism, extracellular Amyloid-beta (Aβ) deposition, and intracellular Tau deposition are cardinal features of AD and display characteristic spatial patterns in the brain. We hypothesize that regional differences in underlying gene expression confer either resistance or susceptibility to AD pathogenic processes and are associated with these spatial patterns. Data-driven methods for the identification of genes involved in AD pathogenesis complement hypothesis-driven approaches that reflect current theories about the disease. Here we present a data driven method for the identification of genes involved in AD pathogenesis based on comparing spatial patterns of normal gene expression to Positron Emission Tomography (PET) images of glucose hypometabolism, Aβ deposition, and Tau deposition.

Methods

We performed correlations between the cerebral cortex microarray samples from the six cognitively normal (CN) post-mortem Allen Human Brain Atlas (AHBA) specimens and PET FDG-18, AV-45, and AV-1451 tracer images from AD and CN participants in the Alzheimer’s Disease and Neuroimaging Initiative (ADNI) database. Correlation coefficients for each gene by each ADNI subject were then entered into a partial least squares discriminant analysis (PLS-DA) to determine sets that best classified the AD and CN groups. Pathway analysis via BioPlanet 2019 was then used to infer the function of implicated genes.

Results

We identified distinct sets of genes strongly associated with each PET modality. Pathway analyses implicated novel genes involved in mitochondrial function, and Notch signaling, as well as genes previously associated with AD.

Conclusion

Using an unbiased approach, we derived sets of genes with expression patterns spatially associated with FDG hypometabolism, Aβ deposition, and Tau deposition in AD. This methodology may complement population-based approaches for identifying the genetic underpinnings of AD.

Abnormal Functional Connectivity of Hippocampal Subdivisions in Obstructive Sleep Apnea: A Resting-State Functional Magnetic Resonance Imaging Study

The hippocampus is involved in various cognitive function, including memory. Hippocampal structural and functional abnormalities have been observed in patients with obstructive sleep apnoea (OSA), but the functional connectivity (FC) patterns among hippocampal subdivisions in OSA patients remain unclear. The purpose of this study was to investigate the changes in FC between hippocampal subdivisions and their relationship with neurocognitive function in male patients with OSA. Resting-state fMRI were obtained from 46 male patients with untreated severe OSA and 46 male good sleepers. The hippocampus was divided into anterior, middle, and posterior parts, and the differences in FC between hippocampal subdivisions and other brain regions were determined. Correlation analysis was used to explore the relationships between abnormal FC of hippocampal subdivisions and clinical characteristics in patients with OSA. Our results revealed increased FC in the OSA group between the left anterior hippocampus and left middle temporal gyrus; between the left middle hippocampus and the left inferior frontal gyrus, right anterior central gyrus, and left anterior central gyrus; between the left posterior hippocampus and right middle frontal gyrus; between the right middle hippocampus and left inferior frontal gyrus; and between the right posterior hippocampus and left middle frontal gyrus. These FC abnormalities predominantly manifested in the sensorimotor network, fronto-parietal network, and semantic/default mode network, which are closely related to the neurocognitive impairment observed in OSA patients. This study advances our understanding of the potential pathophysiological mechanism of neurocognitive dysfunction in OSA.

Transcriptomic and cellular decoding of functional brain connectivity changes reveal regional brain vulnerability to pro- and anti-inflammatory therapies

BACKGROUND

Systemic inflammation induces acute changes in mood, motivation and cognition that closely resemble those observed in depressed individuals. However, the mechanistic pathways linking peripheral inflammation to depression-like psychopathology via intermediate effects on brain function remain incompletely understood.

METHODS

We combined data from 30 patients initiating interferon-α treatment for Hepatitis-C and 20 anti-tumour necrosis factor (TNF) therapy for inflammatory arthritis and used resting-state functional magnetic resonance imaging to investigate acute effects of each treatment on regional global brain connectivity (GBC). We leveraged transcriptomic data from the Allen Human Brain Atlas to uncover potential biological and cellular pathways underpinning regional vulnerability to GBC changes induced by each treatment.

RESULTS

Interferon-α and anti-TNF therapies both produced differential small-to-medium sized decreases in regional GBC. However, these were observed within distinct brain regions and the regional patterns of GBC changes induced by each treatment did not correlate suggesting independent underlying processes. Further, the spatial distribution of these differential GBC decreases could be captured by multivariate patterns of constitutive regional expression of genes respectively related to: i) neuroinflammation and glial cells; and ii) glutamatergic neurotransmission and neurons. The extent to which each participant expressed patterns of GBC changes aligning with these patterns of transcriptomic vulnerability also correlated with both acute treatment-induced changes in interleukin-6 (IL-6) and, for Interferon-α, longer-term treatment-associated changes in depressive symptoms.

CONCLUSIONS

Together, we present two transcriptomic models separately linking regional vulnerability to the acute effects of interferon-α and anti-TNF treatments on brain function to glial neuroinflammation and glutamatergic neurotransmission. These findings generate hypotheses about two potential brain mechanisms through which bidirectional changes in peripheral inflammation may contribute to the development/resolution of psychopathology.

Alternations and Applications of the Structural and Functional Connectome in Gliomas: A Mini-Review

In the central nervous system, gliomas are the most common, but complex primary tumors. Genome-based molecular and clinical studies have revealed different classifications and subtypes of gliomas. Neuroradiological approaches have non-invasively provided a macroscopic view for surgical resection and therapeutic effects. The connectome is a structural map of a physical object, the brain, which raises issues of spatial scale and definition, and it is calculated through diffusion magnetic resonance imaging (MRI) and functional MRI. In this study, we reviewed the basic principles and attributes of the structural and functional connectome, followed by the alternations of connectomes and their influences on glioma. To extend the applications of connectome, we demonstrated that a series of multi-center projects still need to be conducted to systemically investigate the connectome and the structural-functional coupling of glioma. Additionally, the brain-computer interface based on accurate connectome could provide more precise structural and functional data, which are significant for surgery and postoperative recovery. Besides, integrating the data from different sources, including connectome and other omics information, and their processing with artificial intelligence, together with validated biological and clinical findings will be significant for the development of a personalized surgical strategy.

Integrative Functional, Molecular, and Transcriptomic Analyses of Altered Intrinsic Timescale Gradient in Depression

The pathophysiology and pharmacology of depression are hypothesized to be related to the imbalance of excitation–inhibition that gives rise to hierarchical dynamics (or intrinsic timescale gradient), further supporting a hierarchy of cortical functions. On this assumption, intrinsic timescale gradient is theoretically altered in depression. However, it remains unknown. We investigated altered intrinsic timescale gradient recently developed to measure hierarchical brain dynamics gradient and its underlying molecular architecture and brain-wide gene expression in depression. We first presented replicable intrinsic timescale gradient in two independent Chinese Han datasets and then investigated altered intrinsic timescale gradient and its possible underlying molecular and transcriptional bases in patients with depression. As a result, patients with depression showed stage-specifically shorter timescales compared with healthy controls according to illness duration. The shorter timescales were spatially correlated with monoamine receptor/transporter densities, suggesting the underlying molecular basis of timescale aberrance and providing clues to treatment. In addition, we identified that timescale aberrance-related genes ontologically enriched for synapse-related and neurotransmitter (receptor) terms, elaborating the underlying transcriptional basis of timescale aberrance. These findings revealed atypical timescale gradient in depression and built a link between neuroimaging, transcriptome, and neurotransmitter information, facilitating an integrative understanding of depression.

Imaging transcriptomics: Convergent cellular, transcriptomic, and molecular neuroimaging signatures in the healthy adult human brain

The integration of transcriptomic and neuroimaging data, "imaging transcriptomics," has recently emerged to generate hypotheses about potential biological pathways underlying regional variability in neuroimaging features. However, the validity of this approach is yet to be examined in depth. Here, we sought to bridge this gap by performing transcriptomic decoding of the regional distribution of well-known molecular markers spanning different elements of the biology of the healthy human brain. Imaging transcriptomics identifies biological and cell pathways that are consistent with the known biology of a wide range of molecular neuroimaging markers. The extent to which it can capture patterns of gene expression that align well with elements of the biology of the neuroinflammatory axis, at least in healthy controls without a proinflammatory challenge, is inconclusive. Imaging transcriptomics might constitute an interesting approach to improve our understanding of the biological pathways underlying regional variability in a wide range of neuroimaging phenotypes.

Dissect Relationships Between Gene Co-expression and Functional Connectivity in Human Brain

Although recent evidence indicates an association between gene co-expression and functional connectivity in human brain, specific association patterns remain largely unknown. Here, using neuroimaging-based functional connectivity data of living brains and brain-wide gene expression data of postmortem brains, we performed comprehensive analyses to dissect relationships between gene co-expression and functional connectivity. We identified 125 connectivity-related genes (20 novel genes) enriched for dendrite extension, signaling pathway and schizophrenia, and 179 gene-related functional connections mainly connecting intra-network regions, especially homologous cortical regions. In addition, 51 genes were associated with connectivity in all brain functional networks and enriched for action potential and schizophrenia; in contrast, 51 genes showed network-specific modulatory effects and enriched for ion transportation. These results indicate that functional connectivity is unequally affected by gene expression, and connectivity-related genes with different biological functions are involved in connectivity modulation of different networks.

Genetic Architecture Underlying Differential Resting-state Functional Connectivity of Subregions Within the Human Visual Cortex

The human visual cortex is a heterogeneous entity that has multiple subregions showing substantial variability in their functions and connections. We aimed to identify genes associated with resting-state functional connectivity (rsFC) of visual subregions using transcriptome-neuroimaging spatial correlations in discovery and validation datasets. Results showed that rsFC of eight visual subregions were associated with expression measures of eight gene sets, which were specifically expressed in brain tissue and showed the strongest correlations with visual behavioral processes. Moreover, there was a significant divergence in these gene sets and their functional features between medial and lateral visual subregions. Relative to those associated with lateral subregions, more genes associated with medial subregions were found to be enriched for neuropsychiatric diseases and more diverse biological functions and pathways, and to be specifically expressed in multiple types of neurons and immune cells and during the middle and late stages of cortical development. In addition to shared behavioral processes, lateral subregion associated genes were uniquely correlated with high-order cognition. These findings of commonalities and differences in the identified rsFC-related genes and their functional features across visual subregions may improve our understanding of the functional heterogeneity of the visual cortex from the perspective of underlying genetic architecture.

The Vital Role of Central Executive Network in Brain Age: Evidence From Machine Learning and Transcriptional Signatures

Recent studies combining neuroimaging with machine learning methods successfully infer an individual’s brain age, and its discrepancy with the chronological age is used to identify age-related diseases. However, which brain networks play decisive roles in brain age prediction and the underlying biological basis of brain age remain unknown. To answer these questions, we estimated an individual’s brain age in the Southwest University Adult Lifespan Dataset (N = 492) from the gray matter volumes (GMV) derived from T1-weighted MRI scans by means of Gaussian process regression. Computational lesion analysis was performed to determine the importance of each brain network in brain age prediction. Then, we identified brain age-related genes by using prior brain-wide gene expression data, followed by gene enrichment analysis using Metascape. As a result, the prediction model successfully inferred an individual’s brain age and the computational lesion prediction results identified the central executive network as a vital network in brain age prediction (Steiger’s Z = 2.114, p = 0.035). In addition, the brain age-related genes were enriched in Gene Ontology (GO) processes/Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways grouped into numbers of clusters, such as regulation of iron transmembrane transport, synaptic signaling, synapse organization, retrograde endocannabinoid signaling (e.g., dopaminergic synapse), behavior (e.g., memory and associative learning), neurotransmitter secretion, and dendrite development. In all, these results reveal that the GMV of the central executive network played a vital role in predicting brain age and bridged the gap between transcriptome and neuroimaging promoting an integrative understanding of the pathophysiology of brain age.

Correlation between cortical gene expression and resting-state functional network centrality in healthy young adults

Resting‐state functional connectivity in the human brain is heritable, and previous studies have investigated the genetic basis underlying functional connectivity. However, at present, the molecular mechanisms associated with functional network centrality are still largely unknown. In this study, functional networks were constructed, and the graph‐theory method was employed to calculate network centrality in 100 healthy young adults from the Human Connectome Project. Specifically, functional connectivity strength (FCS), also known as the “degree centrality” of weighted networks, is calculated to measure functional network centrality. A multivariate technique of partial least squares regression (PLSR) was then conducted to identify genes whose spatial expression profiles best predicted the FCS distribution. We found that FCS spatial distribution was significantly positively correlated with the expression of genes defined by the first PLSR component. The FCS‐related genes we identified were significantly enriched for ion channels, axon guidance, and synaptic transmission. Moreover, FCS‐related genes were preferentially expressed in cortical neurons and young adulthood and were enriched in numerous neurodegenerative and neuropsychiatric disorders. Furthermore, a series of validation and robustness analyses demonstrated the reliability of the results. Overall, our results suggest that the spatial distribution of FCS is modulated by the expression of a set of genes associated with ion channels, axon guidance, and synaptic transmission.