Abstract
Genomic intratumoral heterogeneity (ITH) is common in cancers, but the extent of phenotypic ITH is uncertain because most subclonal mutations are passengers. Since tumor phenotypes are largely driven by epigenetics, methylomic analyses can provide insights into phenotypic ITH. Following this principle, we determined the extent of epigenetic ITH in 16 human colorectal tumors by comparing the methylomes from spatially separated regions in each tumor. Methylomes from opposite tumor sides were similar (Pearson correlation >0.95) with little evidence of ITH or stepwise selection during growth, suggesting that the epigenome of a sampled tumor largely reflects that of its founder cell. Epigenetic conservation was functional, with higher conservation at promoters and expressed genes compared to non-coding regions. Despite epigenomic conservation, RNA expression varied between individual tumor glands, indicating continued adaption during growth. Because many promoters and enhancers were unmethylated, continued adaptation may be due to phenotypic plasticity. Gene enrichment analyses identified that interferon signaling and antigen-processing and presenting pathways were strongly conserved during tumor growth, suggesting a mechanism for immune evasion. In summary, our findings suggest that epigenomes are preferentially conserved during tumor growth and that early tumor cells are poised for rapid growth, phenotypic adaptation, and immune evasion.
Similar content being viewed by others
Introduction
Multiregional sampling reveals that genomic intratumor heterogeneity (ITH) is common in many human tumors1,2,3, yet its clinical significance remains uncertain4. A possible interpretation of genomic ITH is that of ongoing clonal evolution, where successively more fit subclones with different mutations replace less fit subclones. However, most subclonal mutations are passengers5 and it is unclear how or if subclonal mutations confer selection during growth. For instance, quantitative analyses of human colorectal tumors revealed topographic patterns of genetic alterations that are consistent with single expansions of neutral evolution rather than stepwise selection during growth6,7,8. More generally, mutation allelic frequencies of several solid tumor types have been shown to be consistent with single neutral expansions during growth9.
To date, little is known to what extent genomic ITH is translated into phenotypic heterogeneity. Epigenomic variations can lead to phenotypic ITH even in absence of genetic driver mutations, and phenotypic ITH may in turn drive clinical progression of the tumor. In practice, multiregional epigenomic comparisons can be used to infer the degree of phenotypic heterogeneity because DNA methylation of enhancers, promoters, and gene bodies modulates gene expression10,11,12. Epigenomes should differ between tumor sides if distinct subclones with fitter phenotypes emerge stepwise during growth (Fig. 1A). In contrast, with single neutral expansions, epigenomes should be similar between sides because the tumor is essentially a single population of cells with similar phenotypes (Fig. 1B).
To evaluate the degree of phenotypic ITH in human colorectal tumors, we sampled and compared the epigenomes from opposite sides of 16 specimens (Fig. 1). Epigenetic ITH was found to be low, with little evidence of recent stepwise selection. In addition, epigenetic conservation was nonuniform along the genome, with preferential conservation in genes that are functionally important.
Results
Epigenomes are similar between opposite tumor sides
DNA methylation was measured at ~850,000 CpG sites from bulk samples of six normal colons and opposite sides of 16 human colorectal tumors (Table 1; 4 adenomas and 12 colorectal cancers (CRCs)). As seen in other studies13,14,15, there were substantial differences between the epigenomes of tumor and normal tissue within individual patients, and between tumors of different patients (Fig. 2A,B). In contrast, the epigenomes between normal colons of different patients and between opposite sides of the same tumors were strongly conserved (Pearson > 0.95). These findings suggest that epigenomes evolve during progression from normal colon to tumor, but are more conserved between two sides of the same tumor.
Colorectal tumors consist of older cell populations
Other studies have documented that epigenetic ITH is common in CRCs14,15,16, but the biological significance is uncertain, especially because even small epigenetic changes may cause phenotypic differences. Here we sought to determine whether epigenomes and their phenotypes were sufficiently different to confer selection. Epigenomes can detect recent selection because, like DNA sequences, epigenetic patterns are copied almost exactly between generations (Fig. 1A). Cells of a recent expansion initially have similar epigenomes that they inherit from their founder cell, but they become polymorphic with time (“molecular clock hypothesis”17). DNA methylation replication error rates are up to 10,000 times higher than DNA mutation rates17 and therefore can record very recent selective sweeps.
To benchmark the tempo of epigenetic changes, epigenomes were compared between different parts of experimental tumors – both xenografts and in vitro cell lines – started from single cells. Although these experimental “clonal” tumors do not fully model sporadic tumorigenesis and lack immune systems, they can calibrate how fast epigenomes diversify after the start of a clonal expansion. Clonal xenografts grown on different mouse flanks had minimal pairwise differences (PWD) after 4 months (Fig. 2C). Consistent with a molecular clock, clonal serial parallel cell lines originating from a single cell had low PWDs at 2.5 months, and substantially higher PWDs after 14 and 20 months. Compared to these experimental tumors, PWDs between tumor sides were greater than between single cell xenografts but overlapped with the PWDs between older single cell serial parallel cultures. Therefore, the colorectal tumors appear to be older than 4 months of age by the time of surgery and there is no evidence of very recent selective sweeps across the entire tumor specimens. Indeed, such sweeps would have led to lower PWDs between opposite tumor sides compared to the PWDs in experimental systems.
Multiregional gland sampling reveals little evidence of recent selection
To search for smaller localized subclonal sweeps within tumor regions, we sequenced the methylomes of gland pairs from opposite sides of 11 tumors (four glands per tumor). Neighboring glands within bulk regions (<0.5 cm apart) all had PWDs (Fig. 2D) greater than the four-month old xenografts (Fig. 2C) and were almost as related as glands centimeters apart on opposite tumor sides (Fig. 2D). This lack of correlation between physical and epigenomic distances is more consistent with single expansions (Fig. 1B) rather than active ongoing selection, where neighboring glands should be much more related than distant glands (Fig. 1A). PWDs between glands were generally greater than between bulk tumor sides (Fig. 2E), likely reflecting that population averages drift slower than smaller tumor gland subpopulations. Four of the 22 gland pairs within a tumor side had PWDs smaller than between their bulk tumor sides, which could reflect regional subclonal selective sweeps during growth. Any such selection did not appear recent because PWDs between glands were all greater than the four-month old xenografts (Fig. 2C).
Preferential epigenetic conservation at functional regions
The relatively uniform spatial epigenetic ITH within and between tumor regions indicates drift rather than stepwise selection is the major driver of epigenetic ITH. For a single neutral expansion, the epigenome of the tumor should resemble the epigenome of its founder cell. Furthermore, subsequent drift during growth is expected to be nonuniform along the genome, with preferential conservation at functional regions because DNA methylation modulates expression10,12 (Fig. 1B). More precisely, methylation changes in non-genic regions are expected to be neutral (passengers) whereas methylation changes in expressed genes are expected to be subject to positive or negative selection. Consistent with this drift hypothesis, we found increased conservation at gene-associated CpG sites, with preferential conservation at more highly expressed genes and a subset of genes associated with immune surveillance (Fig. 2E). Within genes, there was greater conservation at promoters (TSS200) and in first exons. Conserved methylation therefore identifies regions in a founder cell epigenome that are selectively maintained during growth, in contrast to non-functional regions that are subject to continued drift (Fig. 2F).
A core set of conserved genes
Preferential epigenetic conservation enabled us to identify genes that were under greater selection during growth. Among the 22,430 annotated genes with four or more CpG sites, we identified four subgroups based on pairwise distances between tumor sides and between normal colons of different patients: super-conserved genes (conserved between tumor sides and between normal colons), tumor-only conserved, normal colon-only conserved and never-conserved (Fig. 3A, Table 2). Super-conserved genes (n = 2,769) were defined as having very similar epigenomes between tumor sides (PWD < 0.05) in 90% or more of tumors, and between normal colons of different patients (PWD < 0.075), which corresponds to little change during progression and growth. Among super-conserved genes, a pathway enrichment analysis with Reactome18 identified significant enrichment (FDR < 0.05) for 612 loci involved in 66 pathways important for general cell metabolism. Many genes involved with core cell functions appear to be under selection in normal colon and tumors because their methylation remains conserved during progression (Table 2).
No loci were consistently conserved only between normal colons but not between tumor sides, indicating that most genes that are highly regulated in the colon are also under selection during growth. A total of 2,977 loci were conserved neither between tumor sides nor between normal colons (PWD > 0.05 between tumor sides in >90% of tumors, PWD > 0.075 between normal colons), with significant nonsensical enrichment in a collagen chain trimerization pathway.
Preferential tumor-specific conservation of immune surveillance genes
Areas of the epigenome that are only conserved between tumor sides and not between normal colons can potentially identify which genes and pathways are uniquely important for tumorigenesis. There were 242 loci consistently conserved only in tumors (PWD < 0.05 between tumor sides in >90% of tumors) but not between normal colons (PWD > 0.075), and pathway analyses with Reactome revealed significant enrichment of 9 genes and 9 immune surveillance (IS) pathways (Table 2).
To facilitate conservation analyses of individual tumors, we introduced an extended tumor-only classification that filtered out genes that are inherently epigenetically stable. More precisely, across all tumors, we identified 1,767 genes that were conserved between tumor sides (average PWD across all tumors <0.05) but not conserved between normal colon and tumor of the same patient (average PWD > 0.075), normal colons of different patients (average PWD > 0.075), and tumors of different patients (average PWD > 0.075). In this extended tumor-only group, there was again significant pathway enrichment of 19 tumor-specific IS genes involved in antigen-presentation and processing, and in interferon signaling (Table 2). Next, among the 1,767 genes in the extended tumor-only group, we focused on individual tumors and identified which genes were conserved between the sides of individual tumors (tumor-specific PWD < 0.05). Through this process, IS gene enrichment was found in 15 of the 16 tumors (Table S1). The only tumor without preferential IS gene conservation enrichment had a POLE mutation, suggesting an alternative mechanism for immune evasion when the mutation burden is exceptionally high.
The preferential conservation of IS genes between different regions of the same tumor is illustrated on the examples of HLA-A, HLA-B, and HLA-C in Fig. 3C,D. IS epigenetic conservation within tumors is greater than between tumors or between multiple samples of individual colons (Fig. 3E). Tumor-specific preferential conservation of antigen processing and presentation genes suggests immune surveillance exerts strong selective pressures on most human colorectal tumors during growth.
Tumor-specific conservation is likely due to selection, but other factors such as greater replication fidelity could also explain a lack of drift. One way to test for selection is to determine whether IS gene conservation persists in cell lines which are not under immune surveillance. Unlike between tumor sides, IS genes were not preferential conserved in the clonal parallel cell lines relative to all genes (Fig. 3E), indicating that methylation at these genes drift more in the absence of an immune system.
Immune surveillance gene methylation differs within healthy adult colons
Patient-specific somatic mutations in human tumors lead to the expression of patient-specific neoantigens, which in turn are presented to the immune system. This can lead to immunoediting, whereby the immune system removes certain neoantigen-presenting cancer cells19. The conservation of IS genes as observed in our study samples may prevent immunoediting and enable the tumor to grow despite active immune surveillance. More precisely, tumor-specific IS gene conservation may indicate that highly specific processing and presentation of neoantigens may enable tumors to avoid immunoediting. To test this hypothesis, we performed multiregional measurements of methylation in healthy colons, where crypts independently accumulate many different mutations with aging20. If IS gene methylation customizes the presentation of specific neoantigens, then conservation should not be observed between different colonic regions because they have acquired different somatic mutations and hence neoantigens over time. Consistent with this hypothesis, while the methylation of highly expressed genes was conserved, IS gene methylation drifted between the four different parts each of nine healthy adult colons (Fig. 3E). In summary, the difference in IS gene conservation in tumors compared to healthy colons (Fig. 3E) suggests that each tumor may be locked into its own IS gene methylation pattern that enables successful evasion of immunoediting.
Phenotypic plasticity despite epigenetic conservation
In most human tumors, many different microscopic phenotypes can be observed. In the selective sweep scenario (Fig. 1A), different phenotypes can arise by stepwise selection. In the single expansion scenario (Fig. 1B), phenotypic heterogeneity can arise thanks to phenotypic plasticity21,22 which allows cells with similar epigenomes to express multiple phenotypes. Such phenotypic plasticity is observed in the normal colon where microscopically different cell types in a single crypt have nearly identical epigenomes23,24. In particular, broadly permissive24 or open chromatin at enhancers and promoters license plasticity because transcription factors induced by the microenvironment determine final cell phenotypes.
In our study samples, the vast majority of promoters were unmethylated in both normal colons and tumors (Fig. 4A), indicating open configurations10 and suggesting phenotypic plasticity. In contrast, enhancer methylation was bimodal, with lower mean methylation in tumors (46%) compared to normal colons (57%). Hence, the open chromatin at promoters in normal colon is largely maintained in tumors, and enhancers are less methylated or more open in tumors, potentially allowing for plasticity and rapid adaptation when cells encounter new microenvironments.
Finally, we sought to determine the relationship between DNA methylation and gene expression by performing RNA sequencing at the tumor gland scale. More precisely, mRNA levels were measured in four glands (two from each side) in 3 tumors (H, R, J). The mRNA expression correlation between glands on the same tumor side was larger compared to glands on opposite sides of the same tumor or glands from different tumors (Fig. 4B). In individual tumors, expression levels differed between glands, with 0.3-, 0.4- and 0.5- fold mean expression differences between glands in tumors J, H and R, respectively. However, there was no notable correlation between the degree of epigenetic conservation and RNA expression variability (Fig. 4C). In summary, although DNA methylation is conserved, gene expression is variable between glands, which may reflect phenotypic plasticity and adaptation of single growing tumor populations to heterogeneous local microenvironments.
Discussion
How human tumors grow is difficult to ascertain because longitudinal observations are generally not feasible. Genome comparisons offer a practical way to reconstruct in vivo tumorigenesis because genomes record ancestries by becoming polymorphic. Multiregional sequencing studies reveal that mutational ITH is common in human tumors1,2,3, but it is uncertain if this genetic ITH is translated into stepwise selective sweeps during growth because most subclonal mutations are neutral passengers. In the case of human CRCs, the patterns of genetic ITH have been found to be more consistent with single neutral expansions, indicating that sequential stepwise selection may be uncommon during growth6,8,11. In this study, we performed multiregional epigenetic sampling and found that spatial epigenetic ITH is very limited, further corroborating that stepwise selection is infrequent or spatially limited to small sections during growth of human CRCs. Notably, DNA passenger methylation is also heterogeneous within individual tumor glands25 indicating that selective sweeps within glands are rare as well.
Although inferential, genome conservation is often used to identify functional regions under selection and is useful when direct observations or experimental manipulations are impractical. Indeed, epigenetic error rates during cell division are high compared to DNA error rates, and the resulting drift can create epigenetic ITH along the genome, with greater changes at non-functional CpG sites. As expected, we found preferential epigenetic conservation of highly expressed genes involved in core cellular functions, indicating that these genes are important in both normal colon and tumors. Interestingly, a search for loci that were specifically conserved in tumors but not in normal tissue revealed several genes that are involved in interferon signaling and antigen-processing and presentation.
Because HLA expression and tumor antigenicity are modulated by DNA methylation26,27,28, the observed conservation of IS genes stipulates a potential mechanism for evasion of immune surveillance during the growth of CRCs29. Whether immune surveillance is a significant extrinsic barrier to human tumorigenesis is unclear19. In the case of CRCs, T-cells reactive to neoantigens are present30 and mutations that could restrict neoantigen presentation to T-cells are rare31. Nevertheless, anti-PD-L1 therapies are ineffective against non-mutator phenotypes CRCs32, suggesting that other mechanisms such as chromatin regulation33,34 may contribute to evasion of immune surveillance. Finally, we note that our findings are consistent with recent CRISPR-Cas9 studies that found multiple IS genes influence immunoediting35,36.
The preferential conservation of IS gene methylation patterns between tumor sides – but not between tumor and normal tissues – implies that the founding tumor cells are under selection right from the start of tumor growth. Immune evasion may not be limited to tumors because mutations accumulate with age in normal colon20, and immune cytolytic activity is higher in the colon than in CRCs, with neoantigens depleted but not absent in CRCs31. Consistent with independent mutations and a customized mechanism where distinct neoantigens must be processed and presented differently to avoid immunoediting, IS gene methylation varied between different parts of normal adult colons. Because every founder cell lineage has a different set of neoantigens and must have previously evaded immunoediting, their progeny may simply perpetuate a successful incumbent strategy, resulting in the observed tumor-specific IS gene conservation.
The tumor microenvironment may play an important role in shaping the phenotypic ITH characterized by cells of different microscopic appearances and variable mRNA expression between regions or cells of the same CRCs2,37,38. If there was a strict correspondence between epigenome and phenotype, the observed epigenetic conservation would strongly limit phenotypic repertoire. However, phenotypic plasticity21,22 allows a single epigenome to confer multiple phenotypes, and broadly permissive chromatin in the colon allows crypt cells with different phenotypes to share similar epigenomes23,24. We found that the unmethylated promoters in normal colon were largely maintained in the tumors, and similar to other studies39, enhancers were more open in the tumors than adjacent normal tissue. Unmethylated promoters and more open enhancer configurations indicate that transcription factors induced by the microenvironment could potentially dictate tumor cell phenotypes. Hence, instead of waiting for stepwise evolution, a conserved but functionally plastic epigenome allows cells to immediately adapt to the many heterogeneous microenvironments encountered during growth. The ability to rapidly adapt rather than evolve through genetic bottlenecks may facilitate growth because tumor microenvironments may vary considerably temporally and spatially, even over a few millimeters40.
In conclusion, multiregional sampling reveals epigenetic ITH is minimal and more consistent with neutral evolution than frequent stepwise selection during primary human colorectal tumor growth. Importantly, this approach can characterize the degree of selection during growth, and epigenetic conservation can indicate which genes are more likely under selection. Our conclusions depend on highly purified tumor specimens because methylation differences can reflect variations in normal contaminating cells41. This epigenetic conservation provides new insights into the epigenome of the founder tumor cell whose stem like phenotypic plasticity can help explain the cellular heterogeneity and drug resistance of human tumors21. Unlike stepwise progression, where progenitors lack a complete set of drivers needed for robust growth, the founder cells of a neutral single expansion are already highly capable and adaptable, equipped with driver mutations and a largely open epigenome needed for rapid growth. Although further progression during growth may be rare, neutral evolution can be highly detrimental to the patient as it maximizes diversity42 and therefore increases the chances for further stepwise evolution (such as metastasis) and for selection of resistant subclones under therapeutic pressure. Taken together, both genetic7,8 and epigenetic multiregional measurements of human CRCs reconstruct portraits of founder cells that are born already poised for growth, phenotypic adaptation, and with an incumbent immune evasion strategy, helping to explain rapid expansions into visible “Big Bang” tumors.
Methods
Tumors
The tumor (n = 16) and matched normal colon (n = 6) samples (Table 1) were obtained fresh at the Norris Comprehensive Cancer Center, and samples (~0.5 cm3) from opposite tumor sides and from different parts of the normal colons were procured, as previously described7. The specimens were procured in accordance with relevant guidelines and regulations, and all experimental protocols were approved by the Institutional Review Board at the University of Southern California. Informed consent was not obtained because the specimens were considered exempt by our Institutional Review Board. Other data from 10 of the tumors have been previously published7 and have the same names (except the tumor originally labeled R). From each bulk sample (tumor or normal colon), individual tumor glands (10,000 to 30,000 cells) or intestinal crypts were isolated by an EDTA washout that results in >90% epithelial cell purity43. Examples of isolated tumor glands from our study and corresponding LUMP purity estimates are provided in Fig. S1. Bulk methylation levels were measured by pooling hundreds of individual tumor glands or normal colon crypts.
DNA methylation
Bulk sample DNA methylation was measured using the EPIC methylation array (Illumina) for opposite sides of 16 tumors, gland pairs from opposite sides of 11 CRCs, and six normal colons. Idat files in each color channel (red and green) were processed with the ‘noob’ function in the minfi R program to extract the methylated and unmethylated probe signal intensities and to calculate the beta value for each data point after within-array normalization for background correction, dye-bias equalization and probe design. Beta values between samples were compared based on pairwise distances (PWDs) defined as the absolute differences in beta values. Only autosomal CpG sites were analyzed.
Methylation of four widely spaced samples from nine normal adult colons (ages 50 to 70 years, mucosal biopsies, not EDTA washout-purified) was also measured with the Infinium methylation arrays (Illumina). Raw intensity data was read, preprocessed and batch corrected using the minfi (v.1.16.0) in R (v3.2.2). Probes were mapped to the genome and probes showing mean intensity p-value > 0.05 across all samples were excluded. Excluded probes also included: probes with proximal SNPs; non-CpG probes; probes on sex chromosomes; probes with observed cross-reactivity between two or more genomic regions. After subsequent filtering, 426,718 probes were left for analysis.
Gene annotation
The CpG sites were classified based on the UCSC annotations provided with the array. Enhancers were defined by FANTOM544. The panel of genes with higher expression in CRCs (RPKM > 1) was obtained from the non-hypermutated data in Supplemental Table 3 in ref.9.
Pathway enrichment analysis
Autosomal genes with at least four annotated CpG sites (N = 22,430) were classified based on methylation differences between normal colon and tumor (average colon methylation versus average tumor methylation) and between sides (PWD between tumor sides). There were 5,988 genes consistently (90% or more of tumors) with the four types of conservation in Fig. 3A and Table 2. The gene lists were analyzed with the Reactome pathway enrichment algorithm18. A false detection rate (FDR) cutoff of <0.05 was used.
Another method to find genes preferentially conserved only during tumor growth identified genes with average PWDs greater than 0.075 between normal colons from different patients, between normal and tumor from the same patient, and between tumors from different patients, but with PWDs less than 0.05 between opposite tumor sides from the same patient. The gene lists from average and individual tumor data were analyzed with the Reactome algorithm.
RNA measurements
RNA was measured using the NanoString GX Human Cancer Reference Kit which measures 230 cancer-related RNAs. Single tumor glands (~10,000 cells each, 2 glands from each side of tumors H, J, R from Table 1) were analyzed directly. Data were normalized for each tumor using the geometric mean with nSolver 3.0 software. To minimize measurement variability, only genes with at least 30 counts in all four glands and with CpG sites in the EPIC array (136 genes for J, 101 genes for R, 84 gene for H) were compared.
Cell line studies
A single cell was isolated from a bulk tissue culture of HCT116 (American Type Culture Collection), grown to confluence, and split into three flasks. Each of the three parallel lines were serially cultured for 18 to 24 passages (52 to 94 days), 70 passages (409 days) and 105 passages (682 days). DNA methylation was measured using EPIC arrays and compared between the samples.
Xenograft studies
Single cells were isolated from bulk tissue cultures and grown to several million cells. About 1 million cells were injected subcutaneously into immune deficient mice (nu/nu) and xenografts (~1 cm3) were harvested within 3 months after injection. DNA was isolated from physically distinct clonal xenografts (separate flanks or different mice) derived from HCT116 (2 xenografts), SW480 (two xenografts) and LoVo (three xenografts). These mouse experiments were performed in accordance with relevant guidelines and regulations, and were approved by the Institutional Animal Care and Use Committee at the University of Southern California. DNA methylation was measured using EPIC arrays.
References
Gerlinger, M. et al. Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. N. Engl. J. Med. 366, 883–892 (2012).
McGranahan, N. & Swanton, C. Clonal heterogeneity and tumor evolution: past, present, and the future. Cell 168, 613–628 (2017).
Shibata, D. Cancer. Heterogeneity and tumor history. Science 336, 304–305 (2012).
Andor, N. et al. Pan-cancer analysis of the extent and consequences of intratumor heterogeneity. Nat. Med. 22, 105–113 (2016).
Vogelstein, B. et al. Cancer genome landscapes. Science 339, 1546–1558 (2013).
Ryser, M. D., Min, B. H., Siegmund, K. D. & Shibata, D. Spatial mutation patterns as markers of early colorectal tumor cell mobility. Proc. Natl. Acad. Sci. USA 115, 5774–5779 (2018).
Sottoriva, A. et al. A Big Bang model of human colorectal tumor growth. Nat. Genet. 47, 209–216 (2015).
Sun, R. et al. Between-region genetic divergence reflects the mode and tempo of tumor evolution. Nat. Genet. 49, 1015–1024 (2017).
Williams, M. J., Werner, B., Barnes, C. P., Graham, T. A. & Sottoriva, A. Identification of neutral tumor evolution across cancer types. Nat. Genet. 48, 238–244 (2016).
Jones, P. A. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat. Rev. Genet. 13, 484–492 (2012).
Shlyueva, D., Stampfel, G. & Stark, A. Transcriptional enhancers: from properties to genome-wide predictions. Nat. Rev. Genet. 15, 272–286 (2014).
Yang, X. et al. Gene body methylation can alter gene expression and is a therapeutic target in cancer. Cancer Cell. 26, 577–590 (2014).
Cancer Genome Atlas Network. Comprehensive molecular characterization of human colon and rectal cancer. Nature 487, 330–337 (2012).
Martínez-Cardús, A. et al. Epigenetic homogeneity within colorectal tumors predicts shorter relapse-free and overall survival times for patients with locoregional cancer. Gastroenterology 151, 961–72 (2016).
Uchi, R. et al. Integrated multiregional analysis proposing a new model of colorectal cancer evolution. PLoS Genet. 12, e1005778, https://doi.org/10.1371/journal.pgen.1005778 (2016).
Roerink, S. F. et al. Intra-tumour diversification in colorectal cancer at the single-cell level. Nature 556, 457–462 (2018).
Shibata, D. & Tavare, S. Counting divisions in a human somatic cell tree: how, what and why? Cell Cycle 5, 610–614 (2006).
Fabregat, A. et al. The Reactome pathway Knowledgebase. Nucleic Acids Res. 44, D481–487 (2016).
Schumacher, T. N. & Schreiber, R. D. Neoantigens in cancer immunotherapy. Science 348, 69–74 (2015).
Blokzijl, F. et al. Tissue-specific mutation accumulation in human adult stem cells during life. Nature 538, 260–264 (2016).
Easwaran, H., Tsai, H. C. & Baylin, S. B. Cancer epigenetics: tumor heterogeneity, plasticity of stem-like states, and drug resistance. Mol. Cell 54, 716–727 (2014).
Meacham, C. E. & Morrison, S. J. Tumour heterogeneity and cancer cell plasticity. Nature 501, 328–337 (2013).
Kaaij, L. T. et al. DNA methylation dynamics during intestinal stem cell differentiation reveals enhancers driving gene expression in the villus. Genome Biol. 14, R50 (2013).
Kim, T. H. et al. Broadly permissive intestinal chromatin underlies lateral inhibition and cell plasticity. Nature 506, 511–515 (2014).
Siegmund, K. D., Marjoram, P., Tavare, S. & Shibata, D. High DNA methylation pattern intratumoral diversity implies weak selection in many human colorectal cancers. PLoS One 6, e21657, https://doi.org/10.1371/journal.pone.0021657 (2011).
Campoli, M. & Ferrone, S. HLA antigen changes in malignant cells: epigenetic mechanisms and biologic significance. Oncogene 27, 5869–5885 (2008).
Fonsatti, E. et al. Functional up-regulation of human leukocyte antigen class I antigens expression by 5-aza-2′-deoxycytidine in cutaneous melanoma: immunotherapeutic implications. Clin. Cancer Res. 13, 3333–3338 (2007).
Siebenkas, C. et al. Inhibiting DNA methylation activates cancer testis antigens and expression of the antigen processing and presentation machinery in colon and ovarian cancer cells. PLoS One 12, e0179501, https://doi.org/10.1371/journal.pone.0179501 (2017).
DuPage, M., Mazumdar, C., Schmidt, L. M., Cheung, A. F. & Jacks, T. Expression of tumour-specific antigens underlies cancer immunoediting. Nature 482, 405–409 (2012).
Tran, E. et al. Immunogenicity of somatic mutations in human gastrointestinal cancers. Science 350, 1387–1390 (2015).
Rooney, M. S., Shukla, S. A., Wu, C. J., Getz, G. & Hacohen, N. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell 160, 48–61 (2015).
Le, D. T. et al. PD-1 Blockade in Tumors with Mismatch-Repair Deficiency. N. Engl. J. Med. 372, 2509–2520 (2015).
Miao, D. et al. Genomic correlates of response to immune checkpoint therapies in clear cell renal cell carcinoma. Science 359, 801–806 (2018).
Pan, D. et al. A major chromatin regulator determines resistance of tumor cells to T cell-mediated killing. Science 359, 770–775 (2018).
Manguso, R. T. et al. In vivo CRISPR screening identifies Ptpn2 as a cancer immunotherapy target. Nature 547, 413–418 (2017).
Patel, S. J. et al. Identification of essential genes for cancer immunotherapy. Nature 548, 537–542 (2017).
Dunne, P. D. et al. Cancer-cell intrinsic gene expression signatures overcome intratumoural heterogeneity bias in colorectal cancer patient classification. Nat. Commun. 8, 15657, https://doi.org/10.1038/ncomms15657 (2017).
Li, H. et al. Reference component analysis of single-cell transcriptomes elucidates cellular heterogeneity in human colorectal tumors. Nat. Genet. 49, 708–718 (2017).
Akhtar-Zaidi, B. et al. Epigenomic enhancer profiling defines a signature of colon cancer. Science 336, 736–739 (2012).
Heindl, A., Nawaz, S. & Yuan, Y. Mapping spatial heterogeneity in the tumor microenvironment: a new era for digital pathology. Lab. Invest. 95, 377–384 (2015).
Aran, D., Sirota, M. & Butte, A. J. Systematic pan-cancer analysis of tumour purity. Nat. Commun. 6, 8971 (2015).
Ling, S. et al. Extremely high genetic diversity in a single tumor points to prevalence of non-Darwinian cell evolution. Proc. Natl. Acad. Sci. USA 112, E6496–505 (2015).
Nakamura, S., Goto, J., Kitayama, M. & Kino, I. Application of the crypt-isolation technique to flow-cytometric analysis of DNA content in colorectal neoplasms. Gastroenterology 106, 100–107 (1994).
Andersson, R. et al. An atlas of active enhancers across human cell types and tissues. Nature 507, 455–461 (2014).
Acknowledgements
Supported by grants from the National Institutes of Health (CA185016, CA196569, P30CA014089, K99CA207872, U54CA217376, CA194663, CA220004, U01CA206110), the National Science Foundation (DMS 1614838), the Swiss National Science Foundation (P300P2-154583), and R.A.C.E. Charities, Cottrell Family Fund, and Rodger Haggitt Endowed Chair to WMG.
Author information
Authors and Affiliations
Contributions
D.S. designed the experiments. D.S., K.S. and M.D.R. analyzed the data and wrote the manuscript. M.Y. and W.M.G. provided the multiregional methylation data on the normal human intestines.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ryser, M.D., Yu, M., Grady, W. et al. Epigenetic Heterogeneity in Human Colorectal Tumors Reveals Preferential Conservation And Evidence of Immune Surveillance. Sci Rep 8, 17292 (2018). https://doi.org/10.1038/s41598-018-35621-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-35621-y
Keywords
This article is cited by
-
Minimal barriers to invasion during human colorectal tumor growth
Nature Communications (2020)
-
Translating insights into tumor evolution to clinical practice: promises and challenges
Genome Medicine (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.