Macrophages at the interface of the co-evolving cancer ecosystem.

Macrophages are versatile and heterogeneous innate immune cells undertaking central functions in balancing immune responses and tissue repair to maintain homeostasis. This plasticity, once co-opted by malignant outgrowth, orchestrates manifold reciprocal interactions within the tumor microenvironment, fueling the evolution of the cancer ecosystem. Here, we review the multilayered sources of influence that jointly underpin and longitudinally shape tumor-associated macrophage (TAM) phenotypic states in solid neoplasms. We discuss how, in response to these signals, TAMs steer tumor evolution in the context of natural selection, biological dispersion, and treatment resistance. A number of research frontiers to be tackled are laid down in this review to therapeutically exploit the complex roles of TAMs in cancer. Building upon knowledge obtained from currently applied TAM-targeting strategies and using next generation technologies, we propose conceptual advances and novel therapeutic avenues to rewire TAM multifaceted regulation of the co-evolving cancer ecosystem.

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

Distinct Classes of Complex Structural Variation Uncovered across Thousands of Cancer Genome Graphs.

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.

Timing the Landmark Events in the Evolution of Clear Cell Renal Cell Cancer: TRACERx Renal.

Clear cell renal cell carcinoma (ccRCC) is characterized by near-universal loss of the short arm of chromosome 3, deleting several tumor suppressor genes. We analyzed whole genomes from 95 biopsies across 33 patients with clear cell renal cell carcinoma. We find hotspots of point mutations in the 5' UTR of TERT, targeting a MYC-MAX-MAD1 repressor associated with telomere lengthening. The most common structural abnormality generates simultaneous 3p loss and 5q gain (36% patients), typically through chromothripsis. This event occurs in childhood or adolescence, generally as the initiating event that precedes emergence of the tumor's most recent common ancestor by years to decades. Similar genomic changes drive inherited ccRCC. Modeling differences in age incidence between inherited and sporadic cancers suggests that the number of cells with 3p loss capable of initiating sporadic tumors is no more than a few hundred. Early development of ccRCC follows well-defined evolutionary trajectories, offering opportunity for early intervention.

Deterministic Evolutionary Trajectories Influence Primary Tumor Growth: TRACERx Renal.

The evolutionary features of clear-cell renal cell carcinoma (ccRCC) have not been systematically studied to date. We analyzed 1,206 primary tumor regions from 101 patients recruited into the multi-center prospective study, TRACERx Renal. We observe up to 30 driver events per tumor and show that subclonal diversification is associated with known prognostic parameters. By resolving the patterns of driver event ordering, co-occurrence, and mutual exclusivity at clone level, we show the deterministic nature of clonal evolution. ccRCC can be grouped into seven evolutionary subtypes, ranging from tumors characterized by early fixation of multiple mutational and copy number drivers and rapid metastases to highly branched tumors with >10 subclonal drivers and extensive parallel evolution associated with attenuated progression. We identify genetic diversity and chromosomal complexity as determinants of patient outcome. Our insights reconcile the variable clinical behavior of ccRCC and suggest evolutionary potential as a biomarker for both intervention and surveillance.

Allele-Specific HLA Loss and Immune Escape in Lung Cancer Evolution.

Immune evasion is a hallmark of cancer. Losing the ability to present neoantigens through human leukocyte antigen (HLA) loss may facilitate immune evasion. However, the polymorphic nature of the locus has precluded accurate HLA copy-number analysis. Here, we present loss of heterozygosity in human leukocyte antigen (LOHHLA), a computational tool to determine HLA allele-specific copy number from sequencing data. Using LOHHLA, we find that HLA LOH occurs in 40% of non-small-cell lung cancers (NSCLCs) and is associated with a high subclonal neoantigen burden, APOBEC-mediated mutagenesis, upregulation of cytolytic activity, and PD-L1 positivity. The focal nature of HLA LOH alterations, their subclonal frequencies, enrichment in metastatic sites, and occurrence as parallel events suggests that HLA LOH is an immune escape mechanism that is subject to strong microenvironmental selection pressures later in tumor evolution. Characterizing HLA LOH with LOHHLA refines neoantigen prediction and may have implications for our understanding of resistance mechanisms and immunotherapeutic approaches targeting neoantigens. VIDEO ABSTRACT.

Boveri and beyond: Chromothripsis and genomic instability from mitotic errors.

Mitotic cell division is tightly monitored by checkpoints that safeguard the genome from instability. Failures in accurate chromosome segregation during mitosis can cause numerical aneuploidy, which was hypothesized by Theodor Boveri over a century ago to promote tumorigenesis. Recent interrogation of pan-cancer genomes has identified unexpected classes of chromosomal abnormalities, including complex rearrangements arising through chromothripsis. This process is driven by mitotic errors that generate abnormal nuclear structures that provoke extensive yet localized shattering of mis-segregated chromosomes. Here, we discuss emerging mechanisms underlying chromothripsis from micronuclei and chromatin bridges, as well as highlight how this mutational cascade converges on the DNA damage response. A fundamental understanding of these catastrophic processes will provide insight into how initial errors in mitosis can precipitate rapid cancer genome evolution.

Clonal selection parallels between normal and cancer tissues.

Clonal selection and drift drive both normal tissue and cancer development. However, the biological mechanisms and environmental conditions underpinning these processes remain to be elucidated. Clonal selection models are centered in Darwinian evolutionary theory, where some clones with the fittest features are selected and populate the tissue or tumor. We suggest that different subclasses of stem cells, each of which is responsible for a distinct feature of the selection process, share common features between normal and cancer conditions. While active stem cells populate the tissue, dormant cells account for tissue replenishment/regeneration in both normal and cancerous tissues. We also discuss potential mechanisms that drive clonal drift, their interactions with clonal selection, and their similarities during normal and cancer tissue development.

Immunosurveillance encounters cancer metabolism.

Tumor cells reprogram nutrient acquisition and metabolic pathways to meet their energetic, biosynthetic, and redox demands. Similarly, metabolic processes in immune cells support host immunity against cancer and determine differentiation and fate of leukocytes. Thus, metabolic deregulation and imbalance in immune cells within the tumor microenvironment have been reported to drive immune evasion and to compromise therapeutic outcomes. Interestingly, emerging evidence indicates that anti-tumor immunity could modulate tumor heterogeneity, aggressiveness, and metabolic reprogramming, suggesting that immunosurveillance can instruct cancer progression in multiple dimensions. This review summarizes our current understanding of how metabolic crosstalk within tumors affects immunogenicity of tumor cells and promotes cancer progression. Furthermore, we explain how defects in the metabolic cascade can contribute to developing dysfunctional immune responses against cancers and discuss the contribution of immunosurveillance to these defects as a feedback mechanism. Finally, we highlight ongoing clinical trials and new therapeutic strategies targeting cellular metabolism in cancer.

Fully accessible fitness landscape of oncogene-negative lung adenocarcinoma.

Cancer genomes are almost invariably complex with genomic alterations cooperating during each step of carcinogenesis. In cancers that lack a single dominant oncogene mutation, cooperation between the inactivation of multiple tumor suppressor genes can drive tumor initiation and growth. Here, we shed light on how the sequential acquisition of genomic alterations generates oncogene-negative lung tumors. We couple tumor barcoding with combinatorial and multiplexed somatic genome editing to characterize the fitness landscapes of three tumor suppressor genes NF1, RASA1, and PTEN, the inactivation of which jointly drives oncogene-negative lung adenocarcinoma initiation and growth. The fitness landscape was surprisingly accessible, with each additional mutation leading to growth advantage. Furthermore, the fitness landscapes remained fully accessible across backgrounds with the inactivation of additional tumor suppressor genes. These results suggest that while predicting cancer evolution will be challenging, acquiring the multiple alterations that drive the growth of oncogene-negative tumors can be facilitated by the lack of constraints on mutational order.

A fitness trade-off explains the early fate of yeast aneuploids with chromosome gains.

The early development of aneuploidy from an accidental chromosome missegregation shows contrasting effects. On the one hand, it is associated with significant cellular stress and decreased fitness. On the other hand, it often carries a beneficial effect and provides a quick (but typically transient) solution to external stress. These apparently controversial trends emerge in several experimental contexts, particularly in the presence of duplicated chromosomes. However, we lack a mathematical evolutionary modeling framework that comprehensively captures these trends from the mutational dynamics and the trade-offs involved in the early stages of aneuploidy. Here, focusing on chromosome gains, we address this point by introducing a fitness model where a fitness cost of chromosome duplications is contrasted by a fitness advantage from the dosage of specific genes. The model successfully captures the experimentally measured probability of emergence of extra chromosomes in a laboratory evolution setup. Additionally, using phenotypic data collected in rich media, we explored the fitness landscape, finding evidence supporting the existence of a per-gene cost of extra chromosomes. Finally, we show that the substitution dynamics of our model, evaluated in the empirical fitness landscape, explains the relative abundance of duplicated chromosomes observed in yeast population genomics data. These findings lay a firm framework for the understanding of the establishment of newly duplicated chromosomes, providing testable quantitative predictions for future observations.

Evolutionary impact of chimeric RNAs on generating phenotypic plasticity in human cells.

Chimeric RNAs are generated by the fusion of the exons or introns of two genes. The generation of chimeric RNAs is important for the functional expansion of cells. Here, we describe the functional implications of chimeric RNAs for generating phenotypic plasticity from an evolutionary perspective.

Cancer evolution: Darwin and beyond.

Clinical and laboratory studies over recent decades have established branched evolution as a feature of cancer. However, while grounded in somatic selection, several lines of evidence suggest a Darwinian model alone is insufficient to fully explain cancer evolution. First, the role of macroevolutionary events in tumour initiation and progression contradicts Darwin's central thesis of gradualism. Whole-genome doubling, chromosomal chromoplexy and chromothripsis represent examples of single catastrophic events which can drive tumour evolution. Second, neutral evolution can play a role in some tumours, indicating that selection is not always driving evolution. Third, increasing appreciation of the role of the ageing soma has led to recent generalised theories of age-dependent carcinogenesis. Here, we review these concepts and others, which collectively argue for a model of cancer evolution which extends beyond Darwin. We also highlight clinical opportunities which can be grasped through targeting cancer vulnerabilities arising from non-Darwinian patterns of evolution.

Genetic and immune landscape evolution in MMR-deficient colorectal cancer.

Mismatch repair-deficient (MMRd) colorectal cancers (CRCs) have high mutation burdens, which make these tumours immunogenic and many respond to immune checkpoint inhibitors. The MMRd hypermutator phenotype may also promote intratumour heterogeneity (ITH) and cancer evolution. We applied multiregion sequencing and CD8 and programmed death ligand 1 (PD-L1) immunostaining to systematically investigate ITH and how genetic and immune landscapes coevolve. All cases had high truncal mutation burdens. Despite pervasive ITH, driver aberrations showed a clear hierarchy. Those in WNT/ß-catenin, mitogen-activated protein kinase, and TGF-ß receptor family genes were almost always truncal. Immune evasion (IE) drivers, such as inactivation of genes involved in antigen presentation or IFN-γ signalling, were predominantly subclonal and showed parallel evolution. These IE drivers have been implicated in immune checkpoint inhibitor resistance or sensitivity. Clonality assessments are therefore important for the development of predictive immunotherapy biomarkers in MMRd CRCs. Phylogenetic analysis identified three distinct patterns of IE driver evolution: pan-tumour evolution, subclonal evolution, and evolutionary stasis. These, but neither mutation burdens nor heterogeneity metrics, significantly correlated with T-cell densities, which were used as a surrogate marker of tumour immunogenicity. Furthermore, this revealed that genetic and T-cell infiltrates coevolve in MMRd CRCs. Low T-cell densities in the subgroup without any known IE drivers may indicate an, as yet unknown, IE mechanism. PD-L1 was expressed in the tumour microenvironment in most samples and correlated with T-cell densities. However, PD-L1 expression in cancer cells was independent of T-cell densities but strongly associated with loss of the intestinal homeobox transcription factor CDX2. This explains infrequent PD-L1 expression by cancer cells and may contribute to a higher recurrence risk of MMRd CRCs with impaired CDX2 expression. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Wild animals as an underused treasure trove for studying the genetics of cancer.

Recent years have seen an emergence of the field of comparative cancer genomics. However, the advancements in this field are held back by the hesitation to use knowledge obtained from human studies to study cancer in other animals, and vice versa. Since cancer is an ancient disease that arose with multicellularity, oncogenes and tumour-suppressor genes are amongst the oldest gene classes, shared by most animal species. Acknowledging that other animals are, in terms of cancer genetics, ecology, and evolution, rather similar to humans, creates huge potential for advancing the fields of human and animal oncology, but also biodiversity conservation. Also see the video abstract here: https://youtu.be/UFqyMx5HETY.

Clonal interactions in cancer: Integrating quantitative models with experimental and clinical data.

Tumors consist of different genotypically distinct subpopulations-or subclones-of cells. These subclones can influence neighboring clones in a process called "clonal interaction." Conventionally, research on driver mutations in cancer has focused on their cell-autonomous effects that lead to an increase in fitness of the cells containing the driver. Recently, with the advent of improved experimental and computational technologies for investigating tumor heterogeneity and clonal dynamics, new studies have shown the importance of clonal interactions in cancer initiation, progression, and metastasis. In this review we provide an overview of clonal interactions in cancer, discussing key discoveries from a diverse range of approaches to cancer biology research. We discuss common types of clonal interactions, such as cooperation and competition, its mechanisms, and the overall effect on tumorigenesis, with important implications for tumor heterogeneity, resistance to treatment, and tumor suppression. Quantitative models-in coordination with cell culture and animal model experiments-have played a vital role in investigating the nature of clonal interactions and the complex clonal dynamics they generate. We present mathematical and computational models that can be used to represent clonal interactions and provide examples of the roles they have played in identifying and quantifying the strength of clonal interactions in experimental systems. Clonal interactions have proved difficult to observe in clinical data; however, several very recent quantitative approaches enable their detection. We conclude by discussing ways in which researchers can further integrate quantitative methods with experimental and clinical data to elucidate the critical-and often surprising-roles of clonal interactions in human cancers.

What can we learn about acid-base transporters in cancer from studying somatic mutations in their genes?

Acidosis is a chemical signature of the tumour microenvironment that challenges intracellular pH homeostasis. The orchestrated activity of acid-base transporters of the solute-linked carrier (SLC) family is critical for removing the end-products of fermentative metabolism (lactate/H+) and maintaining a favourably alkaline cytoplasm. Given the critical role of pH homeostasis in enabling cellular activities, mutations in relevant SLC genes may impact the oncogenic process, emerging as negatively or positively selected, or as driver or passenger mutations. To address this, we performed a pan-cancer analysis of The Cancer Genome Atlas simple nucleotide variation data for acid/base-transporting SLCs (ABT-SLCs). Somatic mutation patterns of monocarboxylate transporters (MCTs) were consistent with their proposed essentiality in facilitating lactate/H+ efflux. Among all cancers, tumours of uterine corpus endometrial cancer carried more ABT-SLC somatic mutations than expected from median tumour mutation burden. Among these, somatic mutations in SLC4A3 had features consistent with meaningful consequences on cellular fitness. Definitive evidence for ABT-SLCs as 'cancer essential' or 'driver genes' will have to consider microenvironmental context in genomic sequencing because bulk approaches are insensitive to pH heterogeneity within tumours. Moreover, genomic analyses must be validated with phenotypic outcomes (i.e. SLC-carried flux) to appreciate the opportunities for targeting acid-base transport in cancers.

The Evolutionary Landscape of Colorectal Tumorigenesis: Recent Paradigms, Models, and Hypotheses.

Using colorectal cancer as a model, we review some of the insights into cancer evolution afforded by cancer sequencing. These include nonlinear and neutral evolution; polyclonality of driver mutations and parallel evolution in adenomas, although these are rare in carcinomas; the ability of mutational processes to shape evolution against the force of selection; the presence of rare driver genes that function in the same signaling pathways as the longstanding canonical drivers; and the existence of selective windows that constrain the functional effects of cancer driver mutations within limits. Many of these nascent evolutionary paradigms are potentially important for treating colorectal cancers as well as understanding their development.

The Ecology and Evolution of Cancer: The Ultra-Microevolutionary Process.

Although tumorigenesis has been accepted as an evolutionary process ( 20 , 102 ), many forces may operate differently in cancers than in organisms, as they evolve at vastly different time scales. Among such forces, natural selection, here defined as differential cellular proliferation among distinct somatic cell genotypes, is particularly interesting because its action might be thwarted in multicellular organisms ( 20 , 29 ). In this review, selection is analyzed in two stages of cancer evolution: Stage I is the evolution between tumors and normal tissues, and Stage II is the evolution within tumors. The Cancer Genome Atlas (TCGA) data show a low degree of convergent evolution in Stage I, where genetic changes are not extensively shared among cases. An equally important, albeit much less highlighted, discovery using TCGA data is that there is almost no net selection in cancer evolution. Both positive and negative selection are evident but they neatly cancel each other out, rendering total selection ineffective in the absence of recombination. The efficacy of selection is even lower in Stage II, where neutral (non-Darwinian) evolution is increasingly supported by high-density sampling studies ( 81 , 123 ). Because natural selection is not a strong deterministic force, cancers usually evolve divergently even in similar tissue environments.

BiTSC 2: Bayesian inference of tumor clonal tree by joint analysis of single-cell SNV and CNA data.

The rapid development of single-cell DNA sequencing (scDNA-seq) technology has greatly enhanced the resolution of tumor cell profiling, providing an unprecedented perspective in characterizing intra-tumoral heterogeneity and understanding tumor progression and metastasis. However, prominent algorithms for constructing tumor phylogeny based on scDNA-seq data usually only take single nucleotide variations (SNVs) as markers, failing to consider the effect caused by copy number alterations (CNAs). Here, we propose BiTSC$^2$, Bayesian inference of Tumor clonal Tree by joint analysis of Single-Cell SNV and CNA data. BiTSC$^2$ takes raw reads from scDNA-seq as input, accounts for the overlapping of CNA and SNV, models allelic dropout rate, sequencing errors and missing rate, as well as assigns single cells into subclones. By applying Markov Chain Monte Carlo sampling, BiTSC$^2$ can simultaneously estimate the subclonal scCNA and scSNV genotype matrices, subclonal assignments and tumor subclonal evolutionary tree. In comparison with existing methods on synthetic and real tumor data, BiTSC$^2$ shows high accuracy in genotype recovery, subclonal assignment and tree reconstruction. BiTSC$^2$ also performs robustly in dealing with scDNA-seq data with low sequencing depth and variant missing rate. BiTSC$^2$ software is available at https://github.com/ucasdp/BiTSC2.