Profiling Numerical and Structural Chromosomal Instability in Different Cancer Types.

Many cancers display whole chromosome instability (W-CIN) and structural chromosomal instability (S-CIN), referring to increased rates of acquiring numerically and structurally abnormal chromosome changes. This protocol provides detailed steps to analyze the W-CIN and S-CIN across cancer types, intending to leverage large-scale bulk sequencing and SNP array data complemented with the computational models to gain a better understanding of W-CIN and S-CIN.

Increased replication origin firing links replication stress to whole chromosomal instability in human cancer.

Chromosomal instability (CIN) is a hallmark of cancer and comprises structural CIN (S-CIN) and numerical or whole chromosomal CIN (W-CIN). Recent work indicated that replication stress (RS), known to contribute to S-CIN, also affects mitotic chromosome segregation, possibly explaining the common co-existence of S-CIN and W-CIN in human cancer. Here, we show that RS-induced increased origin firing is sufficient to trigger W-CIN in human cancer cells. We discovered that overexpression of origin firing genes, including GINS1 and CDC45, correlates with W-CIN in human cancer specimens and causes W-CIN in otherwise chromosomally stable human cells. Furthermore, modulation of the ATR-CDK1-RIF1 axis increases the number of firing origins and leads to W-CIN. Importantly, chromosome missegregation upon additional origin firing is mediated by increased mitotic microtubule growth rates, a mitotic defect prevalent in chromosomally unstable cancer cells. Thus, our study identifies increased replication origin firing as a cancer-relevant trigger for chromosomal instability.

Stemness Correlates Inversely with MHC Class I Expression in Pediatric Small Round Blue Cell Tumors.

Recently, immunotherapeutic approaches have become a feasible option for a subset of pediatric cancer patients. Low MHC class I expression hampers the use of immunotherapies relying on antigen presentation. A well-established stemness score (mRNAsi) was determined using the bulk transcriptomes of 1134 pediatric small round blue cell tumors. Interestingly, MHC class I gene expression (HLA-A/-B/-C) was correlated negatively with mRNAsi throughout all diagnostic entities: neuroblastomas (NB) (n = 88, r = −0.41, p < 0.001), the Ewing's sarcoma family of tumors (ESFT) (n = 117, r = −0.46, p < 0.001), rhabdomyosarcomas (RMS) (n = 158, r = −0.5, p < 0.001), Wilms tumors (WT) (n = 224, r = −0.39, p < 0.001), and central nervous system-primitive neuroectodermal tumors CNS-PNET (r = −0.49, p < 0.001), with the exception of medulloblastoma (MB) (n = 76, r = −0.24, p = 0.06). The negative correlation of MHC class I and mRNAsi was independent of clinical features in NB, RMS, and WT. In NB and WT, increased MHC class I was correlated negatively with tumor stage. RMS patients with a high expression of MHC class I and abundant CD8 T cells showed a prolonged overall survival (n = 148, p = 0.004). Possibly, low MHC class I expression and stemness in pediatric tumors are remnants of prenatal tumorigenesis from multipotent precursor cells. Further studies are needed to assess the usefulness of stemness and MHC class I as predictive markers.

Distinct and Common Features of Numerical and Structural Chromosomal Instability across Different Cancer Types.

A large proportion of tumours is characterised by numerical or structural chromosomal instability (CIN), defined as an increased rate of gaining or losing whole chromosomes (W-CIN) or of accumulating structural aberrations (S-CIN). Both W-CIN and S-CIN are associated with tumourigenesis, cancer progression, treatment resistance and clinical outcome. Although W-CIN and S-CIN can co-occur, they are initiated by different molecular events. By analysing tumour genomic data from 33 cancer types, we show that the majority of tumours with high levels of W-CIN underwent whole genome doubling, whereas S-CIN levels are strongly associated with homologous recombination deficiency. Both CIN phenotypes are prognostic in several cancer types. Most drugs are less efficient in high-CIN cell lines, but we also report compounds and drugs which should be investigated as targets for W-CIN or S-CIN. By analysing associations between CIN and bio-molecular entities with pathway and gene expression levels, we complement gene signatures of CIN and report that the drug resistance gene CKS1B is strongly associated with S-CIN. Finally, we propose a potential copy number-dependent mechanism to activate the PI3K pathway in high-S-CIN tumours.

Loss of USP28 and SPINT2 expression promotes cancer cell survival after whole genome doubling.

BACKGROUND: Whole genome doubling is a frequent event during cancer evolution and shapes the cancer genome due to the occurrence of chromosomal instability. Yet, erroneously arising human tetraploid cells usually do not proliferate due to p53 activation that leads to CDKN1A expression, cell cycle arrest, senescence and/or apoptosis. METHODS: To uncover the barriers that block the proliferation of tetraploids, we performed a RNAi mediated genome-wide screen in a human colorectal cancer cell line (HCT116). RESULTS: We identified 140 genes whose depletion improved the survival of tetraploid cells and characterized in depth two of them: SPINT2 and USP28. We found that SPINT2 is a general regulator of CDKN1A transcription via histone acetylation. Using mass spectrometry and immunoprecipitation, we found that USP28 interacts with NuMA1 and affects centrosome clustering. Tetraploid cells accumulate DNA damage and loss of USP28 reduces checkpoint activation, thus facilitating their proliferation. CONCLUSIONS: Our results indicate three aspects that contribute to the survival of tetraploid cells: (i) increased mitogenic signaling and reduced expression of cell cycle inhibitors, (ii) the ability to establish functional bipolar spindles and (iii) reduced DNA damage signaling.

MFmap: A semi-supervised generative model matching cell lines to tumours and cancer subtypes.

Translating in vitro results from experiments with cancer cell lines to clinical applications requires the selection of appropriate cell line models. Here we present MFmap (model fidelity map), a machine learning model to simultaneously predict the cancer subtype of a cell line and its similarity to an individual tumour sample. The MFmap is a semi-supervised generative model, which compresses high dimensional gene expression, copy number variation and mutation data into cancer subtype informed low dimensional latent representations. The accuracy (test set F1 score >90%) of the MFmap subtype prediction is validated in ten different cancer datasets. We use breast cancer and glioblastoma cohorts as examples to show how subtype specific drug sensitivity can be translated to individual tumour samples. The low dimensional latent representations extracted by MFmap explain known and novel subtype specific features and enable the analysis of cell-state transformations between different subtypes. From a methodological perspective, we report that MFmap is a semi-supervised method which simultaneously achieves good generative and predictive performance and thus opens opportunities in other areas of computational biology.

Systems approaches identify the consequences of monosomy in somatic human cells.

Chromosome loss that results in monosomy is detrimental to viability, yet it is frequently observed in cancers. How cancers survive with monosomy is unknown. Using p53-deficient monosomic cell lines, we find that chromosome loss impairs proliferation and genomic stability. Transcriptome and proteome analysis demonstrates reduced expression of genes encoded on the monosomes, which is partially compensated in some cases. Monosomy also induces global changes in gene expression. Pathway enrichment analysis reveals that genes involved in ribosome biogenesis and translation are downregulated in all monosomic cells analyzed. Consistently, monosomies display defects in protein synthesis and ribosome assembly. We further show that monosomies are incompatible with p53 expression, likely due to defects in ribosome biogenesis. Accordingly, impaired ribosome biogenesis and p53 inactivation are associated with monosomy in cancer. Our systematic study of monosomy in human cells explains why monosomy is so detrimental and reveals the importance of p53 for monosomy occurrence in cancer.

The p53/p73 - p21CIP1 tumor suppressor axis guards against chromosomal instability by restraining CDK1 in human cancer cells.

Whole chromosome instability (W-CIN) is a hallmark of human cancer and contributes to the evolvement of aneuploidy. W-CIN can be induced by abnormally increased microtubule plus end assembly rates during mitosis leading to the generation of lagging chromosomes during anaphase as a major form of mitotic errors in human cancer cells. Here, we show that loss of the tumor suppressor genes TP53 and TP73 can trigger increased mitotic microtubule assembly rates, lagging chromosomes, and W-CIN. CDKN1A, encoding for the CDK inhibitor p21CIP1, represents a critical target gene of p53/p73. Loss of p21CIP1 unleashes CDK1 activity which causes W-CIN in otherwise chromosomally stable cancer cells. Consequently, induction of CDK1 is sufficient to induce abnormal microtubule assembly rates and W-CIN. Vice versa, partial inhibition of CDK1 activity in chromosomally unstable cancer cells corrects abnormal microtubule behavior and suppresses W-CIN. Thus, our study shows that the p53/p73 - p21CIP1 tumor suppressor axis, whose loss is associated with W-CIN in human cancer, safeguards against chromosome missegregation and aneuploidy by preventing abnormally increased CDK1 activity.

Germline ancestry influences the evolutionary disease course in lung adenocarcinomas.

Precision medicine relies on targeting specific somatic alterations present in a patient's tumor. However, the extent to which germline ancestry may influence the somatic burden of disease has received little attention. We estimated the genetic ancestry of non-small-cell lung cancer (NSCLC) patients and performed an in-depth analysis of the influence of genetic ancestry on the evolutionary disease course. Compared with European Americans (EA), African Americans (AA) with lung adenocarcinoma (LUAD) were found to be significantly younger and smoke significantly less. However, LUADs from AAs exhibited a significantly higher somatic mutation burden, with a more pronounced tobacco carcinogen footprint and increased frequencies of alterations affecting cancer genes. Conversely, no significant differences were observed between lung squamous cell carcinomas (LUSC) from EAs and AAs. Our results suggest germline ancestry influences the somatic evolution of LUAD but not LUSC.

Modelling and mathematical analysis of the M$_{2}$ receptor-dependent joint signalling and secondary messenger network in CHO cells.

The muscarinic M$_{2}$ receptor is a prominent member of the GPCR family and strongly involved in heart diseases. Recently published experimental work explored the cellular response to iperoxo-induced M$_{2}$ receptor stimulation in Chinese hamster ovary (CHO) cells. To better understand these responses, we modelled and analysed the muscarinic M$_{2}$ receptor-dependent signalling pathway combined with relevant secondary messenger molecules using mass action. In our literature-based joint signalling and secondary messenger model, all binding and phosphorylation events are explicitly taken into account in order to enable subsequent stoichiometric matrix analysis. We propose constraint flux sampling (CFS) as a method to characterize the expected shift of the steady state reaction flux distribution due to the known amount of cAMP production and PDE4 activation. CFS correctly predicts an experimentally observable influence on the cytoskeleton structure (marked by actin and tubulin) and in consequence a change of the optical density of cells. In a second step, we use CFS to simulate the effect of knock-out experiments within our biological system, and thus to rank the influence of individual molecules on the observed change of the optical cell density. In particular, we confirm the relevance of the protein RGS14, which is supported by current literature. A combination of CFS with Elementary Flux Mode analysis enabled us to determine the possible underlying mechanism. Our analysis suggests that mathematical tools developed for metabolic network analysis can also be applied to mixed secondary messenger and signalling models. This could be very helpful to perform model checking with little effort and to generate hypotheses for further research if parameters are not known.

Correction: The Subclonal Architecture of Metastatic Breast Cancer: Results from a Prospective Community-Based Rapid Autopsy Program "CASCADE".

[This corrects the article DOI: 10.1371/journal.pmed.1002204.].

The Subclonal Architecture of Metastatic Breast Cancer: Results from a Prospective Community-Based Rapid Autopsy Program "CASCADE".

BACKGROUND: Understanding the cancer genome is seen as a key step in improving outcomes for cancer patients. Genomic assays are emerging as a possible avenue to personalised medicine in breast cancer. However, evolution of the cancer genome during the natural history of breast cancer is largely unknown, as is the profile of disease at death. We sought to study in detail these aspects of advanced breast cancers that have resulted in lethal disease. METHODS AND FINDINGS: Three patients with oestrogen-receptor (ER)-positive, human epidermal growth factor receptor 2 (HER2)-negative breast cancer and one patient with triple negative breast cancer underwent rapid autopsy as part of an institutional prospective community-based rapid autopsy program (CASCADE). Cases represented a range of management problems in breast cancer, including late relapse after early stage disease, de novo metastatic disease, discordant disease response, and disease refractory to treatment. Between 5 and 12 metastatic sites were collected at autopsy together with available primary tumours and longitudinal metastatic biopsies taken during life. Samples underwent paired tumour-normal whole exome sequencing and single nucleotide polymorphism (SNP) arrays. Subclonal architectures were inferred by jointly analysing all samples from each patient. Mutations were validated using high depth amplicon sequencing. Between cases, there were significant differences in mutational burden, driver mutations, mutational processes, and copy number variation. Within each case, we found dramatic heterogeneity in subclonal structure from primary to metastatic disease and between metastatic sites, such that no single lesion captured the breadth of disease. Metastatic cross-seeding was found in each case, and treatment drove subclonal diversification. Subclones displayed parallel evolution of treatment resistance in some cases and apparent augmentation of key oncogenic drivers as an alternative resistance mechanism. We also observed the role of mutational processes in subclonal evolution. Limitations of this study include the potential for bias introduced by joint analysis of formalin-fixed archival specimens with fresh specimens and the difficulties in resolving subclones with whole exome sequencing. Other alterations that could define subclones such as structural variants or epigenetic modifications were not assessed. CONCLUSIONS: This study highlights various mechanisms that shape the genome of metastatic breast cancer and the value of studying advanced disease in detail. Treatment drives significant genomic heterogeneity in breast cancers which has implications for disease monitoring and treatment selection in the personalised medicine paradigm.

DNA replication stress mediates APOBEC3 family mutagenesis in breast cancer.

BACKGROUND: The APOBEC3 family of cytidine deaminases mutate the cancer genome in a range of cancer types. Although many studies have documented the downstream effects of APOBEC3 activity through next-generation sequencing, less is known about their upstream regulation. In this study, we sought to identify a molecular basis for APOBEC3 expression and activation. RESULTS: HER2 amplification and PTEN loss promote DNA replication stress and APOBEC3B activity in vitro and correlate with APOBEC3 mutagenesis in vivo. HER2-enriched breast carcinomas display evidence of elevated levels of replication stress-associated DNA damage in vivo. Chemical and cytotoxic induction of replication stress, through aphidicolin, gemcitabine, camptothecin or hydroxyurea exposure, activates transcription of APOBEC3B via an ATR/Chk1-dependent pathway in vitro. APOBEC3B activation can be attenuated through repression of oncogenic signalling, small molecule inhibition of receptor tyrosine kinase signalling and alleviation of replication stress through nucleoside supplementation. CONCLUSION: These data link oncogene, loss of tumour suppressor gene and drug-induced replication stress with APOBEC3B activity, providing new insights into how cytidine deaminase-induced mutagenesis might be activated in tumourigenesis and limited therapeutically.

Expression of regulators of mitotic fidelity are associated with intercellular heterogeneity and chromosomal instability in primary breast cancer.

Regulators of transition through mitosis such as SURVIVIN and Aurora kinase A (AURKA) have been previously implicated in the initiation of chromosomal instability (CIN), a driver of intratumour heterogeneity. We investigate the relationship between protein expression of these genes and directly quantified CIN, and their prognostic utility in breast cancer. The expression of SURVIVIN and AURKA was determined by immunohistochemistry in a cohort of 426 patients with primary breast cancer. The association between protein expression and histopathological characteristics, clinical outcome and CIN status, as determined by centromeric FISH and defined by modal centromere deviation, was analysed. Significantly poorer clinical outcome was observed in patients with high AURKA expression levels. Expression of SURVIVIN was elevated in ER-negative relative to ER-positive breast cancer. Both AURKA and SURVIVIN increased expression were significantly associated with breast cancer grade. There was a significant association between increased CIN and both increased AURKA and SURVIVIN expression. AURKA gene amplification was also associated with increased CIN. To our knowledge this is the largest study assessing CIN status in parallel with the expression of the mitotic regulators AURKA and SURVIVIN. These data suggest that elevated expression of AURKA and SURVIVIN, together with AURKA gene amplification, are associated with increased CIN in breast cancer, and may be used as a proxy for CIN in breast cancer samples in the absence of more advanced molecular measurements.

Chromosomal instability selects gene copy-number variants encoding core regulators of proliferation in ER+ breast cancer.

Chromosomal instability (CIN) is associated with poor outcome in epithelial malignancies, including breast carcinomas. Evidence suggests that prognostic signatures in estrogen receptor-positive (ER(+)) breast cancer define tumors with CIN and high proliferative potential. Intriguingly, CIN induction in lower eukaryotic cells and human cells is context dependent, typically resulting in a proliferation disadvantage but conferring a fitness benefit under strong selection pressures. We hypothesized that CIN permits accelerated genomic evolution through the generation of diverse DNA copy-number events that may be selected during disease development. In support of this hypothesis, we found evidence for selection of gene amplification of core regulators of proliferation in CIN-associated cancer genomes. Stable DNA copy-number amplifications of the core regulators TPX2 and UBE2C were associated with expression of a gene module involved in proliferation. The module genes were enriched within prognostic signature gene sets for ER(+) breast cancer, providing a logical connection between CIN and prognostic signature expression. Our results provide a framework to decipher the impact of intratumor heterogeneity on key cancer phenotypes, and they suggest that CIN provides a permissive landscape for selection of copy-number alterations that drive cancer proliferation.

Role of Saccharomyces cerevisiae Trk1 in stabilization of intracellular potassium content upon changes in external potassium levels.

Saccharomyces cerevisiae cells are able to grow at very different potassium concentrations adapting its intracellular cation levels to changes in the external milieu. Potassium homeostasis in wild type cells resuspended in media with low potassium is an example of non-perfect adaptation since the same intracellular concentration is not approached irrespective of the extracellular levels of the cation. By using yeasts lacking the Trk1,2 system or expressing different versions of the mutated main plasma membrane potassium transporter (Trk1), we show that Trk1 is not essential for adaptation to potassium changes but the dynamics of potassium loss is very different in the wild type and in trk1,2 mutant or in yeasts expressing Trk1 versions with highly impaired transport characteristics. We also show that the pattern here described can be also fulfilled by heterologous expression of NcHAK1, a potassium transporter not belonging to the TRK family. Hyperpolarization and cationic drugs sensitivity in mutants with defective transport capacity provide additional support to the hypothesis of connections between the activity of the Trk system and the plasma membrane H(+) ATPase (Pma1) in the adaptive process.

Replication stress links structural and numerical cancer chromosomal instability.

Cancer chromosomal instability (CIN) results in an increased rate of change of chromosome number and structure and generates intratumour heterogeneity. CIN is observed in most solid tumours and is associated with both poor prognosis and drug resistance. Understanding a mechanistic basis for CIN is therefore paramount. Here we find evidence for impaired replication fork progression and increased DNA replication stress in CIN(+) colorectal cancer (CRC) cells relative to CIN(-) CRC cells, with structural chromosome abnormalities precipitating chromosome missegregation in mitosis. We identify three new CIN-suppressor genes (PIGN (also known as MCD4), MEX3C (RKHD2) and ZNF516 (KIAA0222)) encoded on chromosome 18q that are subject to frequent copy number loss in CIN(+) CRC. Chromosome 18q loss was temporally associated with aneuploidy onset at the adenoma-carcinoma transition. CIN-suppressor gene silencing leads to DNA replication stress, structural chromosome abnormalities and chromosome missegregation. Supplementing cells with nucleosides, to alleviate replication-associated damage, reduces the frequency of chromosome segregation errors after CIN-suppressor gene silencing, and attenuates segregation errors and DNA damage in CIN(+) cells. These data implicate a central role for replication stress in the generation of structural and numerical CIN, which may inform new therapeutic approaches to limit intratumour heterogeneity.

LRIG1 regulates cadherin-dependent contact inhibition directing epithelial homeostasis and pre-invasive squamous cell carcinoma development.

Epidermal growth factor receptor (EGFR) pathway activation is a frequent event in human carcinomas. Mutations in EGFR itself are, however, rare, and the mechanisms regulating EGFR activation remain elusive. Leucine-rich immunoglobulin repeats-1 (LRIG1), an inhibitor of EGFR activity, is one of four genes identified that predict patient survival across solid tumour types including breast, lung, melanoma, glioma, and bladder. We show that deletion of Lrig1 is sufficient to promote murine airway hyperplasia through loss of contact inhibition and that re-expression of LRIG1 in human lung cancer cells inhibits tumourigenesis. LRIG1 regulation of contact inhibition occurs via ternary complex formation with EGFR and E-cadherin with downstream modulation of EGFR activity. We find that LRIG1 LOH is frequent across cancers and its loss is an early event in the development of human squamous carcinomas. Our findings imply that the early stages of squamous carcinoma development are driven by a change in amplitude of EGFR signalling governed by the loss of contact inhibition.

Potassium starvation in yeast: mechanisms of homeostasis revealed by mathematical modeling.

The intrinsic ability of cells to adapt to a wide range of environmental conditions is a fundamental process required for survival. Potassium is the most abundant cation in living cells and is required for essential cellular processes, including the regulation of cell volume, pH and protein synthesis. Yeast cells can grow from low micromolar to molar potassium concentrations and utilize sophisticated control mechanisms to keep the internal potassium concentration in a viable range. We developed a mathematical model for Saccharomyces cerevisiae to explore the complex interplay between biophysical forces and molecular regulation facilitating potassium homeostasis. By using a novel inference method ("the reverse tracking algorithm") we predicted and then verified experimentally that the main regulators under conditions of potassium starvation are proton fluxes responding to changes of potassium concentrations. In contrast to the prevailing view, we show that regulation of the main potassium transport systems (Trk1,2 and Nha1) in the plasma membrane is not sufficient to achieve homeostasis.

CERT depletion predicts chemotherapy benefit and mediates cytotoxic and polyploid-specific cancer cell death through autophagy induction.

Chromosomal instability (CIN) has been implicated in multidrug resistance and the silencing of the ceramide transporter, CERT, promotes sensitization to diverse cytotoxics. An improved understanding of mechanisms governing multidrug sensitization might provide insight into pathways contributing to the death of CIN cancer cells. Using an integrative functional genomics approach, we find that CERT-specific multidrug sensitization is associated with enhanced autophagosome-lysosome flux, resulting from the expression of LAMP2 following CERT silencing in colorectal and HER2(+) breast cancer cell lines. Live cell microscopy analysis revealed that CERT depletion induces LAMP2-dependent death of polyploid cells following exit from mitosis in the presence of paclitaxel. We find that CERT is relatively over-expressed in HER2(+) breast cancer and CERT protein expression acts as an independent prognostic variable and predictor of outcome in adjuvant chemotherapy-treated patients with primary breast cancer. These data suggest that the induction of LAMP2-dependent autophagic flux through CERT targeting may provide a rational approach to enhance multidrug sensitization and potentiate the death of polyploid cells following paclitaxel exposure to limit the acquisition of CIN and intra-tumour heterogeneity.