Prefrontal Coexpression of Schizophrenia Risk Genes Is Associated With Treatment Response in Patients

Network-wide risk convergence in gene co-expression identifies reproducible genetic hubs of schizophrenia risk

The omnigenic model posits that genetic risk for traits with complex heritability involves cumulative effects of peripheral genes on mechanistic "core genes," suggesting that in a network of genes, those closer to clusters including core genes should have higher GWAS signals. In gene co-expression networks, we confirmed that GWAS signals accumulate in genes more connected to risk-enriched gene clusters, highlighting across-network risk convergence. This was strongest in adult psychiatric disorders, especially schizophrenia (SCZ), spanning 70% of network genes, suggestive of super-polygenic architecture. In snRNA-seq cell type networks, SCZ risk convergence was strongest in L2/L3 excitatory neurons. We prioritized genes most connected to SCZ-GWAS genes, which showed robust association to a CRISPRa measure of PGC3 regulation and were consistently identified across several brain regions. Several genes, including dopamine-associated ones, were prioritized specifically in the striatum. This strategy thus retrieves current drug targets and can be used to prioritize other potential drug targets.

Critical reasoning on the co-expression module QTL in the dorsolateral prefrontal cortex

Expression quantitative trait locus (eQTL) analysis is a popular method of gaining insight into the function of regulatory variation. While cis-eQTL resources have been instrumental in linking genome-wide association study variants to gene function, complex trait heritability may be additionally mediated by other forms of gene regulation. Toward this end, novel eQTL methods leverage gene co-expression (module-QTL) to investigate joint regulation of gene modules by single genetic variants. Here we broadly define a "module-QTL" as the association of a genetic variant with a summary measure of gene co-expression. This approach aims to reduce the multiple testing burden of a trans-eQTL search through the consolidation of gene-based testing and provide biological context to eQTLs shared between genes. In this article we provide an in-depth examination of the co-expression module eQTL (module-QTL) through literature review, theoretical investigation, and real-data application of the module-QTL to three large prefrontal cortex genotype-RNA sequencing datasets. We find module-QTLs in our study that are disease associated and reproducible are not additionally informative beyond cis- or trans-eQTLs for module genes. Through comparison to prior studies, we highlight promises and limitations of the module-QTL across study designs and provide recommendations for further investigation of the module-QTL framework.

Dopamine signaling enriched striatal gene set predicts striatal dopamine synthesis and physiological activity in vivo

The polygenic architecture of schizophrenia implicates several molecular pathways involved in synaptic function. However, it is unclear how polygenic risk funnels through these pathways to translate into syndromic illness. Using tensor decomposition, we analyze gene co-expression in the caudate nucleus, hippocampus, and dorsolateral prefrontal cortex of post-mortem brain samples from 358 individuals. We identify a set of genes predominantly expressed in the caudate nucleus and associated with both clinical state and genetic risk for schizophrenia that shows dopaminergic selectivity. A higher polygenic risk score for schizophrenia parsed by this set of genes predicts greater dopamine synthesis in the striatum and greater striatal activation during reward anticipation. These results translate dopamine-linked genetic risk variation into in vivo neurochemical and hemodynamic phenotypes in the striatum that have long been implicated in the pathophysiology of schizophrenia.

Sex dimorphism controls dysbindin-related cognitive dysfunctions in mice and humans with the contribution of COMT

Cognitive dysfunctions are core-enduring symptoms of schizophrenia, with important sex-related differences. Genetic variants of the DTBPN1 gene associated with reduced dysbindin-1 protein (Dys) expression negatively impact cognitive functions in schizophrenia through a functional epistatic interaction with Catechol-O-methyltransferase (COMT). Dys is involved in the trafficking of dopaminergic receptors, crucial for prefrontal cortex (PFC) signaling regulation. Moreover, dopamine signaling is modulated by estrogens via inhibition of COMT expression. We hypothesized a sex dimorphism in Dys-related cognitive functions dependent on COMT and estrogen levels. Our multidisciplinary approach combined behavioral-molecular findings on genetically modified mice, human postmortem Dys expression data, and in vivo fMRI during a working memory task performance. We found cognitive impairments in male mice related to genetic variants characterized by reduced Dys protein expression (pBonferroni = 0.0001), as well as in male humans through a COMT/Dys functional epistatic interaction involving PFC brain activity during working memory (t(23) = -3.21; pFDR = 0.004). Dorsolateral PFC activity was associated with lower working memory performance in males only (p = 0.04). Also, male humans showed decreased Dys expression in dorsolateral PFC during adulthood (pFDR = 0.05). Female Dys mice showed preserved cognitive performances with deficits only with a lack of estrogen tested in an ovariectomy model (pBonferroni = 0.0001), suggesting that genetic variants reducing Dys protein expression could probably become functional in females when the protective effect of estrogens is attenuated, i.e., during menopause. Overall, our results show the differential impact of functional variants of the DTBPN1 gene interacting with COMT on cognitive functions across sexes in mice and humans, underlying the importance of considering sex as a target for patient stratification and precision medicine in schizophrenia.

Post-GWAS functional analyses of CNTNAP5 suggests its role in glaucomatous neurodegeneration

Primary angle closure glaucoma (PACG) affects more than 20 million people worldwide, with an increased prevalence in south-east Asia. In a prior haplotype-based GWAS, we identified a novel CNTNAP5 genic region, significantly associated with PACG. In the current study, we have extended our perception of CNTNAP5 involvement in glaucomatous neurodegeneration in a zebrafish model, through investigating phenotypic consequences pertinent to retinal degeneration upon knockdown of cntnap5 by translation-blocking morpholinos. While cntnap5 knockdown was successfully validated using an antibody, immunofluorescence followed by western blot analyses in cntnap5-morphant (MO) zebrafish revealed increased expression of acetylated tubulin indicative of perturbed cytoarchitecture of retinal layers. Moreover, significant loss of Nissl substance is observed in the neuro-retinal layers of cntnap5-MO zebrafish eye, indicating neurodegeneration. Additionally, in spontaneous movement behavioural analysis, cntnap5-MO zebrafish have a significantly lower average distance traversed in light phase compared to mismatch-controls, whereas no significant difference was observed in the dark phase, corroborating with vision loss in the cntnap5-MO zebrafish. This study provides the first direct functional evidence of a putative role of CNTNAP5 in visual neurodegeneration.

A miR-137-Related Biological Pathway of Risk for Schizophrenia Is Associated With Human Brain Emotion Processing

BACKGROUND

miR-137 is a microRNA involved in brain development, regulating neurogenesis and neuronal maturation. Genome-wide association studies have implicated miR-137 in schizophrenia risk but do not explain its involvement in brain function and underlying biology. Polygenic risk for schizophrenia mediated by miR-137 targets is associated with working memory, although other evidence points to emotion processing. We characterized the functional brain correlates of miR-137 target genes associated with schizophrenia while disentangling previously reported associations of miR-137 targets with working memory and emotion processing.

METHODS

Using RNA sequencing data from postmortem prefrontal cortex (N = 522), we identified a coexpression gene set enriched for miR-137 targets and schizophrenia risk genes. We validated the relationship of this set to miR-137 in vitro by manipulating miR-137 expression in neuroblastoma cells. We translated this gene set into polygenic scores of coexpression prediction and associated them with functional magnetic resonance imaging activation in healthy volunteers (n1 = 214; n2 = 136; n3 = 2075; n4 = 1800) and with short-term treatment response in patients with schizophrenia (N = 427).

RESULTS

In 4652 human participants, we found that 1) schizophrenia risk genes were coexpressed in a biologically validated set enriched for miR-137 targets; 2) increased expression of miR-137 target risk genes was mediated by low prefrontal miR-137 expression; 3) alleles that predict greater gene set coexpression were associated with greater prefrontal activation during emotion processing in 3 independent healthy cohorts (n1, n2, n3) in interaction with age (n4); and 4) these alleles predicted less improvement in negative symptoms following antipsychotic treatment in patients with schizophrenia.

CONCLUSIONS

The functional translation of miR-137 target gene expression linked with schizophrenia involves the neural substrates of emotion processing.

Temporal changes of gene expression in health, schizophrenia, bipolar disorder, and major depressive disorder

The molecular events underlying the development, manifestation, and course of schizophrenia, bipolar disorder, and major depressive disorder span from embryonic life to advanced age. However, little is known about the early dynamics of gene expression in these disorders due to their relatively late manifestation. To address this, we conducted a secondary analysis of post-mortem prefrontal cortex datasets using bioinformatics and machine learning techniques to identify differentially expressed gene modules associated with aging and the diseases, determine their time-perturbation points, and assess enrichment with expression quantitative trait loci (eQTL) genes. Our findings revealed early, mid, and late deregulation of expression of functional gene modules involved in neurodevelopment, plasticity, homeostasis, and immune response. This supports the hypothesis that multiple hits throughout life contribute to disease manifestation rather than a single early-life event. Moreover, the time-perturbed functional gene modules were associated with genetic loci affecting gene expression, highlighting the role of genetic factors in gene expression dynamics and the development of disease phenotypes. Our findings emphasize the importance of investigating time-dependent perturbations in gene expression before the age of onset in elucidating the molecular mechanisms of psychiatric disorders.

Investigating trait variability of gene co-expression network architecture in brain by controlling for genomic risk of schizophrenia

The effect of schizophrenia (SCZ) genetic risk on gene expression in brain remains elusive. A popular approach to this problem has been the application of gene co-expression network algorithms (e.g., WGCNA). To improve reliability with this method it is critical to remove unwanted sources of variance while also preserving biological signals of interest. In this WCGNA study of RNA-Seq data from postmortem prefrontal cortex (78 neurotypical donors, EUR ancestry), we tested the effects of SCZ genetic risk on co-expression networks. Specifically, we implemented a novel design in which gene expression was adjusted by linear regression models to preserve or remove variance explained by biological signal of interest (GWAS genomic scores for SCZ risk-(GS-SCZ), and genomic scores- GS of height (GS-Ht) as a negative control), while removing variance explained by covariates of non-interest. We calculated co-expression networks from adjusted expression (GS-SCZ and GS-Ht preserved or removed), and consensus between them (representative of a "background" network free of genomic scores effects). We then tested the overlap between GS-SCZ preserved modules and background networks reasoning that modules with reduced overlap would be most affected by GS-SCZ biology. Additionally, we tested these modules for convergence of SCZ risk (i.e., enrichment in PGC3 SCZ GWAS priority genes, enrichment in SCZ risk heritability and relevant biological ontologies. Our results highlight key aspects of GS-SCZ effects on brain co-expression networks, specifically: 1) preserving/removing SCZ genetic risk alters the co-expression modules; 2) biological pathways enriched in modules affected by GS-SCZ implicate processes of transcription, translation and metabolism that converge to influence synaptic transmission; 3) priority PGC3 SCZ GWAS genes and SCZ risk heritability are enriched in modules associated with GS-SCZ effects. Overall, our results indicate that gene co-expression networks that selectively integrate information about genetic risk can reveal novel combinations of biological pathways involved in schizophrenia.

Dopamine and schizophrenia from bench to bedside: Discovery of a striatal co-expression risk gene set that predicts in vivo measures of striatal function

Schizophrenia (SCZ) is characterized by a polygenic risk architecture implicating diverse molecular pathways important for synaptic function. However, how polygenic risk funnels through these pathways to translate into syndromic illness is unanswered. To evaluate biologically meaningful pathways of risk, we used tensor decomposition to characterize gene co-expression in post-mortem brain (of neurotypicals: N=154; patients with SCZ: N=84; and GTEX samples N=120) from caudate nucleus (CN), hippocampus (HP), and dorsolateral prefrontal cortex (DLPFC). We identified a CN-predominant gene set showing dopaminergic selectivity that was enriched for genes associated with clinical state and for genes associated with SCZ risk. Parsing polygenic risk score for SCZ based on this specific gene set (parsed-PRS), we found that greater pathway-specific SCZ risk predicted greater in vivo striatal dopamine synthesis capacity measured by [ 18 F]-FDOPA PET in three independent cohorts of neurotypicals and patients (total N=235) and greater fMRI striatal activation during reward anticipation in two additional independent neurotypical cohorts (total N=141). These results reveal a 'bench to bedside' translation of dopamine-linked genetic risk variation in driving in vivo striatal neurochemical and hemodynamic phenotypes that have long been implicated in the pathophysiology of SCZ.

Prioritization of potential causative genes for schizophrenia in placenta

Our earlier work has shown that genomic risk for schizophrenia converges with early life complications in affecting risk for the disorder and sex-biased neurodevelopmental trajectories. Here, we identify specific genes and potential mechanisms that, in placenta, may mediate such outcomes. We performed TWAS in healthy term placentae (N = 147) to derive candidate placental causal genes that we confirmed with SMR; to search for placenta and schizophrenia-specific associations, we performed an analogous analysis in fetal brain (N = 166) and additional placenta TWAS for other disorders/traits. The analyses in the whole sample and stratifying by sex ultimately highlight 139 placenta and schizophrenia-specific risk genes, many being sex-biased; the candidate molecular mechanisms converge on the nutrient-sensing capabilities of placenta and trophoblast invasiveness. These genes also implicate the Coronavirus-pathogenesis pathway and showed increased expression in placentae from a small sample of SARS-CoV-2-positive pregnancies. Investigating placental risk genes for schizophrenia and candidate mechanisms may lead to opportunities for prevention that would not be suggested by study of the brain alone.

Consensus molecular environment of schizophrenia risk genes in coexpression networks shifting across age and brain regions

Schizophrenia is a neurodevelopmental brain disorder whose genetic risk is associated with shifting clinical phenomena across the life span. We investigated the convergence of putative schizophrenia risk genes in brain coexpression networks in postmortem human prefrontal cortex (DLPFC), hippocampus, caudate nucleus, and dentate gyrus granule cells, parsed by specific age periods (total N = 833). The results support an early prefrontal involvement in the biology underlying schizophrenia and reveal a dynamic interplay of regions in which age parsing explains more variance in schizophrenia risk compared to lumping all age periods together. Across multiple data sources and publications, we identify 28 genes that are the most consistently found partners in modules enriched for schizophrenia risk genes in DLPFC; twenty-three are previously unidentified associations with schizophrenia. In iPSC-derived neurons, the relationship of these genes with schizophrenia risk genes is maintained. The genetic architecture of schizophrenia is embedded in shifting coexpression patterns across brain regions and time, potentially underwriting its shifting clinical presentation.

Incorporating multi-stage diagnosis status to mine associations between genetic risk variants and the multi-modality phenotype network in major depressive disorder

Depression (major depressive disorder, MDD) is a common and serious medical illness. Globally, it is estimated that 5% of adults suffer from depression. Recently, imaging genetics receives growing attention and become a powerful strategy for discoverying the associations between genetic variants (e.g., single-nucleotide polymorphisms, SNPs) and multi-modality brain imaging data. However, most of the existing MDD imaging genetic research studies conducted by clinicians usually utilize simple statistical analysis methods and only consider single-modality brain imaging, which are limited in the deeper discovery of the mechanistic understanding of MDD. It is therefore imperative to utilize a powerful and efficient technology to fully explore associations between genetic variants and multi-modality brain imaging. In this study, we developed a novel imaging genetic association framework to mine the multi-modality phenotype network between genetic risk variants and multi-stage diagnosis status. Specifically, the multi-modality phenotype network consists of voxel node features and connectivity edge features from structural magnetic resonance imaging (sMRI) and resting-state functional magnetic resonance imaging (rs-fMRI). Thereafter, an association model based on multi-task learning strategy was adopted to fully explore the relationship between the MDD risk SNP and the multi-modality phenotype network. The multi-stage diagnosis status was introduced to further mine the relation among the multiple modalities of different subjects. A multi-modality brain imaging data and genotype data were collected by us from two hospitals. The experimental results not only demonstrate the effectiveness of our proposed method but also identify some consistent and stable brain regions of interest (ROIs) biomarkers from the node and edge features of multi-modality phenotype network. Moreover, four new and potential risk SNPs associated with MDD were discovered.

Lessons Learned From Parsing Genetic Risk for Schizophrenia Into Biological Pathways

The clinically heterogeneous presentation of schizophrenia is compounded by the heterogeneity of risk factors and neurobiological correlates of the disorder. Genome-wide association studies in schizophrenia have uncovered a remarkably high number of genetic variants, but the biological pathways they impact upon remain largely unidentified. Among the diverse methodological approaches employed to provide a more granular understanding of genetic risk for schizophrenia, the use of biological labels, such as gene ontologies, regulome approaches, and gene coexpression have all provided novel perspectives into how genetic risk translates into the neurobiology of schizophrenia. Here, we review the salient aspects of parsing polygenic risk for schizophrenia into biological pathways. We argue that parsed scores, compared to standard polygenic risk scores, may afford a more biologically plausible and accurate physiological modeling of the different dimensions involved in translating genetic risk into brain mechanisms, including multiple brain regions, cell types, and maturation stages. We discuss caveats, opportunities, and pitfalls inherent in the parsed risk approach.

A polygenic score indexing a DRD2-related co-expression network is associated with striatal dopamine function

The D2 dopamine receptor (D2R) is the primary site of the therapeutic action of antipsychotics and is involved in essential brain functions relevant to schizophrenia, such as attention, memory, motivation, and emotion processing. Moreover, the gene coding for D2R (DRD2) has been associated with schizophrenia at a genome-wide level. Recent studies have shown that a polygenic co-expression index (PCI) predicting the brain-specific expression of a network of genes co-expressed with DRD2 was associated with response to antipsychotics, brain function during working memory in patients with schizophrenia, and with the modulation of prefrontal cortex activity after pharmacological stimulation of D2 receptors. We aimed to investigate the relationship between the DRD2 gene network and in vivo striatal dopaminergic function, which is a phenotype robustly associated with psychosis and schizophrenia. To this aim, a sample of 92 healthy subjects underwent 18F-DOPA PET and was genotyped for genetic variations indexing the co-expression of the DRD2-related genetic network in order to calculate the PCI for each subject. The PCI was significantly associated with whole striatal dopamine synthesis capacity (p = 0.038). Exploratory analyses on the striatal subdivisions revealed a numerically larger effect size of the PCI on dopamine function for the associative striatum, although this was not significantly different than effects in other sub-divisions. These results are in line with a possible relationship between the DRD2-related co-expression network and schizophrenia and extend it by identifying a potential mechanism involving the regulation of dopamine synthesis. Future studies are needed to clarify the molecular mechanisms implicated in this relationship.

Integrative omics of schizophrenia: from genetic determinants to clinical classification and risk prediction

Schizophrenia (SCZ) is a debilitating neuropsychiatric disorder with high heritability and complex inheritance. In the past decade, successful identification of numerous susceptibility loci has provided useful insights into the molecular etiology of SCZ. However, applications of these findings to clinical classification and diagnosis, risk prediction, or intervention for SCZ have been limited, and elucidating the underlying genomic and molecular mechanisms of SCZ is still challenging. More recently, multiple Omics technologies - genomics, transcriptomics, epigenomics, proteomics, metabolomics, connectomics, and gut microbiomics - have all been applied to examine different aspects of SCZ pathogenesis. Integration of multi-Omics data has thus emerged as an approach to provide a more comprehensive view of biological complexity, which is vital to enable translation into assessments and interventions of clinical benefit to individuals with SCZ. In this review, we provide a broad survey of the single-omics studies of SCZ, summarize the advantages and challenges of different Omics technologies, and then focus on studies in which multiple omics data are integrated to unravel the complex pathophysiology of SCZ. We believe that integration of multi-Omics technologies would provide a roadmap to create a more comprehensive picture of interactions involved in the complex pathogenesis of SCZ, constitute a rich resource for elucidating the potential molecular mechanisms of the illness, and eventually improve clinical assessments and interventions of SCZ to address clinical translational questions from bench to bedside.

The Etiology of Auditory Hallucinations in Schizophrenia: From Multidimensional Levels

Enormous efforts have been made to unveil the etiology of auditory hallucinations (AHs), and multiple genetic and neural factors have already been shown to have their own roles. Previous studies have shown that AHs in schizophrenia vary from those in other disorders, suggesting that they have unique features and possibly distinguishable mechanisms worthy of further investigation. In this review, we intend to offer a comprehensive summary of current findings related to AHs in schizophrenia from aspects of genetics and transcriptome, neurophysiology (neurometabolic and electroencephalogram studies), and neuroimaging (structural and functional magnetic resonance imaging studies and transcriptome–neuroimaging association study). Main findings include gene polymorphisms, glutamate level change, electroencephalographic alterations, and abnormalities of white matter fasciculi, cortical structure, and cerebral activities, especially in multiple regions, including auditory and language networks. More solid and comparable research is needed to replicate and integrate ongoing findings from multidimensional levels.

Limited Association between Schizophrenia Genetic Risk Factors and Transcriptomic Features

Schizophrenia is a polygenic disorder with many genomic regions contributing to schizophrenia risk. The majority of genetic variants associated with schizophrenia lie in the non-coding genome and are thought to contribute to transcriptional regulation. Extensive transcriptomic dysregulation has been detected from postmortem brain samples of schizophrenia-affected individuals. However, the relationship between schizophrenia genetic risk factors and transcriptomic features has yet to be explored. Herein, we examined whether varying gene expression features, including differentially expressed genes (DEGs), co-expression networks, and central hubness of genes, contribute to the heritability of schizophrenia. We leveraged quantitative trait loci and chromatin interaction profiles to identify schizophrenia risk variants assigned to the genes that represent different transcriptomic features. We then performed stratified linkage disequilibrium score regression analysis on these variants to estimate schizophrenia heritability enrichment for different gene expression features. Notably, DEGs and co-expression networks showed nominal heritability enrichment. This nominal association can be partly explained by cellular heterogeneity, as DEGs were associated with the genetic risk of schizophrenia in a cell type-specific manner. Moreover, DEGs were enriched for target genes of schizophrenia-associated transcription factors, suggesting that the transcriptomic signatures of schizophrenia are the result of transcriptional regulatory cascades elicited by genetic risk factors.

Brain network dynamics during working memory are modulated by dopamine and diminished in schizophrenia

Dynamical brain state transitions are critical for flexible working memory but the network mechanisms are incompletely understood. Here, we show that working memory performance entails brain-wide switching between activity states using a combination of functional magnetic resonance imaging in healthy controls and individuals with schizophrenia, pharmacological fMRI, genetic analyses and network control theory. The stability of states relates to dopamine D1 receptor gene expression while state transitions are influenced by D2 receptor expression and pharmacological modulation. Individuals with schizophrenia show altered network control properties, including a more diverse energy landscape and decreased stability of working memory representations. Our results demonstrate the relevance of dopamine signaling for the steering of whole-brain network dynamics during working memory and link these processes to schizophrenia pathophysiology.

Working memory requires the brain to switch between cognitive states and activity patterns. Here, the authors show that the steering of these neural network dynamics is influenced by dopamine D1- and D2-receptor function and altered in schizophrenia.

Identification of Reproducible BCL11A Alterations in Schizophrenia Through Individual-Level Prediction of Coexpression

Previous studies have provided evidence for an alteration of genetic coexpression in schizophrenia (SCZ). However, such analyses have thus far lacked biological specificity for individual genes, which may be critical for identifying illness-relevant effects. Therefore, we applied machine learning to identify gene-specific coexpression differences at the individual subject level and compared these between individuals with SCZ, bipolar disorder, major depressive disorder (MDD), autism spectrum disorder (ASD), and healthy controls. Utilizing transcriptome-wide gene expression data from 21 independent datasets, comprising a total of 9509 participants, we identified a reproducible decrease of BCL11A coexpression across 4 SCZ datasets that showed diagnostic specificity for SCZ when compared with ASD and MDD. We further demonstrate that individual-level coexpression differences can be combined in multivariate coexpression scores that show reproducible illness classification across independent datasets in SCZ and ASD. This study demonstrates that machine learning can capture gene-specific coexpression differences at the individual subject level for SCZ and identify novel biomarker candidates.

A Pattern of Cognitive Deficits Stratified for Genetic and Environmental Risk Reliably Classifies Patients With Schizophrenia From Healthy Control Subjects

BACKGROUND

Schizophrenia risk is associated with both genetic and environmental risk factors. Furthermore, cognitive abnormalities are established core characteristics of schizophrenia. We aim to assess whether a classification approach encompassing risk factors, cognition, and their associations can discriminate patients with schizophrenia (SCZs) from healthy control subjects (HCs). We hypothesized that cognition would demonstrate greater HC-SCZ classification accuracy and that combined gene-environment stratification would improve the discrimination performance of cognition.

METHODS

Genome-wide association study-based genetic, environmental, and neurocognitive classifiers were trained to separate 337 HCs from 103 SCZs using support vector classification and repeated nested cross-validation. We validated classifiers on independent datasets using within-diagnostic (SCZ) and cross-diagnostic (clinically isolated syndrome for multiple sclerosis, another condition with cognitive abnormalities) approaches. Then, we tested whether gene-environment multivariate stratification modulated the discrimination performance of the cognitive classifier in iterative subsamples.

RESULTS

The cognitive classifier discriminated SCZs from HCs with a balanced accuracy (BAC) of 88.7%, followed by environmental (BAC = 65.1%) and genetic (BAC = 55.5%) classifiers. Similar classification performance was measured in the within-diagnosis validation sample (HC-SCZ BACs, cognition = 70.5%; environment = 65.8%; genetics = 49.9%). The cognitive classifier was relatively specific to schizophrenia (HC-clinically isolated syndrome for multiple sclerosis BAC = 56.7%). Combined gene-environment stratification allowed cognitive features to classify HCs from SCZs with 89.4% BAC.

CONCLUSIONS

Consistent with cognitive deficits being core features of the phenotype of SCZs, our results suggest that cognitive features alone bear the greatest amount of information for classification of SCZs. Consistent with genes and environment being risk factors, gene-environment stratification modulates HC-SCZ classification performance of cognition, perhaps providing another target for refining early identification and intervention strategies.

Brain-derived neurotrophic factor and schizophrenia

The brain-derived neurotrophic factor (BDNF) is a secretory growth factor that promotes neuronal proliferation and survival, synaptic plasticity and long-term potentiation in the central nervous system. Brain-derived neurotrophic factor biosynthesis and secretion are chrono-topically regulated processes at the cellular level, accounting for specific localizations and functions. Given its role in regulating brain development and activity, BDNF represents a potentially relevant gene for schizophrenia, and indeed BDNF and its non-synonymous functional variant, rs6265 (C → T, Val → Met) have been widely studied in psychiatric genetics. Human and animal studies have indicated that brain-derived neurotrophic factor is relevant for schizophrenia-related phenotypes, and that: (1) fine-tuned regulation of brain-derived neurotrophic factor secretion and activity is necessary to guarantee brain optimal development and functioning; (2) the Val → Met substitution is associated with impaired activity-dependent secretion of brain-derived neurotrophic factor; (3) disruption of brain-derived neurotrophic factor signaling is associated with altered synaptic plasticity and neurodevelopment. However, genome-wide association studies failed to associate the BDNF locus with schizophrenia, even though a sub-threshold association exists. Here, we will review studies focused on the relationship between the genetic variation of BDNF and schizophrenia, trying to fill the gap between genetic risk per se and insights from molecular biology. A deeper understanding of brain-derived neurotrophic factor biology and of the epigenetic regulation of brain-derived neurotrophic factor and its interactome during development may help clarifying the potential role of this gene in schizophrenia, thus informing development of brain-derived neurotrophic factor-based strategies of prevention and treatment of this disorder.

A Novel Schizophrenia Diagnostic Model Based on Statistically Significant Changes in Gene Methylation in Specific Brain Regions

Objective

The present study identified methylation patterns of schizophrenia- (SCZ-) related genes in different brain regions and used them to construct a novel DNA methylation-based SCZ diagnostic model.

Methods

Four DNA methylation datasets representing different brain regions were downloaded from the Gene Expression Omnibus. The common differentially methylated genes (CDMGs) in all datasets were identified to perform functional enrichment analysis. The differential methylation sites of 10 CDMGs involved in the largest numbers of neurological or psychiatric-related biological processes were used to construct a DNA methylation-based diagnostic model for SCZ in the respective datasets.

Results

A total of 849 CDMGs were identified in the four datasets, but the methylation sites as well as degree of methylation differed across the brain regions. Functional enrichment analysis showed CDMGs were significantly involved in biological processes associated with neuronal axon development, intercellular adhesion, and cell morphology changes and, specifically, in PI3K-Akt, AMPK, and MAPK signaling pathways. Four DNA methylation-based classifiers for diagnosing SCZ were constructed in the four datasets, respectively. The sample recognition efficiency of the classifiers showed an area under the receiver operating characteristic curve of 1.00 in three datasets and >0.9 in one dataset.

Conclusion

DNA methylation patterns in SCZ vary across different brain regions, which may be a useful epigenetic characteristic for diagnosing SCZ. Our novel model based on SCZ-gene methylation shows promising diagnostic power.

Integrative network analysis identifies differential regulation of neuroimmune system in Schizophrenia and Bipolar disorder

Background

Neuropsychiatric disorders such as Schizophrenia (SCZ) and Bipolar disorder (BPD) pose a broad range of problems with different symptoms mainly characterized by some combination of abnormal thoughts, emotions, behaviour, etc. However, in depth molecular and pathophysiological mechanisms among different neuropsychiatric disorders have not been clearly understood yet. We have used RNA-seq data to investigate unique and overlapping molecular signatures between SCZ and BPD using an integrative network biology approach.

Methods

RNA-seq count data were collected from NCBI-GEO database generated on post-mortem brain tissues of controls (n = 24) and patients of BPD (n = 24) and SCZ (n = 24). Differentially expressed genes (DEGs) were identified using the consensus of DESeq2 and edgeR tools and used for downstream analysis. Weighted gene correlation networks were constructed to find non-preserved (NP) modules for SCZ, BPD and control conditions. Topological analysis and functional enrichment analysis were performed on NP modules to identify unique and overlapping expression signatures during SCZ and BPD conditions.

Results

We have identified four NP modules from the DEGs of BPD and SCZ. Eleven overlapping genes have been identified between SCZ and BPD networks, and they were found to be highly enriched in inflammatory responses. Among these eleven genes, TNIP2, TNFRSF1A and AC005840.1 had higher sum of connectivity exclusively in BPD network. In addition, we observed that top five genes of NP module from SCZ network were downregulated which may be a key factor for SCZ disorder.

Conclusions

Differential activation of the immune system components and pathways may drive the common and unique pathogenesis of the BPD and SCZ.

NURR1 and ERR1 Modulate the Expression of Genes of a DRD2 Coexpression Network Enriched for Schizophrenia Risk

Multiple schizophrenia (SCZ) risk loci may be involved in gene co-regulation mechanisms, and analysis of coexpressed gene networks may help to clarify SCZ molecular basis. We have previously identified a dopamine D2 receptor (DRD2) coexpression module enriched for SCZ risk genes and associated with cognitive and neuroimaging phenotypes of SCZ, as well as with response to treatment with antipsychotics. Here we aimed to identify regulatory factors modulating this coexpression module and their relevance to SCZ. We performed motif enrichment analysis to identify transcription factor (TF) binding sites in human promoters of genes coexpressed with DRD2. Then, we measured transcript levels of a group of these genes in primary mouse cortical neurons in basal conditions and upon overexpression and knockdown of predicted TFs. Finally, we analyzed expression levels of these TFs in dorsolateral prefrontal cortex (DLPFC) of SCZ patients. Our in silico analysis revealed enrichment for NURR1 and ERR1 binding sites. In neuronal cultures, the expression of genes either relevant to SCZ risk (Drd2, Gatad2a, Slc28a1, Cnr1) or indexing coexpression in our module (Btg4, Chit1, Osr1, Gpld1) was significantly modified by gain and loss of Nurr1 and Err1. Postmortem DLPFC expression data analysis showed decreased expression levels of NURR1 and ERR1 in patients with SCZ. For NURR1 such decreased expression is associated with treatment with antipsychotics. Our results show that NURR1 and ERR1 modulate the transcription of DRD2 coexpression partners and support the hypothesis that NURR1 is involved in the response to SCZ treatment.SIGNIFICANCE STATEMENT In the present study, we provide in silico and experimental evidence for a role of the TFs NURR1 and ERR1 in modulating the expression pattern of genes coexpressed with DRD2 in human DLPFC. Notably, genetic variations in these genes is associated with SCZ risk and behavioral and neuroimaging phenotypes of the disease, as well as with response to treatment. Furthermore, this study presents novel findings on a possible interplay between D2 receptor-mediated dopamine signaling involved in treatment with antipsychotics and the transcriptional regulation mechanisms exerted by NURR1. Our results suggest that coexpression and co-regulation mechanisms may help to explain some of the complex biology of genetic associations with SCZ.

Multiscale Neuroscience of Psychiatric Disorders

The human brain comprises a multiscale network with multiple levels of organization. Neurons with dendritic and axonal connections form the microscale fabric of brain circuitry, and macroscale brain regions and white matter connections form the infrastructure for system-level brain communication and information integration. In this review, we discuss the emerging trend of multiscale neuroscience, the multidisciplinary field that brings together data from these different levels of nervous system organization to form a better understanding of between-scale relationships of brain structure, function, and behavior in health and disease. We provide a broad overview of this developing field and discuss recent findings of exemplary multiscale neuroscience studies that illustrate the importance of studying cross-scale interactions among the genetic, molecular, cellular, and macroscale levels of brain circuitry and connectivity and behavior. We particularly consider a central, overarching goal of these multiscale neuroscience studies of human brain connectivity: to obtain insight into how disease-related alterations at one level of organization may underlie alterations observed at other scales of brain network organization in mental disorders. We conclude by discussing the current limitations, challenges, and future directions of the field.

Functional neuroimaging in treatment resistant schizophrenia: A systematic review

Despite the availability of several drugs, about 30% of patients with schizophrenia still fail to respond properly to a course of appropriate antipsychotic treatment. Functional neuroimaging studies have shown widespread patterns of altered activation and functional connectivity in treatment-resistant schizophrenia (TRS). The aim of the present study was to examine the available functional magnetic resonance imaging studies investigating TRS and to identify common patterns of altered brain function that could predict the lack of response to antipsychotic treatment in this disorder. Alterations of activation and functional connectivity in fronto-temporal, cortico-striatal, default mode network and salience networks, and of their interplay, were associated with TRS. Our findings support the notion that large-scale network alterations present in schizophrenia lie in a continuum within treatment response with the most severe dysfunction in TRS. Few studies with small sample size and without adequate control group limit the generalizability of current literature. Future controlled longitudinal studies are needed to identify neuroimaging biomarkers of pharmacotherapy response to inform individual treatment selection and facilitate early clinical response.