A Markov chain for numerical chromosomal instability in clonally expanding populations

Inference of chromosome selection parameters and missegregation rate in cancer from DNA-sequencing data

Aneuploidy is frequently observed in cancers and has been linked to poor patient outcome. Analysis of aneuploidy in DNA-sequencing (DNA-seq) data necessitates untangling the effects of the Copy Number Aberration (CNA) occurrence rates and the selection coefficients that act upon the resulting karyotypes. We introduce a parameter inference algorithm that takes advantage of both bulk and single-cell DNA-seq cohorts. The method is based on Approximate Bayesian Computation (ABC) and utilizes CINner, our recently introduced simulation algorithm of chromosomal instability in cancer. We examine three groups of statistics to summarize the data in the ABC routine: (A) Copy Number-based measures, (B) phylogeny tip statistics, and (C) phylogeny balance indices. Using these statistics, our method can recover both the CNA probabilities and selection parameters from ground truth data, and performs well even for data cohorts of relatively small sizes. We find that only statistics in groups A and C are well-suited for identifying CNA probabilities, and only group A carries the signals for estimating selection parameters. Moreover, the low number of CNA events at large scale compared to cell counts in single-cell samples means that statistics in group B cannot be estimated accurately using phylogeny reconstruction algorithms at the chromosome level. As data from both bulk and single-cell DNA-sequencing techniques becomes increasingly available, our inference framework promises to facilitate the analysis of distinct cancer types, differentiation between selection and neutral drift, and prediction of cancer clonal dynamics.

YY2/BUB3 Axis promotes SAC Hyperactivation and Inhibits Colorectal Cancer Progression via Regulating Chromosomal Instability

Spindle assembly checkpoint (SAC) is a crucial safeguard mechanism of mitosis fidelity that ensures equal division of duplicated chromosomes to the two progeny cells. Impaired SAC can lead to chromosomal instability (CIN), a well-recognized hallmark of cancer that facilitates tumor progression; paradoxically, high CIN levels are associated with better therapeutic response and prognosis. However, the mechanism by which CIN determines tumor cell survival and therapeutic response remains poorly understood. Here, using a cross-omics approach, YY2 is identified as a mitotic regulator that promotes SAC activity by activating the transcription of budding uninhibited by benzimidazole 3 (BUB3), a component of SAC. While both conditions induce CIN, a defect in YY2/SAC activity enhances mitosis and tumor growth. Meanwhile, hyperactivation of SAC mediated by YY2/BUB3 triggers a delay in mitosis and suppresses growth. Furthermore, it is revealed that YY2/BUB3-mediated excessive CIN causes higher cell death rates and drug sensitivity, whereas residual tumor cells that survived DNA damage-based therapy have moderate CIN and increased drug resistance. These results provide insights into the role of SAC activity and CIN levels in influencing tumor cell survival and drug response, as well as suggest a novel anti-tumor therapeutic strategy that combines SAC activity modulators and DNA-damage agents.

CINner: modeling and simulation of chromosomal instability in cancer at single-cell resolution

Cancer development is characterized by chromosomal instability, manifesting in frequent occurrences of different genomic alteration mechanisms ranging in extent and impact. Mathematical modeling can help evaluate the role of each mutational process during tumor progression, however existing frameworks can only capture certain aspects of chromosomal instability (CIN). We present CINner, a mathematical framework for modeling genomic diversity and selection during tumor evolution. The main advantage of CINner is its flexibility to incorporate many genomic events that directly impact cellular fitness, from driver gene mutations to copy number alterations (CNAs), including focal amplifications and deletions, missegregations and whole-genome duplication (WGD). We apply CINner to find chromosome-arm selection parameters that drive tumorigenesis in the absence of WGD in chromosomally stable cancer types. We found that the selection parameters predict WGD prevalence among different chromosomally unstable tumors, hinting that the selective advantage of WGD cells hinges on their tolerance for aneuploidy and escape from nullisomy. Direct application of CINner to model the WGD proportion and fraction of genome altered (FGA) further uncovers the increase in CNA probabilities associated with WGD in each cancer type. CINner can also be utilized to study chromosomally stable cancer types, by applying a selection model based on driver gene mutations and focal amplifications or deletions. Finally, we used CINner to analyze the impact of CNA probabilities, chromosome selection parameters, tumor growth dynamics and population size on cancer fitness and heterogeneity. We expect that CINner will provide a powerful modeling tool for the oncology community to quantify the impact of newly uncovered genomic alteration mechanisms on shaping tumor progression and adaptation.

The reckoning of chromosomal instability: past, present, future

Quantitative measures of CIN are crucial to our understanding of its role in cancer. Technological advances have changed the way CIN is quantified, offering increased accuracy and insight. Here, we review measures of CIN through its rise as a field, discuss considerations for its measurement, and look forward to future quantification of CIN.

CNETML: maximum likelihood inference of phylogeny from copy number profiles of multiple samples

Phylogenetic trees based on copy number profiles from multiple samples of a patient are helpful to understand cancer evolution. Here, we develop a new maximum likelihood method, CNETML, to infer phylogenies from such data. CNETML is the first program to jointly infer the tree topology, node ages, and mutation rates from total copy numbers of longitudinal samples. Our extensive simulations suggest CNETML performs well on copy numbers relative to ploidy and under slight violation of model assumptions. The application of CNETML to real data generates results consistent with previous discoveries and provides novel early copy number events for further investigation.

Targeting the undruggable

Journey through basic biology reveals a way to treat chromosomally unstable cancers.

Proliferative advantage of specific aneuploid cells drives evolution of tumor karyotypes

Most tumors have abnormal karyotypes, which arise from mistakes during mitotic division of healthy euploid cells and evolve through numerous complex mechanisms. In a recent mouse model with increased chromosome missegregation, chromosome gains dominate over losses both in pretumor and tumor tissues, whereas T-cell lymphomas are characterized by gains of chromosomes 14 and 15. However, the quantitative understanding of clonal selection leading to tumor karyotype evolution remains unknown. Here we show, by introducing a mathematical model based on a concept of a macro-karyotype, that tumor karyotypes can be explained by proliferation-driven evolution of aneuploid cells. In pretumor cells, increased apoptosis and slower proliferation of cells with monosomies lead to predominant chromosome gains over losses. Tumor karyotypes with gain of one chromosome can be explained by karyotype-dependent proliferation, whereas, for those with two chromosomes, an interplay with karyotype-dependent apoptosis is an additional possible pathway. Thus, evolution of tumor-specific karyotypes requires proliferative advantage of specific aneuploid karyotypes.

Intra-tumor heterogeneity, turnover rate and karyotype space shape susceptibility to missegregation-induced extinction

The phenotypic efficacy of somatic copy number alterations (SCNAs) stems from their incidence per base pair of the genome, which is orders of magnitudes greater than that of point mutations. One mitotic event stands out in its potential to significantly change a cell's SCNA burden-a chromosome missegregation. A stochastic model of chromosome mis-segregations has been previously developed to describe the evolution of SCNAs of a single chromosome type. Building upon this work, we derive a general deterministic framework for modeling missegregations of multiple chromosome types. The framework offers flexibility to model intra-tumor heterogeneity in the SCNAs of all chromosomes, as well as in missegregation- and turnover rates. The model can be used to test how selection acts upon coexisting karyotypes over hundreds of generations. We use the model to calculate missegregation-induced population extinction (MIE) curves, that separate viable from non-viable populations as a function of their turnover- and missegregation rates. Turnover- and missegregation rates estimated from scRNA-seq data are then compared to theoretical predictions. We find convergence of theoretical and empirical results in both the location of MIE curves and the necessary conditions for MIE. When a dependency of missegregation rate on karyotype is introduced, karyotypes associated with low missegregation rates act as a stabilizing refuge, rendering MIE impossible unless turnover rates are exceedingly high. Intra-tumor heterogeneity, including heterogeneity in missegregation rates, increases as tumors progress, rendering MIE unlikely.

GPER1 links estrogens to centrosome amplification and chromosomal instability in human colon cells

The role of the alternate G protein-coupled estrogen receptor 1 (GPER1) in colorectal cancer (CRC) development and progression is unclear, not least because of conflicting clinical and experimental evidence for pro- and anti-tumorigenic activities. Here, we show that low concentrations of the estrogenic GPER1 ligands, 17β-estradiol, bisphenol A, and diethylstilbestrol cause the generation of lagging chromosomes in normal colon and CRC cell lines, which manifest in whole chromosomal instability and aneuploidy. Mechanistically, (xeno)estrogens triggered centrosome amplification by inducing centriole overduplication that leads to transient multipolar mitotic spindles, chromosome alignment defects, and mitotic laggards. Remarkably, we could demonstrate a significant role of estrogen-activated GPER1 in centrosome amplification and increased karyotype variability. Indeed, both gene-specific knockdown and inhibition of GPER1 effectively restored normal centrosome numbers and karyotype stability in cells exposed to 17β-estradiol, bisphenol A, or diethylstilbestrol. Thus, our results reveal a novel link between estrogen-activated GPER1 and the induction of key CRC-prone lesions, supporting a pivotal role of the alternate estrogen receptor in colon neoplastic transformation and tumor progression.

MEDICC2: whole-genome doubling aware copy-number phylogenies for cancer evolution

Aneuploidy, chromosomal instability, somatic copy-number alterations, and whole-genome doubling (WGD) play key roles in cancer evolution and provide information for the complex task of phylogenetic inference. We present MEDICC2, a method for inferring evolutionary trees and WGD using haplotype-specific somatic copy-number alterations from single-cell or bulk data. MEDICC2 eschews simplifications such as the infinite sites assumption, allowing multiple mutations and parallel evolution, and does not treat adjacent loci as independent, allowing overlapping copy-number events. Using simulations and multiple data types from 2780 tumors, we use MEDICC2 to demonstrate accurate inference of phylogenies, clonal and subclonal WGD, and ancestral copy-number states.

Quantifying chromosomal instability from intratumoral karyotype diversity using agent-based modeling and Bayesian inference

Chromosomal instability (CIN)-persistent chromosome gain or loss through abnormal mitotic segregation-is a hallmark of cancer that drives aneuploidy. Intrinsic chromosome mis-segregation rate, a measure of CIN, can inform prognosis and is a promising biomarker for response to anti-microtubule agents. However, existing methodologies to measure this rate are labor intensive, indirect, and confounded by selection against aneuploid cells, which reduces observable diversity. We developed a framework to measure CIN, accounting for karyotype selection, using simulations with various levels of CIN and models of selection. To identify the model parameters that best fit karyotype data from single-cell sequencing, we used approximate Bayesian computation to infer mis-segregation rates and karyotype selection. Experimental validation confirmed the extensive chromosome mis-segregation rates caused by the chemotherapy paclitaxel (18.5 ± 0.5/division). Extending this approach to clinical samples revealed that inferred rates fell within direct observations of cancer cell lines. This work provides the necessary framework to quantify CIN in human tumors and develop it as a predictive biomarker.

In silico Prediction on the PI3K/AKT/mTOR Pathway of the Antiproliferative Effect of O. joconostle in Breast Cancer Models

The search for new cancer treatments from traditional medicine involves developing studies to understand at the molecular level different cell signaling pathways involved in cancer development. In this work, we present a model of the PI3K/Akt/mTOR pathway, which plays a key role in cell cycle regulation and is related to cell survival, proliferation, and growth in cancer, as well as resistance to antitumor therapies, so finding drugs that act on this pathway is ideal to propose a new adjuvant treatment. The aim of this work was to model, simulate and predict in silico using the Big Data-Cellulat platform the possible targets in the PI3K/Akt/mTOR pathway on which the Opuntia joconostle extract acts, as well as to indicate the concentration range to be used to find the mean lethal dose in in vitro experiments on breast cancer cells. The in silico results show that, in a cancer cell, the activation of JAK and STAT, as well as PI3K and Akt is related to the effect of cell proliferation, angiogenesis, and inhibition of apoptosis, and that the extract of O. joconostle has an antiproliferative effect on breast cancer cells by inhibiting cell proliferation, regulating the cell cycle and inhibiting apoptosis through this signaling pathway . In vitro it was demonstrated that the extract shows an antiproliferative effect, causing the arrest of cells in the G2/M phase of the cell cycle. Therefore, it is concluded that the use of in silico tools is a valuable method to perform virtual experiments and discover new treatments. The use of this type of model supports in vitro experimentation, reducing the costs and number of experiments in the real laboratory.

Karyotype Aberrations in Action: The Evolution of Cancer Genomes and the Tumor Microenvironment

Cancer is a disease of cellular evolution. For this cellular evolution to take place, a population of cells must contain functional heterogeneity and an assessment of this heterogeneity in the form of natural selection. Cancer cells from advanced malignancies are genomically and functionally very different compared to the healthy cells from which they evolved. Genomic alterations include aneuploidy (numerical and structural changes in chromosome content) and polyploidy (e.g., whole genome doubling), which can have considerable effects on cell physiology and phenotype. Likewise, conditions in the tumor microenvironment are spatially heterogeneous and vastly different than in healthy tissues, resulting in a number of environmental niches that play important roles in driving the evolution of tumor cells. While a number of studies have documented abnormal conditions of the tumor microenvironment and the cellular consequences of aneuploidy and polyploidy, a thorough overview of the interplay between karyotypically abnormal cells and the tissue and tumor microenvironments is not available. Here, we examine the evidence for how this interaction may unfold during tumor evolution. We describe a bidirectional interplay in which aneuploid and polyploid cells alter and shape the microenvironment in which they and their progeny reside; in turn, this microenvironment modulates the rate of genesis for new karyotype aberrations and selects for cells that are most fit under a given condition. We conclude by discussing the importance of this interaction for tumor evolution and the possibility of leveraging our understanding of this interplay for cancer therapy.

Integrating Mathematical Modeling with High-Throughput Imaging Explains How Polyploid Populations Behave in Nutrient-Sparse Environments

Breast cancer progresses in a multistep process from primary tumor growth and stroma invasion to metastasis. Nutrient-limiting environments promote chemotaxis with aggressive morphologies characteristic of invasion. It is unknown how coexisting cells differ in their response to nutrient limitations and how this impacts invasion of the metapopulation as a whole. In this study, we integrate mathematical modeling with microenvironmental perturbation data to investigate invasion in nutrient-limiting environments inhabited by one or two cancer cell subpopulations. Subpopulations were defined by their energy efficiency and chemotactic ability. Invasion distance traveled by a homogeneous population was estimated. For heterogeneous populations, results suggest that an imbalance between nutrient efficacy and chemotactic superiority accelerates invasion. Such imbalance will spatially segregate the two populations and only one type will dominate at the invasion front. Only if these two phenotypes are balanced, the two subpopulations compete for the same space, which decelerates invasion. We investigate ploidy as a candidate biomarker of this phenotypic heterogeneity and discuss its potential to inform the dose of mTOR inhibitors (mTOR-I) that can inhibit chemotaxis just enough to facilitate such competition. SIGNIFICANCE: This study identifies the double-edged sword of high ploidy as a prerequisite to personalize combination therapies with cytotoxic drugs and inhibitors of signal transduction pathways such as mTOR-Is. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/22/5109/F1.large.jpg.

Pervasive chromosomal instability and karyotype order in tumour evolution

Chromosomal instability in cancer consists of dynamic changes to the number and structure of chromosomes1,2. The resulting diversity in somatic copy number alterations (SCNAs) may provide the variation necessary for tumour evolution1,3,4. Here we use multi-sample phasing and SCNA analysis of 1,421 samples from 394 tumours across 22 tumour types to show that continuous chromosomal instability results in pervasive SCNA heterogeneity. Parallel evolutionary events, which cause disruption in the same genes (such as BCL9, MCL1, ARNT (also known as HIF1B), TERT and MYC) within separate subclones, were present in 37% of tumours. Most recurrent losses probably occurred before whole-genome doubling, that was found as a clonal event in 49% of tumours. However, loss of heterozygosity at the human leukocyte antigen (HLA) locus and loss of chromosome 8p to a single haploid copy recurred at substantial subclonal frequencies, even in tumours with whole-genome doubling, indicating ongoing karyotype remodelling. Focal amplifications that affected chromosomes 1q21 (which encompasses BCL9, MCL1 and ARNT), 5p15.33 (TERT), 11q13.3 (CCND1), 19q12 (CCNE1) and 8q24.1 (MYC) were frequently subclonal yet appeared to be clonal within single samples. Analysis of an independent series of 1,024 metastatic samples revealed that 13 focal SCNAs were enriched in metastatic samples, including gains in chromosome 8q24.1 (encompassing MYC) in clear cell renal cell carcinoma and chromosome 11q13.3 (encompassing CCND1) in HER2+ breast cancer. Chromosomal instability may enable the continuous selection of SCNAs, which are established as ordered events that often occur in parallel, throughout tumour evolution.

Short-term in vitro glutamine restriction differentially impacts the chromosomal stability of transformed and non-transformed cells

Glutamine (Gln) is a non-essential amino acid central for generating building blocks and cellular energy in tumours and rapidly proliferating non-transformed cells. However, the influence of Gln on regulating chromosomal stability of transformed and non-transformed cells remain poorly understand. We hypothesised that Gln is required for maintaining a homeostatic level of chromosomal stability. To this end, transformed cells HeLa and A375 and non-transformed cells NCM460 and HUVEC cells were intervened with varying concentrations of Gln (10, 1, 0.1 and 0.01 mM), with or without cisplatin (0.1 µg/ml), for 24 h. The cytokinesis-block micronucleus (MN) assay was used to determine chromosomal instability (CIN), the extent of which is reflected by the frequency of MN, nucleoplasmic bridge (NPB) and nuclear bud (NB). We demonstrated an unexpected decrease in the spontaneous rate of MN, but not NPB and NB, after Gln restriction in HeLa and A375 cells. Gln restriction reduced cisplatin-induced MN, but not NPB and NB, in HeLa and A375 cells. We further revealed that Gln restriction suppressed the proliferation of HeLa cells with high CIN induced by nocodazole, partially explaining why Gln restriction decreased the frequency of spontaneous and cisplatin-induced MN in transformed cells. In contrast, Gln restriction increased MN and NB, but not NPB, in NCM460 cells. In HUVEC cells, Gln restriction increased MN, NPB and NB. Meanwhile, Gln restriction sensitised NCM460 cells to cisplatin-induced genotoxicity. A similar but more pronounced pattern was observed in HUVEC cells. Collectively, these results suggest that the in vitro influences of Gln metabolism on CIN depend on cellular contexts: Transformed cells require high Gln to fine tune their CIN in an optimal rate to maximise genomic heterogeneity and fitness, whereas non-transformed cells need high Gln to prevent CIN.