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.

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.

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.

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.

Acute lymphoblastic leukemia displays a distinct highly methylated genome.

DNA methylation is tightly regulated during development and is stably maintained in healthy cells. In contrast, cancer cells are commonly characterized by a global loss of DNA methylation co-occurring with CpG island hypermethylation. In acute lymphoblastic leukemia (ALL), the commonest childhood cancer, perturbations of CpG methylation have been reported to be associated with genetic disease subtype and outcome, but data from large cohorts at a genome-wide scale are lacking. Here, we performed whole-genome bisulfite sequencing across ALL subtypes, leukemia cell lines and healthy hematopoietic cells, and show that unlike most cancers, ALL samples exhibit CpG island hypermethylation but minimal global loss of methylation. This was most pronounced in T cell ALL and accompanied by an exceptionally broad range of hypermethylation of CpG islands between patients, which is influenced by TET2 and DNMT3B. These findings demonstrate that ALL is characterized by an unusually highly methylated genome and provide further insights into the non-canonical regulation of methylation in cancer.

Systematic Profiling of DNMT3A Variants Reveals Protein Instability Mediated by the DCAF8 E3 Ubiquitin Ligase Adaptor.

Clonal hematopoiesis is a prevalent age-related condition associated with a greatly increased risk of hematologic disease; mutations in DNA methyltransferase 3A (DNMT3A) are the most common driver of this state. DNMT3A variants occur across the gene with some particularly associated with malignancy, but the functional relevance and mechanisms of pathogenesis of the majority of mutations are unknown. Here, we systematically investigated the methyltransferase activity and protein stability of 253 disease-associated DNMT3A mutations, and found that 74% were loss-of-function mutations. Half of these variants exhibited reduced protein stability and, as a class, correlated with greater clonal expansion and acute myeloid leukemia development. We investigated the mechanisms underlying the instability using a CRISPR screen and uncovered regulated destruction of DNMT3A mediated by the DCAF8 E3 ubiquitin ligase adaptor. We establish a new paradigm to classify novel variants that has prognostic and potential therapeutic significance for patients with hematologic disease. SIGNIFICANCE: DNMT3A has emerged as the most important epigenetic regulator and tumor suppressor in the hematopoietic system. Our study represents a systematic and high-throughput method to characterize the molecular impact of DNMT3A missense mutations and the discovery of a regulated destruction mechanism of DNMT3A offering new prognostic and future therapeutic avenues.See related commentary by Ma and Will, p. 23.This article is highlighted in the In This Issue feature, p. 1.

Chromatin-state barriers enforce an irreversible mammalian cell fate decision.

Stem and progenitor cells have the capacity to balance self-renewal and differentiation. Hematopoietic myeloid progenitors replenish more than 25 billion terminally differentiated neutrophils every day under homeostatic conditions and can increase this output in response to stress or infection. At what point along the spectrum of maturation do progenitors lose capacity for self-renewal and become irreversibly committed to differentiation? Using a system of conditional myeloid development that can be toggled between self-renewal and differentiation, we interrogate determinants of this "point of no return" in differentiation commitment. Irreversible commitment is due primarily to loss of open regulatory site access and disruption of a positive feedback transcription factor activation loop. Restoration of the transcription factor feedback loop extends the window of cell plasticity and alters the point of no return. These findings demonstrate how the chromatin state enforces and perpetuates cell fate and identify potential avenues for manipulating cell identity.

Activation of Notch and Myc Signaling via B-cell-Restricted Depletion of Dnmt3a Generates a Consistent Murine Model of Chronic Lymphocytic Leukemia.

Chronic lymphocytic leukemia (CLL) is characterized by disordered DNA methylation, suggesting these epigenetic changes might play a critical role in disease onset and progression. The methyltransferase DNMT3A is a key regulator of DNA methylation. Although DNMT3A somatic mutations in CLL are rare, we found that low DNMT3A expression is associated with more aggressive disease. A conditional knockout mouse model showed that homozygous depletion of Dnmt3a from B cells results in the development of CLL with 100% penetrance at a median age of onset of 5.3 months, and heterozygous Dnmt3a depletion yields a disease penetrance of 89% with a median onset at 18.5 months, confirming its role as a haploinsufficient tumor suppressor. B1a cells were confirmed as the cell of origin of disease in this model, and Dnmt3a depletion resulted in focal hypomethylation and activation of Notch and Myc signaling. Amplification of chromosome 15 containing the Myc gene was detected in all CLL mice tested, and infiltration of high-Myc-expressing CLL cells in the spleen was observed. Notably, hyperactivation of Notch and Myc signaling was exclusively observed in the Dnmt3a CLL mice, but not in three other CLL mouse models tested (Sf3b1-Atm, Ikzf3, and MDR), and Dnmt3a-depleted CLL were sensitive to pharmacologic inhibition of Notch signaling in vitro and in vivo. Consistent with these findings, human CLL samples with lower DNMT3A expression were more sensitive to Notch inhibition than those with higher DNMT3A expression. Altogether, these results suggest that Dnmt3a depletion induces CLL that is highly dependent on activation of Notch and Myc signaling. SIGNIFICANCE: Loss of DNMT3A expression is a driving event in CLL and is associated with aggressive disease, activation of Notch and Myc signaling, and enhanced sensitivity to Notch inhibition.

Smart-RRBS for single-cell methylome and transcriptome analysis.

The integration of DNA methylation and transcriptional state within single cells is of broad interest. Several single-cell dual- and multi-omics approaches have been reported that enable further investigation into cellular heterogeneity, including the discovery and in-depth study of rare cell populations. Such analyses will continue to provide important mechanistic insights into the regulatory consequences of epigenetic modifications. We recently reported a new method for profiling the DNA methylome and transcriptome from the same single cells in a cancer research study. Here, we present details of the protocol and provide guidance on its utility. Our Smart-RRBS (reduced representation bisulfite sequencing) protocol combines Smart-seq2 and RRBS and entails physically separating mRNA from the genomic DNA. It generates paired epigenetic promoter and RNA-expression measurements for ~24% of protein-coding genes in a typical single cell. It also works for micro-dissected tissue samples comprising hundreds of cells. The protocol, excluding flow sorting of cells and sequencing, takes ~3 d to process up to 192 samples manually. It requires basic molecular biology expertise and laboratory equipment, including a PCR workstation with UV sterilization, a DNA fluorometer and a microfluidic electrophoresis system.

Discovery of Candidate DNA Methylation Cancer Driver Genes.

Epigenetic alterations, such as promoter hypermethylation, may drive cancer through tumor suppressor gene inactivation. However, we have limited ability to differentiate driver DNA methylation (DNAme) changes from passenger events. We developed DNAme driver inference-MethSig-accounting for the varying stochastic hypermethylation rate across the genome and between samples. We applied MethSig to bisulfite sequencing data of chronic lymphocytic leukemia (CLL), multiple myeloma, ductal carcinoma in situ, glioblastoma, and to methylation array data across 18 tumor types in TCGA. MethSig resulted in well-calibrated quantile-quantile plots and reproducible inference of likely DNAme drivers with increased sensitivity/specificity compared with benchmarked methods. CRISPR/Cas9 knockout of selected candidate CLL DNAme drivers provided a fitness advantage with and without therapeutic intervention. Notably, DNAme driver risk score was closely associated with adverse outcome in independent CLL cohorts. Collectively, MethSig represents a novel inference framework for DNAme driver discovery to chart the role of aberrant DNAme in cancer. SIGNIFICANCE: MethSig provides a novel statistical framework for the analysis of DNA methylation changes in cancer, to specifically identify candidate DNA methylation driver genes of cancer progression and relapse, empowering the discovery of epigenetic mechanisms that enhance cancer cell fitness.This article is highlighted in the In This Issue feature, p. 2113.

Preneoplastic Alterations Define CLL DNA Methylome and Persist through Disease Progression and Therapy.

Most human cancers converge to a deregulated methylome with reduced global levels and elevated methylation at select CpG islands. To investigate the emergence and dynamics of the cancer methylome, we characterized genome-wide DNA methylation in pre-neoplastic monoclonal B cell lymphocytosis (MBL) and chronic lymphocytic leukemia (CLL), including serial samples collected across disease course. We detected the aberrant tumor-associated methylation landscape at CLL diagnosis and found no significantly differentially methylated regions in the high-count MBL-to-CLL transition. Patient methylomes showed remarkable stability with natural disease and post-therapy progression. Single CLL cells were consistently aberrantly methylated, indicating a homogeneous transition to the altered epigenetic state, and a distinct expression profile together with MBL cells compared to normal B cells. Our longitudinal analysis reveals the cancer methylome to emerge early, which may provide a platform for subsequent genetically-driven growth dynamics and together with its persistent presence suggests a central role in the normal-to-cancer transition.

TETs compete with DNMT3 activity in pluripotent cells at thousands of methylated somatic enhancers.

Mammalian cells stably maintain high levels of DNA methylation despite expressing both positive (DNMT3A/B) and negative (TET1-3) regulators. Here, we analyzed the independent and combined effects of these regulators on the DNA methylation landscape using a panel of knockout human embryonic stem cell (ESC) lines. The greatest impact on global methylation levels was observed in DNMT3-deficient cells, including reproducible focal demethylation at thousands of normally methylated loci. Demethylation depends on TET expression and occurs only when both DNMT3s are absent. Dynamic loci are enriched for hydroxymethylcytosine and overlap with subsets of putative somatic enhancers that are methylated in ESCs and can be activated upon differentiation. We observe similar dynamics in mouse ESCs that were less frequent in epiblast stem cells (EpiSCs) and scarce in somatic tissues, suggesting a conserved pluripotency-linked mechanism. Taken together, our data reveal tightly regulated competition between DNMT3s and TETs at thousands of somatic regulatory sequences within pluripotent cells.

Netrin-1 promotes naive pluripotency through Neo1 and Unc5b co-regulation of Wnt and MAPK signalling.

In mouse embryonic stem cells (mESCs), chemical blockade of Gsk3α/ß and Mek1/2 (2i) instructs a self-renewing ground state whose endogenous inducers are unknown. Here we show that the axon guidance cue Netrin-1 promotes naive pluripotency by triggering profound signalling, transcriptomic and epigenetic changes in mESCs. Furthermore, we demonstrate that Netrin-1 can substitute for blockade of Gsk3α/ß and Mek1/2 to sustain self-renewal of mESCs in combination with leukaemia inhibitory factor and regulates the formation of the mouse pluripotent blastocyst. Mechanistically, we reveal how Netrin-1 and the balance of its receptors Neo1 and Unc5B co-regulate Wnt and MAPK pathways in both mouse and human ESCs. Netrin-1 induces Fak kinase to inactivate Gsk3α/ß and stabilize ß-catenin while increasing the phosphatase activity of a Ppp2r2c-containing Pp2a complex to reduce Erk1/2 activity. Collectively, this work identifies Netrin-1 as a regulator of pluripotency and reveals that it mediates different effects in mESCs depending on its receptor dosage, opening perspectives for balancing self-renewal and lineage commitment.

Circadian Entrainment Triggers Maturation of Human In Vitro Islets.

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.

Author Correction: LKB1 loss links serine metabolism to DNA methylation and tumorigenesis.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Inducible histone K-to-M mutations are dynamic tools to probe the physiological role of site-specific histone methylation in vitro and in vivo.

Development and differentiation are associated with profound changes to histone modifications, yet their in vivo function remains incompletely understood. Here, we generated mouse models expressing inducible histone H3 lysine-to-methionine (K-to-M) mutants, which globally inhibit methylation at specific sites. Mice expressing H3K36M developed severe anaemia with arrested erythropoiesis, a marked haematopoietic stem cell defect, and rapid lethality. By contrast, mice expressing H3K9M survived up to a year and showed expansion of multipotent progenitors, aberrant lymphopoiesis and thrombocytosis. Additionally, some H3K9M mice succumbed to aggressive T cell leukaemia/lymphoma, while H3K36M mice exhibited differentiation defects in testis and intestine. Mechanistically, induction of either mutant reduced corresponding histone trimethylation patterns genome-wide and altered chromatin accessibility as well as gene expression landscapes. Strikingly, discontinuation of transgene expression largely restored differentiation programmes. Our work shows that individual chromatin modifications are required at several specific stages of differentiation and introduces powerful tools to interrogate their roles in vivo.

The RNA Helicase DDX6 Controls Cellular Plasticity by Modulating P-Body Homeostasis.

Post-transcriptional mechanisms have the potential to influence complex changes in gene expression, yet their role in cell fate transitions remains largely unexplored. Here, we show that suppression of the RNA helicase DDX6 endows human and mouse primed embryonic stem cells (ESCs) with a differentiation-resistant, "hyper-pluripotent" state, which readily reprograms to a naive state resembling the preimplantation embryo. We further demonstrate that DDX6 plays a key role in adult progenitors where it controls the balance between self-renewal and differentiation in a context-dependent manner. Mechanistically, DDX6 mediates the translational suppression of target mRNAs in P-bodies. Upon loss of DDX6 activity, P-bodies dissolve and release mRNAs encoding fate-instructive transcription and chromatin factors that re-enter the ribosome pool. Increased translation of these targets impacts cell fate by rewiring the enhancer, heterochromatin, and DNA methylation landscapes of undifferentiated cell types. Collectively, our data establish a link between P-body homeostasis, chromatin organization, and stem cell potency.

Epigenetic evolution and lineage histories of chronic lymphocytic leukaemia.

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.

Corrupted coordination of epigenetic modifications leads to diverging chromatin states and transcriptional heterogeneity in CLL.

Cancer evolution is fueled by epigenetic as well as genetic diversity. In chronic lymphocytic leukemia (CLL), intra-tumoral DNA methylation (DNAme) heterogeneity empowers evolution. Here, to comprehensively study the epigenetic dimension of cancer evolution, we integrate DNAme analysis with histone modification mapping and single cell analyses of RNA expression and DNAme in 22 primary CLL and 13 healthy donor B lymphocyte samples. Our data reveal corrupted coherence across different layers of the CLL epigenome. This manifests in decreased mutual information across epigenetic modifications and gene expression attributed to cell-to-cell heterogeneity. Disrupted epigenetic-transcriptional coordination in CLL is also reflected in the dysregulation of the transcriptional output as a function of the combinatorial chromatin states, including incomplete Polycomb-mediated gene silencing. Notably, we observe unexpected co-mapping of typically mutually exclusive activating and repressing histone modifications, suggestive of intra-tumoral epigenetic diversity. Thus, CLL epigenetic diversification leads to decreased coordination across layers of epigenetic information, likely reflecting an admixture of cells with diverging cellular identities.