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.

Shank3/PROSAP2 gene mutations are associated with cognitive impairment ranging from mental retardation to autism. Shank3 is a large scaffold postsynaptic density protein implicated in dendritic spines and synapse formation; however, its specific functions have not been clearly demonstrated. We have used RNAi to knockdown Shank3 expression in neuronal cultures and showed that this treatment specifically reduced the synaptic expression of the metabotropic glutamate receptor 5 (mGluR5), but did not affect the expression of other major synaptic proteins. The functional consequence of Shank3 RNAi knockdown was impaired signaling via mGluR5, as shown by reduction in ERK1/2 and CREB phosphorylation induced by stimulation with (S)-3,5-dihydroxyphenylglycine (DHPG) as the agonist of mGluR5 receptors, impaired mGluR5-dependent synaptic plasticity (DHPG-induced long-term depression), and impaired mGluR5-dependent modulation of neural network activity. We also found morphological abnormalities in the structure of synapses (spine number, width, and length) and impaired glutamatergic synaptic transmission, as shown by reduction in the frequency of miniature excitatory postsynaptic currents (mEPSC). Notably, pharmacological augmentation of mGluR5 activity using 3-cyano-N-(1,3-diphenyl-1H-pyrazol-5-yl)-benzamide as the positive allosteric modulator of these receptors restored mGluR5-dependent signaling (DHPG-induced phosphorylation of ERK1/2) and normalized the frequency of mEPSCs in Shank3-knocked down neurons. These data demonstrate that a deficit in mGluR5-mediated intracellular signaling in Shank3 knockdown neurons can be compensated by 3-cyano-N-(1,3-diphenyl-1H-pyrazol-5-yl)-benzamide; this raises the possibility that pharmacological augmentation of mGluR5 activity represents a possible new therapeutic approach for patients with Shank3 mutations.

Exome sequencing has become a popular method to evaluate undirected mutagenesis experiments in mice. However, the most suitable mouse strain for the biological model may be relatively distant from the standard mouse reference genome. For pinpointing causative variants, a matching reference with gene annotations is essential, but not always readily available.

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.

Cactophilic Drosophila species provide a valuable model to study gene-environment interactions and ecological adaptation. Drosophila buzzatii and Drosophila mojavensis are two cactophilic species that belong to the repleta group, but have very different geographical distributions and primary host plants. To investigate the genomic basis of ecological adaptation, we sequenced the genome and developmental transcriptome of D. buzzatii and compared its gene content with that of D. mojavensis and two other noncactophilic Drosophila species in the same subgenus. The newly sequenced D. buzzatii genome (161.5 Mb) comprises 826 scaffolds (>3 kb) and contains 13,657 annotated protein-coding genes. Using RNA sequencing data of five life-stages we found expression of 15,026 genes, 80% protein-coding genes, and 20% noncoding RNA genes. In total, we detected 1,294 genes putatively under positive selection. Interestingly, among genes under positive selection in the D. mojavensis lineage, there is an excess of genes involved in metabolism of heterocyclic compounds that are abundant in Stenocereus cacti and toxic to nonresident Drosophila species. We found 117 orphan genes in the shared D. buzzatii-D. mojavensis lineage. In addition, gene duplication analysis identified lineage-specific expanded families with functional annotations associated with proteolysis, zinc ion binding, chitin binding, sensory perception, ethanol tolerance, immunity, physiology, and reproduction. In summary, we identified genetic signatures of adaptation in the shared D. buzzatii-D. mojavensis lineage, and in the two separate D. buzzatii and D. mojavensis lineages. Many of the novel lineage-specific genomic features are promising candidates for explaining the adaptation of these species to their distinct ecological niches.

T-cell acute lymphoblastic leukaemias (T-ALL) are aggressive malignant proliferations characterized by high relapse rates and great genetic heterogeneity. TAL1 is amongst the most frequently deregulated oncogenes. Yet, over half of the TAL1(+) cases lack TAL1 lesions, suggesting unrecognized (epi)genetic deregulation mechanisms. Here we show that TAL1 is normally silenced in the T-cell lineage, and that the polycomb H3K27me3-repressive mark is focally diminished in TAL1(+) T-ALLs. Sequencing reveals that >20% of monoallelic TAL1(+) patients without previously known alterations display microinsertions or RAG1/2-mediated episomal reintegration in a single site 5' to TAL1. Using 'allelic-ChIP' and CrispR assays, we demonstrate that such insertions induce a selective switch from H3K27me3 to H3K27ac at the inserted but not the germline allele. We also show that, despite a considerable mechanistic diversity, the mode of oncogenic TAL1 activation, rather than expression levels, impact on clinical outcome. Altogether, these studies establish site-specific epigenetic desilencing as a mechanism of oncogenic activation.

The molecular mechanisms underlying specification from embryonic stem cells (ESCs) and maintenance of neural progenitor cells (NPCs) are largely unknown. Recently, we reported that the Zuotin-related factor 1 (Zrf1) is necessary for chromatin displacement of the Polycomb-repressive complex 1 (PRC1). We found that Zrf1 is required for NPC specification from ESCs and that it promotes the expression of NPC markers, including the key regulator Pax6. Moreover, Zrf1 is essential to establish and maintain Wnt ligand expression levels, which are necessary for NPC self-renewal. Reactivation of proper Wnt signaling in Zrf1-depleted NPCs restores Pax6 expression and the self-renewal capacity. ESC-derived NPCs in vitro resemble most of the characteristics of the self-renewing NPCs located in the developing embryonic cortex, which are termed radial glial cells (RGCs). Depletion of Zrf1 in vivo impairs the expression of key self-renewal regulators and Wnt ligand genes in RGCs. Thus, we demonstrate that Zrf1 plays an essential role in NPC generation and maintenance.

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.

The REvolutionary Approaches and Devices for Nucleic Acid analysis (READNA) project received funding from the European Commission for 41/2 years. The objectives of the project revolved around technological developments in nucleic acid analysis. The project partners have discovered, created and developed a huge body of insights into nucleic acid analysis, ranging from improvements and implementation of current technologies to the most promising sequencing technologies that constitute a 3(rd) and 4(th) generation of sequencing methods with nanopores and in situ sequencing, respectively.

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.

Here, we report DNA methylation and hydroxymethylation dynamics at nucleotide resolution using C/EBPα-enhanced reprogramming of B cells into induced pluripotent cells (iPSCs). We observed successive waves of hydroxymethylation at enhancers, concomitant with a decrease in DNA methylation, suggesting active demethylation. Consistent with this finding, ablation of the DNA demethylase Tet2 almost completely abolishes reprogramming. C/EBPα, Klf4, and Tfcp2l1 each interact with Tet2 and recruit the enzyme to specific DNA sites. During reprogramming, some of these sites maintain high levels of 5hmC, and enhancers and promoters of key pluripotency factors become demethylated as early as 1 day after Yamanaka factor induction. Surprisingly, methylation changes precede chromatin opening in distinct chromatin regions, including Klf4 bound sites, revealing a pioneer factor activity associated with alternation in DNA methylation. Rapid changes in hydroxymethylation similar to those in B cells were also observed during compound-accelerated reprogramming of fibroblasts into iPSCs, highlighting the generality of our observations.

Six extant species of non-human great apes are currently recognized: Sumatran and Bornean orangutans, eastern and western gorillas, and chimpanzees and bonobos [1]. However, large gaps remain in our knowledge of fine-scale variation in hominoid morphology, behavior, and genetics, and aspects of great ape taxonomy remain in flux. This is particularly true for orangutans (genus: Pongo), the only Asian great apes and phylogenetically our most distant relatives among extant hominids [1]. Designation of Bornean and Sumatran orangutans, P. pygmaeus (Linnaeus 1760) and P. abelii (Lesson 1827), as distinct species occurred in 2001 [1, 2]. Here, we show that an isolated population from Batang Toru, at the southernmost range limit of extant Sumatran orangutans south of Lake Toba, is distinct from other northern Sumatran and Bornean populations. By comparing cranio-mandibular and dental characters of an orangutan killed in a human-animal conflict to those of 33 adult male orangutans of a similar developmental stage, we found consistent differences between the Batang Toru individual and other extant Ponginae. Our analyses of 37 orangutan genomes provided a second line of evidence. Model-based approaches revealed that the deepest split in the evolutionary history of extant orangutans occurred ∼3.38 mya between the Batang Toru population and those to the north of Lake Toba, whereas both currently recognized species separated much later, about 674 kya. Our combined analyses support a new classification of orangutans into three extant species. The new species, Pongo tapanuliensis, encompasses the Batang Toru population, of which fewer than 800 individuals survive. VIDEO ABSTRACT.

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.

Determining the effect of DNA methylation on chromatin structure and function in higher organisms is challenging due to the extreme complexity of epigenetic regulation. We studied a simpler model system, budding yeast, that lacks DNA methylation machinery making it a perfect model system to study the intrinsic role of DNA methylation in chromatin structure and function. We expressed the murine DNA methyltransferases in Saccharomyces cerevisiae and analyzed the correlation between DNA methylation, nucleosome positioning, gene expression and 3D genome organization. Despite lacking the machinery for positioning and reading methylation marks, induced DNA methylation follows a conserved pattern with low methylation levels at the 5' end of the gene increasing gradually toward the 3' end, with concentration of methylated DNA in linkers and nucleosome free regions, and with actively expressed genes showing low and high levels of methylation at transcription start and terminating sites respectively, mimicking the patterns seen in mammals. We also see that DNA methylation increases chromatin condensation in peri-centromeric regions, decreases overall DNA flexibility, and favors the heterochromatin state. Taken together, these results demonstrate that methylation intrinsically modulates chromatin structure and function even in the absence of cellular machinery evolved to recognize and process the methylation signal.

Bringing together cancer genomes from different projects increases power and allows the investigation of pan-cancer, molecular mechanisms. However, working with whole genomes sequenced over several years in different sequencing centres requires a framework to compare the quality of these sequences. We used the Pan-Cancer Analysis of Whole Genomes cohort as a test case to construct such a framework. This cohort contains whole cancer genomes of 2832 donors from 18 sequencing centres. We developed a non-redundant set of five quality control (QC) measurements to establish a star rating system. These QC measures reflect known differences in sequencing protocol and provide a guide to downstream analyses and allow for exclusion of samples of poor quality. We have found that this is an effective framework of quality measures. The implementation of the framework is available at: https://dockstore.org/containers/quay.io/jwerner_dkfz/pancanqc:1.2.2 .