Increased RNA and protein degradation is required for counteracting transcriptional burden and proteotoxic stress in human aneuploid cells.

Aneuploidy results in a stoichiometric imbalance of protein complexes that jeopardizes cellular fitness. Aneuploid cells thus need to compensate for the imbalanced DNA levels by regulating their RNA and protein levels, but the underlying molecular mechanisms remain unknown. Here, we dissected multiple diploid vs. aneuploid cell models. We found that aneuploid cells cope with transcriptional burden by increasing several RNA degradation pathways, and are consequently more sensitive to the perturbation of RNA degradation. At the protein level, aneuploid cells mitigate proteotoxic stress by reducing protein translation and increasing protein degradation, rendering them more sensitive to proteasome inhibition. These findings were recapitulated across hundreds of human cancer cell lines and primary tumors, and aneuploidy levels were significantly associated with the response of multiple myeloma patients to proteasome inhibitors. Aneuploid cells are therefore preferentially dependent on several key nodes along the gene expression process, creating clinically-actionable vulnerabilities in aneuploid cells.

Human aneuploid cells depend on the RAF/MEK/ERK pathway for overcoming increased DNA damage.

Aneuploidy is a hallmark of human cancer, yet the molecular mechanisms to cope with aneuploidy-induced cellular stresses remain largely unknown. Here, we induce chromosome mis-segregation in non-transformed RPE1-hTERT cells and derive multiple stable clones with various degrees of aneuploidy. We perform a systematic genomic, transcriptomic and proteomic profiling of 6 isogenic clones, using whole-exome DNA, mRNA and miRNA sequencing, as well as proteomics. Concomitantly, we functionally interrogate their cellular vulnerabilities, using genome-wide CRISPR/Cas9 and large-scale drug screens. Aneuploid clones activate the DNA damage response and are more resistant to further DNA damage induction. Aneuploid cells also exhibit elevated RAF/MEK/ERK pathway activity and are more sensitive to clinically-relevant drugs targeting this pathway, and in particular to CRAF inhibition. Importantly, CRAF and MEK inhibition sensitize aneuploid cells to DNA damage-inducing chemotherapies and to PARP inhibitors. We validate these results in human cancer cell lines. Moreover, resistance of cancer patients to olaparib is associated with high levels of RAF/MEK/ERK signaling, specifically in highly-aneuploid tumors. Overall, our study provides a comprehensive resource for genetically-matched karyotypically-stable cells of various aneuploidy states, and reveals a therapeutically-relevant cellular dependency of aneuploid cells.

A p62-dependent rheostat dictates micronuclei catastrophe and chromosome rearrangements.

Chromosomal instability (CIN) generates micronuclei-aberrant extranuclear structures that catalyze the acquisition of complex chromosomal rearrangements present in cancer. Micronuclei are characterized by persistent DNA damage and catastrophic nuclear envelope collapse, which exposes DNA to the cytoplasm. We found that the autophagic receptor p62/SQSTM1 modulates micronuclear stability, influencing chromosome fragmentation and rearrangements. Mechanistically, proximity of micronuclei to mitochondria led to oxidation-driven homo-oligomerization of p62, limiting endosomal sorting complex required for transport (ESCRT)-dependent micronuclear envelope repair by triggering autophagic degradation. We also found that p62 levels correlate with increased chromothripsis across human cancer cell lines and with increased CIN in colorectal tumors. Thus, p62 acts as a regulator of micronuclei and may serve as a prognostic marker for tumors with high CIN.

Chromosome 7 gain compensates for chromosome 10 loss in glioma.

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. This phenomenon has been investigated since the late 1980s without resolution. Expanding beyond previous gene-centric studies, we investigated the co-occurrence in a genome-wide manner taking an evolutionary perspective. Mining of large-scale tumor aneuploidy data confirmed the previous finding of a small-scale longitudinal study that the most likely order is chromosome 10 loss followed by chromosome 7 gain. Extensive analysis of genomic and transcriptomic data from both patients and cell lines revealed that this co-occurrence can be explained by functional rescue interactions that are highly enriched on chromosome 7, which could potentially compensate for any detrimental consequences arising from the loss of chromosome 10. Transcriptomic data from various normal, non-cancerous human brain tissues were analyzed to assess which tissues may be most predisposed to tolerate compensation of chromosome 10 loss by chromosome 7 gain. The analysis indicated that the pre-existing transcriptomic states in the cortex and frontal cortex, where gliomas arise, are more favorable than other brain regions for compensation by rescuer genes that are active on chromosome 7. Collectively, these findings suggest that the phenomenon of chromosome 10 loss and chromosome 7 gain in gliomas is orchestrated by a complex interaction of many genes residing within these two chromosomes and provide a plausible reason why this co-occurrence happens preferentially in cancers originating in certain regions of the brain.

Feeder-free culture of human pluripotent stem cells drives MDM4-mediated gain of chromosome 1q.

Culture-acquired variants in human pluripotent stem cells (hPSCs) hinder their applications in research and clinic. However, the mechanisms that underpin selection of variants remain unclear. Here, through analysis of comprehensive karyotyping datasets from over 23,000 hPSC cultures of more than 1,500 lines, we explored how culture conditions shape variant selection. Strikingly, we identified an association of chromosome 1q gains with feeder-free cultures and noted a rise in its prevalence in recent years, coinciding with increased usage of feeder-free regimens. Competition experiments of multiple isogenic lines with and without a chromosome 1q gain confirmed that 1q variants have an advantage in feeder-free (E8/vitronectin), but not feeder-based, culture. Mechanistically, we show that overexpression of MDM4, located on chromosome 1q, drives variants' advantage in E8/vitronectin by alleviating genome damage-induced apoptosis, which is lower in feeder-based conditions. Our study explains condition-dependent patterns of hPSC aberrations and offers insights into the mechanisms of variant selection.

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.