Human forebrain organoid-based multi-omics analyses of PCCB as a schizophrenia associated gene linked to GABAergic pathways.

Identifying genes whose expression is associated with schizophrenia (SCZ) risk by transcriptome-wide association studies (TWAS) facilitates downstream experimental studies. Here, we integrated multiple published datasets of TWAS, gene coexpression, and differential gene expression analysis to prioritize SCZ candidate genes for functional study. Convergent evidence prioritized Propionyl-CoA Carboxylase Subunit Beta (PCCB), a nuclear-encoded mitochondrial gene, as an SCZ risk gene. However, the PCCB's contribution to SCZ risk has not been investigated before. Using dual luciferase reporter assay, we identified that SCZ-associated SNPs rs6791142 and rs35874192, two eQTL SNPs for PCCB, showed differential allelic effects on transcriptional activities. PCCB knockdown in human forebrain organoids (hFOs) followed by RNA sequencing analysis revealed dysregulation of genes enriched with multiple neuronal functions including gamma-aminobutyric acid (GABA)-ergic synapse. The metabolomic and mitochondrial function analyses confirmed the decreased GABA levels resulted from inhibited tricarboxylic acid cycle in PCCB knockdown hFOs. Multielectrode array recording analysis showed that PCCB knockdown in hFOs resulted into SCZ-related phenotypes including hyper-neuroactivities and decreased synchronization of neural network. In summary, this study utilized hFOs-based multi-omics analyses and revealed that PCCB downregulation may contribute to SCZ risk through regulating GABAergic pathways, highlighting the mitochondrial function in SCZ.

Human forebrain organoids-based multi-omics analyses reveal PCCB's regulation on GABAergic system contributing to schizophrenia.

Identifying genes whose expression is associated with schizophrenia (SCZ) risk by transcriptome-wide association studies (TWAS) facilitates downstream experimental studies. Here, we integrated multiple published datasets of TWAS (including FUSION, PrediXcan, summary-data-based Mendelian randomization (SMR), joint-tissue imputation approach with Mendelian randomization (MR-JTI)), gene coexpression, and differential gene expression analysis to prioritize SCZ candidate genes for functional study. Convergent evidence prioritized Propionyl-CoA Carboxylase Subunit Beta ( PCCB ), a nuclear-encoded mitochondrial gene, as an SCZ risk gene. However, the PCCB ’s contribution to SCZ risk has not been investigated before. Using dual luciferase reporter assay, we identified that SCZ-associated SNP rs35874192, an eQTL SNP for PCCB , showed differential allelic effects on transcriptional activities. PCCB knockdown in human forebrain organoids (hFOs) followed by RNA-seq revealed dysregulation of genes enriched with multiple neuronal functions including gamma-aminobutyric acid (GABA)-ergic synapse, as well as genes dysregulated in postmortem brains of SCZ patients or in cerebral organoids derived from SCZ patients. The metabolomic and mitochondrial function analyses confirmed the deceased GABA levels resulted from reduced tricarboxylic acid cycle in PCCB knockdown hFOs. Multielectrode array recording analysis showed that PCCB knockdown in hFOs resulted into SCZ-related phenotypes including hyper-neuroactivities and decreased synchronization of neural network. In summary, this study utilized hFOs-based multi-omics data and revealed that PCCB downregulation may contribute to SCZ risk through regulating GABAergic system, highlighting the mitochondrial function in SCZ.

Microglia-containing human brain organoids for the study of brain development and pathology.

Microglia are resident immune cells in the central nervous system, playing critical roles in brain development and homeostasis. Increasing evidence has implicated microglia dysfunction in the pathogenesis of various brain disorders ranging from psychiatric disorders to neurodegenerative diseases. Using a human cell-based model to illuminate the functional mechanisms of microglia will promote pathological studies and drug development. The recently developed microglia-containing human brain organoids (MC-HBOs), in-vitro three-dimensional cell cultures that recapitulate key features of the human brain, have provided a new avenue to model brain development and pathology. However, MC-HBOs generated from different methods differ in the origin, proportion, and fidelity of microglia within the organoids, and may have produced inconsistent results. To help researchers to develop a robust and reproducible model that recapitulates in-vivo signatures of human microglia to study brain development and pathology, this review summarized the current methods used to generate MC-HBOs and provided opinions on the use of MC-HBOs for disease modeling and functional studies.

Expression of expanded GGC repeats within NOTCH2NLC causes behavioral deficits and neurodegeneration in a mouse model of neuronal intranuclear inclusion disease.

GGC repeat expansions within NOTCH2NLC have been identified as the genetic cause of neuronal intranuclear inclusion disease (NIID). To understand the molecular pathogenesis of NIID, here, we established both a transgenic mouse model and a human neural progenitor cells (hNPCs) model. Expression of the NOTCH2NLC with expanded GGC repeats produced widespread intranuclear and perinuclear polyglycine (polyG), polyalanine (polyA), and polyarginine (polyR) inclusions, leading to behavioral deficits and severe neurodegeneration, which faithfully mimicked the clinical and pathological features associated with NIID. Furthermore, conserved alternative splicing events were identified between the NIID mouse and hNPC models, among which was the enrichment of the binding motifs of hnRNPM, an RNA binding protein known as alternative splicing regulator. Expanded NOTCH2NLC-polyG and NOTCH2NLC-polyA could interact with and sequester hnRNPM, while overexpression of hnRNPM could ameliorate the cellular toxicity. These results together suggested that dysfunction of hnRNPM could play an important role in the molecular pathogenesis of NIID.

Profiling the Genome-Wide Landscape of Short Tandem Repeats by Long-Read Sequencing.

Background: Short tandem repeats (STRs) are highly variable elements that play a pivotal role in multiple genetic diseases and the regulation of gene expression. Long-read sequencing (LRS) offers a potential solution to genome-wide STR analysis. However, characterizing STRs in human genomes using LRS on a large population scale has not been reported. Methods: We conducted the large LRS-based STR analysis in 193 unrelated samples of the Chinese population and performed genome-wide profiling of STR variation in the human genome. The repeat dynamic index (RDI) was introduced to evaluate the variability of STR. We sourced the expression data from the Genotype-Tissue Expression to explore the tissue specificity of highly variable STRs related genes across tissues. Enrichment analyses were also conducted to identify potential functional roles of the high variable STRs. Results: This study reports the large-scale analysis of human STR variation by LRS and offers a reference STR database based on the LRS dataset. We found that the disease-associated STRs (dSTRs) and STRs associated with the expression of nearby genes (eSTRs) were highly variable in the general population. Moreover, tissue-specific expression analysis showed that those highly variable STRs related genes presented the highest expression level in brain tissues, and enrichment pathways analysis found those STRs are involved in synaptic function-related pathways. Conclusion: Our study profiled the genome-wide landscape of STR using LRS and highlighted the highly variable STRs in the human genome, which provide a valuable resource for studying the role of STRs in human disease and complex traits.

Human forebrain organoids reveal connections between valproic acid exposure and autism risk.

Valproic acid (VPA) exposure as an environmental factor that confers risk of autism spectrum disorder (ASD), its functional mechanisms in the human brain remain unclear since relevant studies are currently restricted to two-dimensional cell cultures and animal models. To identify mechanisms by which VPA contribute to ASD risk in human, here we used human forebrain organoids (hFOs), in vitro derived three-dimensional cell cultures that recapitulate key human brain developmental features. We identified that VPA exposure in hFOs affected the expression of genes enriched in neural development, synaptic transmission, oxytocin signaling, calcium, and potassium signaling pathways, which have been implicated in ASD. Genes (e.g., CAMK4, CLCN4, DPP10, GABRB3, KCNB1, PRKCB, SCN1A, and SLC24A2) that affected by VPA were significantly overlapped with those dysregulated in brains or organoids derived from ASD patients, and known ASD risk genes, as well as genes in ASD risk-associated gene coexpression modules. Single-cell RNA sequencing analysis showed that VPA exposure affected the expression of genes in choroid plexus, excitatory neuron, immature neuron, and medial ganglionic eminence cells annotated in hFOs. Microelectrode array further identified that VPA exposure in hFOs disrupted synaptic transmission. Taken together, this study connects VPA exposure to ASD pathogenesis using hFOs, which is valuable for illuminating the etiology of ASD and screening for potential therapeutic targets.

Transcription factor POU3F2 regulates TRIM8 expression contributing to cellular functions implicated in schizophrenia.

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.

De novo mutation of cancer-related genes associates with particular neurodevelopmental disorders.

Epidemiological studies have shown an increased prevalence of cancer in some patients with neurodevelopmental disorder (NDD); however, the genetic mechanisms regarding how cancer-related genes (CRGs) contribute to NDD remain unclear. We performed bioinformatic analyses on 219 CRGs from OMIM and de novo mutations (DNMs) from 16,498 patients with different NDDs and 3391 controls. Our results showed that autism spectrum disorder, undiagnosed neurodevelopmental disorder, congenital heart disease and intellectual disability, but not epileptic encephalopathy and schizophrenia, harboured significantly more putative functional DNMs in CRGs, compared with controls, providing genetic evidence supporting previous epidemiological surveys. We further detected 26 CRGs with recurrent putative functional DNMs that showed high expression in the human brain during the prenatal stage and in non-brain organs in adults. The proteins coded by the 26 CRGs and known NDD candidate genes formed a functional network that is involved in brain development and tumorigenesis. Overall, we proposed 39 cancer-targeting drugs that could be investigated for treating patients with NDD, which would be potentially cost-effective. In conclusion, DNMs contribute to specific NDDs and there may be a shared genetic basis between NDDs and cancer, highlighting the importance of considering cancer-targeting drugs with potential curative effects in patients with NDDs. KEY MESSAGES: • The contribution of DNMs in NDD is consistent with epidemiological surveys. • We highlighted 26 CRGs, including nine genes with more than five functional DNMs. • Specific expression patterns underlie the genetic mechanism of CRGs in NDD. • Specific functional networks underlie the genetic mechanism of CRGs in NDD. • The shared genetic aetiology suggests potential mutual treatment strategies.

Integrative analyses prioritize GNL3 as a risk gene for bipolar disorder.

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.

Brain Banks Spur New Frontiers in Neuropsychiatric Research and Strategies for Analysis and Validation.

Neuropsychiatric disorders affect hundreds of millions of patients and families worldwide. To decode the molecular framework of these diseases, many studies use human postmortem brain samples. These studies reveal brain-specific genetic and epigenetic patterns via high-throughput sequencing technologies. Identifying best practices for the collection of postmortem brain samples, analyzing such large amounts of sequencing data, and interpreting these results are critical to advance neuropsychiatry. We provide an overview of human brain banks worldwide, including progress in China, highlighting some well-known projects using human postmortem brain samples to understand molecular regulation in both normal brains and those with neuropsychiatric disorders. Finally, we discuss future research strategies, as well as state-of-the-art statistical and experimental methods that are drawn upon brain bank resources to improve our understanding of the agents of neuropsychiatric disorders.

BrainEXP: a database featuring with spatiotemporal expression variations and co-expression organizations in human brains.

Summary: Gene expression changes over the lifespan and varies among different tissues or cell types. Gene co-expression also changes by sex, age, different tissues or cell types. However, gene expression under the normal state and gene co-expression in the human brain has not been fully defined and quantified. Here we present a database named Brain EXPression Database (BrainEXP) which provides spatiotemporal expression of individual genes and co-expression in normal human brains. BrainEXP consists of 4567 samples from 2863 healthy individuals gathered from existing public databases and our own data, in either microarray or RNA-Seq library types. We mainly provide two analysis results based on the large dataset: (i) basic gene expression across specific brain regions, age ranges and sexes; (ii) co-expression analysis from different platforms. Availability and implementation: http://www.brainexp.org/. Supplementary information: Supplementary data are available at Bioinformatics online.

The transcription factor POU3F2 regulates a gene coexpression network in brain tissue from patients with psychiatric disorders.

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.

The DGCR5 long noncoding RNA may regulate expression of several schizophrenia-related genes.

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.