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

Cetuximab chemotherapy resistance: Insight into the homeostatic evolution of head and neck cancer (Review)

The complex evolution of genetic alterations in cancer that occurs in vivo is a selective process involving numerous factors and mechanisms. Chemotherapeutic agents that prevent the growth and spread of cancer cells induce selective pressure, leading to rapid artificial selection of resistant subclones. This rapid evolution is possible because antineoplastic drugs promote alterations in tumor‑cell metabolism, thus creating a bottleneck event. The few resistant cells that survive in this new environment obtain differential reproductive success that enables them to pass down the newly selected resistant gene pool. The present review aims to summarize key findings of tumor evolution, epithelial‑mesenchymal transition and resistance to cetuximab therapy in head and neck squamous cell carcinoma.

Time-resolved, integrated analysis of clonally evolving genomes

Clonal genome evolution is a key feature of asexually reproducing species and human cancer development. While many studies have described the landscapes of clonal genome evolution in cancer, few determine the underlying evolutionary parameters from molecular data, and even fewer integrate theory with data. We derived theoretical results linking mutation rate, time, expansion dynamics, and biological/clinical parameters. Subsequently, we inferred time-resolved estimates of evolutionary parameters from mutation accumulation, mutational signatures and selection. We then applied this framework to predict the time of speciation of the marbled crayfish, an enigmatic, globally invasive parthenogenetic freshwater crayfish. The results predict that speciation occurred between 1986 and 1990, which is consistent with biological records. We also used our framework to analyze whole-genome sequencing datasets from primary and relapsed glioblastoma, an aggressive brain tumor. The results identified evolutionary subgroups and showed that tumor cell survival could be inferred from genomic data that was generated during the resection of the primary tumor. In conclusion, our framework allowed a time-resolved, integrated analysis of key parameters in clonally evolving genomes, and provided novel insights into the evolutionary age of marbled crayfish and the progression of glioblastoma.

Mutation divergence over space in tumour expansion

Mutation accumulation in tumour evolution is one major cause of intra-tumour heterogeneity (ITH), which often leads to drug resistance during treatment. Previous studies with multi-region sequencing have shown that mutation divergence among samples within the patient is common, and the importance of spatial sampling to obtain a complete picture in tumour measurements. However, quantitative comparisons of the relationship between mutation heterogeneity and tumour expansion modes, sampling distances as well as the sampling methods are still few. Here, we investigate how mutations diverge over space by varying the sampling distance and tumour expansion modes using individual-based simulations. We measure ITH by the Jaccard index between samples and quantify how ITH increases with sampling distance, the pattern of which holds in various sampling methods and sizes. We also compare the inferred mutation rates based on the distributions of variant allele frequencies under different tumour expansion modes and sampling sizes. In exponentially fast expanding tumours, a mutation rate can always be inferred for any sampling size. However, the accuracy compared with the true value decreases when the sampling size decreases, where small sampling sizes result in a high estimate of the mutation rate. In addition, such an inference becomes unreliable when the tumour expansion is slow, such as in surface growth.

Tumour Genetic Heterogeneity in Relation to Oral Squamous Cell Carcinoma and Anti-Cancer Treatment

Oral squamous cell carcinoma (SCC) represents more than 90% of all oral cancers and is the most frequent SCC of the head and neck region. It may affect any oral mucosal subsite but most frequently the tongue, followed by the floor of the mouth. The use of tobacco and betel nut, either smoked or chewed, and abuse of alcohol are the main risk factors for oral SCC. Oral SCC is characterized by considerable genetic heterogeneity and diversity, which together have a significant impact on the biological behaviour, clinical course, and response to treatment and on the generally poor prognosis of this carcinoma. Characterization of spatial and temporal tumour-specific molecular profiles and of person-specific resource availability and environmental and biological selective pressures could assist in personalizing anti-cancer treatment for individual patients, with the aim of improving treatment outcomes. In this narrative review, we discuss some of the events in cancer evolution and the functional significance of driver-mutations in carcinoma-related genes in general and elaborate on mechanisms mediating resistance to anti-cancer treatment.

Intracellular and extracellular factors of colorectal cancer liver metastasis: a pivotal perplex to be fully elucidated

Metastasis is the leading cause of death in colorectal cancer (CRC) patients, and the liver is the most common site of metastasis. Tumor cell metastasis can be thought of as an invasion-metastasis cascade and metastatic organotropism is thought to be a process that relies on the intrinsic properties of tumor cells and their interactions with molecules and cells in the microenvironment. Many studies have provided new insights into the molecular mechanism and contributing factors involved in CRC liver metastasis for a better understanding of the organ-specific metastasis process. The purpose of this review is to summarize the theories that explain CRC liver metastasis at multiple molecular dimensions (including genetic and non-genetic factors), as well as the main factors that cause CRC liver metastasis. Many findings suggest that metastasis may occur earlier than expected and with specific organ-anchoring property. The emergence of potential metastatic clones, the timing of dissemination, and the distinct routes of metastasis have been explained by genomic studies. The main force of CRC liver metastasis is also thought to be epigenetic alterations and dynamic phenotypic traits. Furthermore, we review key extrinsic factors that influence CRC cell metastasis and liver tropisms, such as pre-niches, tumor stromal cells, adhesion molecules, and immune/inflammatory responses in the tumor microenvironment. In addition, biomarkers associated with early diagnosis, prognosis, and recurrence of liver metastasis from CRC are summarized to enlighten potential clinical practice, including some markers that can be used as therapeutic targets to provide new perspectives for the treatment strategies of CRC liver metastasis.

Longitudinal Monitoring of Intra-Tumoural Heterogeneity Using Optical Barcoding of Patient-Derived Colorectal Tumour Models

Simple Summary

Colorectal cancer (CRC) is the second most common cancer worldwide. Despite improvements in the clinical management of CRC, outcomes for those with metastatic disease remain extremely poor. One reason for this is tumour heterogeneity, which refers to the observation that each cell within complex tumour cell populations displays different genetic features and biological behaviours. Such tumour heterogeneity is known to impact treatment efficacy and promote tumour recurrence. Here, we present a multi-colour barcoding methodology that allows for different lineages of colorectal cancer cells to be identified and monitored, thus allowing for tumour heterogeneity to be quantified in real-time. We show that discrete cell lineages can be quantified by both fluorescence microscopy and flow cytometry. Using this approach, we show that the cell culture models that are traditionally used in cancer research display limited heterogeneity, whereas patient-derived organoids—which are generated from fresh tumour resections—more faithfully represent the heterogeneity observed in cancer patients.

Abstract

Geno- and phenotypic heterogeneity amongst cancer cell subpopulations are established drivers of treatment resistance and tumour recurrence. However, due to the technical difficulty associated with studying such intra-tumoural heterogeneity, this phenomenon is seldom interrogated in conventional cell culture models. Here, we employ a fluorescent lineage technique termed “optical barcoding” (OBC) to perform simultaneous longitudinal tracking of spatio-temporal fate in 64 patient-derived colorectal cancer subclones. To do so, patient-derived cancer cell lines and organoids were labelled with discrete combinations of reporter constructs, stably integrated into the genome and thus passed on from the founder cell to all its clonal descendants. This strategy enables the longitudinal monitoring of individual cell lineages based upon their unique optical barcodes. By designing a novel panel of six fluorescent proteins, the maximum theoretical subpopulation resolution of 64 discriminable subpopulations was achieved, greatly improving throughput compared with previous studies. We demonstrate that all subpopulations can be purified from complex clonal mixtures via flow cytometry, permitting the downstream isolation and analysis of any lineages of interest. Moreover, we outline an optimized imaging protocol that can be used to image optical barcodes in real-time, allowing for clonal dynamics to be resolved in live cells. In contrast with the limited intra-tumour heterogeneity observed in conventional 2D cell lines, the OBC technique was successfully used to quantify dynamic clonal expansions and contractions in 3D patient-derived organoids, which were previously demonstrated to better recapitulate the heterogeneity of their parental tumour material. In summary, we present OBC as a user-friendly, inexpensive, and high-throughput technique for monitoring intra-tumoural heterogeneity in in vitro cell culture models.

The Genomic Processes of Biological Invasions: From Invasive Species to Cancer Metastases and Back Again

The concept of invasion is useful across a broad range of contexts, spanning from the fine scale landscape of cancer tumors up to the broader landscape of ecosystems. Invasion biology provides extraordinary opportunities for studying the mechanistic basis of contemporary evolution at the molecular level. Although the field of invasion genetics was established in ecology and evolution more than 50 years ago, there is still a limited understanding of how genomic level processes translate into invasive phenotypes across different taxa in response to complex environmental conditions. This is largely because the study of most invasive species is limited by information about complex genome level processes. We lack good reference genomes for most species. Rigorous studies to examine genomic processes are generally too costly. On the contrary, cancer studies are fortified with extensive resources for studying genome level dynamics and the interactions among genetic and non-genetic mechanisms. Extensive analysis of primary tumors and metastatic samples have revealed the importance of several genomic mechanisms including higher mutation rates, specific types of mutations, aneuploidy or whole genome doubling and non-genetic effects. Metastatic sites can be directly compared to primary tumor cell counterparts. At the same time, clonal dynamics shape the genomics and evolution of metastatic cancers. Clonal diversity varies by cancer type, and the tumors’ donor and recipient tissues. Still, the cancer research community has been unable to identify any common events that provide a universal predictor of “metastatic potential” which parallels findings in evolutionary ecology. Instead, invasion in cancer studies depends strongly on context, including order of events and clonal composition. The detailed studies of the behavior of a variety of human cancers promises to inform our understanding of genome level dynamics in the diversity of invasive species and provide novel insights for management.

Clonal Evolution of Multiple Myeloma-Clinical and Diagnostic Implications

Plasma cell dyscrasias are a heterogeneous group of diseases characterized by the expansion of bone marrow plasma cells. Malignant transformation of plasma cells depends on the continuity of events resulting in a sequence of well-defined disease stages, from monoclonal gammopathy of undetermined significance (MGUS) through smoldering myeloma (SMM) to symptomatic multiple myeloma (MM). Evolution of a pre-malignant cell into a malignant cell, as well as further tumor progression, dissemination, and relapse, require development of multiple driver lesions conferring selective advantage of the dominant clone and allowing subsequent evolution under selective pressure of microenvironment and treatment. This process of natural selection facilitates tumor plasticity leading to the formation of genetically complex and heterogenous tumors that are notoriously difficult to treat. Better understanding of the mechanisms underlying tumor evolution in MM and identification of lesions driving the evolution from the premalignant clone is therefore a key to development of effective treatment and long-term disease control. Here, we review recent advances in clonal evolution patterns and genomic landscape dynamics of MM, focusing on their clinical implications.

To portray clonal evolution in blood cancer, count your stem cells

Clonal evolution, the process of expansion and diversification of mutated cells, plays an important role in cancer development, resistance, and relapse. Although clonal evolution is most often conceived of as driven by natural selection, recent studies uncovered that neutral evolution shapes clonal evolution in a significant proportion of solid cancers. In hematological malignancies, the interplay between neutral evolution and natural selection is also disputed. Because natural selection selects cells with a greater fitness, providing a growth advantage to some cells relative to others, the architecture of clonal evolution serves as indirect evidence to distinguish natural selection from neutral evolution and has been associated with different prognoses for the patient. Linear architecture, when the new mutant clone grows within the previous one, is characteristic of hematological malignancies and is typically interpreted as being driven by natural selection. Here, we discuss the role of natural selection and neutral evolution in the production of linear clonal architectures in hematological malignancies. Although it is tempting to attribute linear evolution to natural selection, we argue that a lower number of contributing stem cells accompanied by genetic drift can also result in a linear pattern of evolution, as illustrated by simulations of clonal evolution in hematopoietic stem cells. The number of stem cells contributing to long-term clonal evolution is not known in the pathological context, and we advocate that estimating these numbers in the context of cancer and aging is crucial to parsing out neutral evolution from natural selection, 2 processes that require different therapeutic strategies.

Genomic landscape and clonal architecture of mouse oral squamous cell carcinomas dictate tumour ecology

To establish whether 4-nitroquinoline N-oxide-induced carcinogenesis mirrors the heterogeneity of human oral squamous cell carcinoma (OSCC), we have performed genomic analysis of mouse tongue lesions. The mutational signatures of human and mouse OSCC overlap extensively. Mutational burden is higher in moderate dysplasias and invasive SCCs than in hyperplasias and mild dysplasias, although mutations in p53, Notch1 and Fat1 occur in early lesions. Laminin-α3 mutations are associated with tumour invasiveness and Notch1 mutant tumours have an increased immune infiltrate. Computational modelling of clonal dynamics indicates that high genetic heterogeneity may be a feature of those mild dysplasias that are likely to progress to more aggressive tumours. These studies provide a foundation for exploring OSCC evolution, heterogeneity and progression.

Genetic Alterations of Metastatic Colorectal Cancer

Genome sequencing studies have characterized the genetic alterations of different tumor types, highlighting the diversity of the molecular processes driving tumor development. Comprehensive sequencing studies have defined molecular subtypes of colorectal cancers (CRCs) through the identification of genetic events associated with microsatellite stability (MSS), microsatellite-instability-high (MSI-H), and hypermutation. Most of these studies characterized primary tumors. Only recent studies have addressed the characterization of the genetic and clinical heterogeneity of metastatic CRC. Metastatic CRC genomes were found to be not fundamentally different from primary CRCs in terms of the mutational landscape or of genes that drive tumorigenesis, and a genomic heterogeneity associated with tumor location of primary tumors helps to define different clinical behaviors of metastatic CRCs. Although CRC metastatic spreading was traditionally seen as a late-occurring event, growing evidence suggests that this process can begin early during tumor development and the clonal architecture of these tumors is consistently influenced by cancer treatment. Although the survival rate of patients with metastatic CRC patients improved in the last years, the response to current treatments and prognosis of many of these patients remain still poor, indicating the need to discover new improvements for therapeutic vulnerabilities and to formulate a rational prospective of personalized therapies.

Intratumor heterogeneity: A new perspective on colorectal cancer research

Colorectal cancers generally consist of multiple subclones. These subclones have their own unique characteristics, resulting in intratumor heterogeneity (ITH). As the discussion of ITH has advanced, a model describing the relationship of ITH to the tumor has gradually emerged. ITH can be divided into two types of intraprimary tumor heterogeneity and intraindividual tumor heterogeneity, the former for further understanding of tumor composition, and the latter for providing more information about evolutionary patterns. With the rapid development of new methods, such as next‐generation, polyguanine region sequencing, and Image detection, researchers may unravel the secrets underlying ITH. The higher the ITH of the tumor, the richer the interaction between the subclones maybe, or the greater the chance of the tumor getting more powerful subclones may be, thus increasing the malignant potential of the tumor. Existing evidence suggests that ITH may increase the ability of tumors to resist treatment and can be used as an independent influence on the prognosis of colorectal cancer. We reviewed 80 recent studies to give researchers a new perspective on colorectal cancer. There is still a limited amount of research in this area. Further study of the relationship between ITH and clinical endpoints may lead to the development of new treatment strategies.

The Intratumoral Heterogeneity Reflects the Intertumoral Subtypes of Glioblastoma Multiforme: A Regional Immunohistochemistry Analysis

Glioblastoma multiforme (GBM) is the most frequent and aggressive primary brain tumor in adults. Despite extensive therapy the prognosis for GBM patients remains poor and the extraordinary therapy resistance has been attributed to intertumoral heterogeneity of glioblastoma. Different prognostic relevant GBM tumor subtypes have been identified based on their molecular profile. This approach, however, neglects the heterogeneity within individual tumors, that is, the intratumoral heterogeneity. Here, we detected the regional immunoreactivity by immunohistochemistry and immunofluorescence using nine different markers on resected GBM specimens (IDH wildtype, WHO grade IV). We found repetitive expression profiles, that could be classified into clusters. These clusters could then be assigned to five pathophysiologically relevant groups that reflect the previously described subclasses of GBM, including mesenchymal, classical, and proneural subtype. Our data indicate the presence of tumor differentiations and tumor subclasses that occur within individual tumors, and might therefore contribute to develop adapted, individual-based therapies.

Tumour heterogeneity and the evolutionary trade-offs of cancer

Tumours vary in gene expression programmes and genetic alterations. Understanding this diversity and its biological meaning requires a theoretical framework, which could in turn guide the development of more accurate prognosis and therapy. Here, we review the theory of multi-task evolution of cancer, which is based upon the premise that tumours evolve in the host and face selection trade-offs between multiple biological functions. This theory can help identify the major biological tasks that cancer cells perform and the trade-offs between these tasks. It introduces the concept of specialist tumours, which focus on one task, and generalist tumours, which perform several tasks. Specialist tumours are suggested to be sensitive to therapy targeting their main task. Driver mutations tune gene expression towards specific tasks in a tissue-dependent manner and thus help to determine whether a tumour is specialist or generalist. We discuss potential applications of the theory of multi-task evolution to interpret the spatial organization of tumours and intratumour heterogeneity.

Assessing reliability of intra-tumor heterogeneity estimates from single sample whole exome sequencing data

Tumors are made of evolving and heterogeneous populations of cells which arise from successive appearance and expansion of subclonal populations, following acquisition of mutations conferring them a selective advantage. Those subclonal populations can be sensitive or resistant to different treatments, and provide information about tumor aetiology and future evolution. Hence, it is important to be able to assess the level of heterogeneity of tumors with high reliability for clinical applications. In the past few years, a large number of methods have been proposed to estimate intra-tumor heterogeneity from whole exome sequencing (WES) data, but the accuracy and robustness of these methods on real data remains elusive. Here we systematically apply and compare 6 computational methods to estimate tumor heterogeneity on 1,697 WES samples from the cancer genome atlas (TCGA) covering 3 cancer types (breast invasive carcinoma, bladder urothelial carcinoma, and head and neck squamous cell carcinoma), and two distinct input mutation sets. We observe significant differences between the estimates produced by different methods, and identify several likely confounding factors in heterogeneity assessment for the different methods. We further show that the prognostic value of tumor heterogeneity for survival prediction is limited in those datasets, and find no evidence that it improves over prognosis based on other clinical variables. In conclusion, heterogeneity inference from WES data on a single sample, and its use in cancer prognosis, should be considered with caution. Other approaches to assess intra-tumoral heterogeneity such as those based on multiple samples may be preferable for clinical applications.

An analysis of genetic heterogeneity in untreated cancers

Genetic intratumoural heterogeneity is a natural consequence of imperfect DNA replication. Any two randomly selected cells, whether normal or cancerous, are therefore genetically different. Here, we review the different forms of genetic heterogeneity in cancer and re-analyse the extent of genetic heterogeneity within seven types of untreated epithelial cancers, with particular regard to its clinical relevance. We find that the homogeneity of predicted functional mutations in driver genes is the rule rather than the exception. In primary tumours with multiple samples, 97% of driver-gene mutations in 38 patients were homogeneous. Moreover, among metastases from the same primary tumour, 100% of the driver mutations in 17 patients were homogeneous. With a single biopsy of a primary tumour in 14 patients, the likelihood of missing a functional driver-gene mutation that was present in all metastases was 2.6%. Furthermore, all functional driver-gene mutations detected in these 14 primary tumours were present among all their metastases. Finally, we found that individual metastatic lesions responded concordantly to targeted therapies in 91% of 44 patients. These analyses indicate that the cells within the primary tumours that gave rise to metastases are genetically homogeneous with respect to functional driver-gene mutations, and we suggest that future efforts to develop combination therapies have the potential to be curative.

Big Bang Tumor Growth and Clonal Evolution

The advent and application of next-generation sequencing (NGS) technologies to tumor genomes has reinvigorated efforts to understand clonal evolution. Although tumor progression has traditionally been viewed as a gradual stepwise process, recent studies suggest that evolutionary rates in tumors can be variable with periods of punctuated mutational bursts and relative stasis. For example, Big Bang dynamics have been reported, wherein after transformation, growth occurs in the absence of stringent selection, consistent with effectively neutral evolution. Although first noted in colorectal tumors, effective neutrality may be relatively common. Additionally, punctuated evolution resulting from mutational bursts and cataclysmic genomic alterations have been described. In this review, we contrast these findings with the conventional gradualist view of clonal evolution and describe potential clinical and therapeutic implications of different evolutionary modes and tempos.