1.1. Varying Degrees of Tumour Heterogeneity
Tumour heterogeneity harbours multiple layers of complexity in human malignancies. It has long been known that tumours of the same histopathological subtype commonly differ from one patient to another (inter-tumour heterogeneity) (
Figure 1). Exacerbating the complexity even further, sizable variations have been reported within a single tumour (intra-tumour heterogeneity, ITH). ITH can be detected between the different geographic regions of the same primary tumour or even between the primary tumour and the metastastic lesions (spatial intra-tumour heterogeneity). Moreover, the analysis of serial tumour samples demonstrated that the cell features may evolve during the course of the disease progression (temporal heterogeneity) under environmental or therapeutic stress
[1][2]. ITH has been observed in most (nearly all) types of cancers, including both haematological malignancies (chronic lymphoblastic leukemia and acute lymphoblastic leukemia), and solid tumours (lung, breast, ovarian, pancreatic, kidney, colorectal, brain and prostate cancers)
[3].
Figure 1. A multifaceted heterogeneity in cancers. (A) Inter-tumour heterogeneity refers to the variability observed in tumours of the same histological subtypes between different patients. (B) Intra-tumour heterogeneity (ITH) is observed across different regions of the primary tumour site and/or metastatic sites (spatial ITH) and can evolve over time (temporal ITH). Colours represent the different characteristics between tumours or tumour cells.
1.1.1. Phenotypic Heterogeneity
The first demonstration of tumour heterogeneity has been made by histopathologists who are familiar with morphological divergence (differentiation status, necrosis, fibrosis and so on) across the tumours or between the different areas of the tumour (
Figure 2)
[1][4][5][6][7][8]. This notion has led to the very basis of tumour classification systems based on histopathological features
[9]. Tumour grading systems notably include the pathological examination of multiple microscopy fields in order to avoid tumour misclassification due to ITH
[10]. Increasing evidence indicates that tumour foci are heterogeneous at other phenotypic levels than merely morphologic, including differential capabilities in terms of proliferation, metabolism, motility, migration, invasiveness, metastasis and stemness, as well as varied sensitivity to therapies
[11][12][13]. The morphological and other phenotypic cell features co-vary in the different tumour regions, notably between the core and the external borders of the tumour.
Figure 2. Sources of intra-tumour heterogeneity (ITH). Genetic, epigenetic, microenvironmental, transcriptomic and proteomic heterogeneities coexist in tumours and are linked with phenotypic diversity. Abbreviations: CAF: Cancer-associated fibroblasts; ECM: extracellular matrix; EMT: epithelial-to-mesenchymal transition; indels: small insertions and deletions; LOH: loss of heterozygoty; MET: mesenchymal-to-epithelial transition; SCNA: somatic copy number alterations; SNV: single nucleotide variants; SV: structural variants.
1.1.2. Molecular Heterogeneity
Advances in next-generation sequencing (NGS) revealed the extent and prevalence of molecular diversity in tumours
[14][15]. The sequencing of multiple regions in space and time demonstrated the various repertoires of genetic events that can occur in cancers, including single nucleotide variants (SNVs), small insertions and deletions (indels), structural variants and somatic copy number alterations (SCNA)
[15]. Large-scale studies indicated that genetic ITH occurs in almost all cancer types, albeit at varying degrees
[16][17]. Melanoma and lung adenocarcinomas notably account for cancers with high mutational tumour burden and the establishment of specific mutational signatures as a result of exposure to exogenous mutagens (UV radiation and tobacco smoke)
[18]. Dietz et al. demonstrated that the frequencies of driver gene mutations in regionally distinct areas of lung adenocarcinomas were correlated with the spatial distribution of histological patterns, highlighting an interplay between histologic and genetic features in a tumour
[19][20].
However, the genetic perspective is insufficient to fully explain the range of phenotypic diversity in solid malignancies, given the fact that cell populations with identical genetic background can lead to distinct morphological patterns and differential responses to treatment or environmental stimuli
[10][21]. Emerging evidence demonstrated that ITH also take place at other levels, such as epigenetics, transcriptomics and proteomics
[22].
The epigenome is defined as a connection between the genome and the environment. Alterations of the epigenetic machinery has been recognised as a hallmark of cancer
[23] and may appear early during carcinogenesis
[24]. Epigenetic marks induce heritable changes in gene expression without any modification in the underlying DNA sequence that allows cells to adapt to microenvironment stimuli (oxygen, nutrient deprivation, acidity and so on) and develop resistance mechanisms against anticancer therapies
[25]. Bidirectional communications between genetics and epigenetics have been reported in cancers, with the detection of somatic mutations in genes encoding epigenome regulators (such as
DNMT3A, IDH1, H3F3A) and inversely the identification of DNA hypomethylation or epigenetic silencing of DNA repair genes (such as
MLH1 or
BRCA1) that can cause genomic instability in cancer cells
[24]. It was assessing histone modifications, chromatin accessibility and DNA methylation profiles demonstrated a high epigenetic variability in cancers
[21][26][27]. Considering the major implications of epigenetics in the development of cancers and their response to anticancer treatments, a better understanding of epigenetic heterogeneity could help to identify novel epigenetic therapies and consider them for a combination with other anticancer treatments (genotoxic/cytotoxic agents, hormone therapy, immunotherapy, targeted therapy) to improve their efficacy or reverse drug resistance
[25].
Transcriptome refers to all RNA species that can be found in cells; however, mRNAs are frequently the most studied. Their composition varies between cell types and tumour types and continuously evolves depending on the local conditions that are applied to cells over time. They can be explored through targeted (RT-PCR) or high throughput approaches (gene expression arrays, RNA sequencing (RNASeq)). A plethora of gene expression signatures have been developed in oncology for tumour classification
[28], prognosis establishment
[29][30][31][32], therapeutic and surveillance decision making
[33] but only a few are already implemented for routine practice
[34].
Because proteins directly reveal the functional mechanisms that occur in cancers and account for most of the therapeutic targets, it appears important to assess tumour heterogeneity at the protein level, which has shown growing interest. Proteomic approaches have long lagged behind those for transcriptome and genome due to technical limitations, high amounts of proteins generated from a single gene (with different isoforms and modification states) and a complex regulation of protein expression at both translational and post-translational levels
[35]. Immunohistochemistry appears as one of the most standard approaches to assess protein abundance changes; however, it provides only semi-quantitative information, interrogates a limited number of proteins and is limited by the availability of appropriate antibodies. The development of reverse-phase protein array (RPPA) and mass spectrometry (MS)-based methods enabled the assessment of the proteomic landscape on a larger scale
[35]. Transcriptomic approaches cannot substitute proteomic investigations, as the analysis of datasets from The Cancer Proteome Atlas (TCPA) found a poor correlation between protein and gene expression in cancer tissues, with Spearman correlation oscillating from 0.1 to 0.3, depending on the cancer type
[36].
Recognition of all these cancer-specific molecular processes as major elements in the evolution of cancers and their considerable diversity has led to the launch of cancer genomics programs by international consortia, such as The Cancer Genome Atlas consortium (TCGA) and the International Cancer Genome Consortium (ICGC). All aim to integrate the analysis of multi-omics datasets (genetic, epigenetic, transcriptomic and proteomic data) to provide a comprehensive overview of the tumour landscape
[16][37][38][39][40][41][42][43][44][45].
1.1.3. Tumour Micro-Environment (TME) Heterogeneity
The influence of the complex ecosystem in which cancer cells evolve has long been overlooked. In the last decade, cancer biology progressively shifted from a cancer cell-centric model to a more ample view, where cancer cells and their near environment are highly interrelated. The tumour microenvironment (TME) is made up of non-transformed cells (endothelial cells, fibroblasts, pericytes, adipocytes, immune cells and so on) and non-cellular constituents (such as the extracellular matrix) which are shaped by cancer cells through the modification of local environmental conditions and the secretion of oncogenic signals
[46][47]. As a consequence, the phenotypic traits and behaviours of TME components are highly heterogeneous, depending on the tumour context
[46][48]. In return, TME can assist in the development of the tumour niche by contributing to cancer progression, metastasis and drug resistance
[48][49][50]. In this way, TME represents an emerging target for treatments (such as immune checkpoint inhibitors or antiangiogenic therapies) and should be taken into consideration for clinical decisions. Recently, Garattini et al. demonstrated that heterogeneity also extends to the drug distribution in tumours, which depends on many aspects of the patient, the tumour and its microenvironment and influences tumour response
[51].
1.2. Unravelling Evolutionary Processes behind Tumour Heterogeneity
Two major and paradoxical theories have been developed to explain the installation of high degree of diversity in tumours. In 1976, Peter Nowell first described the cancer development as a continuous evolutionary process originating from a single renagade mutant cell and driven by the accumulation of stepwise somatic mutations during proliferation processes that give rise to various clones and subclones
[52]. The development of multiple cell groups with distinct genomic profiles is amplified by genomic instability that arises in most solid tumours and haematopoietic malignancies as a result from both exposure to exogenous mutagens and defects in DNA repair pathways
[53]. A subclone is characterised as a set of cells that diverge from the cell ancestor lineage (clone) by the presence of additional genetic alterations. Equivalent to Darwinian natural selection, most stochastic events that appears during the evolution process probably do not confer any selective benefit to the cancer cells (passenger mutations)
[54]. In contrast, certain mutations can provide a fitness advantage over adjacent cells (driver mutations) and enables them to become predominant and outcompete other ones
[55]. Most driver mutations are clonal. They appear early during cancer progression under a given microenvironmental context and foster cancer progression but they seem not essential for cancer maintenance once installed
[54]. The clonal genomic architecture is distinct from a tumour to another considering that the emergence of subclones strongly depends on specific environmental stresses (local hypoxia or inflammation, treatment exposure, etc.) applied in each tumour over time. More recently, epigenetics and genetics were shown to follow convergent evolutionary trajectories in the development of cancers, highlighting the potential interest of combining epigenetic agents with other anticancer therapies
[56][57][58].
In contrast to the Darwinian clonal evolution theory, where all subclones possess tumorigenic potential, a second model proposed that only a small subgroup of cancer cells (named cancer stem cells (CSC) or tumour-initiating cells) has the capacity to generate new tumours
[59][60]. In this model, tumours are structured in a unidirectional hierarchy fashion, whereby CSC can either indefinitely self-renew (symmetric division) or differentiate into multiple cancer cell types (asymmetric division). CSC with stem cell-like characteristics have been observed in several cancers, including leukemia, breast, colon, head and neck and oesophageal cancers
[55]. CSC are thought to be more drug-resistant than non-CSC and in such ways, they may be responsible for recurrence and therapeutic evasion
[59]. Increasing evidence, however, indicated that non-CSC can readily convert to a CSC state through cell plasticity programmes, such as epithelial-to-mesenchymal transition (EMT), indicating that the hierarchy seems less rigid than previously thought
[61]. In the same manner, different subsets of CSC with variable EMT phenotypes can coexist in tumours and can switch from one to another
[61]. Stemness and CSC plasticity may be modulated by internal (genetic and epigenetic) and external (TME) factors that can work apart or simultaneously
[62]. Van Niekerk and colleagues show that certain stem cell features can be acquired by cancer cells through clonal selection, highlighting the fact that clonal evolution and the CSC theories are not necessarily mutually exclusive and can intertwine
[63].
In these first models of tumour evolution, ITH was thought to gain gradually over time as the tumour grew (
Figure 3). Although this concept of continuous clonal evolution is still applicable to describe most cancer evolutionary processes, increasing evidence supported the idea that this model cannot explain the full spectrum of observed evolutionary behaviours
[64]. Notably, single catastrophic events, such as whole-genome doubling, chromosomal chromoplexy and chromothripsis, can arise suddenly as single macroevolutionary jumps over long periods of relative stasis. In some extreme cases of punctuated tumour evolution, the development of colorectal cancers and other tumour types has been modelled as “Big-bang” dynamics, whereby a single or few mutational bursts occur early during carcinogenesis and result in a large number of intermixed subclones that are not subjected to selective pressure and coexist during growth (neutral evolution)
[65][66][67].
Figure 3. Models of tumour evolution described by Muller plots, which represent the tumour clonal dynamics over time. Colours indicate the different genotypes of the tumour cell clones.
Branched evolutionary trajectories have been extensively described in a wide range of tumour types, such as childhood acute lymphoblastic leukaemia, clear cell renal carcinomas, pancreatic, colorectal, breast and prostate cancers
[68]. Evolutionary pathways can then be represented as a phylogenetic tree, where truncal mutations (clonal) represent the alterations occurring early in cancer development progenitors, while nontruncal mutations (subclonal) emerge during cancer progression and are shared by only a small group but not all cancer cells. In a branched evolution pattern, several distinct subclones co-exist and can be either intermingled in the same area or regionally separated, depending on the presence of physical barriers, such as blood vessels or microenvironment specificities
[53].
ITH has also been described in cancer cases with linear evolutionary trajectories whereby a predominant subclone outgrows at the expense of its predecessor(s) followed by incomplete selective sweeps
[53][68]. Although most studies described a single model of evolution in cancers, emerging data suggest that tumours may follow different models of evolution (linear, branched, punctuated or neutral) sequentially or simultaneously during the course of the disease
[69]. The full context of tumour evolution is still to be explored in detail in order to better define effective therapeutic strategies.