Long-term Multimodal Recording Reveals Epigenetic Adaptation Routes in Dormant Breast Cancer Cells.

Patients with estrogen receptor-positive breast cancer receive adjuvant endocrine therapies (ET) that delay relapse by targeting clinically undetectable micrometastatic deposits. Yet, up to 50% of patients relapse even decades after surgery through unknown mechanisms likely involving dormancy. To investigate genetic and transcriptional changes underlying tumor awakening, we analyzed late relapse patients and longitudinally profiled a rare cohort treated with long-term neoadjuvant ETs until progression. Next, we developed an in vitro evolutionary study to record the adaptive strategies of individual lineages in unperturbed parallel experiments. Our data demonstrate that ETs induce nongenetic cell state transitions into dormancy in a stochastic subset of cells via epigenetic reprogramming. Single lineages with divergent phenotypes awaken unpredictably in the absence of recurrent genetic alterations. Targeting the dormant epigenome shows promising activity against adapting cancer cells. Overall, this study uncovers the contribution of epigenetic adaptation to the evolution of resistance to ETs. SIGNIFICANCE: This study advances the understanding of therapy-induced dormancy with potential clinical implications for breast cancer. Estrogen receptor-positive breast cancer cells adapt to endocrine treatment by entering a dormant state characterized by strong heterochromatinization with no recurrent genetic changes. Targeting the epigenetic rewiring impairs the adaptation of cancer cells to ETs. See related commentary by Llinas-Bertran et al., p. 704. This article is featured in Selected Articles from This Issue, p. 695.

Epigenome and early selection determine the tumour-immune evolutionary trajectory of colorectal cancer.

Immune system control is a major hurdle that cancer evolution must circumvent. The relative timing and evolutionary dynamics of subclones that have escaped immune control remain incompletely characterized, and how immune-mediated selection shapes the epigenome has received little attention. Here, we infer the genome- and epigenome-driven evolutionary dynamics of tumour-immune coevolution within primary colorectal cancers (CRCs). We utilise our existing CRC multi-region multi-omic dataset that we supplement with high-resolution spatially-resolved neoantigen sequencing data and highly multiplexed imaging of the tumour microenvironment (TME). Analysis of somatic chromatin accessibility alterations (SCAAs) reveals frequent somatic loss of accessibility at antigen presenting genes, and that SCAAs contribute to silencing of neoantigens. We observe that strong immune escape and exclusion occur at the outset of CRC formation, and that within tumours, including at the microscopic level of individual tumour glands, additional immune escape alterations have negligible consequences for the immunophenotype of cancer cells. Further minor immuno-editing occurs during local invasion and is associated with TME reorganisation, but that evolutionary bottleneck is relatively weak. Collectively, we show that immune evasion in CRC follows a "Big Bang" evolutionary pattern, whereby genetic, epigenetic and TME-driven immune evasion acquired by the time of transformation defines subsequent cancer-immune evolution.

The co-evolution of the genome and epigenome in colorectal cancer.

Colorectal malignancies are a leading cause of cancer-related death1 and have undergone extensive genomic study2,3. However, DNA mutations alone do not fully explain malignant transformation4-7. Here we investigate the co-evolution of the genome and epigenome of colorectal tumours at single-clone resolution using spatial multi-omic profiling of individual glands. We collected 1,370 samples from 30 primary cancers and 8 concomitant adenomas and generated 1,207 chromatin accessibility profiles, 527 whole genomes and 297 whole transcriptomes. We found positive selection for DNA mutations in chromatin modifier genes and recurrent somatic chromatin accessibility alterations, including in regulatory regions of cancer driver genes that were otherwise devoid of genetic mutations. Genome-wide alterations in accessibility for transcription factor binding involved CTCF, downregulation of interferon and increased accessibility for SOX and HOX transcription factor families, suggesting the involvement of developmental genes during tumourigenesis. Somatic chromatin accessibility alterations were heritable and distinguished adenomas from cancers. Mutational signature analysis showed that the epigenome in turn influences the accumulation of DNA mutations. This study provides a map of genetic and epigenetic tumour heterogeneity, with fundamental implications for understanding colorectal cancer biology.

Phenotypic plasticity and genetic control in colorectal cancer evolution.

Genetic and epigenetic variation, together with transcriptional plasticity, contribute to intratumour heterogeneity1. The interplay of these biological processes and their respective contributions to tumour evolution remain unknown. Here we show that intratumour genetic ancestry only infrequently affects gene expression traits and subclonal evolution in colorectal cancer (CRC). Using spatially resolved paired whole-genome and transcriptome sequencing, we find that the majority of intratumour variation in gene expression is not strongly heritable but rather 'plastic'. Somatic expression quantitative trait loci analysis identified a number of putative genetic controls of expression by cis-acting coding and non-coding mutations, the majority of which were clonal within a tumour, alongside frequent structural alterations. Consistently, computational inference on the spatial patterning of tumour phylogenies finds that a considerable proportion of CRCs did not show evidence of subclonal selection, with only a subset of putative genetic drivers associated with subclone expansions. Spatial intermixing of clones is common, with some tumours growing exponentially and others only at the periphery. Together, our data suggest that most genetic intratumour variation in CRC has no major phenotypic consequence and that transcriptional plasticity is, instead, widespread within a tumour.

Circulating tumour DNA sequencing to determine therapeutic response and identify tumour heterogeneity in patients with paediatric solid tumours.

OBJECTIVE: Clinical diagnostic sequencing of circulating tumour DNA (ctDNA) is well advanced for adult patients, but application to paediatric cancer patients lags behind. METHODS: To address this, we have developed a clinically relevant (67 gene) NGS capture panel and accompanying workflow that enables sensitive and reliable detection of low-frequency genetic variants in cell-free DNA (cfDNA) from children with solid tumours. We combined gene panel sequencing with low pass whole-genome sequencing of the same library to inform on genome-wide copy number changes in the blood. RESULTS: Analytical validity was evaluated using control materials, and the method was found to be highly sensitive (0.96 for SNVs and 0.97 for INDEL), specific (0.82 for SNVs and 0.978 for INDEL), repeatable (>0.93 [95% CI: 0.89-0.95]) and reproducible (>0.87 [95% CI: 0.87-0.95]). Potential for clinical application was demonstrated in 39 childhood cancer patients with a spectrum of solid tumours in which the single nucleotide variants expected from tumour sequencing were detected in cfDNA in 94.4% (17/18) of cases with active extracranial disease. In 13 patients, where serial samples were available, we show a close correlation between events detected in cfDNA and treatment response, demonstrate that cfDNA analysis could be a useful tool to monitor disease progression, and show cfDNA sequencing has the potential to identify targetable variants that were not detected in tumour samples. CONCLUSIONS: This is the first pan-cancer DNA sequencing panel that we know to be optimised for cfDNA in children for blood-based molecular diagnostics in paediatric solid tumours.

Functional reconstruction of human AML reveals stem cell origin and vulnerability of treatment-resistant MLL-rearranged leukemia.

Chemoresistance remains the major challenge for successful treatment of acute myeloid leukemia (AML). Although recent mouse studies suggest that treatment response of genetically and immunophenotypically indistinguishable AML can be influenced by their different cells of origin, corresponding evidence in human disease is still largely lacking. By combining prospective disease modeling using highly purified human hematopoietic stem or progenitor cells with retrospective deconvolution study of leukemia stem cells (LSCs) from primary patient samples, we identified human hematopoietic stem cells (HSCs) and common myeloid progenitors (CMPs) as two distinctive origins of human AML driven by Mixed Lineage Leukemia (MLL) gene fusions (MLL-AML). Despite LSCs from either MLL-rearranged HSCs or MLL-rearranged CMPs having a mature CD34-/lo/CD38+ immunophenotype in both a humanized mouse model and primary patient samples, the resulting AML cells exhibited contrasting responses to chemotherapy. HSC-derived MLL-AML was highly resistant to chemotherapy and expressed elevated amounts of the multispecific anion transporter ABCC3. Inhibition of ABCC3 by shRNA-mediated knockdown or with small-molecule inhibitor fidaxomicin, currently used for diarrhea associated with Clostridium difficile infection, effectively resensitized HSC-derived MLL-AML toward standard chemotherapeutic drugs. This study not only functionally established two distinctive origins of human LSCs for MLL-AML and their role in mediating chemoresistance but also identified a potential therapeutic avenue for stem cell-associated treatment resistance by repurposing a well-tolerated antidiarrhea drug already used in the clinic.

Transcriptional memory of cells of origin overrides ß-catenin requirement of MLL cancer stem cells.

While ß-catenin has been demonstrated as an essential molecule and therapeutic target for various cancer stem cells (CSCs) including those driven by MLL fusions, here we show that transcriptional memory from cells of origin predicts AML patient survival and allows ß-catenin-independent transformation in MLL-CSCs derived from hematopoietic stem cell (HSC)-enriched LSK population but not myeloid-granulocyte progenitors. Mechanistically, ß-catenin regulates expression of downstream targets of a key transcriptional memory gene, Hoxa9 that is highly enriched in LSK-derived MLL-CSCs and helps sustain leukemic self-renewal. Suppression of Hoxa9 sensitizes LSK-derived MLL-CSCs to ß-catenin inhibition resulting in abolishment of CSC transcriptional program and transformation ability. In addition, further molecular and functional analyses identified Prmt1 as a key common downstream mediator for ß-catenin/Hoxa9 functions in LSK-derived MLL-CSCs. Together, these findings not only uncover an unexpectedly important role of cells of origin transcriptional memory in regulating CSC self-renewal, but also reveal a novel molecular network mediated by ß-catenin/Hoxa9/Prmt1 in governing leukemic self-renewal.

Molecular synergy underlies the co-occurrence patterns and phenotype of NPM1-mutant acute myeloid leukemia.

NPM1 mutations define the commonest subgroup of acute myeloid leukemia (AML) and frequently co-occur with FLT3 internal tandem duplications (ITD) or, less commonly, NRAS or KRAS mutations. Co-occurrence of mutant NPM1 with FLT3-ITD carries a significantly worse prognosis than NPM1-RAS combinations. To understand the molecular basis of these observations, we compare the effects of the 2 combinations on hematopoiesis and leukemogenesis in knock-in mice. Early effects of these mutations on hematopoiesis show that compound Npm1cA/+;NrasG12D/+ or Npm1cA;Flt3ITD share a number of features: Hox gene overexpression, enhanced self-renewal, expansion of hematopoietic progenitors, and myeloid differentiation bias. However, Npm1cA;Flt3ITD mutants displayed significantly higher peripheral leukocyte counts, early depletion of common lymphoid progenitors, and a monocytic bias in comparison with the granulocytic bias in Npm1cA/+;NrasG12D/+ mutants. Underlying this was a striking molecular synergy manifested as a dramatically altered gene expression profile in Npm1cA;Flt3ITD , but not Npm1cA/+;NrasG12D/+ , progenitors compared with wild-type. Both double-mutant models developed high-penetrance AML, although latency was significantly longer with Npm1cA/+;NrasG12D/+ During AML evolution, both models acquired additional copies of the mutant Flt3 or Nras alleles, but only Npm1cA/+;NrasG12D/+ mice showed acquisition of other human AML mutations, including IDH1 R132Q. We also find, using primary Cas9-expressing AMLs, that Hoxa genes and selected interactors or downstream targets are required for survival of both types of double-mutant AML. Our results show that molecular complementarity underlies the higher frequency and significantly worse prognosis associated with NPM1c/FLT3-ITD vs NPM1/NRAS-G12D-mutant AML and functionally confirm the role of HOXA genes in NPM1c-driven AML.

Novel Denisovan and Neanderthal retroviruses.

Following the recent availability of high-coverage genomes for Denisovan and Neanderthal hominids, we conducted a screen for endogenized retroviruses, identifying six novel, previously unreported HERV-K(HML2) elements (HERV-K is human endogenous retrovirus K). These elements are absent from the human genome (hg38) and appear to be unique to archaic hominids. These findings provide further evidence supporting the recent activity of the HERV-K(HML2) group, which has been implicated in human disease. They will also provide insights into the evolution of archaic hominids.