# Group phenotypic composition in cancer

Jean-Pascal Capp, James Degregori, Aurora M Nedelcu, Antoine M Dujon, Justine Boutry, Pascal Pujol, Catherine Alix-Panabières, Rodrigo Hamede, Benjamin Roche, Beata Ujvari, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Jean-Pascal Capp, James Degregori, Aurora M Nedelcu, Antoine M Dujon, Justine Boutry, et al. Group phenotypic composition in cancer. eLife, 2021, 10, pp.e63518. 10.7554/eLife.63518. hal-03261102

HAL Id: hal-03261102

https://hal.inrae.fr/hal-03261102

Submitted on 15 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

# Group phenotypic composition in cancer

Jean-Pascal Capp¹†, James DeGregori²†, Aurora M Nedelcu³†, Antoine M Dujon⁴,⁵, Justine Boutry⁴, Pascal Pujol⁴, Catherine Alix-Panabières⁴,⁶, Rodrigo Hamede⁷, Benjamin Roche⁴, Beata Ujvari⁵,⁷†, Andriy Marusyk⁸†, Robert Gatenby⁸†, Frédéric Thomas⁴*

¹Toulouse Biotechnology Institute, University of Toulouse, INSA, CNRS, INRAE, Toulouse, France; ²Department of Biochemistry and Molecular Genetics, University of Colorado Anschutz Medical Campus, Aurora, United States; ³Department of Biology, University of New Brunswick, Fredericton, New Brunswick, Canada; ⁴CREEC/CANECEV, MIVEGEC (CREES), University of Montpellier, CNRS, IRD, Montpellier, France; ⁵Centre for Integrative Ecology, School of Life and Environmental Sciences, Deakin University, Geelong, Australia; ⁶Laboratory of Rare Human Circulating Cells (LCCRH), University Medical Centre of Montpellier, Montpellier, France; ⁷School of Natural Sciences, University of Tasmania, Hobart, Australia; ⁸Department of Cancer Physiology, H Lee Moffitt Cancer Center and Research Institute, Tampa, United States

**Abstract** Although individual cancer cells are generally considered the Darwinian units of selection in malignant populations, they frequently act as members of groups where fitness of the group cannot be reduced to the average fitness of individual group members. A growing body of studies reveals limitations of reductionist approaches to explaining biological and clinical observations. For example, induction of angiogenesis, inhibition of the immune system, and niche engineering through environmental acidification and/or remodeling of extracellular matrix cannot be achieved by single tumor cells and require collective actions of groups of cells. Success or failure of such group activities depends on the phenotypic makeup of the individual group members. Conversely, these group activities affect the fitness of individual members of the group, ultimately affecting the composition of the group. This phenomenon, where phenotypic makeup of individual group members impacts the fitness of both members and groups, has been captured in the term 'group phenotypic composition' (GPC). We provide examples where considerations of GPC could help in understanding the evolution and clinical progression of cancers and argue that use of the GPC framework can facilitate new insights into cancer biology and assist with the development of new therapeutic strategies.

## Introduction

Primary cancers often display high levels of phenotypic and genetic intratumor heterogeneity (ITH). Whereas the mutational and epigenetic mechanisms that generate this heterogeneity are relatively well understood, mechanisms responsible for maintaining this heterogeneity in cancer cell populations, as well as the impact of ITH on clinically important properties of the disease, are less clear. Yet, high levels of ITH have been linked to more aggressive tumor behavior, resistance to therapies, and overall poor prognosis (Marusyk and Polyak, 2010; Meacham and Morrison, 2013; McGranahan and Swanton, 2015; Morris et al., 2016; McGranahan et al., 2017; Dagogo-Jack and Shaw, 2018; Hausser and Alon, 2020; Marusyk et al., 2020). Thus, a deeper understanding of ITH at different stages during cancer progression might result in novel therapeutic strategies directed at altering tumor composition towards decreased malignancy.

<footnote>*For correspondence: frederic.thomas2@ird.fr</footnote>

<footnote>†These authors contributed equally to this work</footnote>

<footnote>**Competing interests:** The authors declare that no competing interests exist.</footnote>

<footnote>**Funding:** See page 14</footnote>

<footnote>**Received:** 02 October 2020</footnote>

<footnote>**Accepted:** 09 March 2021</footnote>

<footnote>**Published:** 30 March 2021</footnote>

<footnote>**Reviewing editor:** Antonis Rokas, Vanderbilt University, United States</footnote>

<footnote>© Copyright Capp et al. This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.</footnote>
ITH can be subdivided into micro-heterogeneity (i.e., cell-cell variability) as well as macro-heterogeneity (i.e., heterogeneous subgroups), yielding a co-existence of either independent genetic subclonal lineages or metastable phenotypic subgroups (Merlo et al., 2006). Despite the rapidly growing interest in assessing ITH, the complex eco-evolutionary dynamics between cancer cells and their environment, the balance between heritable and non-heritable cellular traits, as well as the resulting effects on global tumor behavior remain poorly understood.

In particular, stable co-existence of distinct subpopulations or subclones within tumor cell populations that follow Darwinian dynamics appears to be counterintuitive. Indeed, the maintenance of stable subpopulations with heritable phenotypes within the same tumor indicates the presence of multiple niches producing complex patterns of co-existence, although weak stabilizing selection and drift may also contribute to phenotypic variability. In turn, the co-existence of different subpopulations can profoundly influence global tumor behavior such that elimination of one subpopulation can induce tumor regression or necrotic collapse (Inda et al., 2010; Marusyk et al., 2014). This co-existence can be based on 'cooperation' (or 'mutualism' in ecological terms) between distinct subpopulations (Merlo et al., 2006; Axelrod et al., 2006; Tabassum and Polyak, 2015; Zhou et al., 2017; Martín-Pardillos et al., 2019). Synergistic interactions between two subpopulations of cancer cells may also permit metastatic spread (Marusyk et al., 2014; Calbo et al., 2011; Neelakantan et al., 2017; Janiszewska et al., 2019). Importantly, reversible phenotype switching seems to have a prominent role in the dynamics of interactions between diverse subclones (Chapman et al., 2014), suggesting that the intratumor phenotypic heterogeneity arising from non-genetic variability (i.e., phenotypic plasticity) permits a rapid, reversible strategy for cancer cells to optimally respond to changing environmental conditions and population fluctuations. In addition, populations of tumor cells also exhibit Allee effects (i.e., positive correlation between population density and individual fitness) in which the proliferation rate paradoxically increases with population size, which promotes aggressive tumor growth (see for instance Korolev et al., 2014; Böttger et al., 2015; Sewalt et al., 2016; Johnson et al., 2019). However, at the moment, investigations on these interactions are limited to few subpopulations (Marusyk et al., 2014; Calbo et al., 2011), molecular interactions (Inda et al., 2010; Vinci et al., 2018; Wu et al., 2010; Cleary et al., 2014), and cancer hallmarks (Chapman et al., 2014; Wagenblast et al., 2015). Moreover, these examples represent a snapshot of tumor progression, while the interactions between subpopulations can be highly dynamic and change over time. To capture these dynamics, we suggest applying an ecological concept termed 'group phenotypic composition' (GPC), to understand the ecology/evolutionary dynamics of cancer.

In nature, the function and stability of a group is dependent on the characteristics of the constituent individuals (and/or subgroups composed of individuals sharing common features) and their interactions, which are always in an evolutionary state of flux. Tumors have already been described as composed of functional compartments with a division of labor among these compartments (Grunewald et al., 2011; Hausser et al., 2019) and conceived as a whole consortium of cooperating malignant clones (Ramón Y Cajal et al., 2017). Other studies suggested that phenotypic groups can be defined among tumoral cells on the basis of common expression of phenotypic markers, even if these phenotypic clusters can display genetic and epigenetic heterogeneity (see Box 1). Here, we propose to extend the concept of GPC toward tumor cell populations and investigate how the phenotypic composition of tumors can impact cancer progression and treatment.

## GPC in the animal kingdom and its relevance to cancer

Group living is a widespread phenomenon within the animal kingdom and has attracted considerable attention among evolutionary ecologists (Farine et al., 2015; Koenig et al., 2020). Recent advances in this area have highlighted that GPC is a crucial parameter to consider. Farine et al., 2015 define GPC as "any descriptor of the types of phenotypes found within a group. Examples include the average body size of individuals, the variation in male color, or the aggressiveness level of the most aggressive individual in the group. In this framework, 'groups' can represent many aspects of the social environments of individuals, including breeding units, social networks, neighborhoods, populations, and communities." Members in any group, be it individuals in populations or species in communities, express some phenotypic variation at several levels (e.g., activity, behavior,
Box 1. Intratumoral phenotypic clustering despite genetic and epigenetic diversity.

Cancer cells can reversibly and stochastically transit between states that differ in their competence to contribute to tumor pro-
gression (Gupta et al., 2011). These stochastic cell-state transitions can not only generate cells with the properties of cancer
stem cells (CSC), but also make cells transdifferentiate without acquisition of a stem-like state. This cancer cell plasticity can con-
tribute to tumor initiation, progression, and therapeutic resistance (Capp, 2019; Gupta et al., 2011; Pisco and Huang, 2015)
and is expected to greatly contribute to cell-to-cell heterogeneity in cancer (Meacham and Morrison, 2013; Marusyk et al.,
2020). For instance, transdifferentiation can be a major driver of drug resistance and metastasis through a mechanism that ulti-
mately relies on epigenetic mechanisms (Sharma et al., 2018). The ability to transit, especially toward a CSC state, is dependent
on either intracellular stochastic events or microenvironmental perturbations (Nguyen et al., 2012; Beck and Blanpain, 2013).
Indeed, various phenomena such as epigenetic regulation, stochastic gene expression, or variability in the microenvironment can
contribute to phenotypic diversification (Marusyk et al., 2012). Among them, histone modifications were shown to confer line-
age plasticity in cancer (Pastore et al., 2019; Liau et al., 2017; Flavahan et al., 2017). For instance, intratumoral epigenetic
diversity allows leukemic cells to stochastically activate alternative gene regulatory programs, facilitating the emergence of novel
cell states, ultimately enabling tumor growth and drug tolerance (Pastore et al., 2019). The strong increase in entropy of the
cancer epigenome appears to increase stochasticity of gene expression at higher levels than normal stem cells because of a less
organized and less stable chromatin structure (Jenkinson et al., 2017), thus greatly impacting phenotypic variability.
Interestingly, when a subpopulation is purified for a given phenotypic state, a return toward equilibrium proportions is then
observed over time (Gupta et al., 2011), suggesting that such phenotypic transitions occur regardless of the cell's genetic con-
tent. A pioneering study already showed that the distinct cell populations and the underlying transcriptional heterogeneity
observed among single human colon cancer cells were not due to underlying genetic heterogeneity, as injection of single CSC
into immune-deficient mice gave rise to monoclonal tumors as heterogeneous as the parental one (Dalerba et al., 2011). To
generalize these observations, it is now recognized that phenotypic clusters (groups of cells expressing the same phenotypic
markers) exist in vivo among tumor cells (Wagner et al., 2019; Jackson et al., 2020), despite divergent genomes, epigenomes,
and clonal origins. The same functional clustering has been suggested for breast cancers despite high intertumoral heterogene-
ity: diverse genetic and epigenetic alterations converge phenotypically into four main breast cancer classes (Cancer Genome
Atlas Network, 2012).

Innovative methods such as single-cell mass cytometry allow the detection of phenotypic clusters among single tumoral cells despite extensive heterogeneity in genetic and epigenetic content. For instance, a consensus clustering approach allowed hierarchical clustering of 45 epithelial clusters in breast cancers based on normalized epithelial cell marker expression (Wagner et al., 2019). Recently, by designing an imaging mass cytometry panel specific to breast histology, it was possible to quantify spatial inter- and intratumor phenotypic heterogeneity at the single-cell level among 281 tumors using 35 antibodies (Jackson et al., 2020). This strategy allowed mapping the cellular spatial organization of these tumors and to observe variable structures and cellular densities. It also revealed relationships between cellular phenotype, tissue organization, and clinical outcome. In summary, both heterogeneous tumors with multiple phenotypically pure communities and homogeneous tumors with one epithelial sheet or with similar communities of different sizes were observed, the former being associated with poorer outcomes. The authors concluded that studying how spatial and phenotypic tissue features influence disease outcome might be medically relevant.

morphology, life-history traits, etc.). The heterogeneous composition of groups provides them with different properties, due to the different average phenotypes, the presence of keystone individuals, and/or the variability in phenotypes etc., that can in return influence group-level outcomes (e.g., foraging success, mating system, predation risk, evolvability). While theory on the implications of GPC is still in its infancy, the concept proved relevant to understand many aspects of group-level dynamics and self-organization (Farine et al., 2015). In addition, selective pressures arising from GPC properties can influence the evolutionary trajectory of individual phenotypes. For instance, large groups composed of individuals with similar phenotypes may benefit from a reduced predation pressure because of the confusion effect that the lack of prey 'oddity' generates on predators (Landeau and Terborgh, 1986; Krakauer, 2004). In prey species such as the cladoceran Daphnia magna for which predation risk (by the three-spine stickleback Gasterosteus aculeatus) is higher on larger individuals (i.e., more profitable preys) and confusion effects operate, some conflicts between individuals with different phenotypes might emerge: both large and small individuals benefit from
<figure><img src="image_4.png" /><figcaption>**Figure 1.** GPC and its ecological/evolutionary implications. Because individuals show significant phenotypic variation (I), groups may vary in their GPC depending on which individuals constitute the group (II). This may in return influence group-level outcomes (III), which as a consequence differentially impacts individual fitness (IV). (V) GPC then drives different evolutionary responses, like for instance decisions to join or leave particular groups or phenotypic plasticity in response to GPC. These phenomena can then influence the distribution of phenotypes in subsequent generations (figure modified from Farine et al., 2015).</figcaption></figure>

joining groups composed of large individuals, but the presence of small individuals enhances the predation risk experienced by large ones (Rodgers et al., 2015). Selection from GPC can thus alter the covariance between individual and group phenotypes through the removal of particular phenotypes within generations. In the *Daphnia/stickleback* system, excluding or avoiding small individuals would be beneficial for large individuals, but it also comes with (energetic) costs, leading to mixed groups in the wild. This example illustrates the kind of complex and reciprocal links that may exist between individuals’ phenotypes and the GPC’s properties. More generally, individual heterogeneity and resulting GPCs result in multiple hierarchical effects that affect collective behavior, from within-group positioning, group coherence, leadership, and collective decision making to group functioning, fission-fusion dynamics, and among group assortment (see Jolles et al., 2020 for a recent review).

Farine et al., 2015 provided a review with many examples as well as a conceptual framework (Figure 1) allowing one to understand and predict when GPC should have important selective consequences. GPC in the animal kingdom has both pervasive consequences for individual fitness and fascinating evolutionary implications, with potential to contribute to many classical questions in ecology and evolution.

In so far as tumors are heterogeneous subgroups of cells that interact among each other and with other living entities, for example host normal cells with the microbiota, it is theoretically possible to examine them within the framework of the GPC. Noticeable ecological and evolutionary singularities however exist, given that most cancers are not transmissible and hence cannot evolve for longer periods than the life expectancy of their host, and because tumors are usually not in competition inside the hosts (except in transmissible tumors like DFTD and DFT2 in Tasmanian devils, see
<figure><img src="image_5.png" /><figcaption>**Figure 2.** The tumoral GPC framework. Cancer cell proliferation and mutation in a tumor can produce different possibilities of tumoral GPCs, depending on the relative fitness of cancer cells at a given time (different colored cells represent distinct evolutionary lineages). Depending on the resulting tumoral GPC, the tumor, viewed as the habitat in which malignant cells live and evolve, possesses specific group properties (e.g., quality of the vascular network, level of immunogenicity, etc.). These properties can, in return, affect (positively or negatively) cell fitness, and hence tumor growth. In the absence of selection at the group level, or of an encoded tumorigenesis program, it is potentially frequent that conditions that increase cell-level fitness of one clonal lineage can result in a non-optimal or even detrimental tumoral GPCs, which can slow down or stop tumor growth, and/ or even induce its size reduction. Since tumors of different sizes have different requirements and interactions with their changing microenvironment, tumoral GPC varies with the tumor stages and the microenvironment; that is, there is no single optimal tumoral GPC that is maintained throughout cancer progression. Only tumors that achieve a successful/adequate tumoral GPC at each step of tumorigenesis will evolve into metastatic tumors. The tumors that fail to generate an adequate tumoral GPC at a given step do not necessarily disappear, they just do not continue to expand. Those that produce an inadequate tumoral GPC, for instance leading to higher immunogenicity, may become reduced in size and even disappear. This hypothesis can explain why we can develop many neoplasms in the body, but the majority of them never grow until the metastatic stage or even regress (Folkman and Kalluri, 2004). Circulating tumor cells (CTCs), especially clusters that can be either homogeneous (organ 3) or heterogeneous (organ 2), can disseminate and initiate metastasis where a novel process of diversification is required so as to harbor the right GPC in a given organ and develop into advanced metastasis.</figcaption></figure>

Thomas et al., 2018). Despite these differences for the majority of other cancers, there is still preda-
tion of malignant cells by the immune system and competition between cells and/or between clus-
ters of cells within tumors. The extrapolation of Farine et al.'s definition of GPC to cancer cannot
totally rely on the same descriptors because groups are composed here of cells, not animals, and
thus cells' specificities need to be considered. For example, phenotypes in animal models that
depend on a high level of mobility and rapid communication within the group might not translate
well because cancer cells move slowly and communicate slowly by diffusible factors. Similarly, the
high mutagenic rate of cancer cells is likely to impact group behavior, while the slow rate of mutation
seen in whole organisms is unlikely to significantly alter group functioning. The extent to which
mechanisms driving evolution or heterogeneity in cancers might have potentially useful analogs in
other group situations still needs to be explored. For instance, promising aspects could include (1)
the importance of epigenetic modulation of gene expression in cancer cells vs the Baldwin effect
(process by which plasticity facilitates evolution) in whole organisms, (2) how might the asexual
nature of cancer cell reproduction influence individual-group evolutionary relationships relative to sexually reproductive individuals in other ecosystems, and (3) what phenotypes in cancer need to be quantified to be able to produce group selection models that are useful for describing progression or influencing clinical practice. Thus, more investigations are required at the moment to fully appreciate what is unique about cancer that will define its individual-group relationships differently from other ecosystems (see also Box 1).

## Group phenotypic composition in the context of cancer

It is traditionally assumed that individual cells are Darwinian units of selection in cancers (Greaves, 2013), as they possess the required prerequisites that include (1) phenotypic variation, (2) differential fitness co-variant with a phenotypic trait, and (3) heritability (Lewontin, 1970). However, some aspects of cancer biology, such as promotion of angiogenesis, evasion of immune predation, or environmental engineering involving remodeling of extracellular matrix require group action. Yet, the evolutionary dynamics that govern these groups (including the level that selection acts on) have been poorly explored (Greaves, 2013; Bertolaso and Dieli, 2017; Lean and Plutynski, 2016). Furthermore, the ecological factors and evolutionary forces shaping these interactions as well as the outcomes of these processes are likely to change during cancer progression. Nevertheless, (1) how the phenotypic composition of a group (i.e., its average phenotype or phenotypic variance) can affect ecological and social processes (e.g., niche construction, competition, cooperation) and (2) how multi-level selection can drive phenotypic covariance among interacting cancer cells are not understood. Following Farine et al., 2015 framework, we propose that addressing the role of GPC as both an agent of selection shaping individual fitness and an emergent property of the individual phenotypes can help understand cancer progression and predict ways that can alter GPC to affect both individual- and group-level outcomes.

How do we define GPC in cancers? We envision a group to (1) be limited to cells that share a set of driver (epi)mutations (but which could subsequently diverge genetically and epigenetically and produce cooperating subgroups each composed of cells sharing at least a common phenotype, see also Box 1), (2) be spatially restricted, given that there are physical limitations to cooperation and competition, and (3) have interactions that are highly dependent on microenvironmental context (considering both other competing normal and malignant cells). Note that we are not including normal cells as part of the malignant group, although its GPC and collective properties will be highly influenced by the interactions with normal cells within the tumor microenvironment.

Soon after initiation, a tumor constitutes a small group, but as it grows, new subgroups can be distinguished, based on new spatial structures and interactions, and its phenotypic composition is likely to change. Whereas specific characteristics, sufficient to define individual subgroups within the GPC still need to be rigorously defined, we envision useful GPC descriptors to include average proliferation rate, degree of genetic or phenotypic heterogeneity, proportion of cells with distinct differentiated states (such as mesenchymal and epithelial for carcinomas), type (luminal and basal, for breast tumors), amount and composition of extracellular matrix, type and numbers of non-malignant cells, proportion of drug-resistant cells, etc.

GPC can be dynamic both in terms of number, types, and proportion of subgroups, and codependency levels between subgroups. For instance, the emergence of macroscopic tumors is dependent on the presence of a threshold number of tumor cells that produce angiogenic factors (Folkman and Klagsbrun, 1987). These factors represent ‘public goods’ that can increase the fitness of all cells, which will ultimately alter the tumor composition. Another example of a public good is the enzyme indoleamine-pyrrole 2,3-dioxygenase, which catalyzes the degradation of tryptophan in the cell’s environment leading to the inhibition of cytotoxic immune responses thereby improving survival of both producing and non-producing cells (Munn and Mellor, 2016). In groups that expand due to sufficiently high proportion of producers, cells that do not produce angiogenic factors (i.e., ‘cheaters’ that rip the benefits of the common goods without incurring the potential cost associated with their production) can potentially outcompete the producers leading to a collapse of the vascular network (Epstein et al., 2017; Gravenmier et al., 2018).

Tumor growth and progression are associated with changes in microenvironmental contexts in space and time (such as acidification, inflammatory response, hypoxia etc.). In turn, these contextual changes would be expected to impact fitness of subgroups of tumor cells. Thus, a GPC that is highly
## Box 2. Examples of GPC in cancer: CTCs.

CTC clusters are defined as groups of two or more CTCs that have stable cell-cell linkages and travel together through the bloodstream (Castro-Giner and Aceto, 2020). Aceto and others reported that CTC clusters have a higher potential for metastasis compared to individual CTCs, suggesting that group dynamics is an important component of metastatic success, at least for a cell's capacity to initiate metastasis (Castro-Giner and Aceto, 2020). Even though characterization of CTC clusters is still in its infancy, different studies have suggested that these groups of cells collectively comprise key and complementary biological features that provide the clustering cells with novel characteristics compared to single CTCs.

For instance, the presence of clusters as well as their size are correlated with a poorer clinical outcome compared to the presence of single CTCs in cancer patients (Aceto et al., 2014). Interestingly, CTC clusters in breast cancer initially detach as multi-cellular aggregates from the primary tumor, precluding the occurrence of intravascular CTC aggregation (Aceto et al., 2014; Cheung et al., 2016).

A recent study might explain the advantage of clustered CTCs. In response to detachment from the extracellular matrix, reactive oxygen species are produced. The cell clusters could be protected from this phenomenon by the production of Hif1α and, subsequently, by a metabolic switch to glycolysis as well as increased survival and metastatic capacity (Labuschagne et al., 2019). This in vitro study remains to be validated in human CTC clusters under physiological conditions. Additional data demonstrated that CTC clusters show different gene expression profiles and dissemination modes compared to single CTCs. Indeed, the epithelial features (cell-cell adhesion via plakoglobin, keratin 14, and E-cadherin) observed in CTC clusters may be important in their formation and maintenance in the circulatory system and may also promote their collective migration and survival beyond maintaining contact between cells (Padmanaban et al., 2019).

Epigenetics also seems to play a key role since DNA methylation is clearly different between single CTCs and CTC clusters. CTC aggregation results in the hypomethylation of binding sites for stemness and proliferation regulators (e.g., OCT4, NANOG, SOX2, and SIN3A). Interestingly, this methylation profile is reversible by dissociating CTC clusters into single cells (Gkountela et al., 2019).

CTC cluster dissemination might also have an impact on a tumor's developmental dynamics in the metastatic setting. In fact, heterogeneous clusters might seed polyclonal metastases, suggesting that this modality of cancer spread may increase the likelihood that a tumor will colonize distant sites successfully and might eventually show increased resistance to anti-cancer therapies (Castro-Giner and Aceto, 2020).

CTCs can form heterotypic interactions with non-tumor cells such as immune cells and cancer-associated fibroblasts (Giuliano et al., 2018; Herath et al., 2020). Szczerba et al., 2019 reported that neutrophils are the most abundant type of immune cell accompanying CTCs; such heterotypic clusters enhance the metastatic potential of CTCs by boosting their proliferative abilities. Alternatively, CTCs can circulate by interacting with fibroblasts, favoring metastasis formation by 'bringing their own soil' (Duda et al., 2010). Moreover, Heeke et al., 2019 mentioned that the CTC clusters have a clear advantage because they 'never travel alone' and they have close crosstalk with the circulating microenvironment. Indeed, cancer-associated fibroblasts and CTCs are considered as 'dangerous liaisons' (Hurtado et al., 2020), since cancer-associated fibroblasts (CAFs) secrete different proteins (cytokines, chemokines, growth factors) that promote invasion and angiogenesis, induce migration of tumor cells, and even confer evasion of the immune system. In conclusion, different types of clusters exist; the very rare and most aggressive ones are the heterotypic CTC clusters compared to the homotypic clusters and single CTCs (Castro-Giner and Aceto, 2020; Duda et al., 2010). The GPC framework could eventually help decipher which CTC clusters are most likely to metastasize (and where). This would provide predictive and/or therapeutic opportunities regarding metastatic success.

'optimal' (reflected in proliferation and death of subgroup members, as well as in expansion of the group) during the early stages of tumorigenesis (e.g., composed of highly proliferative cells) may not be optimal in a large tumor (which requires vasculature) (Losi et al., 2005; Lee et al., 2015; Nguyen et al., 2016; Saito et al., 2018; Figure 2; steps 2 and 3). Similarly, the transition to a metastatic tumor relies on a different 'optimal' GPC (e.g., including a higher proportion of 'aggressive' cells) that reflects changes in the tumor microenvironment (Figure 2; step 4). This change in GPC might involve epigenetic changes as well as phenotypic plastic responses in response to microenvironmental and/or paracrine signals, from both cancer and non-cancer cells. In addition to these spatio-temporal contextual changes within dynamic tumors, the differences in the overall context prior to cancer initiation, such as genetic makeup, sex, hormonal environment, aging, are expected to alter the match between GPC and its impact on the fitness of the group (Laconi et al., 2020; Vibishan and Watve, 2020).
Thus, tumoral GPC, and its relationships with the fitness of individual malignant cells in different ecological contexts, may represent a crucial, previously unexplored aspect that influences evolutionary dynamics within tumor cell populations. In this context, GPC will influence selection pressures experienced by individual tumor cells, which will influence individual cell-level outcomes. As the relative proportion and type of interactions among subgroups of tumor cells vary in response to changing microenvironments, tumoral GPCs will also change resulting in either new optimal tumoral GPCs or, potentially, catastrophic imbalances producing tumor extinction (Figure 2; in this case, tumoral GPC can be seen as the target of selection). As discussed further below, it might be possible to affect GPC to achieve desirable therapeutic outcomes, including tumor stasis or extinction.

The GPC concept might also be applicable to circulating tumor cells (CTCs) (see also Box 2). Individual metastatic tumors can be viewed as reproducing groups (Lean and Plutynski, 2016), where CTCs are seeds of successful groups that travel through blood and lymphatic vasculature to reach distant fertile 'soils' (tissues, organs) (Pantel and Alix-Panabières, 2019). In many cases, successful metastatic colonization is achieved by clusters of CTCs rather than by individual tumor cells (Heeke et al., 2019). While these clusters represent a minority of the CTCs found in the bloodstream, they have been shown to have a metastatic potential 23–50 times greater than that of their single-cell counterparts (Aceto et al., 2014). Moreover, the presence and size of CTC clusters are associated with worse clinical outcome when compared to single CTCs in multiple cancer types (Murlidhar et al., 2017).

CTC clusters can consist of similar or different tumor cells (e.g., expressing epithelial or mesenchymal markers). Clusters can also include non-tumor cells (e.g., blood cells, immune cells, and/or stromal cells) that can affect the fitness of cancer cells (Szczerba et al., 2019). According to our terminology, the single CTCs and CTC clusters composed of homogeneous cells would correspond to the simplest expression of a subgroup that, if successful, would proliferate and diversify to create larger groups with new specific GPCs. On the other hand, CTC clusters consisting of different cells can be viewed as small groups with a GPC that can influence their success during the dispersal, dissemination, and colonization steps. Importantly, the GPC of metastasis resulting from CTC clusters can also change over time, and these changes might be essential for the success of metastatic colonization and progression. For instance, dynamic changes in the epithelial and mesenchymal composition of CTC clusters in breast cancer have been associated with disease progression (Yu et al., 2013). Nevertheless, it is less clear whether this heterogeneity reflects potential cooperative interactions among cells within a CTC cluster and whether this change in GPC impacts metastatic efficiency. The extent to which CTC cluster composition should recapitulate the parental GPC (perhaps with modifications) to be successful remains to be determined. We cannot exclude that specific GPCs, different from the parental one are required to target given organs, with convergence between metastatic tumors independent of the localization of their primary tumors (Cunningham et al., 2015). Box 2 discusses CTCs as an example of the GPC dynamics in cancer.

## GPC in transmissible cancers

The GPC concept can also be extended to transmissible cancers (see also Box 3 that uses transmissible cancers as another example of the GPC in cancer). Previous work has argued that the emergence of transmissible cancers required a 'perfect storm' with the convergence of multiple host (micro- and macroenvironmental factors) and tumor cell traits (Ujvari et al., 2016a). Here, we hypothesize that the tumoral GPC is also an important component of the 'perfect storm' necessary for the evolution of a transmissible cancer. The majority of cancer cells that emerge in multicellular organisms die with their host, but there are at least nine well-described, independent exceptions of transmissible cancers in the wild: one in dogs (called canine transmissible venereal tumor, CTVT), two in Tasmanian devils (called devil facial tumor diseases, DFTD and DFT2), and six in six bivalve species (called bivalve transmissible neoplasias [BTN], see Ujvari et al., 2017; Yonemitsu et al., 2019). In these diseases, malignant lineages behave like infectious agents, with cancer cells spreading horizontally as allografts and/or xenografts from one individual to conspecifics, or to individuals of (closely) related species. These cancer cells are currently viewed as new parasitic species, with their own evolution and dynamics (Dingli and Nowak, 2006).

The first steps during the emergence of a transmissible cancer are likely to be similar to non-transmissible ones (i.e., as described in the sections above), with similar constraints on the GPC composition required for a tumor to grow and metastasize. However, to become transmissible, a cancer
must also gain the ability to shed a very large number of cells that can both survive in the environment outside the host's body (Sunila and Farley, 1989), but also be able to infect new hosts (Ujvari et al., 2016b). Since we have thus far only observed transmissible cancers that successfully spread to a large number of individuals, there is at the moment no evidence of what would be a favorable tumoral GPC for an emerging cancer to become transmissible (at least for the infectious stage). Both CTVT and DFTD manifest as undifferentiated cancers, showing strikingly uniform cell morphology. However, the tumors are not completely homogeneous, and variability in antigen expression in and between tumors has been observed (O'Neill, 2011). Since they are found in both host species, those tumoral GPC characteristics are likely to be important ones for a cancer to be transmissible.

Once inside the body of a new host, most transmissible cancer cells fail to establish a new tumor (10⁵–10⁷ cells were required in transmission experiments [House, 1997a; Kreiss et al., 2011]). This high failure rate can be partly explained by the fact that a large proportion of cells are eliminated by the host's immune system or fail to find a suitable environment to establish. In addition, it is also possible that most of the newly initiated transmissible cancer tumors regress because of a tumoral GPC mismatch (in a similar way many of the de novo tumors can fail) (see further details in Box 3).

It is also possible that different types of tumoral GPC are required for every step of the intra- and inter-individual transmission process. As transmissible cancers show intra- and inter-individual metastasis, every successful transmission requires establishment in a novel microenvironmental niche and the avoidance of different levels of immune recognitions. For example, a particular GPC may be beneficial when spreading within a host via initiating immune tolerance to the tumoral GPC, but the same GPC may also present as additional immune target, and thus not be 'optimum', when transmitting across hosts. For example, the newly emerged devil facial tumor disease, DFT2, avoids immune recognition by various mechanisms: some tumors express MHC class I molecules that are shared with host devils, some express non-classical MHC molecules, while in some MHC class I, expression has been completely lost (Caldwell et al., 2018). Further work would be needed to determine whether this also occurs within tumors and to establish the extent of intratumoral variability in these features.

A major difference between non-communicable cancers and transmissible cancers is that classical cancer cells have to 'reinvent the wheel' at every step – from initiation and progression to metastasis, but ultimately, they all die with the death of the host. In contrast, in transmissible cancers, clonal cell lineages are passed from host to host, and thus evolve as novel parasitic life forms, as seen in CTVT and DFTD, where these cancer lineages have continued to evolve since their initial emergence (Pearse et al., 2012; Dujon et al., 2020a; Kwon et al., 2020). Thus, investigating whether the tumoral GPC has been adjusted by selection in a way that optimizes the parasitic lifestyle will also provide insight into the evolution of parasitism. Focus on BTN, where host switch has been observed, would be particularly interesting as different tumoral RPCs might have been required at every stage of transmission (within host, between hosts, and between species). There are, however, no data available concerning the GPC of transmissible cancers, especially at infection time. Therefore, all the hypotheses discussed above remain to be tested.

## Therapeutic implications

Currently, the major focus of therapeutic development is on directly targeting viability or proliferation of tumor cells. Consideration of the GPC perspective could suggest new therapeutic approaches to achieve desirable clinical outcomes. While, at this point, prerequisite knowledge is still lacking, below we provide several suggestions on the new therapeutic angles inspired from the GPC framework.

First, we should attempt to identify keystone tumor cell subgroups/phenotypes – those critical for the maintenance of tumors and transitions (i.e., those that contribute to the 'optimum' GPC, Figure 2) at key steps of their progression (invasion, metastasis, acquisition of therapy resistance), and either target them directly or focus on the mechanisms underlying cross talk among phenotypes. This strategy, compared to those that target the entire heterogeneous tumor, will be more efficient because it would target smaller and more homogenous and less diverse subgroups. This should decrease the risk of emergence of resistant clones. The emergence of therapy resistance is also accompanied by changes in GPC; understanding the specific RPCs that promote resistance should
**Box 3. Examples of GPC in cancer: transmissible cancers.**

Compared to human cancer research, the study of transmissible cancers is in its relative infancy, with fewer and smaller teams working on the topic (Dujon et al., 2020b). Finding empirical evidence for GPC to exist in these (so far) rare cancers (Dujon et al., 2021) is less straightforward. Nevertheless, based on emerging transcriptomic, genomic, and experimental data, predictions could be made for their existence. For example, experimental transmission of all three types of transmissible cancers indicates the potential relevance of GPC. Immunocompromised xenograft murine models have been developed for both DFTD (Kreiss et al., 2011) and CTVT (Holmes, 1981; Harmelin et al., 2001). In the DFTD model, NOD/SCID mice inoculated with either 10⁵ or 10⁶ viable DFTD tumor cells developed tumors after 7 and 4 weeks, respectively, and X-DFTD tumors reinoculated into NOD/SCID mice formed visible tumors after 5 weeks (Kreiss et al., 2011). In CTVT murine models, 10⁶ cells were sufficient to cause tumor development in NOD/SCID mice (Harmelin et al., 2001). A study from the early 1980s that followed the shading of radioisotope-labeled CTVT cells in urine found that after the initial rapid phase of cell death, approximately 10% of the injected label was detectable, and subcutaneous CTVT nodules developed by day 7 (Holmes, 1981). Similarly, in the transmissible cancer system of bivalves, experimental transmission of 10⁶ BTN cells induced cancer in 33% of new hosts, and transfection of 10⁷ BTN cells resulted in a 50% infection rate (House, 1997b). These studies indicate that although most likely a large proportion of cancer cells die during transmission, a small population of viable cells could be sufficient to establish a tumor mass. While there is currently no evidence available, it could be proposed that—similar to CTCs—the success and survival of transmissible cancer cells may be higher when ‘traveling’ in groups or having a particular GPC at a given time point of transmission (e.g., extravasation, transmission, intravasation). Collective transmission of a group of cells with diverse and/or particular phenotypes may modulate the emergent properties of transmissible cancer cells that facilitate the breakdown of tissues as well as individual and environmental barriers, eventually leading to persistent group-level differences in fitness. This could potentially be most relevant for transmissible cancers in bivalves, where transmission occurs by release and uptake of hemocytes into and from the aquatic environment (Metzger et al., 2016). Disseminated neoplasia cells have been shown to have relatively broad tolerance to environmental conditions (10–15% mortality below 20°C that rapidly increases to 80% at 30°C, salinity tolerance between 0.5 psu and <35 psu, and pH preference between pH 4 and 9.3 [Bramwell et al., 2021; Sunila, 1991]); this plasticity might be the result of and facilitated by the production of different GPC released by BTN cells at a given point in time.

Although empirical evidence is currently lacking to support the hypotheses outlined above, future transfection experiments could prove or disprove these theories. For example, mussels have a natural tendency to open their valves and to release and uptake hemocytes (Caza et al., 2020). Placing mussels under different environmental conditions and exposing them to stained BTN cells (e.g., CFSE, a cytosol dye commonly used for tracking cells in vivo by flow cytometry), followed by transcriptome and epitome analyses of both established BTN and host cells may reveal the phenotypic composition of BTN cells that were successfully taken up by and survived in the new host (Caza et al., 2020).

Another direction for investigating the importance of GPC in transmissible cancers could be the Tasmanian devil and its two independently emerged contagious cell lineages. The first lineage, DFT1, was initially observed in 1996 in the northeastern corner of Tasmania; it has quickly spread across Tasmania, including the southern regions. The second independently risen lineage, DFT2, was discovered in 2014 (McCallum et al., 2009; Pye et al., 2016) in the southern corner of Tasmania. The distributions of the two lineages now overlap in the d'Entrecasteaux Peninsula (James et al., 2019). DFT1 tumors mostly occur on the face, with non-facial tumors more commonly found in DFT2 (130), indicating potential niche partitioning of the host by the two cancer lineages. Although co-infection by DFT1 and DFT2 occurs, it has been rare thus far, and DFT2 has seemingly started replacing DFT1 on the local scale (James et al., 2019). The current rapid expansion of the DFT2 lineage may be occurring because its high levels of GPC facilitate evasion of immune recognition. As described above, the two lineages employ different immune system evasion strategies: DFT1 cells lack MHC expression on their cell surface while DFT2 cells express classical and non-classical MHC alleles (Caldwell et al., 2018). The expression of MHC molecules in DFT2 is not uniform – it is highly heterogenic in vitro and in vivo. For example, some DFT2 and DFT1 cells share MHC class I alleles, some DFT2 cells also express high levels of an MHC class I allele common to its hosts, and some DFT2 cells do not express MHC at all (Caldwell et al., 2018). Although only conjectural, it is possible that during transmission of a cluster of tumor cells, the heterogeneity in immune recognition evasion may be beneficial for these contagious cancers since it may affect phenotype-dependent predation risk via confusion and oddity effects (i.e., a complete lack of MHC expression could trigger natural killer cell responses, but these are being silenced due to certain cells expressing alleles similar to the host).

While conducting co-infection experiments on Tasmanian devils would be nearly impossible, in vitro cell competition assays could potentially be employed to investigate whether the GPC of the two lineages and co-infection by 'mixed-species' (the two cell lines representing species in this case) would benefit one or both lineages via resource allocation and niche partitioning.
help develop strategies that prevent the acquisition of resistant GPCs. Below we expand on these
two complementary approaches and provide specific experimental directions.

In order to gain the ability to consider GPC in therapeutic decision making, we first need to develop approaches to identify keystone subgroups/phenotypes and specific combinations of interacting phenotypes that result in optimum GPCs. Examples of such keystone subgroups and interactions have been identified in multiple experimental contexts (Tabassum and Polyak, 2015). For example, in experimental models of breast cancer heterogeneity, interleukin 11 expressing subpopulations can be responsible for driving tumor growth and metastatic dissemination via mechanisms involving microenvironmental changes (Marusyk et al., 2014). Similarly, aggressive subtypes in glioma and breast cancer can affect, via paracrine signaling, the phenotype and metastatic potential of less aggressive subtypes (Neelakantan et al., 2017; Espinoza-Sánchez et al., 2017). However, transitioning from proof of principle studies toward characterizing GPC in primary tumors and identifying key subpopulations would necessitate stepping out of the dominant focus on molecular mechanisms, following 'Gene X is a critical regulator of hallmark H' algorithm, and developing system-level understanding of tumor cell subpopulations and their interactions. Recent advances in single-cell analyses, such as single-cell RNA sequencing and high multiplexing histological analyses, are already providing us with detailed characterization of phenotypes of individual tumor cells and/or subgroups that constitute the GPC. However, in order to transition from simply describing the GPC into identification of keystone subgroups and key biological interactions, we would need to develop network analysis pipelines that explicitly focus on the analyses of expression of receptors and corresponding ligands (including both secreted ligands and extracellular matrix (ECM)), and identify or engineer experimental systems that recapitulate individual interaction nodes, enabling the interrogation of the consequences of their disruption. Moreover, rather than focusing on a single interaction at a time, we will need to develop strategies for co-disrupting redundant interactions responsible for robustness and antifragility of tumor ecosystems.

Second, in many cases of systemic treatment of multisite metastatic disease, individual tumors display marked variability in responses to therapy, ranging from complete elimination to lack of detectable therapeutic responses (e.g., *Patel et al., 2021*). Comparative analyses of the GPC of 'successful' versus 'failed' tumors represent another potential approach to interrogate the role of GPC in the development of resistance to therapies. While, traditionally, the explanation sought is in the identification of cell-intrinsic 'drivers' of resistance, comprehensive single-cell characterization of pretreatment biopsies could discriminate between GPCs associated with tumors that succumb to therapy and those that develop resistance. Even if metastases are rarely resected, which hampers their characterization at the cellular/molecular level, such analyses might sometimes be feasible in clinical presentations of multiple skin metastases, such as in *Wagle et al., 2011*, or in animal models of heterogeneous multi-metastatic disease. A limitation to this suggestion is, however, the fact that certain cells and/or subgroups important for the GPC could be missed for technical reasons.

Third, the phenotypic composition of successful/resistant tumors is also influenced by the ability of tumor cells to switch phenotypes. Therefore, explicitly targeting phenotypic plasticity (independently of cytotoxic/cytostatic effects of treatment) might be able to prevent plastic adaptive changes in GPC and improve therapeutic outcomes. Phenotypic plasticity should not be considered strictly as a means to increase individual/cell-level fitness. In our framework, phenotypic plasticity is also a way for cancer cells to both rapidly change the tumoral GPC to maintain the combination of phenotypes required for progression or overcome the therapeutic suppression of a dominant aggressive phenotype. This means that optimal or suboptimal tumoral GPCs could be restored after such interventions simply as a consequence of these non-genetic events. To circumvent these phenomena, targeting phenotypic plasticity should be relevant, either alone or in combination with targeted therapies (Doherty et al., 2016). This is essential when considering combinations with other therapies targeting one or several particular phenotype(s)/subgroup(s) because alteration of the possibility to change the tumoral GPC should avoid reconstituting the missing subpopulation.

Phenotypic plasticity is extremely relevant especially during the epithelial-to-mesenchymal transition (EMT). EMT generates mesenchymal-like cells and a variety of intermediate cell states between the epithelial and the mesenchymal state and changes the GPC toward tumor progression and metastasis (Brooks et al., 2015; Dongre and Weinberg, 2019; Rios et al., 2019). Thus, targeting molecular machinery involved in phenotypic rewiring (such as inhibitors of HDAC, KDM, BETs) should limit the evolution of the GPC both during tumor growth and after a therapeutic intervention. Many
other molecular mechanisms underlie cell plasticity, for instance the Notch and Wnt development pathways, or those involving transcription factors from the Snail, Zeb, and Twist families, whose role in the generation of plasticity and cell-to-cell heterogeneity in cancer partially overlaps with their role during development and wound healing (Gupta et al., 2019). Deregulation of these transcriptional regulators may allow cancer cells to access resistant phenotypes by rewiring of gene regulatory networks (Marusyk et al., 2020).

More globally, insights can be gained on the role of phenotypic plasticity in affecting tumor GPC using dynamical systems theory (Huang et al., 2009; Jia et al., 2017) where reshaping of the epigenetic landscape in cancer cells allows acquisition of novel abnormal cell states corresponding to 'cancer attractors' in the gene regulatory network (Huang et al., 2009). The enhanced cellular stochasticity in cancer cells would facilitate these phenotypic conversions toward novel states (Huang et al., 2009; Li et al., 2016). Stochasticity of gene expression appears to be increased in cancer cells (Jenkinson et al., 2017), and diverging chromatin states and transcriptional heterogeneity produce corrupted coordination of epigenetic modifications (Pastore et al., 2019). Also, the stochastic nature of gene expression fosters a transient resistant phenotype by chance, which can then be stabilized by epigenetic mechanisms (Shaffer et al., 2017). As such, modulating stochasticity, especially through tuning of gene expression variability, might represent a good strategy to control tumoral GPC evolution (reduced phenotypic heterogeneity will limit cooperative interactions). The global epigenetic changes generally observed in cancer cell epigenomes could find their origin in 'tumor reprogramming' due to genetic alterations initiating cancer that would contribute to cancer growth less by inducing cell proliferation than by causing 'developmental reprogramming' of the epigenome and allowing cells to aberrantly and pathologically differentiate (Vicente-Dueñas et al., 2015; Xiong et al., 2019). Another possibility is that this global destabilization of gene expression could originate from tissue disruption-induced stochasticity (Capp, 2005; Capp, 2017; Capp, 2019; Capp and Bataille, 2020), suggesting that restoring cellular interactions present in the initial tissue for instance with molecules that would mimic healthy cellular interactions would reduce cell stochasticity and restabilize cellular phenotypes (Capp and Bataille, 2020; Capp, 2012; Capp and Bataille, 2018). In any case, disrupting the tumoral GPC will require innovative strategies based on the control of phenotypic plasticity.

Fourth, while currently relevant knowledge is almost entirely lacking, it will be important to understand whether and how the tumoral GPC is reestablished in a secondary tumor, and the role the tumoral GPC plays in metastatic colonization and outgrowth. We can speculate that since metastases can be seeded by individual cells or small groups of cells (Cheung and Ewald, 2016), the odds that the 'right combination of cells' is initially present will be very small. Consequently, if the optimum GPC needs to be reestablished at a metastatic site, the failure to do so could contribute (along with the many other hurdles faced by prospective metastases) to the extremely low colonization rate of circulating tumor cells. In addition, the establishment of a new tumor in a completely different tissue microenvironment (e.g., breast cancer metastasis to the bone) should require a GPC that is quite distinct from that in the native site. The growing metastatic tumor would need to adapt to a very different environment, which could be facilitated by phenotypes that secrete cytokines or chemokines that both recruit regional cells into the tumor's service and impede antagonistic cells like cytotoxic T cells, by other phenotypes that promote degradation of the extracellular matrix and/or acidification, and so forth (Quail and Joyce, 2013). The assembly of such a GPC may (fortunately) represent a substantial hurdle to metastatic success. On the other hand, the seeding metastatic cells may already possess greater phenotypic plasticity and an enhanced capacity for generating genetic and epigenetic heterogeneities.

Fifth, in addition to targeting keystone subgroups, tackling specific interactions (such as clonal cooperation) may be of equal or higher therapeutic benefit due to altering tumoral GPC quality. For instance, disrupting clonal cooperation by targeting key factors was already suggested (Zhou et al., 2017), as well as their shared external products (i.e., 'public goods') (Pepper, 2012). A list of potential targets includes members of signaling pathways, such as the proto-oncogenes Wnt1 or Hedgehog, and microenvironment-related factors, such as matrix metalloproteinases (Zhou et al., 2017), although the exact useful target is likely to be specific to each individual tumor type. Inhibition of tumor progression can be achieved by disruption of these clonal cooperative interactions, which impair the tumor ecosystem and the tumoral GPC. While only a few of the many pathways/molecules
involved in clonal cooperativity are characterized, it should be expected that the coming years will
bring new opportunities to inhibit this aspect of tumor biology.

Concluding remarks

Tumors can be seen as examples of social heterogeneous groupings. A proper mechanistic under-
standing of the impact of ITH and tumoral GPC on cancer progression requires bridging oncological
sciences with evolutionary ecology. What are the potential advantages of GPC consideration over
traditional, reductionistic approach of viewing tumors as collection of specific genetic clones or phe-
notypic subpopulations? A first obvious benefit of the GPC concept resides in the fact that it permits
assessment of some of the functional aspects of cancer dynamics, which are impossible to under-
stand by simply examining cellular collections only. In the same way that we cannot fully understand
the biology of any living organisms when examined in isolation from its abiotic and biotic environ-
ment (including its conspecific peers), the biology of malignant cells must consider the social context
(i.e., the GPC) to which it is exposed. The GPC concept should also permit us to understand pheno-
typic convergences, when cells differ at the genetic or epigenetic level but accomplish more or less
the same tasks (see Box 1).

On the other hand, the utility of considering tumor GPC is currently limited by the lack of criteria to define relevant groups of tumor cells. Strong interdisciplinary research that integrates observational, experimental, and theoretical studies will be crucial to develop methodology to define and categorize GPCs, identify functionally relevant subpopulations within heterogeneous tumors, identify critical interactions between these subpopulations, and identify the parameter values of GPC descriptors that are associated with an 'optimum' GPC at each cancer progression step. Such studies will ultimately allow us to understand the impact of GPCs (and the interactions between metastable components of the GPCs and their microenvironments) on the collective behavior of tumors and their evolutionary and clinical trajectories.

In this paper, we did not consider the potential existence of metphenotypes of structured
groups of cells, that is functional clusters, within tumors. Such a situation would yield to different
and more complex GPC-related dynamics on tumor progression. While whole tumors are (most
often) not in competition within tissues, this can conversely be the case for functional clusters within
tumors, suggesting that the GPC of tumors could display a 'nested structure'. That is, the GPCs of
functional clusters within the tumor are also influenced by the competition between groups, but
these groups also contribute together to a whole tumoral GPC that is not in a competitive context
with other tumors. In addition, if a group within the tumor has a better GPC than others, it can
increase in prevalence, with tumor growth then seeding the formation of new GPCs. More research
is needed at the moment to explore the extent to which tumor dynamics are influenced by the whole
tumor GPC, by GPCs of functional groups that are within tumors, and/or by a mix of both. Also, we
need to determine whether there are temporal patterns in these processes, exploring for instance
the importance of these different dynamics throughout tumorigenesis. Still, the formation of 'the
right GPC', appearing in the right condition at the right time, is heavily influenced by chance, sto-
chasticity and spatial context, and thus remains challenging to investigate and act on. For instance,
cells that would form a viable GPC for host-to-host transmission may exist but are not co-localizing
or cannot be transmitted as a unit during shedding/infection. Similarly, some phenotypes may bene-
fit from an advantage simply because they are close to angiogenic cells, while other GPC-formation
candidates in other parts of the tumor will die due to nutrient scarcity. Other information is also
needed to determine whether the propensity for developing the 'right' GPC (see Figure 2) during
tumorigenesis may rely on the individual's genetic predisposition. It is noteworthy that even for
genetic predisposition syndromes that confer the highest probability of developing cancer such as
those resulting from inherited dysfunction of BRCA1/2 (predisposing to breast and ovarian cancers)
and mismatch repair genes (e.g., Lynch syndrome), a significant proportion of patients (30-40%) do
not develop cancer by the age of 70 (Antoniou and Easton, 2006). The GPC might well be an
important modulating factor of the risk over time.

Overall, examining tumor GPC and addressing the effects of GPC on individual tumor cell fitness (selective consequences) and the responses of individual cells to the selective consequences of GPC (evolutionary consequences) should improve our understanding of tumor biology and cancer progression. At the moment, it remains mostly an abstract concept because we need to develop novel
tools to describe, understand, and monitor the tumoral GPC functional dynamics and its resilience to local disturbances. Also, criteria used to describe groups remain by definition subjective, limiting the possibility at the moment to provide a unique definition of GPC in cancers (see Box 1 for a tentative example). This knowledge will also need to move beyond this theoretical phase to become useful for diagnostics or therapeutic design, providing opportunities for therapies that will specifically and durably disrupt key components and interactions within tumors.

## Acknowledgements

This work was supported by the MAVA Foundation and by the following grants: ANR TRANSCAN (ANR-18-CE35-0009), ARC Linkage (LP170101105), Deakin SEBE_RGS_2019, a grant from AACR/ Johnson and Johnson grant, and a CNRS International Associated Laboratory Grant. The funders had no role in study design, data collection and analysis, decision to publish, or manuscript preparation. We would like to thank Guy Lyons, Pierre Martinez, and one anonymous referee for their constructive comments on an earlier version of this manuscript.

## Additional information

### Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>ANR</td><td>ANR-18-CE35-0009</td><td>Frédéric Thomas</td></tr><tr><td>MAVA Foundation</td><td>2020-2023</td><td>Frédéric Thomas</td></tr><tr><td>Australian Research Council</td><td>LP170101105</td><td>Rodrigo Hamede</td></tr><tr><td>Deakin University</td><td>SEBE_RGS_2019</td><td>Beata Ujvari</td></tr><tr><td>AACR/Johnson&Johnson</td><td></td><td>James DeGregori</td></tr><tr><td>CNRS International Associated Laboratory</td><td></td><td>Rodrigo Hamede<br>Beata Ujvari<br>Frédéric Thomas</td></tr></tbody></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

### Author contributions

Jean-Pascal Capp, James DeGregori, Aurora M Nedelcu, Beata Ujvari, Andriy Marusyk, Robert Gatenby, Frédéric Thomas, Conceptualization, Supervision, Writing - original draft, Writing - review and editing; Antoine M Dujon, Pascal Pujol, Catherine Alix-Panabières, Rodrigo Hamede, Writing - original draft, Writing - review and editing; Justine Boutry, Writing - original draft; Benjamin Roche, Conceptualization, Writing - original draft, Writing - review and editing

### Author ORCIDs

Jean-Pascal Capp https://orcid.org/0000-0002-6470-079X

Aurora M Nedelcu https://orcid.org/0000-0002-7517-2419

Catherine Alix-Panabières https://orcid.org/0000-0002-6401-2903

Beata Ujvari https://orcid.org/0000-0003-2391-2988

Robert Gatenby https://orcid.org/0000-0002-1621-1510

Frédéric Thomas https://orcid.org/0000-0003-2238-1978

## References

* Aceto N, Bardia A, Miyamoto DT, Donaldson MC, Wittner BS, Spencer JA, Yu M, Pely A, Engstrom A, Zhu H, Brannigan BW, Kapur R, Stott SL, Shioda T, Ramaswamy S, Ting DT, Lin CP, Toner M, Haber DA, Maheswaran S. 2014. Circulating tumor cell clusters are oligoclonal precursors of breast Cancer metastasis. *Cell* **158**:1110-1122. DOI: https://doi.org/10.1016/j.cell.2014.07.013, PMID: 25171411
Antoniou AC, Easton DF. 2006. Models of genetic susceptibility to breast Cancer. *Oncogene* **25**:5898-5905.
DOI: https://doi.org/10.1038/sj.onc.1209879, PMID: 16998504

Axelrod R, Axelrod DE, Pienta KJ. 2006. Evolution of cooperation among tumor cells. *PNAS* **103**:13474-13479.
DOI: https://doi.org/10.1073/pnas.0606053103

Beck B, Blanpain C. 2013. Unravelling Cancer stem cell potential. *Nature Reviews Cancer* **13**:727-738.
DOI: https://doi.org/10.1038/nrc3597, PMID: 24060864

Bertolaso M, Dieli AM. 2017. Cancer and intercellular cooperation. *Royal Society Open Science* **4**:170470.
DOI: https://doi.org/10.1098/rsos.170470, PMID: 29134064

Böttger K, Hatzikirou H, Voss-Böhme A, Cavalcanti-Adam EA, Herrero MA, Deutsch A. 2015. An emerging allee effect is critical for tumor initiation and persistence. *PLOS Computational Biology* **11**:e1004366. DOI: https://doi.org/10.1371/journal.pcbi.1004366, PMID: 26335202

Bramwell G, Schultz AG, Sherman CDH, Giraudeau M, Thomas F, Ujvari B, Dujon AM. 2021. A review of the potential effects of climate change on disseminated neoplasia with an emphasis on efficient detection in marine bivalve populations. *Science of the Total Environment* **775**:145134. DOI: https://doi.org/10.1016/j.scitotenv.2021.145134

Brooks MD, Burness ML, Wicha MS. 2015. Therapeutic implications of cellular heterogeneity and plasticity in breast Cancer. *Cell Stem Cell* **17**:260-271. DOI: https://doi.org/10.1016/j.stem.2015.08.014, PMID: 26340526

Calbo J, van Montfort E, Proost N, van Drunen E, Beverloo HB, Meuwissen R, Berns A. 2011. A functional role for tumor cell heterogeneity in a mouse model of small cell lung Cancer. *Cancer Cell* **19**:244-256. DOI: https://doi.org/10.1016/j.ccr.2010.12.021, PMID: 21316603

Caldwell A, Coleby R, Tovar C, Stammnitz MR, Kwon YM, Owen RS, Tringides M, Murchison EP, Skjødt K, Thomas GJ, Kaufman J, Elliott T, Woods GM, Siddle HV. 2018. The newly-arisen Devil facial tumour disease 2 (DFT2) reveals a mechanism for the emergence of a contagious Cancer. *eLife* **7**:e35314. DOI: https://doi.org/10.7554/eLife.35314, PMID: 30103855

Cancer Genome Atlas Network. 2012. Comprehensive molecular portraits of human breast tumours. *Nature* **490**:61-70. DOI: https://doi.org/10.1038/nature11412, PMID: 23000897

Capp JP. 2005. Stochastic gene expression, disruption of tissue averaging effects and Cancer as a disease of development. *BioEssays* **27**:1277-1285. DOI: https://doi.org/10.1002/bies.20326, PMID: 16299757

Capp JP. 2012. Stochastic gene expression stabilization as a new therapeutic strategy for Cancer. *BioEssays* **34**:170-173. DOI: https://doi.org/10.1002/bies.201100149, PMID: 22231918

Capp JP. 2017. Tissue disruption increases stochastic gene expression thus producing tumors: cancer initiation without driver mutation. *International Journal of Cancer* **140**:2408-2413. DOI: https://doi.org/10.1002/ijc.30596, PMID: 28052327

Capp JP. 2019. Cancer stem cells: from historical roots to a new perspective. *Journal of Oncology* **2019**:1-10.
DOI: https://doi.org/10.1155/2019/5189232, PMID: 31308849

Capp JP, Bataille R. 2018. Multiple myeloma exemplifies a model of Cancer based on tissue disruption as the initiator event. *Frontiers in Oncology* **8**:355. DOI: https://doi.org/10.3389/fonc.2018.00355, PMID: 30250824

Capp JP, Bataille R. 2020. Multiple myeloma as a bone disease? The tissue Disruption-Induced cell stochasticity (TiDiS) Theory. *Cancers* **12**:2158. DOI: https://doi.org/10.3390/cancers12082158

Castro-Giner F, Aceto N. 2020. Tracking Cancer progression: from circulating tumor cells to metastasis. *Genome Medicine* **12**:31. DOI: https://doi.org/10.1186/s13073-020-00728-3, PMID: 32192534

Caza F, Bernet E, Veyrier FJ, Betoulle S, St-Pierre Y. 2020. Hemocytes released in seawater act as trojan horses for spreading of bacterial infections in mussels. *Scientific Reports* **10**:19696. DOI: https://doi.org/10.1038/s41598-020-76677-z, PMID: 33184419

Chapman A, Fernandez del Ama L, Ferguson J, Kamarashev J, Wellbrock C, Hurlstone A. 2014. Heterogeneous tumor subpopulations cooperate to drive invasion. *Cell Reports* **8**:688-695. DOI: https://doi.org/10.1016/j.celrep.2014.06.045, PMID: 25066122

Cheung KJ, Padmanaban V, Silvestri V, Schipper K, Cohen JD, Fairchild AN, Gorin MA, Verdone JE, Pienta KJ, Bader JS, Ewald AJ. 2016. Polyclonal breast Cancer metastases arise from collective dissemination of keratin 14-expressing tumor cell clusters. *PNAS* **113**:E854-E863. DOI: https://doi.org/10.1073/pnas.1508541113, PMID: 26831077

Cheung KJ, Ewald AJ. 2016. A collective route to metastasis: seeding by tumor cell clusters. *Science* **352**:167-169. DOI: https://doi.org/10.1126/science.aaf6546, PMID: 27124449

Cleary AS, Leonard TL, Gestl SA, Gunther EJ. 2014. Tumour cell heterogeneity maintained by cooperating subclones in Wnt-driven mammary cancers. *Nature* **508**:113-117. DOI: https://doi.org/10.1038/nature13187, PMID: 24695311

Cunningham JJ, Brown JS, Vincent TL, Gatenby RA. 2015. Divergent and convergent evolution in metastases suggest treatment strategies based on specific metastatic sites. *Evolution, Medicine, and Public Health* **2015**:76-87. DOI: https://doi.org/10.1093/emph/eov006, PMID: 25794501

Dagogo-Jack I, Shaw AT. 2018. Tumour heterogeneity and resistance to Cancer therapies. *Nature Reviews Clinical Oncology* **15**:81-94. DOI: https://doi.org/10.1038/nrclinonc.2017.166, PMID: 29115304

Dalerba P, Kalisky T, Sahoo D, Rajendran PS, Rothenberg ME, Leyrat AA, Sim S, Okamoto J, Johnston DM, Qian D, Zabala M, Bueno J, Neff NF, Wang J, Shelton AA, Visser B, Hisamori S, Shimono Y, van de Wetering M, Clevers H, et al. 2011. Single-cell dissection of transcriptional heterogeneity in human Colon tumors. *Nature Biotechnology* **29**:1120-1127. DOI: https://doi.org/10.1038/nbt.2038, PMID: 22081019

Dingli D, Nowak MA. 2006. Cancer biology: infectious tumour cells. *Nature* **443**:35-36. DOI: https://doi.org/10.1038/443035a, PMID: 16957717
Doherty M, Smigiel J, Junk D, Jackson M. 2016. Cancer stem cell plasticity drives therapeutic resistance. *Cancers* **8**:8. DOI: https://doi.org/10.3390/cancers8010008

Dongre A, Weinberg RA. 2019. New insights into the mechanisms of epithelial-mesenchymal transition and implications for Cancer. *Nature Reviews Molecular Cell Biology* **20**:69–84. DOI: https://doi.org/10.1038/s41580-018-0080-4, PMID: 30459476

Duda DG, Duyverman AM, Kohno M, Snuderl M, Steller EJ, Fukumura D, Jain RK. 2010. Malignant cells facilitate lung metastasis by bringing their own soil. *PNAS* **107**:21677–21682. DOI: https://doi.org/10.1073/pnas.1016234107, PMID: 21098274

Dujon AM, Gatenby RA, Bramwell G, MacDonald N, Dohrmann E, Raven N, Schultz A, Hamede R, Gérard AL, Giraudeau M, Thomas F, Ujvari B. 2020a. Transmissible cancers in an evolutionary perspective. *iScience* **23**: 101269. DOI: https://doi.org/10.1016/j.isci.2020.101269, PMID: 32592998

Dujon AM, Schofield G, Bramwell G, Raven N, Hamede R, Thomas F, Ujvari B. 2020b. Global meta-analysis of over 50 years of multidisciplinary and international collaborations on transmissible cancers. *Evolutionary Applications* **13**:1745–1755. DOI: https://doi.org/10.1111/eva.12938

Dujon AM, Bramwell G, Roche B, Thomas F, Ujvari B. 2021. Transmissible cancers in mammals and bivalves: how many examples are there?: predictions indicate widespread occurrence. *BioEssays: News and Reviews in Molecular, Cellular and Developmental Biology* **43**:e2000222. DOI: https://doi.org/10.1002/bies.202000222, PMID: 33210313

Epstein T, Gatenby RA, Brown JS. 2017. The warburg effect as an adaptation of Cancer cells to rapid fluctuations in energy demand. *PLOS ONE* **12**:e0185085. DOI: https://doi.org/10.1371/journal.pone.0185085, PMID: 28922380

Espinoza-Sánchez NA, Vadillo E, Balandrán JC, Monroy-García A, Pelayo R, Fuentes-Paná EM. 2017. Evidence of lateral transmission of aggressive features between different types of breast Cancer cells. *International Journal of Oncology* **51**:1482–1496. DOI: https://doi.org/10.3892/ijo.2017.4128, PMID: 29048610

Farine DR, Montiglio PO, Spiegel O. 2015. From individuals to groups and back: the evolutionary implications of group phenotypic composition. *Trends in Ecology & Evolution* **30**:609–621. DOI: https://doi.org/10.1016/j.tree.2015.07.005, PMID: 26411618

Flavahan WA, Gaskell E, Bernstein BE. 2017. Epigenetic plasticity and the hallmarks of Cancer. *Science* **357**: eaal2380. DOI: https://doi.org/10.1126/science.aal2380, PMID: 28729483

Folkman J, Kalluri R. 2004. Cancer without disease. *Nature* **427**:787. DOI: https://doi.org/10.1038/427787a, PMID: 14985739

Folkman J, Klagsbrun M. 1987. Angiogenic factors. *Science* **235**:442–447. DOI: https://doi.org/10.1126/science.2432664, PMID: 2432664

Giuliano M, Shaikh A, Lo HC, Arpino G, De Placido S, Zhang XH, Cristofanilli M, Schiff R, Trivedi MV. 2018. Perspective on circulating tumor cell clusters: why it takes a village to metastasize. *Cancer Research* **78**:845–852. DOI: https://doi.org/10.1158/0008-5472.CAN-17-2748, PMID: 29437766

Gkountela S, Castro-Giner F, Szczerba BM, Vetter M, Landin J, Scherrer R, Krol I, Scheidmann MC, Beisel C, Stirnimann CU, Kurzeder C, Heinzelmann-Schwarz V, Rochlitz C, Weber WP, Aceto N. 2019. Circulating tumor cell clustering shapes DNA methylation to enable metastasis seeding. *Cell* **176**:98–112. DOI: https://doi.org/10.1016/j.cell.2018.11.046, PMID: 30633912

Gravenmier CA, Siddique M, Gatenby RA. 2018. Adaptation to stochastic temporal variations in intratumoral blood flow: the warburg effect as a bet hedging strategy. *Bulletin of Mathematical Biology* **80**:954–970. DOI: https://doi.org/10.1007/s11538-017-0261-x, PMID: 28508297

Greaves M. 2013. Cancer stem cells as 'units of selection'. *Evolutionary Applications* **6**:102–108. DOI: https://doi.org/10.1111/eva.12017, PMID: 23396760

Grunewald TG, Herbst SM, Heinze J, Burdach S. 2011. Understanding tumor heterogeneity as functional compartments-superorganisms revisited. *Journal of Translational Medicine* **9**:79. DOI: https://doi.org/10.1186/1479-5876-9-79, PMID: 21619636

Gupta PB, Fillmore CM, Jiang G, Shapira SD, Tao K, Kuperwasser C, Lander ES. 2011. Stochastic state transitions give rise to phenotypic equilibrium in populations of Cancer cells. *Cell* **146**:633–644. DOI: https://doi.org/10.1016/j.cell.2011.07.026, PMID: 21854987

Gupta PB, Pastushenko I, Skibinski A, Blanpain C, Kuperwasser C. 2019. Phenotypic plasticity: driver of Cancer initiation, progression, and therapy resistance. *Cell Stem Cell* **24**:65–78. DOI: https://doi.org/10.1016/j.stem.2018.11.011, PMID: 30554963

Harmelin A, Pinthus JH, Katzir N, Kapon A, Volcani Y, Amariglio EN, Rehavi G. 2001. Use of a murine xenograft model for canine transmissible venereal tumor. *American Journal of Veterinary Research* **62**:907–911. DOI: https://doi.org/10.2460/ajvr.2001.62.907, PMID: 11400849

Hausser J, Szekely P, Bar N, Zimmer A, Sheftel H, Caldas C, Alon U. 2019. Tumor diversity and the trade-off between universal Cancer tasks. *Nature Communications* **10**:5423. DOI: https://doi.org/10.1038/s41467-019-13195-1, PMID: 31780652

Hausser J, Alon U. 2020. Tumour heterogeneity and the evolutionary trade-offs of Cancer. *Nature Reviews Cancer* **20**:247–257. DOI: https://doi.org/10.1038/s41568-020-0241-6, PMID: 32094544

Heeke S, Mograbi B, Alix-Panabières C, Hofman P. 2019. Never travel alone: the crosstalk of circulating tumor cells and the blood microenvironment. *Cells* **8**:714. DOI: https://doi.org/10.3390/cells8070714

Herath S, Razavi Bazaz S, Monkman J, Ebrahimi Warkiani M, Richard D, O'Byrne K, Kulasinghe A. 2020. Circulating tumor cell clusters: insights into tumour dissemination and metastasis. *Expert Review of Molecular Diagnostics* **20**:1139–1147. DOI: https://doi.org/10.1080/14737159.2020.1846523, PMID: 33140979
Holmes JM. 1981. Measurement of the rate of death of canine transmissible venereal tumour cells transplanted into dogs and nude mice. *Research in Veterinary Science* **30**:248-250. DOI: https://doi.org/10.1016/S0034-5288(18)32592-X, PMID: 7255917

House M. 1997a. Transmission of Disseminated Neoplasia in the Soft Shell Clam, Mya Arenaria [Dissertation]. Oregon State University.

House M. 1997b. Transmission of Disseminated Neoplasia in the Soft Shell Clam Mya Arenaria. Oregon State University.

Huang S, Ernberg I, Kauffman S. 2009. Cancer attractors: a systems view of tumors from a gene network dynamics and developmental perspective. *Seminars in Cell & Developmental Biology* **20**:869-876. DOI: https://doi.org/10.1016/j.semcdb.2009.07.003, PMID: 19595782

Hurtado P, Martínez-Pena I, Piñeiro R. 2020. Dangerous liaisons: circulating tumor cells (CTCs) and Cancer-Associated fibroblasts (CAFs). *Cancers* **12**:2861. DOI: https://doi.org/10.3390/cancers12102861

Inda MM, Bonavia R, Mukasa A, Narita Y, Sah DW, Vandenberg S, Brennan C, Johns TG, Bachoo R, Hadwiger P, Tan P, Depinho RA, Cavenee W, Furnari F. 2010. Tumor heterogeneity is an active process maintained by a mutant EGFR-induced cytokine circuit in glioblastoma. *Genes & Development* **24**:1731-1745. DOI: https://doi.org/10.1101/gad.1890510, PMID: 20713517

Jackson HW, Fischer JR, Zanotelli VRT, Ali HR, Mechera R, Soysal SD, Moch H, Muenst S, Varga Z, Weber WP, Bodenmiller B. 2020. The single-cell pathology landscape of breast Cancer. *Nature* **578**:615-620. DOI: https://doi.org/10.1038/s41586-019-1876-x, PMID: 31959985

James S, Jennings G, Kwon YM, Stammnitz M, Fraik A, Storfer A, Comte S, Pemberton D, Fox S, Brown B, Pye R, Woods G, Lyons B, Hohenlohe PA, McCallum H, Siddle H, Thomas F, Ujvari B, Murchison EP, Jones M, et al. 2019. Tracing the rise of malignant cell lines: distribution, epidemiology and evolutionary interactions of two transmissible cancers in tasmanian devils. *Evolutionary Applications* **12**:1772-1780. DOI: https://doi.org/10.1111/eva.12831, PMID: 31548856

Janiszewska M, Tabassum DP, Castaño Z, Cristea S, Yamamoto KN, Kingston NL, Murphy KC, Shu S, Harper NW, Del Alcazar CG, Alečković M, Ekram MB, Cohen O, Kwak M, Qin Y, Laszewski T, Luoma A, Marusyk A, Wucherpfennig KW, Wagle N, et al. 2019. Subclonal cooperation drives metastasis by modulating local and systemic immune microenvironments. *Nature Cell Biology* **21**:879-888. DOI: https://doi.org/10.1038/s41556-019-0346-x, PMID: 31263265

Jenkinson G, Pujadas E, Goutsias J, Feinberg AP. 2017. Potential energy landscapes identify the information-theoretic nature of the epigenome. *Nature Genetics* **49**:719-729. DOI: https://doi.org/10.1038/ng.3811, PMID: 28346445

Jia D, Jolly MK, Kulkarni P, Levine H. 2017. Phenotypic plasticity and cell fate decisions in Cancer: insights from dynamical systems theory. *Cancers* **9**:70. DOI: https://doi.org/10.3390/cancers9070070

Johnson KE, Howard G, Mo W, Strasser MK, Lima E, Huang S, Brock A. 2019. Cancer cell population growth kinetics at low densities deviate from the exponential growth model and suggest an allee effect. *PLOS Biology* **17**:e3000399. DOI: https://doi.org/10.1371/journal.pbio.3000399, PMID: 31381560

Jolles JW, King AJ, Killen SS. 2020. The role of individual heterogeneity in collective animal behaviour. *Trends in Ecology & Evolution* **35**:278-291. DOI: https://doi.org/10.1016/j.tree.2019.11.001, PMID: 31879039

Koenig WD, Mumme RL, Groups LI. 2020. Population Ecology of the Cooperatively Breeding Acorn Woodpecker (MPB-24). **24** Princeton University Press. DOI: https://doi.org/10.2307/j.ctvx5wc1m

Korolev KS, Xavier JB, Gore J. 2014. Turning ecology and evolution against Cancer. *Nature Reviews Cancer* **14**: 371-380. DOI: https://doi.org/10.1038/nrc3712, PMID: 24739582

Krakauer DC. 2004. Groups confuse predators by exploiting perceptual bottlenecks: a connectionist model of the confusion effect. *Behavioral Ecology and Sociobiology* **36**:421-429. DOI: https://doi.org/10.1007/BF00177338

Kreiss A, Tovar C, Obendorf DL, Dun K, Woods GM. 2011. A murine xenograft model for a transmissible Cancer in tasmanian devils. *Veterinary Pathology* **48**:475-481. DOI: https://doi.org/10.1177/0300985810380398, PMID: 20861503

Kwon YM, Gori K, Park N, Potts N, Swift K, Wang J, Stammnitz MR, Cannell N, Baez-Ortega A, Comte S, Fox S, Harmsen C, Huxtable S, Jones M, Kreiss A, Lawrence C, Lazenby B, Peck S, Pye R, Woods G, et al. 2020. Evolution and lineage dynamics of a transmissible Cancer in tasmanian devils. *PLOS Biology* **18**:e3000926. DOI: https://doi.org/10.1371/journal.pbio.3000926, PMID: 33232318

Labuschagne CF, Cheung EC, Blagih J, Domart MC, Vousden KH. 2019. Cell clustering promotes a metabolic switch that supports metastatic colonization. *Cell Metabolism* **30**:720-734. DOI: https://doi.org/10.1016/j.cmet.2019.07.014, PMID: 31447323

Laconi E, Marongiu F, DeGregori J. 2020. Cancer as a disease of old age: changing mutational and microenvironmental landscapes. *British Journal of Cancer* **122**:943-952. DOI: https://doi.org/10.1038/s41416-019-0721-1, PMID: 32042067

Landeau L, Terborgh J. 1986. Oddity and the 'confusion effect' in predation. *Animal Behaviour* **34**:1372-1380. DOI: https://doi.org/10.1016/S0003-3472(86)80208-1

Lean C, Plutynski A. 2016. The evolution of failure: explaining Cancer as an evolutionary process. *Biology & Philosophy* **31**:39-57. DOI: https://doi.org/10.1007/s10539-015-9511-1

Lee JY, Yoon JK, Kim B, Kim S, Kim MA, Lim H, Bang D, Song YS. 2015. Tumor evolution and intratumor heterogeneity of an epithelial ovarian Cancer investigated using next-generation sequencing. *BMC Cancer* **15**: 85. DOI: https://doi.org/10.1186/s12885-015-1077-4, PMID: 25881093
Lewontin RC. 1970. The units of selection. *Annual Review of Ecology and Systematics* **1**:1–18. DOI: https://doi.org/10.1146/annurev.es.01.110170.000245

Li Q, Wennborg A, Aurell E, Dekel E, Zou JZ, Xu Y, Huang S, Ernberg I. 2016. Dynamics inside the Cancer cell attractor reveal cell heterogeneity, limits of stability, and escape. *PNAS* **113**:2672–2677. DOI: https://doi.org/10.1073/pnas.1519210113, PMID: 26929366

Liau BB, Sievers C, Donohue LK, Gillespie SM, Flavahan WA, Miller TE, Venteicher AS, Hebert CH, Carey CD, Rodig SJ, Shareef SJ, Najm FJ, van Galen P, Wakimoto H, Cahill DP, Rich JN, Aster JC, Suvà ML, Patel AP, Bernstein BE. 2017. Adaptive chromatin remodeling drives glioblastoma stem cell plasticity and drug tolerance. *Cell Stem Cell* **20**:233–246. DOI: https://doi.org/10.1016/j.stem.2016.11.003, PMID: 27989769

Losi L, Baisse B, Bouzourene H, Benhattar J. 2005. Evolution of intratumoral genetic heterogeneity during colorectal Cancer progression. *Carcinogenesis* **26**:916–922. DOI: https://doi.org/10.1093/carcin/bgi044, PMID: 15731168

Martín-Pardillos A, Valls Chiva Á, Bande Vargas G, Hurtado Blanco P, Piñeiro Cid R, Guijarro PJ, Hümmer S, Bejar Serrano E, Rodriguez-Casanova A, Diaz-Lagares Á, Castellvi J, Miravet-Verde S, Serrano L, Lluch-Senar M, Sebastian V, Bribian A, López-Mascaraque L, López-López R, Ramón Y Cajal S. 2019. The role of clonal communication and heterogeneity in breast Cancer. *BMC Cancer* **19**:666. DOI: https://doi.org/10.1186/s12885-019-5883-y, PMID: 31277602

Marusyk A, Almendro V, Polyak K. 2012. Intra-tumour heterogeneity: a looking glass for Cancer? *Nature Reviews Cancer* **12**:323–334. DOI: https://doi.org/10.1038/nrc3261, PMID: 22513401

Marusyk A, Tabassum DP, Altrock PM, Almendro V, Michor F, Polyak K. 2014. Non-cell-autonomous driving of tumour growth supports sub-clonal heterogeneity. *Nature* **514**:54–58. DOI: https://doi.org/10.1038/nature13556, PMID: 25079331

Marusyk A, Janiszewska M, Polyak K. 2020. Intratumor heterogeneity: the rosetta stone of therapy resistance. *Cancer Cell* **37**:471–484. DOI: https://doi.org/10.1016/j.ccell.2020.03.007, PMID: 32289271

Marusyk A, Polyak K. 2010. Tumor heterogeneity: causes and consequences. *Biochimica Et Biophysica Acta (BBA) - Reviews on Cancer* **1805**:105–117. DOI: https://doi.org/10.1016/j.bbcan.2009.11.002, PMID: 19931353

McCallum H, Jones M, Hawkins C, Hamede R, Lachish S, Sinn DL, Beeton N, Lazenby B. 2009. Transmission dynamics of tasmanian Devil facial tumor disease may lead to disease-induced extinction. *Ecology* **90**:3379–3392. DOI: https://doi.org/10.1890/08-1763.1, PMID: 20120807

McGranahan N, Swanton C, Heterogeneity C. 2017. Tumor evolution: past, present, and the future. *Cell* **168**: 613–628. DOI: https://doi.org/10.1016/j.cell.2017.01.018

McGranahan N, Swanton C. 2015. Biological and therapeutic impact of intratumor heterogeneity in Cancer evolution. *Cancer Cell* **27**:15–26. DOI: https://doi.org/10.1016/j.ccell.2014.12.001, PMID: 25584892

Meacham CE, Morrison SJ. 2013. Tumour heterogeneity and Cancer cell plasticity. *Nature* **501**:328–337. DOI: https://doi.org/10.1038/nature12624, PMID: 24048065

Merlo LM, Pepper JW, Reid BJ, Maley CC. 2006. Cancer as an evolutionary and ecological process. *Nature Reviews Cancer* **6**:924–935. DOI: https://doi.org/10.1038/nrc2013, PMID: 17109012

Metzger MJ, Villalba A, Carballal MJ, Iglesias D, Sherry J, Reinisch C, Muttray AF, Baldwin SA, Goff SP. 2016. Widespread transmission of independent Cancer lineages within multiple bivalve species. *Nature* **534**:705–709. DOI: https://doi.org/10.1038/nature18599, PMID: 27338791

Morris LG, Riaz N, Desrichard A, Şenbabaoğlu Y, Hakimi AA, Makarov V, Reis-Filho JS, Chan TA. 2016. Pan-cancer analysis of intratumor heterogeneity as a prognostic determinant of survival. *Oncotarget* **7**:10051–10063. DOI: https://doi.org/10.18632/oncotarget.7067, PMID: 26840267

Munn DH, Mellor AL. 2016. IDO in the tumor microenvironment: inflammation, Counter-Regulation, and tolerance. *Trends in Immunology* **37**:193–207. DOI: https://doi.org/10.1016/j.it.2016.01.002, PMID: 26839260

Murlidhar V, Reddy RM, Fouladdel S, Zhao L, Ishikawa MK, Grabauskiene S, Zhang Z, Lin J, Chang AC, Carrott P, Lynch WR, Orringer MB, Kumar-Sinha C, Palanisamy N, Beer DG, Wicha MS, Ramnath N, Azizi E, Nagrath S. 2017. Poor prognosis indicated by venous circulating tumor cell clusters in Early-Stage lung cancers. *Cancer Research* **77**:5194–5206. DOI: https://doi.org/10.1158/0008-5472.CAN-16-2072, PMID: 28716896

Neelakantan D, Zhou H, Oliphant MUJ, Zhang X, Simon LM, Henke DM, Shaw CA, Wu MF, Hilsenbeck SG, White LD, Lewis MT, Ford HL. 2017. EMT cells increase breast Cancer metastasis via paracrine GLI activation in neighbouring tumour cells. *Nature Communications* **8**:15773. DOI: https://doi.org/10.1038/ncomms15773, PMID: 28604738

Nguyen LV, Vanner R, Dirks P, Eaves CJ. 2012. Cancer stem cells: an evolving concept. *Nature Reviews Cancer* **12**:133–143. DOI: https://doi.org/10.1038/nrc3184, PMID: 22237392

Nguyen A, Yoshida M, Goodarzi H, Tavazoie SF. 2016. Highly variable Cancer subpopulations that exhibit enhanced transcriptome variability and metastatic fitness. *Nature Communications* **7**:11246. DOI: https://doi.org/10.1038/ncomms11246, PMID: 27138336

O'Neill ID. 2011. Concise review: transmissible animal tumors as models of the Cancer stem-cell process. *Stem Cells* **29**:1909–1914. DOI: https://doi.org/10.1002/stem.751, PMID: 21956952

Padmanaban V, Krol I, Suhail Y, Szczerba BM, Aceto N, Bader JS, Ewald AJ. 2019. E-cadherin is required for metastasis in multiple models of breast Cancer. *Nature* **573**:439–444. DOI: https://doi.org/10.1038/s41586-019-1526-3, PMID: 31485072

Pantel K, Alix-Panabières C. 2019. Liquid biopsy and minimal residual disease - latest advances and implications for cure. *Nature Reviews Clinical Oncology* **16**:409–424. DOI: https://doi.org/10.1038/s41571-019-0187-3, PMID: 30796368
Pastore A, Gaiti F, Lu SX, Brand RM, Kulm S, Chaligne R, Gu H, Huang KY, Stamenova EK, Béguelin W, Jiang Y, Schulman RC, Kim KT, Alonso A, Allan JN, Furman RR, Gnrirke A, Wu CJ, Melnick AM, Meissner A, et al. 2019. Corrupted coordination of epigenetic modifications leads to diverging chromatin states and transcriptional heterogeneity in CLL. *Nature Communications* **10**:1874. DOI: https://doi.org/10.1038/s41467-019-09645-5, PMID: 31015400

Patel RR, Verma V, Barsoumian HB, Ning MS, Chun SG, Tang C, Chang JY, Lee PP, Gandhi S, Balter P, Dunn JD, Chen D, Puebla-Osorio N, Cortez MA, Welsh JW. 2021. Use of Multi-Site radiation therapy for systemic disease control. *International Journal of Radiation Oncology*Biology*Physics **109**:352–364. DOI: https://doi.org/10.1016/j.ijrobp.2020.08.025, PMID: 32798606

Pearse AM, Swift K, Hodson P, Hua B, McCallum H, Pyecroft S, Taylor R, Eldridge MD, Belov K. 2012. Evolution in a transmissible Cancer: a study of the chromosomal changes in Devil facial tumor (DFT) as it spreads through the wild tasmanian Devil population. *Cancer Genetics* **205**:101–112. DOI: https://doi.org/10.1016/j.cancergen.2011.12.001, PMID: 22469509

Pepper JW. 2012. Drugs that target pathogen public goods are robust against evolved drug resistance. *Evolutionary Applications* **5**:757–761. DOI: https://doi.org/10.1111/j.1752-4571.2012.00254.x, PMID: 23144661

Pisco AO, Huang S. 2015. Non-genetic Cancer cell plasticity and therapy-induced stemness in tumour relapse: 'What does not kill me strengthens me'. *British Journal of Cancer* **112**:1725-1732. DOI: https://doi.org/10.1038/bjc.2015.146, PMID: 25965164

Pye RJ, Pemberton D, Tovar C, Tubio JM, Dun KA, Fox S, Darby J, Hayes D, Knowles GW, Kreiss A, Siddle HV, Swift K, Lyons AB, Murchison EP, Woods GM. 2016. A second transmissible Cancer in tasmanian devils. *PNAS* **113**:374–379. DOI: https://doi.org/10.1073/pnas.1519691113, PMID: 26711993

Quail DF, Joyce JA. 2013. Microenvironmental regulation of tumor progression and metastasis. *Nature Medicine* **19**:1423–1437. DOI: https://doi.org/10.1038/nm.3394, PMID: 24202395

Ramón Y Cajal S, Capdevila C, Hernandez-Losa J, De Mattos-Arruda L, Ghosh A, Lorent J, Larsson O, Aasen T, Postovit LM, Topisirovic I. 2017. Cancer as an ecomolecular disease and a neoplastic consortium. *Biochimica Et Biophysica Acta (BBA) - Reviews on Cancer* **1868**:484–499. DOI: https://doi.org/10.1016/j.bbcan.2017.09.004, PMID: 28947238

Rios AC, Capaldo BD, Vaillant F, Pal B, van Ineveld R, Dawson CA, Chen Y, Nolan E, Fu NY, Jackling FC, Devi S, Clouston D, Whitehead L, Smyth GK, Mueller SN, Lindeman GJ, Visvader JE, 3DTCLSM Group. 2019. Intraclonal plasticity in mammary tumors revealed through Large-Scale Single-Cell resolution 3D imaging. *Cancer Cell* **35**:618–632. DOI: https://doi.org/10.1016/j.ccell.2019.02.010, PMID: 30930118

Rodgers GM, Downing B, Morrell LJ. 2015. Prey body size mediates the predation risk associated with being "odd". *Behavioral Ecology* **26**:242–246. DOI: https://doi.org/10.1093/beheco/aru185

Saito T, Niida A, Uchi R, Hirata H, Komatsu H, Sakimura S, Hayashi S, Nambara S, Kuroda Y, Ito S, Eguchi H, Masuda T, Sugimachi K, Tobo T, Nishida H, Daa T, Chiba K, Shiraishi Y, Yoshizato T, Kodama M, et al. 2018. A temporal shift of the evolutionary principle shaping intratumor heterogeneity in colorectal Cancer. *Nature Communications* **9**:2884. DOI: https://doi.org/10.1038/s41467-018-05226-0, PMID: 30038269

Sewalt L, Harley K, van Heijster P, Balasuriya S. 2016. Influences of allee effects in the spreading of malignant tumours. *Journal of Theoretical Biology* **394**:77–92. DOI: https://doi.org/10.1016/j.jtbi.2015.12.024, PMID: 26802481

Shaffer SM, Dunagin MC, Torborg SR, Torre EA, Emert B, Krepler C, Beqiri M, Sproesser K, Brafford PA, Xiao M, Eggan E, Anastopoulos IN, Vargas-Garcia CA, Singh A, Nathanson KL, Herlyn M, Raj A. 2017. Rare cell variability and drug-induced reprogramming as a mode of Cancer drug resistance. *Nature* **546**:431–435. DOI: https://doi.org/10.1038/nature22794, PMID: 28607484

Sharma A, Cao EY, Kumar V, Zhang X, Leong HS, Wong AML, Ramakrishnan N, Hakimullah M, Teo HMV, Chong FT, Chia S, Thangavelu MT, Kwang XL, Gupta R, Clark JR, Periyasamy G, Iyer NG, DasGupta R. 2018. Longitudinal single-cell RNA sequencing of patient-derived primary cells reveals drug-induced infidelity in stem cell hierarchy. *Nature Communications* **9**:4931. DOI: https://doi.org/10.1038/s41467-018-07261-3, PMID: 30467425

Sunila I. 1991. Respiration of sarcoma cells from the soft-shell clam Mya arenaria L. under various conditions. *Journal of Experimental Marine Biology and Ecology* **150**:19–29. DOI: https://doi.org/10.1016/0022-0981(91)90103-4

Sunila I, Farley CA. 1989. Environmental limits for survival of sarcoma cells from the soft-shell clam Mya arenaria. *Diseases of Aquatic Organisms* **7**:111–115. DOI: https://doi.org/10.3354/dao007111

Szczerba BM, Castro-Giner F, Vetter M, Krol I, Gkountela S, Landin J, Scheidmann MC, Donato C, Scherrer R, Singer J, Beisel C, Kurzeder C, Heinzelmann-Schwarz V, Rochlitz C, Weber WP, Beerenwinkel N, Aceto N. 2019. Neutrophils escort circulating tumour cells to enable cell cycle progression. *Nature* **566**:553–557. DOI: https://doi.org/10.1038/s41586-019-0915-y, PMID: 30728496

Tabassum DP, Polyak K. 2015. Tumorigenesis: it takes a village. *Nature Reviews Cancer* **15**:473–483. DOI: https://doi.org/10.1038/nrc3971, PMID: 26156638

Thomas F, Kareva I, Raven N, Hamede R, Pujol P, Roche B, Ujvari B. 2018. Evolved dependence in response to Cancer. *Trends in Ecology & Evolution* **33**:269–276. DOI: https://doi.org/10.1016/j.tree.2018.01.012, PMID: 29475741

Ujvari B, Gatenby RA, Thomas F. 2016a. The evolutionary ecology of transmissible cancers. *Infection, Genetics and Evolution* **39**:293–303. DOI: https://doi.org/10.1016/j.meegid.2016.02.005, PMID: 26861618

Ujvari B, Gatenby RA, Thomas F. 2016b. Transmissible cancers, are they more common than thought? *Evolutionary Applications* **9**:633–634. DOI: https://doi.org/10.1111/eva.12372, PMID: 27398093
Ujvari B, Gatenby RA, Thomas F. 2017. Chapter 12 - Transmissible Cancer: The Evolution of Interindividual Metastasis. In: Ujvari B, Roche B, Thomas F (Eds). *Ecology and Evolution of Cancer: Academic Press*. Wiley. p. 167–179. DOI: https://doi.org/10.1002/icl3.1020

Vibishan B, Watve MG. 2020. Context-dependent selection as the keystone in the somatic evolution of Cancer. *Scientific Reports* **10**:4223. DOI: https://doi.org/10.1038/s41598-020-61046-7, PMID: 32144283

Vicente-Dueñas C, Hauer J, Ruiz-Roca L, Ingenhag D, Rodríguez-Meira A, Auer F, Borkhardt A, Sánchez-García I. 2015. Tumoral stem cell reprogramming as a driver of Cancer: theory, biological models, implications in Cancer therapy. *Seminars in Cancer Biology* **32**:3–9. DOI: https://doi.org/10.1016/j.semcancer.2014.02.001, PMID: 24530939

Vinci M, Burford A, Molinari V, Kessler K, Popov S, Clarke M, Taylor KR, Pemberton HN, Lord CJ, Gutteridge A, Forshew T, Carvalho D, Marshall LV, Qin EY, Ingram WJ, Moore AS, Ng HK, Trabelsi S, H'mida-Ben Brahim D, Entz-Werle N, et al. 2018. Functional diversity and cooperativity between subclonal populations of pediatric glioblastoma and diffuse intrinsic pontine glioma cells. *Nature Medicine* **24**:1204–1215. DOI: https://doi.org/10.1038/s41591-018-0086-7, PMID: 29967352

Wagenblast E, Soto M, Gutiérrez-Ángel S, Hartl CA, Gable AL, Maceli AR, Erard N, Williams AM, Kim SY, Dickopf S, Harrell JC, Smith AD, Perou CM, Wilkinson JE, Hannon GJ, Knott SR. 2015. A model of breast Cancer heterogeneity reveals vascular mimicry as a driver of metastasis. *Nature* **520**:358–362. DOI: https://doi.org/10.1038/nature14403, PMID: 25855289

Wagle N, Emery C, Berger MF, Davis MJ, Sawyer A, Pochanard P, Kehoe SM, Johannessen CM, Macconaill LE, Hahn WC, Meyerson M, Garraway LA. 2011. Dissecting therapeutic resistance to RAF inhibition in melanoma by tumor genomic profiling. *Journal of Clinical Oncology* **29**:3085–3096. DOI: https://doi.org/10.1200/JCO.2010.33.2312, PMID: 21383288

Wagner J, Rapsomaniki MA, Chevrier S, Anzeneder T, Langwieder C, Dykgers A, Rees M, Ramaswamy A, Muenst S, Soysal SD, Jacobs A, Windhager J, Silina K, van den Broek M, Dedes KJ, Rodríguez Martínez M, Weber WP, Bodenmiller B. 2019. A Single-Cell atlas of the tumor and immune ecosystem of human breast Cancer. *Cell* **177**:1330–1345. DOI: https://doi.org/10.1016/j.cell.2019.03.005, PMID: 30982598

Wu M, Pastor-Pareja JC, Xu T. 2010. Interaction between ras(V12) and scribbled clones induces tumour growth and invasion. *Nature* **463**:545–548. DOI: https://doi.org/10.1038/nature08702, PMID: 20072127

Xiong S, Feng Y, Cheng L. 2019. Cellular reprogramming as a therapeutic target in Cancer. *Trends in Cell Biology* **29**:623–634. DOI: https://doi.org/10.1016/j.tcb.2019.05.001, PMID: 31153655

Yonemitsu MA, Giersch RM, Polo-Prieto M, Hammel M, Simon A, Cremonte F, Avilés FT, Merino-Véliz N, Burioli EA, Muttray AF, Sherry J, Reinisch C, Baldwin SA, Goff SP, Houssin M, Arriagada G, Vázquez N, Bierne N, Metzger MJ. 2019. A single clonal lineage of transmissible Cancer identified in two marine mussel species in South America and Europe. *eLife* **8**:e47788. DOI: https://doi.org/10.7554/eLife.47788, PMID: 31686650

Yu M, Bardia A, Wittner BS, Stott SL, Smas ME, Ting DT, Isakoff SJ, Ciciliano JC, Wells MN, Shah AM, Concannon KF, Donaldson MC, Sequist LV, Brachtel E, Sgroi D, Baselga J, Ramaswamy S, Toner M, Haber DA, Maheswaran S. 2013. Circulating breast tumor cells exhibit dynamic changes in epithelial and mesenchymal composition. *Science* **339**:580–584. DOI: https://doi.org/10.1126/science.1228522, PMID: 23372014

Zhou H, Neelakantan D, Ford HL. 2017. Clonal cooperativity in heterogenous cancers. *Seminars in Cell & Developmental Biology* **64**:79–89. DOI: https://doi.org/10.1016/j.semcdb.2016.08.028, PMID: 27582427