Despite rapid progress in characterizing transcription factor-driven reprogramming of somatic cells to an induced pluripotent stem cell (iPSC) state, many mechanistic questions still remain. To gain insight into the earliest events in the reprogramming process, we systematically analyzed the transcriptional and epigenetic changes that occur during early factor induction after discrete numbers of divisions. We observed rapid, genome-wide changes in the euchromatic histone modification, H3K4me2, at more than a thousand loci including large subsets of pluripotency-related or developmentally regulated gene promoters and enhancers. In contrast, patterns of the repressive H3K27me3 modification remained largely unchanged except for focused depletion specifically at positions where H3K4 methylation is gained. These chromatin regulatory events precede transcriptional changes within the corresponding loci. Our data provide evidence for an early, organized, and population-wide epigenetic response to ectopic reprogramming factors that clarify the temporal order through which somatic identity is reset during reprogramming.

Reprogramming to pluripotency is a low-efficiency process at the population level. Despite notable advances to molecularly characterize key steps, several fundamental aspects remain poorly understood, including when the potential to reprogram is first established. Here, we apply live-cell imaging combined with a novel statistical approach to infer when somatic cells become fated to generate downstream pluripotent progeny. By tracing cell lineages from several divisions before factor induction through to pluripotent colony formation, we find that pre-induction sister cells acquire similar outcomes. Namely, if one daughter cell contributes to a lineage that generates induced pluripotent stem cells (iPSCs), its paired sibling will as well. This result suggests that the potential to reprogram is predetermined within a select subpopulation of cells and heritable, at least over the short term. We also find that expanding cells over several divisions prior to factor induction does not increase the per-lineage likelihood of successful reprogramming, nor is reprogramming fate correlated to neighboring cell identity or cell-specific reprogramming factor levels. By perturbing the epigenetic state of somatic populations with Ezh2 inhibitors prior to factor induction, we successfully modulate the fraction of iPSC-forming lineages. Our results therefore suggest that reprogramming potential may in part reflect preexisting epigenetic heterogeneity that can be tuned to alter the cellular response to factor induction.

Differentiation of human embryonic stem cells (hESCs) provides a unique opportunity to study the regulatory mechanisms that facilitate cellular transitions in a human context. To that end, we performed comprehensive transcriptional and epigenetic profiling of populations derived through directed differentiation of hESCs representing each of the three embryonic germ layers. Integration of whole-genome bisulfite sequencing, chromatin immunoprecipitation sequencing, and RNA sequencing reveals unique events associated with specification toward each lineage. Lineage-specific dynamic alterations in DNA methylation and H3K4me1 are evident at putative distal regulatory elements that are frequently bound by pluripotency factors in the undifferentiated hESCs. In addition, we identified germ-layer-specific H3K27me3 enrichment at sites exhibiting high DNA methylation in the undifferentiated state. A better understanding of these initial specification events will facilitate identification of deficiencies in current approaches, leading to more faithful differentiation strategies as well as providing insights into the rewiring of human regulatory programs during cellular transitions.

Research on human pluripotent stem cells has been hampered by the lack of a standardized, quantitative, scalable assay of pluripotency. We previously described an assay called ScoreCard that used gene expression signatures to quantify differentiation efficiency. Here we report an improved version of the assay based on qPCR that enables faster, more quantitative assessment of functional pluripotency. We provide an in-depth characterization of the revised signature panel (commercially available as the TaqMan hPSC Scorecard Assay) through embryoid body and directed differentiation experiments as well as a detailed comparison to the teratoma assay. We further show that the improved ScoreCard enables a wider range of applications, such as screening of small molecules, genetic perturbations and assessment of culture conditions. Our approach can be extended beyond stem cell applications to characterize and assess the utility of other cell types and lineages.

Polycomb Repressive Complex 2 (PRC2) function and DNA methylation (DNAme) are typically correlated with gene repression. Here, we show that PRC2 is required to maintain expression of maternal microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) from the Gtl2-Rian-Mirg locus, which is essential for full pluripotency of iPSCs. In the absence of PRC2, the entire locus becomes transcriptionally repressed due to gain of DNAme at the intergenic differentially methylated regions (IG-DMRs). Furthermore, we demonstrate that the IG-DMR serves as an enhancer of the maternal Gtl2-Rian-Mirg locus. Further analysis reveals that PRC2 interacts physically with Dnmt3 methyltransferases and reduces recruitment to and subsequent DNAme at the IG-DMR, thereby allowing for proper expression of the maternal Gtl2-Rian-Mirg locus. Our observations are consistent with a mechanism through which PRC2 counteracts the action of Dnmt3 methyltransferases at an imprinted locus required for full pluripotency.

The cost of whole-genome bisulfite sequencing (WGBS) remains a bottleneck for many studies and it is therefore imperative to extract as much information as possible from a given dataset. This is particularly important because even at the recommend 30X coverage for reference methylomes, up to 50% of high-resolution features such as differentially methylated positions (DMPs) cannot be called with current methods as determined by saturation analysis. To address this limitation, we have developed a tool that dynamically segments WGBS methylomes into blocks of comethylation (COMETs) from which lost information can be recovered in the form of differentially methylated COMETs (DMCs). Using this tool, we demonstrate recovery of ∼30% of the lost DMP information content as DMCs even at very low (5X) coverage. This constitutes twice the amount that can be recovered using an existing method based on differentially methylated regions (DMRs). In addition, we explored the relationship between COMETs and haplotypes in lymphoblastoid cell lines of African and European origin. Using best fit analysis, we show COMETs to be correlated in a population-specific manner, suggesting that this type of dynamic segmentation may be useful for integrated (epi)genome-wide association studies in the future.

Despite their central importance in mammalian development, the mechanisms that regulate the DNA methylation machinery and thereby the generation of genomic methylation patterns are still poorly understood. Here, we identify the 5mC-binding protein MeCP2 as a direct and strong interactor of DNA methyltransferase 3 (DNMT3) proteins. We mapped the interaction interface to the transcriptional repression domain of MeCP2 and the ADD domain of DNMT3A and find that binding of MeCP2 strongly inhibits the activity of DNMT3A in vitro. This effect was reinforced by cellular studies where a global reduction of DNA methylation levels was observed after overexpression of MeCP2 in human cells. By engineering conformationally locked DNMT3A variants as novel tools to study the allosteric regulation of this enzyme, we show that MeCP2 stabilizes the closed, autoinhibitory conformation of DNMT3A. Interestingly, the interaction with MeCP2 and its resulting inhibition were relieved by the binding of K4 unmodified histone H3 N-terminal tail to the DNMT3A-ADD domain. Taken together, our data indicate that the localization and activity of DNMT3A are under the combined control of MeCP2 and H3 tail modifications where, depending on the modification status of the H3 tail at the binding sites, MeCP2 can act as either a repressor or activator of DNA methylation.

Mutations in KCNQ2, which encodes a pore-forming K+ channel subunit responsible for neuronal M-current, cause neonatal epileptic encephalopathy, a complex disorder presenting with severe early-onset seizures and impaired neurodevelopment. The condition is exceptionally difficult to treat, partially because the effects of KCNQ2 mutations on the development and function of human neurons are unknown. Here, we used induced pluripotent stem cells (iPSCs) and gene editing to establish a disease model and measured the functional properties of differentiated excitatory neurons. We find that patient iPSC-derived neurons exhibit faster action potential repolarization, larger post-burst afterhyperpolarization and a functional enhancement of Ca2+-activated K+ channels. These properties, which can be recapitulated by chronic inhibition of M-current in control neurons, facilitate a burst-suppression firing pattern that is reminiscent of the interictal electroencephalography pattern in patients. Our findings suggest that dyshomeostatic mechanisms compound KCNQ2 loss-of-function leading to alterations in the neurodevelopmental trajectory of patient iPSC-derived neurons.

Stem-cell-derived tissues could transform disease research and therapy, yet most methods generate functionally immature products. We investigate how human pluripotent stem cells (hPSCs) differentiate into pancreatic islets in vitro by profiling DNA methylation, chromatin accessibility, and histone modification changes. We find that enhancer potential is reset upon lineage commitment and show how pervasive epigenetic priming steers endocrine cell fates. Modeling islet differentiation and maturation regulatory circuits reveals genes critical for generating endocrine cells and identifies circadian control as limiting for in vitro islet function. Entrainment to circadian feeding/fasting cycles triggers islet metabolic maturation by inducing cyclic synthesis of energy metabolism and insulin secretion effectors, including antiphasic insulin and glucagon pulses. Following entrainment, hPSC-derived islets gain persistent chromatin changes and rhythmic insulin responses with a raised glucose threshold, a hallmark of functional maturity, and function within days of transplantation. Thus, hPSC-derived tissues are amenable to functional improvement by circadian modulation.

Genetic and epigenetic intra-tumoral heterogeneity cooperate to shape the evolutionary course of cancer1. Chronic lymphocytic leukaemia (CLL) is a highly informative model for cancer evolution as it undergoes substantial genetic diversification and evolution after therapy2,3. The CLL epigenome is also an important disease-defining feature4,5, and growing populations of cells in CLL diversify by stochastic changes in DNA methylation known as epimutations6. However, previous studies using bulk sequencing methods to analyse the patterns of DNA methylation were unable to determine whether epimutations affect CLL populations homogeneously. Here, to measure the epimutation rate at single-cell resolution, we applied multiplexed single-cell reduced-representation bisulfite sequencing to B cells from healthy donors and patients with CLL. We observed that the common clonal origin of CLL results in a consistently increased epimutation rate, with low variability in the cell-to-cell epimutation rate. By contrast, variable epimutation rates across healthy B cells reflect diverse evolutionary ages across the trajectory of B cell differentiation, consistent with epimutations serving as a molecular clock. Heritable epimutation information allowed us to reconstruct lineages at high-resolution with single-cell data, and to apply this directly to patient samples. The CLL lineage tree shape revealed earlier branching and longer branch lengths than in normal B cells, reflecting rapid drift after the initial malignant transformation and a greater proliferative history. Integration of single-cell bisulfite sequencing analysis with single-cell transcriptomes and genotyping confirmed that genetic subclones mapped to distinct clades, as inferred solely on the basis of epimutation information. Finally, to examine potential lineage biases during therapy, we profiled serial samples during ibrutinib-associated lymphocytosis, and identified clades of cells that were preferentially expelled from the lymph node after treatment, marked by distinct transcriptional profiles. The single-cell integration of genetic, epigenetic and transcriptional information thus charts the lineage history of CLL and its evolution with therapy.

De novo establishment of DNA methylation is accomplished by DNMT3A and DNMT3B. Here, we analyze de novo DNA methylation in mouse embryonic fibroblasts (2i-MEFs) derived from DNA-hypomethylated 2i/L ES cells with genetic ablation of Dnmt3a or Dnmt3b. We identify 355 and 333 uniquely unmethylated genes in Dnmt3a and Dnmt3b knockout (KO) 2i-MEFs, respectively. We find that Dnmt3a is exclusively required for de novo methylation at both TSS regions and gene bodies of Polycomb group (PcG) target developmental genes, while Dnmt3b has a dominant role on the X chromosome. Consistent with this, tissue-specific DNA methylation at PcG target genes is substantially reduced in Dnmt3a KO embryos. Finally, we find that human patients with DNMT3 mutations exhibit reduced DNA methylation at regions that are hypomethylated in Dnmt3 KO 2i-MEFs. In conclusion, here we report a set of unique de novo DNA methylation target sites for both DNMT3 enzymes during mammalian development that overlap with hypomethylated sites in human patients.

Schizophrenia (SCZ) is a debilitating psychiatric disorder with a large genetic contribution; however, its neurodevelopmental substrates remain largely unknown. Modeling pathogenic processes in SCZ using human induced pluripotent stem cell-derived neurons (iNs) has emerged as a promising strategy. Copy number variants confer high genetic risk for SCZ, with duplication of the 16p11.2 locus increasing the risk 14.5-fold.

Long non-coding RNAs (lncRNAs) have emerged as fundamental regulators in various biological processes, including embryonic development and cellular differentiation. Despite much progress over the past decade, the genome-wide annotation of lncRNAs remains incomplete and many known non-coding loci are still poorly characterized. Here, we report the discovery of a previously unannotated lncRNA that is transcribed 230 kb upstream of the SOX17 gene and located within the same topologically associating domain. We termed it T-REX17 (Transcript Regulating Endoderm and activated by soX17) and show that it is induced following SOX17 activation but its expression is more tightly restricted to early definitive endoderm. Loss of T-REX17 affects crucial functions independent of SOX17 and leads to an aberrant endodermal transcriptome, signaling pathway deregulation and epithelial to mesenchymal transition defects. Consequently, cells lacking the lncRNA cannot further differentiate into more mature endodermal cell types. Taken together, our study identified and characterized T-REX17 as a transiently expressed and essential non-coding regulator in early human endoderm differentiation.

The segregation of definitive endoderm (DE) from bipotent mesendoderm progenitors leads to the formation of two distinct germ layers. Dissecting DE commitment and onset has been challenging as it occurs within a narrow spatiotemporal window in the embryo. Here, we employ a dual Bra/Sox17 reporter cell line to study DE onset dynamics. We find Sox17 expression initiates in vivo in isolated cells within a temporally restricted window. In 2D and 3D in vitro models, DE cells emerge from mesendoderm progenitors at a temporally regular, but spatially stochastic pattern, which is subsequently arranged by self-sorting of Sox17 + cells. A subpopulation of Bra-high cells commits to a Sox17+ fate independent of external Wnt signal. Self-sorting coincides with upregulation of E-cadherin but is not necessary for DE differentiation or proliferation. Our in vivo and in vitro results highlight basic rules governing DE onset and patterning through the commonalities and differences between these systems.

DNA and Histone 3 Lysine 27 methylation typically function as repressive modifications and operate within distinct genomic compartments. In mammals, the majority of the genome is kept in a DNA methylated state, whereas the Polycomb repressive complexes regulate the unmethylated CpG-rich promoters of developmental genes. In contrast to this general framework, the extra-embryonic lineages display non-canonical, globally intermediate DNA methylation levels, including disruption of local Polycomb domains. Here, to better understand this unusual landscape's molecular properties, we genetically and chemically perturbed major epigenetic pathways in mouse trophoblast stem cells. We find that the extra-embryonic epigenome reflects ongoing and dynamic de novo methyltransferase recruitment, which is continuously antagonized by Polycomb to maintain intermediate, locally disordered methylation. Despite its disorganized molecular appearance, our data point to a highly controlled equilibrium between counteracting repressors within extra-embryonic cells, one that can seemingly persist indefinitely without bistable features typically seen for embryonic forms of epigenetic regulation.

X-chromosome inactivation (XCI) balances gene expression between the sexes in female mammals. Shortly after fertilization, upregulation of Xist RNA from one X chromosome initiates XCI, leading to chromosome-wide gene silencing. XCI is maintained in all cell types, except the germ line and the pluripotent state where XCI is reversed. The mechanisms triggering Xist upregulation have remained elusive. Here we identify GATA transcription factors as potent activators of Xist. Through a pooled CRISPR activation screen in murine embryonic stem cells, we demonstrate that GATA1, as well as other GATA transcription factors can drive ectopic Xist expression. Moreover, we describe GATA-responsive regulatory elements in the Xist locus bound by different GATA factors. Finally, we show that GATA factors are essential for XCI induction in mouse preimplantation embryos. Deletion of GATA1/4/6 or GATA-responsive Xist enhancers in mouse zygotes effectively prevents Xist upregulation. We propose that the activity or complete absence of various GATA family members controls initial Xist upregulation, XCI maintenance in extra-embryonic lineages and XCI reversal in the epiblast.

DNA methylation contributes to the maintenance of genomic integrity in somatic cells, in part through the silencing of transposable elements. In this study, we use CRISPR-Cas9 technology to delete DNMT1, the DNA methyltransferase key for DNA methylation maintenance, in human neural progenitor cells (hNPCs). We observe that inactivation of DNMT1 in hNPCs results in viable, proliferating cells despite a global loss of DNA CpG-methylation. DNA demethylation leads to specific transcriptional activation and chromatin remodeling of evolutionarily young, hominoid-specific LINE-1 elements (L1s), while older L1s and other classes of transposable elements remain silent. The activated L1s act as alternative promoters for many protein-coding genes involved in neuronal functions, revealing a hominoid-specific L1-based transcriptional network controlled by DNA methylation that influences neuronal protein-coding genes. Our results provide mechanistic insight into the role of DNA methylation in silencing transposable elements in somatic human cells, as well as further implicating L1s in human brain development and disease.

Although thousands of large intergenic non-coding RNAs (lincRNAs) have been identified in mammals, few have been functionally characterized, leading to debate about their biological role. To address this, we performed loss-of-function studies on most lincRNAs expressed in mouse embryonic stem (ES) cells and characterized the effects on gene expression. Here we show that knockdown of lincRNAs has major consequences on gene expression patterns, comparable to knockdown of well-known ES cell regulators. Notably, lincRNAs primarily affect gene expression in trans. Knockdown of dozens of lincRNAs causes either exit from the pluripotent state or upregulation of lineage commitment programs. We integrate lincRNAs into the molecular circuitry of ES cells and show that lincRNA genes are regulated by key transcription factors and that lincRNA transcripts bind to multiple chromatin regulatory proteins to affect shared gene expression programs. Together, the results demonstrate that lincRNAs have key roles in the circuitry controlling ES cell state.

The developmental potential of human pluripotent stem cells suggests that they can produce disease-relevant cell types for biomedical research. However, substantial variation has been reported among pluripotent cell lines, which could affect their utility and clinical safety. Such cell-line-specific differences must be better understood before one can confidently use embryonic stem (ES) or induced pluripotent stem (iPS) cells in translational research. Toward this goal we have established genome-wide reference maps of DNA methylation and gene expression for 20 previously derived human ES lines and 12 human iPS cell lines, and we have measured the in vitro differentiation propensity of these cell lines. This resource enabled us to assess the epigenetic and transcriptional similarity of ES and iPS cells and to predict the differentiation efficiency of individual cell lines. The combination of assays yields a scorecard for quick and comprehensive characterization of pluripotent cell lines.

Over the past years we have witnessed an explosion in the generation of freely available genome-wide data sets, including maps of various histone modifications, transcription factor binding, DNase hypersensitivity, and DNA methylation, which provide valuable resources for data validation, exploration, and hypothesis generation. The goal of this review is to provide the reader with information on where to find many of the data sets and how to utilize the various (epi)genome browsers for display and initial analysis. We provide selected examples to highlight key features and demonstrate the application of these browsers to stem cell biology.