Estimating growth patterns and driver effects in tumor evolution from individual samples

Parameter estimation from single patient, single time-point sequencing data of recurrent tumors

In this study, we develop consistent estimators for key parameters that govern the dynamics of tumor cell populations when subjected to pharmacological treatments. While these treatments often lead to an initial reduction in the abundance of drug-sensitive cells, a population of drug-resistant cells frequently emerges over time, resulting in cancer recurrence. Samples from recurrent tumors present as an invaluable data source that can offer crucial insights into the ability of cancer cells to adapt and withstand treatment interventions. To effectively utilize the data obtained from recurrent tumors, we derive several large number limit theorems, specifically focusing on the metrics that quantify the clonal diversity of cancer cell populations at the time of cancer recurrence. These theorems then serve as the foundation for constructing our estimators. A distinguishing feature of our approach is that our estimators only require a single time-point sequencing data from a single tumor, thereby enhancing the practicality of our approach and enabling the understanding of cancer recurrence at the individual level.

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.

Latent evolutionary signatures: a general framework for analysing music and cultural evolution

Cultural processes of change bear many resemblances to biological evolution. The underlying units of non-biological evolution have, however, remained elusive, especially in the domain of music. Here, we introduce a general framework to jointly identify underlying units and their associated evolutionary processes. We model musical styles and principles of organization in dimensions such as harmony and form as following an evolutionary process. Furthermore, we propose that such processes can be identified by extracting latent evolutionary signatures from musical corpora, analogously to identifying mutational signatures in genomics. These signatures provide a latent embedding for each song or musical piece. We develop a deep generative architecture for our model, which can be viewed as a type of variational autoencoder with an evolutionary prior constraining the latent space; specifically, the embeddings for each song are tied together via an energy-based prior, which encourages songs close in evolutionary space to share similar representations. As illustration, we analyse songs from the McGill Billboard dataset. We find frequent chord transitions and formal repetition schemes and identify latent evolutionary signatures related to these features. Finally, we show that the latent evolutionary representations learned by our model outperform non-evolutionary representations in such tasks as period and genre prediction.

Quantification of rare somatic single nucleotide variants by droplet digital PCR using SuperSelective primers

Somatic single-nucleotide variants (SNVs) occur every time a cell divides, appearing even in healthy tissues at low frequencies. These mutations may accumulate as neutral variants during aging, or eventually, promote the development of neoplasia. Here, we present the SP-ddPCR, a droplet digital PCR (ddPCR) based approach that utilizes customized SuperSelective primers aiming at quantifying the proportion of rare SNVs. For that purpose, we selected five potentially pathogenic variants identified by whole-exome sequencing (WES) occurring at low variant allele frequency (VAF) in at-risk colon healthy mucosa of patients diagnosed with colorectal cancer or advanced adenoma. Additionally, two APC SNVs detected in two cancer lesions were added to the study for WES-VAF validation. SuperSelective primers were designed to quantify SNVs at low VAFs both in silico and in clinical samples. In addition to the two APC SNVs in colonic lesions, SP-ddPCR confirmed the presence of three out of five selected SNVs in the normal colonic mucosa with allelic frequencies ≤ 5%. Moreover, SP-ddPCR showed the presence of two potentially pathogenic variants in the distal normal mucosa of patients with colorectal carcinoma. In summary, SP-ddPCR offers a rapid and feasible methodology to validate next-generation sequencing data and accurately quantify rare SNVs, thus providing a potential tool for diagnosis and stratification of at-risk patients based on their mutational profiling.

Netie: inferring the evolution of neoantigen-T cell interactions in tumors

Neoantigens are the key targets of antitumor immune responses from cytotoxic T cells and play a critical role in affecting tumor progressions and immunotherapy treatment responses. However, little is known about how the interaction between neoantigens and T cells ultimately affects the evolution of cancerous masses. Here, we develop a hierarchical Bayesian model, named neoantigen-T cell interaction estimation (netie) to infer the history of neoantigen-CD8⁺ T cell interactions in tumors. Netie was systematically validated and applied to examine the molecular patterns of 3,219 tumors, compiled from a panel of 18 cancer types. We showed that tumors with an increase in immune selection pressure over time are associated with T cells that have an activation-related expression signature. We also identified a subset of exhausted cytotoxic T cells postimmunotherapy associated with tumor clones that newly arise after treatment. These analyses demonstrate how netie enables the interrogation of the relationship between individual neoantigen repertoires and the tumor molecular profiles. We found that a T cell inflammation gene expression profile (TIGEP) is more predictive of patient outcomes in the tumors with an increase in immune pressure over time, which reveals a curious synergy between T cells and neoantigen distributions. Overall, we provide a new tool that is capable of revealing the imprints left by neoantigens during each tumor's developmental process and of predicting how tumors will progress under further pressure of the host's immune system.

Inferring parameters of cancer evolution in chronic lymphocytic leukemia

As a cancer develops, its cells accrue new mutations, resulting in a heterogeneous, complex genomic profile. We make use of this heterogeneity to derive simple, analytic estimates of parameters driving carcinogenesis and reconstruct the timeline of selective events following initiation of an individual cancer, where two longitudinal samples are available for sequencing. Using stochastic computer simulations of cancer growth, we show that we can accurately estimate mutation rate, time before and after a driver event occurred, and growth rates of both initiated cancer cells and subsequently appearing subclones. We demonstrate that in order to obtain accurate estimates of mutation rate and timing of events, observed mutation counts should be corrected to account for clonal mutations that occurred after the founding of the tumor, as well as sequencing coverage. Chronic lymphocytic leukemia (CLL), which often does not require treatment for years after diagnosis, presents an optimal system to study the untreated, natural evolution of cancer cell populations. When we apply our methodology to reconstruct the individual evolutionary histories of CLL patients, we find that the parental leukemic clone typically appears within the first fifteen years of life.

Cancer Stem Cells (CSCs), Circulating Tumor Cells (CTCs) and Their Interplay with Cancer Associated Fibroblasts (CAFs): A New World of Targets and Treatments

Simple Summary

The world of small molecules in solid tumors as cancer stem cells (CSCs), circulating tumor cells (CTCs) and cancer-associated fibroblasts (CAFs) continues to be under-debated, but not of minor interest in recent decades. One of the main problems in regard to cancer is the development of tumor recurrence, even in the early stages, in addition to drug resistance and, consequently, ineffective or an incomplete response against the tumor. The findings behind this resistance are probably justified by the presence of small molecules such as CSCs, CTCs and CAFs connected with the tumor microenvironment, which may influence the aggressiveness and the metastatic process. The mechanisms, connections, and molecular pathways behind them are still unknown. Our review would like to represent an important step forward to highlight the roles of these molecules and the possible connections among them.

Abstract

The importance of defining new molecules to fight cancer is of significant interest to the scientific community. In particular, it has been shown that cancer stem cells (CSCs) are a small subpopulation of cells within tumors with capabilities of self-renewal, differentiation, and tumorigenicity; on the other side, circulating tumor cells (CTCs) seem to split away from the primary tumor and appear in the circulatory system as singular units or clusters. It is becoming more and more important to discover new biomarkers related to these populations of cells in combination to define the network among them and the tumor microenvironment. In particular, cancer-associated fibroblasts (CAFs) are a key component of the tumor microenvironment with different functions, including matrix deposition and remodeling, extensive reciprocal signaling interactions with cancer cells and crosstalk with immunity. The settings of new markers and the definition of the molecular connections may present new avenues, not only for fighting cancer but also for the definition of more tailored therapies.

Inferring ongoing cancer evolution from single tumour biopsies using synthetic supervised learning

Variant allele frequencies (VAF) encode ongoing evolution and subclonal selection in growing tumours. However, existing methods that utilize VAF information for cancer evolutionary inference are compressive, slow, or incorrectly specify the underlying cancer evolutionary dynamics. Here, we provide a proof-of-principle synthetic supervised learning method, TumE, that integrates simulated models of cancer evolution with Bayesian neural networks, to infer ongoing selection in bulk-sequenced single tumour biopsies. Analyses in synthetic and patient tumours show that TumE significantly improves both accuracy and inference time per sample when detecting positive selection, deconvoluting selected subclonal populations, and estimating subclone frequency. Importantly, we show how transfer learning can leverage stored knowledge within TumE models for related evolutionary inference tasks-substantially reducing data and computational time for further model development and providing a library of recyclable deep learning models for the cancer evolution community. This extensible framework provides a foundation and future directions for harnessing progressive computational methods for the benefit of cancer genomics and, in turn, the cancer patient.

Integrative analysis of genomic and transcriptomic data of normal, tumour and co-occurring leukoplakia tissue triads drawn from patients with gingivobuccal oral cancer identifies signatures of tumour initiation and progression

A thickened, white patch — leukoplakia — in the oral cavity is usually benign, but sometimes (in ~9% of individuals) it progresses to malignant tumour. Because the genomic basis of this progression is poorly understood, we undertook this study and collected samples of four tissues — leukoplakia, tumour, adjacent normal, and blood — from each of 28 patients suffering from gingivobuccal oral cancer. We performed multiomics analysis of the 112 collected tissues (four tissues per patient from 28 patients) and integrated information on progressive changes in the mutational and transcriptional profiles of each patient to create this genomic narrative. Additionally, we generated and analysed whole‐exome sequence data from leukoplakia tissues collected from 11 individuals not suffering from oral cancer. Nonsynonymous somatic mutations in the CASP8 gene were identified as the likely events to initiate malignant transformation, since these were frequently shared between tumour and co‐occurring leukoplakia. CASP8 alterations were also shown to enhance expressions of genes that favour lateral spread of mutant cells. During malignant transformation, additional pathogenic mutations are acquired in key genes (TP53, NOTCH1, HRAS) (41% of patients); chromosomal‐instability (arm‐level deletions of 19p and q, focal‐deletion of DNA‐repair pathway genes and NOTCH1, amplification of EGFR) (77%), and increased APOBEC‐activity (23%) are also observed. These additional alterations were present singly (18% of patients) or in combination (68%). Some of these alterations likely impact immune‐dynamics of the evolving transformed tissue; progression to malignancy is associated with immune suppression through infiltration of regulatory T‐cells (56%), depletion of cytotoxic T‐cells (68%), and antigen‐presenting dendritic cells (72%), with a concomitant increase in inflammation (92%). Patients can be grouped into three clusters by the estimated time to development of cancer from precancer by acquiring additional mutations (range: 4–10 years). Our findings provide deep molecular insights into the evolutionary processes and trajectories of oral cancer initiation and progression. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Computational Reconstruction of Clonal Hierarchies From Bulk Sequencing Data of Acute Myeloid Leukemia Samples

Acute myeloid leukemia is an aggressive cancer of the blood forming system. The malignant cell population is composed of multiple clones that evolve over time. Clonal data reflect the mechanisms governing treatment response and relapse. Single cell sequencing provides most direct insights into the clonal composition of the leukemic cells, however it is still not routinely available in clinical practice. In this work we develop a computational algorithm that allows identifying all clonal hierarchies that are compatible with bulk variant allele frequencies measured in a patient sample. The clonal hierarchies represent descendance relations between the different clones and reveal the order in which mutations have been acquired. The proposed computational approach is tested using single cell sequencing data that allow comparing the outcome of the algorithm with the true structure of the clonal hierarchy. We investigate which problems occur during reconstruction of clonal hierarchies from bulk sequencing data. Our results suggest that in many cases only a small number of possible hierarchies fits the bulk data. This implies that bulk sequencing data can be used to obtain insights in clonal evolution.

Topographic correlates of driver mutations and endogenous gene expression in pediatric diffuse midline gliomas and hemispheric high-grade gliomas

We evaluate the topographic distribution of diffuse midline gliomas and hemispheric high-grade gliomas in children with respect to their normal gene expression patterns and pathologic driver mutation patterns. We identified 19 pediatric patients with diffuse midline or high-grade glioma with preoperative MRI from tumor board review. 7 of these had 500 gene panel mutation testing, 11 patients had 50 gene panel mutation testing and one 343 gene panel testing from a separate institution were included as validation set. Tumor imaging features and gene expression patterns were analyzed using Allen Brain Atlas. Twelve patients had diffuse midline gliomas and seven had hemispheric high-grade gliomas. Three diffuse midline gliomas had the K27M mutation in the tail of histone H3 protein. All patients undergoing 500 gene panel testing had additional mutations, the most common being in ACVR1, PPM1D, and p53. Hemispheric high-grade gliomas had either TP53 or IDH1 mutation and diffuse midline gliomas had H3 K27M-mutation. Gene expression analysis in normal brains demonstrated that genes mutated in diffuse midline gliomas had higher expression along midline structures as compared to the cerebral hemispheres. Our study suggests that topographic location of pediatric diffuse midline gliomas and hemispheric high-grade gliomas correlates with driver mutations of tumor to the endogenous gene expression in that location. This correlation suggests that cellular state that is required for increased gene expression predisposes that location to mutations and defines the driver mutations within tumors that arise from that region.

Clonal fitness inferred from time-series modelling of single-cell cancer genomes

Progress in defining genomic fitness landscapes in cancer, especially those defined by copy number alterations (CNAs), has been impeded by lack of time-series single-cell sampling of polyclonal populations and temporal statistical models1-7. Here we generated 42,000 genomes from multi-year time-series single-cell whole-genome sequencing of breast epithelium and primary triple-negative breast cancer (TNBC) patient-derived xenografts (PDXs), revealing the nature of CNA-defined clonal fitness dynamics induced by TP53 mutation and cisplatin chemotherapy. Using a new Wright-Fisher population genetics model8,9 to infer clonal fitness, we found that TP53 mutation alters the fitness landscape, reproducibly distributing fitness over a larger number of clones associated with distinct CNAs. Furthermore, in TNBC PDX models with mutated TP53, inferred fitness coefficients from CNA-based genotypes accurately forecast experimentally enforced clonal competition dynamics. Drug treatment in three long-term serially passaged TNBC PDXs resulted in cisplatin-resistant clones emerging from low-fitness phylogenetic lineages in the untreated setting. Conversely, high-fitness clones from treatment-naive controls were eradicated, signalling an inversion of the fitness landscape. Finally, upon release of drug, selection pressure dynamics were reversed, indicating a fitness cost of treatment resistance. Together, our findings define clonal fitness linked to both CNA and therapeutic resistance in polyclonal tumours.

Clone Wars: Quantitatively Understanding Cancer Drug Resistance

A key aim of early clinical development for new cancer treatments is to detect the potential for efficacy early and to identify a safe therapeutic dose to take forward to phase II. Because of this need, researchers have sought to build mathematical models linking initial radiologic tumor response, often assessed after 6 to 8 weeks of treatment, with overall survival. However, there has been mixed success of this approach in the literature. We argue that evolutionary selection pressure should be considered to interpret these early efficacy signals and so optimize cancer therapy.

New Insights into Therapy-Induced Progression of Cancer

The malignant tumor is a complex heterogeneous set of cells functioning in a no less heterogeneous microenvironment. Like any dynamic system, cancerous tumors evolve and undergo changes in response to external influences, including therapy. Initially, most tumors are susceptible to treatment. However, remaining cancer cells may rapidly reestablish the tumor after a temporary remission. These new populations of malignant cells usually have increased resistance not only to the first-line agent, but also to the second- and third-line drugs, leading to a significant decrease in patient survival. Multiple studies describe the mechanism of acquired therapy resistance. In past decades, it became clear that, in addition to the simple selection of pre-existing resistant clones, therapy induces a highly complicated and tightly regulated molecular response that allows tumors to adapt to current and even subsequent therapeutic interventions. This review summarizes mechanisms of acquired resistance, such as secondary genetic alterations, impaired function of drug transporters, and autophagy. Moreover, we describe less obvious molecular aspects of therapy resistance in cancers, including epithelial-to-mesenchymal transition, cell cycle alterations, and the role of intercellular communication. Understanding these molecular mechanisms will be beneficial in finding novel therapeutic approaches for cancer therapy.

Multifunctional Supramolecular Assemblies with Aggregation-Induced Emission (AIE) for Cell Line Identification, Cell Contamination Evaluation, and Cancer Cell Discrimination

It is undoubted the important role of cells in biology and medicine, but worldwide misidentified and cross-contaminated cell lines have caused much trouble in related fields. Herein, three kinds of supramolecular AIE (aggregation-induced emission) nanoassemblies were constructed by the host-guest interaction between tetraphenylethene (TPE) derivatives and cucurbit[8]uril (CB[8]). Based on the recognized mechanism of AIE, the TPE derivatives could achieve stronger fluorescence emission and higher fluorescence quantum yield after assembling with CB[8]. Moreover, the constructed supramolecular AIE complexes obtained well-confirmed nanostructures and exhibited different sizes and shapes. Correspondingly, they generated characteristic biological properties and fluorescence enhancement of cells. Inspired by the concept of Big Data Analysis, these fluorescence signals were further transformed into a unique fingerprint of cells via linear discriminant analysis. Immediately, we realized the veracious identification between a normal cell line, two cancer cell lines, and two metastasized cancer cell lines in a qualitative analysis. More importantly, it was well used to monitor the evaluation of cross-contaminated cells and the discrimination of cancer cells. As a proper bioapplication of ideal supramolecular nanomaterials, this system was easy to learn and apply, and the whole procedure was kept to 20 min, without cell disruption, centrifugation, or washing steps.