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.

Author Correction: Resolving genetic heterogeneity in cancer.

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

Resolving genetic heterogeneity in cancer.

To a large extent, cancer conforms to evolutionary rules defined by the rates at which clones mutate, adapt and grow. Next-generation sequencing has provided a snapshot of the genetic landscape of most cancer types, and cancer genomics approaches are driving new insights into cancer evolutionary patterns in time and space. In contrast to species evolution, cancer is a particular case owing to the vast size of tumour cell populations, chromosomal instability and its potential for phenotypic plasticity. Nevertheless, an evolutionary framework is a powerful aid to understand cancer progression and therapy failure. Indeed, such a framework could be applied to predict individual tumour behaviour and support treatment strategies.

Measuring Clonal Evolution in Cancer with Genomics.

Cancers originate from somatic cells in the human body that have accumulated genetic alterations. These mutations modify the phenotype of the cells, allowing them to escape the homeostatic regulation that maintains normal cell number. Viewed through the lens of evolutionary biology, the transformation of normal cells into malignant cells is evolution in action. Evolution continues throughout cancer growth, progression, treatment resistance, and disease relapse, driven by adaptation to changes in the cancer's environment, and intratumor heterogeneity is an inevitable consequence of this evolutionary process. Genomics provides a powerful means to characterize tumor evolution, enabling quantitative measurement of evolving clones across space and time. In this review, we discuss concepts and approaches to quantify and measure this evolutionary process in cancer using genomics.

EGFR amplification and outcome in a randomised phase III trial of chemotherapy alone or chemotherapy plus panitumumab for advanced gastro-oesophageal cancers.

OBJECTIVE: Epidermal growth factor receptor (EGFR) inhibition may be effective in biomarker-selected populations of advanced gastro-oesophageal adenocarcinoma (aGEA) patients. Here, we tested the association between outcome and EGFR copy number (CN) in pretreatment tissue and plasma cell-free DNA (cfDNA) of patients enrolled in a randomised first-line phase III clinical trial of chemotherapy or chemotherapy plus the anti-EGFR monoclonal antibody panitumumab in aGEA (NCT00824785). DESIGN: EGFR CN by either fluorescence in situ hybridisation (n=114) or digital-droplet PCR in tissues (n=250) and plasma cfDNAs (n=354) was available for 474 (86%) patients in the intention-to-treat (ITT) population. Tissue and plasma low-pass whole-genome sequencing was used to screen for coamplifications in receptor tyrosine kinases. Interaction between chemotherapy and EGFR inhibitors was modelled in patient-derived organoids (PDOs) from aGEA patients. RESULTS: EGFR amplification in cfDNA correlated with poor survival in the ITT population and similar trends were observed when the analysis was conducted in tissue and plasma by treatment arm. EGFR inhibition in combination with chemotherapy did not correlate with improved survival, even in patients with significant EGFR CN gains. Addition of anti-EGFR inhibitors to the chemotherapy agent epirubicin in PDOs, resulted in a paradoxical increase in viability and accelerated progression through the cell cycle, associated with p21 and cyclin B1 downregulation and cyclin E1 upregulation, selectively in organoids from EGFR-amplified aGEA. CONCLUSION: EGFR CN can be accurately measured in tissue and liquid biopsies and may be used for the selection of aGEA patients. EGFR inhibitors may antagonise the antitumour effect of anthracyclines with important implications for the design of future combinatorial trials.

Detecting repeated cancer evolution from multi-region tumor sequencing data.

Recurrent successions of genomic changes, both within and between patients, reflect repeated evolutionary processes that are valuable for the anticipation of cancer progression. Multi-region sequencing allows the temporal order of some genomic changes in a tumor to be inferred, but the robust identification of repeated evolution across patients remains a challenge. We developed a machine-learning method based on transfer learning that allowed us to overcome the stochastic effects of cancer evolution and noise in data and identified hidden evolutionary patterns in cancer cohorts. When applied to multi-region sequencing datasets from lung, breast, renal, and colorectal cancer (768 samples from 178 patients), our method detected repeated evolutionary trajectories in subgroups of patients, which were reproduced in single-sample cohorts (n = 2,935). Our method provides a means of classifying patients on the basis of how their tumor evolved, with implications for the anticipation of disease progression.

The MOBSTER R package for tumour subclonal deconvolution from bulk DNA whole-genome sequencing data.

BACKGROUND: The large-scale availability of whole-genome sequencing profiles from bulk DNA sequencing of cancer tissues is fueling the application of evolutionary theory to cancer. From a bulk biopsy, subclonal deconvolution methods are used to determine the composition of cancer subpopulations in the biopsy sample, a fundamental step to determine clonal expansions and their evolutionary trajectories. RESULTS: In a recent work we have developed a new model-based approach to carry out subclonal deconvolution from the site frequency spectrum of somatic mutations. This new method integrates, for the first time, an explicit model for neutral evolutionary forces that participate in clonal expansions; in that work we have also shown that our method improves largely over competing data-driven methods. In this Software paper we present mobster, an open source R package built around our new deconvolution approach, which provides several functions to plot data and fit models, assess their confidence and compute further evolutionary analyses that relate to subclonal deconvolution. CONCLUSIONS: We present the mobster package for tumour subclonal deconvolution from bulk sequencing, the first approach to integrate Machine Learning and Population Genetics which can explicitly model co-existing neutral and positive selection in cancer. We showcase the analysis of two datasets, one simulated and one from a breast cancer patient, and overview all package functionalities.

Prediction of Benefit from Checkpoint Inhibitors in Mismatch Repair Deficient Metastatic Colorectal Cancer: Role of Tumor Infiltrating Lymphocytes.

BACKGROUND: Immunotherapy with immune checkpoint inhibitors (ICIs) is highly effective in microsatellite instability-high (MSI-H) metastatic colorectal cancer (mCRC); however, specific predictive biomarkers are lacking. PATIENTS AND METHODS: Data and samples from 85 patients with MSI-H mCRC treated with ICIs were gathered. Tumor infiltrating lymphocytes (TILs) and tumor mutational burden (TMB) were analyzed in an exploratory cohort of "super" responders and "clearly" refractory patients; TILs were then evaluated in the whole cohort of patients. Primary objectives were the correlation between the number of TILs and TMB and their role as biomarkers of ICI efficacy. Main endpoints included response rate (RR), progression-free survival (PFS), and overall survival (OS). RESULTS: In the exploratory cohort, an increasing number of TILs correlated to higher TMB (Pearson's test, p = .0429). In the whole cohort, median number of TILs was 3.6 in responders compared with 1.8 in nonresponders (Mann-Whitney test, p = .0448). RR was 70.6% in patients with high number of TILs (TILs-H) compared with 42.9% in patients with low number of TILs (odds ratio = 3.20, p = .0291). Survival outcomes differed significantly in favor of TILs-H (PFS: hazard ratio [HR] = 0.42, p = .0278; OS: HR = 0.41, p = .0463). CONCLUSION: A significant correlation between higher TMB and increased number of TILs was shown. A significantly higher activity and better PFS and OS with ICI in MSI-H mCRC were reported in cases with high number of TILs, thus supporting further studies of TIL count as predictive biomarker of ICI efficacy. IMPLICATIONS FOR PRACTICE: Microsatellite instability is the result of mismatch repair protein deficiency, caused by germline mutations or somatic modifications in mismatch repair genes. In metastatic colorectal cancer (mCRC), immunotherapy (with immune checkpoint inhibitors [ICIs]) demonstrated remarkable clinical benefit in microsatellite instability-high (MSI-H) patients. ICI primary resistance has been observed in approximately 25% of patients with MSI-H mCRC, underlining the need for predictive biomarkers. In this study, tumor mutational burden (TMB) and tumor infiltrating lymphocyte (TIL) analyses were performed in an exploratory cohort of patients with MSI-H mCRC treated with ICIs, demonstrating a significant correlation between higher TMB and increased number of TILs. Results also demonstrated a significant correlation between high number of TILs and clinical responses and survival benefit in a large data set of patients with MSI-H mCRC treated with ICI. TMB and TILs could represent predictive biomarkers of ICI efficacy in MSI-H mCRC and should be incorporated in future trials testing checkpoint inhibitors in colorectal cancer.

Spatially constrained tumour growth affects the patterns of clonal selection and neutral drift in cancer genomic data.

Quantification of the effect of spatial tumour sampling on the patterns of mutations detected in next-generation sequencing data is largely lacking. Here we use a spatial stochastic cellular automaton model of tumour growth that accounts for somatic mutations, selection, drift and spatial constraints, to simulate multi-region sequencing data derived from spatial sampling of a neoplasm. We show that the spatial structure of a solid cancer has a major impact on the detection of clonal selection and genetic drift from both bulk and single-cell sequencing data. Our results indicate that spatial constrains can introduce significant sampling biases when performing multi-region bulk sampling and that such bias becomes a major confounding factor for the measurement of the evolutionary dynamics of human tumours. We also propose a statistical inference framework that incorporates spatial effects within a growing tumour and so represents a further step forwards in the inference of evolutionary dynamics from genomic data. Our analysis shows that measuring cancer evolution using next-generation sequencing while accounting for the numerous confounding factors remains challenging. However, mechanistic model-based approaches have the potential to capture the sources of noise and better interpret the data.

Catch my drift? Making sense of genomic intra-tumour heterogeneity.

The cancer genome is shaped by three components of the evolutionary process: mutation, selection and drift. While many studies have focused on the first two components, the role of drift in cancer evolution has received little attention. Drift occurs when all individuals in the population have the same likelihood of producing surviving offspring, and so by definition a drifting population is one that is evolving neutrally. Here we focus on how neutral evolution is manifested in the cancer genome. We discuss how neutral passenger mutations provide a magnifying glass that reveals the evolutionary dynamics underpinning cancer development, and outline how statistical inference can be used to quantify these dynamics from sequencing data. We argue that only after we understand the impact of neutral drift on the genome can we begin to make full sense of clonal selection. This article is part of a Special Issue entitled: Evolutionary principles - heterogeneity in cancer? Edited by Dr. Robert A. Gatenby.

Variation of mutational burden in healthy human tissues suggests non-random strand segregation and allows measuring somatic mutation rates.

The immortal strand hypothesis poses that stem cells could produce differentiated progeny while conserving the original template strand, thus avoiding accumulating somatic mutations. However, quantitating the extent of non-random DNA strand segregation in human stem cells remains difficult in vivo. Here we show that the change of the mean and variance of the mutational burden with age in healthy human tissues allows estimating strand segregation probabilities and somatic mutation rates. We analysed deep sequencing data from healthy human colon, small intestine, liver, skin and brain. We found highly effective non-random DNA strand segregation in all adult tissues (mean strand segregation probability: 0.98, standard error bounds (0.97,0.99)). In contrast, non-random strand segregation efficiency is reduced to 0.87 (0.78,0.88) in neural tissue during early development, suggesting stem cell pool expansions due to symmetric self-renewal. Healthy somatic mutation rates differed across tissue types, ranging from 3.5 × 10-9/bp/division in small intestine to 1.6 × 10-7/bp/division in skin.

The shaping and functional consequences of the microRNA landscape in breast cancer.

MicroRNAs (miRNAs) show differential expression across breast cancer subtypes, and have both oncogenic and tumour-suppressive roles. Here we report the miRNA expression profiles of 1,302 breast tumours with matching detailed clinical annotation, long-term follow-up and genomic and messenger RNA expression data. This provides a comprehensive overview of the quantity, distribution and variation of the miRNA population and provides information on the extent to which genomic, transcriptional and post-transcriptional events contribute to miRNA expression architecture, suggesting an important role for post-transcriptional regulation. The key clinical parameters and cellular pathways related to the miRNA landscape are characterized, revealing context-dependent interactions, for example with regards to cell adhesion and Wnt signalling. Notably, only prognostic miRNA signatures derived from breast tumours devoid of somatic copy-number aberrations (CNA-devoid) are consistently prognostic across several other subtypes and can be validated in external cohorts. We then use a data-driven approach to seek the effects of miRNAs associated with differential co-expression of mRNAs, and find that miRNAs act as modulators of mRNA-mRNA interactions rather than as on-off molecular switches. We demonstrate such an important modulatory role for miRNAs in the biology of CNA-devoid breast cancers, a common subtype in which the immune response is prominent. These findings represent a new framework for studying the biology of miRNAs in human breast cancer.

Measuring cancer evolution from the genome.

The temporal dynamics of cancer evolution remain elusive, because it is impractical to longitudinally observe cancers unperturbed by treatment. Consequently, our knowledge of how cancers grow largely derives from inferences made from a single point in time - the endpoint in the cancer's evolution, when it is removed from the body and studied in the laboratory. Fortuitously however, the cancer genome, by virtue of ongoing mutations that uniquely mark clonal lineages within the tumour, provides a rich, yet surreptitious, record of cancer development. In this review, we describe how a cancer's genome can be analysed to reveal the temporal history of mutation and selection, and discuss why both selective and neutral evolution feature prominently in carcinogenesis. We argue that selection in cancer can only be properly studied once we have some understanding of what the absence of selection looks like. We review the data describing punctuated evolution in cancer, and reason that punctuated phenotype evolution is consistent with both gradual and punctuated genome evolution. We conclude that, to map and predict evolutionary trajectories during carcinogenesis, it is critical to better understand the relationship between genotype change and phenotype change. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Epigenetic profile of human adventitial progenitor cells correlates with therapeutic outcomes in a mouse model of limb ischemia.

OBJECTIVE: We investigated the association between the functional, epigenetic, and expressional profile of human adventitial progenitor cells (APCs) and therapeutic activity in a model of limb ischemia. APPROACH AND RESULTS: Antigenic and functional features were analyzed throughout passaging in 15 saphenous vein (SV)-derived APC lines, of which 10 from SV leftovers of coronary artery bypass graft surgery and 5 from varicose SV removal. Moreover, 5 SV-APC lines were transplanted (8×10(5) cells, IM) in mice with limb ischemia. Blood flow and capillary and arteriole density were correlated with functional characteristics and DNA methylation/expressional markers of transplanted cells. We report successful expansion of tested lines, which reached the therapeutic target of 30 to 50 million cells in ≈10 weeks. Typical antigenic profile, viability, and migratory and proangiogenic activities were conserved through passaging, with low levels of replicative senescence. In vivo, SV-APC transplantation improved blood flow recovery and revascularization of ischemic limbs. Whole genome screening showed an association between DNA methylation at the promoter or gene body level and microvascular density and to a lesser extent with blood flow recovery. Expressional studies highlighted the implication of an angiogenic network centered on the vascular endothelial growth factor receptor as a predictor of microvascular outcomes. FLT-1 gene silencing in SV-APCs remarkably reduced their ability to form tubes in vitro and support tube formation by human umbilical vein endothelial cells, thus confirming the importance of this signaling in SV-APC angiogenic function. CONCLUSIONS: DNA methylation landscape illustrates different therapeutic activities of human APCs. Epigenetic screening may help identify determinants of therapeutic vasculogenesis in ischemic disease.

Many private mutations originate from the first few divisions of a human colorectal adenoma.

Intratumoural mutational heterogeneity (ITH) or the presence of different private mutations in different parts of the same tumour is commonly observed in human tumours. The mechanisms generating such ITH are uncertain. Here we find that ITH can be remarkably well structured by measuring point mutations, chromosome copy numbers, and DNA passenger methylation from opposite sides and individual glands of a 6 cm human colorectal adenoma. ITH was present between tumour sides and individual glands, but the private mutations were side-specific and subdivided the adenoma into two major subclones. Furthermore, ITH disappeared within individual glands because the glands were clonal populations composed of cells with identical mutant genotypes. Despite mutation clonality, the glands were relatively old, diverse populations when their individual cells were compared for passenger methylation and by FISH. These observations can be organized into an expanding star-like ancestral tree with co-clonal expansion, where many private mutations and multiple related clones arise during the first few divisions. As a consequence, most detectable mutational ITH in the final tumour originates from the first few divisions. Much of the early history of a tumour, especially the first few divisions, may be embedded within the detectable ITH of tumour genomes.

Intratumor heterogeneity in human glioblastoma reflects cancer evolutionary dynamics.

Glioblastoma (GB) is the most common and aggressive primary brain malignancy, with poor prognosis and a lack of effective therapeutic options. Accumulating evidence suggests that intratumor heterogeneity likely is the key to understanding treatment failure. However, the extent of intratumor heterogeneity as a result of tumor evolution is still poorly understood. To address this, we developed a unique surgical multisampling scheme to collect spatially distinct tumor fragments from 11 GB patients. We present an integrated genomic analysis that uncovers extensive intratumor heterogeneity, with most patients displaying different GB subtypes within the same tumor. Moreover, we reconstructed the phylogeny of the fragments for each patient, identifying copy number alterations in EGFR and CDKN2A/B/p14ARF as early events, and aberrations in PDGFRA and PTEN as later events during cancer progression. We also characterized the clonal organization of each tumor fragment at the single-molecule level, detecting multiple coexisting cell lineages. Our results reveal the genome-wide architecture of intratumor variability in GB across multiple spatial scales and patient-specific patterns of cancer evolution, with consequences for treatment design.

Immune evasion impacts the landscape of driver genes during cancer evolution.

BACKGROUND: Carcinogenesis is driven by interactions between genetic mutations and the local tumor microenvironment. Recent research has identified hundreds of cancer driver genes; however, these studies often include a mixture of different molecular subtypes and ecological niches and ignore the impact of the immune system. RESULTS: In this study, we compare the landscape of driver genes in tumors that escaped the immune system (escape +) versus those that did not (escape -). We analyze 9896 primary tumors from The Cancer Genome Atlas using the ratio of non-synonymous to synonymous mutations (dN/dS) and find 85 driver genes, including 27 and 16 novel genes, in escape - and escape + tumors, respectively. The dN/dS of driver genes in immune escaped tumors is significantly lower and closer to neutrality than in non-escaped tumors, suggesting selection buffering in driver genes fueled by immune escape. Additionally, we find that immune evasion leads to more mutated sites, a diverse array of mutational signatures and is linked to tumor prognosis. CONCLUSIONS: Our findings highlight the need for improved patient stratification to identify new therapeutic targets for cancer treatment.

Cancer Evolution: A Multifaceted Affair.

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery. SIGNIFICANCE: Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

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.