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.

Chromosome 7 to the rescue: overcoming chromosome 10 loss in gliomas.

The co-occurrence of chromosome 10 loss and chromosome 7 gain in gliomas is the most frequent loss-gain co-aneuploidy pair in human cancers, a phenomenon that has been investigated without resolution since the late 1980s. Expanding beyond previous gene-centric studies, we investigate the co-occurrence in a genome-wide manner taking an evolutionary perspective. First, by mining large tumor aneuploidy data, we predict that the more likely order is 10 loss followed by 7 gain. Second, by analyzing extensive genomic and transcriptomic data from both patients and cell lines, we find that this co-occurrence can be explained by functional rescue interactions that are highly enriched on 7, which can possibly compensate for any detrimental consequences arising from the loss of 10. Finally, by analyzing transcriptomic data from normal, non-cancerous, human brain tissues, we provide a plausible reason why this co-occurrence happens preferentially in cancers originating in certain regions of the brain.

CKAP5 stabilizes CENP-E at kinetochores by regulating microtubule-chromosome attachments.

Stabilization of microtubule plus end-directed kinesin CENP-E at the metaphase kinetochores is important for chromosome alignment, but its mechanism remains unclear. Here, we show that CKAP5, a conserved microtubule plus tip protein, regulates CENP-E at kinetochores in human cells. Depletion of CKAP5 impairs CENP-E localization at kinetochores at the metaphase plate and results in increased kinetochore-microtubule stability and attachment errors. Erroneous attachments are also supported by computational modeling. Analysis of CKAP5 knockout cancer cells of multiple tissue origins shows that CKAP5 is preferentially essential in aneuploid, chromosomally unstable cells, and the sensitivity to CKAP5 depletion is correlated to that of CENP-E depletion. CKAP5 depletion leads to reduction in CENP-E-BubR1 interaction and the interaction is rescued by TOG4-TOG5 domain of CKAP5. The same domain can rescue CKAP5 depletion-induced CENP-E removal from the kinetochores. Interestingly, CKAP5 depletion facilitates recruitment of PP1 to the kinetochores and furthermore, a PP1 target site-specific CENP-E phospho-mimicking mutant gets stabilized at kinetochores in the CKAP5-depleted cells. Together, the results support a model in which CKAP5 controls mitotic chromosome attachment errors by stabilizing CENP-E at kinetochores and by regulating stability of the kinetochore-attached microtubules.

Acute expression of human APOBEC3B in mice results in RNA editing and lethality.

BACKGROUND: RNA editing has been described as promoting genetic heterogeneity, leading to the development of multiple disorders, including cancer. The cytosine deaminase APOBEC3B is implicated in tumor evolution through DNA mutation, but whether it also functions as an RNA editing enzyme has not been studied. RESULTS: Here, we engineer a novel doxycycline-inducible mouse model of human APOBEC3B-overexpression to understand the impact of this enzyme in tissue homeostasis and address a potential role in C-to-U RNA editing. Elevated and sustained levels of APOBEC3B lead to rapid alteration of cellular fitness, major organ dysfunction, and ultimately lethality in mice. Importantly, RNA-sequencing of mouse tissues expressing high levels of APOBEC3B identifies frequent UCC-to-UUC RNA editing events that are not evident in the corresponding genomic DNA. CONCLUSIONS: This work identifies, for the first time, a new deaminase-dependent function for APOBEC3B in RNA editing and presents a preclinical tool to help understand the emerging role of APOBEC3B as a driver of carcinogenesis.

Immunosignatures associated with TP53 status and co-mutations classify prognostically head and neck cancer patients.

BACKGROUND: Immune checkpoint inhibitors (ICIs) are a therapeutic strategy for various cancers although only a subset of patients respond to the therapy. Identifying patients more prone to respond to ICIs may increase the therapeutic benefit and allow studying new approaches for resistant patients. METHODS: We analyzed the TCGA cohort of HNSCC patients in relation to their activation of 26 immune gene expression signatures, as well as their cell type composition, in order to define signaling pathways associated with resistance to ICIs. Results were validated on two cohorts of 102 HNSCC patients and 139 HNSCC patients under treatment with PD-L1 inhibitors, respectively, and a cohort of 108 HNSCC HPV negative patients and by in vitro experiments in HNSCC cell lines. RESULTS: We observed a significant association between the gene set and TP53 gene status and OS and PFS of HNSCC patients. Surprisingly, the presence of a TP53 mutation together with another co-driver mutation was associated with significantly higher levels of the immune gene expression, in comparison to tumors in which the TP53 gene was mutated alone. In addition, the higher level of TP53 mutated-dependent MYC signature was associated with lower levels of the immune gene expression signature. In vitro and three different patient cohorts validation analyses corroborated these findings. CONCLUSIONS: Immune gene signature sets associated with TP53 status and co-mutations classify with more accuracy HNSCC patients. These biomarkers may be easily implemented in clinical setting.

Spectrum of Response to Platinum and PARP Inhibitors in Germline BRCA-Associated Pancreatic Cancer in the Clinical and Preclinical Setting.

Germline BRCA-associated pancreatic ductal adenocarcinoma (glBRCA PDAC) tumors are susceptible to platinum and PARP inhibition. The clinical outcomes of 125 patients with glBRCA PDAC were stratified based on the spectrum of response to platinum/PARP inhibition: (i) refractory [overall survival (OS) <6 months], (ii) durable response followed by acquired resistance (OS <36 months), and (iii) long-term responders (OS >36 months). Patient-derived xenografts (PDX) were generated from 25 patients with glBRCA PDAC at different clinical time points. Response to platinum/PARP inhibition in vivo and ex vivo culture (EVOC) correlated with clinical response. We deciphered the mechanisms of resistance in glBRCA PDAC and identified homologous recombination (HR) proficiency and secondary mutations restoring partial functionality as the most dominant resistant mechanism. Yet, a subset of HR-deficient (HRD) patients demonstrated clinical resistance. Their tumors displayed basal-like molecular subtype and were more aneuploid. Tumor mutational burden was high in HRD PDAC and significantly higher in tumors with secondary mutations. Anti-PD-1 attenuated tumor growth in a novel humanized glBRCA PDAC PDX model. This work demonstrates the utility of preclinical models, including EVOC, to predict the response of glBRCA PDAC to treatment, which has the potential to inform time-sensitive medical decisions. SIGNIFICANCE: glBRCA PDAC has a favorable response to platinum/PARP inhibition. However, most patients develop resistance. Additional treatment options for this unique subpopulation are needed. We generated model systems in PDXs and an ex vivo system (EVOC) that faithfully recapitulate these specific clinical scenarios as a platform to investigate the mechanisms of resistance for further drug development. This article is highlighted in the In This Issue feature, p. 1749.

Optimizing cancer immunotherapy response prediction by tumor aneuploidy score and fraction of copy number alterations.

Identifying patients that are likely to respond to cancer immunotherapy is an important, yet highly challenging clinical need. Using 3139 patients across 17 different cancer types, we comprehensively studied the ability of two common copy-number alteration (CNA) scores-the tumor aneuploidy score (AS) and the fraction of genome single nucleotide polymorphism encompassed by copy-number alterations (FGA)-to predict survival following immunotherapy in both pan-cancer and individual cancer types. First, we show that choice of cutoff during CNA calling significantly influences the predictive power of AS and FGA for patient survival following immunotherapy. Remarkably, by using proper cutoff during CNA calling, AS and FGA can predict pan-cancer survival following immunotherapy for both high-TMB and low-TMB patients. However, at the individual cancer level, our data suggest that the use of AS and FGA for predicting immunotherapy response is currently limited to only a few cancer types. Therefore, larger sample sizes are needed to evaluate the clinical utility of these measures for patient stratification in other cancer types. Finally, we propose a simple, non-parameterized, elbow-point-based method to help determine the cutoff used for calling CNAs.

Short-term molecular consequences of chromosome mis-segregation for genome stability.

Chromosome instability (CIN) is the most common form of genome instability and is a hallmark of cancer. CIN invariably leads to aneuploidy, a state of karyotype imbalance. Here, we show that aneuploidy can also trigger CIN. We found that aneuploid cells experience DNA replication stress in their first S-phase and precipitate in a state of continuous CIN. This generates a repertoire of genetically diverse cells with structural chromosomal abnormalities that can either continue proliferating or stop dividing. Cycling aneuploid cells display lower karyotype complexity compared to the arrested ones and increased expression of DNA repair signatures. Interestingly, the same signatures are upregulated in highly-proliferative cancer cells, which might enable them to proliferate despite the disadvantage conferred by aneuploidy-induced CIN. Altogether, our study reveals the short-term origins of CIN following aneuploidy and indicates the aneuploid state of cancer cells as a point mutation-independent source of genome instability, providing an explanation for aneuploidy occurrence in tumors.

A palmitate-rich metastatic niche enables metastasis growth via p65 acetylation resulting in pro-metastatic NF-κB signaling.

Metabolic rewiring is often considered an adaptive pressure limiting metastasis formation; however, some nutrients available at distant organs may inherently promote metastatic growth. We find that the lung and liver are lipid-rich environments. Moreover, we observe that pre-metastatic niche formation increases palmitate availability only in the lung, whereas a high-fat diet increases it in both organs. In line with this, targeting palmitate processing inhibits breast cancer-derived lung metastasis formation. Mechanistically, breast cancer cells use palmitate to synthesize acetyl-CoA in a carnitine palmitoyltransferase 1a-dependent manner. Concomitantly, lysine acetyltransferase 2a expression is promoted by palmitate, linking the available acetyl-CoA to the acetylation of the nuclear factor-kappaB subunit p65. Deletion of lysine acetyltransferase 2a or carnitine palmitoyltransferase 1a reduces metastasis formation in lean and high-fat diet mice, and lung and liver metastases from patients with breast cancer show coexpression of both proteins. In conclusion, palmitate-rich environments foster metastases growth by increasing p65 acetylation, resulting in a pro-metastatic nuclear factor-kappaB signaling.

Frequent aneuploidy in primary human T cells after CRISPR-Cas9 cleavage.

Multiple clinical trials of allogeneic T cell therapy use site-specific nucleases to disrupt T cell receptor (TCR) and other genes1-6. In this study, using single-cell RNA sequencing, we investigated genome editing outcomes in primary human T cells transfected with CRISPR-Cas9 and guide RNAs targeting genes for TCR chains and programmed cell death protein 1. Four days after transfection, we found a loss of chromosome 14, harboring the TCRα locus, in up to 9% of the cells and a chromosome 14 gain in up to 1.4% of the cells. Chromosome 7, harboring the TCRß locus, was truncated in 9.9% of the cells. Aberrations were validated using fluorescence in situ hybridization and digital droplet PCR. Aneuploidy was associated with reduced proliferation, induced p53 activation and cell death. However, at 11 days after transfection, 0.9% of T cells still had a chromosome 14 loss. Aneuploidy and chromosomal truncations are, thus, frequent outcomes of CRISPR-Cas9 cleavage that should be monitored and minimized in clinical protocols.

Whole-Genome Duplication Shapes the Aneuploidy Landscape of Human Cancers.

Aneuploidy is a hallmark of cancer with tissue-specific prevalence patterns that suggest it plays a driving role in cancer initiation and progression. However, the contribution of aneuploidy to tumorigenesis depends on both cellular and genomic contexts. Whole-genome duplication (WGD) is a common macroevolutionary event that occurs in more than 30% of human tumors early in tumorigenesis. Although tumors that have undergone WGD are reported to be more permissive to aneuploidy, it remains unknown whether WGD also affects aneuploidy prevalence patterns. Here we analyzed clinical tumor samples from 5,586 WGD- tumors and 3,435 WGD+ tumors across 22 tumor types and found distinct patterns of aneuploidy in WGD- and WGD+ tumors. WGD+ tumors were characterized by more promiscuous aneuploidy patterns, in line with increased aneuploidy tolerance. Moreover, the genetic interactions between chromosome arms differed between WGD- and WGD+ tumors, giving rise to distinct cooccurrence and mutual exclusivity aneuploidy patterns. The proportion of whole-chromosome aneuploidy compared with arm-level aneuploidy was significantly higher in WGD+ tumors, indicating distinct dominant mechanisms for aneuploidy formation. Human cancer cell lines successfully reproduced these WGD/aneuploidy interactions, confirming the relevance of studying this phenomenon in culture. Finally, induction of WGD and assessment of aneuploidy in isogenic WGD-/WGD+ human colon cancer cell lines under standard or selective conditions validated key findings from the clinical tumor analysis, supporting a causal link between WGD and altered aneuploidy landscapes. We conclude that WGD shapes the aneuploidy landscape of human tumors and propose that this interaction contributes to tumor evolution. SIGNIFICANCE: These findings suggest that the interactions between whole-genome duplication and aneuploidy are important for tumor evolution, highlighting the need to consider genome status in the analysis and modeling of cancer aneuploidy.

DNA-based copy number analysis confirms genomic evolution of PDX models.

Genomic evolution of patient-derived xenografts (PDXs) may lead to their gradual divergence away of their tumors of origin. We previously reported the genomic evolution of the copy number (CN) landscapes of PDXs during their engraftment and passaging1. However, whether PDX models are highly stable throughout passaging2, or can evolve CNAs rapidly1,3, remains controversial. Here, we reassess the genomic evolution of PDXs using DNA-based CN profiles. We find strong evidence for genomic evolution in the DNA-based PDX data: a median of ~10% of the genome is differentially altered between matched primary tumors (PTs) and PDXs across cohorts (range, 0% to 73% across all models). In 24% of the matched PT-PDX samples, over a quarter of the genome is differentially affected by CN alterations. Moreover, in matched analyses of PTs and their derived PDXs at multiple passages, later-passage PDXs are significantly less similar to their parental PTs than earlier-passage PDXs, indicative of genomic divergence. We conclude that PDX models indeed evolve throughout their derivation and propagation, and that the phenotypic consequences of this evolution ought to be assessed in order to determine its relevance to the proper application of these valuable cancer models.

Gene copy-number changes and chromosomal instability induced by aneuploidy confer resistance to chemotherapy.

Mitotic errors lead to aneuploidy, a condition of karyotype imbalance, frequently found in cancer cells. Alterations in chromosome copy number induce a wide variety of cellular stresses, including genome instability. Here, we show that cancer cells might exploit aneuploidy-induced genome instability and the resulting gene copy-number changes to survive under conditions of selective pressure, such as chemotherapy. Resistance to chemotherapeutic drugs was dictated by the acquisition of recurrent karyotypes, indicating that gene dosage might play a role in driving chemoresistance. Thus, our study establishes a causal link between aneuploidy-driven changes in gene copy number and chemoresistance and might explain why some chemotherapies fail to succeed.

Aneuploid senescent cells activate NF-κB to promote their immune clearance by NK cells.

The immune system plays a major role in the protection against cancer. Identifying and characterizing the pathways mediating this immune surveillance are thus critical for understanding how cancer cells are recognized and eliminated. Aneuploidy is a hallmark of cancer, and we previously found that untransformed cells that had undergone senescence due to highly abnormal karyotypes are eliminated by natural killer (NK) cells in vitro. However, the mechanisms underlying this process remained elusive. Here, using an in vitro NK cell killing system, we show that non-cell-autonomous mechanisms in aneuploid cells predominantly mediate their clearance by NK cells. Our data indicate that in untransformed aneuploid cells, NF-κB signaling upregulation is central to elicit this immune response. Inactivating NF-κB abolishes NK cell-mediated clearance of untransformed aneuploid cells. In cancer cell lines, NF-κB upregulation also correlates with the degree of aneuploidy. However, such upregulation in cancer cells is not sufficient to trigger NK cell-mediated clearance, suggesting that additional mechanisms might be at play during cancer evolution to counteract NF-κB-mediated immunogenicity.

Relevance of aneuploidy for cancer therapies targeting the spindle assembly checkpoint and KIF18A.

Aneuploidy, a common feature of cancer cells, results in increased sensitivity to the inhibition of the spindle assembly checkpoint (SAC) and the mitotic motor protein Kinesin Family Member 18A (KIF18A). We discuss the importance of drugs targeting SAC core members and KIF18A. We stress the need to assess the sensitivity to this class of drugs at appropriate time points, and propose that aneuploidy could serve as a biomarker to stratify patients for SAC-targeting treatments.

FOXM1 repression increases mitotic death upon antimitotic chemotherapy through BMF upregulation.

Inhibition of spindle microtubule (MT) dynamics has been effectively used in cancer treatment. Although the mechanisms by which MT poisons elicit mitotic arrest are fairly understood, efforts are still needed towards elucidating how cancer cells respond to antimitotic drugs owing to cytotoxicity and resistance side effects. Here, we identified the critical G2/M transcription factor Forkhead box M1 (FOXM1) as a molecular determinant of cell response to antimitotics. We found FOXM1 repression to increase death in mitosis (DiM) due to upregulation of the BCL-2 modifying factor (BMF) gene involved in anoikis, an apoptotic process induced upon cell detachment from the extracellular matrix. FOXM1 binds to a BMF intronic cis-regulatory element that interacts with both the BMF and the neighbor gene BUB1B promoter regions, to oppositely regulate their expression. This mechanism ensures that cells treated with antimitotics repress BMF and avoid DiM when FOXM1 levels are high. In addition, we show that this mechanism is partly disrupted in anoikis/antimitotics-resistant tumor cells, with resistance correlating with lower BMF expression but in a FOXM1-independent manner. These findings provide a stratification biomarker for antimitotic chemotherapy response.

Aneuploidy renders cancer cells vulnerable to mitotic checkpoint inhibition.

Selective targeting of aneuploid cells is an attractive strategy for cancer treatment1. However, it is unclear whether aneuploidy generates any clinically relevant vulnerabilities in cancer cells. Here we mapped the aneuploidy landscapes of about 1,000 human cancer cell lines, and analysed genetic and chemical perturbation screens2-9 to identify cellular vulnerabilities associated with aneuploidy. We found that aneuploid cancer cells show increased sensitivity to genetic perturbation of core components of the spindle assembly checkpoint (SAC), which ensures the proper segregation of chromosomes during mitosis10. Unexpectedly, we also found that aneuploid cancer cells were less sensitive than diploid cells to short-term exposure to multiple SAC inhibitors. Indeed, aneuploid cancer cells became increasingly sensitive to inhibition of SAC over time. Aneuploid cells exhibited aberrant spindle geometry and dynamics, and kept dividing when the SAC was inhibited, resulting in the accumulation of mitotic defects, and in unstable and less-fit karyotypes. Therefore, although aneuploid cancer cells could overcome inhibition of SAC more readily than diploid cells, their long-term proliferation was jeopardized. We identified a specific mitotic kinesin, KIF18A, whose activity was perturbed in aneuploid cancer cells. Aneuploid cancer cells were particularly vulnerable to depletion of KIF18A, and KIF18A overexpression restored their response to SAC inhibition. Our results identify a therapeutically relevant, synthetic lethal interaction between aneuploidy and the SAC.