An increasing amount of evidence indicates that developmental programs are tightly regulated by the complex interplay between signaling pathways, as well as transcriptional and epigenetic processes. Here, we have uncovered coordination between transcriptional and morphogen cues to specify Merkel cells, poorly understood skin cells that mediate light touch sensations. In murine dorsal skin, Merkel cells are part of touch domes, which are skin structures consisting of specialized keratinocytes, Merkel cells, and afferent neurons, and are located exclusively around primary hair follicles. We show that the developing primary hair follicle functions as a niche required for Merkel cell specification. We find that intraepidermal Sonic hedgehog (Shh) signaling, initiated by the production of Shh ligand in the developing hair follicles, is required for Merkel cell specification. The importance of Shh for Merkel cell formation is further reinforced by the fact that Shh overexpression in embryonic epidermal progenitors leads to ectopic Merkel cells. Interestingly, Shh signaling is common to primary, secondary, and tertiary hair follicles, raising the possibility that there are restrictive mechanisms that regulate Merkel cell specification exclusively around primary hair follicles. Indeed, we find that loss of Polycomb repressive complex 2 (PRC2) in the epidermis results in the formation of ectopic Merkel cells that are associated with all hair types. We show that PRC2 loss expands the field of epidermal cells competent to differentiate into Merkel cells through the upregulation of key Merkel-differentiation genes, which are known PRC2 targets. Importantly, PRC2-mediated repression of the Merkel cell differentiation program requires inductive Shh signaling to form mature Merkel cells. Our study exemplifies how the interplay between epigenetic and morphogen cues regulates the complex patterning and formation of the mammalian skin structures.

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.

Multiple endocrine neoplasia type 1 (MEN1) syndrome results from mutations in the MEN1 gene and causes tumor formation via largely unknown mechanisms. Using a novel genome-wide methylation analysis, we studied tissues from MEN1-parathyroid tumors, Men1 knockout (KO) mice, and Men1 null mouse embryonic fibroblast (MEF) cell lines. We demonstrated that inactivation of menin (the protein product of MEN1) increases activity of DNA (cytosine-5)-methyltransferase 1 (DNMT1) by activating retinoblastoma-binding protein 5 (Rbbp5). The increased activity of DNMT1 mediates global DNA hypermethylation, which results in aberrant activation of the Wnt/β-catenin signaling pathway through inactivation of Sox regulatory genes. Our study provides important insights into the role of menin in DNA methylation and its impact on the pathogenesis of MEN1 tumor development.

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.

Gene expression correlates with local chromatin structure. Our studies have mapped histone post-translational modifications, RNA polymerase II (pol II), and transcription factor Pax6 in lens chromatin. These data represent the first genome-wide insights into the relationship between lens chromatin structure and lens transcriptomes and serve as an excellent source for additional data analysis and refinement. The principal lens proteins, the crystallins, are encoded by predominantly expressed mRNAs; however, the regulatory mechanisms underlying their high expression in the lens remain poorly understood.

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.

The Tet family of enzymes (Tet1/2/3) converts 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC). Mouse embryonic stem cells (mESCs) highly express Tet1 and have an elevated level of 5hmC. Tet1 has been implicated in ESC maintenance and lineage specification in vitro but its precise function in development is not well defined. To establish the role of Tet1 in pluripotency and development, we have generated Tet1 mutant mESCs and mice. Tet1(-/-) ESCs have reduced levels of 5hmC and subtle changes in global gene expression, and are pluripotent and support development of live-born mice in tetraploid complementation assay, but display skewed differentiation toward trophectoderm in vitro. Tet1 mutant mice are viable, fertile, and grossly normal, though some mutant mice have a slightly smaller body size at birth. Our data suggest that Tet1 loss leading to a partial reduction in 5hmC levels does not affect pluripotency in ESCs and is compatible with embryonic and postnatal development.

Pax6 encodes a specific DNA-binding transcription factor that regulates the development of multiple organs, including the eye, brain and pancreas. Previous studies have shown that Pax6 regulates the entire process of ocular lens development. In the developing forebrain, Pax6 is expressed in ventricular zone precursor cells and in specific populations of neurons; absence of Pax6 results in disrupted cell proliferation and cell fate specification in telencephalon. In the pancreas, Pax6 is essential for the differentiation of α-, β- and δ-islet cells. To elucidate molecular roles of Pax6, chromatin immunoprecipitation experiments combined with high-density oligonucleotide array hybridizations (ChIP-chip) were performed using three distinct sources of chromatin (lens, forebrain and β-cells). ChIP-chip studies, performed as biological triplicates, identified a total of 5,260 promoters occupied by Pax6. 1,001 (133) of these promoter regions were shared between at least two (three) distinct chromatin sources, respectively. In lens chromatin, 2,335 promoters were bound by Pax6. RNA expression profiling from Pax6⁺/⁻ lenses combined with in vivo Pax6-binding data yielded 76 putative Pax6-direct targets, including the Gaa, Isl1, Kif1b, Mtmr2, Pcsk1n, and Snca genes. RNA and ChIP data were validated for all these genes. In lens cells, reporter assays established Kib1b and Snca as Pax6 activated and repressed genes, respectively. In situ hybridization revealed reduced expression of these genes in E14 cerebral cortex. Moreover, we examined differentially expressed transcripts between E9.5 wild type and Pax6⁻/⁻ lens placodes that suggested Efnb2, Fat4, Has2, Nav1, and Trpm3 as novel Pax6-direct targets. Collectively, the present studies, through the identification of Pax6-direct target genes, provide novel insights into the molecular mechanisms of Pax6 gene control during mouse embryonic development. In addition, the present data demonstrate that Pax6 interacts preferentially with promoter regions in a tissue-specific fashion. Nevertheless, nearly 20% of the regions identified are accessible to Pax6 in multiple tissues.

The transcriptional repressor REST (RE1 silencing transcription factor, also called NRSF for neuron-restrictive silencing factor) binds to a conserved RE1 motif and represses many neuronal genes in non-neuronal cells. This transcriptional regulation is transacted by several nucleosome-modifying enzymes recruited by REST to RE1 sites, including histone deacetylases (for example, HDAC1/2), demethylases (for example, LSD1), and methyltransferases (for example, G9a).

We report the generation and analysis of functional data from multiple, diverse experiments performed on a targeted 1% of the human genome as part of the pilot phase of the ENCODE Project. These data have been further integrated and augmented by a number of evolutionary and computational analyses. Together, our results advance the collective knowledge about human genome function in several major areas. First, our studies provide convincing evidence that the genome is pervasively transcribed, such that the majority of its bases can be found in primary transcripts, including non-protein-coding transcripts, and those that extensively overlap one another. Second, systematic examination of transcriptional regulation has yielded new understanding about transcription start sites, including their relationship to specific regulatory sequences and features of chromatin accessibility and histone modification. Third, a more sophisticated view of chromatin structure has emerged, including its inter-relationship with DNA replication and transcriptional regulation. Finally, integration of these new sources of information, in particular with respect to mammalian evolution based on inter- and intra-species sequence comparisons, has yielded new mechanistic and evolutionary insights concerning the functional landscape of the human genome. Together, these studies are defining a path for pursuit of a more comprehensive characterization of human genome function.

Tissue regeneration relies on resident stem cells (SCs), whose activity and lineage choices are influenced by the microenvironment. Exploiting the synchronized, cyclical bouts of tissue regeneration in hair follicles (HFs), we investigate how microenvironment dynamics shape the emergence of stem cell lineages. Employing epigenetic and ChIP-seq profiling, we uncover how signal-dependent transcription factors couple spatiotemporal cues to chromatin dynamics, thereby choreographing stem cell lineages. Using enhancer-driven reporters, mutagenesis, and genetics, we show that simultaneous BMP-inhibitory and WNT signals set the stage for lineage choices by establishing chromatin platforms permissive for diversification. Mechanistically, when binding of BMP effector pSMAD1 is relieved, enhancers driving HF-stem cell master regulators are silenced. Concomitantly, multipotent, lineage-fated enhancers silent in HF-stem cells become activated by exchanging WNT effectors TCF3/4 for LEF1. Throughout regeneration, lineage enhancers continue reliance upon LEF1 but then achieve specificity by accommodating additional incoming signaling effectors. Barriers to progenitor plasticity increase when diverse, signal-sensitive transcription factors shape LEF1-regulated enhancer dynamics.

Self-renewal and differentiation of adult stem cells are tightly regulated partly through configuration of chromatin structure by chromatin remodelers. Using knockout mice, we here demonstrate that bromodomain PHD finger transcription factor (BPTF), a component of the nucleosome remodeling factor (NURF) chromatin-remodeling complex, is essential for maintaining the population size of hematopoietic stem/progenitor cells (HSPCs), including long-term hematopoietic stem cells (HSCs). Bptf-deficient HSCs are defective in reconstituted hematopoiesis, and hematopoietic-specific knockout of Bptf caused profound defects including bone marrow failure and anemia. Genome-wide transcriptome profiling revealed that BPTF loss caused downregulation of HSC-specific gene-expression programs, which contain several master transcription factors (Meis1, Pbx1, Mn1, and Lmo2) required for HSC maintenance and self-renewal. Furthermore, we show that BPTF potentiates the chromatin accessibility of key HSC "stemness" genes. These results demonstrate an essential requirement of the chromatin remodeler BPTF and NURF for activation of "stemness" gene-expression programs and proper function of adult HSCs.

Gata3 is a DNA-binding transcription factor involved in cellular differentiation in a variety of tissues including inner ear, hair follicle, kidney, mammary gland and T-cells. In a previous study in 2009, Maeda et al. (Dev. Dyn.238, 2280-2291; doi:10.1002/dvdy.22035) found that Gata3 mutants could be rescued from midgestational lethality by the expression of a Gata3 transgene in sympathoadrenal neuroendocrine cells. The rescued embryos clearly showed multiple defects in lens fibre cell differentiation. To determine whether these defects were truly due to the loss of Gata3 expression in the lens, we generated a lens-specific Gata3 loss-of-function model. Analogous to the previous findings, our Gata3 null embryos showed abnormal regulation of cell cycle exit during lens fibre cell differentiation, marked by reduction in the expression of the cyclin-dependent kinase inhibitors Cdkn1b/p27 and Cdkn1c/p57, and the retention of nuclei accompanied by downregulation of Dnase IIβ. Comparisons of transcriptomes between control and mutated lenses by RNA-Seq revealed dysregulation of lens-specific crystallin genes and intermediate filament protein Bfsp2. Both Cdkn1b/p27 and Cdkn1c/p57 loci are occupied in vivo by Gata3, as well as Prox1 and c-Jun, in lens chromatin. Collectively, our studies suggest that Gata3 regulates lens differentiation through the direct regulation of the Cdkn1b/p27and Cdkn1c/p57 expression, and the direct/or indirect transcriptional control of Bfsp2 and Dnase IIβ.

Tet enzymes (Tet1/2/3) oxidize 5-methylcytosine to promote DNA demethylation and partner with chromatin modifiers to regulate gene expression. Tet1 is highly expressed in embryonic stem cells (ESCs), but its enzymatic and non-enzymatic roles in gene regulation are not dissected. We have generated Tet1 catalytically inactive (Tet1m/m) and knockout (Tet1-/-) ESCs and mice to study these functions. Loss of Tet1, but not loss of its catalytic activity, caused aberrant upregulation of bivalent (H3K4me3+; H3K27me3+) developmental genes, leading to defects in differentiation. Wild-type and catalytic-mutant Tet1 occupied similar genomic loci which overlapped with H3K27 tri-methyltransferase PRC2 and the deacetylase complex Sin3a at promoters of bivalent genes and with the helicase Chd4 at active genes. Loss of Tet1, but not loss of its catalytic activity, impaired enrichment of PRC2 and Sin3a at bivalent promoters leading to reduced H3K27 trimethylation and deacetylation, respectively, in absence of any changes in DNA methylation. Tet1-/-, but not Tet1m/m, embryos expressed higher levels of Gata6 and were developmentally delayed. Thus, the critical functions of Tet1 in ESCs and early development are mediated through its non-catalytic roles in regulating H3K27 modifications to silence developmental genes, and are more important than its catalytic functions in DNA demethylation.