Impact of risk factors on early cancer evolution.

Recent identification of oncogenic cells within healthy tissues and the prevalence of indolent cancers found incidentally at autopsies reveal a greater complexity in tumor initiation than previously appreciated. The human body contains roughly 40 trillion cells of 200 different types that are organized within a complex three-dimensional matrix, necessitating exquisite mechanisms to restrain aberrant outgrowth of malignant cells that have the capacity to kill the host. Understanding how this defense is overcome to trigger tumorigenesis and why cancer is so extraordinarily rare at the cellular level is vital to future prevention therapies. In this review, we discuss how early initiated cells are protected from further tumorigenesis and the non-mutagenic pathways by which cancer risk factors promote tumor growth. By nature, the absence of permanent genomic alterations potentially renders these tumor-promoting mechanisms clinically targetable. Finally, we consider existing strategies for early cancer interception with perspectives on the next steps for molecular cancer prevention.

Deterministic Evolutionary Trajectories Influence Primary Tumor Growth: TRACERx Renal.

The evolutionary features of clear-cell renal cell carcinoma (ccRCC) have not been systematically studied to date. We analyzed 1,206 primary tumor regions from 101 patients recruited into the multi-center prospective study, TRACERx Renal. We observe up to 30 driver events per tumor and show that subclonal diversification is associated with known prognostic parameters. By resolving the patterns of driver event ordering, co-occurrence, and mutual exclusivity at clone level, we show the deterministic nature of clonal evolution. ccRCC can be grouped into seven evolutionary subtypes, ranging from tumors characterized by early fixation of multiple mutational and copy number drivers and rapid metastases to highly branched tumors with >10 subclonal drivers and extensive parallel evolution associated with attenuated progression. We identify genetic diversity and chromosomal complexity as determinants of patient outcome. Our insights reconcile the variable clinical behavior of ccRCC and suggest evolutionary potential as a biomarker for both intervention and surveillance.

Pervasive chromosomal instability and karyotype order in tumour evolution.

Chromosomal instability in cancer consists of dynamic changes to the number and structure of chromosomes1,2. The resulting diversity in somatic copy number alterations (SCNAs) may provide the variation necessary for tumour evolution1,3,4. Here we use multi-sample phasing and SCNA analysis of 1,421 samples from 394 tumours across 22 tumour types to show that continuous chromosomal instability results in pervasive SCNA heterogeneity. Parallel evolutionary events, which cause disruption in the same genes (such as BCL9, MCL1, ARNT (also known as HIF1B), TERT and MYC) within separate subclones, were present in 37% of tumours. Most recurrent losses probably occurred before whole-genome doubling, that was found as a clonal event in 49% of tumours. However, loss of heterozygosity at the human leukocyte antigen (HLA) locus and loss of chromosome 8p to a single haploid copy recurred at substantial subclonal frequencies, even in tumours with whole-genome doubling, indicating ongoing karyotype remodelling. Focal amplifications that affected chromosomes 1q21 (which encompasses BCL9, MCL1 and ARNT), 5p15.33 (TERT), 11q13.3 (CCND1), 19q12 (CCNE1) and 8q24.1 (MYC) were frequently subclonal yet appeared to be clonal within single samples. Analysis of an independent series of 1,024 metastatic samples revealed that 13 focal SCNAs were enriched in metastatic samples, including gains in chromosome 8q24.1 (encompassing MYC) in clear cell renal cell carcinoma and chromosome 11q13.3 (encompassing CCND1) in HER2+ breast cancer. Chromosomal instability may enable the continuous selection of SCNAs, which are established as ordered events that often occur in parallel, throughout tumour evolution.

Corrigendum: Phylogenetic ctDNA analysis depicts early-stage lung cancer evolution.

This corrects the article DOI: 10.1038/nature22364.

Phylogenetic ctDNA analysis depicts early-stage lung cancer evolution.

The early detection of relapse following primary surgery for non-small-cell lung cancer and the characterization of emerging subclones, which seed metastatic sites, might offer new therapeutic approaches for limiting tumour recurrence. The ability to track the evolutionary dynamics of early-stage lung cancer non-invasively in circulating tumour DNA (ctDNA) has not yet been demonstrated. Here we use a tumour-specific phylogenetic approach to profile the ctDNA of the first 100 TRACERx (Tracking Non-Small-Cell Lung Cancer Evolution Through Therapy (Rx)) study participants, including one patient who was also recruited to the PEACE (Posthumous Evaluation of Advanced Cancer Environment) post-mortem study. We identify independent predictors of ctDNA release and analyse the tumour-volume detection limit. Through blinded profiling of postoperative plasma, we observe evidence of adjuvant chemotherapy resistance and identify patients who are very likely to experience recurrence of their lung cancer. Finally, we show that phylogenetic ctDNA profiling tracks the subclonal nature of lung cancer relapse and metastasis, providing a new approach for ctDNA-driven therapeutic studies.

Replication stress links structural and numerical cancer chromosomal instability.

Cancer chromosomal instability (CIN) results in an increased rate of change of chromosome number and structure and generates intratumour heterogeneity. CIN is observed in most solid tumours and is associated with both poor prognosis and drug resistance. Understanding a mechanistic basis for CIN is therefore paramount. Here we find evidence for impaired replication fork progression and increased DNA replication stress in CIN(+) colorectal cancer (CRC) cells relative to CIN(-) CRC cells, with structural chromosome abnormalities precipitating chromosome missegregation in mitosis. We identify three new CIN-suppressor genes (PIGN (also known as MCD4), MEX3C (RKHD2) and ZNF516 (KIAA0222)) encoded on chromosome 18q that are subject to frequent copy number loss in CIN(+) CRC. Chromosome 18q loss was temporally associated with aneuploidy onset at the adenoma-carcinoma transition. CIN-suppressor gene silencing leads to DNA replication stress, structural chromosome abnormalities and chromosome missegregation. Supplementing cells with nucleosides, to alleviate replication-associated damage, reduces the frequency of chromosome segregation errors after CIN-suppressor gene silencing, and attenuates segregation errors and DNA damage in CIN(+) cells. These data implicate a central role for replication stress in the generation of structural and numerical CIN, which may inform new therapeutic approaches to limit intratumour heterogeneity.

Intratumor heterogeneity and branched evolution revealed by multiregion sequencing.

BACKGROUND: Intratumor heterogeneity may foster tumor evolution and adaptation and hinder personalized-medicine strategies that depend on results from single tumor-biopsy samples. METHODS: To examine intratumor heterogeneity, we performed exome sequencing, chromosome aberration analysis, and ploidy profiling on multiple spatially separated samples obtained from primary renal carcinomas and associated metastatic sites. We characterized the consequences of intratumor heterogeneity using immunohistochemical analysis, mutation functional analysis, and profiling of messenger RNA expression. RESULTS: Phylogenetic reconstruction revealed branched evolutionary tumor growth, with 63 to 69% of all somatic mutations not detectable across every tumor region. Intratumor heterogeneity was observed for a mutation within an autoinhibitory domain of the mammalian target of rapamycin (mTOR) kinase, correlating with S6 and 4EBP phosphorylation in vivo and constitutive activation of mTOR kinase activity in vitro. Mutational intratumor heterogeneity was seen for multiple tumor-suppressor genes converging on loss of function; SETD2, PTEN, and KDM5C underwent multiple distinct and spatially separated inactivating mutations within a single tumor, suggesting convergent phenotypic evolution. Gene-expression signatures of good and poor prognosis were detected in different regions of the same tumor. Allelic composition and ploidy profiling analysis revealed extensive intratumor heterogeneity, with 26 of 30 tumor samples from four tumors harboring divergent allelic-imbalance profiles and with ploidy heterogeneity in two of four tumors. CONCLUSIONS: Intratumor heterogeneity can lead to underestimation of the tumor genomics landscape portrayed from single tumor-biopsy samples and may present major challenges to personalized-medicine and biomarker development. Intratumor heterogeneity, associated with heterogeneous protein function, may foster tumor adaptation and therapeutic failure through Darwinian selection. (Funded by the Medical Research Council and others.).

Extrachromosomal DNA driven oncogene spatial heterogeneity and evolution in glioblastoma.

Oncogene amplification on extrachromosomal DNA (ecDNA) is strongly associated with treatment resistance and shorter survival for patients with cancer, including patients with glioblastoma. The non-chromosomal inheritance of ecDNA during cell division is a major contributor to intratumoral genetic heterogeneity. At present, the spatial dynamics of ecDNA, and the impact on tumor evolutionary trajectories, are not well understood. Here, we investigate the spatial-temporal evolution of ecDNA and its clinical impact by analyzing tumor samples from 94 treatment-naive human IDH -wildtype glioblastoma patients. We developed a spatial-temporal computational model of ecDNA positive tumors ('SPECIES') that integrates whole-genome sequencing, multi-region DNA FISH, and nascent RNAscope, to provide unique insight into the spatial dynamics of ecDNA evolution. Random segregation in combination with positive selection of ecDNAs induce large, predictable spatial patterns of cell-to-cell ecDNA copy number variation that are highly dependent on the oncogene encoded on the circular DNA. EGFR ecDNAs often reach high mean copy number (mean of 50 copies per tumor cell), are under strong positive selection (mean selection coefficient, s > 2) and do not co-amplify other oncogenes on the same ecDNA particles. In contrast, PDGFRA ecDNAs have lower mean copy number (mean of 15 copies per cell), are under weaker positive selection and frequently co-amplify other oncogenes on the same ecDNA. Evolutionary modeling suggests that EGFR ecDNAs often accumulate prior to clonal expansion. EGFR structural variants, including vIII and c-terminal deletions are under strong positive selection, are found exclusively on ecDNA, and are intermixed with wild-type EGFR ecDNAs. Simulations show EGFRvIII ecDNA likely arises after ecDNA formation in a cell with high wild-type EGFR copy number (> 10) before the onset of the most recent clonal expansion. This remains true even in cases of co-selection and co-amplification of multiple oncogenic ecDNA species in a subset of patients. Overall, our results suggest a potential time window in which early ecDNA detection may provide an opportunity for more effective intervention. Highlights: ecDNA is the most common mechanism of focal oncogene amplification in IDH wt glioblastoma. EGFR and its variants on ecDNA are particularly potent, likely arising early in tumor development, providing a strong oncogenic stimulus to drive tumorigenesis. Wild-type and variant EGFR ecDNA heteroplasmy (co-occurrence) is common with EGFR vIII or c-terminal deletions being derived from EGFR wild-type ecDNA prior to the most recent clonal expansion. Tumors with ecDNA amplified EGFR versus PDGFRA exhibit different evolutionary trajectories. SPECIES model can infer spatial evolutionary dynamics of ecDNA in cancer.A delay between ecDNA accumulation and subsequent oncogenic mutation may give a therapeutic window for early intervention.

Representation of genomic intratumor heterogeneity in multi-region non-small cell lung cancer patient-derived xenograft models.

Patient-derived xenograft (PDX) models are widely used in cancer research. To investigate the genomic fidelity of non-small cell lung cancer PDX models, we established 48 PDX models from 22 patients enrolled in the TRACERx study. Multi-region tumor sampling increased successful PDX engraftment and most models were histologically similar to their parent tumor. Whole-exome sequencing enabled comparison of tumors and PDX models and we provide an adapted mouse reference genome for improved removal of NOD scid gamma (NSG) mouse-derived reads from sequencing data. PDX model establishment caused a genomic bottleneck, with models often representing a single tumor subclone. While distinct tumor subclones were represented in independent models from the same tumor, individual PDX models did not fully recapitulate intratumor heterogeneity. On-going genomic evolution in mice contributed modestly to the genomic distance between tumors and PDX models. Our study highlights the importance of considering primary tumor heterogeneity when using PDX models and emphasizes the benefit of comprehensive tumor sampling.

The role of APOBEC3B in lung tumor evolution and targeted cancer therapy resistance.

In this study, the impact of the apolipoprotein B mRNA-editing catalytic subunit-like (APOBEC) enzyme APOBEC3B (A3B) on epidermal growth factor receptor (EGFR)-driven lung cancer was assessed. A3B expression in EGFR mutant (EGFRmut) non-small-cell lung cancer (NSCLC) mouse models constrained tumorigenesis, while A3B expression in tumors treated with EGFR-targeted cancer therapy was associated with treatment resistance. Analyses of human NSCLC models treated with EGFR-targeted therapy showed upregulation of A3B and revealed therapy-induced activation of nuclear factor kappa B (NF-κB) as an inducer of A3B expression. Significantly reduced viability was observed with A3B deficiency, and A3B was required for the enrichment of APOBEC mutation signatures, in targeted therapy-treated human NSCLC preclinical models. Upregulation of A3B was confirmed in patients with NSCLC treated with EGFR-targeted therapy. This study uncovers the multifaceted roles of A3B in NSCLC and identifies A3B as a potential target for more durable responses to targeted cancer therapy.

Mixed responses to targeted therapy driven by chromosomal instability through p53 dysfunction and genome doubling.

The phenomenon of mixed/heterogenous treatment responses to cancer therapies within an individual patient presents a challenging clinical scenario. Furthermore, the molecular basis of mixed intra-patient tumor responses remains unclear. Here, we show that patients with metastatic lung adenocarcinoma harbouring co-mutations of EGFR and TP53, are more likely to have mixed intra-patient tumor responses to EGFR tyrosine kinase inhibition (TKI), compared to those with an EGFR mutation alone. The combined presence of whole genome doubling (WGD) and TP53 co-mutations leads to increased genome instability and genomic copy number aberrations in genes implicated in EGFR TKI resistance. Using mouse models and an in vitro isogenic p53-mutant model system, we provide evidence that WGD provides diverse routes to drug resistance by increasing the probability of acquiring copy-number gains or losses relative to non-WGD cells. These data provide a molecular basis for mixed tumor responses to targeted therapy, within an individual patient, with implications for therapeutic strategies.

ALDH1L2 regulation of formate, formyl-methionine, and ROS controls cancer cell migration and metastasis.

Mitochondrial 10-formyltetrahydrofolate (10-formyl-THF) is utilized by three mitochondrial enzymes to produce formate for nucleotide synthesis, NADPH for antioxidant defense, and formyl-methionine (fMet) to initiate mitochondrial mRNA translation. One of these enzymes-aldehyde dehydrogenase 1 family member 2 (ALDH1L2)-produces NADPH by catabolizing 10-formyl-THF into CO2 and THF. Using breast cancer cell lines, we show that reduction of ALDH1L2 expression increases ROS levels and the production of both formate and fMet. Both depletion of ALDH1L2 and direct exposure to formate result in enhanced cancer cell migration that is dependent on the expression of the formyl-peptide receptor (FPR). In various tumor models, increased ALDH1L2 expression lowers formate and fMet accumulation and limits metastatic capacity, while human breast cancer samples show a consistent reduction of ALDH1L2 expression in metastases. Together, our data suggest that loss of ALDH1L2 can support metastatic progression by promoting formate and fMet production, resulting in enhanced FPR-dependent signaling.

CKS1 inhibition depletes leukemic stem cells and protects healthy hematopoietic stem cells in acute myeloid leukemia.

Acute myeloid leukemia (AML) is an aggressive hematological disorder comprising a hierarchy of quiescent leukemic stem cells (LSCs) and proliferating blasts with limited self-renewal ability. AML has a dismal prognosis, with extremely low 2-year survival rates in the poorest cytogenetic risk patients, primarily due to the failure of intensive chemotherapy protocols to deplete LSCs and toxicity of therapy toward healthy hematopoietic cells. We studied the role of cyclin-dependent kinase regulatory subunit 1 (CKS1)-dependent protein degradation in primary human AML and healthy hematopoiesis xenograft models in vivo. Using a small-molecule inhibitor (CKS1i), we demonstrate a dual role for CKS1-dependent protein degradation in reducing patient-derived AML blasts in vivo and, importantly, depleting LSCs, whereas inhibition of CKS1 has the opposite effect on normal hematopoiesis, protecting normal hematopoietic stem cells from chemotherapeutic toxicity. Proteomic analysis of responses to CKS1i in our patient-derived xenograft mouse model demonstrate that inhibition of CKS1 in AML leads to hyperactivation of RAC1 and accumulation of lethal reactive oxygen species, whereas healthy hematopoietic cells enter quiescence in response to CKS1i, protecting hematopoietic stem cells. Together, these findings demonstrate that CKS1-dependent proteostasis is a key vulnerability in malignant stem cell biology.

Progress towards non-small-cell lung cancer models that represent clinical evolutionary trajectories.

Non-small-cell lung cancer (NSCLC) is the leading cause of cancer-related deaths worldwide. Although advances are being made towards earlier detection and the development of impactful targeted therapies and immunotherapies, the 5-year survival of patients with advanced disease is still below 20%. Effective cancer research relies on pre-clinical model systems that accurately reflect the evolutionary course of disease progression and mimic patient responses to therapy. Here, we review pre-clinical models, including genetically engineered mouse models and patient-derived materials, such as cell lines, primary cell cultures, explant cultures and xenografts, that are currently being used to interrogate NSCLC evolution from pre-invasive disease through locally invasive cancer to the metastatic colonization of distant organ sites.

Induction of APOBEC3 Exacerbates DNA Replication Stress and Chromosomal Instability in Early Breast and Lung Cancer Evolution.

APOBEC3 enzymes are cytosine deaminases implicated in cancer. Precisely when APOBEC3 expression is induced during cancer development remains to be defined. Here we show that specific APOBEC3 genes are upregulated in breast ductal carcinoma in situ, and in preinvasive lung cancer lesions coincident with cellular proliferation. We observe evidence of APOBEC3-mediated subclonal mutagenesis propagated from TRACERx preinvasive to invasive non-small cell lung cancer (NSCLC) lesions. We find that APOBEC3B exacerbates DNA replication stress and chromosomal instability through incomplete replication of genomic DNA, manifested by accumulation of mitotic ultrafine bridges and 53BP1 nuclear bodies in the G1 phase of the cell cycle. Analysis of TRACERx NSCLC clinical samples and mouse lung cancer models revealed APOBEC3B expression driving replication stress and chromosome missegregation. We propose that APOBEC3 is functionally implicated in the onset of chromosomal instability and somatic mutational heterogeneity in preinvasive disease, providing fuel for selection early in cancer evolution. SIGNIFICANCE: This study reveals the dynamics and drivers of APOBEC3 gene expression in preinvasive disease and the exacerbation of cellular diversity by APOBEC3B through DNA replication stress to promote chromosomal instability early in cancer evolution.This article is highlighted in the In This Issue feature, p. 2355.

Cancer-Specific Loss of p53 Leads to a Modulation of Myeloid and T Cell Responses.

Loss of p53 function contributes to the development of many cancers. While cell-autonomous consequences of p53 mutation have been studied extensively, the role of p53 in regulating the anti-tumor immune response is still poorly understood. Here, we show that loss of p53 in cancer cells modulates the tumor-immune landscape to circumvent immune destruction. Deletion of p53 promotes the recruitment and instruction of suppressive myeloid CD11b+ cells, in part through increased expression of CXCR3/CCR2-associated chemokines and macrophage colony-stimulating factor (M-CSF), and attenuates the CD4+ T helper 1 (Th1) and CD8+ T cell responses in vivo. p53-null tumors also show an accumulation of suppressive regulatory T (Treg) cells. Finally, we show that two key drivers of tumorigenesis, activation of KRAS and deletion of p53, cooperate to promote immune tolerance.

Selective inhibition of cancer cell self-renewal through a Quisinostat-histone H1.0 axis.

Continuous cancer growth is driven by subsets of self-renewing malignant cells. Targeting of uncontrolled self-renewal through inhibition of stem cell-related signaling pathways has proven challenging. Here, we show that cancer cells can be selectively deprived of self-renewal ability by interfering with their epigenetic state. Re-expression of histone H1.0, a tumor-suppressive factor that inhibits cancer cell self-renewal in many cancer types, can be broadly induced by the clinically well-tolerated compound Quisinostat. Through H1.0, Quisinostat inhibits cancer cell self-renewal and halts tumor maintenance without affecting normal stem cell function. Quisinostat also hinders expansion of cells surviving targeted therapy, independently of the cancer types and the resistance mechanism, and inhibits disease relapse in mouse models of lung cancer. Our results identify H1.0 as a major mediator of Quisinostat's antitumor effect and suggest that sequential administration of targeted therapy and Quisinostat may be a broadly applicable strategy to induce a prolonged response in patients.

Tolerance of Chromosomal Instability in Cancer: Mechanisms and Therapeutic Opportunities.

Chromosomal instability (CIN) is the result of ongoing changes in the number (aneuploidy) and structure of chromosomes. CIN is induced by chromosome missegregation in mitosis and leads to karyotypic diversity within the cancer cell population, thereby adding to intratumor heterogeneity. Regardless of the overall pro-oncogenic function of CIN, its onset is typically detrimental for cell fitness and thus tumors must develop CIN-tolerance mechanisms in order to propagate. There is overwhelming genetic and functional evidence linking mutations in the tumor suppressor TP53 with CIN-tolerance. However, the pathways leading to p53 activation following chromosome missegregation remain controversial. Recently, additional mechanisms have been identified in CIN-surveillance, resulting in a more complex network of pathways acting independently or in cooperation with p53. Tolerance might also be achieved by modifying aspects of the cancer cell physiology in order to attenuate CIN or by adaptation to the consequences of aneuploid karyotypes. In this review, we summarize the current knowledge about p53-dependent and -independent mechanisms of CIN-tolerance in cancer, the adaptations observed in CIN cells buffering CIN levels, its consequences for cellular homeostasis, and the potential of exploiting these adaptations in order to design new cancer therapies.

Coactivator-dependent acetylation stabilizes members of the SREBP family of transcription factors.

Members of the SREBP family of transcription factors control cholesterol and lipid homeostasis and play important roles during adipocyte differentiation. The transcriptional activity of SREBPs is dependent on the coactivators p300 and CBP. We now present evidence that SREBPs are acetylated by the intrinsic acetyltransferase activity of p300 and CBP. In SREBP1a, the acetylated lysine residue resides in the DNA-binding domain of the protein. Coexpression with p300 dramatically increases the expression of both SREBP1a and SREBP2, and this effect is dependent on the acetyltransferase activity of p300, indicating that acetylation of SREBPs regulates their stability. Indeed, acetylation or mutation of the acetylated lysine residue in SREBP1a stabilizes the protein. We demonstrate that the acetylated residue in SREBP1a is also targeted by ubiquitination and that acetylation inhibits this process. Thus, our studies define acetylation-dependent stabilization of transcription factors as a novel mechanism for coactivators to regulate gene expression.

BCL9L Dysfunction Impairs Caspase-2 Expression Permitting Aneuploidy Tolerance in Colorectal Cancer.

Chromosomal instability (CIN) contributes to cancer evolution, intratumor heterogeneity, and drug resistance. CIN is driven by chromosome segregation errors and a tolerance phenotype that permits the propagation of aneuploid genomes. Through genomic analysis of colorectal cancers and cell lines, we find frequent loss of heterozygosity and mutations in BCL9L in aneuploid tumors. BCL9L deficiency promoted tolerance of chromosome missegregation events, propagation of aneuploidy, and genetic heterogeneity in xenograft models likely through modulation of Wnt signaling. We find that BCL9L dysfunction contributes to aneuploidy tolerance in both TP53-WT and mutant cells by reducing basal caspase-2 levels and preventing cleavage of MDM2 and BID. Efforts to exploit aneuploidy tolerance mechanisms and the BCL9L/caspase-2/BID axis may limit cancer diversity and evolution.