scAbsolute: measuring single-cell ploidy and replication status.

Cancer cells often exhibit DNA copy number aberrations and can vary widely in their ploidy. Correct estimation of the ploidy of single-cell genomes is paramount for downstream analysis. Based only on single-cell DNA sequencing information, scAbsolute achieves accurate and unbiased measurement of single-cell ploidy and replication status, including whole-genome duplications. We demonstrate scAbsolute's capabilities using experimental cell multiplets, a FUCCI cell cycle expression system, and a benchmark against state-of-the-art methods. scAbsolute provides a robust foundation for single-cell DNA sequencing analysis across different technologies and has the potential to enable improvements in a number of downstream analyses.

The copy number and mutational landscape of recurrent ovarian high-grade serous carcinoma.

The drivers of recurrence and resistance in ovarian high grade serous carcinoma remain unclear. We investigate the acquisition of resistance by collecting tumour biopsies from a cohort of 276 women with relapsed ovarian high grade serous carcinoma in the BriTROC-1 study. Panel sequencing shows close concordance between diagnosis and relapse, with only four discordant cases. There is also very strong concordance in copy number between diagnosis and relapse, with no significant difference in purity, ploidy or focal somatic copy number alterations, even when stratified by platinum sensitivity or prior chemotherapy lines. Copy number signatures are strongly correlated with immune cell infiltration, whilst diagnosis samples from patients with primary platinum resistance have increased rates of CCNE1 and KRAS amplification and copy number signature 1 exposure. Our data show that the ovarian high grade serous carcinoma genome is remarkably stable between diagnosis and relapse and acquired chemotherapy resistance does not select for common copy number drivers.

High-grade serous ovarian carcinoma organoids as models of chromosomal instability.

High-grade serous ovarian carcinoma (HGSOC) is the most genomically complex cancer, characterized by ubiquitous TP53 mutation, profound chromosomal instability, and heterogeneity. The mutational processes driving chromosomal instability in HGSOC can be distinguished by specific copy number signatures. To develop clinically relevant models of these mutational processes we derived 15 continuous HGSOC patient-derived organoids (PDOs) and characterized them using bulk transcriptomic, bulk genomic, single-cell genomic, and drug sensitivity assays. We show that HGSOC PDOs comprise communities of different clonal populations and represent models of different causes of chromosomal instability including homologous recombination deficiency, chromothripsis, tandem-duplicator phenotype, and whole genome duplication. We also show that these PDOs can be used as exploratory tools to study transcriptional effects of copy number alterations as well as compound-sensitivity tests. In summary, HGSOC PDO cultures provide validated genomic models for studies of specific mutational processes and precision therapeutics.

A pan-cancer compendium of chromosomal instability.

Chromosomal instability (CIN) results in the accumulation of large-scale losses, gains and rearrangements of DNA1. The broad genomic complexity caused by CIN is a hallmark of cancer2; however, there is no systematic framework to measure different types of CIN and their effect on clinical phenotypes pan-cancer. Here we evaluate the extent, diversity and origin of CIN across 7,880 tumours representing 33 cancer types. We present a compendium of 17 copy number signatures that characterize specific types of CIN, with putative aetiologies supported by multiple independent data sources. The signatures predict drug response and identify new drug targets. Our framework refines the understanding of impaired homologous recombination, which is one of the most therapeutically targetable types of CIN. Our results illuminate a fundamental structure underlying genomic complexity in human cancers and provide a resource to guide future CIN research.

CNpare: matching DNA copy number profiles.

SUMMARY: Selecting the optimal cancer cell line for an experiment can be challenging given the diversity of lines available. Here, we present CNpare, which identifies similar cell line models based on genome-wide DNA copy number. AVAILABILITY AND IMPLEMENTATION: CNpare is available as an R package at https://github.com/macintyrelab/CNpare. All analysis performed in the manuscript can be reproduced via the code found at https://github.com/macintyrelab/CNpare_analyses. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Prognostic markers of inflammation in endometrioid and clear cell ovarian cancer.

OBJECTIVES: Cancer-related systemic inflammation has been associated with prognosis in multiple cancer types. Conversely, local inflammation, which is characterized by dense intratumoral immune infiltrates, is a favorable predictor of survival outcome. However, these survival associations are not well established in ovarian cancer, particularly in the less frequent endometrioid and clear cell endometriosis associated histotypes. METHODS: This retrospective study included 119 patients (63 endometrioid and 56 clear cell ovarian carcinomas). We performed a comprehensive survival association analysis of both systemic (neutrophil-to-lymphocyte ratio or presence of endometriosis) and local inflammation markers (CD3+ and CD8+ tumor infiltrating lymphocytes) using multivariate Cox proportional hazards models that account for confounding factors. RESULTS: Medium to high levels of intraepithelial CD8+ tumor infiltrating lymphocytes are associated with longer survival in endometrioid ovarian cancer (p=0.04). In addition, we found that intraepithelial CD8+ tumor infiltrating lymphocytes are prognostic in clear cell ovarian cancer (p=0.02), and that intraepithelial CD3+ tumor infiltrating lymphocytes are also associated with improved outcome (p=0.02). Furthermore, intratumoral CD3+ and CD8+ tumor infiltrating lymphocytes showed improved prognosis in the endometrioid subtype (p<0.1). No prognostic value was observed for systemic immune markers. CONCLUSIONS: In this study, patients with endometrioid and clear cell ovarian cancer with moderate to high CD8+ and CD3+ intraepithelial tumor infiltrating lymphocytes had longer overall survival. Higher expression of intratumoral CD3+ and CD8+ tumor infiltrating lymphocytes also showed an improved outcome in endometrioid ovarian cancer. In contrast, systemic inflammation, evaluated by neutrophil-to-lymphocyte ratio or presence of endometriosis, did not have a prognostic impact in these histologic subtypes.

Loss of SNAI2 in Prostate Cancer Correlates With Clinical Response to Androgen Deprivation Therapy.

PURPOSE: Androgen receptor (AR) signaling is important in prostate cancer progression, and therapies that target this pathway have been the mainstay of treatment for advanced disease for over 70 years. Tumors eventually progress despite castration through a number of well-characterized mechanisms; however, little is known about what determines the magnitude of response to short-term pathway inhibition. METHODS: We evaluated a novel combination of AR-targeting therapies (degarelix, abiraterone, and bicalutamide) and noted that the objective patient response to therapy was highly variable. To investigate what was driving treatment resistance in poorly responding patients, as a secondary outcome we comprehensively characterized pre- and post-treatment samples using both whole-genome and RNA sequencing. RESULTS: We find that resistance following short-term treatment differs molecularly from typical progressive castration-resistant disease, associated with transcriptional reprogramming, to a transitional epithelial-to-mesenchymal transition (EMT) phenotype rather than an upregulation of AR signaling. Unexpectedly, tolerance to therapy appears to be the default state, with treatment response correlating with the prevalence of tumor cells deficient for SNAI2, a key regulator of EMT reprogramming. CONCLUSION: We show that EMT characterizes acutely resistant prostate tumors and that deletion of SNAI2, a key transcriptional regulator of EMT, correlates with clinical response.

MSH2-deficient prostate tumours have a distinct immune response and clinical outcome compared to MSH2-deficient colorectal or endometrial cancer.

BACKGROUND: Recent publications have shown patients with defects in the DNA mismatch repair (MMR) pathway driven by either MSH2 or MSH6 loss experience a significant increase in the incidence of prostate cancer. Moreover, this increased incidence of prostate cancer is accompanied by rapid disease progression and poor clinical outcomes. METHODS AND RESULTS: We show that androgen-receptor activation, a key driver of prostate carcinogenesis, can disrupt the MSH2 gene in prostate cancer. We screened tumours from two cohorts (recurrent/non-recurrent) of prostate cancer patients to confirm the loss of MSH2 protein expression and identified decreased MSH2 expression in recurrent cases. Stratifying the independent TCGA prostate cancer cohort for MSH2/6 expression revealed that patients with lower levels of MSH2/6 had significant worse outcomes, in contrast, endometrial and colorectal cancer patients with lower MSH2/6 levels. MMRd endometrial and colorectal tumours showed the expected increase in mutational burden, microsatellite instability and enhanced immune cell mobilisation but this was not evident in prostate tumours. CONCLUSIONS: We have shown that loss or reduced levels of MSH2/MSH6 protein in prostate cancer is associated with poor outcome. However, our data indicate that this is not associated with a statistically significant increase in mutational burden, microsatellite instability or immune cell mobilisation in a cohort of primary prostate cancers.

Characterizing genetic intra-tumor heterogeneity across 2,658 human cancer genomes.

Intra-tumor heterogeneity (ITH) is a mechanism of therapeutic resistance and therefore an important clinical challenge. However, the extent, origin, and drivers of ITH across cancer types are poorly understood. To address this, we extensively characterize ITH across whole-genome sequences of 2,658 cancer samples spanning 38 cancer types. Nearly all informative samples (95.1%) contain evidence of distinct subclonal expansions with frequent branching relationships between subclones. We observe positive selection of subclonal driver mutations across most cancer types and identify cancer type-specific subclonal patterns of driver gene mutations, fusions, structural variants, and copy number alterations as well as dynamic changes in mutational processes between subclonal expansions. Our results underline the importance of ITH and its drivers in tumor evolution and provide a pan-cancer resource of comprehensively annotated subclonal events from whole-genome sequencing data.

FrenchFISH: Poisson Models for Quantifying DNA Copy Number From Fluorescence In Situ Hybridization of Tissue Sections.

PURPOSE: Chromosomal aberration and DNA copy number change are robust hallmarks of cancer. The gold standard for detecting copy number changes in tumor cells is fluorescence in situ hybridization (FISH) using locus-specific probes that are imaged as fluorescent spots. However, spot counting often does not perform well on solid tumor tissue sections due to partially represented or overlapping nuclei. MATERIALS AND METHODS: To overcome these challenges, we have developed a computational approach called FrenchFISH, which comprises a nuclear volume correction method coupled with two types of Poisson models: either a Poisson model for improved manual spot counting without the need for control probes or a homogeneous Poisson point process model for automated spot counting. RESULTS: We benchmarked the performance of FrenchFISH against previous approaches using a controlled simulation scenario and tested it experimentally in 12 ovarian carcinoma FFPE-tissue sections for copy number alterations at three loci (c-Myc, hTERC, and SE7). FrenchFISH outperformed standard spot counting with 74% of the automated counts having < 1 copy number difference from the manual counts and 17% having < 2 copy number differences, while taking less than one third of the time of manual counting. CONCLUSION: FrenchFISH is a general approach that can be used to enhance clinical diagnosis on sections of any tissue by both speeding up and improving the accuracy of spot count estimates.

Unraveling tumor-immune heterogeneity in advanced ovarian cancer uncovers immunogenic effect of chemotherapy.

In metastatic cancer, the degree of heterogeneity of the tumor microenvironment (TME) and its molecular underpinnings remain largely unstudied. To characterize the tumor-immune interface at baseline and during neoadjuvant chemotherapy (NACT) in high-grade serous ovarian cancer (HGSOC), we performed immunogenomic analysis of treatment-naive and paired samples from before and after treatment with chemotherapy. In treatment-naive HGSOC, we found that immune-cell-excluded and inflammatory microenvironments coexist within the same individuals and within the same tumor sites, indicating ubiquitous variability in immune cell infiltration. Analysis of TME cell composition, DNA copy number, mutations and gene expression showed that immune cell exclusion was associated with amplification of Myc target genes and increased expression of canonical Wnt signaling in treatment-naive HGSOC. Following NACT, increased natural killer (NK) cell infiltration and oligoclonal expansion of T cells were detected. We demonstrate that the tumor-immune microenvironment of advanced HGSOC is intrinsically heterogeneous and that chemotherapy induces local immune activation, suggesting that chemotherapy can potentiate the immunogenicity of immune-excluded HGSOC tumors.

Genomic landscape of platinum resistant and sensitive testicular cancers.

While most testicular germ cell tumours (TGCTs) exhibit exquisite sensitivity to platinum chemotherapy, ~10% are platinum resistant. To gain insight into the underlying mechanisms, we undertake whole exome sequencing and copy number analysis in 40 tumours from 26 cases with platinum-resistant TGCT, and combine this with published genomic data on an additional 624 TGCTs. We integrate analyses for driver mutations, mutational burden, global, arm-level and focal copy number (CN) events, and SNV and CN signatures. Albeit preliminary and observational in nature, these analyses provide support for a possible mechanistic link between early driver mutations in RAS and KIT and the widespread copy number events by which TGCT is characterised.

Inferring structural variant cancer cell fraction.

We present SVclone, a computational method for inferring the cancer cell fraction of structural variant (SV) breakpoints from whole-genome sequencing data. SVclone accurately determines the variant allele frequencies of both SV breakends, then simultaneously estimates the cancer cell fraction and SV copy number. We assess performance using in silico mixtures of real samples, at known proportions, created from two clonal metastases from the same patient. We find that SVclone's performance is comparable to single-nucleotide variant-based methods, despite having an order of magnitude fewer data points. As part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) consortium, which aggregated whole-genome sequencing data from 2658 cancers across 38 tumour types, we use SVclone to reveal a subset of liver, ovarian and pancreatic cancers with subclonally enriched copy-number neutral rearrangements that show decreased overall survival. SVclone enables improved characterisation of SV intra-tumour heterogeneity.

The evolutionary history of 2,658 cancers.

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

Copy number signatures and mutational processes in ovarian carcinoma.

The genomic complexity of profound copy number aberrations has prevented effective molecular stratification of ovarian cancers. Here, to decode this complexity, we derived copy number signatures from shallow whole-genome sequencing of 117 high-grade serous ovarian cancer (HGSOC) cases, which were validated on 527 independent cases. We show that HGSOC comprises a continuum of genomes shaped by multiple mutational processes that result in known patterns of genomic aberration. Copy number signature exposures at diagnosis predict both overall survival and the probability of platinum-resistant relapse. Measurement of signature exposures provides a rational framework to choose combination treatments that target multiple mutational processes.

Periprostatic fat tissue transcriptome reveals a signature diagnostic for high-risk prostate cancer.

Evidence suggests that altered adipose tissue homeostasis may be an important contributor to the development and/or progression of prostate cancer. In this study, we investigated the adipose transcriptional profiles of low- and high-risk disease to determine both prognostic potential and possible biological drivers of aggressive disease. RNA was extracted from periprostatic adipose tissue from patients categorised as having prostate cancer with either a low or high risk of progression based on tumour characteristics at prostatectomy and profiled by RNA sequencing. The expression of selected genes was then quantified by qRT-PCR in a cross-validation cohort. In the first phase, a total of 677 differentially transcribed genes were identified, from which a subset of 14 genes was shortlisted. In the second phase, a 3 gene (IGHA1, OLFM4, RERGL) signature was refined and evaluated using recursive feature selection and cross-validation, obtaining a promising discriminatory utility (area under curve 0.72) at predicting the presence of high-risk disease. Genes implicated in immune and/or inflammatory responses predominated. Periprostatic adipose tissue from patients with high-risk prostate cancer has a distinct transcriptional signature that may be useful for detecting its occult presence. Differential expression appears to be driven by a local immune/inflammatory reaction to more advanced tumours, than any specific adipose tissue-specific tumour-promoting mechanism. This signature is transferable into a clinically usable PCR-based assay, which in a cross-validation cohort shows diagnostic potential.