Main

Neoplasms evolve1,2,3. This evolution has been recognized since 1976 (Ref. 4), and it explains the processes of both carcinogenesis and acquired therapeutic resistance1. The evolution of neoplasms is shaped by the selective pressures of their microenvironmental ecology. But between and within cancer types, tumours probably display differences in the dynamics of cancer evolution and ecology, including the rates at which new clones appear and go extinct, how different those clones are from one another and whether they appear in bursts or at a more regular pace. Many of the evolutionary and ecological properties of a neoplasm are clinically relevant5,6,7,8,9,10,11,12,13,14,15,16, though this is not always true6,16,17, and in most cases their clinical relevance has not yet been tested. There is a need for a common language and conceptual categories for drawing clinical distinctions that capture the relevant genetic, environmental and kinetic parameters that impact tumour adaptation and progression, as well as response to therapy. A classification system for the evolution and ecology of neoplasms would provide clinicians and researchers with a foundation for developing better prognostic and predictive assessments of tumour behaviour, such as response to an intervention.

The ultimate purpose of a classification system for the evolution and ecology of neoplasms is to provide a descriptive tool by which to improve clinical management with respect to the overall survival and quality of life of the patient. It would also help to drive research and discovery in cancer biology and oncology.

Below, we discuss the methods by which we reached consensus as well as the goals and guiding principles we aspired to in the development of a framework for classifying neoplasms. We then discuss each of the components of the classification system as well as methods for measuring them and for dividing tumours into an initial set of 16 classes. We discuss how such a classification system could be developed, improved and used clinically in the future.

Methods

We convened a consensus conference of experts in the fields of cancer evolution and cancer ecology to lay the groundwork for the development of an evolutionary and ecological classification system. The initial participants (Maley, Aktipis, Graham, Sottoriva, Boddy, Janiszewska, Silva, Gerlinger, Anderson, Brown and Shibata) were among the faculty for the Evolution and Ecology of Cancer summer school funded by Wellcome and held at the Wellcome Genome Campus in Hinxton, UK, in July of 2016. Input from all participants was solicited, and after discussion, we identified areas of consensus. Afterwards, other leaders in the field were invited to join the effort by co-editing and discussing the developing statement. All authors reviewed and approved the final statement. Wellcome Genome Campus Advanced Courses and Scientific Conferences provided financial support for the consensus meeting. We have named the classification system, with their permission, in appreciation of Wellcome's support. Please note that the statement reflects the opinions of the authors and not necessarily those of Wellcome.

Goals and guiding principles

Our development of this framework has been guided by several goals and principles. We agreed that an ideal classification system should have the following properties. First, it must be able to alter a clinical decision point. Second, it should be simple enough to be easily remembered and applied. Third, it should also align with our current understanding of the dynamics of neoplasms. Fourth, the classification system should be general enough to be applied across different types of neoplasm, recognizing that the types of measurement may need to be individualized to a given type of cancer.

This framework is based on fundamental theoretical principles underlying evolutionary and ecological dynamics. It is not based on any particular assay or parameter but rather captures the fundamental drivers of tumour evolution. This is a necessary first step that we hope will lead to many methodological and measurement innovations to quantify the key components of tumour evolution and ecology that we identify here. Because the evolution of cancer is still a relatively new field, there is still uncertainty about the best ways to measure and describe the evolution and ecology of a tumour.

There are also practical considerations in the construction of a classification system. If a tumour could be classified based on a single biopsy from standard assays such as those that can be done on formalin-fixed paraffin-embedded (FFPE) tissue or standard radiological images, translation to the clinic would be relatively easy. However, studies have not yet been done to test whether measures of the evolvability of a tumour from a single biopsy sample are sufficient or whether multiple samples substantially improve predictions of clinical outcomes15. We hypothesize that we will need to extensively sample neoplasms over both space and time in order to accurately quantify their evolvability, but this remains an open question. It is clear, however, that evolutionary analyses are limited if the clonal structure of the primary tumour is unknown18. The use of cell-free DNA (cfDNA) from liquid biopsy samples should facilitate longitudinal studies19, although deconvoluting the clones within such a mixed sample remains a challenge20.

Framework for classifying tumour evolution

There are many well-established ways to classify tumours, largely based on extent of spread and morphological appearances (for example, stage and grade). An evolutionary classification system would augment current schemes by further capturing the evolvability of a tumour. How much intrinsic genetic instability does it have? How likely is it to respond quickly to a new selective pressure such as a therapeutic intervention? For example, rapid progression after chemotherapy is probably driven by pre-existing resistant variants, and therefore, failure is more likely in tumours with more subclonal diversity (intratumoural heterogeneity)6. Moreover, it would be useful to classify evolution through time. For example, a second biopsy from the same patient after therapy may reveal minimal diversity, indicating a recurrent tumour derived from a single clone, or substantial diversity, suggesting intrinsic resistance by the majority of tumour cells. There was widespread agreement at the consensus conference that both the evolutionary dynamics of the neoplastic cells themselves (cancer cell intrinsic factors) and the microenvironment that defines the ecology of those cells (cancer cell extrinsic factors) are important in predicting the future behaviour and response of a tumour. To capture this, we have developed a framework for both an evolutionary index (Evo-index) that describes the intrinsic evolvability of the neoplastic cell population and an ecological index (Eco-index) that describes potential selective pressures imposed by the surrounding microenvironment.

The Evo-index

The Evo-index (D#Δ#) is a combination of two fundamental components: the diversity (D) or intratumoural heterogeneity of the neoplasm and how it changes over time (Δ). In other words, the Evo-index quantifies heterogeneity in both space and time (Fig. 1a). Both diversity and changes in the clonal structure of a tumour over time are objective measures and may be assessed as part of preclinical studies or clinical trials.

Figure 1: The Evo-index and how it changes.
figure 1

a | The evolutionary index (Evo-index) is composed of two factors corresponding to heterogeneity over space (diversity, D) and heterogeneity over time (change over time, Δ). By 'change', we mean both change in the genetic, epigenetic and phenotypic alterations present in the population and change in the frequencies of those alterations in the neoplastic cell population. What measures of D and Δ are best is an open question. In addition, how these factors should be stratified into two, three or more classes is also an open question. Here, for simplicity, we provide examples of the kinds of dynamics that could be categorized into a simple 2 × 2 classification. b | The genetic composition of a tumour may change either slowly (Δ1) or rapidly (Δ2) in a variety of ways. On the left, a tumour may have low diversity (D1) at time 0 because it is a new tumour or there has been a recent homogenizing clonal expansion. That tumour may be quiescent and so appear substantially the same at time 1 (D1Δ1), or it may accumulate clones, some of which expand, to generate a diverse tumour by time 1 (D2Δ2). Alternatively, a tumour may be diverse (D2) at time 0 because it is old or has a high mutation rate and is evolving neutrally. At time 1, that tumour may have been homogenized by a selective sweep (D1Δ2) or may continue on its current trajectory with gradual turnover of its clones (D2Δ1).

PowerPoint slide

Diversity. The heterogeneity that is currently present in a population defines its capacity to respond, at a population level, to selective pressures. This diversity is the fuel for the engine of natural selection. There are different forms of diversity, including genetic diversity, epigenetic diversity, phenotypic diversity and functional diversity. Genetic diversity can predict progression to invasive cancers12,13 as well as recurrence and survival5,6,7,8,9,16. The relationship between diversity and clinical outcomes is not universally consistent across different cancer types6,16 and can be complicated (Box 1).

Diversity can be a proxy for the likelihood that a resistant clone is present in a neoplasm. We currently do not know all the mutations and epigenetic alterations that make a neoplastic cell resistant to a particular therapy, and even those we do know are difficult to detect if they are present in only a small region of the tumour. Compared with homogeneous neoplasms, diverse neoplasms are more likely to harbour resistant clones and are also probably more likely to evolve resistance in the future.

Multiple forms of diversity within a neoplasm may be clinically important, not only as fuel for natural selection but also as biomarkers of clinically targetable dynamics. For example, if high levels of genetic diversity are indicative of high levels of moderately deleterious passenger mutations21,22, then suppressing mechanisms in the cell that buffer against those deleterious effects, such as chaperone proteins, should preferentially harm neoplastic cells21. Alternatively, diversity may be indicative of cooperation between clones, through mechanisms such as cross feeding23,24,25,26,27. These mechanisms of cooperation are themselves potential therapeutic targets. Theory suggests that targeting cancer cell cooperation should provide weaker selection for resistance than cytotoxic therapies28.

It is likely that not all forms of diversity are equal, and future work must test which are clinically relevant. It may be the case that measures of functional diversity or even phenotypic diversity are better predictors of clinical outcomes than measures of genetic diversity (as many genetic mutations will have no phenotypic consequence), and the ideal measures may vary between tumour types.

Measuring diversity. Of the four components of the classification framework, the largest number of methods has been developed for measuring diversity (intratumoural heterogeneity)13,25,29 (Table 1). There is a large literature in ecology on the quantification of diversity30. The overall diversity of a large area, or landscape (gamma diversity), can be broken down into the diversity within local regions (alpha diversity) and the differences between regions (beta diversity)31. Inherent in this definition is the concept that measuring diversity requires defining the spatial scale that one is examining. One might define within-region diversity as the diversity measured within a biopsy sample, while between-region diversity would account for differences between biopsy samples in multi-region sampling studies. Alternatively, one could take a sample across an entire tumour, perhaps using cfDNA, and estimate the diversity of the entire population. Most of the studies to date have focused on within-region diversity5,6,32,33 or the diversity of the entire tumour12,13,14,25. The use of ecological statistics for measuring between-region diversity in tumours remains relatively unexplored. Established measures of differences between microbial communities34 could possibly be applied to measuring differences between biopsy samples.

Table 1 Measures and assays for the factors that go into the Evo- and Eco-indices

There are many ways to measure diversity30 and a number of challenges to measuring diversity in neoplasms, as discussed in Box 1. In Barrett oesophagus, Merlo and colleagues tested many of those measures of diversity and found that high levels of diversity were predictive of progression to cancer, regardless of the measure13,14. Because evolution is driven by the fitness outliers35, and it may take only one resistant cell at diagnosis to eventually cause drug resistance or relapse after therapy, much of the predictive value of measuring diversity may lie in the long tail of rare clones. Because of this, we recommend using either a count of the number of clones ('species richness') or Shannon index, which equally weights number and relative abundance of clones, to quantify diversity30.

The feasibility of obtaining a complete picture of the diversity of a neoplasm, through multi-region sampling or cfDNA, varies across tumour types. In Barrett oesophagus, bladder cancer and prostate cancer, multi-region sampling is part of the current standard of care36,37,38. In a well-mixed neoplasm, such as a blood cancer, a single sample may be sufficient, but it requires single-cell assays, which have their own challenges (Box 1). In other tumours that are difficult to sample, such as pancreatic cancers, we are lucky to get more than one biopsy sample. The main challenge in using cfDNA is detecting it in serum for cancers that have not yet metastasized, although the level of tumour cfDNA in serum varies across cancer types. A recent study was able to detect tumour cfDNA in 97% of early-stage lung squamous cell carcinomas but only 19% of early-stage lung adenocarcinomas39.

The interpretation of the diversity of a neoplasm depends on the context of its history. A neoplasm that has just been homogenized by a therapy that killed most of the clones in that neoplasm is different from a neoplasm that is homogeneous because it has a very low mutation rate and has not had enough time to accumulate many clones. By contrast, a high level of diversity in a neoplasm that has just passed through a therapeutic bottleneck may be a sign that therapy selected for a mutator phenotype40. Because of this complication, we agreed that we must measure how neoplasms are changing over time as well as diversity.

Change over time. There are a variety of ways that a neoplastic cell population changes over time. These include mutations, natural selection and genetic drift. One important parameter of change over time is the mutation rate, which describes how fast a lineage accumulates new mutations. Of course, there are different mutation rates induced by each mechanism for genetic and epigenetic alteration, including mutation signatures induced by specific agents41 as well as telomere erosion, non-homologous recombination, other forms of chromosomal instability, CpG methylation and histone modifications. Which mechanisms are relevant will depend on individual tumours and may vary across the different clones within the same tumour.

When we talk about and measure mutation rates, we are implicitly assuming that mutations happen at a regular rate. Evolutionary biologists call these 'molecular clocks' (Ref. 42). However, a catastrophic mitosis can generate chromosomal alterations across the genome in a single event43,44. There is a continuum from regular, gradual, clock-like small alterations to sporadic, punctuated, large alterations. For example, a lineage may evolve different mutation rates across its history, as happens with the evolution of a mutator phenotype45,46. If a cell lineage can change suddenly, in what used to be called a 'macromutation' generating a 'hopeful monster' (Ref. 47), then that tumour may have a different capacity for evolution compared with a tumour that is constrained to evolve through the slow accumulation of mutations with small phenotypic effects. There is a large cancer literature on genetic instability that is relevant here48,49, and evidence has shown that tumours with extremely high mutation rates may have a better prognosis than tumours with moderate rates6,11,21,22,50. High levels of genomic instability may make it difficult for cell lineages to maintain the adaptive information encoded in their genomes, generating non-viable daughter cells, and may also produce an abundance of neo-antigens that stimulate an antitumour immune response6. Furthermore, high mutation rates of single nucleotide variants can generate deleterious mutations, leading to the fitness decline of neoplastic cell lineages in a form of Muller's ratchet21,51. This may even cause tumour regression in some cases21,22.

The genetic composition of a population changes over time not only through the rate at which mutations arise and the genetic drift of those alleles but also through the action of natural selection. Natural selection leads to adaptations, such as drug resistance52, that are clinically relevant. Detecting and measuring natural selection is likely to be an important component of our future clinical management of cancers.

The classification of a neoplasm's change over time (Δ) will probably need to take into account both the speed at which a tumour acquires genetic or epigenetic alterations, or changes phenotypically, including how fast clones spread by natural selection, as well as the tempo of that change (from gradual to punctuated). The appropriate intervals for longitudinal sampling will depend on the rate of change over time53. Note that neutral, or 'passenger', mutations should not be ignored in these calculations because selective pressures change over time, particularly with the onset of therapy. Thus, resistance mutations, which may be deleterious or neutral in the absence of therapy, can become selectively advantageous for neoplastic cells exposed to therapy54.

Measuring change over time. Measuring change over time is complicated, whether it is genetic or phenotypic change (Table 1). Figure 1b illustrates a simple version of how the Evo-index can describe evolutionary changes in tumour cell populations. It is possible for there to be change over time but for diversity to remain stable, with a dynamic equilibrium of clones appearing and going extinct14. For single samples, past genetic changes over time can be indirectly inferred based on mutation frequencies17,55. Sottoriva and Graham have pioneered methods to infer the mutation rate and to distinguish between tumours that are dominated by genetic drift versus those with evidence of natural selection after transformation. In the absence of selection, mutations that occur in the first cell division after transformation should appear in approximately one-half of all cancer cells, mutations that occur in the second round of cell division should appear in one-quarter of all cancer cells, and so on17,56.

There are a number of measures of genetic change over time from population genetics that might be used on neoplasms, including Nei's standard genetic distance57,58 and the Jaccard similarity coefficient59, as well as measures of beta diversity that can also quantify changes in a community over time, such as UniFrac34 or the fixation index60. The degree of genetic divergence between samples (called 'nucleotide diversity' in molecular population genetics) provides indirect information on the degree of change over time. Genetic divergence is often defined as the percentage of the genome that is different between pairs of samples12,13,14. This statistic provides predictive power independent of the number of clones for predicting progression12,13, supporting the framework of including both diversity and change over time in the Evo-index. Note that the same clonal structure can have radically different degrees of genetic divergence (Fig. 2). Maley and colleagues have calculated a mean pairwise divergence score between all pairs of samples from a neoplasm12,13,14. As the chance that two samples come from the same clone (and so have minimal divergence) depends on the size of the clone, the mean pairwise divergence blends the degree of divergence with clone size measures (and so blends D with Δ).

Figure 2: Clonal divergence is independent of clonal structure.
figure 2

The cell lineages from two tumours may have the exact same clonal structure when they are sampled at the far right but have radically different degrees of genetic divergence. If one tumour (part a) has a higher mutation rate or has been accumulating genetic alterations for a longer period of time because those cells had a common ancestor, it will have a higher level of genetic divergence than another tumour (part b).

PowerPoint slide

One of the primary tools for measuring change over time in evolutionary biology is phylogenetic inference, which reconstructs the history of a neoplasm61,62. Phylogenetic methods can be used to describe and quantify diversity patterns as well as rates of evolution across both space and time. Multiple phylogenetic approaches have been developed in recent years to study tumour evolution within a patient, both for bulk and single-cell data and from a variety of data types20,63. These methods depend on evolutionary models for the likelihood of molecular alterations occurring in neoplastic cell lineages, although the development of these models is still in its infancy.

All of the measures discussed so far can be calculated from a single timepoint. Of course, the degree and nature of change over time can be better measured directly with longitudinal samples. Minimally invasive assays, such as sequencing cfDNA from longitudinal blood samples, could reveal the action of natural or artificial selection in patients.

Incorporation of the Evo-index into clinical trials can better describe, in evolutionary terms, why interventions fail. Most human tumours at the time of clinical presentation contain multiple large clones6,16 and probably many more small clones64,65, and relapse without a reduction in diversity would probably imply intrinsic resistance or perhaps that an intervention resulted in increased mutagenesis. By contrast, relapse with less diversity (D1) implies a bottleneck effect where only a minority of tumour cells survived the intervention, probably indicating selection for one or a few resistant clones.

The Eco-index

From the perspective of an organism or a neoplastic cell, its ecology can be broadly described by two characteristics: hazards (H) and resources (R)66,67,68,69 (Fig. 3). Hazards, here, are the things that can kill a cell. The relevant resources required for cell maintenance and growth are many and varied; whatever may potentially limit the growth of the neoplastic cell population66. Note that hazards and resources here are understood from the perspective of the neoplastic cell, not the patient. This is an important point from ecology — we can understand the evolution and responses of a population best when we take the perspective of an organism in that population70.

Figure 3: The Eco-index.
figure 3

The ecological index (Eco-index) is composed of two factors corresponding to the hazards (H) and resources (R) available to the neoplastic cells. These capture the broad categories of selective pressures on a population. We have included example phenomena in this figure that might be observed in the different combinations of the degrees of hazards and resources. For example, a tumour with low hazards (H1) and low resources (R1) might be relatively barren, with few infiltrating lymphocytes but also poor perfusion and few supporting cells. Such an environment would select for cells that can either survive on few resources or move to locate more resources. High levels of hazards (H2) should, according to life history theory71, select for rapid proliferation, evasion of predation, migration away from the hazards67 and little investment in cell (and DNA) maintenance. High levels of resources allow neoplastic cells to rapidly proliferate. Thus, an H2R2 tumour would probably undergo massive cell turnover as cells are killed by the hazards and replaced by their rapidly proliferating sisters.

PowerPoint slide

From an ecological perspective, the hazard and resource profiles for a species select for the particular life history strategies of that species. Aktipis and colleagues argued that the same principles are true for neoplastic cells71. Species that are exposed to high levels of hazard tend to evolve fast life history strategies, reproducing quickly and investing little in maintenance and survival. Organisms subjected to hazards generally leave behind higher levels of unexploited resources. Ecosystems with high or fluctuating resource supplies favour organisms that can rapidly reproduce to exploit those opportunities. This selects for speed over efficiency and can result in very high population densities but also fluctuating levels of unexploited resources. By contrast, populations that have few hazards and a steady supply of resources will tend to expand to the carrying capacity of the habitat, at which point natural selection favours organisms that can best compete for and efficiently utilize the limiting resources72. The heterogeneity of resources and hazards across space also has important impacts on the future evolution of cancer cell populations and prognosis for patients73,74.

Hazards. There are multiple sources of hazards for neoplastic cells, including immune cells, toxins, waste products, microorganisms and anticancer therapies. There is good evidence that immune predation is associated with improved cancer prognosis73,75,76,77,78,79,80,81,82,83. Furthermore, there is emerging evidence linking high mutation loads that result in the formation of neo-antigens with immune predation and better survival in patients treated with immune checkpoint blockade therapies84,85,86. In addition, a high subclonal neo-antigenic burden is associated with worse outcomes in lung cancer when patients are treated with checkpoint inhibitors87. These data suggest that subclonal neo-antigens might impede cytotoxic immune responses against neo-antigens that are present in every tumour cell.

Other hazards faced by neoplastic cells include the accumulation of waste products in their microenvironments67,69,88,89. This may include lactic acid and lactate build-up from glycolysis88,90 as well as reactive oxygen species from excessive cellular proliferation91. Methylglyoxal92,93, nitric oxide94,95 and advanced glycation end products96,97 have also been implicated as toxic waste products in cancer microenvironments.

The role of the microbiome in cancer is complicated and largely unknown. While some microorganisms may promote tumours98,99, others have antitumour effects98, enhancing the efficacy of chemotherapy100. Thus, microorganisms may act as both resources and hazards for neoplastic cells.

Measuring hazards. The current best measures of hazards for a neoplastic cell depend on measures of immune predation (Table 1). There is a large literature on the association between infiltrating lymphocytes and favourable prognosis in cancer73,75,76,77,78,79,80,81,82,83. In addition, a pan-cancer analysis revealed T cell signatures to be broadly favourable prognostic markers across 25 cancer types101. Galon and colleagues have found that a signature of activated T cells from bulk tumour samples is also strongly predictive of favourable survival76,77,78,83. Yuan and colleagues have shown that haematoxylin and eosin images can be computationally analysed to identify neoplastic cells, fibroblasts and lymphocytes and, furthermore, that patients with breast cancer who show colocalization of neoplastic cells with lymphocytes in the tumour have a better prognosis than patients with tumours in which the lymphocytes are separated from the neoplastic cells75. This is based on a standard ecological statistic, the Morisita–Horn index102, for quantifying statistically significant colocalization in order to detect ecological interactions (in this case, predation). These results suggest that immune predation is a major form of hazard for a neoplastic cell, and measures of that predation should be a central component of the ecological index.

While much research has investigated the potentially toxic effects of low pH (Refs 103,104), fewer studies have examined the fitness consequences to cancer cells from various metabolites. Future research should determine the effects of different concentrations of putative toxic metabolites on cancer cell survival and proliferation in both cell culture experiments and mouse models. Measurements of anticancer drug concentrations in the tumour are also likely to quantify important hazards for the neoplastic cells. In addition, the microbiome (including the virome) of tumours can be surveyed to reveal microbial hazards for the neoplastic cells105.

Resources. Resources, including oxygen, glucose, micronutrients, survival signals, growth signals and space, are also critical to the future behaviour of a tumour. Surprisingly little is known about the interactions between cell metabolism and the availability of key resources, which ecologists term the organism's 'foraging ecology'. Almost all cancers rely on glycolytic as opposed to aerobic metabolism, suggesting that resources can select for tumour phenotypes106,107. From nature, we know that selection favours feeding behaviours that balance speed, efficiency and safety108. There must be strong selection for cancer cells to do the same (for example, through upregulation of transporters such as glucose transporter type 1, erythrocyte/brain (GLUT1, also known as SLC2A1)109). Measuring which resources limit the population size and proliferation of neoplastic cells would allow researchers to identify some of the strongest selective pressures on the tumour and to predict how it will change in the future. This approach would also provide targets for further reducing the evolvability of the neoplasm by lowering the carrying capacity of its microenvironment.

In the broader ecological literature, consumer–resource theory110 shows that resource supply, depletion and availability affect population growth rates, population sizes and competition between different species (that is, distinct clonal lineages). Resource supply represents the rate at which new resources enter the system (in this case, the tumour) and the rate at which resources become available through nutrient cycling within the system. The aggregate consumption of all cells depletes the resources, typically to levels much lower than experienced by normal tissues111. In fact, glucose becomes depleted below levels detectable by most analyses112. However, in some cases, immune predation and fluctuations in resource supply can prevent the complete exploitation of resources113,114, leaving patches of residual resources available for future exploitation115.

The potential resources for a tumour include the contents of plasma and the metabolites synthesized and secreted by the normal cells of the tumour and its microenvironment. Hence, the list includes proteins (albumins, globulins and fibrinogens), glucose, amino acids, fatty acids, hormones, electrolytes, oxygen and trace elements. The functional response and the value of the resources to the consumer are dictated by nutritional relationships116. In some cases, lack of a resource may trigger stasis, but in others, it may lead to cell death or dispersal117. At the moment, there are many open questions about the intratumoural cycle of critical nutrients other than carbon and nitrogen (that is, phosphate, iron, copper, etc.)118. These nutrient cycles may contain valuable therapeutic targets.

Some resources, particularly growth and survival signals, may be provided by the neighbouring stromal cells119,120. Nutrients may also be provided by the stroma. Pyruvate and lactate can be supplied to cancer cells by activated fibroblasts121,122, and fatty acids may be supplied by activated adipocytes123,124. Tumour and stroma only come into physical contact when the basement membrane is breached by malignant neoplastic cells. At this stage, cancer cells can directly interact with cancer-associated fibroblasts, which are known to play a key role in the regulation and development of tumours, especially solid tumours120,125. In this secretory reactive state, fibroblasts facilitate not only cancer growth and progression126,127 but also treatment resistance128. In addition, their presence in a tumour has been correlated with poor outcomes129.

Other resources must be delivered through the vasculature. Folkman made the crucial link between angiogenesis and tumour invasion and metastasis, realizing that preventing new vessels from forming could be a simple way to inhibit further tumour growth130,131. The presence in many tumours of necrosis and hypoxia, which are major drivers of angiogenesis, attests to the importance of resource limitation in tumours. Furthermore, there is evidence that necrosis is a prognostic factor in many cancers132.

The effects of resources on the evolution of a tumour are not defined simply by their supply, depletion and availability. Resource diversity may also be important. Whether resources are uniform across space or heterogeneous ('patchy' or exhibiting gradients) makes a difference67,133. Patchy resources (and hazards) create multiple habitats (for example, rich and sparse regions) that may select for different clones that can survive in those regions and may be differentially responsive to (and differentially exposed to) therapies. Furthermore, we and others have shown that if those patchy resources change over time, then there is selective pressure on cells to move to escape regions of scarce resources and exploit transient regions of more plentiful resources67,113,114,134,135,136. Thus, ecological theory predicts that heterogeneous resources should select for invasion and metastasis134,135, and there is evidence to support that prediction in cancer137,138,139,140,141,142,143. Verduzco and colleagues found that intermittent exposure of some cell lines to hypoxia selected for increased resistance to a variety of chemotherapies, including etoposide, docetaxel and methotrexate, compared with unselected controls144. In addition, resource gradients often lead to rapid evolution, as organisms that are able to invade more stressful environments can escape competition and flourish145. Much needs to be learned about resource heterogeneity, consumer–resource dynamics and the foraging ecology of neoplastic cells.

Measuring resources. Measuring resources (and hazards) requires the consideration of relevant spatial and temporal scales. It is not yet clear how to combine measures of the level of resources, their spatial variance and their stability over time into a single statistic.

There are various resources and methods to measure them that may be prognostically relevant (Table 1). The proportion of a tumour that is necrotic or poorly perfused may be read from standard positron emission tomography and computed tomography (PET–CT) images146 and through other measures of blood vessel density147,148. The degree and patchiness of hypoxia can also be assayed in FFPE samples with antibodies against carbonic anhydrase 9 (CA9) or hypoxia-inducible factor 1α (HIF1α)115 or via intravenous introduction of 2-(2-nitro-1-H-imidazol-1-yl)-N-(2,2,3,3,3-pentafluoropropyl) acetamide (EF5) and the subsequent measurement of its binding in the tumour tissue149. EF5 binding and related techniques have proved useful in the clinic for detecting regions of hypoxia, determining prognosis and measuring response to therapy150. While it is difficult to measure glucose concentration directly, an indirect measure may be made via immunohistochemistry staining for expression of GLUT1115. Measures of ATP may also be a good indirect measure of the amount of resources available to neoplastic cells151. Glutamine, pyruvate, lactate, fatty acids, calcium, potassium, phosphorus and various trace metals may also be limiting and important to measure, but this appears to be unexplored. Most of these measures will be limited to biopsy samples analysed ex vivo and thus will suffer the problems of spatial heterogeneity and sampling error.

In some cases, the problem of spatial heterogeneity and sampling error can be avoided through gross measures of resources from radiological images152,153,154 Radiographic images such as those obtained using PET–CT and magnetic resonance imaging (MRI) can provide valuable habitat data. In natural systems, there is usually a tight correlation between habitat and the types and characteristics of species inhabiting the habitat. Similarly, simply knowing the different habitat types within a tumour may be prognostic of the community of cancer cells and therapeutic outcomes. For instance, in glioblastoma, measures of fluid-attenuated inversion recovery (FLAIR), T1 and T2 from MRI examinations after gadolinium administration identified distinct habitats that correlated with therapeutic outcome, independent of tumour size153. Texture analysis of MRI scans has been used to identify spatial heterogeneity and regional variations that are associated with microenvironmental conditions, including cell density, tissue stiffness, blood flow and nutrient dispersion152,154. These may also be used to measure functional diversity (D) in tumours. Geographic information systems (GIS)155,156,157 and ecology158 provide a rich literature and a source of tools for analysing spatial resource information, but these are rarely utilized in cancer research73,74.

Standard histopathology can provide measures of T cell infiltration and vascular and lymphatic density77. Using digital pathology, Lloyd et al. investigated the spatial distributions of oestrogen receptor (ER) expression in relation to vascular density and tissue necrosis in breast cancer histology specimens, revealing considerable regional variations in cancer proliferation phenotypes accompanied by vascularity and immune response115,159. Yuan and colleagues also used digital pathology to analyse the spatial relationships between fibroblasts and neoplastic cells160. We have summarized the statistics and assays for measuring diversity, change over time, hazards and resources in Table 1.

Categories of tumours

The future behaviour of a tumour depends on both its evolutionary potential (the Evo-index) and the selective pressures on the tumour (the Eco-index). A highly evolvable tumour may or may not evolve immune evasion depending on whether the immune system is imposing a strong selective pressure on the tumour. By contrast, an immune response may or may not lead to immune evasion depending on the evolvability of the tumour. Thus, both the evolution and ecology of a tumour must be considered in predicting cancer outcomes. We therefore propose to combine the Evo- and Eco-indices to classify tumours. Dichotomizing each evolutionary and ecological factor of the Evo- and Eco-indices into high and low values would produce 16 possible types of tumour (Table 2).

Table 2 An initial classification scheme

In order to classify a tumour, investigators will first need to define and validate clinically relevant thresholds for dichotomizing diversity, change over time, hazards and resources (Table 1). For example, in Barrett oesophagus, Maley and colleagues found that the upper quartile of diversity statistics distinguished patients who are likely to progress to oesophageal adenocarcinoma12,13,14. Once those thresholds are validated, a tumour would be measured for each of the four evolutionary and ecological factors to determine which of the 16 types it falls into. For example, if a tumour was below the thresholds for all four factors (that is, a D1Δ1H1R1 tumour), it would be a type 1 tumour.

A roadmap for improvements. We are not yet in a position to specify which measures and thresholds should be used to determine the D#Δ# or H#R# type of a tumour. Initial studies should test if these classifications significantly predict clinical outcomes and which evolutionary and ecological measures provide independent predictive value. They should also test if there are measures that can apply across cancer types or if they have to be uniquely defined for specific organs or tumour subtypes. Future studies should test alternative measures of diversity, change over time (Box 2), hazards and resources to help standardize useful metrics for the classifications. They should also quantify the improvements to prognosis gained by sampling multiple regions at multiple timepoints.

The ecology of a tumour affects its evolution, and the evolution of the cells in a tumour change their ecology. Neoplastic cells evolve genomic instability161, generating neo-antigens as well as adaptations, such as recruitment of resources, through activating fibroblasts162 and neo-angiogenesis161. Evolution of neo-antigens triggers immune predation, which may reduce diversity and select for immune evasion163. High levels of extrinsic mortality and resources select for rapid proliferation with little investment in somatic maintenance71. These interactions imply that not all possible combinations of ecological and evolutionary measurements are equally likely. We will probably be able to drop some of the 16 possible tumour types in Table 2 and focus on the subset of classes that present in the clinic.

The framework for a classification system that we have proposed could be incorporated into clinical trials, which could allow us to gather data on how the different types of evolving tumour respond to different types of intervention (Fig. 4). Clinical trials could then be developed to stratify treatment of patients based on the Evo- and Eco-indices of their tumours. We could use the results to develop guidelines for best practice in managing cancers.

Figure 4: Changing the evolutionary class of a tumour through interventions.
figure 4

With the classification system outlined in Table 2, we could examine how different interventions move tumours between categories. a | In this example, chemotherapy can be mutagenic and can select for hypermutator clones, generating new clones and more diversity40,186,187. It can also kill endothelial cells and thus have an anti-angiogenic effect188, resulting in a tumour (type 13) with one of the worst predicted prognoses. This may partly explain why tumours that recur after chemotherapy are so difficult to control. b | Immunotherapy, if successful, may increase the predation hazards to the tumour and perhaps select for a subclone, reducing diversity. Targeted therapy, unlike chemotherapy, probably does not cause significant DNA damage and may further genetically homogenize the tumour. Anti-angiogenic therapy is designed to restrict the resources of the tumour. At the end of this example sequence, the tumour is in the most manageable, least evolvable category (type 3 in Table 2). Of course, chemotherapy, immunotherapy and targeted therapy may have different effects depending on the details of those therapies and their interaction with the clones in the tumour and their ecosystem.

PowerPoint slide

Vision of the future

In the future, the pathology report for a neoplasm could include its Evo-index and Eco-index classifications. Ideally, these classifications would provide 'chessboard'-like scenarios where, based on the current evolutionary class of a tumour, one could anticipate how the tumour type will change with different possible therapeutic moves (Fig. 4). Clinicians would then be able to choose appropriate interventions for the evolvability of those neoplasms and would also be able to track whether the neoplasms change substantially in response to interventions. A D1Δ1 tumour or even a D1Δ2 tumour would be a prime candidate for aggressive therapy with curative intent. In fact, a D1Δ1 tumour may be so evolutionarily indolent as to not require any form of intervention. On the other hand, a D2Δ2 tumour is likely to have multiple resistant subclones present at diagnosis, and future clinical trials should test if such a tumour can be managed through strategies that minimize the expansion of resistant subclones by exploiting their disadvantage in competition with sensitive subclones164. A legitimate clinical strategy might be to down-stage a tumour from a highly evolvable one to a much more clinically manageable class that could be contained in a non-lethal state indefinitely (Fig. 4b). If validated, the Evo- and Eco-indices could be used as surrogate measures for overall survival or disease-free survival.

Conclusions

The evolutionary biology of cancer is, clinically, in a similar state to psychiatry in the nineteenth century. At that time, there was no standard classification system for mental illness used by practitioners. Without such a classification system, it was difficult to even talk about the illness, let alone make progress, as a common language was lacking. With the American Medical Association's Standard Classified Nomenclature of Disease published in 1933 (Ref. 165) and the first Diagnostic and Statistical Manual of Mental Disorders published in 1952 (Ref. 166), no matter how flawed they were, diagnoses of mental disorders became standardized, which facilitated studies to refine both the classifications as well as the treatment of those disorders. Studies based on the same classification system were then comparable, which further facilitated meta-analyses and overall progress in the field.

We have diagnostic categories for types of tumour based on their tissue of origin and staging, as well as some molecular markers, but we have lacked a system for classifying the evolvability and ecology of a tumour, which help determine how it will respond to interventions and how it might best be managed. Evolutionary oncology requires a shared lexicon upon which to base discovery. We reached consensus on the proposed framework for a classification system to characterize evolutionary differences between tumours that is applicable across all cancer types. Importantly, an evolutionary classification system will facilitate future efforts to study this fundamental property of tumours to reveal implications for treatment.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.