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.

Human GABAergic interneurons (GIN) are implicated in normal brain function and in numerous mental disorders. However, the generation of functional human GIN subtypes from human pluripotent stem cells (hPSCs) has not been established. By expressing LHX6, a transcriptional factor that is critical for GIN development, we induced hPSCs to form GINs, including somatostatin (SST, 29%) and parvalbumin (PV, 21%) neurons. Our RNAseq results also confirmed the alteration of GIN identity with the overexpression of LHX6. Five months after transplantation into the mouse brain, the human GABA precursors generated increased population of SST and PV neurons by overexpressing LHX6. Importantly, the grafted human GINs exhibited functional electrophysiological properties and even fast-spiking-like action potentials. Thus, expression of the single transcription factor LHX6 under our GIN differentiation condition is sufficient to robustly induce human PV and SST subtypes.

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.

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.

Mammalian aging is a highly complex process, a full mechanistic understanding of which is still lacking. One way to help understand the molecular changes underlying aging is through a comprehensive analysis of the transcriptome, the primary determinant of age-related phenotypic diversity. Previous studies have relied on microarray analysis to examine gene expression profiles in different tissues of aging organisms. However, studies have shown microarray-based transcriptional profiling is less accurate and not fully capable of capturing certain intricacies of the global transcriptome.

The 22q11.2 deletion syndrome (22q11.2DS; velo-cardio-facial syndrome; DiGeorge syndrome) is a congenital anomaly disorder in which haploinsufficiency of TBX1, encoding a T-box transcription factor, is the major candidate for cardiac outflow tract (OFT) malformations. Inactivation of Tbx1 in the anterior heart field (AHF) mesoderm in the mouse results in premature expression of pro-differentiation genes and a persistent truncus arteriosus (PTA) in which septation does not form between the aorta and pulmonary trunk. Canonical Wnt/β-catenin has major roles in cardiac OFT development that may act upstream of Tbx1. Consistent with an antagonistic relationship, we found the opposite gene expression changes occurred in the AHF in β-catenin loss of function embryos compared to Tbx1 loss of function embryos, providing an opportunity to test for genetic rescue. When both alleles of Tbx1 and one allele of β-catenin were inactivated in the Mef2c-AHF-Cre domain, 61% of them (n = 34) showed partial or complete rescue of the PTA defect. Upregulated genes that were oppositely changed in expression in individual mutant embryos were normalized in significantly rescued embryos. Further, β-catenin was increased in expression when Tbx1 was inactivated, suggesting that there may be a negative feedback loop between canonical Wnt and Tbx1 in the AHF to allow the formation of the OFT. We suggest that alteration of this balance may contribute to variable expressivity in 22q11.2DS.

Meiosis entry during spermatogenesis requires reprogramming from mitotic to meiotic gene expression profiles. Transcriptional regulation has been extensively studied in meiosis entry, but gain of function for master transcription factors is insufficient to down-regulate mitotic genes. RNA helicase YTHDC2 and its partner MEIOC emerge as essential posttranscriptional regulators of meiotic entry. However, it is unclear what governs the RNA binding specificity of YTHDC2/MEIOC. Here, we identified RNA binding protein RBM46 as a component of the YTHDC2/MEIOC complex. Testis-specific Rbm46 knockout in mice causes infertility with defective mitotic-to-meiotic transition, phenocopying global Ythdc2 or Meioc knockout. RBM46 binds to 3' UTR of mitotic transcripts within 100 nucleotides from YTHDC2 U-rich motifs and targets these transcripts for degradation. Dysregulated RBM46 expression is associated with human male fertility disorders. These findings establish the RBM46/YTHDC2/MEIOC complex as the major posttranscriptional regulator responsible for down-regulating mitotic transcripts during meiosis entry in mammalian spermatogenesis, with implications for understanding meiosis-related fertility disorders.

The transition from meiotic spermatocytes to postmeiotic haploid germ cells constitutes an essential step in spermatogenesis. The epigenomic regulatory mechanisms underlying this transition remain unclear. Here, we find a prominent transcriptomic switch from the late spermatocytes to the early round spermatids during the meiotic-to-postmeiotic transition, which is associated with robust histone acetylation changes across the genome. Among histone deacetylases (HDACs) and acetyltransferases, we find that HDAC3 is selectively expressed in the late meiotic and early haploid stages. Three independent mouse lines with the testis-specific knockout of HDAC3 show infertility and defects in meiotic exit with an arrest at the late stage of meiosis or early stage of round spermatids. Stage-specific RNA-seq and histone acetylation ChIP-seq analyses reveal that HDAC3 represses meiotic/spermatogonial genes and activates postmeiotic haploid gene programs during meiotic exit, with associated histone acetylation alterations. Unexpectedly, abolishing HDAC3 catalytic activity by missense mutations in the nuclear receptor corepressor (NCOR or SMRT) does not cause infertility, despite causing histone hyperacetylation as HDAC3 knockout, demonstrating that HDAC3 enzyme activity is not required for spermatogenesis. Motif analysis of the HDAC3 cistrome in the testes identified SOX30, which has a similar spatiotemporal expression pattern as HDAC3 during spermatogenesis. Depletion of SOX30 in the testes abolishes the genomic recruitment of the HDAC3 to the binding sites. Collectively, these results establish the SOX30/HDAC3 signaling as a key regulator of the transcriptional program in a deacetylase-independent manner during the meiotic-to-postmeiotic transition in spermatogenesis.

Genetic biomarkers of lung cancer (LC) susceptibility may provide a basis for treatment and prevention. This study analyzed an association between SNPs (single nucleotide polymorphisms) in the complementary region of the 3'-UTR (3' untranslated region) of microRNAs of the gene RIPK1 (receptor-interacting serine/threonine-protein kinase 1) and LC among an adult Han Chinese population aged younger than 60 years. Also explored the effect of regulation of the RIPK1 gene via rs17548629 and microRNA-1197 on the occurrence of LC.

DNA methylation is a major epigenetic mechanism altering gene expression in development and disease. However, its role in the regulation of gene expression during heart development is incompletely understood. The aim of this study is to reveal DNA methylation in mouse embryonic hearts and its role in regulating gene expression during heart development.

RNA regulation by RNA-binding proteins (RBPs) involve extremely complicated mechanisms. MOV10 and MOV10L1 are two homologous RNA helicases implicated in distinct intracellular pathways. MOV10L1 participates specifically in Piwi-interacting RNA (piRNA) biogenesis and protects mouse male fertility. In contrast, the functional complexity of MOV10 remains incompletely understood, and its role in the mammalian germline is unknown. Here, we report a study of the biological and molecular functions of the RNA helicase MOV10 in mammalian male germ cells.

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.

Malignant peripheral nerve sheath tumors (MPNSTs) represent a group of highly aggressive soft-tissue sarcomas that may occur sporadically, in association with neurofibromatosis type I (NF1 associated) or after radiotherapy. Using comprehensive genomic approaches, we identified loss-of-function somatic alterations of the Polycomb repressive complex 2 (PRC2) components (EED or SUZ12) in 92% of sporadic, 70% of NF1-associated and 90% of radiotherapy-associated MPNSTs. MPNSTs with PRC2 loss showed complete loss of trimethylation at lysine 27 of histone H3 (H3K27me3) and aberrant transcriptional activation of multiple PRC2-repressed homeobox master regulators and their regulated developmental pathways. Introduction of the lost PRC2 component in a PRC2-deficient MPNST cell line restored H3K27me3 levels and decreased cell growth. Additionally, we identified frequent somatic alterations of CDKN2A (81% of all MPNSTs) and NF1 (72% of non-NF1-associated MPNSTs), both of which significantly co-occur with PRC2 alterations. The highly recurrent and specific inactivation of PRC2 components, NF1 and CDKN2A highlights their critical and potentially cooperative roles in MPNST pathogenesis.