Integrated structural variation and point mutation signatures in cancer genomes using correlated topic models

Predictive biomarkers for the efficacy of PARP inhibitors in ovarian cancer: an updated systematic review

BACKGROUND

PARP inhibitors are effective in treating ovarian cancer, especially for BRCA1/2 pathogenic variant carriers and those with HRD (homologous recombination deficiency). Concerns over toxicity and costs have led to the search for predictive biomarkers. We present an updated systematic review, expanding on a previous ESMO review on PARP inhibitor biomarkers.

METHODS

Following ESMO's 2020 review protocol, we extended our search to March 31, 2023, including PubMed and clinical trial data. We also reviewed the reference lists of review articles. We conducted a meta-analysis using a random-effects model to evaluate hazard ratios and assess the predictive potential of biomarkers and the effectiveness of PARP inhibitors in survival.

RESULTS

We found 375 articles, 103 of which were included after screening (62 primary research, 41 reviews). HRD remained the primary biomarker (95%), particularly BRCA1/2 variants (77%). In the non-HRD category, six articles (10%) introduced innovative biomarkers, including ADP-ribosylation, HOXA9 promoter methylation, patient-derived organoids, KELIM, and SLFN11.

DISCUSSION

Prospective assessment of real-time homologous recombination repair via nuclear RAD51 levels shows promise but needs validation. Emerging biomarkers like ADP-ribosylation, HOXA9 promoter methylation, patient-derived organoids, KELIM, and SLFN11 offer potential but require large-scale validation.

Realizing the promise of machine learning in precision oncology: expert perspectives on opportunities and challenges

BACKGROUND

The ability of machine learning (ML) to process and learn from large quantities of heterogeneous patient data is gaining attention in the precision oncology community. Some remarkable developments have taken place in the domain of image classification tasks in areas such as digital pathology and diagnostic radiology. The application of ML approaches to the analysis of DNA data, including tumor-derived genomic profiles, microRNAs, and cancer epigenetic signatures, while relatively more recent, has demonstrated some utility in identifying driver variants and molecular signatures with possible prognostic and therapeutic applications.

METHODS

We conducted semi-structured interviews with academic and clinical experts to capture the status quo, challenges, opportunities, ethical implications, and future directions.

RESULTS

Our participants agreed that machine learning in precision oncology is in infant stages, with clinical integration still rare. Overall, participants equated ongoing developments with better clinical workflows and improved treatment decisions for more cancer patients. They underscored the ability of machine learning to tackle the dynamic nature of cancer, break down the complexity of molecular data, and support decision-making. Our participants emphasized obstacles related to molecular data access, clinical utility, and guidelines. The availability of reliable and well-curated data to train and validate machine learning algorithms and integrate multiple data sources were described as constraints yet necessary for future clinical implementation. Frequently mentioned ethical challenges included privacy risks, equity, explainability, trust, and incidental findings, with privacy being the most polarizing. While participants recognized the issue of hype surrounding machine learning in precision oncology, they agreed that, in an assistive role, it represents the future of precision oncology.

CONCLUSIONS

Given the unique nature of medical AI, our findings highlight the field's potential and remaining challenges. ML will continue to advance cancer research and provide opportunities for patient-centric, personalized, and efficient precision oncology. Yet, the field must move beyond hype and toward concrete efforts to overcome key obstacles, such as ensuring access to molecular data, establishing clinical utility, developing guidelines and regulations, and meaningfully addressing ethical challenges.

Identifying somatic fingerprints of cancers defined by germline and environmental risk factors

Numerous studies over the past generation have identified germline variants that increase specific cancer risks. Simultaneously, a revolution in sequencing technology has permitted high-throughput annotations of somatic genomes characterizing individual tumors. However, examining the relationship between germline variants and somatic alteration patterns is hugely challenged by the large numbers of variants in a typical tumor, the rarity of most individual variants, and the heterogeneity of tumor somatic fingerprints. In this article, we propose statistical methodology that frames the investigation of germline-somatic relationships in an interpretable manner. The method uses meta-features embodying biological contexts of individual somatic alterations to implicitly group rare mutations. Our team has used this technique previously through a multilevel regression model to diagnose with high accuracy tumor site of origin. Herein, we further leverage topic models from computational linguistics to achieve interpretable lower-dimensional embeddings of the meta-features. We demonstrate how the method can identify distinctive somatic profiles linked to specific germline variants or environmental risk factors. We illustrate the method using The Cancer Genome Atlas whole-exome sequencing data to characterize somatic tumor fingerprints in breast cancer patients with germline BRCA1/2 mutations and in head and neck cancer patients exposed to human papillomavirus.

Unlocking the potential of immunotherapy in platinum-resistant ovarian cancer: rationale, challenges, and novel strategies

Ovarian cancer is a significant global health challenge, with cytoreductive surgery and platinum-based chemotherapy serving as established primary treatments. Unfortunately, most patients relapse and ultimately become platinum-resistant, at which point there are limited effective treatment options. Given the success of immunotherapy in inducing durable treatment responses in several other cancers, its potential in platinum-resistant ovarian cancer (PROC) is currently being investigated. However, in unselected advanced ovarian cancer populations, researchers have reported low response rates to immune checkpoint inhibition, and thus far, no validated biomarkers are predictive of response. Understanding the intricate interplay between platinum resistance, immune recognition, and the tumour microenvironment (TME) is crucial. In this review, we examine the research challenges encountered thus far, the biological rationale for immunotherapy, the underlying mechanisms of immune resistance, and new strategies to overcome resistance.

Inferring replication timing and proliferation dynamics from single-cell DNA sequencing data

Dysregulated DNA replication is a cause and a consequence of aneuploidy in cancer, yet the interplay between copy number alterations (CNAs), replication timing (RT) and cell cycle dynamics remain understudied in aneuploid tumors. We developed a probabilistic method, PERT, for simultaneous inference of cell-specific replication and copy number states from single-cell whole genome sequencing (scWGS) data. We used PERT to investigate clone-specific RT and proliferation dynamics in  >50,000 cells obtained from aneuploid and clonally heterogeneous cell lines, xenografts and primary cancers. We observed bidirectional relationships between RT and CNAs, with CNAs affecting X-inactivation producing the largest RT shifts. Additionally, we found that clone-specific S-phase enrichment positively correlated with ground-truth proliferation rates in genomically stable but not unstable cells. Together, these results demonstrate robust computational identification of S-phase cells from scWGS data, and highlight the importance of RT and cell cycle properties in studying the genomic evolution of aneuploid tumors.

Integrating PARP Inhibitors in mCRPC Therapy: Current Strategies and Emerging Trends

Metastatic castrate-resistant prostate cancer (mCRPC) is associated with poor prognosis. DNA damage response (DDR) genes are commonly altered in mCRPC rendering them as promising therapeutic targets. Poly (ADP ribose) polymerase inhibitors (PARPi) demonstrated antitumor activity in mCRPC patients with DDR gene mutations through synthetic lethality. Multiple clinical trials with PARPi monotherapy exhibited encouraging clinical outcomes in selected patients with mCRPC. More recently, three Phase III randomized clinical trials (RCTs) combining PARPi with androgen receptor signaling inhibitors (ARSIs) demonstrated improved antitumor activity compared to ARSI monotherapy in mCRPC patients as the first-line therapy. Clinical benefit was more pronounced in patients harboring DDR alterations, specifically BRCA1/2. Interestingly, antitumor activity was also observed irrespective of DDR gene mutations, highlighting BRCAness phenotype with androgen receptor blockade resulting in synergistic activity between ARSIs and PARPi. In this review, we discuss the clinical efficacy and safety data of the combination of PARPi plus ARSI in all Phase 3 randomized controlled trials (RCTs), emphasizing strategies for patient selection and highlighting emerging trends based on clinical trial data.

Tracking clonal evolution of drug resistance in ovarian cancer patients by exploiting structural variants in cfDNA

Drug resistance is the major cause of therapeutic failure in high-grade serous ovarian cancer (HGSOC). Yet, the mechanisms by which tumors evolve to drug resistant states remains largely unknown. To address this, we aimed to exploit clone-specific genomic structural variations by combining scaled single-cell whole genome sequencing with longitudinally collected cell-free DNA (cfDNA), enabling clonal tracking before, during and after treatment. We developed a cfDNA hybrid capture, deep sequencing approach based on leveraging clone-specific structural variants as endogenous barcodes, with orders of magnitude lower error rates than single nucleotide variants in ctDNA (circulating tumor DNA) detection, demonstrated on 19 patients at baseline. We then applied this to monitor and model clonal evolution over several years in ten HGSOC patients treated with systemic therapy from diagnosis through recurrence. We found drug resistance to be polyclonal in most cases, but frequently dominated by a single high-fitness and expanding clone, reducing clonal diversity in the relapsed disease state in most patients. Drug-resistant clones frequently displayed notable genomic features, including high-level amplifications of oncogenes such as CCNE1, RAB25, NOTCH3, and ERBB2. Using a population genetics Wright-Fisher model, we found evolutionary trajectories of these features were consistent with drug-induced positive selection. In select cases, these alterations impacted selection of secondary lines of therapy with positive patient outcomes. For cases with matched single-cell RNA sequencing data, pre-existing and genomically encoded phenotypic states such as upregulation of EMT and VEGF were linked to drug resistance. Together, our findings indicate that drug resistant states in HGSOC pre-exist at diagnosis and lead to dramatic clonal expansions that alter clonal composition at the time of relapse. We suggest that combining tumor single cell sequencing with cfDNA enables clonal tracking in patients and harbors potential for evolution-informed adaptive treatment decisions.

Ongoing genome doubling promotes evolvability and immune dysregulation in ovarian cancer

Whole-genome doubling (WGD) is a critical driver of tumor development and is linked to drug resistance and metastasis in solid malignancies. Here, we demonstrate that WGD is an ongoing mutational process in tumor evolution. Using single-cell whole-genome sequencing, we measured and modeled how WGD events are distributed across cellular populations within tumors and associated WGD dynamics with properties of genome diversification and phenotypic consequences of innate immunity. We studied WGD evolution in 65 high-grade serous ovarian cancer (HGSOC) tissue samples from 40 patients, yielding 29,481 tumor cell genomes. We found near-ubiquitous evidence of WGD as an ongoing mutational process promoting cell-cell diversity, high rates of chromosomal missegregation, and consequent micronucleation. Using a novel mutation-based WGD timing method, doubleTime , we delineated specific modes by which WGD can drive tumor evolution: (i) unitary evolutionary origin followed by significant diversification, (ii) independent WGD events on a pre-existing background of copy number diversity, and (iii) evolutionarily late clonal expansions of WGD populations. Additionally, through integrated single-cell RNA sequencing and high-resolution immunofluorescence microscopy, we found that inflammatory signaling and cGAS-STING pathway activation result from ongoing chromosomal instability and are restricted to tumors that remain predominantly diploid. This contrasted with predominantly WGD tumors, which exhibited significant quiescent and immunosuppressive phenotypic states. Together, these findings establish WGD as an evolutionarily 'active' mutational process that promotes evolvability and dysregulated immunity in late stage ovarian cancer.

A text analytics approach for mining public discussions in online cancer forum: Analysis of multi-intent lung cancer treatment dataset

BACKGROUND

Online cancer forums (OCF) are increasingly popular platforms for patients and caregivers to discuss, seek information on, and share opinions about diseases and treatments. This interaction generates a substantial amount of unstructured text data, necessitating deeper exploration. Using time series data, our study exploits topic modeling in the novel domain of online cancer forums (OCFs) to identify meaningful topics and changing dynamics of online discussion across different lung cancer treatment intent groups.

METHODS

For this purpose, a dataset comprising 27,998 forum posts about lung cancer was collected from three OCFs: lungcancer.net, lungevity.org, and reddit.com, spanning the years 2016 to 2018.

RESULTS

The analysis reflects the public discussion on multi-intent lung cancer treatment over time, taking into account seasonal variations. Discussions on cancer symptoms and prevention garnered the most attention, dominating both curative and palliative care discussions. There were distinct seasonal peaks: curative care topics surged from winter to late spring, while palliative care topics peaked from late summer to mid-autumn. Keyword analysis highlighted that lung cancer diagnosis and treatment were primary topics, whereas cancer prevention and treatment outcomes were predominant across multi-care contexts. For the study period, curative care discussions predominantly revolved around informational support and disease syndromes. In contrast, social support and cancer prevention prevailed in the palliative care context. Notably, topics such as cancer screening and cancer treatment exhibit pronounced seasonal variations in curative care, peaking in frequency during the summers (May to August) of the study period. Meanwhile, the topic of tumor control within palliative care showed significant seasonal influence during the winters and summers of 2017 and 2018.

CONCLUSION

Our text analysis approach using OCF data shows potential for computational methods in this novel domain to gain insights into trends in public cancer communication and seasonal variations for a better understanding of improving personalized care, decision support, treatment outcomes, and quality of life.

Topical hidden genome: discovering latent cancer mutational topics using a Bayesian multilevel context-learning approach

Inferring the cancer-type specificities of ultra-rare, genome-wide somatic mutations is an open problem. Traditional statistical methods cannot handle such data due to their ultra-high dimensionality and extreme data sparsity. To harness information in rare mutations, we have recently proposed a formal multilevel multilogistic "hidden genome" model. Through its hierarchical layers, the model condenses information in ultra-rare mutations through meta-features embodying mutation contexts to characterize cancer types. Consistent, scalable point estimation of the model can incorporate 10s of millions of variants across thousands of tumors and permit impressive prediction and attribution. However, principled statistical inference is infeasible due to the volume, correlation, and noninterpretability of mutation contexts. In this paper, we propose a novel framework that leverages topic models from computational linguistics to effectuate dimension reduction of mutation contexts producing interpretable, decorrelated meta-feature topics. We propose an efficient MCMC algorithm for implementation that permits rigorous full Bayesian inference at a scale that is orders of magnitude beyond the capability of existing out-of-the-box inferential high-dimensional multi-class regression methods and software. Applying our model to the Pan Cancer Analysis of Whole Genomes dataset reveals interesting biological insights including somatic mutational topics associated with UV exposure in skin cancer, aging in colorectal cancer, and strong influence of epigenome organization in liver cancer. Under cross-validation, our model demonstrates highly competitive predictive performance against blackbox methods of random forest and deep learning.

Proteogenomic analysis of enriched HGSOC tumor epithelium identifies prognostic signatures and therapeutic vulnerabilities

We performed a deep proteogenomic analysis of bulk tumor and laser microdissection enriched tumor cell populations from high-grade serous ovarian cancer (HGSOC) tissue specimens spanning a broad spectrum of purity. We identified patients with longer progression-free survival had increased immune-related signatures and validated proteins correlating with tumor-infiltrating lymphocytes in 65 tumors from an independent cohort of HGSOC patients, as well as with overall survival in an additional 126 HGSOC patient cohort. We identified that homologous recombination deficient (HRD) tumors are enriched in pathways associated with metabolism and oxidative phosphorylation that we validated in independent patient cohorts. We further identified that polycomb complex protein BMI-1 is elevated in HR proficient (HRP) tumors, that elevated BMI-1 correlates with poor overall survival in HRP but not HRD HGSOC patients, and that HRP HGSOC cells are uniquely sensitive to BMI-1 inhibition.

Accurate and sensitive mutational signature analysis with MuSiCal

Mutational signature analysis is a recent computational approach for interpreting somatic mutations in the genome. Its application to cancer data has enhanced our understanding of mutational forces driving tumorigenesis and demonstrated its potential to inform prognosis and treatment decisions. However, methodological challenges remain for discovering new signatures and assigning proper weights to existing signatures, thereby hindering broader clinical applications. Here we present Mutational Signature Calculator (MuSiCal), a rigorous analytical framework with algorithms that solve major problems in the standard workflow. Our simulation studies demonstrate that MuSiCal outperforms state-of-the-art algorithms for both signature discovery and assignment. By reanalyzing more than 2,700 cancer genomes, we provide an improved catalog of signatures and their assignments, discover nine indel signatures absent in the current catalog, resolve long-standing issues with the ambiguous 'flat' signatures and give insights into signatures with unknown etiologies. We expect MuSiCal and the improved catalog to be a step towards establishing best practices for mutational signature analysis.

Unsupervised learning for medical data: A review of probabilistic factorization methods

We review popular unsupervised learning methods for the analysis of high-dimensional data encountered in, for example, genomics, medical imaging, cohort studies, and biobanks. We show that four commonly used methods, principal component analysis, K-means clustering, nonnegative matrix factorization, and latent Dirichlet allocation, can be written as probabilistic models underpinned by a low-rank matrix factorization. In addition to highlighting their similarities, this formulation clarifies the various assumptions and restrictions of each approach, which eases identifying the appropriate method for specific applications for applied medical researchers. We also touch upon the most important aspects of inference and model selection for the application of these methods to health data.

A topic modeling approach reveals the dynamic T cell composition of peripheral blood during cancer immunotherapy

We present TopicFlow, a computational framework for flow cytometry data analysis of patient blood samples for the identification of functional and dynamic topics in circulating T cell population. This framework applies a Latent Dirichlet Allocation (LDA) model, adapting the concept of topic modeling in text mining to flow cytometry. To demonstrate the utility of our method, we conducted an analysis of ∼17 million T cells collected from 138 peripheral blood samples in 51 patients with melanoma undergoing treatment with immune checkpoint inhibitors (ICIs). Our study highlights three latent dynamic topics identified by LDA: a T cell exhaustion topic that independently recapitulates the previously identified LAG-3+ immunotype associated with ICI resistance, a naive topic and its association with immune-related toxicity, and a T cell activation topic that emerges upon ICI treatment. Our approach can be broadly applied to mine high-parameter flow cytometry data for insights into mechanisms of treatment response and toxicity.

A mutation-level covariate model for mutational signatures

Mutational processes and their exposures in particular genomes are key to our understanding of how these genomes are shaped. However, current analyses assume that these processes are uniformly active across the genome without accounting for potential covariates such as strand or genomic region that could impact such activities. Here we suggest the first mutation-covariate models that explicitly model the effect of different covariates on the exposures of mutational processes. We apply these models to test the impact of replication strand on these processes and compare them to strand-oblivious models across a range of data sets. Our models capture replication strand specificity, point to signatures affected by it, and score better on held-out data compared to standard models that do not account for mutation-level covariate information.

Uncovering the hidden structure of dynamic T cell composition in peripheral blood during cancer immunotherapy: a topic modeling approach

Immune checkpoint inhibitors (ICIs), now mainstays in the treatment of cancer treatment, show great potential but only benefit a subset of patients. A more complete understanding of the immunological mechanisms and pharmacodynamics of ICI in cancer patients will help identify the patients most likely to benefit and will generate knowledge for the development of next-generation ICI regimens. We set out to interrogate the early temporal evolution of T cell populations from longitudinal single-cell flow cytometry data. We developed an innovative statistical and computational approach using a Latent Dirichlet Allocation (LDA) model that extends the concept of topic modeling used in text mining. This powerful unsupervised learning tool allows us to discover compositional topics within immune cell populations that have distinct functional and differentiation states and are biologically and clinically relevant. To illustrate the model's utility, we analyzed ∼17 million T cells obtained from 138 pre- and on-treatment peripheral blood samples from a cohort of melanoma patients treated with ICIs. We identified three latent dynamic topics: a T-cell exhaustion topic that recapitulates a LAG3+ predominant patient subgroup with poor clinical outcome; a naive topic that shows association with immune-related toxicity; and an immune activation topic that emerges upon ICI treatment. We identified that a patient subgroup with a high baseline of the naïve topic has a higher toxicity grade. While the current application is demonstrated using flow cytometry data, our approach has broader utility and creates a new direction for translating single-cell data into biological and clinical insights.

Cancer driver mutations: predictions and reality

Cancer cells accumulate many genetic alterations throughout their lifetime, but only a few of them drive cancer progression, termed driver mutations. Driver mutations may vary between cancer types and patients, can remain latent for a long time and become drivers at particular cancer stages, or may drive oncogenesis only in conjunction with other mutations. The high mutational, biochemical, and histological tumor heterogeneity makes driver mutation identification very challenging. In this review we summarize recent efforts to identify driver mutations in cancer and annotate their effects. We underline the success of computational methods to predict driver mutations in finding novel cancer biomarkers, including in circulating tumor DNA (ctDNA). We also report on the boundaries of their applicability in clinical research.

A generalizable machine learning framework for classifying DNA repair defects using ctDNA exomes

Specific classes of DNA damage repair (DDR) defect can drive sensitivity to emerging therapies for metastatic prostate cancer. However, biomarker approaches based on DDR gene sequencing do not accurately predict DDR deficiency or treatment benefit. Somatic alteration signatures may identify DDR deficiency but historically require whole-genome sequencing of tumour tissue. We assembled whole-exome sequencing data for 155 high ctDNA fraction plasma cell-free DNA and matched leukocyte DNA samples from patients with metastatic prostate or bladder cancer. Labels for DDR gene alterations were established using deep targeted sequencing. Per sample mutation and copy number features were used to train XGBoost ensemble models. Naive somatic features and trinucleotide signatures were associated with specific DDR gene alterations but insufficient to resolve each class. Conversely, XGBoost-derived models showed strong performance including an area under the curve of 0.99, 0.99 and 1.00 for identifying BRCA2, CDK12, and mismatch repair deficiency in metastatic prostate cancer. Our machine learning approach re-classified several samples exhibiting genomic features inconsistent with original labels, identified a metastatic bladder cancer sample with a homozygous BRCA2 copy loss, and outperformed an existing exome-based classifier for BRCA2 deficiency. We present DARC Sign (DnA Repair Classification SIGNatures); a public machine learning tool leveraging clinically-practical liquid biopsy specimens for simultaneously identifying multiple types of metastatic prostate cancer DDR deficiencies. We posit that it will be useful for understanding differential responses to DDR-directed therapies in ongoing clinical trials and may ultimately enable prospective identification of prostate cancers with phenotypic evidence of DDR deficiency.

mSigHdp: hierarchical Dirichlet process mixture modeling for mutational signature discovery

Mutational signatures are characteristic patterns of mutations caused by endogenous or exogenous mutational processes. These signatures can be discovered by analyzing mutations in large sets of samples-usually somatic mutations in tumor samples. Most programs for discovering mutational signatures are based on non-negative matrix factorization (NMF). Alternatively, signatures can be discovered using hierarchical Dirichlet process (HDP) mixture models, an approach that has been less explored. These models assign mutations to clusters and view each cluster as being generated from the signature of a particular mutational process. Here, we describe mSigHdp, an improved approach to using HDP mixture models to discover mutational signatures. We benchmarked mSigHdp and state-of-the-art NMF-based approaches on four realistic synthetic data sets. These data sets encompassed 18 cancer types. In total, they contained 3.5 × 10⁷ single-base-substitution mutations representing 32 signatures and 6.1 × 10⁶ small insertion and deletion mutations representing 13 signatures. For three of the four data sets, mSigHdp had the best positive predictive value for discovering mutational signatures, and for all four data sets, it had the best true positive rate. Its CPU usage was similar to that of the NMF-based approaches. Thus, mSigHdp is an important and practical addition to the set of tools available for discovering mutational signatures.

Ovarian cancer mutational processes drive site-specific immune evasion

High-grade serous ovarian cancer (HGSOC) is an archetypal cancer of genomic instability1-4 patterned by distinct mutational processes5,6, tumour heterogeneity7-9 and intraperitoneal spread7,8,10. Immunotherapies have had limited efficacy in HGSOC11-13, highlighting an unmet need to assess how mutational processes and the anatomical sites of tumour foci determine the immunological states of the tumour microenvironment. Here we carried out an integrative analysis of whole-genome sequencing, single-cell RNA sequencing, digital histopathology and multiplexed immunofluorescence of 160 tumour sites from 42 treatment-naive patients with HGSOC. Homologous recombination-deficient HRD-Dup (BRCA1 mutant-like) and HRD-Del (BRCA2 mutant-like) tumours harboured inflammatory signalling and ongoing immunoediting, reflected in loss of HLA diversity and tumour infiltration with highly differentiated dysfunctional CD8+ T cells. By contrast, foldback-inversion-bearing tumours exhibited elevated immunosuppressive TGFβ signalling and immune exclusion, with predominantly naive/stem-like and memory T cells. Phenotypic state associations were specific to anatomical sites, highlighting compositional, topological and functional differences between adnexal tumours and distal peritoneal foci. Our findings implicate anatomical sites and mutational processes as determinants of evolutionary phenotypic divergence and immune resistance mechanisms in HGSOC. Our study provides a multi-omic cellular phenotype data substrate from which to develop and interpret future personalized immunotherapeutic approaches and early detection research.

The genomic and immune landscape of long-term survivors of high-grade serous ovarian cancer

Fewer than half of all patients with advanced-stage high-grade serous ovarian cancers (HGSCs) survive more than five years after diagnosis, but those who have an exceptionally long survival could provide insights into tumor biology and therapeutic approaches. We analyzed 60 patients with advanced-stage HGSC who survived more than 10 years after diagnosis using whole-genome sequencing, transcriptome and methylome profiling of their primary tumor samples, comparing this data to 66 short- or moderate-term survivors. Tumors of long-term survivors were more likely to have multiple alterations in genes associated with DNA repair and more frequent somatic variants resulting in an increased predicted neoantigen load. Patients clustered into survival groups based on genomic and immune cell signatures, including three subsets of patients with BRCA1 alterations with distinctly different outcomes. Specific combinations of germline and somatic gene alterations, tumor cell phenotypes and differential immune responses appear to contribute to long-term survival in HGSC.

Single-cell genomic variation induced by mutational processes in cancer

How cell-to-cell copy number alterations that underpin genomic instability1 in human cancers drive genomic and phenotypic variation, and consequently the evolution of cancer2, remains understudied. Here, by applying scaled single-cell whole-genome sequencing3 to wild-type, TP53-deficient and TP53-deficient;BRCA1-deficient or TP53-deficient;BRCA2-deficient mammary epithelial cells (13,818 genomes), and to primary triple-negative breast cancer (TNBC) and high-grade serous ovarian cancer (HGSC) cells (22,057 genomes), we identify three distinct 'foreground' mutational patterns that are defined by cell-to-cell structural variation. Cell- and clone-specific high-level amplifications, parallel haplotype-specific copy number alterations and copy number segment length variation (serrate structural variations) had measurable phenotypic and evolutionary consequences. In TNBC and HGSC, clone-specific high-level amplifications in known oncogenes were highly prevalent in tumours bearing fold-back inversions, relative to tumours with homologous recombination deficiency, and were associated with increased clone-to-clone phenotypic variation. Parallel haplotype-specific alterations were also commonly observed, leading to phylogenetic evolutionary diversity and clone-specific mono-allelic expression. Serrate variants were increased in tumours with fold-back inversions and were highly correlated with increased genomic diversity of cellular populations. Together, our findings show that cell-to-cell structural variation contributes to the origins of phenotypic and evolutionary diversity in TNBC and HGSC, and provide insight into the genomic and mutational states of individual cancer cells.

Proteogenomic analysis of lung adenocarcinoma reveals tumor heterogeneity, survival determinants, and therapeutically relevant pathways

We present a deep proteogenomic profiling study of 87 lung adenocarcinoma (LUAD) tumors from the United States, integrating whole-genome sequencing, transcriptome sequencing, proteomics and phosphoproteomics by mass spectrometry, and reverse-phase protein arrays. We identify three subtypes from somatic genome signature analysis, including a transition-high subtype enriched with never smokers, a transversion-high subtype enriched with current smokers, and a structurally altered subtype enriched with former smokers, TP53 alterations, and genome-wide structural alterations. We show that within-tumor correlations of RNA and protein expression associate with tumor purity and immune cell profiles. We detect and independently validate expression signatures of RNA and protein that predict patient survival. Additionally, among co-measured genes, we found that protein expression is more often associated with patient survival than RNA. Finally, integrative analysis characterizes three expression subtypes with divergent mutations, proteomic regulatory networks, and therapeutic vulnerabilities. This proteogenomic characterization provides a foundation for molecularly informed medicine in LUAD.

Variant calling enhances the identification of cancer cells in single-cell RNA sequencing data

Single-cell RNA-sequencing is an invaluable research tool that allows for the investigation of gene expression in heterogeneous cancer cell populations in ways that bulk RNA-seq cannot. However, normal (i.e., non tumor) cells in cancer samples have the potential to confound the downstream analysis of single-cell RNA-seq data. Existing methods for identifying cancer and normal cells include copy number variation inference, marker-gene expression analysis, and expression-based clustering. This work aims to extend the existing approaches for identifying cancer cells in single-cell RNA-seq samples by incorporating variant calling and the identification of putative driver alterations. We found that putative driver alterations can be detected in single-cell RNA-seq data obtained with full-length transcript technologies and noticed that a subset of cells in tumor samples are enriched for putative driver alterations as compared to normal cells. Furthermore, we show that the number of putative driver alterations and inferred copy number variation are not correlated in all samples. Taken together, our findings suggest that augmenting existing cancer-cell filtering methods with variant calling and analysis can increase the number of tumor cells that can be confidently included in downstream analyses of single-cell full-length transcript RNA-seq datasets.

Accurate determination of CRISPR-mediated gene fitness in transplantable tumours

Assessing tumour gene fitness in physiologically-relevant model systems is challenging due to biological features of in vivo tumour regeneration, including extreme variations in single cell lineage progeny. Here we develop a reproducible, quantitative approach to pooled genetic perturbation in patient-derived xenografts (PDXs), by encoding single cell output from transplanted CRISPR-transduced cells in combination with a Bayesian hierarchical model. We apply this to 181 PDX transplants from 21 breast cancer patients. We show that uncertainty in fitness estimates depends critically on the number of transplant cell clones and the variability in clone sizes. We use a pathway-directed allelic series to characterize Notch signaling, and quantify TP53 / MDM2 drug-gene conditional fitness in outlier patients. We show that fitness outlier identification can be mirrored by pharmacological perturbation. Overall, we demonstrate that the gene fitness landscape in breast PDXs is dominated by inter-patient differences.

Comprehensive mutations analyses of FTO (fat mass and obesity-associated gene) and their effects on FTO’s substrate binding implicated in obesity

An excessive amount of fat deposition in the body leads to obesity which is a complex disease and poses a generic threat to human health. It increases the risk of various other diseases like diabetes, cardiovascular disease, and multiple types of cancer. Genomic studies have shown that the expression of the fat mass obesity (FTO) gene was highly altered and identified as one of the key biomarkers for obesity. This study has been undertaken to investigate the mutational profile of the FTO gene and elucidates its effect on the protein structure and function. Harmful effects of various missense mutations were predicted using different independent tools and it was observed that all mutations were highly pathogenic. Molecular dynamics (MD) simulations were performed to study the structure and function of FTO protein upon different mutations and it was found that mutations decreased the structure stability and affected protein conformation. Furthermore, a protein residue network analysis suggested that the mutations affected the overall residues bonding and topology. Finally, molecular docking coupled with MD simulation suggested that mutations affected FTO substrate binding by changing the protein-ligand affinity. Hence, the results of this finding would help in an in-depth understanding of the molecular biology of the FTO gene and its variants and lead to the development of effective therapeutics against associated diseases and disorders.

Multimodal data integration using machine learning improves risk stratification of high-grade serous ovarian cancer

Patients with high-grade serous ovarian cancer suffer poor prognosis and variable response to treatment. Known prognostic factors for this disease include homologous recombination deficiency status, age, pathological stage and residual disease status after debulking surgery. Recent work has highlighted important prognostic information captured in computed tomography and histopathological specimens, which can be exploited through machine learning. However, little is known about the capacity of combining features from these disparate sources to improve prediction of treatment response. Here, we assembled a multimodal dataset of 444 patients with primarily late-stage high-grade serous ovarian cancer and discovered quantitative features, such as tumor nuclear size on staining with hematoxylin and eosin and omental texture on contrast-enhanced computed tomography, associated with prognosis. We found that these features contributed complementary prognostic information relative to one another and clinicogenomic features. By fusing histopathological, radiologic and clinicogenomic machine-learning models, we demonstrate a promising path toward improved risk stratification of patients with cancer through multimodal data integration.

Computational analysis of cancer genome sequencing data

Distilling biologically meaningful information from cancer genome sequencing data requires comprehensive identification of somatic alterations using rigorous computational methods. As the amount and complexity of sequencing data have increased, so has the number of tools for analysing them. Here, we describe the main steps involved in the bioinformatic analysis of cancer genomes, review key algorithmic developments and highlight popular tools and emerging technologies. These tools include those that identify point mutations, copy number alterations, structural variations and mutational signatures in cancer genomes. We also discuss issues in experimental design, the strengths and limitations of sequencing modalities and methodological challenges for the future.

Spectrum of DNA mismatch repair failures viewed through the lens of cancer genomics and implications for therapy

Genome sequencing can be used to detect DNA repair failures in tumors and learn about underlying mechanisms. Here, we synthesize findings from genomic studies that examined deficiencies of the DNA mismatch repair (MMR) pathway. The impairment of MMR results in genome-wide hypermutation and in the 'microsatellite instability' (MSI) phenotype-occurrence of indel mutations at short tandem repeat (microsatellite) loci. The MSI status of tumors was traditionally assessed by molecular testing of a selected set of MS loci or by measuring MMR protein expression levels. Today, genomic data can provide a more complete picture of the consequences on genomic instability. Multiple computational studies examined somatic mutation distributions that result from failed DNA repair pathways in tumors. These include analyzing the commonly studied trinucleotide mutational spectra of single-nucleotide variants (SNVs), as well as of other features such as indels, structural variants, mutation clusters and regional mutation rate redistribution. The identified mutation patterns can be used to rigorously measure prevalence of MMR failures across cancer types, and potentially to subcategorize the MMR deficiencies. Diverse data sources, genomic and pre-genomic, from human and from experimental models, suggest there are different ways in which MMR can fail, and/or that the cell-type or genetic background may result in different types of MMR mutational patterns. The spectrum of MMR failures may direct cancer evolution, generating particular sets of driver mutations. Moreover, MMR affects outcomes of therapy by DNA damaging drugs, antimetabolites, nonsense-mediated mRNA decay (NMD) inhibitors, and immunotherapy by promoting either resistance or sensitivity, depending on the type of therapy.

Harnessing multimodal data integration to advance precision oncology

Advances in quantitative biomarker development have accelerated new forms of data-driven insights for patients with cancer. However, most approaches are limited to a single mode of data, leaving integrated approaches across modalities relatively underdeveloped. Multimodal integration of advanced molecular diagnostics, radiological and histological imaging, and codified clinical data presents opportunities to advance precision oncology beyond genomics and standard molecular techniques. However, most medical datasets are still too sparse to be useful for the training of modern machine learning techniques, and significant challenges remain before this is remedied. Combined efforts of data engineering, computational methods for analysis of heterogeneous data and instantiation of synergistic data models in biomedical research are required for success. In this Perspective, we offer our opinions on synthesizing complementary modalities of data with emerging multimodal artificial intelligence methods. Advancing along this direction will result in a reimagined class of multimodal biomarkers to propel the field of precision oncology in the coming decade.

Accuracy of mutational signature software on correlated signatures

Mutational signatures are characteristic patterns of mutations generated by exogenous mutagens or by endogenous mutational processes. Mutational signatures are important for research into DNA damage and repair, aging, cancer biology, genetic toxicology, and epidemiology. Unsupervised learning can infer mutational signatures from the somatic mutations in large numbers of tumors, and separating correlated signatures is a notable challenge for this task. To investigate which methods can best meet this challenge, we assessed 18 computational methods for inferring mutational signatures on 20 synthetic data sets that incorporated varying degrees of correlated activity of two common mutational signatures. Performance varied widely, and four methods noticeably outperformed the others: hdp (based on hierarchical Dirichlet processes), SigProExtractor (based on multiple non-negative matrix factorizations over resampled data), TCSM (based on an approach used in document topic analysis), and mutSpec.NMF (also based on non-negative matrix factorization). The results underscored the complexities of mutational signature extraction, including the importance and difficulty of determining the correct number of signatures and the importance of hyperparameters. Our findings indicate directions for improvement of the software and show a need for care when interpreting results from any of these methods, including the need for assessing sensitivity of the results to input parameters.

Gaining insights into relevance across cancers based on mutation features of TP53 gene

The tumor suppressor gene TP53, one of the most frequently mutated genes, is recognized as the guardian of genome and can provide a significant barrier to neoplastic transformation and tumor progression. Traditional theory believes that TP53 mutations are equal among cancer types. However, to date, no study has explored the TP53 mutation profile from a holistic and systematic standpoint to discovery its relevance and feature with cancers. Mutation signature, an unbiased approach to identify the mutational processes, can be a potent indicator for exploring mutation-driven tumor occurrence and progression. In this research, several features such as hotspots, mutability and mutation signature of somatic TP53 mutations derived from 18 types of cancer tissues from cBioPortal were analyzed and manifested the organizational preference among cancers. Mutation signatures found in almost all cancer types were Signature 6 related to mismatch repair deficiency, and Signature 1 that reflects the natural decomposition of 5-methylcytosine into thymine associated with aging. Meanwhile, several signatures of TP53 mutations displayed tissue-selective. Mutations enriched in bladder, skin, lung cancer were associated with signatures of APOBEC activity (Signature 2 and 13), alkylating agents (Signature 11), and tobacco smoke (Signature 4), respectively. Moreover, Signature 4 and 29 associated with tobacco smoking or chewing found in lung, sarcoma, esophageal, and head and neck cancer may be related to their smoking history. In addition, several digestive cancers, including colorectal, stomach, pancreatic and esophageal cancers, showed the high correlation in context and mutation signature profiles. Our study suggests that the tissue-selective activity of mutational processes would reflect the tissue-specific enrichment of TP53 mutations and provides a new perspective to understand the relevance of diverse diseases based on the spectrum of TP53 mutations.

•

TP53 mutations display the tissue-selective among cancers.

•

Mutation signature of TP53 reflects the tissue-selective activity of mutational processes among cancers.

•

Digestive cancers show the high correlation based on TP53 mutation signature.

Homologous Recombination Deficiency Assays in Epithelial Ovarian Cancer: Current Status and Future Direction

Epithelial ovarian cancer (EOC) patients are generally diagnosed at an advanced stage, usually relapse after initial treatments, which include debulking surgery and adjuvant platinum-based chemotherapy, and eventually have poor 5-year survival of less than 50%. In recent years, promising survival benefits from maintenance therapy with poly(ADP-ribose) polymerase (PARP) inhibitor (PARPi) has changed the management of EOC in newly diagnosed and recurrent disease. Identification of BRCA mutations and/or homologous recombination deficiency (HRD) is critical for selecting patients for PARPi treatment. However, the currently available HRD assays are not perfect predictors of the clinical response to PARPis in EOC patients. In this review, we introduce the concept of synthetic lethality, the rationale of using PARPi when HRD is present in tumor cells, the clinical trials of PARPi incorporating the HRD assays for EOC, the current HRD assays, and other HRD assays in development.

A mixture model for signature discovery from sparse mutation data

Mutational signatures are key to understanding the processes that shape cancer genomes, yet their analysis requires relatively rich whole-genome or whole-exome mutation data. Recently, orders-of-magnitude sparser gene-panel-sequencing data have become increasingly available in the clinic. To deal with such sparse data, we suggest a novel mixture model, Mix. In application to simulated and real gene-panel sequences, Mix is shown to outperform current approaches and yield mutational signatures and patient stratifications that are in higher agreement with the literature. We further demonstrate its utility in several clinical settings, successfully predicting therapy benefit and patient groupings from MSK-IMPACT pan-cancer data. Availability: https://github.com/itaysason/Mix-MMM .

Mutational Signatures: From Methods to Mechanisms

Mutations are the driving force of evolution, yet they underlie many diseases, in particular, cancer. They are thought to arise from a combination of stochastic errors in DNA processing, naturally occurring DNA damage (e.g., the spontaneous deamination of methylated CpG sites), replication errors, and dysregulation of DNA repair mechanisms. High-throughput sequencing has made it possible to generate large datasets to study mutational processes in health and disease. Since the emergence of the first mutational process studies in 2012, this field is gaining increasing attention and has already accumulated a host of computational approaches and biomedical applications.

Learning mutational signatures and their multidimensional genomic properties with TensorSignatures

We present TensorSignatures, an algorithm to learn mutational signatures jointly across different variant categories and their genomic localisation and properties. The analysis of 2778 primary and 3824 metastatic cancer genomes of the PCAWG consortium and the HMF cohort shows that all signatures operate dynamically in response to genomic states. The analysis pins differential spectra of UV mutagenesis found in active and inactive chromatin to global genome nucleotide excision repair. TensorSignatures accurately characterises transcription-associated mutagenesis in 7 different cancer types. The algorithm also extracts distinct signatures of replication- and double strand break repair-driven mutagenesis by APOBEC3A and 3B with differential numbers and length of mutation clusters. Finally, TensorSignatures reproduces a signature of somatic hypermutation generating highly clustered variants at transcription start sites of active genes in lymphoid leukaemia, distinct from a general and less clustered signature of Polη-driven translesion synthesis found in a broad range of cancer types. In summary, TensorSignatures elucidates complex mutational footprints by characterising their underlying processes with respect to a multitude of genomic variables.

ESMO recommendations on predictive biomarker testing for homologous recombination deficiency and PARP inhibitor benefit in ovarian cancer

BACKGROUND

Homologous recombination repair deficiency (HRD) is a frequent feature of high-grade serous ovarian, fallopian tube and peritoneal carcinoma (HGSC) and is associated with sensitivity to PARP inhibitor (PARPi) therapy. HRD testing provides an opportunity to optimise PARPi use in HGSC but methodologies are diverse and clinical application remains controversial.

MATERIALS AND METHODS

To define best practice for HRD testing in HGSC the ESMO Translational Research and Precision Medicine Working Group launched a collaborative project that incorporated a systematic review approach. The main aims were to (i) define the term 'HRD test'; (ii) provide an overview of the biological rationale and the level of evidence supporting currently available HRD tests; (iii) provide recommendations on the clinical utility of HRD tests in clinical management of HGSC.

RESULTS

A broad range of repair genes, genomic scars, mutational signatures and functional assays are associated with a history of HRD. Currently, the clinical validity of HRD tests in ovarian cancer is best assessed, not in terms of biological HRD status per se, but in terms of PARPi benefit. Clinical trials evidence supports the use of BRCA mutation testing and two commercially available assays that also incorporate genomic instability for identifying subgroups of HGSCs that derive different magnitudes of benefit from PARPi therapy, albeit with some variation by clinical scenario. These tests can be used to inform treatment selection and scheduling but their use is limited by a failure to consistently identify a subgroup of patients who derive no benefit from PARPis in most studies. Existing tests lack negative predictive value and inadequately address the complex and dynamic nature of the HRD phenotype.

CONCLUSIONS

Currently available HRD tests are useful for predicting likely magnitude of benefit from PARPis but better biomarkers are urgently needed to better identify current homologous recombination proficiency status and stratify HGSC management.

Mechanisms underlying genome instability mediated by formation of foldback inversions in Saccharomyces cerevisiae

Foldback inversions, also called inverted duplications, have been observed in human genetic diseases and cancers. Here, we used a Saccharomyces cerevisiae genetic system that generates gross chromosomal rearrangements (GCRs) mediated by foldback inversions combined with whole-genome sequencing to study their formation. Foldback inversions were mediated by formation of single-stranded DNA hairpins. Two types of hairpins were identified: small-loop hairpins that were suppressed by MRE11, SAE2, SLX1, and YKU80 and large-loop hairpins that were suppressed by YEN1, TEL1, SWR1, and MRC1. Analysis of CRISPR/Cas9-induced double strand breaks (DSBs) revealed that long-stem hairpin-forming sequences could form foldback inversions when proximal or distal to the DSB, whereas short-stem hairpin-forming sequences formed foldback inversions when proximal to the DSB. Finally, we found that foldback inversion GCRs were stabilized by secondary rearrangements, mostly mediated by different homologous recombination mechanisms including single-strand annealing; however, POL32-dependent break-induced replication did not appear to be involved forming secondary rearrangements.

The repertoire of mutational signatures in human cancer

Somatic mutations in cancer genomes are caused by multiple mutational processes, each of which generates a characteristic mutational signature1. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium2 of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we characterized mutational signatures using 84,729,690 somatic mutations from 4,645 whole-genome and 19,184 exome sequences that encompass most types of cancer. We identified 49 single-base-substitution, 11 doublet-base-substitution, 4 clustered-base-substitution and 17 small insertion-and-deletion signatures. The substantial size of our dataset, compared with previous analyses3-15, enabled the discovery of new signatures, the separation of overlapping signatures and the decomposition of signatures into components that may represent associated-but distinct-DNA damage, repair and/or replication mechanisms. By estimating the contribution of each signature to the mutational catalogues of individual cancer genomes, we revealed associations of signatures to exogenous or endogenous exposures, as well as to defective DNA-maintenance processes. However, many signatures are of unknown cause. This analysis provides a systematic perspective on the repertoire of mutational processes that contribute to the development of human cancer.

Why Are Some Driver Mutations Rare?

Understanding why driver mutations that promote cancer are sometimes rare is important for precision medicine since it would help in their identification. Driver mutations are largely discovered through their frequencies. Thus, rare mutations often escape detection. Unlike high-frequency drivers, low-frequency drivers can be tissue specific; rare drivers have extremely low frequencies. Here, we discuss rare drivers and strategies to discover them. We suggest that allosteric driver mutations shift the protein ensemble from the inactive to the active state. Rare allosteric drivers are statistically rare since, to switch the protein functional state, they cooperate with additional mutations, and these are not considered in the patient cancer-specific protein sequence analysis. A complete landscape of mutations that drive cancer will reveal tumor-specific therapeutic vulnerabilities.

Computed Tomography-Derived Radiomic Metrics Can Identify Responders to Immunotherapy in Ovarian Cancer

PURPOSE

To determine if radiomic measures of tumor heterogeneity derived from baseline contrast-enhanced computed tomography (CE-CT) are associated with durable clinical benefit and time to off-treatment in patients with recurrent ovarian cancer (OC) enrolled in prospective immunotherapeutic trials.

MATERIALS AND METHODS

This retrospective study included 75 patients with recurrent OC who were enrolled in prospective immunotherapeutic trials (n = 74) or treated off-label (n = 1) and had baseline CE-CT scans. Disease burden (total tumor volume, number of disease sites), radiomic measures of intertumor heterogeneity (cluster-site entropy, cluster-site dissimilarity), and intratumor heterogeneity of the largest lesion (Haralick texture features) were computed. Associations of clinical, conventional imaging, and radiomic measures with durable clinical benefit and time to off-treatment were examined.

RESULTS

In univariable analysis, fewer disease sites, lower intertumor heterogeneity (lower cluster-site entropy, lower cluster-site dissimilarity), and lower intratumor heterogeneity of the largest lesion (higher energy) were significantly associated with durable clinical benefit (P ≤ .031). More disease sites, presence of pleural disease and/or distant metastases, higher intertumor heterogeneity (higher cluster-site entropy, higher cluster-site dissimilarity), and higher intratumor heterogeneity of the largest lesion (higher Contrast_{largest-lesion}) were significantly associated with shorter time to off-treatment (P ≤ .034). In multivariable analysis, higher Energy_{largest-lesion} (indicator of lower intratumor heterogeneity; P = .006; odds ratio, 1.41) and fewer disease sites (P = .003; odds ratio, 1.64) remained significant indicators of durable clinical benefit (multivariable model C-index, 0.821). Higher cluster-site dissimilarity (indicator of higher intertumor heterogeneity) was a modest but single independent indicator of shorter time to off-treatment (P = .004; hazard ratio, 1.19; C-index, 0.6).

CONCLUSION

Fewer disease sites and lower intra- and intertumor heterogeneity modeled from the baseline CE-CT may indicate better response of OC to immunotherapy.