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.

Regional homogeneity alterations in patients with functional constipation and their associations with gene expression profiles.

Functional constipation, a highly prevalent functional gastrointestinal disorder, often accompanies by mental and psychological disorders. Previous neuroimaging studies have demonstrated brain functional and structural alterations in patients with functional constipation. However, little is known about whether and how regional homogeneity is altered in these patients. Moreover, the potential genetic mechanisms associated with these alterations remain largely unknown. The study included 73 patients with functional constipation and 68 healthy controls, and regional homogeneity comparison was conducted to identify the abnormal spontaneous brain activities in patients with functional constipation. Using Allen Human Brain Atlas, we further investigated gene expression profiles associated with regional homogeneity alterations in functional constipation patients with partial least squares regression analysis applied. Compared with healthy controls, functional constipation patients demonstrated significantly decreased regional homogeneity in both bilateral caudate nucleus, putamen, anterior insula, thalamus and right middle cingulate cortex, supplementary motor area, and increased regional homogeneity in the bilateral orbitofrontal cortex. Genes related to synaptic signaling, central nervous system development, fatty acid metabolism, and immunity were spatially correlated with abnormal regional homogeneity patterns. Our findings showed significant regional homogeneity alterations in functional constipation patients, and the changes may be caused by complex polygenetic and poly-pathway mechanisms, which provides a new perspective on functional constipation's pathophysiology.

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.

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.

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.

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.

Genetic mechanisms underlying gray matter volume changes in patients with drug-naive first-episode schizophrenia.

Brain structural damage is a typical feature of schizophrenia. Investigating such disease phenotype in patients with drug-naive first-episode schizophrenia (DFSZ) may exclude the confounds of antipsychotics and illness chronicity. However, small sample sizes and marked clinical heterogeneity have precluded definitive identification of gray matter volume (GMV) changes in DFSZ as well as their underlying genetic mechanisms. Here, GMV changes in DFSZ were assessed using a neuroimaging meta-analysis of 19 original studies, including 605 patients and 637 controls. Gene expression data were derived from the Allen Human Brain Atlas and processed with a newly proposed standardized pipeline. Then, we used transcriptome-neuroimaging spatial correlations to identify genes associated with GMV changes in DFSZ, followed by a set of gene functional feature analyses. Meta-analysis revealed consistent GMV reduction in the right superior temporal gyrus, right insula and left inferior temporal gyrus in DFSZ. Moreover, we found that these GMV changes were spatially correlated with expression levels of 1,201 genes, which exhibited a wide range of functional features. Our findings may provide important insights into the genetic mechanisms underlying brain morphological abnormality in schizophrenia.

Cortical gradient perturbation in attention deficit hyperactivity disorder correlates with neurotransmitter-, cell type-specific and chromosome- transcriptomic signatures.

AIMS: This study aimed to illuminate the neuropathological landscape of attention deficit hyperactivity disorder (ADHD) by a multiscale macro-micro-molecular perspective from in vivo neuroimaging data. METHODS: The "ADHD-200 initiative" repository provided multi-site high-quality resting-state functional connectivity (rsfc-) neuroimaging for ADHD children and matched typically developing (TD) cohort. Diffusion mapping embedding model to derive the functional connectome gradient detecting biologically plausible neural pattern was built, and the multivariate partial least square method to uncover the enrichment of neurotransmitomic, cellular and chromosomal gradient-transcriptional signatures of AHBA enrichment and meta-analytic decoding. RESULTS: Compared to TD, ADHD children presented connectopic cortical gradient perturbations in almost all the cognition-involved brain macroscale networks (all pBH <0.001), but not in the brain global topology. As an intermediate phenotypic variant, such gradient perturbation was spatially enriched into distributions of GABAA/BZ and 5-HT2A receptors (all pBH <0.01) and co-varied with genetic transcriptional expressions (e.g. DYDC2, ATOH7, all pBH <0.01), associated with phenotypic variants in episodic memory and emotional regulations. Enrichment models demonstrated such gradient-transcriptional variants indicated the risk of both cell-specific and chromosome- dysfunctions, especially in enriched expression of oligodendrocyte precursors and endothelial cells (all pperm <0.05) as well enrichment into chromosome 18, 19 and X (pperm <0.05). CONCLUSIONS: Our findings bridged brain macroscale neuropathological patterns to microscale/cellular biological architectures for ADHD children, demonstrating the neurobiologically pathological mechanism of ADHD into the genetic and molecular variants in GABA and 5-HT systems as well brain-derived enrichment of specific cellular/chromosomal expressions.

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

BACKGROUND: Neuroimaging techniques provide insights into the brain abnormalities secondary to degenerative cervical myelopathy (DCM) and their association with neurological deficits. However, the neural correlates underlying the discrepancy between symptom severity and the degree of spinal cord compression, as well as the transcriptional correlates of these cortical abnormalities, remain unknown in DCM patients. METHODS: In this cross-sectional study, which collected resting-state functional MRI (rs-fMRI) images and the Japanese Orthopedic Association (JOA) score, enrolled 104 participants (54 patients and 50 healthy controls). The frequency-dependent amplitude of low-frequency fluctuation (ALFF) was obtained for all participants. We investigated the ALFF differences between mild-symptom DCM patients and severe-symptom DCM patients while carefully matching the degree of compression between these two groups via both univariate comparison and searchlight classification for three frequency bands (e.g., Slow-4, Slow-5, and Full-band). Additionally, we identified genes associated with symptom severity in DCM patients by linking the spatial patterns of gene expression of Allen Human Brain Atlas and brain functional differences between mild symptom and severe symptom groups. RESULTS: (1) We found that the frequency-specific brain activities within the sensorimotor network (SMN), visual network (VN), and default mode network (DMN) were associated with the varying degrees of functional impairment in DCM patients; (2) the frequency-specific brain activity within the SMN correlated with the functional recovery in patients with DCM; (3) a spatial correlation between the brain-wide expression of genes involved in neuronal migration and the brain functional activities associated with symptom severity was identified in DCM patients. CONCLUSION: In conclusion, our study bridges gaps between genes, cell classes, biological processes, and brain functional correlates of DCM. While our findings are correlational in nature, they suggest that the neural activities of sensorimotor cortices in DCM are associated with the severity of symptoms and might be associated with neuronal migration within the brain.

Association of protein distribution and gene expression revealed by positron emission tomography and postmortem gene expression in the dopaminergic system of the human brain.

PURPOSE: The topological distribution of dopamine-related proteins is determined by gene transcription and subsequent regulations. Recent research strategies integrating positron emission tomography with a transcriptome atlas have opened new opportunities to understand the influence of regulation after transcription on protein distribution. Previous studies have reported that messenger (m)-RNA expression levels spatially correlate with the density maps of serotonin receptors but not with those of transporters. This discrepancy may be due to differences in regulation after transcription between presynaptic and postsynaptic proteins, which have not been studied in the dopaminergic system. Here, we focused on dopamine D1 and D2/D3 receptors and dopamine transporters and investigated their region-wise relationship between mRNA expression and protein distribution. METHODS: We examined the region-wise correlation between regional binding potentials of the target region relative to that of non-displaceable tissue (BPND) values of 11C-SCH-23390 and mRNA expression levels of dopamine D1 receptors (D1R); regional BPND values of 11C-FLB-457 and mRNA expression levels of dopamine D2/D3 receptors (D2/D3R); and regional total distribution volume (VT) values of 18F-FE-PE2I and mRNA expression levels of dopamine transporters (DAT) using Spearman's rank correlation. RESULTS: We found significant positive correlations between regional BPND values of 11C-SCH-23390 and the mRNA expression levels of D1R (r = 0.769, p = 0.0021). Similar to D1R, regional BPND values of 11C-FLB-457 positively correlated with the mRNA expression levels of D2R (r = 0.809, p = 0.0151) but not with those of D3R (r = 0.413, p = 0.3095). In contrast to D1R and D2R, no significant correlation between VT values of 18F-FE-PE2I and mRNA expression levels of DAT was observed (r = -0.5934, p = 0.140). CONCLUSION: We found a region-wise correlation between the mRNA expression levels of dopamine D1 and D2 receptors and their respective protein distributions. However, we found no region-wise correlation between the mRNA expression levels of dopamine transporters and their protein distributions, indicating different regulatory mechanisms for the localization of pre- and postsynaptic proteins. These results provide a broader understanding of the application of the transcriptome atlas to neuroimaging studies of the dopaminergic nervous system.

Differences in whole-brain metabolism are associated with the expression of genes related to neurovascular unit integrity and synaptic plasticity in temporal lobe epilepsy.

PURPOSE: Temporal lobe epilepsy (TLE) is a common, polygenic epilepsy syndrome that involves glucose hypometabolism in the epileptogenic zone. However, the transcriptional and cellular signatures underlying the metabolism in TLE remain unclear. METHODS: In this retrospective study, 2-[18F]-fluoro-2-deoxy-D-glucose ([18F]FDG) positron emission tomography (PET) scans of TLE patients (n = 104) who underwent anterior temporal lobectomy were consecutively collected between 2016 and 2021. The transcriptional profiles of TLE risk genes across the brain were identified by the gene expression analyses from six TLE patients and twelve postmortem donors (six from the Allen Human Brain Atlas). Integrating the neuroimaging and transcriptomic data, we examined the relationship between the expression of TLE-associated genes and metabolic alterations in TLE. Furthermore, we performed functional enrichment analyses of the genes with higher weight in partial least squares regression using Metascape. RESULTS: A total of 104 patients with TLE (mean age 29 ± 9 years, 50% male) and 30 healthy controls (HCs) (mean age 31 ± 6 years, 53% male) were enrolled. Compared to that of HCs, patients with TLE showed hypometabolism in the temporal lobes and adjacent structures but hypermetabolism in the thalamus and basal ganglia. The cortical map of inter-group differences in cerebral metabolism was spatially correlated with the expression of a weighted combination of genes enriched in ontology terms and pathways related to neurovascular unit (NVU) integrity and synaptic plasticity. DISCUSSION: Our findings, combined with the analysis of neuroimaging and transcriptional data, suggest that genes related to NVU integrity and synaptic plasticity may drive alterations to brain metabolism that mediate the genetic risk of TLE.

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.

Coexpression of Gene Transcripts with Monoamine Oxidase A Quantified by Human In Vivo Positron Emission Tomography.

The monoamine oxidase A (MAO-A) is integral to monoamine metabolism and is thus relevant to the pathophysiology of various neuropsychiatric disorders; however, associated gene-enzyme relations are not well understood. This study aimed to unveil genes coexpressed with MAO-A. Therefore, 18 179 mRNA expression maps (based on the Allen Human Brain Atlas) were correlated with the cerebral distribution volume (VT) of MAO-A assessed in 36 healthy subjects (mean age ± standard deviation: 32.9 ± 8.8 years, 18 female) using [11C]harmine positron emission tomography scans. Coexpression analysis was based on Spearman's ρ, over-representation tests on Fisher's exact test with false discovery rate (FDR) correction. The analysis revealed 35 genes in cortex (including B-cell translocation gene family, member 3, implicated in neuroinflammation) and 247 genes in subcortex (including kallikrein-related peptidase 10, implicated in Alzheimer's disease). Significantly over-represented Gene Ontology terms included "neuron development", "neuron differentiation", and "cell-cell signaling" as well as "axon" and "neuron projection". In vivo MAO-A enzyme distribution and MAOA expression did not correlate in cortical areas (ρ = 0.08) while correlation was found in subcortical areas (ρ = 0.52), suggesting influences of region-specific post-transcriptional and -translational modifications. The herein reported information could contribute to guide future genetic studies, deepen the understanding of associated pathomechanisms and assist in the pursuit of novel therapeutic targets.

Frequency-dependent genetic modulation of neuronal oscillations: a combined transcriptome and resting-state functional MRI study.

Neuronal oscillations within certain frequency bands are assumed to associate with specific neural processes and cognitive functions. To examine this hypothesis, transcriptome-neuroimaging spatial correlation analysis was applied to resting-state functional magnetic resonance imaging data from 793 healthy individuals and gene expression data from the Allen Human Brain Atlas. We found that expression measures of 336 genes were correlated with fractional amplitude of low-frequency fluctuations (fALFF) in the slow-4 band (0.027-0.073 Hz), whereas there were no expression-fALFF correlations for the other frequency bands. Furthermore, functional enrichment analyses showed that these slow-4 fALFF-related genes were mainly enriched for ion channel, synaptic function, and neuronal system as well as many neuropsychiatric disorders. Specific expression analyses demonstrated that these genes were specifically expressed in brain tissue, in neurons, and during the late stage of cortical development. Concurrently, the fALFF-related genes were linked to multiple behavioral domains, including dementia, attention, and emotion. In addition, these genes could construct a protein-protein interaction network supported by 30 hub genes. Our findings of a frequency-dependent genetic modulation of spontaneous neuronal activity may support the concept that neuronal oscillations within different frequency bands capture distinct neurobiological processes from the perspective of underlying molecular mechanisms.

Conservative and disruptive modes of adolescent change in human brain functional connectivity.

Adolescent changes in human brain function are not entirely understood. Here, we used multiecho functional MRI (fMRI) to measure developmental change in functional connectivity (FC) of resting-state oscillations between pairs of 330 cortical regions and 16 subcortical regions in 298 healthy adolescents scanned 520 times. Participants were aged 14 to 26 y and were scanned on 1 to 3 occasions at least 6 mo apart. We found 2 distinct modes of age-related change in FC: "conservative" and "disruptive." Conservative development was characteristic of primary cortex, which was strongly connected at 14 y and became even more connected in the period from 14 to 26 y. Disruptive development was characteristic of association cortex and subcortical regions, where connectivity was remodeled: connections that were weak at 14 y became stronger during adolescence, and connections that were strong at 14 y became weaker. These modes of development were quantified using the maturational index (MI), estimated as Spearman's correlation between edgewise baseline FC (at 14 y, [Formula: see text]) and adolescent change in FC ([Formula: see text]), at each region. Disruptive systems (with negative MI) were activated by social cognition and autobiographical memory tasks in prior fMRI data and significantly colocated with prior maps of aerobic glycolysis (AG), AG-related gene expression, postnatal cortical surface expansion, and adolescent shrinkage of cortical thickness. The presence of these 2 modes of development was robust to numerous sensitivity analyses. We conclude that human brain organization is disrupted during adolescence by remodeling of FC between association cortical and subcortical areas.

Correlation of receptor density and mRNA expression patterns in the human cerebral cortex.

Changes in distribution of associated molecular targets have been reported across several neuropsychiatric disorders. However, the high-resolution topology of most proteins is unknown and simultaneous in vivo measurement in multi-receptor systems is complicated. To account for the missing proteomic information, messenger ribonucleic acid (mRNA) transcripts are typically used as a surrogate. Nonetheless, post-transcriptional and post-translational processes might cause the discrepancy between the final distribution of proteins and gene expression patterns. Therefore, this study aims to investigate ex vivo links between mRNA expression and corresponding receptor density in the human cerebral cortex. To this end, autoradiography data on the density of 15 different receptors in 38 brain regions were correlated with the expression patterns of 50 associated genes derived from microarray data (mA), RNA sequencing data (RNA-Seq) provided by the Allen Human Brain Atlas and predicted mRNA expression patterns (pred-mRNA). Spearman's rank correlation was used to evaluate the possible links between proteomic data and mRNA expression patterns. Correlations between mRNA and protein density varied greatly between targets: Positive associations were found for e.g. the serotonin 1A (pred-mRNA: rs = 0.708; mA: rs = 0.601) or kainate receptor (pred-mRNA: rs = 0.655; mA: rs = 0.601; RNA-Seq: rs = 0.575) as well as a few negative associations e.g. γ-Aminobutyric acid (GABA) A receptor subunit α3 (pred-mRNA: rs = -0.638; mA: rs = -0.619) or subunit α5 (pred-mRNA: rs = -0.565; mA: rs = -0.563), while most of the other investigated target receptors showed low correlations. The high variability in the correspondence of mRNA expression and receptor spatial distribution warrants caution when inferring the topology of molecular targets in the brain from transcriptome data. This not only highlights the longstanding value of molecular imaging but also indicates a need for comprehensive proteomic studies.

Combined analysis of cytoarchitectonic, molecular and transcriptomic patterns reveal differences in brain organization across human functional brain systems.

Brain areas show specific cellular, molecular, and gene expression patterns that are linked to function, but their precise relationships are largely unknown. To unravel these structure-function relationships, a combined analysis of 53 neurotransmitter receptor genes, receptor densities of six transmitter systems and cytoarchitectonic data of the auditory, somatosensory, visual, motor systems was conducted. Besides covariation of areal gene expression with receptor density, the study reveals specific gene expression patterns in functional systems, which are most prominent for the inhibitory GABAA and excitatory glutamatergic NMDA receptors. Furthermore, gene expression-receptor relationships changed in a systematic manner according to information flow from primary to higher associative areas. The findings shed new light on the relationship of anatomical, functional, and molecular and transcriptomic principles of cortical segregation towards a more comprehensive understanding of human brain organization.

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.

Brain Gene Expression Pattern Correlated with the Differential Brain Activation by Pain and Touch in Humans.

Genes involved in pain and touch sensations have been studied extensively, but very few studies have tried to link them with neural activities in the brain. Here, we aimed to identify genes preferentially correlated to painful activation patterns by linking the spatial patterns of gene expression of Allen Human Brain Atlas with the pain-elicited neural responses in the human brain, with a parallel, control analysis for identification of genes preferentially correlated to tactile activation patterns. We identified 1828 genes whose expression patterns preferentially correlated to painful activation patterns and 411 genes whose expression patterns preferentially correlated to tactile activation pattern at the cortical level. In contrast to the enrichment for astrocyte and inhibitory synaptic transmission of genes preferentially correlated to tactile activation, the genes preferentially correlated to painful activation were mainly enriched for neuron and opioid- and addiction-related pathways and showed significant overlap with pain-related genes identified in previous studies. These findings not only provide important evidence for the differential genetic architectures of specific brain activation patterns elicited by painful and tactile stimuli but also validate a new approach to studying pain- and touch-related genes more directly from the perspective of neural responses in the human brain.

Cortical patterning of abnormal morphometric similarity in psychosis is associated with brain expression of schizophrenia-related genes.

Schizophrenia has been conceived as a disorder of brain connectivity, but it is unclear how this network phenotype is related to the underlying genetics. We used morphometric similarity analysis of MRI data as a marker of interareal cortical connectivity in three prior case-control studies of psychosis: in total, n = 185 cases and n = 227 controls. Psychosis was associated with globally reduced morphometric similarity in all three studies. There was also a replicable pattern of case-control differences in regional morphometric similarity, which was significantly reduced in patients in frontal and temporal cortical areas but increased in parietal cortex. Using prior brain-wide gene expression data, we found that the cortical map of case-control differences in morphometric similarity was spatially correlated with cortical expression of a weighted combination of genes enriched for neurobiologically relevant ontology terms and pathways. In addition, genes that were normally overexpressed in cortical areas with reduced morphometric similarity were significantly up-regulated in three prior post mortem studies of schizophrenia. We propose that this combined analysis of neuroimaging and transcriptional data provides insight into how previously implicated genes and proteins as well as a number of unreported genes in their topological vicinity on the protein interaction network may drive structural brain network changes mediating the genetic risk of schizophrenia.