Negative selection in tumor genome evolution acts on essential cellular functions and the immunopeptidome

Copy-number alterations in cell-free DNA can be transient or harbingers of clonal hematopoiesis

Genome-wide plasma cfDNA pan-cancer screening of 1002 healthy elderly identified 15 individuals with CNAs of unknown origin. Nine participants were reassessed over 3-5 years through health questionnaires, WB-MRI, and cfDNA and blood analyses. CNAs resolved in two cases but persisted in seven mainly associated with low-grade clonal mosaicism. These findings suggest cfDNA CNAs may be transient or serve as early markers of clonal mosaicism, preceding clinical detection by years.

Comprehensive assessment of computational methods for cancer immunoediting

Cancer immunoediting reflects the role of the immune system in eliminating tumor cells and shaping tumor immunogenicity, which leaves marks in the genome. In this study, we systematically evaluate four methods for quantifying immunoediting. In colorectal cancer samples from The Cancer Genome Atlas, we found that these methods identified 78.41%, 46.17%, 36.61%, and 4.92% of immunoedited samples, respectively, covering 92.90% of all colorectal cancer samples. Comparison of 10 patient-derived xenografts (PDXs) with their original tumors showed that different methods identified reduced immune selection in PDXs ranging from 44.44% to 60.0%. The proportion of such PDX-tumor pairs increases to 77.78% when considering the union of results from multiple methods, indicating the complementarity of these methods. We find that observed-to-expected ratios highly rely on neoantigen selection criteria and reference datasets. In contrast, HLA-binding mutation ratio, immune dN/dS, and enrichment score of cancer cell fraction were less affected by these factors. Our findings suggest integration of multiple methods may benefit future immunoediting analyses.

The theory of massively repeated evolution and full identifications of cancer-driving nucleotides (CDNs)

Tumorigenesis, like most complex genetic traits, is driven by the joint actions of many mutations. At the nucleotide level, such mutations are cancer-driving nucleotides (CDNs). The full sets of CDNs are necessary, and perhaps even sufficient, for the understanding and treatment of each cancer patient. Currently, only a small fraction of CDNs is known as most mutations accrued in tumors are not drivers. We now develop the theory of CDNs on the basis that cancer evolution is massively repeated in millions of individuals. Hence, any advantageous mutation should recur frequently and, conversely, any mutation that does not is either a passenger or deleterious mutation. In the TCGA cancer database (sample size n=300-1000), point mutations may recur in i out of n patients. This study explores a wide range of mutation characteristics to determine the limit of recurrences (i*) driven solely by neutral evolution. Since no neutral mutation can reach i*=3, all mutations recurring at i≥3 are CDNs. The theory shows the feasibility of identifying almost all CDNs if n increases to 100,000 for each cancer type. At present, only <10% of CDNs have been identified. When the full sets of CDNs are identified, the evolutionary mechanism of tumorigenesis in each case can be known and, importantly, gene targeted therapy will be far more effective in treatment and robust against drug resistance.

Cancers adapt to their mutational load by buffering protein misfolding stress

In asexual populations that don't undergo recombination, such as cancer, deleterious mutations are expected to accrue readily due to genome-wide linkage between mutations. Despite this mutational load of often thousands of deleterious mutations, many tumors thrive. How tumors survive the damaging consequences of this mutational load is not well understood. Here, we investigate the functional consequences of mutational load in 10,295 human tumors by quantifying their phenotypic response through changes in gene expression. Using a generalized linear mixed model (GLMM), we find that high mutational load tumors up-regulate proteostasis machinery related to the mitigation and prevention of protein misfolding. We replicate these expression responses in cancer cell lines and show that the viability in high mutational load cancer cells is strongly dependent on complexes that degrade and refold proteins. This indicates that the upregulation of proteostasis machinery is causally important for high mutational burden tumors and uncovers new therapeutic vulnerabilities.

Evolutionary dependency of cancer mutations in gene pairs inferred by nonsynonymous-synonymous mutation ratios

BACKGROUND

Determining the impact of somatic mutations requires understanding the functional relationship of genes acquiring mutations; however, it is largely unknown how mutations in functionally related genes influence each other.

METHODS

We employed non-synonymous-to-synonymous or dNdS ratios to evaluate the evolutionary dependency (ED) of gene pairs, assuming a mutation in one gene of a gene pair can affect the evolutionary fitness of mutations in its partner genes as mutation context. We employed PanCancer- and tumor type-specific mutational profiles to infer the ED of gene pairs and evaluated their biological relevance with respect to gene dependency and drug sensitivity.

RESULTS

We propose that dNdS ratios of gene pairs and their derived cdNS (context-dependent dNdS) scores as measure of ED distinguishing gene pairs either as synergistic (SYN) or antagonistic (ANT). Mutation contexts can induce substantial changes in the evolutionary fitness of mutations in the paired genes, e.g., IDH1 and IDH2 mutation contexts lead to substantial increase and decrease of dNdS ratios of ATRX indels and IDH1 missense mutations corresponding to SYN and ANT relationship with positive and negative cdNS scores, respectively. The impact of gene silencing or knock-outs on cell viability (genetic dependencies) often depends on ED, suggesting that ED can guide the selection of candidates for synthetic lethality such as TCF7L2-KRAS mutations. Using cell line-based drug sensitivity data, the effects of targeted agents on cell lines are often associated with mutations of genes exhibiting ED with the target genes, informing drug sensitizing or resistant mutations for targeted inhibitors, e.g., PRSS1 and CTCF mutations as resistant mutations to EGFR and BRAF inhibitors for lung adenocarcinomas and melanomas, respectively.

CONCLUSIONS

We propose that the ED of gene pairs evaluated by dNdS ratios can advance our understanding of the functional relationship of genes with potential biological and clinical implications.

Copy number losses of oncogenes and gains of tumor suppressor genes generate common driver mutations

Cancer driver genes can undergo positive selection for various types of genetic alterations, including gain-of-function or loss-of-function mutations and copy number alterations (CNA). We investigated the landscape of different types of alterations affecting driver genes in 17,644 cancer exomes and genomes. We find that oncogenes may simultaneously exhibit signatures of positive selection and also negative selection in different gene segments, suggesting a method to identify additional tumor types where an oncogene is a driver or a vulnerability. Next, we characterize the landscape of CNA-dependent selection effects, revealing a general trend of increased positive selection on oncogene mutations not only upon CNA gains but also upon CNA deletions. Similarly, we observe a positive interaction between mutations and CNA gains in tumor suppressor genes. Thus, two-hit events involving point mutations and CNA are universally observed regardless of the type of CNA and may signal new therapeutic opportunities. An analysis with focus on the somatic CNA two-hit events can help identify additional driver genes relevant to a tumor type. By a global inference of point mutation and CNA selection signatures and interactions thereof across genes and tissues, we identify 9 evolutionary archetypes of driver genes, representing different mechanisms of (in)activation by genetic alterations.

Neoadjuvant sintilimab plus chemotherapy in EGFR-mutant NSCLC: Phase 2 trial interim results (NEOTIDE/CTONG2104)

The clinical efficacy of neoadjuvant immunotherapy plus chemotherapy remains elusive in localized epidermal growth factor receptor (EGFR)-mutant non-small cell lung cancer (NSCLC). Here, we report interim results of a Simon's two-stage design, phase 2 trial using neoadjuvant sintilimab with carboplatin and nab-paclitaxel in resectable EGFR-mutant NSCLC. All 18 patients undergo radical surgery, with one patient experiencing surgery delay. Fourteen patients exhibit confirmed radiological response, with 44% achieving major pathological response (MPR) and no pathological complete response (pCR). Similar genomic alterations are observed before and after treatment without influencing the efficacy of subsequent EGFR-tyrosine kinase inhibitors (TKIs) in vitro. Infiltration and T cell receptor (TCR) clonal expansion of CCR8+ regulatory T (Treg)hi/CXCL13+ exhausted T (Tex)lo cells define a subtype of EGFR-mutant NSCLC highly resistant to immunotherapy, with the phenotype potentially serving as a promising signature to predict immunotherapy efficacy. Informed circulating tumor DNA (ctDNA) detection in EGFR-mutant NSCLC could help identify patients nonresponsive to neoadjuvant immunochemotherapy. These findings provide supportive data for the utilization of neoadjuvant immunochemotherapy and insight into immune resistance in EGFR-mutant NSCLC.

Immune evasion impacts the landscape of driver genes during cancer evolution

BACKGROUND

Carcinogenesis is driven by interactions between genetic mutations and the local tumor microenvironment. Recent research has identified hundreds of cancer driver genes; however, these studies often include a mixture of different molecular subtypes and ecological niches and ignore the impact of the immune system.

RESULTS

In this study, we compare the landscape of driver genes in tumors that escaped the immune system (escape +) versus those that did not (escape -). We analyze 9896 primary tumors from The Cancer Genome Atlas using the ratio of non-synonymous to synonymous mutations (dN/dS) and find 85 driver genes, including 27 and 16 novel genes, in escape - and escape + tumors, respectively. The dN/dS of driver genes in immune escaped tumors is significantly lower and closer to neutrality than in non-escaped tumors, suggesting selection buffering in driver genes fueled by immune escape. Additionally, we find that immune evasion leads to more mutated sites, a diverse array of mutational signatures and is linked to tumor prognosis.

CONCLUSIONS

Our findings highlight the need for improved patient stratification to identify new therapeutic targets for cancer treatment.

Frequency-dependent selection of neoantigens fosters tumor immune escape and predicts immunotherapy response

Cancer is an evolutionary process shaped by selective pressure from the microenvironments. However, recent studies reveal that certain tumors undergo neutral evolution where there is no detectable fitness difference amongst the cells following malignant transformation. Here, through computational modeling, we demonstrate that negative frequency-dependent selection (or NFDS), where the immune response against cancer cells depends on the clonality of neoantigens, can lead to an immunogenic landscape that is highly similar to neutral evolution. Crucially, NFDS promotes high antigenic heterogeneity and early immune evasion in hypermutable tumors, leading to poor responses to immune checkpoint blockade (ICB) therapy. Our model also reveals that NFDS is characterized by a negative association between average clonality and total burden of neoantigens. Indeed, this unique feature of NFDS is common in the whole-exome sequencing (WES) datasets (357 tumor samples from 275 patients) from four melanoma cohorts with ICB therapy and a non-small cell lung cancer (NSCLC) WES dataset (327 tumor samples from 100 patients). Altogether, our study provides quantitative evidence supporting the theory of NFDS in cancer, explaining the high prevalence of neutral-looking tumors. These findings also highlight the critical role of frequency-dependent selection in devising more efficient and predictive immunotherapies.

Strength of selection in lung tumors correlates with clinical features better than tumor mutation burden

Single nucleotide substitutions are the most common type of somatic mutations in cancer genome. The goal of this study was to use publicly available somatic mutation data to quantify negative and positive selection in individual lung tumors and test how strength of directional and absolute selection is associated with clinical features. The analysis found a significant variation in strength of selection (both negative and positive) among tumors, with median selection tending to be negative even though tumors with strong positive selection also exist. Strength of selection estimated as the density of missense mutations relative to the density of silent mutations showed only a weak correlation with tumor mutation burden. In the "all histology together" analysis we found that absolute strength of selection was strongly correlated with all clinically relevant features analyzed. In histology-stratified analysis selection was strongest in small cell lung cancer. Selection in adenocarcinoma was somewhat higher compared to squamous cell carcinoma. The study suggests that somatic mutation- based quantifying of directional and absolute selection in individual tumors can be a useful biomarker of tumor aggressiveness.

Machine-learning analysis reveals an important role for negative selection in shaping cancer aneuploidy landscapes

BACKGROUND

Aneuploidy, an abnormal number of chromosomes within a cell, is a hallmark of cancer. Patterns of aneuploidy differ across cancers, yet are similar in cancers affecting closely related tissues. The selection pressures underlying aneuploidy patterns are not fully understood, hindering our understanding of cancer development and progression.

RESULTS

Here, we apply interpretable machine learning methods to study tissue-selective aneuploidy patterns. We define 20 types of features corresponding to genomic attributes of chromosome-arms, normal tissues, primary tumors, and cancer cell lines (CCLs), and use them to model gains and losses of chromosome arms in 24 cancer types. To reveal the factors that shape the tissue-specific cancer aneuploidy landscapes, we interpret the machine learning models by estimating the relative contribution of each feature to the models. While confirming known drivers of positive selection, our quantitative analysis highlights the importance of negative selection for shaping aneuploidy landscapes. This is exemplified by tumor suppressor gene density being a better predictor of gain patterns than oncogene density, and vice versa for loss patterns. We also identify the importance of tissue-selective features and demonstrate them experimentally, revealing KLF5 as an important driver for chr13q gain in colon cancer. Further supporting an important role for negative selection in shaping the aneuploidy landscapes, we find compensation by paralogs to be among the top predictors of chromosome arm loss prevalence and demonstrate this relationship for one paralog interaction. Similar factors shape aneuploidy patterns in human CCLs, demonstrating their relevance for aneuploidy research.

CONCLUSIONS

Our quantitative, interpretable machine learning models improve the understanding of the genomic properties that shape cancer aneuploidy landscapes.

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.

Incorporating mutational heterogeneity to identify genes that are enriched for synonymous mutations in cancer

BACKGROUND

Synonymous mutations, which change the DNA sequence but not the encoded protein sequence, can affect protein structure and function, mRNA maturation, and mRNA half-lives. The possibility that synonymous mutations might be enriched in cancer has been explored in several recent studies. However, none of these studies control for all three types of mutational heterogeneity (patient, histology, and gene) that are known to affect the accurate identification of non-synonymous cancer-associated genes. Our goal is to adopt the current standard for non-synonymous mutations in an investigation of synonymous mutations.

RESULTS

Here, we create an algorithm, MutSigCVsyn, an adaptation of MutSigCV, to identify cancer-associated genes that are enriched for synonymous mutations based on a non-coding background model that takes into account the mutational heterogeneity across these levels. Using MutSigCVsyn, we first analyzed 2572 cancer whole-genome samples from the Pan-cancer Analysis of Whole Genomes (PCAWG) to identify non-synonymous cancer drivers as a quality control. Indicative of the algorithm accuracy we find that 58.6% of these candidate genes were also found in Cancer Census Gene (CGC) list, and 66.2% were found within the PCAWG cancer driver list. We then applied it to identify 30 putative cancer-associated genes that are enriched for synonymous mutations within the same samples. One of the promising gene candidates is the B cell lymphoma 2 (BCL-2) gene. BCL-2 regulates apoptosis by antagonizing the action of proapoptotic BCL-2 family member proteins. The synonymous mutations in BCL2 are enriched in its anti-apoptotic domain and likely play a role in cancer cell proliferation.

CONCLUSION

Our study introduces MutSigCVsyn, an algorithm that accounts for mutational heterogeneity at patient, histology, and gene levels, to identify cancer-associated genes that are enriched for synonymous mutations using whole genome sequencing data. We identified 30 putative candidate genes that will benefit from future experimental studies on the role of synonymous mutations in cancer biology.

Paralog dispensability shapes homozygous deletion patterns in tumor genomes

Genomic instability is a hallmark of cancer, resulting in tumor genomes having large numbers of genetic aberrations, including homozygous deletions of protein coding genes. That tumor cells remain viable in the presence of such gene loss suggests high robustness to genetic perturbation. In model organisms and cancer cell lines, paralogs have been shown to contribute substantially to genetic robustness-they are generally more dispensable for growth than singletons. Here, by analyzing copy number profiles of > 10,000 tumors, we test the hypothesis that the increased dispensability of paralogs shapes tumor genome evolution. We find that genes with paralogs are more likely to be homozygously deleted and that this cannot be explained by other factors known to influence copy number variation. Furthermore, features that influence paralog dispensability in cancer cell lines correlate with paralog deletion frequency in tumors. Finally, paralogs that are broadly essential in cancer cell lines are less frequently deleted in tumors than non-essential paralogs. Overall, our results suggest that homozygous deletions of paralogs are more frequently observed in tumor genomes because paralogs are more dispensable.

Developmental Selection and the Perception of Mutation Bias

The notion that mutations are random relative to their fitness effects is central to the Neo-Darwinian view of evolution. However, a recent interpretation of the patterns of mutation accumulation in the genome of Arabidopsis thaliana has challenged this notion, arguing for the presence of a targeted DNA repair mechanism that causes a nonrandom association of mutation rates and fitness effects. Specifically, this mechanism was suggested to cause a reduction in the rates of mutations on essential genes, thus lowering the rates of deleterious mutations. Central to this argument were attempts to rule out selection at the population level. Here, we offer an alternative and parsimonious interpretation of the patterns of mutation accumulation previously attributed to mutation bias, showing how they can instead or additionally be caused by developmental selection, that is selection occurring at the cellular level during the development of a multicellular organism. Thus, the depletion of deleterious mutations in A. thaliana may indeed be the result of a selective process, rather than a bias in mutation. More broadly, our work highlights the importance of considering development in the interpretation of population-genetic analyses of multicellular organisms, and it emphasizes that efforts to identify mechanisms involved in mutational biases should explicitly account for developmental selection.

Precise microdissection of gastric mixed adeno-neuroendocrine carcinoma dissects its genomic landscape and evolutionary clonal origins

Gastric mixed adenoneuroendocrine carcinoma (MANEC) is a clinically aggressive and heterogeneous tumor composed of adenocarcinoma (ACA) and neuroendocrine carcinoma (NEC). The genomic properties and evolutionary clonal origins of MANEC remain unclear. We conduct whole-exome and multiregional sequencing on 101 samples from 33 patients to elucidate their evolutionary paths. We identify four significantly mutated genes, TP53, RB1, APC, and CTNNB1. MANEC resembles chromosomal instability stomach adenocarcinoma in that whole-genome doubling in MANEC is predominant and occurs earlier than most copy-number losses. All tumors are of monoclonal origin, and NEC components show more aggressive genomic properties than their ACA counterparts. The phylogenetic trees show two tumor divergence patterns, including sequential and parallel divergence. Furthermore, ACA-to-NEC rather than NEC-to-ACA transition is confirmed by immunohistochemistry on 6 biomarkers in ACA- and NEC-dominant regions. These results provide insights into the clonal origin and tumor differentiation of MANEC.

Cancer genomes tolerate deleterious coding mutations through somatic copy number amplifications of wild-type regions

Cancers evolve under the accumulation of thousands of somatic mutations and chromosomal aberrations. While most coding mutations are deleterious, almost all protein-coding genes lack detectable signals of negative selection. This raises the question of how tumors tolerate such large amounts of deleterious mutations. Using 8,690 tumor samples from The Cancer Genome Atlas, we demonstrate that copy number amplifications frequently cover haploinsufficient genes in mutation-prone regions. This could increase tolerance towards the deleterious impact of mutations by creating safe copies of wild-type regions and, hence, protecting the genes therein. Our findings demonstrate that these potential buffering events are highly influenced by gene functions, essentiality, and mutation impact and that they occur early during tumor evolution. We show how cancer type-specific mutation landscapes drive copy number alteration patterns across cancer types. Ultimately, our work paves the way for the detection of novel cancer vulnerabilities by revealing genes that fall within amplifications likely selected during evolution to mitigate the effect of mutations.

Unraveling the Drivers of Tumorigenesis in the Context of Evolution: Theoretical Models and Bioinformatics Tools

Cancer originates from somatic cells that have accumulated mutations. These mutations alter the phenotype of the cells, allowing them to escape homeostatic regulation that maintains normal cell numbers. The emergence of malignancies is an evolutionary process in which the random accumulation of somatic mutations and sequential selection of dominant clones cause cancer cells to proliferate. The development of technologies such as high-throughput sequencing has provided a powerful means to measure subclonal evolutionary dynamics across space and time. Here, we review the patterns that may be observed in cancer evolution and the methods available for quantifying the evolutionary dynamics of cancer. An improved understanding of the evolutionary trajectories of cancer will enable us to explore the molecular mechanism of tumorigenesis and to design tailored treatment strategies.

Distinct Traits of Structural and Regulatory Evolutional Conservation of Human Genes with Specific Focus on Major Cancer Molecular Pathways

The evolution of protein-coding genes has both structural and regulatory components. The first can be assessed by measuring the ratio of non-synonymous to synonymous nucleotide substitutions. The second component can be measured as the normalized proportion of transposable elements that are used as regulatory elements. For the first time, we characterized in parallel the regulatory and structural evolutionary profiles for 10,890 human genes and 2972 molecular pathways. We observed a ~0.1 correlation between the structural and regulatory metrics at the gene level, which appeared much higher (~0.4) at the pathway level. We deposited the data in the publicly available database RetroSpect. We also analyzed the evolutionary dynamics of six cancer pathways of two major axes: Notch/WNT/Hedgehog and AKT/mTOR/EGFR. The Hedgehog pathway had both components slower, whereas the Akt pathway had clearly accelerated structural evolution. In particular, the major hub nodes Akt and beta-catenin showed both components strongly decreased, whereas two major regulators of Akt TCL1 and CTMP had outstandingly high evolutionary rates. We also noticed structural conservation of serine/threonine kinases and the genes related to guanosine metabolism in cancer signaling: GPCRs, G proteins, and small regulatory GTPases (Src, Rac, Ras); however, this was compensated by the accelerated regulatory evolution.

Cancer driver mutations: predictions and reality

Cancer cells accumulate many genetic alterations throughout their lifetime, but only a few of them drive cancer progression, termed driver mutations. Driver mutations may vary between cancer types and patients, can remain latent for a long time and become drivers at particular cancer stages, or may drive oncogenesis only in conjunction with other mutations. The high mutational, biochemical, and histological tumor heterogeneity makes driver mutation identification very challenging. In this review we summarize recent efforts to identify driver mutations in cancer and annotate their effects. We underline the success of computational methods to predict driver mutations in finding novel cancer biomarkers, including in circulating tumor DNA (ctDNA). We also report on the boundaries of their applicability in clinical research.

The role of physics in multiomics and cancer evolution

Complex interactions between the physical environment and phenotype of a tumour, and genomics, transcriptomics, proteomics and epigenomics, are increasingly known to have a significant influence on cancer development, progression and evolution. For example, mechanical stress can alter both genome maintenance and histone modifications, which consequently affect transcription and the epigenome. Increased stiffness has been linked to genetic heterogeneity and is responsible for heterochromatin accumulations. Stiffness thereby leads to deregulation in gene expression, disrupts the proteome and can impact angiogenesis. Several studies have shown how the physics of cancer can influence diverse cancer hallmarks such as resistance to cell death, angiogenesis and evasion from immune destruction. In this review, we will explain the role that physics of cancer plays in cancer evolution and explore how multiomics are being used to elucidate the mechanisms underpinning them.

Immune selection determines tumor antigenicity and influences response to checkpoint inhibitors

In cancer, evolutionary forces select for clones that evade the immune system. Here we analyzed >10,000 primary tumors and 356 immune-checkpoint-treated metastases using immune dN/dS, the ratio of nonsynonymous to synonymous mutations in the immunopeptidome, to measure immune selection in cohorts and individuals. We classified tumors as immune edited when antigenic mutations were removed by negative selection and immune escaped when antigenicity was covered up by aberrant immune modulation. Only in immune-edited tumors was immune predation linked to CD8 T cell infiltration. Immune-escaped metastases experienced the best response to immunotherapy, whereas immune-edited patients did not benefit, suggesting a preexisting resistance mechanism. Similarly, in a longitudinal cohort, nivolumab treatment removes neoantigens exclusively in the immunopeptidome of nonimmune-edited patients, the group with the best overall survival response. Our work uses dN/dS to differentiate between immune-edited and immune-escaped tumors, measuring potential antigenicity and ultimately helping predict response to treatment.

First passage time analysis of spatial mutation patterns reveals sub-clonal evolutionary dynamics in colorectal cancer

The signature of early cancer dynamics on the spatial arrangement of tumour cells is poorly understood, and yet could encode information about how sub-clones grew within the expanding tumour. Novel methods of quantifying spatial tumour data at the cellular scale are required to link evolutionary dynamics to the resulting spatial architecture of the tumour. Here, we propose a framework using first passage times of random walks to quantify the complex spatial patterns of tumour cell population mixing. First, using a simple model of cell mixing we demonstrate how first passage time statistics can distinguish between different pattern structures. We then apply our method to simulated patterns of mutated and non-mutated tumour cell population mixing, generated using an agent-based model of expanding tumours, to explore how first passage times reflect mutant cell replicative advantage, time of emergence and strength of cell pushing. Finally, we explore applications to experimentally measured human colorectal cancer, and estimate parameters of early sub-clonal dynamics using our spatial computational model. We infer a wide range of sub-clonal dynamics, with mutant cell division rates varying between 1 and 4 times the rate of non-mutated cells across our sample set. Some mutated sub-clones emerged after as few as 100 non-mutant cell divisions, and others only after 50,000 divisions. The majority were consistent with boundary driven growth or short-range cell pushing. By analysing multiple sub-sampled regions in a small number of samples, we explore how the distribution of inferred dynamics could inform about the initial mutational event. Our results demonstrate the efficacy of first passage time analysis as a new methodology in spatial analysis of solid tumour tissue, and suggest that patterns of sub-clonal mixing can provide insights into early cancer dynamics.

MHC II immunogenicity shapes the neoepitope landscape in human tumors

Despite advances in predicting physical peptide-major histocompatibility complex I (pMHC I) binding, it remains challenging to identify functionally immunogenic neoepitopes, especially for MHC II. By using the results of >36,000 immunogenicity assay, we developed a method to identify pMHC whose structural alignment facilitates T cell reaction. Our method predicted neoepitopes for MHC II and MHC I that were responsive to checkpoint blockade when applied to >1,200 samples of various tumor types. To investigate selection by spontaneous immunity at the single epitope level, we analyzed the frequency spectrum of >25 million mutations in >9,000 treatment-naive tumors with >100 immune phenotypes. MHC II immunogenicity specifically lowered variant frequencies in tumors under high immune pressure, particularly with high TCR clonality and MHC II expression. A similar trend was shown for MHC I neoepitopes, but only in particular tissue types. In summary, we report immune selection imposed by MHC II-restricted natural or therapeutic T cell reactivity.

In Vitro and In Vivo Models for Metastatic Intestinal Tumors Using Genotype-Defined Organoids

It has been established that the accumulation of driver gene mutations causes malignant progression of colorectal cancer (CRC) through positive selection and clonal expansion, similar to Darwin's evolution. Following this multistep tumorigenesis concept, we previously showed the specific mutation patterns for each process of malignant progression, including submucosal invasion, epithelial mesenchymal transition (EMT), intravasation, and metastasis, using genetically engineered mouse and organoid models. However, we also found that certain populations of cancer-derived organoid cells lost malignant characteristics of metastatic ability, although driver mutations were not impaired, and such subpopulations were eliminated from the tumor tissues by negative selection. These organoid model studies have contributed to our understanding of the cancer evolution mechanism. We herein report the in vitro and in vivo experimental protocols to investigate the survival, growth, and metastatic ability of intestinal tumor-derived organoids. The model system will be useful for basic research as well as the development of clinical strategies.

Frequent loss of metastatic ability in subclones of Apc, Kras, Tgfbr2, and Trp53 mutant intestinal tumor organoids

Cancer evolution is explained by the accumulation of driver mutations and subsequent positive selection by acquired growth advantages, like Darwin's evolution theory. However, whether the negative selection of cells that have lost malignant properties contributes to cancer progression has not yet been fully investigated. Using intestinal metastatic tumor-derived organoids carrying Apc, Kras, Tgfbr2, and Trp53 quadruple mutations, we demonstrate here that approximately 30% of subclones of the organoids show loss of metastatic ability to the liver while keeping the driver mutations and oncogenic pathways. Notably, highly metastatic subclones also showed a gradual loss of metastatic ability during further passages. Such non-metastatic subclones revealed significantly decreased survival and proliferation ability in Matrigel and collagen gel culture conditions, which may cause elimination from the tumor tissues in vivo. RNA sequencing indicated that stemness-related genes, including Lgr5 and Myb, were significantly downregulated in non-metastatic subclones as well as subclones that lost metastatic ability during additional passages. Furthermore, a CGH analysis showed that non-metastatic subclones were derived from a minor population of parental organoid cells. These results indicate that metastatic ability is continuously lost with decreased stem cell property in certain subpopulations of malignant tumors, and such subpopulations are eliminated by negative selection. Therefore, it is possible that cancer evolution is regulated not only by positive selection but also by negative selection. The mechanism underlying the loss of metastatic ability will be important for the future development of therapeutic strategies against metastasis.

Evolutionary Analysis of TCGA Data Using Over- and Under- Mutated Genes Identify Key Molecular Pathways and Cellular Functions in Lung Cancer Subtypes

We identify critical conserved and mutated genes through a theoretical model linking a gene’s fitness contribution to its observed mutational frequency in a clinical cohort. “Passenger” gene mutations do not alter fitness and have mutational frequencies determined by gene size and the mutation rate. Driver mutations, which increase fitness (and proliferation), are observed more frequently than expected. Non-synonymous mutations in essential genes reduce fitness and are eliminated by natural selection resulting in lower prevalence than expected. We apply this “evolutionary triage” principle to TCGA data from EGFR-mutant, KRAS-mutant, and NEK (non-EGFR/KRAS) lung adenocarcinomas. We find frequent overlap of evolutionarily selected non-synonymous gene mutations among the subtypes suggesting enrichment for adaptations to common local tissue selection forces. Overlap of conserved genes in the LUAD subtypes is rare suggesting negative evolutionary selection is strongly dependent on initiating mutational events during carcinogenesis. Highly expressed genes are more likely to be conserved and significant changes in expression (>20% increased/decreased) are common in genes with evolutionarily selected mutations but not in conserved genes. EGFR-mut cancers have fewer average mutations (89) than KRAS-mut (228) and NEK (313). Subtype-specific variation in conserved and mutated genes identify critical molecular components in cell signaling, extracellular matrix remodeling, and membrane transporters. These findings demonstrate subtype-specific patterns of co-adaptations between the defining driver mutation and somatically conserved genes as well as novel insights into epigenetic versus genetic contributions to cancer evolution.

The co-evolution of the genome and epigenome in colorectal cancer

Colorectal malignancies are a leading cause of cancer-related death1 and have undergone extensive genomic study2,3. However, DNA mutations alone do not fully explain malignant transformation4-7. Here we investigate the co-evolution of the genome and epigenome of colorectal tumours at single-clone resolution using spatial multi-omic profiling of individual glands. We collected 1,370 samples from 30 primary cancers and 8 concomitant adenomas and generated 1,207 chromatin accessibility profiles, 527 whole genomes and 297 whole transcriptomes. We found positive selection for DNA mutations in chromatin modifier genes and recurrent somatic chromatin accessibility alterations, including in regulatory regions of cancer driver genes that were otherwise devoid of genetic mutations. Genome-wide alterations in accessibility for transcription factor binding involved CTCF, downregulation of interferon and increased accessibility for SOX and HOX transcription factor families, suggesting the involvement of developmental genes during tumourigenesis. Somatic chromatin accessibility alterations were heritable and distinguished adenomas from cancers. Mutational signature analysis showed that the epigenome in turn influences the accumulation of DNA mutations. This study provides a map of genetic and epigenetic tumour heterogeneity, with fundamental implications for understanding colorectal cancer biology.

Most cancers carry a substantial deleterious load due to Hill-Robertson interference

Cancer genomes exhibit surprisingly weak signatures of negative selection^1,2. This may be because selective pressures are relaxed or because genome-wide linkage prevents deleterious mutations from being removed (Hill-Robertson interference)³. By stratifying tumors by their genome-wide mutational burden, we observe negative selection (dN/dS ~ 0.56) in low mutational burden tumors, while remaining cancers exhibit dN/dS ratios ~1. This suggests that most tumors do not remove deleterious passengers. To buffer against deleterious passengers, tumors upregulate heat shock pathways as their mutational burden increases. Finally, evolutionary modeling finds that Hill-Robertson interference alone can reproduce patterns of attenuated selection and estimates the total fitness cost of passengers to be 46% per cell on average. Collectively, our findings suggest that the lack of observed negative selection in most tumors is not due to relaxed selective pressures, but rather the inability of selection to remove deleterious mutations in the presence of genome-wide linkage.

Evolutionary selection identifies critical immune-relevant genes in lung cancer subtypes

In an evolving population, proliferation is dependent on fitness so that a numerically dominant population typically possesses the most well adapted phenotype. In contrast, the evolutionary "losers" typically disappear from the population so that their genetic record is lost. Historically, cancer research has focused on observed genetic mutations in the dominant tumor cell populations which presumably increase fitness. Negative selection, i.e., removal of deleterious mutations from a population, is not observable but can provide critical information regarding genes involved in essential cellular processes. Similar to immunoediting, "evolutionary triage" eliminates mutations in tumor cells that increase susceptibility to the host immune response while mutations that shield them from immune attack increase proliferation and are readily observable (e.g., B2M mutations). These dynamics permit an "inverse problem" analysis linking the fitness consequences of a mutation to its prevalence in a tumor cohort. This is evident in "driver mutations" but, equally important, can identify essential genes in which mutations are seen significantly less than expected by chance. Here we utilized this new approach to investigate evolutionary triage in immune-related genes from TCGA lung adenocarcinoma cohorts. Negative selection differs between the two cohorts and is observed in endoplasmic reticulum aminopeptidase genes, ERAP1 and ERAP2 genes, and DNAM-1/TIGIT ligands. Targeting genes or molecular pathways under positive or negative evolutionary selection may permit new treatment options and increase the efficacy of current immunotherapy.

Cancer evolution: special focus on the immune aspect of cancer

Cancer is an evolutionary disease. Intra-tumor heterogeneity (ITH), which describes the diversity within individual tumors, sets the foundation for evolution. The fitness of tumor cells is determined by their microenvironment, which exerts intense selection pressure that generally favors cells with survival and proliferation advantages. It has been revealed that host immunity dramatically influences the evolutionary trajectory of cancer. As technologies advance, a refined map of the immune system's involvement in cancer evolution has gradually come to our knowledge. Here we specifically view cancer through the lens of evolutionary immunological biology. We will cover the neoplastic evolution under immunosurveillance, including how the host immunity shapes the tumor evolutionary trajectory and how progressive tumors modulate the host immunity to survive. A comprehensive understanding of the interplay between cancer evolution and cancer immunity provides clues to combating cancer strategically.

Quantification of neoantigen-mediated immunoediting in cancer evolution

Immunoediting includes three temporally distinct stages, termed elimination, equilibrium, and escape, and has been proposed to explain the interactions between cancer cells and the immune system during the evolution of cancer. However, the status of immunoediting in cancer remains unclear, and the existence of neoantigen depletion in untreated cancer has been debated. Here we developed a distribution pattern-based method for quantifying neoantigen-mediated negative selection in cancer evolution. The method can provide a robust and reliable quantification for immunoediting signal in individual cancer patients. Moreover, this method demonstrated the prevalence of immunoediting in immunotherapy-naive cancer genome. The elimination and escape stages of immunoediting can be quantified separately, where tumor types with strong immunoediting-elimination exhibit a weak immunoediting-escape signal, and vice versa. The quantified immunoediting-elimination signal was predictive of clinical response to cancer immunotherapy. Collectively, immunoediting quantification provides an evolutionary perspective for evaluating the antigenicity of neoantigens and reveals a potential biomarker for precision immunotherapy in cancer.

Identification of neoantigens for individualized therapeutic cancer vaccines

Somatic mutations in cancer cells can generate tumour-specific neoepitopes, which are recognized by autologous T cells in the host. As neoepitopes are not subject to central immune tolerance and are not expressed in healthy tissues, they are attractive targets for therapeutic cancer vaccines. Because the vast majority of cancer mutations are unique to the individual patient, harnessing the full potential of this rich source of targets requires individualized treatment approaches. Many computational algorithms and machine-learning tools have been developed to identify mutations in sequence data, to prioritize those that are more likely to be recognized by T cells and to design tailored vaccines for every patient. In this Review, we fill the gaps between the understanding of basic mechanisms of T cell recognition of neoantigens and the computational approaches for discovery of somatic mutations and neoantigen prediction for cancer immunotherapy. We present a new classification of neoantigens, distinguishing between guarding, restrained and ignored neoantigens, based on how they confer proficient antitumour immunity in a given clinical context. Such context-based differentiation will contribute to a framework that connects neoantigen biology to the clinical setting and medical peculiarities of cancer, and will enable future neoantigen-based therapies to provide greater clinical benefit.

Influence of T Cell-Mediated Immune Surveillance on Somatic Mutation Occurrences in Melanoma

Background

Neoantigens are presented on the cancer cell surface by peptide-restricted human leukocyte antigen (HLA) proteins and can subsequently activate cognate T cells. It has been hypothesized that the observed somatic mutations in tumors are shaped by immunosurveillance.

Methods

We investigated all somatic mutations identified in The Cancer Genome Atlas (TCGA) Skin Cutaneous Melanoma (SKCM) samples. By applying a computational algorithm, we calculated the binding affinity of the resulting neo-peptides and their corresponding wild-type peptides with the major histocompatibility complex (MHC) Class I complex. We then examined the relationship between binding affinity alterations and mutation frequency.

Results

Our results show that neoantigens derived from recurrent mutations tend to have lower binding affinities with the MHC Class I complex compared to peptides from non-recurrent mutations. Tumor samples harboring recurrent SKCM mutations exhibited lower immune infiltration levels, indicating a relatively colder immune microenvironment.

Conclusions

These results suggested that the occurrences of somatic mutations in melanoma have been shaped by immunosurveillance. Mutations that lead to neoantigens with high MHC class I binding affinity are more likely to be eliminated and thus are less likely to be present in tumors.

Extreme value theory as a framework for understanding mutation frequency distribution in cancer genomes

Currently, the population dynamics of preclonal cancer cells before clonal expansion of tumors has not been sufficiently addressed thus far. By focusing on preclonal cancer cell population as a Darwinian evolutionary system, we formulated and analyzed the observed mutation frequency among tumors (MFaT) as a proxy for the hypothesized sequence read frequency and beneficial fitness effect of a cancer driver mutation. Analogous to intestinal crypts, we assumed that sample donor patients are separate culture tanks where proliferating cells follow certain population dynamics described by extreme value theory (EVT). To validate this, we analyzed three large-scale cancer genome datasets, each harboring > 10000 tumor samples and in total involving > 177898 observed mutation sites. We clarified the necessary premises for the application of EVT in the strong selection and weak mutation (SSWM) regime in relation to cancer genome sequences at scale. We also confirmed that the stochastic distribution of MFaT is likely of the Fréchet type, which challenges the well-known Gumbel hypothesis of beneficial fitness effects. Based on statistical data analysis, we demonstrated the potential of EVT as a population genetics framework to understand and explain the stochastic behavior of driver-mutation frequency in cancer genomes as well as its applicability in real cancer genome sequence data.

No evidence that HLA genotype influences the driver mutations that occur in cancer patients

The major histocompatibility (MHC) molecules are capable of presenting neoantigens resulting from somatic mutations on cell surfaces, potentially directing immune responses against cancer. This led to the hypothesis that cancer driver mutations may occur in gaps in the capacity to present neoantigens that are dependent on MHC genotype. If this is correct, it has important implications for understanding oncogenesis and may help to predict driver mutations based on genotype data. In support of this hypothesis, it has been reported that driver mutations that occur frequently tend to be poorly presented by common MHC alleles and that the capacity of a patient’s MHC alleles to present the resulting neoantigens is predictive of the driver mutations that are observed in their tumor. Here we show that these reports of a strong relationship between driver mutation occurrence and patient MHC alleles are a consequence of unjustified statistical assumptions. Our reanalysis of the data provides no evidence of an effect of MHC genotype on the oncogenic mutation landscape.

Improvement of Enzymatic Saccharification of Cellulose-Containing Raw Materials Using Aspergillus niger

Enzymatic hydrolysis of cellulose-containing raw materials, using Aspergillus niger, were studied. Filter paper, secondary cellulose-containing or starch-containing raw materials, miscanthus cellulose after alkaline or acid pretreatment, and wood chip cellulose, were used as substrates. The study focused on a wild A. niger strain, treated, or not (control), by ultraviolet (UV) irradiations for 45, 60, or 120 min (UV45, UV60, or UV120), or by UV irradiation for 120 min followed by a chemical treatment with NaN3 + ItBr for 30 min or 80 min (UV120 + CH30 or UV120 + CH80). A mixture of all the A. niger strains (MIX) was also tested. A citrate buffer, at 50 mM, wasthe most suitable for enzymatic hydrolysis. As the UV exposure time increased to 2 h, the cellulase activity of the surviving culturewas increased (r = 0.706; p < 0.05). The enzymatic activities of the obtained strains, towards miscanthus cellulose, wood chips, and filter paper, were inferior to those obtained with commercial enzymes (8.6 versus 9.1 IU), in some cases. Under stationary hydrolysis at 37 °C, pH = 4.7, the enzymatic activity of A. niger UV120 + CH30 was 24.9 IU. The enzymatic hydrolysis of secondary raw materials, using treated A. niger strains, was themost effective at 37 °C. Similarly, the most effective treatment of miscanthus cellulose and wood chips occurred at 50 °C. The maximum conversion of cellulose to glucose was observed using miscanthus cellulose (with alkaline pretreatment), and the minimum conversion was observed when using wood chips. The greatest value of cellulase activity was evidenced in the starch-containing raw materials, indicating that A. niger can ferment not only through cellulase activity, but also via an amylolytic one.

Genome-wide profiles of UV lesion susceptibility, repair, and mutagenic potential in melanoma

Exposure to the ultraviolet (UV) radiation in sunlight creates DNA lesions, which if left unrepaired can induce mutations and contribute to skin cancer. The two most common UV-induced DNA lesions are the cis-syn cyclobutane pyrimidine dimers (CPDs) and pyrimidine (6-4) pyrimidone photoproducts (6-4PPs), both of which can initiate mutations. Interestingly, mutation frequency across the genomes of many cancers is heterogenous with significant increases in heterochromatin. Corresponding increases in UV lesion susceptibility and decreases in repair are observed in heterochromatin versus euchromatin. However, the individual contributions of CPDs and 6-4PPs to mutagenesis have not been systematically examined in specific genomic and epigenomic contexts. In this study, we compared genome-wide maps of 6-4PP and CPD lesion abundances in primary cells and conducted comprehensive analyses to determine the genetic and epigenetic features associated with susceptibility. Overall, we found a high degree of similarity between 6-4PP and CPD formation, with an enrichment of both in heterochromatin regions. However, when examining the relative levels of the two UV lesions, we found that bivalent and Polycomb-repressed chromatin states were uniquely more susceptible to 6-4PPs. Interestingly, when comparing UV susceptibility and repair with melanoma mutation frequency in these regions, disparate patterns were observed in that susceptibility was not always inversely associated with repair and mutation frequency. Functional enrichment analysis hint at mechanisms of negative selection for these regions that are essential for cell viability, immune function and induce cell death when mutated. Ultimately, these results reveal both the similarities and differences between UV-induced lesions that contribute to melanoma.

Identifying Cancer Drivers Using DRIVE: A Feature-Based Machine Learning Model for a Pan-Cancer Assessment of Somatic Missense Mutations

Sporadic cancer develops from the accrual of somatic mutations. Out of all small-scale somatic aberrations in coding regions, 95% are base substitutions, with 90% being missense mutations. While multiple studies focused on the importance of this mutation type, a machine learning method based on the number of protein-protein interactions (PPIs) has not been fully explored. This study aims to develop an improved computational method for driver identification, validation and evaluation (DRIVE), which is compared to other methods for assessing its performance. DRIVE aims at distinguishing between driver and passenger mutations using a feature-based learning approach comprising two levels of biological classification for a pan-cancer assessment of somatic mutations. Gene-level features include the maximum number of protein-protein interactions, the biological process and the type of post-translational modifications (PTMs) while mutation-level features are based on pathogenicity scores. Multiple supervised classification algorithms were trained on Genomics Evidence Neoplasia Information Exchange (GENIE) project data and then tested on an independent dataset from The Cancer Genome Atlas (TCGA) study. Finally, the most powerful classifier using DRIVE was evaluated on a benchmark dataset, which showed a better overall performance compared to other state-of-the-art methodologies, however, considerable care must be taken due to the reduced size of the dataset. DRIVE outlines the outstanding potential that multiple levels of a feature-based learning model will play in the future of oncology-based precision medicine.

Genetic and non-genetic clonal diversity in cancer evolution

The observation and analysis of intra-tumour heterogeneity (ITH), particularly in genomic studies, has advanced our understanding of the evolutionary forces that shape cancer growth and development. However, only a subset of the variation observed in a single tumour will have an impact on cancer evolution, highlighting the need to distinguish between functional and non-functional ITH. Emerging studies highlight a role for the cancer epigenome, transcriptome and immune microenvironment in functional ITH. Here, we consider the importance of both genetic and non-genetic ITH and their role in tumour evolution, and present the rationale for a broad research focus beyond the cancer genome. Systems-biology analytical approaches will be necessary to outline the scale and importance of functional ITH. By allowing a deeper understanding of tumour evolution this will, in time, encourage development of novel therapies and improve outcomes for patients.

The Burden of Post-Translational Modification (PTM)-Disrupting Mutations in the Tumor Matrisome

BACKGROUND

To evaluate the occurrence of mutations affecting post-translational modification (PTM) sites in matrisome genes across different tumor types, in light of their genomic and functional contexts and in comparison with the rest of the genome.

METHODS

This study spans 9075 tumor samples and 32 tumor types from The Cancer Genome Atlas (TCGA) Pan-Cancer cohort and identifies 151,088 non-silent mutations in the coding regions of the matrisome, of which 1811 affecting known sites of hydroxylation, phosphorylation, N- and O-glycosylation, acetylation, ubiquitylation, sumoylation and methylation PTM.

RESULTS

PTM-disruptive mutations (PTMmut) in the matrisome are less frequent than in the rest of the genome, seem independent of cell-of-origin patterns but show dependence on the nature of the matrisome protein affected and the background PTM types it generally harbors. Also, matrisome PTMmut are often found among structural and functional protein regions and in proteins involved in homo- and heterotypic interactions, suggesting potential disruption of matrisome functions.

CONCLUSIONS

Though quantitatively minoritarian in the spectrum of matrisome mutations, PTMmut show distinctive features and damaging potential which might concur to deregulated structural, functional, and signaling networks in the tumor microenvironment.

Low immunogenicity of common cancer hot spot mutations resulting in false immunogenic selection signals

Cancer is driven by somatic mutations that result in a cellular fitness advantage. This selective advantage is expected to be counterbalanced by the immune system when these driver mutations simultaneously lead to the generation of neoantigens, novel peptides that are presented at the cancer cell membrane via HLA molecules from the MHC complex. The presentability of these peptides is determined by a patient’s MHC genotype and it has been suggested that this results in MHC genotype-specific restrictions of the oncogenic mutational landscape. Here, we generated a set of virtual patients, each with an identical and prototypical MHC genotype, and show that the earlier reported HLA affinity differences between observed and unobserved mutations are unrelated to MHC genotype variation. We demonstrate how these differences are secondary to high frequencies of 13 hot spot driver mutations in 6 different genes. Several oncogenic mechanisms were identified that lower the peptides’ HLA affinity, including phospho-mimicking substitutions in BRAF, destabilizing tyrosine mutations in TP53 and glycine-rich mutational contexts in the GTP-binding KRAS domain. In line with our earlier findings, our results emphasize that HLA affinity predictions are easily misinterpreted when studying immunogenic selection processes.

The diagnosis of a malignant tumor is preceded by a long process of tumor evolution, characterized by the gradual accumulation of driver mutations, genomic alterations that give cancer cells their typical growth advantage. This process is controlled by the immune system and understanding tumor-immune interactions is critical for the development of new anti-cancer therapies. Immune cells mainly respond to neoantigens, small mutated peptides that are presented at the cancer cell membrane by binding to the Major Histocompatibility Complex (MHC). MHC genes are highly variable between individuals and it was recently suggested that an individual’s MHC genotype determines cancer susceptibility. In this study we used a computational approach and simulated a set of virtual cancer patients, based on the expected driver mutation frequencies and each with a similar prototypical MHC genotype. Using these simulations, we show that the earlier perceived signals are unrelated to the underlying MHC genotypes, but rather secondary to high frequencies of 13 driver mutations in 6 cancer genes.

Prediction of neo-epitope immunogenicity reveals TCR recognition determinants and provides insight into immunoediting

CD8+ T cell recognition of peptide epitopes plays a central role in immune responses against pathogens and tumors. However, the rules that govern which peptides are truly recognized by existing T cell receptors (TCRs) remain poorly understood, precluding accurate predictions of neo-epitopes for cancer immunotherapy. Here, we capitalize on recent (neo-)epitope data to train a predictor of immunogenic epitopes (PRIME), which captures molecular properties of both antigen presentation and TCR recognition. PRIME not only improves prioritization of neo-epitopes but also correlates with T cell potency and unravels biophysical determinants of TCR recognition that we experimentally validate. Analysis of cancer genomics data reveals that recurrent mutations tend to be less frequent in patients where they are predicted to be immunogenic, providing further evidence for immunoediting in human cancer. PRIME will facilitate identification of pathogen epitopes in infectious diseases and neo-epitopes in cancer immunotherapy.

Targeting Loss of Heterozygosity: A Novel Paradigm for Cancer Therapy

Loss of heterozygosity (LOH) is a common genetic event in the development of cancer. In certain tumor types, LOH can affect more than 20% of the genome, entailing loss of allelic variation in thousands of genes. This reduction of heterozygosity creates genetic differences between tumor and normal cells, providing opportunities for development of novel cancer therapies. Here, we review and summarize (1) mutations associated with LOH on chromosomes which have been shown to be promising biomarkers of cancer risk or the prediction of clinical outcomes in certain types of tumors; (2) loci undergoing LOH that can be targeted for development of novel anticancer drugs as well as (3) LOH in tumors provides up-and-coming possibilities to understand the underlying mechanisms of cancer evolution and to discover novel cancer vulnerabilities which are worth a further investigation in the near future.

Use of signals of positive and negative selection to distinguish cancer genes and passenger genes

A major goal of cancer genomics is to identify all genes that play critical roles in carcinogenesis. Most approaches focused on genes positively selected for mutations that drive carcinogenesis and neglected the role of negative selection. Some studies have actually concluded that negative selection has no role in cancer evolution. We have re-examined the role of negative selection in tumor evolution through the analysis of the patterns of somatic mutations affecting the coding sequences of human genes. Our analyses have confirmed that tumor suppressor genes are positively selected for inactivating mutations, oncogenes, however, were found to display signals of both negative selection for inactivating mutations and positive selection for activating mutations. Significantly, we have identified numerous human genes that show signs of strong negative selection during tumor evolution, suggesting that their functional integrity is essential for the growth and survival of tumor cells.

Extensive heterogeneity in somatic mutation and selection in the human bladder

The extent of somatic mutation and clonal selection in the human bladder remains unknown. We sequenced 2097 bladder microbiopsies from 20 individuals using targeted (n = 1914 microbiopsies), whole-exome (n = 655), and whole-genome (n = 88) sequencing. We found widespread positive selection in 17 genes. Chromatin remodeling genes were frequently mutated, whereas mutations were absent in several major bladder cancer genes. There was extensive interindividual variation in selection, with different driver genes dominating the clonal landscape across individuals. Mutational signatures were heterogeneous across clones and individuals, which suggests differential exposure to mutagens in the urine. Evidence of APOBEC mutagenesis was found in 22% of the microbiopsies. Sequencing multiple microbiopsies from five patients with bladder cancer enabled comparisons with cancer-free individuals and across histological features. This study reveals a rich landscape of mutational processes and selection in normal urothelium with large heterogeneity across clones and individuals.

Generation of multiepitope cancer vaccines based on large combinatorial libraries of survivin-derived mutant epitopes

Immune tolerance is the main challenge in the field of cancer vaccines, so modified peptide sequences or naturally occurring mutated versions of cancer-related wild-type (WT) antigens represent a promising pathway. However, the low immunogenicity of mutation-induced neoantigens and, particularly, their incapacity to activate CD8+ T cells are generating doubts on the success of neoantigen-based cancer vaccines in clinical trials. We developed a novel vaccine approach based on a new class of vaccine immunogens, called variable epitope libraries (VELs). We used three regions of survivin (SVN), composed of 40, 49 and 51 amino acids, along with the complete SVN protein to generate the VELs as multiepitope vaccines. BALB/c mice, challenged with the aggressive and highly metastatic 4T1 cell line, were vaccinated in a therapeutic setting. We showed significant tumor growth inhibition and, most importantly, strong suppression of lung metastasis after a single immunization using VEL vaccines. We demonstrated vaccine-induced broad cellular immune responses concomitant with extensive tumor infiltration of T cells, the activation of CD107a+  IFN-γ+ T cells in the spleen and a significant increase in the number of CD3+  CD8+  Ly6C+ effector T cells. In addition, we observed the presence of interferon-γ-, granzyme B- and perforin-producing lymphocytes along with modifications in the amount of CD11b+  Ly6Cint/low  Ly6G+ granulocytic myeloid-derived suppressor cells and CD4+  CD25+  FoxP3+ regulatory T cells in the lungs and tumors of mice. In summary, we showed that the VELs represent a potent new class of cancer immunotherapy and propose the application of the VEL vaccine concept as a true alternative to currently available vaccine platforms.

Evolutionary dynamics of neoantigens in growing tumors

Cancers accumulate mutations that lead to neoantigens, novel peptides that elicit an immune response, and consequently undergo evolutionary selection. Here we establish how negative selection shapes the clonality of neoantigens in a growing cancer by constructing a mathematical model of neoantigen evolution. The model predicts that, without immune escape, tumor neoantigens are either clonal or at low frequency; hypermutated tumors can only establish after the evolution of immune escape. Moreover, the site frequency spectrum of somatic variants under negative selection appears more neutral as the strength of negative selection increases, which is consistent with classical neutral theory. These predictions are corroborated by the analysis of neoantigen frequencies and immune escape in exome and RNA sequencing data from 879 colon, stomach and endometrial cancers.

Subclonal reconstruction of tumors by using machine learning and population genetics

Most cancer genomic data are generated from bulk samples composed of mixtures of cancer subpopulations, as well as normal cells. Subclonal reconstruction methods based on machine learning aim to separate those subpopulations in a sample and infer their evolutionary history. However, current approaches are entirely data driven and agnostic to evolutionary theory. We demonstrate that systematic errors occur in the analysis if evolution is not accounted for, and this is exacerbated with multi-sampling of the same tumor. We present a novel approach for model-based tumor subclonal reconstruction, called MOBSTER, which combines machine learning with theoretical population genetics. Using public whole-genome sequencing data from 2,606 samples from different cohorts, new data and synthetic validation, we show that this method is more robust and accurate than current techniques in single-sample, multiregion and longitudinal data. This approach minimizes the confounding factors of nonevolutionary methods, thus leading to more accurate recovery of the evolutionary history of human cancers.