Search

Schizophrenia (SCZ) is a neuropsychiatric disorder with aberrant expression of multiple genes. However, identifying its exact causal genes remains a considerable challenge. The brain-specific transcription factor POU3F2 (POU domain, class 3, transcription factor 2) has been recognized as a risk factor for SCZ, but our understanding of its target genes and pathogenic mechanisms are still limited. Here we report that POU3F2 regulates 42 SCZ-related genes in knockdown and RNA-sequencing experiments of human neural progenitor cells (NPCs). Among those SCZ-related genes, TRIM8 (Tripartite motif containing 8) is located in SCZ-associated genetic locus and is aberrantly expressed in patients with SCZ. Luciferase reporter and electrophoretic mobility shift assays (EMSA) showed that POU3F2 induces TRIM8 expression by binding to the SCZ-associated SNP (single nucleotide polymorphism) rs5011218, which affects POU3F2-binding efficiency at the promoter region of TRIM8. We investigated the cellular functions of POU3F2 and TRIM8 as they co-regulate several pathways related to neural development and synaptic function. Knocking down either POU3F2 or TRIM8 promoted the proliferation of NPCs, inhibited their neuronal differentiation, and impaired the excitatory synaptic transmission of NPC-derived neurons. These results indicate that POU3F2 regulates TRIM8 expression through the SCZ-associated SNP rs5011218, and both genes may be involved in the etiology of SCZ by regulating neural development and synaptic function.

Genome-wide association studies (GWASs) have identified numerous single nucleotide polymorphisms (SNPs) associated with bipolar disorder (BD), but what the causal variants are and how they contribute to BD is largely unknown. In this study, we used FUMA, a GWAS annotation tool, to pinpoint potential causal variants and genes from the latest BD GWAS findings, and performed integrative analyses, including brain expression quantitative trait loci (eQTL), gene coexpression network, differential gene expression, protein-protein interaction, and brain intermediate phenotype association analysis to identify the functions of a prioritized gene and its connection to BD. Convergent lines of evidence prioritized protein-coding gene G Protein Nucleolar 3 (GNL3) as a BD risk gene, with integrative analyses revealing GNL3's roles in cell proliferation, neuronal functions, and brain phenotypes. We experimentally revealed that BD-related eQTL SNPs rs10865973, rs12635140, and rs4687644 regulate GNL3 expression using dual luciferase reporter assay and CRISPR interference experiment in human neural progenitor cells. We further identified that GNL3 knockdown and overexpression led to aberrant neuronal proliferation and differentiation, using two-dimensional human neural cell cultures and three-dimensional forebrain organoid model. This study gathers evidence that BD-related genetic variants regulate GNL3 expression which subsequently affects neuronal proliferation and differentiation.

A number of studies indicate that rare copy number variations (CNVs) contribute to the risk of schizophrenia (SCZ). Most of these studies have focused on protein-coding genes residing in the CNVs. Here, we investigated long noncoding RNAs (lncRNAs) within 10 SCZ risk-associated CNV deletion regions (CNV-lncRNAs) and examined their potential contribution to SCZ risk. We used RNA sequencing transcriptome data derived from postmortem brain tissue from control individuals without psychiatric disease as part of the PsychENCODE BrainGVEX and Developmental Capstone projects. We carried out weighted gene coexpression network analysis to identify protein-coding genes coexpressed with CNV-lncRNAs in the human brain. We identified one neuronal function-related coexpression module shared by both datasets. This module contained a lncRNA called DGCR5 within the 22q11.2 CNV region, which was identified as a hub gene. Protein-coding genes associated with SCZ genome-wide association study signals, de novo mutations, or differential expression were also contained in this neuronal module. Using DGCR5 knockdown and overexpression experiments in human neural progenitor cells derived from human induced pluripotent stem cells, we identified a potential role for DGCR5 in regulating certain SCZ-related genes.

Schizophrenia and bipolar disorder are complex psychiatric diseases with risks contributed by multiple genes. Dysregulation of gene expression has been implicated in these disorders, but little is known about such dysregulation in the human brain. We analyzed three transcriptome datasets from 394 postmortem brain tissue samples from patients with schizophrenia or bipolar disorder or from healthy control individuals without a known history of psychiatric disease. We built genome-wide coexpression networks that included microRNAs (miRNAs). We identified a coexpression network module that was differentially expressed in the brain tissue from patients compared to healthy control individuals. This module contained genes that were principally involved in glial and neural cell genesis and glial cell differentiation, and included schizophrenia risk genes carrying rare variants. This module included five miRNAs and 545 mRNAs, with six transcription factors serving as hub genes in this module. We found that the most connected transcription factor gene POU3F2, also identified on a genome-wide association study for bipolar disorder, could regulate the miRNA hsa-miR-320e and other putative target mRNAs. These regulatory relationships were replicated using PsychENCODE/BrainGVEX datasets and validated by knockdown and overexpression experiments in SH-SY5Y cells and human neural progenitor cells in vitro. Thus, we identified a brain gene expression module that was enriched for rare coding variants in genes associated with schizophrenia and that contained the putative bipolar disorder risk gene POU3F2 The transcription factor POU3F2 may be a key regulator of gene expression in this disease-associated gene coexpression module.

Schizophrenia is a complex psychiatric disorder with strong genetic background. Studies demonstrated that rare Copy Number Variations (CNVs) contribute to the risk of schizophrenia. More attention has been put to protein-coding genes residing in those CNVs. However, the exact mechanisms still remain unclear. In this study, we explored the potential roles of long non-coding RNAs (lncRNAs) inside CNV deletions (CNV-lncRNAs) in the risk of schizophrenia. We retrieved lncRNAs mapped to the CNV deletion regions known to increase risks of developing schizophrenia. Seven such regions were repeatedly reported by large case-control studies, including 1q21.1, 3q29, 15q11.2, 15q13.3, 17p12, 17q12, and 22q11.2. We carried out weighted gene co-expression network analysis (WGCNA) using RNA-seq data from Genotype-Tissue Expression (GTEx) and BrainSpan projects to look for co-expression modules that harbor CNV-lncRNAs. We identified one male reproduction-related module in male individuals and one neuronal functions-related module both in male and female individuals. In addition, the neuronal functions module comprised of protein-coding genes that were involved in several ion channel activities including calcium and potassium channel activities, which were known to be related to schizophrenia. Pathway analysis of these two modules further suggested that CNV-lncRNAs were involved in the olfactory transduction, neuroactive ligand-receptor and calcium signaling pathways. The CNV-lncRNA co-expression patterns were preserved through the different development and aging time points, and across different tissue types of human body. Our results suggested that lncRNAs inside those rare CNVs might play significant temporal and spatial roles in regulating other protein-coding genes and subsequently contribute to schizophrenia risk.