Gene Expression in Patient-Derived Neural Progenitors Implicates WNT5A Signaling in the Etiology of Schizophrenia.

BACKGROUND: Genome-wide association studies of schizophrenia have demonstrated that variations in noncoding regions are responsible for most of the common variation heritability of the disease. It is hypothesized that these risk variants alter gene expression. Therefore, studying alterations in gene expression in schizophrenia may provide a direct approach to understanding the etiology of the disease. In this study we use cultured neural progenitor cells derived from olfactory neuroepithelium (CNON cells) as a genetically unaltered cellular model to elucidate the neurodevelopmental aspects of schizophrenia. METHODS: We performed a gene expression study using RNA sequencing of CNON cells from 111 control subjects and 144 individuals with schizophrenia. Differentially expressed genes were identified with DESeq2 software, using covariates to correct for sex, age, library batches, and 1 surrogate variable component. RESULTS: A total of 80 genes were differentially expressed (false discovery rate < 10%), showing enrichment in cell migration, cell adhesion, developmental process, synapse assembly, cell proliferation, and related Gene Ontology categories. Cadherin and Wnt signaling pathways were positive in overrepresentation test, and, in addition, many genes were specifically involved in WNT5A signaling. The differentially expressed genes were modestly, but significantly, enriched in the genes overlapping single nucleotide polymorphisms with genome-wide significant association from the Psychiatric Genomics Consortium genome-wide association study of schizophrenia. We also found substantial overlap with genes associated with other psychiatric disorders or brain development, enrichment in the same Gene Ontology categories as genes with mutations de novo in schizophrenia, and studies of induced pluripotent stem cell-derived neural progenitor cells. CONCLUSIONS: CNON cells are a good model of the neurodevelopmental aspects of schizophrenia and can be used to elucidate the etiology of the disorder.

Deconvolution of transcriptional networks identifies TCF4 as a master regulator in schizophrenia.

Applying tissue-specific deconvolution of transcriptional networks to identify their master regulators (MRs) in neuropsychiatric disorders has been largely unexplored. Here, using two schizophrenia (SCZ) case-control RNA-seq datasets, one on postmortem dorsolateral prefrontal cortex (DLPFC) and another on cultured olfactory neuroepithelium, we deconvolved the transcriptional networks and identified TCF4 as a top candidate MR that may be dysregulated in SCZ. We validated TCF4 as a MR through enrichment analysis of TCF4-binding sites in induced pluripotent stem cell (hiPSC)-derived neurons and in neuroblastoma cells. We further validated the predicted TCF4 targets by knocking down TCF4 in hiPSC-derived neural progenitor cells (NPCs) and glutamatergic neurons (Glut_Ns). The perturbed TCF4 gene network in NPCs was more enriched for pathways involved in neuronal activity and SCZ-associated risk genes, compared to Glut_Ns. Our results suggest that TCF4 may serve as a MR of a gene network dysregulated in SCZ at early stages of neurodevelopment.

Androgenic regulation of sexually dimorphic expression of RNA binding motif protein 48 in the developing mouse cortex and hippocampus.

To further reveal the molecular mechanism underlying sexual differentiation of the mouse cerebral cortex and hippocampus, we reanalyzed our previous microarray study with Gene Ontology (GO) term enrichment and found that the GO term "RNA binding" was over-represented among the 89 sexually dimorphic candidate genes. Thus, we selected 16 autosomal genes annotated to the term RNA binding and profiled their mRNA expression in the developing male and female mouse cortex/hippocampus. During the first three weeks after birth, sex differences in mRNA levels of Khdrbs2, Nanos2, Rbm48, and Tdrd3 were observed in the mouse cortex/hippocampus. Of these genes, only the female-biased expression of Rbm48 in neonates was abolished by prenatal exposure to testosterone propionate (TP), while postnatal treatment of TP three weeks after birth increased Rbm48 and Tdrd3 mRNA levels in both sexes. Regardless of sex, the postnatal cortex/hippocampus also showed a marked increase in the content of androgen receptor (Ar) and estrogen receptor ß (Esr2), but a decrease in estrogen receptor α (Esr1) and aromatase (Cyp19a1), which might confer the different responses of Rbm48 to prenatal and postnatal TP. Our results suggest that androgen-regulated, sexually dimorphic Rbm48 expression might present a novel molecular mechanism by which perinatal androgens control development of sexual dimorphism in cortical and hippocampal structure and function.

Using 3D epigenomic maps of primary olfactory neuronal cells from living individuals to understand gene regulation.

As part of PsychENCODE, we developed a three-dimensional (3D) epigenomic map of primary cultured neuronal cells derived from olfactory neuroepithelium (CNON). We mapped topologically associating domains and high-resolution chromatin interactions using Hi-C and identified regulatory elements using chromatin immunoprecipitation and nucleosome positioning assays. Using epigenomic datasets from biopsies of 63 living individuals, we found that epigenetic marks at distal regulatory elements are more variable than marks at proximal regulatory elements. By integrating genotype and metadata, we identified enhancers that have different levels corresponding to differences in genetic variation, gender, smoking, and schizophrenia. Motif searches revealed that many CNON enhancers are bound by neuronal-related transcription factors. Last, we combined 3D epigenomic maps and gene expression profiles to predict enhancer-target gene interactions on a genome-wide scale. This study not only provides a framework for understanding individual epigenetic variation using a primary cell model system but also contributes valuable data resources for epigenomic studies of neuronal epithelium.

Analysis of Gene Expression Variance in Schizophrenia Using Structural Equation Modeling.

Schizophrenia (SCZ) is a psychiatric disorder of unknown etiology. There is evidence suggesting that aberrations in neurodevelopment are a significant attribute of schizophrenia pathogenesis and progression. To identify biologically relevant molecular abnormalities affecting neurodevelopment in SCZ we used cultured neural progenitor cells derived from olfactory neuroepithelium (CNON cells). Here, we tested the hypothesis that variance in gene expression differs between individuals from SCZ and control groups. In CNON cells, variance in gene expression was significantly higher in SCZ samples in comparison with control samples. Variance in gene expression was enriched in five molecular pathways: serine biosynthesis, PI3K-Akt, MAPK, neurotrophin and focal adhesion. More than 14% of variance in disease status was explained within the logistic regression model (C-value = 0.70) by predictors accounting for gene expression in 69 genes from these five pathways. Structural equation modeling (SEM) was applied to explore how the structure of these five pathways was altered between SCZ patients and controls. Four out of five pathways showed differences in the estimated relationships among genes: between KRAS and NF1, and KRAS and SOS1 in the MAPK pathway; between PSPH and SHMT2 in serine biosynthesis; between AKT3 and TSC2 in the PI3K-Akt signaling pathway; and between CRK and RAPGEF1 in the focal adhesion pathway. Our analysis provides evidence that variance in gene expression is an important characteristic of SCZ, and SEM is a promising method for uncovering altered relationships between specific genes thus suggesting affected gene regulation associated with the disease. We identified altered gene-gene interactions in pathways enriched for genes with increased variance in expression in SCZ. These pathways and loci were previously implicated in SCZ, providing further support for the hypothesis that gene expression variance plays important role in the etiology of SCZ.

Genetic variants specific to aging-related verbal memory: Insights from GWASs in a population-based cohort.

Verbal memory is typically studied using immediate recall (IR) and delayed recall (DR) scores, although DR is dependent on IR capability. Separating these components may be useful for deciphering the genetic variation in age-related memory abilities. This study was conducted to (a) construct individual trajectories in IR and independent aspects of delayed recall, or residualized-DR (rDR), across older adulthood; and (b) identify genetic markers that contribute to four estimated phenotypes: IR and rDR levels and changes after age 60. A cognitively intact sample (N = 20,650 with 125,164 observations) was drawn from the U.S. Health and Retirement Study, a nationally representative study of adults aged 50 and older. Mixed effects regression models were constructed using repeated measures from data collected every two years (1996-2012) to estimate level at age 60 and change in memory post-60 in IR and rDR. Genome-wide association scans (GWAS) were conducted in the genotypic subsample (N = 7,486) using ~1.2 million single nucleotide polymorphisms (SNPs). One SNP (rs2075650) in TOMM40 associated with rDR level at the genome-wide level (p = 5.0x10-08), an effect that replicated in an independent sample from the English Longitudinal Study on Ageing (N = 6,898 with 41,328 observations). Meta-analysis of rDR level confirmed the association (p = 5.0x10-11) and identified two others in TOMM40 (rs71352238 p = 1.0x10-10; rs157582 p = 7.0x10-09), and one in APOE (rs769449 p = 3.1 x10-12). Meta-analysis of IR change identified associations with three of the same SNPs in TOMM40 (rs157582 p = 8.3x10-10; rs71352238 p = 1.9x10-09) and APOE (rs769449 p = 2.2x10-08). Conditional analyses indicate GWAS signals on rDR level were driven by APOE, whereas signals on IR change were driven by TOMM40. Additionally, we found that TOMM40 had effects independent of APOE e4 on both phenotypes. Findings from this first U.S. population-based GWAS study conducted on both age-related immediate and delayed verbal memory merit continued examination in other samples and additional measures of verbal memory.

Effects of Prenatal Testosterone Exposure on Sexually Dimorphic Gene Expression in the Neonatal Mouse Cortex and Hippocampus.

OBJECTIVE: Using gene expression microarrays and reverse transcription with quantitative polymerase chain reaction (RT-qPCR), we have recently identified several novel genes that are differentially expressed in the neonatal male versus female mouse cortex/hippocampus (Armoskus et al.). Since perinatal testosterone (T) secreted by the developing testes masculinizes cortical and hippocampal structures and the behaviors regulated by these brain regions, we hypothesized that sexually dimorphic expression of specific selected genes in these areas might be regulated by T during early development. METHODS: To test our hypothesis, we treated timed pregnant female mice daily with vehicle or testosterone propionate (TP) starting on embryonic day 16 until the day of birth. The cortex/hippocampus was collected from vehicle- and TP-treated, male and female neonatal pups. Total RNA was extracted from these brain tissues, followed by RT-qPCR to measure relative mRNA levels of seven sex chromosome genes and three autosomal genes that have previously showed sex differences. RESULTS: The effect of prenatal TP was confirmed as it stimulated Dhcr24 expression in the neonatal mouse cortex/hippocampus and increased the anogenital distance in females. We found a significant effect of sex, but not TP, on expression of three Y-linked (Ddx3y, Eif2s3y, and Kdm5d), four X-linked (Eif2s3x, Kdm6a, Mid1, and Xist), and one autosomal (Klk8) genes in the neonatal mouse cortex/hippocampus. CONCLUSION: Although most of the selected genes are not directly regulated by prenatal T, their sexually dimorphic expression might play an important role in the control of sexually differentiated cognitive and social behaviors as well as in the etiology of sex-biased neurological disorders and mental illnesses.

Identification of sexually dimorphic genes in the neonatal mouse cortex and hippocampus.

The cerebral cortex and hippocampus are important for the control of cognitive functions and social behaviors, many of which are sexually dimorphic and tightly regulated by gonadal steroid hormones via activation of their respective nuclear receptors. As different levels of sex steroid hormones are present between the sexes during early development and their receptors act as transcription factors to regulate gene expression, we hypothesize that sexually dimorphic gene expression in the developing mouse cortex and hippocampus might result in sex differences in brain structures and neural circuits governing distinct behaviors between the sexes as adults. To test our hypothesis, we used gene expression microarrays to identify 90 candidate genes differentially expressed in the neonatal cortex/hippocampus between male and female mice, including 55 male-biased and 35 female-biased genes. Among these genes, sexually dimorphic expression of eight sex chromosome genes was confirmed by reverse transcription with quantitative PCR (RT-qPCR), including three located on the X chromosome (Xist, Eif2s3x, and Kdm6a), three on the Y chromosome (Ddx3y, Eif2s3y, and Kdm5d), and two in the pseudoautosomal region of the X and Y chromosomes (Erdr1 and Mid1). In addition, five autosomal genes (Cd151, Dab2, Klk8, Meg3, and Prkdc) were also validated for their sexually dimorphic expression in the neonatal mouse cortex/hippocampus. Gene Ontology annotation analysis suggests that many of these sexually dimorphic genes are involved in histone modifications, cell proliferation/death, androgen/estrogen signaling pathways, and synaptic organization, and these biological processes have been implicated in differential neural development, cognitive function, and neurological diseases between the sexes.