DNA methylation in mammalian development and disease.

The DNA methylation field has matured from a phase of discovery and genomic characterization to one seeking deeper functional understanding of how this modification contributes to development, ageing and disease. In particular, the past decade has seen many exciting mechanistic discoveries that have substantially expanded our appreciation for how this generic, evolutionarily ancient modification can be incorporated into robust epigenetic codes. Here, we summarize the current understanding of the distinct DNA methylation landscapes that emerge over the mammalian lifespan and discuss how they interact with other regulatory layers to support diverse genomic functions. We then review the rising interest in alternative patterns found during senescence and the somatic transition to cancer. Alongside advancements in single-cell and long-read sequencing technologies, the collective insights made across these fields offer new opportunities to connect the biochemical and genetic features of DNA methylation to cell physiology, developmental potential and phenotype.

Cohesin-mediated 3D contacts tune enhancer-promoter regulation.

Enhancers are key drivers of gene regulation thought to act via 3D physical interactions with the promoters of their target genes. However, genome-wide depletions of architectural proteins such as cohesin result in only limited changes in gene expression, despite a loss of contact domains and loops. Consequently, the role of cohesin and 3D contacts in enhancer function remains debated. Here, we developed CRISPRi of regulatory elements upon degron operation (CRUDO), a novel approach to measure how changes in contact frequency impact enhancer effects on target genes by perturbing enhancers with CRISPRi and measuring gene expression in the presence or absence of cohesin. We systematically perturbed all 1,039 candidate enhancers near five cohesin-dependent genes and identified 34 enhancer-gene regulatory interactions. Of 26 regulatory interactions with sufficient statistical power to evaluate cohesin dependence, 18 show cohesin-dependent effects. A decrease in enhancer-promoter contact frequency upon removal of cohesin is frequently accompanied by a decrease in the regulatory effect of the enhancer on gene expression, consistent with a contact-based model for enhancer function. However, changes in contact frequency and regulatory effects on gene expression vary as a function of distance, with distal enhancers (e.g., >50Kb) experiencing much larger changes than proximal ones (e.g., <50Kb). Because most enhancers are located close to their target genes, these observations can explain how only a small subset of genes - those with strong distal enhancers - are sensitive to cohesin. Together, our results illuminate how 3D contacts, influenced by both cohesin and genomic distance, tune enhancer effects on gene expression.

Extraembryonic gut endoderm cells undergo programmed cell death during development.

Despite a distinct developmental origin, extraembryonic cells in mice contribute to gut endoderm and converge to transcriptionally resemble their embryonic counterparts. Notably, all extraembryonic progenitors share a non-canonical epigenome, raising several pertinent questions, including whether this landscape is reset to match the embryonic regulation and if extraembryonic cells persist into later development. Here we developed a two-colour lineage-tracing strategy to track and isolate extraembryonic cells over time. We find that extraembryonic gut cells display substantial memory of their developmental origin including retention of the original DNA methylation landscape and resulting transcriptional signatures. Furthermore, we show that extraembryonic gut cells undergo programmed cell death and neighbouring embryonic cells clear their remnants via non-professional phagocytosis. By midgestation, we no longer detect extraembryonic cells in the wild-type gut, whereas they persist and differentiate further in p53-mutant embryos. Our study provides key insights into the molecular and developmental fate of extraembryonic cells inside the embryo.

TET activity safeguards pluripotency throughout embryonic dormancy.

Dormancy is an essential biological process for the propagation of many life forms through generations and stressful conditions. Early embryos of many mammals are preservable for weeks to months within the uterus in a dormant state called diapause, which can be induced in vitro through mTOR inhibition. Cellular strategies that safeguard original cell identity within the silent genomic landscape of dormancy are not known. Here we show that the protection of cis-regulatory elements from silencing is key to maintaining pluripotency in the dormant state. We reveal a TET-transcription factor axis, in which TET-mediated DNA demethylation and recruitment of methylation-sensitive transcription factor TFE3 drive transcriptionally inert chromatin adaptations during dormancy transition. Perturbation of TET activity compromises pluripotency and survival of mouse embryos under dormancy, whereas its enhancement improves survival rates. Our results reveal an essential mechanism for propagating the cellular identity of dormant cells, with implications for regeneration and disease.

Reconstructing axial progenitor field dynamics in mouse stem cell-derived embryoids.

Embryogenesis requires substantial coordination to translate genetic programs to the collective behavior of differentiating cells, but understanding how cellular decisions control tissue morphology remains conceptually and technically challenging. Here, we combine continuous Cas9-based molecular recording with a mouse embryonic stem cell-based model of the embryonic trunk to build single-cell phylogenies that describe the behavior of transient, multipotent neuro-mesodermal progenitors (NMPs) as they commit into neural and somitic cell types. We find that NMPs show subtle transcriptional signatures related to their recent differentiation and contribute to downstream lineages through a surprisingly broad distribution of individual fate outcomes. Although decision-making can be heavily influenced by environmental cues to induce morphological phenotypes, axial progenitors intrinsically mature over developmental time to favor the neural lineage. Using these data, we present an experimental and analytical framework for exploring the non-homeostatic dynamics of transient progenitor populations as they shape complex tissues during critical developmental windows.

Basal delamination during mouse gastrulation primes pluripotent cells for differentiation.

The blueprint of the mammalian body plan is laid out during gastrulation, when a trilaminar embryo is formed. This process entails a burst of proliferation, the ingression of embryonic epiblast cells at the primitive streak, and their priming toward primitive streak fates. How these different events are coordinated remains unknown. Here, we developed and characterized a 3D culture of self-renewing mouse embryonic cells that captures the main transcriptional and architectural features of the early gastrulating mouse epiblast. Using this system in combination with microfabrication and in vivo experiments, we found that proliferation-induced crowding triggers delamination of cells that express high levels of the apical polarity protein aPKC. Upon delamination, cells become more sensitive to Wnt signaling and upregulate the expression of primitive streak markers such as Brachyury. This mechanistic coupling between ingression and differentiation ensures that the right cell types become specified at the right place during embryonic development.

Lambda3: homology search for protein, nucleotide, and bisulfite-converted sequences.

MOTIVATION: Local alignments of query sequences in large databases represent a core part of metagenomic studies and facilitate homology search. Following the development of NCBI Blast, many applications aimed to provide faster and equally sensitive local alignment frameworks. Most applications focus on protein alignments, while only few also facilitate DNA-based searches. None of the established programs allow searching DNA sequences from bisulfite sequencing experiments commonly used for DNA methylation profiling, for which specific alignment strategies need to be implemented. RESULTS: Here, we introduce Lambda3, a new version of the local alignment application Lambda. Lambda3 is the first solution that enables the search of protein, nucleotide as well as bisulfite-converted nucleotide query sequences. Its protein mode achieves comparable performance to that of the highly optimized protein alignment application Diamond, while the nucleotide mode consistently outperforms established local nucleotide aligners. Combined, Lambda3 presents a universal local alignment framework that enables fast and sensitive homology searches for a wide range of use-cases. AVAILABILITY AND IMPLEMENTATION: Lambda3 is free and open-source software publicly available at https://github.com/seqan/lambda/.

GATA transcription factors drive initial Xist upregulation after fertilization through direct activation of long-range enhancers.

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.

AML with complex karyotype: extreme genomic complexity revealed by combined long-read sequencing and Hi-C technology.

Acute myeloid leukemia with complex karyotype (CK-AML) is associated with poor prognosis, which is only in part explained by underlying TP53 mutations. Especially in the presence of complex chromosomal rearrangements, such as chromothripsis, the outcome of CK-AML is dismal. However, this degree of complexity of genomic rearrangements contributes to the leukemogenic phenotype and treatment resistance of CK-AML remains largely unknown. Applying an integrative workflow for the detection of structural variants (SVs) based on Oxford Nanopore (ONT) genomic DNA long-read sequencing (gDNA-LRS) and high-throughput chromosome confirmation capture (Hi-C) in a well-defined cohort of CK-AML identified regions with an extreme density of SVs. These rearrangements consisted to a large degree of focal amplifications enriched in the proximity of mammalian-wide interspersed repeat elements, which often result in oncogenic fusion transcripts, such as USP7::MVD, or the deregulation of oncogenic driver genes as confirmed by RNA-seq and ONT direct complementary DNA sequencing. We termed this novel phenomenon chromocataclysm. Thus, our integrative SV detection workflow combing gDNA-LRS and Hi-C enables to unravel complex genomic rearrangements at a very high resolution in regions hard to analyze by conventional sequencing technology, thereby providing an important tool to identify novel important drivers underlying cancer with complex karyotypic changes.

H3K36 methylation maintains cell identity by regulating opposing lineage programmes.

The epigenetic mechanisms that maintain differentiated cell states remain incompletely understood. Here we employed histone mutants to uncover a crucial role for H3K36 methylation in the maintenance of cell identities across diverse developmental contexts. Focusing on the experimental induction of pluripotency, we show that H3K36M-mediated depletion of H3K36 methylation endows fibroblasts with a plastic state poised to acquire pluripotency in nearly all cells. At a cellular level, H3K36M facilitates epithelial plasticity by rendering fibroblasts insensitive to TGFß signals. At a molecular level, H3K36M enables the decommissioning of mesenchymal enhancers and the parallel activation of epithelial/stem cell enhancers. This enhancer rewiring is Tet dependent and redirects Sox2 from promiscuous somatic to pluripotency targets. Our findings reveal a previously unappreciated dual role for H3K36 methylation in the maintenance of cell identity by integrating a crucial developmental pathway into sustained expression of cell-type-specific programmes, and by opposing the expression of alternative lineage programmes through enhancer methylation.

Spatiotemporal transcriptomic maps of whole mouse embryos at the onset of organogenesis.

Spatiotemporal orchestration of gene expression is required for proper embryonic development. The use of single-cell technologies has begun to provide improved resolution of early regulatory dynamics, including detailed molecular definitions of most cell states during mouse embryogenesis. Here we used Slide-seq to build spatial transcriptomic maps of complete embryonic day (E) 8.5 and E9.0, and partial E9.5 embryos. To support their utility, we developed sc3D, a tool for reconstructing and exploring three-dimensional 'virtual embryos', which enables the quantitative investigation of regionalized gene expression patterns. Our measurements along the main embryonic axes of the developing neural tube revealed several previously unannotated genes with distinct spatial patterns. We also characterized the conflicting transcriptional identity of 'ectopic' neural tubes that emerge in Tbx6 mutant embryos. Taken together, we present an experimental and computational framework for the spatiotemporal investigation of whole embryonic structures and mutant phenotypes.

DNA hypomethylation silences anti-tumor immune genes in early prostate cancer and CTCs.

Cancer is characterized by hypomethylation-associated silencing of large chromatin domains, whose contribution to tumorigenesis is uncertain. Through high-resolution genome-wide single-cell DNA methylation sequencing, we identify 40 core domains that are uniformly hypomethylated from the earliest detectable stages of prostate malignancy through metastatic circulating tumor cells (CTCs). Nested among these repressive domains are smaller loci with preserved methylation that escape silencing and are enriched for cell proliferation genes. Transcriptionally silenced genes within the core hypomethylated domains are enriched for immune-related genes; prominent among these is a single gene cluster harboring all five CD1 genes that present lipid antigens to NKT cells and four IFI16-related interferon-inducible genes implicated in innate immunity. The re-expression of CD1 or IFI16 murine orthologs in immuno-competent mice abrogates tumorigenesis, accompanied by the activation of anti-tumor immunity. Thus, early epigenetic changes may shape tumorigenesis, targeting co-located genes within defined chromosomal loci. Hypomethylation domains are detectable in blood specimens enriched for CTCs.

Dynamic antagonism between key repressive pathways maintains the placental epigenome.

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.

T-REX17 is a transiently expressed non-coding RNA essential for human endoderm formation.

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.

JAK-STAT Signaling in Inflammatory Breast Cancer Enables Chemotherapy-Resistant Cell States.

Inflammatory breast cancer (IBC) is a difficult-to-treat disease with poor clinical outcomes due to high risk of metastasis and resistance to treatment. In breast cancer, CD44+CD24- cells possess stem cell-like features and contribute to disease progression, and we previously described a CD44+CD24-pSTAT3+ breast cancer cell subpopulation that is dependent on JAK2/STAT3 signaling. Here we report that CD44+CD24- cells are the most frequent cell type in IBC and are commonly pSTAT3+. Combination of JAK2/STAT3 inhibition with paclitaxel decreased IBC xenograft growth more than either agent alone. IBC cell lines resistant to paclitaxel and doxorubicin were developed and characterized to mimic therapeutic resistance in patients. Multi-omic profiling of parental and resistant cells revealed enrichment of genes associated with lineage identity and inflammation in chemotherapy-resistant derivatives. Integrated pSTAT3 chromatin immunoprecipitation sequencing and RNA sequencing (RNA-seq) analyses showed pSTAT3 regulates genes related to inflammation and epithelial-to-mesenchymal transition (EMT) in resistant cells, as well as PDE4A, a cAMP-specific phosphodiesterase. Metabolomic characterization identified elevated cAMP signaling and CREB as a candidate therapeutic target in IBC. Investigation of cellular dynamics and heterogeneity at the single cell level during chemotherapy and acquired resistance by CyTOF and single cell RNA-seq identified mechanisms of resistance including a shift from luminal to basal/mesenchymal cell states through selection for rare preexisting subpopulations or an acquired change. Finally, combination treatment with paclitaxel and JAK2/STAT3 inhibition prevented the emergence of the mesenchymal chemo-resistant subpopulation. These results provide mechanistic rational for combination of chemotherapy with inhibition of JAK2/STAT3 signaling as a more effective therapeutic strategy in IBC. SIGNIFICANCE: Chemotherapy resistance in inflammatory breast cancer is driven by the JAK2/STAT3 pathway, in part via cAMP/PKA signaling and a cell state switch, which can be overcome using paclitaxel combined with JAK2 inhibitors.

Sensitivity of Human Induced Pluripotent Stem Cells and Thereof Differentiated Kidney Proximal Tubular Cells towards Selected Nephrotoxins.

Proximal tubular epithelial cells (PTEC) are constantly exposed to potentially toxic metabolites and xenobiotics. The regenerative potential of the kidney enables the replacement of damaged cells either via the differentiation of stem cells or the re-acquisition of proliferative properties of the PTEC. Nevertheless, it is known that renal function declines, suggesting that the deteriorated cells are not replaced by fully functional cells. To understand the possible causes of this loss of kidney cell function, it is crucial to understand the role of toxins during the regeneration process. Therefore, we investigated the sensitivity and function of human induced pluripotent stem cells (hiPSC), hiPSC differentiating, and hiPSC differentiated into proximal tubular epithelial-like cells (PTELC) to known nephrotoxins. hiPSC were differentiated into PTELC, which exhibited similar morphology to PTEC, expressed prototypical PTEC markers, and were able to undergo albumin endocytosis. When treated with two nephrotoxins, hiPSC and differentiating hiPSC were more sensitive to cisplatin than differentiated PTELC, whereas all stages were equally sensitive to cyclosporin A. Both toxins also had an inhibitory effect on albumin uptake. Our results suggest a high sensitivity of differentiating cells towards toxins, which could have an unfavorable effect on regenerative processes. To study this, our model of hiPSC differentiating into PTELC appears suitable.

Repression and 3D-restructuring resolves regulatory conflicts in evolutionarily rearranged genomes.

Regulatory landscapes drive complex developmental gene expression, but it remains unclear how their integrity is maintained when incorporating novel genes and functions during evolution. Here, we investigated how a placental mammal-specific gene, Zfp42, emerged in an ancient vertebrate topologically associated domain (TAD) without adopting or disrupting the conserved expression of its gene, Fat1. In ESCs, physical TAD partitioning separates Zfp42 and Fat1 with distinct local enhancers that drive their independent expression. This separation is driven by chromatin activity and not CTCF/cohesin. In contrast, in embryonic limbs, inactive Zfp42 shares Fat1's intact TAD without responding to active Fat1 enhancers. However, neither Fat1 enhancer-incompatibility nor nuclear envelope-attachment account for Zfp42's unresponsiveness. Rather, Zfp42's promoter is rendered inert to enhancers by context-dependent DNA methylation. Thus, diverse mechanisms enabled the integration of independent Zfp42 regulation in the Fat1 locus. Critically, such regulatory complexity appears common in evolution as, genome wide, most TADs contain multiple independently expressed genes.

Molecular map of chronic lymphocytic leukemia and its impact on outcome.

Recent advances in cancer characterization have consistently revealed marked heterogeneity, impeding the completion of integrated molecular and clinical maps for each malignancy. Here, we focus on chronic lymphocytic leukemia (CLL), a B cell neoplasm with variable natural history that is conventionally categorized into two subtypes distinguished by extent of somatic mutations in the heavy-chain variable region of immunoglobulin genes (IGHV). To build the 'CLL map,' we integrated genomic, transcriptomic and epigenomic data from 1,148 patients. We identified 202 candidate genetic drivers of CLL (109 new) and refined the characterization of IGHV subtypes, which revealed distinct genomic landscapes and leukemogenic trajectories. Discovery of new gene expression subtypes further subcategorized this neoplasm and proved to be independent prognostic factors. Clinical outcomes were associated with a combination of genetic, epigenetic and gene expression features, further advancing our prognostic paradigm. Overall, this work reveals fresh insights into CLL oncogenesis and prognostication.

Temporally divergent regulatory mechanisms govern neuronal diversification and maturation in the mouse and marmoset neocortex.

Mammalian neocortical neurons span one of the most diverse cell type spectra of any tissue. Cortical neurons are born during embryonic development, and their maturation extends into postnatal life. The regulatory strategies underlying progressive neuronal development and maturation remain unclear. Here we present an integrated single-cell epigenomic and transcriptional analysis of individual mouse and marmoset cortical neuron classes, spanning both early postmitotic stages of identity acquisition and later stages of neuronal plasticity and circuit integration. We found that, in both species, the regulatory strategies controlling early and late stages of pan-neuronal development diverge. Early postmitotic neurons use more widely shared and evolutionarily conserved molecular regulatory programs. In contrast, programs active during later neuronal maturation are more brain- and neuron-specific and more evolutionarily divergent. Our work uncovers a temporal shift in regulatory choices during neuronal diversification and maturation in both mice and marmosets, which likely reflects unique evolutionary constraints on distinct events of neuronal development in the neocortex.

Hijacking of transcriptional condensates by endogenous retroviruses.

Most endogenous retroviruses (ERVs) in mammals are incapable of retrotransposition; therefore, why ERV derepression is associated with lethality during early development has been a mystery. Here, we report that rapid and selective degradation of the heterochromatin adapter protein TRIM28 triggers dissociation of transcriptional condensates from loci encoding super-enhancer (SE)-driven pluripotency genes and their association with transcribed ERV loci in murine embryonic stem cells. Knockdown of ERV RNAs or forced expression of SE-enriched transcription factors rescued condensate localization at SEs in TRIM28-degraded cells. In a biochemical reconstitution system, ERV RNA facilitated partitioning of RNA polymerase II and the Mediator coactivator into phase-separated droplets. In TRIM28 knockout mouse embryos, single-cell RNA-seq analysis revealed specific depletion of pluripotent lineages. We propose that coding and noncoding nascent RNAs, including those produced by retrotransposons, may facilitate 'hijacking' of transcriptional condensates in various developmental and disease contexts.