Disruptive mutation in the CHD8 gene is one of the top genetic risk factors in autism spectrum disorders (ASDs). Previous analyses of genome-wide CHD8 occupancy and reduced expression of CHD8 by shRNA knockdown in committed neural cells showed that CHD8 regulates multiple cell processes critical for neural functions, and its targets are enriched with ASD-associated genes.

ZNF804A (Zinc Finger Protein 804A) has been identified as a candidate gene for schizophrenia (SZ), autism spectrum disorders (ASD), and bipolar disorder (BD) in replicated genome wide association studies (GWAS) and by copy number variation (CNV) analysis. Although its function has not been well-characterized, ZNF804A contains a C2H2-type zinc-finger domain, suggesting that it has DNA binding properties, and consequently, a role in regulating gene expression. To further explore the role of ZNF804A on gene expression and its downstream targets, we used a gene knockdown (KD) approach to reduce its expression in neural progenitor cells (NPCs) derived from induced pluripotent stem cells (iPSCs). KD was accomplished by RNA interference (RNAi) using lentiviral particles containing shRNAs that target ZNF804A mRNA. Stable transduced NPC lines were generated after puromycin selection. A control cell line expressing a random (scrambled) shRNA was also generated. Neuronal differentiation was induced, RNA was harvested after 14 days and transcriptome analysis was carried out using RNA-seq. 1815 genes were found to be differentially expressed at a nominally significant level (p<0.05); 809 decreased in expression in the KD samples, while 1106 increased. Of these, 370 achieved genome wide significance (FDR<0.05); 125 were lower in the KD samples, 245 were higher. Pathway analysis showed that genes involved in interferon-signaling were enriched among those that were down-regulated in the KD samples. Correspondingly, ZNF804A KD was found to affect interferon-alpha 2 (IFNA2)-mediated gene expression. The findings suggest that ZNF804A may affect a differentiating neuron's response to inflammatory cytokines, which is consistent with models of SZ and ASD that support a role for infectious disease, and/or autoimmunity in a subgroup of patients.

Genome-wide expression analysis using next generation sequencing (RNA-Seq) provides an opportunity for in-depth molecular profiling of fundamental biological processes, such as cellular differentiation and malignant transformation. Differentiating human neurons derived from induced pluripotent stem cells (iPSCs) provide an ideal system for RNA-Seq since defective neurogenesis caused by abnormalities in transcription factors, DNA methylation, and chromatin modifiers lie at the heart of some neuropsychiatric disorders. As a preliminary step towards applying next generation sequencing using neurons derived from patient-specific iPSCs, we have carried out an RNA-Seq analysis on control human neurons. Dramatic changes in the expression of coding genes, long non-coding RNAs (lncRNAs), pseudogenes, and splice isoforms were seen during the transition from pluripotent stem cells to early differentiating neurons. A number of genes that undergo radical changes in expression during this transition include candidates for schizophrenia (SZ), bipolar disorder (BD) and autism spectrum disorders (ASD) that function as transcription factors and chromatin modifiers, such as POU3F2 and ZNF804A, and genes coding for cell adhesion proteins implicated in these conditions including NRXN1 and NLGN1. In addition, a number of novel lncRNAs were found to undergo dramatic changes in expression, one of which is HOTAIRM1, a regulator of several HOXA genes during myelopoiesis. The increase we observed in differentiating neurons suggests a role in neurogenesis as well. Finally, several lncRNAs that map near SNPs associated with SZ in genome wide association studies also increase during neuronal differentiation, suggesting that these novel transcripts may be abnormally regulated in a subgroup of patients.

We are using induced pluripotent stem cell (iPSC) technology to study neuropsychiatric disorders associated with 22q11.2 microdeletions (del), the most common known schizophrenia (SZ)-associated genetic factor. Several genes in the region have been implicated; a promising candidate is DGCR8, which codes for a protein involved in microRNA (miRNA) biogenesis. We carried out miRNA expression profiling (miRNA-seq) on neurons generated from iPSCs derived from controls and SZ patients with 22q11.2 del. Using thresholds of p<0.01 for nominal significance and 1.5-fold differences in expression, 45 differentially expressed miRNAs were detected (13 lower in SZ and 32 higher). Of these, 6 were significantly down-regulated in patients after correcting for genome wide significance (FDR<0.05), including 4 miRNAs that map to the 22q11.2 del region. In addition, a nominally significant increase in the expression of several miRNAs was found in the 22q11.2 neurons that were previously found to be differentially expressed in autopsy samples and peripheral blood in SZ and autism spectrum disorders (e.g., miR-34, miR-4449, miR-146b-3p, and miR-23a-5p). Pathway and function analysis of predicted mRNA targets of the differentially expressed miRNAs showed enrichment for genes involved in neurological disease and psychological disorders for both up and down regulated miRNAs. Our findings suggest that: i. neurons with 22q11.2 del recapitulate the miRNA expression patterns expected of 22q11.2 haploinsufficiency, ii. differentially expressed miRNAs previously identified using autopsy samples and peripheral cells, both of which have significant methodological problems, are indeed disrupted in neuropsychiatric disorders and likely have an underlying genetic basis.

Induced pluripotent stem cell (iPSC) technology is providing an opportunity to study neuropsychiatric disorders through the capacity to grow patient-specific neurons in vitro. Skin fibroblasts obtained by biopsy have been the most reliable source of cells for reprogramming. However, using other somatic cells obtained by less invasive means would be ideal, especially in children with autism spectrum disorders (ASD) and other neurodevelopmental conditions. In addition to fibroblasts, iPSCs have been developed from cord blood, lymphocytes, hair keratinocytes, and dental pulp from deciduous teeth. Of these, dental pulp would be a good source for neurodevelopmental disorders in children because obtaining material is non-invasive. We investigated its suitability for disease modeling by carrying out gene expression profiling, using RNA-seq, on differentiated neurons derived from iPSCs made from dental pulp extracted from deciduous teeth (T-iPSCs) and fibroblasts (F-iPSCs). This is the first RNA-seq analysis comparing gene expression profiles in neurons derived from iPSCs made from different somatic cells. For the most part, gene expression profiles were quite similar with only 329 genes showing differential expression at a nominally significant p-value (p<0.05), of which 63 remained significant after correcting for genome-wide analysis (FDR <0.05). The most striking difference was the lower level of expression detected for numerous members of the all four HOX gene families in neurons derived from T-iPSCs. In addition, an increased level of expression was seen for several transcription factors expressed in the developing forebrain (FOXP2, OTX1, and LHX2, for example). Overall, pathway analysis revealed that differentially expressed genes that showed higher levels of expression in neurons derived from T-iPSCs were enriched for genes implicated in schizophrenia (SZ). The findings suggest that neurons derived from T-iPSCs are suitable for disease-modeling neuropsychiatric disorder and may have some advantages over those derived from F-iPSCs.

Recent studies have shown that the transcriptional functions of REST are much broader than repressing neuronal genes in non-neuronal systems. Whether REST occupies similar chromatin regions in different cell types and how it interacts with other transcriptional regulators to execute its functions in a context-dependent manner has not been adequately investigated. We have applied ChIP-seq analysis to identify the REST cistrome in human CD4+ T cells and compared it with published data from 15 other cell types. We found that REST cistromes were distinct among cell types, with REST binding to several tumor suppressors specifically in cancer cells, whereas 7% of the REST peaks in non-neuronal cells were ubiquitously called and <25% were identified for ≥ 5 cell types. Nevertheless, using a quantitative metric directly comparing raw ChIP-seq signals, we found the majority (∼80%) was shared by ≥ 2 cell types. Integration with RNA-seq data showed that REST binding was generally correlated with low gene expression. Close examination revealed that multiple contexts were correlated with reduced expression of REST targets, e.g., the presence of a cognate RE1 motif and cellular specificity of REST binding. These contexts were shown to play a role in differential corepressor recruitment. Furthermore, transcriptional outcome was highly influenced by REST cofactors, e.g., SIN3 and EZH2 co-occupancy marked higher and lower expression of REST targets, respectively. Unexpectedly, the REST cistrome in differentiated neurons exhibited unique features not observed in non-neuronal cells, e.g., the lack of RE1 motifs and an association with active gene expression. Finally, our analysis demonstrated how REST could differentially regulate a transcription network constituted of miRNAs, REST complex and neuronal factors. Overall, our findings of contexts playing critical roles in REST occupancy and regulatory outcome provide insights into the molecular interactions underlying REST's diverse functions, and point to novel roles of REST in differentiated neurons.

Thousands of pseudogenes exist in the human genome and many are transcribed, but their functional potential remains elusive and understudied. To explore these issues systematically, we first developed a computational pipeline to identify transcribed pseudogenes from RNA-Seq data. Applying the pipeline to datasets from 16 distinct normal human tissues identified ∼ 3,000 pseudogenes that could produce non-coding RNAs in a manner of low abundance but high tissue specificity under normal physiological conditions. Cross-tissue comparison revealed that the transcriptional profiles of pseudogenes and their parent genes showed mostly positive correlations, suggesting that pseudogene transcription could have a positive effect on the expression of their parent genes, perhaps by functioning as competing endogenous RNAs (ceRNAs), as previously suggested and demonstrated with the PTEN pseudogene, PTENP1. Our analysis of the ENCODE project data also found many transcriptionally active pseudogenes in the GM12878 and K562 cell lines; moreover, it showed that many human pseudogenes produced small RNAs (sRNAs) and some pseudogene-derived sRNAs, especially those from antisense strands, exhibited evidence of interfering with gene expression. Further integrated analysis of transcriptomics and epigenomics data, however, demonstrated that trimethylation of histone 3 at lysine 9 (H3K9me3), a posttranslational modification typically associated with gene repression and heterochromatin, was enriched at many transcribed pseudogenes in a transcription-level dependent manner in the two cell lines. The H3K9me3 enrichment was more prominent in pseudogenes that produced sRNAs at pseudogene loci and their adjacent regions, an observation further supported by the co-enrichment of SETDB1 (a H3K9 methyltransferase), suggesting that pseudogene sRNAs may have a role in regional chromatin repression. Taken together, our comprehensive and systematic characterization of pseudogene transcription uncovers a complex picture of how pseudogene ncRNAs could influence gene and pseudogene expression, at both epigenetic and post-transcriptional levels.

Individuals with 22q11.2 Deletion Syndrome (22q11.2 DS) are a specific high-risk group for developing schizophrenia (SZ), schizoaffective disorder (SAD) and autism spectrum disorders (ASD). Several genes in the deleted region have been implicated in the development of SZ, e.g., PRODH and DGCR8. However, the mechanistic connection between these genes and the neuropsychiatric phenotype remains unclear. To elucidate the molecular consequences of 22q11.2 deletion in early neural development, we carried out RNA-seq analysis to investigate gene expression in early differentiating human neurons derived from induced pluripotent stem cells (iPSCs) of 22q11.2 DS SZ and SAD patients.

Rett syndrome (RTT) is a severe, neurodevelopmental disorder primarily affecting girls, characterized by progressive loss of cognitive, social, and motor skills after a relatively brief period of typical development. It is usually due to de novo loss of function mutations in the X-linked gene, MeCP2, which codes for the gene expression and chromatin regulator, methyl-CpG binding protein 2. Although the behavioral phenotype appears to be primarily due to neuronal Mecp2 deficiency in mice, other cell types, including astrocytes and oligodendrocytes, also appear to contribute to some aspects of the RTT phenotype. In addition, microglia may also play a role. However, the effect of Mecp2 deficiency in microglia on RTT pathogenesis is controversial.

Schizophrenia (SZ) and autism spectrum disorders (ASD) are highly heritable neuropsychiatric disorders, although environmental factors, such as maternal immune activation (MIA), play a role as well. Cytokines mediate the effects of MIA on neurogenesis and behavior in animal models. However, MIA stimulators can also induce a febrile reaction, which could have independent effects on neurogenesis through heat shock (HS)-regulated cellular stress pathways. However, this has not been well-studied. To help understand the role of fever in MIA, we used a recently described model of human brain development in which induced pluripotent stem cells (iPSCs) differentiate into 3-dimensional neuronal aggregates that resemble a first trimester telencephalon. RNA-seq was carried out on aggregates that were heat shocked at 39°C for 24 hours, along with their control partners maintained at 37°C. 186 genes showed significant differences in expression following HS (p<0.05), including known HS-inducible genes, as expected, as well as those coding for NGFR and a number of SZ and ASD candidates, including SMARCA2, DPP10, ARNT2, AHI1 and ZNF804A. The degree to which the expression of these genes decrease or increase during HS is similar to that found in copy loss and copy gain copy number variants (CNVs), although the effects of HS are likely to be transient. The dramatic effect on the expression of some SZ and ASD genes places HS, and perhaps other cellular stressors, into a common conceptual framework with disease-causing genetic variants. The findings also suggest that some candidate genes that are assumed to have a relatively limited impact on SZ and ASD pathogenesis based on a small number of positive genetic findings, such as SMARCA2 and ARNT2, may in fact have a much more substantial role in these disorders - as targets of common environmental stressors.

Autism spectrum disorder (ASD) is highly heritable but genetically heterogeneous. The affected neural circuits and cell types remain unclear and may vary at different developmental stages. By analyzing multiple sets of human single cell transcriptome profiles, we found that ASD candidates showed relatively enriched gene expression in neurons, especially in inhibitory neurons. ASD candidates were also more likely to be the hubs of the co-expression gene module that is highly expressed in inhibitory neurons, a feature not detected for excitatory neurons. In addition, we found that upregulated genes in multiple ASD cortex samples were enriched with genes highly expressed in inhibitory neurons, suggesting a potential increase of inhibitory neurons and an imbalance in the ratio between excitatory and inhibitory neurons in ASD brains. Furthermore, the downstream targets of several ASD candidates, such as CHD8, EHMT1 and SATB2, also displayed enriched expression in inhibitory neurons. Taken together, our analyses of single cell transcriptomic data suggest that inhibitory neurons may be a major neuron subtype affected by the disruption of ASD gene networks, providing single cell functional evidence to support the excitatory/inhibitory (E/I) imbalance hypothesis.

Lowe syndrome (LS) is a rare genetic disorder caused by loss of function mutations in the X-linked gene, OCRL, which codes for inositol polyphosphate 5-phosphatase. LS is characterized by the triad of congenital cataracts, neurodevelopmental impairment (primarily intellectual and developmental disabilities [IDD]), and renal proximal tubular dysfunction. Studies carried out over the years have shown that hypomorphic mutations in OCRL adversely affect endosome recycling and actin polymerization in kidney cells and patient-derived fibroblasts. The renal problem has been traced to an impaired recycling of megalin, a multi-ligand receptor that plays a key role in the reuptake of lipoproteins, amino acids, vitamin-binding proteins, and hormones. However, the neurodevelopmental aspects of the disorder have been difficult to study because the mouse knockout (KO) model does not display LS-related phenotypes. Fortunately, the discovery of induced pluripotent stem (iPS) cells has provided an opportunity to grow patient-specific neurons, which can be used to model neurodevelopmental disorders in vitro, as demonstrated in the many studies that have been published in the past few years in autism spectrum disorders (ASD), schizophrenia (SZ), bipolar disorder (BD), and IDD.

Deletions encompassing the BP1-2 region at 15q11.2 increase schizophrenia and epilepsy risk, but only some carriers have either disorder. To investigate the role of CYFIP1, a gene within the region, we performed knockdown experiments in human neural progenitors derived from donors with 2 copies of each gene at the BP1-2 locus. RNA-seq and cellular assays determined that knockdown of CYFIP1 compromised cytoskeletal remodeling. FMRP targets and postsynaptic density genes, each implicated in schizophrenia, were significantly overrepresented among differentially expressed genes (DEGs). Schizophrenia and/or epilepsy genes, but not those associated with randomly selected disorders, were likewise significantly overrepresented. Mirroring the variable expressivity seen in deletion carriers, marked between-line differences were observed for dysregulation of disease genes. Finally, a subset of DEGs showed a striking similarity to known epilepsy genes and represents novel disease candidates. Results support a role for CYFIP1 in disease and demonstrate that disease-related biological signatures are apparent prior to neuronal differentiation.

Molecular genetic studies are typically performed on homogenized biological samples, resulting in contamination from non-neuronal cells. To improve expression profiling of neurons we combined patch recordings with single-cell PCR. Two iPSC lines (healthy subject and 22q11.2 deletion) were differentiated into neurons. Patch electrode recordings were performed on 229 human cells from Day-13 to Day-88, followed by capture and single-cell PCR for 13 genes: ACTB, HPRT, vGLUT1, βTUBIII, COMT, DISC1, GAD1, PAX6, DTNBP1, ERBB4, FOXP1, FOXP2, and GIRK2. Neurons derived from both iPSC lines expressed βTUBIII, fired action potentials, and experienced spontaneous depolarizations (UP states) ~2 weeks before vGLUT1, GAD1 and GIRK2 appeared. Multisite calcium imaging revealed that these UP states were not synchronized among hESC-H9-derived neurons. The expression of FOXP1, FOXP2 and vGLUT1 was lost after 50 days in culture, in contrast to other continuously expressed genes. When gene expression was combined with electrophysiology, two subsets of genes were apparent; those irrelevant to spontaneous depolarizations (including vGLUT1, GIRK2, FOXP2 and DISC1) and those associated with spontaneous depolarizations (GAD1 and ERBB4). The results demonstrate that in the earliest stages of neuron development, it is useful to combine genetic analysis with physiological characterizations, on a cell-to-cell basis.

Jansen de Vries Syndrome (JdVS) is a rare neurodevelopmental disorder (NDD) caused by gain-of-function (GOF) truncating mutations in PPM1D exons 5 or 6. PPM1D is a serine/threonine phosphatase that plays an important role in the DNA damage response (DDR) by negatively regulating TP53 (P53). JdVS-associated mutations lead to the formation of a truncated PPM1D protein that retains catalytic activity and has a GOF effect because of reduced degradation. Somatic PPM1D exons 5 and 6 truncating mutations are well-established factors in a number of cancers, due to excessive dephosphorylation and reduced function of P53 and other substrates involved in DDR. Children with JdVS have a variety of neurodevelopmental, psychiatric, and physical problems. In addition, a small fraction has acute neuropsychiatric decompensation apparently triggered by infection or severe non-infectious environmental stress factors.

CHD8 (chromodomain helicase DNA-binding protein 8), which codes for a member of the CHD family of ATP-dependent chromatin-remodeling factors, is one of the most commonly mutated genes in autism spectrum disorders (ASD) identified in exome-sequencing studies. Loss of function mutations in the gene have also been found in schizophrenia (SZ) and intellectual disabilities and influence cancer cell proliferation. We previously reported an RNA-seq analysis carried out on neural progenitor cells (NPCs) and monolayer neurons derived from induced pluripotent stem (iPS) cells that were heterozygous for CHD8 knockout (KO) alleles generated using CRISPR-Cas9 gene editing. A significant number of ASD and SZ candidate genes were among those that were differentially expressed in a comparison of heterozygous KO lines (CHD8+/-) vs isogenic controls (CHD8+/-), including the SZ and bipolar disorder (BD) candidate gene TCF4, which was markedly upregulated in CHD8+/- neuronal cells.

Schizophrenia (SZ) is a highly heterogeneous disorder in both its symptoms and risk factors. One of the most prevalent genetic risk factors for SZ is the hemizygous microdeletion at chromosome 22q11.2 (22q11DS) that confers a 25-fold increased risk. Six of the genes directly disrupted in 22qDS encode for mitochondrial-localizing proteins. Here, we test the hypothesis that stem cell-derived neurons from subjects with the 22q11DS and SZ have mitochondrial deficits relative to typically developing controls. Human iPSCs from four lines of affected subjects and five lines of controls were differentiated into forebrain-like excitatory neurons. In the patient group, we find significant reductions of ATP levels that appear to be secondary to reduced activity in oxidative phosphorylation complexes I and IV. Protein products of mitochondrial-encoded genes are also reduced. As one of the genes deleted in the 22q11.2 region is MRPL40, a component of the mitochondrial ribosome, we generated a heterozygous mutation of MRPL40 in a healthy control iPSC line. Relative to its isogenic control, this line shows similar deficits in mitochondrial DNA-encoded proteins, ATP level, and complex I and IV activity. These results suggest that in the 22q11DS MRPL40 heterozygosity leads to reduced mitochondria ATP production secondary to altered mitochondrial protein levels. Such defects could have profound effects on neuronal function in vivo.

Stochastic processes and imprinting, along with genetic factors, lead to monoallelic or allele-biased gene expression. Stochastic monoallelic expression fine-tunes information processing in immune cells and the olfactory system, and imprinting plays an important role in development. Recent studies suggest that both stochastic events and imprinting may be more widespread than previously considered. We are interested in allele-biased gene expression occurring in the brain because parent-of-origin effects suggestive of imprinting appear to play a role in the transmission of schizophrenia (SZ) and autism spectrum disorders (ASD) in some families. In addition, allele-biased expression could help explain monozygotic (MZ) twin discordance and reduced penetrance. The ability to study allele-biased expression in human neurons has been transformed with the advent of induced pluripotent stem cell (iPSC) technology and next generation sequencing. Using transcriptome sequencing (RNA-Seq) we identified 801 genes in differentiating neurons that were expressed in an allele-biased manner. These included a number of putative SZ and ASD candidates, such as A2BP1 (RBFOX1), ERBB4, NLGN4X, NRG1, NRG3, NRXN1, and NLGN1. Overall, there was a modest enrichment for SZ and ASD candidate genes among those that showed evidence for allele-biased expression (chi-square, p = 0.02). In addition to helping explain MZ twin discordance and reduced penetrance, the capacity to group many candidate genes affecting a variety of molecular and cellular pathways under a common regulatory process - allele-biased expression - could have therapeutic implications.

Monoallelic expression of autosomal genes has been implicated in human psychiatric disorders. However, there is a paucity of allelic expression studies in human brain cells at the single cell and genome wide levels.

Lowe syndrome (LS) is caused by loss-of-function mutations in the X-linked gene OCRL, which codes for an inositol polyphosphate 5-phosphatase that plays a key role in endosome recycling, clathrin-coated pit formation, and actin polymerization. It is characterized by congenital cataracts, intellectual and developmental disability, and renal proximal tubular dysfunction. Patients are also at high risk for developing glaucoma and seizures. We recently developed induced pluripotent stem cell (iPSC) lines from three patients with LS who have hypomorphic variants affecting the 3' end of the gene, and their neurotypical brothers to serve as controls.