Canonical genetic signatures of the adult human brain

MET receptor tyrosine kinase promotes the generation of functional synapses in adult cortical circuits

JOURNAL/nrgr/04.03/01300535-202505000-00026/figure1/v/2024-07-28T173839Z/r/image-tiff Loss of synapse and functional connectivity in brain circuits is associated with aging and neurodegeneration, however, few molecular mechanisms are known to intrinsically promote synaptogenesis or enhance synapse function. We have previously shown that MET receptor tyrosine kinase in the developing cortical circuits promotes dendritic growth and dendritic spine morphogenesis. To investigate whether enhancing MET in adult cortex has synapse regenerating potential, we created a knockin mouse line, in which the human MET gene expression and signaling can be turned on in adult (10-12 months) cortical neurons through doxycycline-containing chow. We found that similar to the developing brain, turning on MET signaling in the adult cortex activates small GTPases and increases spine density in prefrontal projection neurons. These findings are further corroborated by increased synaptic activity and transient generation of immature silent synapses. Prolonged MET signaling resulted in an increased α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid/N-methyl-D-aspartate (AMPA/NMDA) receptor current ratio, indicative of enhanced synaptic function and connectivity. Our data reveal that enhancing MET signaling could be an interventional approach to promote synaptogenesis and preserve functional connectivity in the adult brain. These findings may have implications for regenerative therapy in aging and neurodegeneration conditions.

Molecular mechanisms and behavioral relevance underlying neural correlates of childhood neglect

BACKGROUND

Childhood neglect is associated with brain changes, yet the molecular mechanisms and behavioral relevance underlying such associations remain elusive.

METHODS

We calculated fractional amplitude of low-frequency fluctuations (fALFF) using resting-state functional MRI and tested their correlation with childhood neglect across a large sample of 510 healthy young adults. Then, we investigated the spatial relationships of the identified neural correlates of childhood neglect with gene expression, neurotransmitter, and behavioral domain atlases.

RESULTS

We found that more severe childhood neglect was correlated with higher fALFF in the bilateral anterior cingulate cortex. Remarkably, the identified neural correlates of childhood neglect were spatially correlated with expression of gene categories primarily involving neuron, synapse, ion channel, cognitive and perceptual processes, and physiological response and regulation. Concurrently, there were significant associations between the neural correlates and specific neurotransmitter systems including serotonin and GABA. Finally, neural correlates of childhood neglect were associated with diverse behavioral domains implicating mental disorders, emotion, cognition, and sensory perception.

LIMITATIONS

The cross-sectional study design cannot unequivocally establish causality.

CONCLUSIONS

Our findings may not only add to the current knowledge regarding the relationship between childhood neglect and mental health, but also have clinical implications for developing preventive strategies for individuals exposed to childhood neglect who are at risk for mental disorders.

Neuro-immune communication at the core of craving-associated brain structural network reconfiguration in methamphetamine users

Methamphetamine (MA) use disorder is a chronic neurotoxic brain disease characterized by a high risk of relapse driven by intense cravings. However, the neurobiological signatures of cravings remain unclear, limiting the effectiveness of various treatment methods. Diffusion MRI (dMRI) scans from 62 MA users and 57 healthy controls (HC) were used in this study. The MA users were longitudinally followed up during their period of long-term abstinence (duration of long-term abstinence: 347.52±99.25 days). We systematically quantified the control ability of each brain region for craving-associated state transitions using network control theory from a causal perspective. Craving-associated structural alterations (CSA) were investigated through multivariate group comparisons and biological relevance analysis. The neural mechanisms underlying CSA were elucidated using transcriptomic and neurochemical analyses. We observed that long-term abstinence-induced structural alterations significantly influenced the state transition energy involved in the cognitive control response to external information, which correlated with changes in craving scores (r ∼ 0.35, P <0.01). Our causal network analysis further supported the crucial role of the prefrontal cortex (PFC) in craving mechanisms. Notably, while the PFC is central to the craving, the CSAs were distributed widely across multiple brain regions (PFDR<0.05), with strong alterations in somatomotor regions (PFDR<0.05) and moderate alterations in high-level association networks (PFDR<0.05). Additionally, transcriptomic, chemical compounds, cell-type analyses, and molecular imaging collectively highlight the influence of neuro-immune communication on human craving modulation. Our results offer an integrative, multi-scale perspective on unraveling the neural underpinnings of craving and suggest that neuro-immune signaling may be a promising target for future human addiction therapeutics.

Multiscale and multimodal signatures of structure-function coupling variability across the human neocortex

The relationship between the brain's structural wiring and its dynamic activity is thought to vary regionally, implying that the mechanisms underlying structure-function coupling may differ depending on a region's position within the brain's hierarchy. To better bridge the gap between structure and function, it is crucial to identify the factors shaping this regionality, not only in terms of how static functional connectivity aligns with structure, but also regarding the time-domain variability of this interplay. Here we map structure - function coupling and its time-domain variability and relate them to the heterogeneity of the cortex. We show that these two properties split the cortical landscape into two districts anchored to the opposite ends of the brain's hierarchy. By looking at statistical relationships with layer-specific gene transcription, T1w/T2 w ratio, and synaptic density, we show that macro-scale structure-function coupling may be rooted in the brain's microstructure and meso‑scale laminar specialization. Finally, we demonstrate that a lower and more variable alignment of function and structure may bestow the emergence of unique functional dynamics.

Genomic structural equation modeling reveals latent phenotypes in the human cortex with distinct genetic architecture

Genetic contributions to human cortical structure manifest pervasive pleiotropy. This pleiotropy may be harnessed to identify unique genetically-informed parcellations of the cortex that are neurobiologically distinct from functional, cytoarchitectural, or other cortical parcellation schemes. We investigated genetic pleiotropy by applying genomic structural equation modeling (SEM) to map the genetic architecture of cortical surface area (SA) and cortical thickness (CT) for 34 brain regions recently reported in the ENIGMA cortical GWAS. Genomic SEM uses the empirical genetic covariance estimated from GWAS summary statistics with LD score regression (LDSC) to discover factors underlying genetic covariance, which we are denoting genetically informed brain networks (GIBNs). Genomic SEM can fit a multivariate GWAS from summary statistics for each of the GIBNs, which can subsequently be used for LD score regression (LDSC). We found the best-fitting model of cortical SA identified 6 GIBNs and CT identified 4 GIBNs, although sensitivity analyses indicated that other structures were plausible. The multivariate GWASs of the GIBNs identified 74 genome-wide significant (GWS) loci (p < 5 × 10-8), including many previously implicated in neuroimaging phenotypes, behavioral traits, and psychiatric conditions. LDSC of GIBN GWASs found that SA-derived GIBNs had a positive genetic correlation with bipolar disorder (BPD), and cannabis use disorder, indicating genetic predisposition to a larger SA in the specific GIBN is associated with greater genetic risk of these disorders. A negative genetic correlation was observed between attention deficit hyperactivity disorder (ADHD) and major depressive disorder (MDD). CT GIBNs displayed a negative genetic correlation with alcohol dependence. Even though we observed model instability in our application of genomic SEM to high-dimensional data, jointly modeling the genetic architecture of complex traits and investigating multivariate genetic links across neuroimaging phenotypes offers new insights into the genetics of cortical structure and relationships to psychopathology.

Molecular mechanisms underlying the neural correlates of working memory

BACKGROUND

Working memory (WM), a core component of executive functions, relies on a dedicated brain system that maintains and stores information in the short term. While extensive neuroimaging research has identified a distributed set of neural substrates relevant to WM, their underlying molecular mechanisms remain enigmatic. This study investigated the neural correlates of WM as well as their underlying molecular mechanisms.

RESULTS

Our voxel-wise analyses of resting-state functional MRI data from 502 healthy young adults showed that better WM performance (higher accuracy and shorter reaction time of the 3-back task) was associated with lower functional connectivity density (FCD) in the left inferior temporal gyrus and higher FCD in the left anterior cingulate cortex. A combination of transcriptome-neuroimaging spatial correlation and the ensemble-based gene category enrichment analysis revealed that the identified neural correlates of WM were associated with expression of diverse gene categories involving important cortical components and their biological processes as well as sodium channels. Cross-region spatial correlation analyses demonstrated significant associations between the neural correlates of WM and a range of neurotransmitters including dopamine, glutamate, serotonin, and acetylcholine.

CONCLUSIONS

These findings may help to shed light on the molecular mechanisms underlying the neural correlates of WM.

Genetic and molecular correlates of cortical thickness alterations in adults with obsessive-compulsive disorder: a transcription-neuroimaging association analysis

BACKGROUND

Although numerous neuroimaging studies have depicted neural alterations in individuals with obsessive-compulsive disorder (OCD), a psychiatric disorder characterized by intrusive cognitions and repetitive behaviors, the molecular mechanisms connecting brain structural changes and gene expression remain poorly understood.

METHODS

This study combined the Allen Human Brain Atlas dataset with neuroimaging data from the Meta-Analysis (ENIGMA) consortium and independent cohorts. Later, partial least squares regression and enrichment analysis were performed to probe the correlation between transcription and cortical thickness variation among adults with OCD.

RESULTS

The cortical map of case-control differences in cortical thickness was spatially correlated with cortical expression of a weighted combination of genes enriched for neurobiologically relevant ontology terms preferentially expressed across different cell types and cortical layers. These genes were specifically expressed in brain tissue, spanning all cortical developmental stages. Protein-protein interaction analysis revealed that these genes coded a network of proteins encompassing various highly interactive hubs.

CONCLUSIONS

The study findings bridge the gap between neural structure and transcriptome data in OCD, fostering an integrative understanding of the potential biological mechanisms.

VEGF-A-mediated venous endothelial cell proliferation results in neoangiogenesis during neuroinflammation

Newly formed leaky vessels and blood-brain barrier (BBB) damage are present in demyelinating acute and chronic lesions in multiple sclerosis (MS) and experimental autoimmune encephalomyelitis (EAE). However, the endothelial cell subtypes and signaling pathways contributing to these leaky neovessels are unclear. Here, using single-cell transcriptional profiling and in vivo validation studies, we show that venous endothelial cells express neoangiogenesis gene signatures and show increased proliferation resulting in enlarged veins and higher venous coverage in acute and chronic EAE lesions in female adult mice. These changes correlate with the upregulation of vascular endothelial growth factor A (VEGF-A) signaling. We also confirmed increased expression of neoangiogenic markers in acute and chronic human MS lesions. Treatment with a VEGF-A blocking antibody diminishes the neoangiogenic transcriptomic signatures and vascular proliferation in female adult mice with EAE, but it does not restore BBB function or ameliorate EAE pathology. Our data demonstrate that venous endothelial cells contribute to neoangiogenesis in demyelinating neuroinflammatory conditions.

Local molecular and connectomic contributions of tau-related neurodegeneration

Neurodegeneration in Alzheimer's disease (AD) is known to be mostly driven by tau neurofibrillary tangles. However, both tau and neurodegeneration exhibit variability in their distribution across the brain and among individuals, and the relationship between tau and neurodegeneration might be influenced by several factors. I aimed to map local molecular and connectivity characteristics that affect the association between tau pathology and neurodegeneration. The current study was conducted on the cross-sectional tau-PET and longitudinal T1-weighted MRI scan data of 186 participants from the ADNI dataset including 71 cognitively unimpaired (CU) and 115 mild cognitive impairment (MCI) individuals. Furthermore, the normative molecular profile of a region was defined using neurotransmitter receptor densities, gene expression, T1w/T2w ratio (myelination), FDG-PET (glycolytic index, glucose metabolism, and oxygen metabolism), and synaptic density. I found that the excitatory-inhibitory (E:I) ratio, myelination, synaptic density, glycolytic index, and functional connectivity are linked with deviation in the relationship between tau and neurodegeneration. Furthermore, there was spatial similarity between tau pathology and glycolytic index, synaptic density, and functional connectivity across brain regions. The current study demonstrates that the regional susceptibility to tau-related neurodegeneration is associated with specific molecular and connectomic characteristics of the affected neural systems. I found that the molecular and connectivity architecture of the human brain is linked to the different effects of tau pathology on downstream neurodegeneration.

Abnormal structural‒functional coupling patterning in progressive supranuclear palsy is associated with diverse gradients and histological features

The anatomy of the brain supports inherent processes, fostering mental abilities and eventually facilitating adaptive behavior. Recent studies have shown that progressive supranuclear palsy (PSP) is accompanied by alterations in functional and structural networks. However, how the structure and function of PSP coordinates change is not clear, and the relationships between structural‒functional coupling (SFC) and the gradient of hierarchical structure and cellular histology remain largely unknown. Here, we use neuroimaging data from two independent cohorts and a public histological dataset to investigate the relationships among the cellular histology, hierarchical structure, and SFC of PSP patients. We find that the SFC of the entire cortex in PSP is severely disrupted, with higher coupling in the visual network (VN). Moreover, coupling differences in PSP follow a macroscopic organizational principle from unimodal to transmodal gradients. Finally, we elucidate greater laminar differentiation in VN regions sensitive to SFC changes in PSP, which is related mainly to the higher cellular density and smaller size of the internal-granular layer. In conclusion, our findings provide an interpretable framework for understanding SFC changes in PSP and provide new insights into the consistency of structural and functional changes in PSP regarding hierarchical structure and cellular histology.

Associations between the multitrajectory neuroplasticity of neuronavigated rTMS-mediated angular gyrus networks and brain gene expression in AD spectrum patients with sleep disorders

INTRODUCTION

The multifactorial influence of repetitive transcranial magnetic stimulation (rTMS) on neuroplasticity in neural networks is associated with improvements in cognitive dysfunction and sleep disorders. The mechanisms of rTMS and the transcriptional-neuronal correlation in Alzheimer's disease (AD) patients with sleep disorders have not been fully elucidated.

METHODS

Forty-six elderly participants with cognitive impairment (23 patients with low sleep quality and 23 patients with high sleep quality) underwent 4-week periods of neuronavigated rTMS of the angular gyrus and neuroimaging tests, and gene expression data for six post mortem brains were collected from another database. Transcription-neuroimaging association analysis was used to evaluate the effects on cognitive dysfunction and the underlying biological mechanisms involved.

RESULTS

Distinct variable neuroplasticity in the anterior and posterior angular gyrus networks was detected in the low sleep quality group. These interactions were associated with multiple gene pathways, and the comprehensive effects were associated with improvements in episodic memory.

DISCUSSION

Multitrajectory neuroplasticity is associated with complex biological mechanisms in AD-spectrum patients with sleep disorders.

HIGHLIGHTS

This was the first transcription-neuroimaging study to demonstrate that multitrajectory neuroplasticity in neural circuits was induced via neuronavigated rTMS, which was associated with complex gene expression in AD-spectrum patients with sleep disorders. The interactions between sleep quality and neuronavigated rTMS were coupled with multiple gene pathways and improvements in episodic memory. The present strategy for integrating neuroimaging, rTMS intervention, and genetic data provide a new approach to comprehending the biological mechanisms involved in AD.

Structural lesions and transcriptomic specializations shape gradient perturbations in Wilson disease

Functional dysregulations in multiple regions are caused by excessive copper deposition in the brain in Wilson disease (WD) patients. The genetic mechanism of WD is thought to involve the abnormal expression of ATP7B in the liver, whereas the biological and molecular processes involved in functional dysregulation within the brain remain unexplored. The objective of this study was to unravel the underpinnings of functional gradient perturbations underlying structural lesions and transcriptomic specializations in WD. In this study, we included 105 WD patients and 93 healthy controls who underwent structural and functional MRI assessments. We used the diffusion mapping embedding model to derive the functional connectome gradient and further employed gray matter volume to uncover structure-function decoupling for WD. Then, we used Neurosynth, clinical data, and whole-brain gene expression data to examine the meta-analytic cognitive function, clinical phenotypes, and transcriptomic specializations related to WD gradient alterations. Compared with controls, WD patients exhibited global topographic changes in the principal pramary-to-transmodal gradient. Meta-analytic terms and clinical characteristics were correlated with these gradient alterations in motor-related processing, higher-order cognition, neurological symptoms, and age. Spatial correlations revealed structure-function decoupling in multiple networks, especially in subcortical and visual networks. Within the cortex, the spatial association between gradient alterations and gene expression profiles has revealed transcriptomic specilizations in WD that display properties indicative of ion homeostasis, neural development, and motor control. Furthermore, for the first time, we characterized the role of the ATP7B gene in impacting subcortical function. The transcriptomic specializations of WD were also associated with other neurological and psychiatric disorders. Finally, we revealed that structural lesions and gradient perturbations may share similar transcriptomic specializations in WD. In conclusion, these findings bridged functional gradient perturbations to structural lesions and gene expression profiles in WD patients, possibly promoting our understanding of the neurobiological mechanisms underlying the emergence of complex neurological and psychiatric phenotypes.

The brain, rapid eye movement sleep, and major depressive disorder: A multimodal neuroimaging study

BACKGROUND

Evidence has established the prominent involvement of rapid eye movement (REM) sleep disturbance in major depressive disorder (MDD). However, the neural correlates of REM sleep in MDD and their clinical significance are less clear.

METHODS

Cross-sectional and longitudinal polysomnography and resting-state functional MRI data were collected from 131 MDD patients and 71 healthy controls to measure REM sleep and voxel-mirrored homotopic connectivity (VMHC). Correlation and mediation analyses were performed to examine the associations between REM sleep, VMHC, and clinical variables. Moreover, we conducted spatial correlations between the neural correlates of REM sleep and a multimodal collection of reference brain maps to facilitate genetic, structural and functional annotations.

RESULTS

MDD patients exhibited REM sleep abnormalities manifesting as higher REM sleep latency and lower REM sleep duration, which were correlated with decreased VMHC of the precentral gyrus and inferior parietal lobe and mediated their associations with more severe anxiety symptoms. Longitudinal data showed that VMHC increase of the inferior parietal lobe was related to improvement of depression symptoms in MDD patients. Spatial correlation analyses revealed that the neural correlates of REM sleep in MDD were linked to gene categories primarily involving cellular metabolic process, signal pathway, and ion channel activity as well as linked to cortical microstructure, metabolism, electrophysiology, and cannabinoid receptor.

CONCLUSION

These findings may add important context to the growing literature on the complex interplay between sleep and MDD, and more broadly may inform future treatment for depression via regulating sleep.

Beyond the usual suspects: multi-factorial computational models in the search for neurodegenerative disease mechanisms

From Alzheimer's disease to amyotrophic lateral sclerosis, the molecular cascades underlying neurodegenerative disorders remain poorly understood. The clinical view of neurodegeneration is confounded by symptomatic heterogeneity and mixed pathology in almost every patient. While the underlying physiological alterations originate, proliferate, and propagate potentially decades before symptomatic onset, the complexity and inaccessibility of the living brain limit direct observation over a patient's lifespan. Consequently, there is a critical need for robust computational methods to support the search for causal mechanisms of neurodegeneration by distinguishing pathogenic processes from consequential alterations, and inter-individual variability from intra-individual progression. Recently, promising advances have been made by data-driven spatiotemporal modeling of the brain, based on in vivo neuroimaging and biospecimen markers. These methods include disease progression models comparing the temporal evolution of various biomarkers, causal models linking interacting biological processes, network propagation models reproducing the spatial spreading of pathology, and biophysical models spanning cellular- to network-scale phenomena. In this review, we discuss various computational approaches for integrating cross-sectional, longitudinal, and multi-modal data, primarily from large observational neuroimaging studies, to understand (i) the temporal ordering of physiological alterations, i(i) their spatial relationships to the brain's molecular and cellular architecture, (iii) mechanistic interactions between biological processes, and (iv) the macroscopic effects of microscopic factors. We consider the extents to which computational models can evaluate mechanistic hypotheses, explore applications such as improving treatment selection, and discuss how model-informed insights can lay the groundwork for a pathobiological redefinition of neurodegenerative disorders.

CYP1B1-RMDN2 Alzheimer's disease endophenotype locus identified for cerebral tau PET

Determining the genetic architecture of Alzheimer's disease pathologies can enhance mechanistic understanding and inform precision medicine strategies. Here, we perform a genome-wide association study of cortical tau quantified by positron emission tomography in 3046 participants from 12 independent studies. The CYP1B1-RMDN2 locus is associated with tau deposition. The most significant signal is at rs2113389, explaining 4.3% of the variation in cortical tau, while APOE4 rs429358 accounts for 3.6%. rs2113389 is associated with higher tau and faster cognitive decline. Additive effects, but no interactions, are observed between rs2113389 and diagnosis, APOE4, and amyloid beta positivity. CYP1B1 expression is upregulated in AD. rs2113389 is associated with higher CYP1B1 expression and methylation levels. Mouse model studies provide additional functional evidence for a relationship between CYP1B1 and tau deposition but not amyloid beta. These results provide insight into the genetic basis of cerebral tau deposition and support novel pathways for therapeutic development in AD.

The macroscale routing mechanism of structural brain connectivity related to body mass index

Understanding the brain's mechanisms in individuals with obesity is important for managing body weight. Prior neuroimaging studies extensively investigated alterations in brain structure and function related to body mass index (BMI). However, how the network communication among the large-scale brain networks differs across BMI is underinvestigated. This study used diffusion magnetic resonance imaging of 290 young adults to identify links between BMI and brain network mechanisms. Navigation efficiency, a measure of network routing, was calculated from the structural connectivity computed using diffusion tractography. The sensory and frontoparietal networks indicated positive associations between navigation efficiency and BMI. The neurotransmitter association analysis identified that serotonergic and dopaminergic receptors, as well as opioid and norepinephrine systems, were related to BMI-related alterations in navigation efficiency. The transcriptomic analysis found that genes associated with network routing across BMI overlapped with genes enriched in excitatory and inhibitory neurons, specifically, gene enrichments related to synaptic transmission and neuron projection. Our findings suggest a valuable insight into understanding BMI-related alterations in brain network routing mechanisms and the potential underlying cellular biology, which might be used as a foundation for BMI-based weight management.

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

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

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.

Accelerated Cortical Thinning in Schizophrenia Is Associated With Rare and Common Predisposing Variation to Schizophrenia and Neurodevelopmental Disorders

BACKGROUND

Schizophrenia is a highly heritable disorder characterized by increased cortical thinning throughout the life span. Studies have reported a shared genetic basis between schizophrenia and cortical thickness. However, no genes whose expression is related to abnormal cortical thinning in schizophrenia have been identified.

METHODS

We conducted linear mixed models to estimate the rates of accelerated cortical thinning across 68 regions from the Desikan-Killiany atlas in individuals with schizophrenia compared with healthy control participants from a large longitudinal sample (ncases = 169 and ncontrols = 298, ages 16-70 years). We studied the correlation between gene expression data from the Allen Human Brain Atlas and accelerated thinning estimates across cortical regions. Finally, we explored the functional and genetic underpinnings of the genes that contribute most to accelerated thinning.

RESULTS

We found a global pattern of accelerated cortical thinning in individuals with schizophrenia compared with healthy control participants. Genes underexpressed in cortical regions that exhibit this accelerated thinning were downregulated in several psychiatric disorders and were enriched for both common and rare disrupting variation for schizophrenia and neurodevelopmental disorders. In contrast, none of these enrichments were observed for baseline cross-sectional cortical thickness differences.

CONCLUSIONS

Our findings suggest that accelerated cortical thinning, rather than cortical thickness alone, serves as an informative phenotype for neurodevelopmental disruptions in schizophrenia. We highlight the genetic and transcriptomic correlates of this accelerated cortical thinning, emphasizing the need for future longitudinal studies to elucidate the role of genetic variation and the temporal-spatial dynamics of gene expression in brain development and aging in schizophrenia.

Universality and language specificity of brain reading networks: A developmental perspective

The development of reading networks across different languages and cultures provides an important window to address gene-culture interactions in brain functionality development. Previous meta-analyses have explored the neural correlates of reading in different languages with diverse orthographic transparencies. However, it remains unknown whether the neural topographic relationship of different languages varies when taking development into account. To address this issue, we conducted meta-analyses of neuroimaging studies with approaches of activation likelihood estimation and seed-based effect size mapping and focused on two highly contrasting languages, Chinese and English. The meta-analyses covered 61 studies of Chinese reading and 64 studies of English reading by native speakers. The brain reading networks of child and adult readers were further separately analyzed and compared to explore the developmental effects. The results revealed that the commonalities and differences in reading networks for Chinese and English were inconsistent between children and adults. In addition, the reading networks converged with development, and the effects of writing systems on brain function organizations were more salient in the initial stages of reading. An interesting finding was that the left inferior parietal lobule demonstrated increased effect sizes in adults compared with children in both Chinese and English reading, indicating a common developmental feature of reading mechanisms across the two languages. These findings provide new insights into the functional evolution and cultural modulation of brain reading networks. RESEARCH HIGHLIGHTS: Meta-analyses with approaches of activation likelihood estimation and seed-based effect size mapping were performed to evaluate the developmental characteristics of brain reading networks. The engagement of universal and language-specific reading networks was different between children and adults, and the reading networks converged with increased reading experience. Overall the middle/inferior occipital and inferior/middle frontal gyrus were specific to Chinese and the middle temporal, right inferior frontal gyrus were specific to English. The left inferior parietal lobule was engaged more in adults than children in Chinese and English reading, demonstrating a common developmental feature of reading mechanisms.

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.

Unique Functional Neuroimaging Signatures of Genetic Versus Clinical High Risk for Psychosis

BACKGROUND

22q11.2 deletion syndrome (22qDel) is a copy number variant that is associated with psychosis and other neurodevelopmental disorders. Adolescents who are at clinical high risk for psychosis (CHR) are identified based on the presence of subthreshold psychosis symptoms. Whether common neural substrates underlie these distinct high-risk populations is unknown. We compared functional brain measures in 22qDel and CHR cohorts and mapped the results to biological pathways.

METHODS

We analyzed 2 large multisite cohorts with resting-state functional magnetic resonance imaging data: 1) a 22qDel cohort (n = 164, 47% female) and typically developing (TD) control participants (n = 134, 56% female); and 2) a cohort of CHR individuals (n = 240, 41% female) and TD control participants (n = 149, 46% female) from the NAPLS-2 (North American Prodrome Longitudinal Study-2). We computed global brain connectivity (GBC), local connectivity (LC), and brain signal variability (BSV) across cortical regions and tested case-control differences for 22qDel and CHR separately. Group difference maps were related to published brain maps using autocorrelation-preserving permutation.

RESULTS

BSV, LC, and GBC were significantly disrupted in individuals with 22qDel compared with TD control participants (false discovery rate-corrected q < .05). Spatial maps of BSV and LC differences were highly correlated with each other, unlike GBC. In the CHR group, only LC was significantly altered versus the control group, with a different spatial pattern than the 22qDel group. Group differences mapped onto biological gradients, with 22qDel effects being strongest in regions with high predicted blood flow and metabolism.

CONCLUSIONS

22qDel carriers and CHR individuals exhibited different effects on functional magnetic resonance imaging temporal variability and multiscale functional connectivity. In 22qDel carriers, strong and convergent disruptions in BSV and LC that were not seen in CHR individuals suggest distinct functional brain alterations.

Voxel-based texture similarity networks reveal individual variability and correlate with biological ontologies

The human brain is organized as a complex, hierarchical network. However, the structural covariance patterns among brain regions and the underlying biological substrates of such covariance networks remain to be clarified. The present study proposed a novel individualized structural covariance network termed voxel-based texture similarity networks (vTSNs) based on 76 refined voxel-based textural features derived from structural magnetic resonance images. Validated in three independent longitudinal healthy cohorts (40, 23, and 60 healthy participants, respectively) with two common brain atlases, we found that the vTSN could robustly resolve inter-subject variability with high test-retest reliability. In contrast to the regional-based texture similarity networks (rTSNs) that calculate radiomic features based on region-of-interest information, vTSNs had higher inter- and intra-subject variability ratios and test-retest reliability in connectivity strength and network topological properties. Moreover, the Spearman correlation indicated a stronger association of the gene expression similarity network (GESN) with vTSNs than with rTSNs (vTSN: r = 0.600, rTSN: r = 0.433, z = 39.784, P < 0.001). Hierarchical clustering identified 3 vTSN subnets with differential association patterns with 13 coexpression modules, 16 neurotransmitters, 7 electrophysiology, 4 metabolism, and 2 large-scale structural and 4 functional organization maps. Moreover, these subnets had unique biological hierarchical organization from the subcortex-limbic system to the ventral neocortex and then to the dorsal neocortex. Based on 424 unrelated, qualified healthy subjects from the Human Connectome Project, we found that vTSNs could sensitively represent sex differences, especially for connections in the subcortex-limbic system and between the subcortex-limbic system and the ventral neocortex. Moreover, a multivariate variance component model revealed that vTSNs could explain a significant proportion of inter-subject behavioral variance in cognition (80.0 %) and motor functions (63.4 %). Finally, using 494 healthy adults (aged 19-80 years old) from the Southwest University Adult Lifespan Dataset, the Spearman correlation identified a significant association between aging and vTSN strength, especially within the subcortex-limbic system and between the subcortex-limbic system and the dorsal neocortex. In summary, our proposed vTSN is robust in uncovering individual variability and neurobiological brain processes, which can serve as biologically plausible measures for linking biological processes and human behavior.

Enhanced cerebral blood flow similarity of the somatomotor network in chronic insomnia: Transcriptomic decoding, gut microbial signatures and phenotypic roles

Chronic insomnia (CI) is a complex disease involving multiple factors including genetics, gut microbiota, and brain structure and function. However, there lacks a unified framework to elucidate how these factors interact in CI. By combining data of clinical assessment, sleep behavior recording, cognitive test, multimodal MRI (structural, functional, and perfusion), gene, and gut microbiota, this study demonstrated that enhanced cerebral blood flow (CBF) similarities of the somatomotor network (SMN) acted as a key mediator to link multiple factors in CI. Specifically, we first demonstrated that only CBF but not morphological or functional networks exhibited alterations in patients with CI, characterized by increases within the SMN and between the SMN and higher-order associative networks. Moreover, these findings were highly reproducible and the CBF similarity method was test-retest reliable. Further, we showed that transcriptional profiles explained 60.4 % variance of the pattern of the increased CBF similarities with the most correlated genes enriched in regulation of cellular and protein localization and material transport, and gut microbiota explained 69.7 % inter-individual variance in the increased CBF similarities with the most contributions from Negativicutes and Lactobacillales. Finally, we found that the increased CBF similarities were correlated with clinical variables, accounted for sleep behaviors and cognitive deficits, and contributed the most to the patient-control classification (accuracy = 84.4 %). Altogether, our findings have important implications for understanding the neuropathology of CI and may inform ways of developing new therapeutic strategies for the disease.

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.

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.

Unraveling multi-scale neuroimaging biomarkers and molecular foundations for schizophrenia: A combined multivariate pattern analysis and transcriptome-neuroimaging association study

AIMS

Schizophrenia is characterized by alterations in resting-state spontaneous brain activity; however, it remains uncertain whether variations at diverse spatial scales are capable of effectively distinguishing patients from healthy controls. Additionally, the genetic underpinnings of these alterations remain poorly elucidated. We aimed to address these questions in this study to gain better understanding of brain alterations and their underlying genetic factors in schizophrenia.

METHODS

A cohort of 103 individuals with diagnosed schizophrenia and 110 healthy controls underwent resting-state functional MRI scans. Spontaneous brain activity was assessed using the regional homogeneity (ReHo) metric at four spatial scales: voxel-level (Scale 1) and regional-level (Scales 2-4: 272, 53, 17 regions, respectively). For each spatial scale, multivariate pattern analysis was performed to classify schizophrenia patients from healthy controls, and a transcriptome-neuroimaging association analysis was performed to establish connections between gene expression data and ReHo alterations in schizophrenia.

RESULTS

The ReHo metrics at all spatial scales effectively discriminated schizophrenia from healthy controls. Scale 2 showed the highest classification accuracy at 84.6%, followed by Scale 1 (83.1%) and Scale 3 (78.5%), while Scale 4 exhibited the lowest accuracy (74.2%). Furthermore, the transcriptome-neuroimaging association analysis showed that there were not only shared but also unique enriched biological processes across the four spatial scales. These related biological processes were mainly linked to immune responses, inflammation, synaptic signaling, ion channels, cellular development, myelination, and transporter activity.

CONCLUSIONS

This study highlights the potential of multi-scale ReHo as a valuable neuroimaging biomarker in the diagnosis of schizophrenia. By elucidating the complex molecular basis underlying the ReHo alterations of this disorder, this study not only enhances our understanding of its pathophysiology, but also pave the way for future advancements in genetic diagnosis and treatment of schizophrenia.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Relationship between brain white matter damage and grey matter atrophy in hereditary spastic paraplegia types 4 and 5

BACKGROUND AND PURPOSE

White matter (WM) damage is the main target of hereditary spastic paraplegia (HSP), but mounting evidence indicates that genotype-specific grey matter (GM) damage is not uncommon. Our aim was to identify and compare brain GM and WM damage patterns in HSP subtypes and investigate how gene expression contributes to these patterns, and explore the relationship between GM and WM damage.

METHODS

In this prospective single-centre cohort study from 2019 to 2022, HSP patients and controls underwent magnetic resonance imaging evaluations. The alterations of GM and WM patterns were compared between groups by applying a source-based morphometry approach. Spearman rank correlation was used to explore the associations between gene expression and GM atrophy patterns in HSP subtypes. Mediation analysis was conducted to investigate the interplay between GM and WM damage.

RESULTS

Twenty-one spastic paraplegia type 4 (SPG4) patients (mean age 50.7 years ± 12.0 SD, 15 men), 21 spastic paraplegia type 5 (SPG5) patients (mean age 29.1 years ± 12.8 SD, 14 men) and 42 controls (sex- and age-matched) were evaluated. Compared to controls, SPG4 and SPG5 showed similar WM damage but different GM atrophy patterns. GM atrophy patterns in SPG4 and SPG5 were correlated with corresponding gene expression (ρ = 0.30, p = 0.008, ρ = 0.40, p < 0.001, respectively). Mediation analysis indicated that GM atrophy patterns were mediated by WM damage in HSP.

CONCLUSIONS

Grey matter atrophy patterns were distinct between SPG4 and SPG5 and were not only secondary to WM damage but also associated with disease-related gene expression.

CLINICAL TRIAL REGISTRATION NO

NCT04006418.

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.

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.

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.

MR Imaging Findings in Anti-Leucine-Rich Glioma Inactivated Protein 1 Encephalitis: A Systematic Review and Meta-analysis

BACKGROUND

Antibodies against leucine-rich glioma inactivated protein 1 (LGI1) constitute a common form of autoimmune encephalitis. On MR imaging, it may show T2 FLAIR hyperintensities of the medial temporal lobe (T2 FLAIR-MTL), involve the basal ganglia, or be unremarkable.

PURPOSE

We performed a systematic review and meta-analysis to obtain prevalence estimates of abnormal findings on MR imaging in anti-LGI1 encephalitis. A human brain map of the LGI1 microarray gene expression was derived from the Allen Human Brain Atlas.

DATA SOURCES

PubMed and Web of Science were searched with the terms "LGI1" and "encephalitis" from inception to April 7, 2022.

STUDY SELECTION

Thirty-one research publications, encompassing case series and retrospective cohort and case-control studies, with >10 patients with anti-LGI1 encephalitis and MR imaging data were included.

DATA ANALYSIS

Pooled prevalence estimates were calculated using Freeman-Tukey double-arcsine transformation. Meta-analysis used DerSimonian and Laird random effects models.

DATA SYNTHESIS

Of 1318 patients in 30 studies, T2 FLAIR-MTL hyperintensities were present in 54% (95% CI, 0.48-0.60; I2 = 76%). Of 394 patients in 13 studies, 27% showed bilateral (95% CI, 0.19-0.36; I2 = 71%) and 24% unilateral T2 FLAIR-MTL abnormalities (95% CI, 0.17-0.32; I2 = 61%). Of 612 patients in 15 studies, basal ganglia abnormalities were present in 10% (95% CI, 0.06-0.15; I2 = 67%). LGI1 expression was highest in the amygdala, hippocampus, and caudate nucleus.

LIMITATIONS

Only part of the spectrum of MR imaging abnormalities in anti-LGI1 encephalitis could be included in a meta-analysis. MR imaging findings were not the main outcomes in most studies, limiting available information. I2 values ranged from 62% to 76%, representing moderate-to-large heterogeneity.

CONCLUSIONS

T2 FLAIR-MTL hyperintensities were present in around one-half of patients with anti-LGI1. The prevalence of unilateral and bilateral presentations was similar, suggesting unilaterality should raise the suspicion of this disease in the appropriate clinical context. Around 10% of patients showed basal ganglia abnormalities, indicating that special attention should be given to this region. LGI1 regional expression coincided with the most frequently reported abnormal findings on MR imaging. Regional specificity might be partially determined by expression levels of the target protein.

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.

Deciphering the functional specialization of whole-brain spatiomolecular gradients in the adult brain

Cortical arealization arises during neurodevelopment from the confluence of molecular gradients representing patterned expression of morphogens and transcription factors. However, whether similar gradients are maintained in the adult brain remains unknown. Here, we uncover three axes of topographic variation in gene expression in the adult human brain that specifically capture previously identified rostral-caudal, dorsal-ventral, and medial-lateral axes of early developmental patterning. The interaction of these spatiomolecular gradients i) accurately reconstructs the position of brain tissue samples, ii) delineates known functional territories, and iii) can model the topographical variation of diverse cortical features. The spatiomolecular gradients are distinct from canonical cortical axes differentiating the primary sensory cortex from the association cortex, but radiate in parallel with the axes traversed by local field potentials along the cortex. We replicate all three molecular gradients in three independent human datasets as well as two nonhuman primate datasets and find that each gradient shows a distinct developmental trajectory across the lifespan. The gradients are composed of several well-known transcription factors (e.g., PAX6 and SIX3), and a small set of genes shared across gradients are strongly enriched for multiple diseases. Together, these results provide insight into the developmental sculpting of functionally distinct brain regions, governed by three robust transcriptomic axes embedded within brain parenchyma.

Irritable Bowel Syndrome Is Associated With Brain Health by Neuroimaging, Behavioral, Biochemical, and Genetic Analyses

BACKGROUND

Irritable bowel syndrome (IBS) interacts with psychopathology in a complex way; however, little is known about the underlying brain, biochemical, and genetic mechanisms.

METHODS

To clarify the phenotypic and genetic associations between IBS and brain health, we performed a comprehensive retrospective cohort study on a large population. Our study included 171,104 participants from the UK Biobank who underwent a thorough assessment of IBS, with the majority also providing neuroimaging, behavioral, biochemical, and genetic information. Multistage linked analyses were conducted, including phenome-wide association analysis, polygenic risk score calculation, and 2-sample Mendelian randomization analysis.

RESULTS

The phenome-wide association analysis showed that IBS was linked to brain health problems, including anxiety and depression, and poor cognitive performance. Significantly lower brain volumes associated with more severe IBS were found in key areas related to emotional regulation and higher-order cognition, including the medial orbitofrontal cortex/ventromedial prefrontal cortex, anterior insula, anterior and mid-cingulate cortices, dorsolateral prefrontal cortex, and hippocampus. Higher triglycerides, lower high-intensity lipoprotein, and lower platelets were also related (p < 1 × 10-10) to more severe IBS. Finally, Mendelian randomization analyses demonstrated potential causal relationships between IBS and brain health and indicated possible mediating effects of dyslipidemia and inflammation.

CONCLUSIONS

For the first time, this study provides a comprehensive understanding of the relationship between IBS and brain health phenotypes, integrating perspectives from neuroimaging, behavioral performance, biochemical factors, and genetics, which is of great significance for clinical applications to potentially address brain health impairments in patients with IBS.

Quantitative susceptibility mapping in the brain reflects spatial expression of genes involved in iron homeostasis and myelination

Quantitative susceptibility mapping (QSM) is an MRI modality used to non-invasively measure iron content in the brain. Iron exhibits a specific anatomically varying pattern of accumulation in the brain across individuals. The highest regions of accumulation are the deep grey nuclei, where iron is stored in paramagnetic molecule ferritin. This form of iron is considered to be what largely contributes to the signal measured by QSM in the deep grey nuclei. It is also known that QSM is affected by diamagnetic myelin contents. Here, we investigate spatial gene expression of iron and myelin related genes, as measured by the Allen Human Brain Atlas, in relation to QSM images of age-matched subjects. We performed multiple linear regressions between gene expression and the average QSM signal within 34 distinct deep grey nuclei regions. Our results show a positive correlation (p < .05, corrected) between expression of ferritin and the QSM signal in deep grey nuclei regions. We repeated the analysis for other genes that encode proteins thought to be involved in the transport and storage of iron in the brain, as well as myelination. In addition to ferritin, our findings demonstrate a positive correlation (p < .05, corrected) between the expression of ferroportin, transferrin, divalent metal transporter 1, several gene markers of myelinating oligodendrocytes, and the QSM signal in deep grey nuclei regions. Our results suggest that the QSM signal reflects both the storage and active transport of iron in the deep grey nuclei regions of the brain.

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

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

Disrupted topological properties of structural brain networks present a glutamatergic neuropathophysiology in people with narcolepsy

STUDY OBJECTIVES

Growing evidences have documented various abnormalities of the white matter bundles in people with narcolepsy. We sought to evaluate topological properties of brain structural networks, and their association with symptoms and neuropathophysiological features in people with narcolepsy.

METHODS

Diffusion tensor imaging was conducted for people with narcolepsy (n = 30) and matched healthy controls as well as symptoms assessment. Structural connectivity for each participant was generated to analyze global and regional topological properties and their correlations with narcoleptic features. Further human brain transcriptome was extracted and spatially registered for connectivity vulnerability. Genetic functional enrichment analysis was performed and further clarified using in vivo emission computed tomography data.

RESULTS

A wide and dramatic decrease in structural connectivities was observed in people with narcolepsy, with descending network degree and global efficiency. These metrics were not only correlated with sleep latency and awakening features, but also reflected alterations of sleep macrostructure in people with narcolepsy. Network-based statistics identified a small hyperenhanced subnetwork of cingulate gyrus that was closely related to rapid eye movement sleep behavior disorder (RBD) in narcolepsy. Further imaging genetics analysis suggested glutamatergic signatures were responsible for the preferential vulnerability of connectivity alterations in people with narcolepsy, while additional PET/SPECT data verified that structural alteration was significantly correlated with metabotropic glutamate receptor 5 (mGlutR5) and N-methyl-D-aspartate receptor (NMDA).

CONCLUSIONS

People with narcolepsy endured a remarkable decrease in the structural architecture, which was not only closely related to narcolepsy symptoms but also glutamatergic signatures.

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

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

Lipidome atlas of the adult human brain

Lipids are the most abundant but poorly explored components of the human brain. Here, we present a lipidome map of the human brain comprising 75 regions, including 52 neocortical ones. The lipidome composition varies greatly among the brain regions, affecting 93% of the 419 analyzed lipids. These differences reflect the brain's structural characteristics, such as myelin content (345 lipids) and cell type composition (353 lipids), but also functional traits: functional connectivity (76 lipids) and information processing hierarchy (60 lipids). Combining lipid composition and mRNA expression data further enhances functional connectivity association. Biochemically, lipids linked with structural and functional brain features display distinct lipid class distribution, unsaturation extent, and prevalence of omega-3 and omega-6 fatty acid residues. We verified our conclusions by parallel analysis of three adult macaque brains, targeted analysis of 216 lipids, mass spectrometry imaging, and lipidome assessment of sorted murine neurons.

TRIM37 is a primate-specific E3 ligase for Huntingtin and accounts for the striatal degeneration in Huntington's disease

Huntington's disease (HD) is an autosomal dominant neurodegenerative disease characterized by preferential neuronal loss in the striatum. The mechanism underlying striatal selective neurodegeneration remains unclear, making it difficult to develop effective treatments for HD. In the brains of nonhuman primates, we examined the expression of Huntingtin (HTT), the gene responsible for HD. We found that HTT protein is highly expressed in striatal neurons due to its slow degradation in the striatum. We also identified tripartite motif-containing 37 (TRIM37) as a primate-specific protein that interacts with HTT and is selectively reduced in the primate striatum. TRIM37 promotes the ubiquitination and degradation of mutant HTT (mHTT) in vitro and modulates mHTT aggregation in mouse and monkey brains. Our findings suggest that nonhuman primates are crucial for understanding the mechanisms of human diseases such as HD and support TRIM37 as a potential therapeutic target for treating HD.

Similar imaging changes and their relations to genetic profiles in bipolar disorder across different clinical stages

Bipolar disorder (BD) across different clinical stages may present shared and distinct changes in brain activity. We aimed to reveal the neuroimaging homogeneity and heterogeneity of BD and its relationship with clinical variables and genetic variations. In present study, we conducted fractional amplitude of low-frequency fluctuations (fALFF), functional connectivity (FC) and genetic neuroimaging association analyses with 32 depressed, 26 manic, 35 euthymic BD patients and 87 healthy controls (HCs). Significant differences were found in the bilateral pre/subgenual anterior cingulate cortex (ACC) across the four groups, and all bipolar patients exhibited decreased fALFF values in the ACC when compared to HCs. Furthermore, positive associations were significantly observed between fALFF values in the pre/subgenual ACC and participants' cognitive functioning. No significant changes were found in ACC-based FC. We identified fALFF-alteration-related genes in BD, with enrichment in biological progress including synaptic and ion transmission. Taken together, abnormal activity in ACC is a characteristic change associated with BD, regardless of specific mood stages, serving as a potential neuroimaging feature in BD patients. Our genetic neuroimaging association analysis highlights possible heterogeneity in biological processes that could be responsible for different clinical stages in BD.

An evolutionary timeline of the oxytocin signaling pathway

Oxytocin is a neuropeptide associated with both psychological and somatic processes like parturition and social bonding. Although oxytocin homologs have been identified in many species, the evolutionary timeline of the entire oxytocin signaling gene pathway has yet to be described. Using protein sequence similarity searches, microsynteny, and phylostratigraphy, we assigned the genes supporting the oxytocin pathway to different phylostrata based on when we found they likely arose in evolution. We show that the majority (64%) of genes in the pathway are 'modern'. Most of the modern genes evolved around the emergence of vertebrates or jawed vertebrates (540 - 530 million years ago, 'mya'), including OXTR, OXT and CD38. Of those, 45% were under positive selection at some point during vertebrate evolution. We also found that 18% of the genes in the oxytocin pathway are 'ancient', meaning their emergence dates back to cellular organisms and opisthokonta (3500-1100 mya). The remaining genes (18%) that evolved after ancient and before modern genes were classified as 'medium-aged'. Functional analyses revealed that, in humans, medium-aged oxytocin pathway genes are highly expressed in contractile organs, while modern genes in the oxytocin pathway are primarily expressed in the brain and muscle tissue.

Multimodal Nature of the Single-cell Primate Brain Atlas: Morphology, Transcriptome, Electrophysiology, and Connectivity

Primates exhibit complex brain structures that augment cognitive function. The neocortex fulfills high-cognitive functions through billions of connected neurons. These neurons have distinct transcriptomic, morphological, and electrophysiological properties, and their connectivity principles vary. These features endow the primate brain atlas with a multimodal nature. The recent integration of next-generation sequencing with modified patch-clamp techniques is revolutionizing the way to census the primate neocortex, enabling a multimodal neuronal atlas to be established in great detail: (1) single-cell/single-nucleus RNA-seq technology establishes high-throughput transcriptomic references, covering all major transcriptomic cell types; (2) patch-seq links the morphological and electrophysiological features to the transcriptomic reference; (3) multicell patch-clamp delineates the principles of local connectivity. Here, we review the applications of these technologies in the primate neocortex and discuss the current advances and tentative gaps for a comprehensive understanding of the primate neocortex.

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.

Functional brain network controllability dysfunction in Alzheimer's disease and its relationship with cognition and gene expression profiling

Objective. In recent studies, network control theory has been applied to clarify transitions between brain states, emphasizing the significance of assessing the controllability of brain networks in facilitating transitions from one state to another. Despite these advancements, the potential alterations in functional network controllability associated with Alzheimer's disease (AD), along with the underlying genetic mechanisms responsible for these alterations, remain unclear.Approach. We conducted a comparative analysis of functional network controllability measures between patients with AD (n= 64) and matched normal controls (NCs,n= 64). We investigated the association between altered controllability measures and cognitive function in AD. Additionally, we conducted correlation analyses in conjunction with the Allen Human Brain Atlas to identify genes whose expression was correlated with changes in functional network controllability in AD, followed by a set of analyses on the functional features of the identified genes.Main results. In comparison to NCs, patients with AD exhibited a reduction in average controllability, predominantly within the default mode network (DMN) (63% of parcellations), and an increase in average controllability within the limbic (LIM) network (33% of parcellations). Conversely, AD patients displayed a decrease in modal controllability within the LIM network (27% of parcellations) and an increase in modal controllability within the DMN (80% of parcellations). In AD patients, a significant positive correlation was found between the average controllability of the salience network and the mini-mental state examination scores. The changes in controllability measures exhibited spatial correlation with transcriptome profiles. The significant genes identified exhibited enrichment in neurobiologically relevant pathways and demonstrated preferential expression in various tissues, cell types, and developmental periods.Significance. Our findings have the potential to offer new insights into the genetic mechanisms underlying alterations in the controllability of functional networks in AD. Additionally, these results offered perspectives for a deeper understanding of the pathogenesis and the development of therapeutic strategies for AD.

Co-representation of Functional Brain Networks Is Shaped by Cortical Myeloarchitecture and Reveals Individual Behavioral Ability

Large-scale functional networks are spatially distributed in the human brain. Despite recent progress in differentiating their functional roles, how the brain navigates the spatial coordination among them and the biological relevance of this coordination is still not fully understood. Capitalizing on canonical individualized networks derived from functional MRI data, we proposed a new concept, that is, co-representation of functional brain networks, to delineate the spatial coordination among them. To further quantify the co-representation pattern, we defined two indexes, that is, the co-representation specificity (CoRS) and intensity (CoRI), for separately measuring the extent of specific and average expression of functional networks at each brain location by using the data from both sexes. We found that the identified pattern of co-representation was anchored by cortical regions with three types of cytoarchitectural classes along a sensory-fugal axis, including, at the first end, primary (idiotypic) regions showing high CoRS, at the second end, heteromodal regions showing low CoRS and high CoRI, at the third end, paralimbic regions showing low CoRI. Importantly, we demonstrated the critical role of myeloarchitecture in sculpting the spatial distribution of co-representation by assessing the association with the myelin-related neuroanatomical and transcriptomic profiles. Furthermore, the significance of manifesting the co-representation was revealed in its prediction of individual behavioral ability. Our findings indicated that the spatial coordination among functional networks was built upon an anatomically configured blueprint to facilitate neural information processing, while advancing our understanding of the topographical organization of the brain by emphasizing the assembly of functional networks.

Abundant pleiotropy across neuroimaging modalities identified through a multivariate genome-wide association study

Genetic pleiotropy is abundant across spatially distributed brain characteristics derived from one neuroimaging modality (e.g. structural, functional or diffusion magnetic resonance imaging [MRI]). A better understanding of pleiotropy across modalities could inform us on the integration of brain function, micro- and macrostructure. Here we show extensive genetic overlap across neuroimaging modalities at a locus and gene level in the UK Biobank (N = 34,029) and ABCD Study (N = 8607). When jointly analysing phenotypes derived from structural, functional and diffusion MRI in a genome-wide association study (GWAS) with the Multivariate Omnibus Statistical Test (MOSTest), we boost the discovery of loci and genes beyond previously identified effects for each modality individually. Cross-modality genes are involved in fundamental biological processes and predominantly expressed during prenatal brain development. We additionally boost prediction of psychiatric disorders by conditioning independent GWAS on our multimodal multivariate GWAS. These findings shed light on the shared genetic mechanisms underlying variation in brain morphology, functional connectivity, and tissue composition.

Exploring functional dysconnectivity in schizophrenia: alterations in eigenvector centrality mapping and insights into related genes from transcriptional profiles

Schizophrenia is a mental health disorder characterized by functional dysconnectivity. Eigenvector centrality mapping (ECM) has been employed to investigate alterations in functional connectivity in schizophrenia, yet the results lack consistency, and the genetic mechanisms underlying these changes remain unclear. In this study, whole-brain voxel-wise ECM analyses were conducted on resting-state functional magnetic resonance imaging data. A cohort of 91 patients with schizophrenia and 91 matched healthy controls were included during the discovery stage. Additionally, in the replication stage, 153 individuals with schizophrenia and 182 healthy individuals participated. Subsequently, a comprehensive analysis was performed using an independent transcriptional database derived from six postmortem healthy adult brains to explore potential genetic factors influencing the observed functional dysconnectivity, and to investigate the roles of identified genes in neural processes and pathways. The results revealed significant and reliable alterations in the ECM across multiple brain regions in schizophrenia. Specifically, there was a significant decrease in ECM in the bilateral superior and middle temporal gyrus, and an increase in the bilateral thalamus in both the discovery and replication stages. Furthermore, transcriptional analysis revealed 420 genes whose expression patterns were related to changes in ECM, and these genes were enriched mainly in biological processes associated with synaptic signaling and transmission. Together, this study enhances our knowledge of the neural processes and pathways involved in schizophrenia, shedding light on the genetic factors that may be linked to functional dysconnectivity in this disorder.