A mutational signature reveals alterations underlying deficient homologous recombination repair in breast cancer

A review on trends in development and translation of omics signatures in cancer

The field of cancer genomics and transcriptomics has evolved from targeted profiling to swift sequencing of individual tumor genome and transcriptome. The steady growth in genome, epigenome, and transcriptome datasets on a genome-wide scale has significantly increased our capability in capturing signatures that represent both the intrinsic and extrinsic biological features of tumors. These biological differences can help in precise molecular subtyping of cancer, predicting tumor progression, metastatic potential, and resistance to therapeutic agents. In this review, we summarized the current development of genomic, methylomic, transcriptomic, proteomic and metabolic signatures in the field of cancer research and highlighted their potentials in clinical applications to improve diagnosis, prognosis, and treatment decision in cancer patients.

(Single-stranded DNA) gaps in understanding BRCAness

The tumour-suppressive roles of BRCA1 and 2 have been attributed to three seemingly distinct functions - homologous recombination, replication fork protection, and single-stranded (ss)DNA gap suppression - and their relative importance is under debate. In this review, we examine the origin and resolution of ssDNA gaps and discuss the recent advances in understanding the role of BRCA1/2 in gap suppression. There are ample data showing that gap accumulation in BRCA1/2-deficient cells is linked to genomic instability and chemosensitivity. However, it remains unclear whether there is a causative role and the function of BRCA1/2 in gap suppression cannot unambiguously be dissected from their other functions. We therefore conclude that the three functions of BRCA1 and 2 are closely intertwined and not mutually exclusive.

The mutagenic consequences of defective DNA repair

Multiple separate repair mechanisms safeguard the genome against various types of DNA damage, and their failure can increase the rate of spontaneous mutagenesis. The malfunction of distinct repair mechanisms leads to genomic instability through different mutagenic processes. For example, defective mismatch repair causes high base substitution rates and microsatellite instability, whereas homologous recombination deficiency is characteristically associated with deletions and chromosome instability. This review presents a comprehensive collection of all mutagenic phenotypes associated with the loss of each DNA repair mechanism, drawing on data from a variety of model organisms and mutagenesis assays, and placing greatest emphasis on systematic analyses of human cancer datasets. We describe the latest theories on the mechanism of each mutagenic process, often explained by reliance on an alternative repair pathway or the error-prone replication of unrepaired, damaged DNA. Aided by the concept of mutational signatures, the genomic phenotypes can be used in cancer diagnosis to identify defective DNA repair pathways.

Luminal androgen receptor subtype and tumor-infiltrating lymphocytes groups based on triple-negative breast cancer molecular subclassification

In our previous study, we developed a triple-negative breast cancer (TNBC) subtype classification that correlated with the TNBC molecular subclassification. In this study, we aimed to evaluate the predictor variables of this subtype classification on the whole slide and to validate the model's performance by using an external test set. We explored the characteristics of this subtype classification and investigated genomic alterations, including genomic scar signature scores. First, TNBC was classified into the luminal androgen receptor (LAR) and non-luminal androgen receptor (non-LAR) subtypes based on the AR Allred score (≥ 6 and < 6, respectively). Then, the non-LAR subtype was further classified into the lymphocyte-predominant (LP), lymphocyte-intermediate (LI), and lymphocyte-depleted (LD) groups based on stromal tumor-infiltrating lymphocytes (TILs) (< 20%, > 20% but < 60%, and ≥ 60%, respectively). This classification showed fair agreement with the molecular classification in the test set. The LAR subtype was characterized by a high rate of PIK3CA mutation, CD274 (encodes PD-L1) and PDCD1LG2 (encodes PD-L2) deletion, and a low homologous recombination deficiency (HRD) score. The non-LAR LD TIL group was characterized by a high frequency of NOTCH2 and MYC amplification and a high HRD score.

Ataxia-Telangiectasia Mutated Loss-of-Function Displays Variant and Tissue-Specific Differences across Tumor Types

PURPOSE

Mutations in the ATM gene are common in multiple cancers, but clinical studies of therapies targeting ATM-aberrant cancers have yielded mixed results. Refinement of ATM loss of function (LOF) as a predictive biomarker of response is urgently needed.

EXPERIMENTAL DESIGN

We present the first disclosure and preclinical development of a novel, selective ATR inhibitor, ART0380, and test its antitumor activity in multiple preclinical cancer models. To refine ATM LOF as a predictive biomarker, we performed a comprehensive pan-cancer analysis of ATM variants in patient tumors and then assessed the ATM variant-to-protein relationship. Finally, we assessed a novel ATM LOF biomarker approach in retrospective clinical data sets of patients treated with platinum-based chemotherapy or ATR inhibition.

RESULTS

ART0380 had potent, selective antitumor activity in a range of preclinical cancer models with differing degrees of ATM LOF. Pan-cancer analysis identified 10,609 ATM variants in 8,587 patient tumors. Cancer lineage-specific differences were seen in the prevalence of deleterious (Tier 1) versus unknown/benign (Tier 2) variants, selective pressure for loss of heterozygosity, and concordance between a deleterious variant and ATM loss of protein (LOP). A novel ATM LOF biomarker approach that accounts for variant classification, relationship to ATM LOP, and tissue-specific penetrance significantly enriched for patients who benefited from platinum-based chemotherapy or ATR inhibition.

CONCLUSIONS

These data help to better define ATM LOF across tumor types in order to optimize patient selection and improve molecularly targeted therapeutic approaches for patients with ATM LOF cancers.

On the Hunt for the Missed Genetic Causes of Multiple Primary Tumors

Improved cancer screening and treatment programs have led to an increased survivorship of patients with cancer, but consequently also to the rise in number of individuals with multiple primary tumors (MPT). Germline testing is the first approach investigating the cause of MPT, as a positive result provides a diagnosis and proper clinical management to the affected individual and their family. Negative or inconclusive genetic results could suggest non-genetic causes, but are negative genetic results truly negative? Herein, we discuss the potential sources of missed genetic causes and highlight the trove of knowledge MPT can provide. See related article by Borja et al., p. 209.

Clinical evaluation of a low-coverage whole-genome test for detecting homologous recombination deficiency in ovarian cancer

BACKGROUND

The PAOLA-1/ENGOT-ov25 trial showed that maintenance olaparib plus bevacizumab increases survival of advanced ovarian cancer patients with homologous recombination deficiency (HRD). However, decentralized solutions to test for HRD in clinical routine are scarce. The goal of this study was to retrospectively validate on tumor samples from the PAOLA-1 trial, the decentralized SeqOne assay, which relies on shallow Whole Genome Sequencing (sWGS) to capture genomic instability and targeted sequencing to determine BRCA status.

METHODS

The study comprised 368 patients from the PAOLA-1 trial. The SeqOne assay was compared to the Myriad MyChoice HRD test (Myriad Genetics), and results were analyzed with respect to Progression-Free Survival (PFS).

RESULTS

We found a 95% concordance between the HRD status of the two tests (95% Confidence Interval (CI); 92%-97%). The Positive Percentage Agreement (PPA) of the sWGS test was 95% (95% CI; 91%-97%) like its Negative Percentage Agreement (NPA) (95% CI; 89%-98%). In patients with HRD-positive tumors treated with olaparib plus bevacizumab, the PFS Hazard Ratio (HR) was 0.38 (95% CI; 0.26-0.54) with SeqOne assay and 0.32 (95% CI; 0.22-0.45) with the Myriad assay. In patients with HRD-negative tumors, HR was 0.99 (95% CI; 0.68-1.42) and 1.05 (95% CI; 0.70-1.57) with SeqOne and Myriad assays. Among patients with BRCA-wildtype tumors, those with HRD-positive tumors, benefited from olaparib plus bevacizumab maintenance, with HR of 0.48 (95% CI: 0.29-0.79) and of 0.38 (95% CI: 0.23 to 0.63) with the SeqOne and Myriad assay.

CONCLUSION

The SeqOne assay offers a clinically validated approach to detect HRD.

Diffsig: Associating Risk Factors with Mutational Signatures

BACKGROUND

Somatic mutational signatures elucidate molecular vulnerabilities to therapy, and therefore detecting signatures and classifying tumors with respect to signatures has clinical value. However, identifying the etiology of the mutational signatures remains a statistical challenge, with both small sample sizes and high variability in classification algorithms posing barriers. As a result, few signatures have been strongly linked to particular risk factors.

METHODS

Here, we develop a statistical model, Diffsig, for estimating the association of one or more continuous or categorical risk factors with DNA mutational signatures. Diffsig takes into account the uncertainty associated with assigning signatures to samples as well as multiple risk factors' simultaneous effect on observed DNA mutations.

RESULTS

We applied Diffsig to breast cancer data to assess relationships between five established breast-relevant mutational signatures and etiologic variables, confirming known mechanisms of cancer development. In simulation, our model was capable of accurately estimating expected associations in a variety of contexts.

CONCLUSIONS

Diffsig allows researchers to quantify and perform inference on the associations of risk factors with mutational signatures.

IMPACT

We expect Diffsig to provide more robust associations of risk factors with signatures to lead to better understanding of the tumor development process and improved models of tumorigenesis.

Exploration of poly (ADP-ribose) polymerase inhibitor resistance in the treatment of BRCA1/2-mutated cancer

Breast cancer susceptibility 1/2 (BRCA1/2) genes play a crucial role in DNA damage repair, yet mutations in these genes increase the susceptibility to tumorigenesis. Exploiting the synthetic lethality mechanism between BRCA1/2 mutations and poly(ADP-ribose) polymerase (PARP) inhibition has led to the development and clinical approval of PARP inhibitor (PARPi), representing a milestone in targeted therapy for BRCA1/2 mutant tumors. This approach has paved the way for leveraging synthetic lethality in tumor treatment strategies. Despite the initial success of PARPis, resistance to these agents diminishes their efficacy in BRCA1/2-mutant tumors. Investigations into PARPi resistance have identified replication fork stability and homologous recombination repair as key factors sensitive to PARPis. Additionally, studies suggest that replication gaps may also confer sensitivity to PARPis. Moreover, emerging evidence indicates a correlation between PARPi resistance and cisplatin resistance, suggesting a potential overlap in the mechanisms underlying resistance to both agents. Given these findings, it is imperative to explore the interplay between replication gaps and PARPi resistance, particularly in the context of platinum resistance. Understanding the impact of replication gaps on PARPi resistance may offer insights into novel therapeutic strategies to overcome resistance mechanisms and enhance the efficacy of targeted therapies in BRCA1/2-mutant tumors.

PARP Inhibitors for Breast Cancer Treatment: A Review

Importance

Poly(adenosine diphosphate-ribose) polymerase (PARP) inhibitors have revolutionized the treatment of patients with germline BRCA1/2-associated breast cancer, representing the first targeted therapy capable of improving outcomes in patients with hereditary tumors. However, resistance to PARP inhibitors occurs in almost all patients.

Observations

This narrative review summarizes the biological rationale behind the use of PARP inhibitors in breast cancer, as well as the available evidence, recent progress, and potential future applications of these agents. Recent studies have shown that the benefit of PARP inhibitors extends beyond patients with germline BRCA1/2-associated metastatic breast cancer to patients with somatic BRCA1/2 variants and to those with germline PALB2 alterations. Moreover, these agents proved to be effective both in the metastatic and adjuvant settings. However, patients with metastatic breast cancer usually do not achieve the long-term benefit from PARP inhibitors observed in other tumor types. Mechanisms of resistance have been identified, but how to effectively target them is largely unknown. Ongoing research is investigating both novel therapeutics and new combination strategies to overcome resistance. PARP1-selective inhibitors, by sparing the hematological toxic effects induced by the PARP2 blockade, are promising agents to be combined with chemotherapy, antibody-drug conjugates, and other targeted therapies.

Conclusions and Relevance

Although the efficacy of PARP inhibitors is well established, many questions persist. Future research should focus on identifying predictive biomarkers and therapeutic strategies to overcome resistance. Integrating well-designed translational efforts into all clinical studies is thereby crucial to laying the groundwork for future insights from ongoing research.

A glycolytic metabolite bypasses "two-hit" tumor suppression by BRCA2

Knudson's "two-hit" paradigm posits that carcinogenesis requires inactivation of both copies of an autosomal tumor suppressor gene. Here, we report that the glycolytic metabolite methylglyoxal (MGO) transiently bypasses Knudson's paradigm by inactivating the breast cancer suppressor protein BRCA2 to elicit a cancer-associated, mutational single-base substitution (SBS) signature in nonmalignant mammary cells or patient-derived organoids. Germline monoallelic BRCA2 mutations predispose to these changes. An analogous SBS signature, again without biallelic BRCA2 inactivation, accompanies MGO accumulation and DNA damage in Kras-driven, Brca2-mutant murine pancreatic cancers and human breast cancers. MGO triggers BRCA2 proteolysis, temporarily disabling BRCA2's tumor suppressive functions in DNA repair and replication, causing functional haploinsufficiency. Intermittent MGO exposure incites episodic SBS mutations without permanent BRCA2 inactivation. Thus, a metabolic mechanism wherein MGO-induced BRCA2 haploinsufficiency transiently bypasses Knudson's two-hit requirement could link glycolysis activation by oncogenes, metabolic disorders, or dietary challenges to mutational signatures implicated in cancer evolution.

Human epidermal growth factor receptor-2 expression and subsequent dynamic changes in patients with ovarian cancer

Human epidermal growth factor receptor-2 (HER2)-targeting drugs are increasingly being incorporated into therapeutic paradigms for non-breast cancers, yet studies on HER2 expression in ovarian cancer (OC) are inadequate. Here, we studied the HER2 status and dynamic changes in OC by reviewing the records of patients who underwent HER2 testing at a single institution. Clinical parameters, including histology, BRCA status, and immunohistochemistry (IHC), were evaluated alongside HER2 expression, timing, and anatomical location. Among 200 patients, 28% and 6% exhibited expression scores of 2+ and 3+, respectively. HER2 3+ scores were observed in 23%, 11%, 9%, and 5% of mucinous, endometrioid, clear cell, and high-grade serous tumors, respectively, and were exclusively identified in BRCA-wildtype, mismatch repair-proficient, or PD-L1-low-expressing tumors. The TP53 mutation rate was low, whereas ARID1A, KRAS, and PIK3CA mutations were relatively more prevalent with HER2 scores of 2+ or 3+ than with 0 or 1+. Four of the five tumors with an HER2 3+ score exhibited ERBB2 amplification. Among 19 patients who underwent multiple time-lagged biopsies, 11 showed increased HER2 expression in subsequent biopsies. Patients with HER2-overexpressing OC exhibited distinct histological, IHC, and genomic profiles. HER2-targeting agents are potential options for BRCA-wildtype patients, particularly as later lines of treatment.

Post-therapy emergence of an NBN reversion mutation in a patient with pancreatic acinar cell carcinoma

Pancreatic acinar cell carcinoma (PACC) is a rare form of pancreatic cancer that commonly harbors targetable alterations, including activating fusions in the MAPK pathway and loss-of-function (LOF) alterations in DNA damage response/homologous recombination DNA repair-related genes. Here, we describe a patient with PACC harboring both somatic biallelic LOF of NBN and an activating NTRK1 fusion. Upon disease progression following 13 months of treatment with folinic acid, fluorouracil, irinotecan, and oxaliplatin (FOLFIRINOX), genomic analysis of a metastatic liver biopsy revealed the emergence of a novel reversion mutation restoring the reading frame of NBN. To our knowledge, genomic reversion of NBN has not been previously reported as a resistance mechanism in any tumor type. The patient was treated with, but did not respond to, targeted treatment with a selective NTRK inhibitor. This case highlights the complex but highly actionable genomic landscape of PACC and underlines the value of genomic profiling of rare tumor types such as PACC.

Free drug and ROS-responsive nanoparticle delivery of synergistic doxorubicin and olaparib combinations to triple negative breast cancer models

Combinations of the topoisomerase II inhibitor doxorubicin and the poly (ADP-ribose) polymerase inhibitor olaparib offer potential drug-drug synergy for the treatment of triple negative breast cancers (TNBC). In this study we performed in vitro screening of combinations of these drugs, administered directly or encapsulated within polymer nanoparticles, in both 2D and in 3D spheroid models of breast cancer. A variety of assays were used to evaluate drug potency, and calculations of combination index (CI) values indicated that synergistic effects of drug combinations occurred in a molar-ratio dependent manner. It is suggested that the mechanisms of synergy were related to enhancement of DNA damage as shown by the level of double-strand DNA breaks, and mechanisms of antagonism associated with mitochondrial mediated cell survival, as indicated by reactive oxygen species (ROS) generation. Enhanced drug delivery and potency was observed with nanoparticle formulations, with a greater extent of doxorubicin localised to cell nuclei as evidenced by microscopy, and higher cytotoxicity at the same time points compared to free drugs. Together, the work presented identifies specific combinations of doxorubicin and olaparib which were most effective in a panel of TNBC cell lines, explores the mechanisms by which these combined agents might act, and shows that formulation of these drug combinations into polymeric nanoparticles at specific ratios conserves synergistic action and enhanced potency in vitro compared to the free drugs.

A deep-learning-based genomic status estimating framework for homologous recombination deficiency detection from low-pass whole genome sequencing

Genome-wide sequencing allows for prediction of clinical treatment responses and outcomes by estimating genomic status. Here, we developed Genomic Status scan (GSscan), a long short-term memory (LSTM)-based deep-learning framework, which utilizes low-pass whole genome sequencing (WGS) data to capture genomic instability-related features. In this study, GSscan directly surveys homologous recombination deficiency (HRD) status independent of other existing biomarkers. In breast cancer, GSscan achieved an AUC of 0.980 in simulated low-pass WGS data, and obtained a higher HRD risk score in clinical BRCA-deficient breast cancer samples (p = 1.3 × 10-4, compared with BRCA-intact samples). In ovarian cancer, GSscan obtained higher HRD risk scores in BRCA-deficient samples in both simulated data and clinical samples (p = 2.3 × 10-5 and p = 0.039, respectively, compared with BRCA-intact samples). Moreover, HRD-positive patients predicted by GSscan showed longer progression-free intervals in TCGA datasets (p = 0.0011) treated with platinum-based adjuvant chemotherapy, outperforming existing low-pass WGS-based methods. Furthermore, GSscan can accurately predict HRD status using only 1 ng of input DNA and a minimum sequencing coverage of 0.02 × , providing a reliable, accessible, and cost-effective approach. In summary, GSscan effectively and accurately detected HRD status, and provide a broadly applicable framework for disease diagnosis and selecting appropriate disease treatment.

Functional Homologous Recombination (HR) Screening Shows the Majority of BRCA1/2-Mutant Breast and Ovarian Cancer Cell Lines Are HR-Proficient

Tumors with a pathogenic BRCA1/2 mutation are homologous recombination (HR)-deficient (HRD) and consequently sensitive to platinum-based chemotherapy and Poly-[ADP-Ribose]-Polymerase inhibitors (PARPi). We hypothesized that functional HR status better reflects real-time HR status than BRCA1/2 mutation status. Therefore, we determined the functional HR status of 53 breast cancer (BC) and 38 ovarian cancer (OC) cell lines by measuring the formation of RAD51 foci after irradiation. Discrepancies between functional HR and BRCA1/2 mutation status were investigated using exome sequencing, methylation and gene expression data from 50 HR-related genes. A pathogenic BRCA1/2 mutation was found in 10/53 (18.9%) of BC and 7/38 (18.4%) of OC cell lines. Among BRCA1/2-mutant cell lines, 14/17 (82.4%) were HR-proficient (HRP), while 1/74 (1.4%) wild-type cell lines was HRD. For most (80%) cell lines, we explained the discrepancy between functional HR and BRCA1/2 mutation status. Importantly, 12/14 (85.7%) BRCA1/2-mutant HRP cell lines were explained by mechanisms directly acting on BRCA1/2. Finally, functional HR status was strongly associated with COSMIC single base substitution signature 3, but not BRCA1/2 mutation status. Thus, the majority of BRCA1/2-mutant cell lines do not represent a suitable model for HRD. Moreover, exclusively determining BRCA1/2 mutation status may not suffice for platinum-based chemotherapy or PARPi patient selection.

Positive and negative regulators of RAD51/DMC1 in homologous recombination and DNA replication

RAD51 recombinase plays a central role in homologous recombination (HR) by forming a nucleoprotein filament on single-stranded DNA (ssDNA) to catalyze homology search and strand exchange between the ssDNA and a homologous double-stranded DNA (dsDNA). The catalytic activity of RAD51 assembled on ssDNA is critical for the DNA-homology-mediated repair of DNA double-strand breaks in somatic and meiotic cells and restarting stalled replication forks during DNA replication. The RAD51-ssDNA complex also plays a structural role in protecting the regressed/reversed replication fork. Two types of regulators control RAD51 filament formation, stability, and dynamics, namely positive regulators, including mediators, and negative regulators, so-called remodelers. The appropriate balance of action by the two regulators assures genome stability. This review describes the roles of positive and negative RAD51 regulators in HR and DNA replication and its meiosis-specific homolog DMC1 in meiotic recombination. We also provide future study directions for a comprehensive understanding of RAD51/DMC1-mediated regulation in maintaining and inheriting genome integrity.

Drug repositioning with adaptive graph convolutional networks

MOTIVATION

Drug repositioning is an effective strategy to identify new indications for existing drugs, providing the quickest possible transition from bench to bedside. With the rapid development of deep learning, graph convolutional networks (GCNs) have been widely adopted for drug repositioning tasks. However, prior GCNs based methods exist limitations in deeply integrating node features and topological structures, which may hinder the capability of GCNs.

RESULTS

In this study, we propose an adaptive GCNs approach, termed AdaDR, for drug repositioning by deeply integrating node features and topological structures. Distinct from conventional graph convolution networks, AdaDR models interactive information between them with adaptive graph convolution operation, which enhances the expression of model. Concretely, AdaDR simultaneously extracts embeddings from node features and topological structures and then uses the attention mechanism to learn adaptive importance weights of the embeddings. Experimental results show that AdaDR achieves better performance than multiple baselines for drug repositioning. Moreover, in the case study, exploratory analyses are offered for finding novel drug-disease associations.

AVAILABILITY AND IMPLEMENTATION

The soure code of AdaDR is available at: https://github.com/xinliangSun/AdaDR.

HRD related signature 3 predicts clinical outcome in advanced tubo-ovarian high-grade serous carcinoma

OBJECTIVES

We evaluated usability of single base substitution signature 3 (Sig3) as a biomarker for homologous recombination deficiency (HRD) in tubo-ovarian high-grade serous carcinoma (HGSC).

MATERIALS AND METHODS

This prospective observational trial includes 165 patients with advanced HGSC. Fresh tissue samples (n = 456) from multiple intra-abdominal areas at diagnosis and after neoadjuvant chemotherapy (NACT) were collected for whole-genome sequencing. Sig3 was assessed by fitting samples independently with COSMIC v3.2 reference signatures. An HR scar assay was applied for comparison. Progression-free survival (PFS) and overall survival (OS) were studied using Kaplan-Meier and Cox regression analysis.

RESULTS

Sig3 has a bimodal distribution, eliminating the need for an arbitrary cutoff typical in HR scar tests. Sig3 could be assessed from samples with low (10%) cancer cell proportion and was consistent between multiple samples and stable during NACT. At diagnosis, 74 (45%) patients were HRD (Sig3+), while 91 (55%) were HR proficient (HRP, Sig3-). Sig3+ patients had longer PFS and OS than Sig3- patients (22 vs. 13 months and 51 vs. 34 months respectively, both p < 0.001). Sig3 successfully distinguished the poor prognostic HRP group among BRCAwt patients (PFS 19 months for Sig3+ and 13 months for Sig3- patients, p < 0.001). However, Sig3 at diagnosis did not predict chemoresponse anymore in the first relapse. The patient-level concordance between Sig3 and HR scar assay was 87%, and patients with HRD according to both tests had the longest median PFS.

CONCLUSIONS

Sig3 is a prognostic marker in advanced HGSC and useful tool in patient stratification for HRD.

BRCAness of brain lesions reflects a worse outcome for patients with metastatic breast cancer

PURPOSE

Breast cancer often metastasizes to the central nervous system. Although the prognosis of brain metastases from breast cancer has been considered poor, and systemic therapy has not contributed to an improved prognosis, newer agents are expected to be more effective. BRCAness is defined as the status of homologous recombination deficiency (HRD) in tumor tissue, regardless of the presence of pathogenic germline BRCA1/2 variants. A study employing next-generation sequencing analysis showed that HRD was found relatively frequently in brain metastases of breast cancer patients. However, there have been no studies evaluating BRCAness in brain metastases of breast cancer with more efficient, rapid, and cost-effective methods.

METHODS

We retrospectively investigated 17 brain metastases of breast cancer that were surgically resected at our hospital from January 2007 to December 2022. Of these, samples from 15 patients were evaluable for BRCAness by employing multiplex ligation-dependent probe amplification (MLPA) assay.

RESULTS

Of the 15 patients, five patients (33%) had tumors with BRCAness. Clinicopathological factors of patients with brain metastases with BRCAness were not statistically different from those of patients who possessed tumors without BRCAness. Patients with brain metastases with BRCAness had shorter overall survival compared to those without BRCAness (BRCAness, median 15 months (95% CI 2-30) vs. non-BRCAness, median 28.5 months (95% CI 10-60); P = 0.013).

CONCLUSION

In this study, we evaluated BRCAness in brain metastases of breast cancer with the MLPA method, and found that about one-third of patients had BRCAness-positive tumors. The analysis of BRCAness using MLPA has the potential for practical clinical use.

Performance of a RAD51-based functional HRD test on paraffin-embedded breast cancer tissue

PURPOSE

BRCA-deficient breast cancers (BC) are highly sensitive to platinum-based chemotherapy and PARP inhibitors due to their deficiency in the homologous recombination (HR) pathway. However, HR deficiency (HRD) extends beyond BRCA-associated BC, highlighting the need for a sensitive method to enrich for HRD tumors in an alternative way. A promising approach is the use of functional HRD tests which evaluate the HR capability of tumor cells by measuring RAD51 protein accumulation at DNA damage sites. This study aims to evaluate the performance of a functional RAD51-based HRD test for the identification of HRD BC.

METHODS

The functional HR status of 63 diagnostic formalin-fixed paraffin-embedded (FFPE) BC samples was determined by applying the RAD51-FFPE test. Samples were screened for the presence of (epi)genetic defects in HR and matching tumor samples were analyzed with the RECAP test, which requires ex vivo irradiated fresh tumor tissue on the premise that the HRD status as determined by the RECAP test faithfully represented the functional HR status.

RESULTS

The RAD51-FFPE test identified 23 (37%) of the tumors as HRD, including three tumors with pathogenic variants in BRCA1/2. The RAD51-FFPE test showed a sensitivity of 88% and a specificity of 76% in determining the HR-class as defined by the RECAP test.

CONCLUSION

Given its high sensitivity and compatibility with FFPE samples, the RAD51-FFPE test holds great potential to enrich for HRD tumors, including those associated with BRCA-deficiency. This potential extends to situations where DNA-based testing may be challenging or not easily accessible in routine clinical practice. This is particularly important considering the potential implications for treatment decisions and patient stratification.

Genomic Profiling Reveals Germline Predisposition and Homologous Recombination Deficiency in Pancreatic Acinar Cell Carcinoma

PURPOSE

To determine the genetic predisposition underlying pancreatic acinar cell carcinoma (PACC) and characterize its genomic features.

METHODS

Both somatic and germline analyses were performed using an Food and Drug Administration-authorized matched tumor/normal sequencing assay on a clinical cohort of 28,780 patients with cancer, 49 of whom were diagnosed with PACC. For a subset of PACCs, whole-genome sequencing (WGS; n = 12) and RNA sequencing (n = 6) were performed.

RESULTS

Eighteen of 49 (36.7%) PACCs harbored germline pathogenic variants in homologous recombination (HR) and DNA damage response (DDR) genes, including BRCA1 (n = 1), BRCA2 (n = 12), PALB2 (n = 2), ATM (n = 2), and CHEK2 (n = 1). Thirty-one PACCs displayed pure, and 18 PACCs harbored mixed acinar cell histology. Fifteen of 31 (48%) pure PACCs harbored a germline pathogenic variant affecting HR-/DDR-related genes. BRCA2 germline pathogenic variants (11 of 31, 35%) were significantly more frequent in pure PACCs than in pancreatic adenocarcinoma (86 of 2,739, 3.1%; P < .001), high-grade serous ovarian carcinoma (67 of 1,318, 5.1%; P < .001), prostate cancer (116 of 3,401, 3.4%; P < .001), and breast cancer (79 of 3,196, 2.5%; P < .001). Genomic features of HR deficiency (HRD) were detected in 7 of 12 PACCs undergoing WGS, including 100% (n = 6) of PACCs with germline HR-related pathogenic mutations and 1 of 6 PACCs lacking known pathogenic alterations in HR-related genes. Exploratory analyses revealed that in PACCs, the repertoire of somatic driver genetic alterations and the load of neoantigens with high binding affinity varied according to the presence of germline pathogenic alterations affecting HR-/DDR-related genes and/or HRD.

CONCLUSION

In a large pan-cancer cohort, PACC was identified as the cancer type with the highest prevalence of both BRCA2 germline pathogenic variants and genomic features of HRD, suggesting that PACC should be considered as part of the spectrum of BRCA-related malignancies.

Clinical utility of genomic signatures for the management of early and metastatic triple-negative breast cancer

PURPOSE OF REVIEW

This comprehensive review aims to provide timely and relevant insights into the current therapeutic landscape for triple-negative breast cancer (TNBC) and the molecular features underlying this subtype. It emphasizes the need for more reliable biomarkers to refine prognostication and optimize therapy, considering the aggressive nature of TNBC and its limited targeted treatment options.

RECENT FINDINGS

The review explores the multidisciplinary management of early TNBC, which typically involves systemic chemotherapy, surgery, and radiotherapy. It highlights the emergence of immune checkpoint inhibitors (ICIs), poly(ADP-ribose) polymerase (PARP) inhibitors, and antibody-drug conjugates (ADCs) as promising therapeutic strategies for TNBC. Recent clinical trials investigating the use of ICIs in combination with chemotherapy and the approval of pembrolizumab and atezolizumab for PD-L1-positive metastatic TNBC are discussed. The efficacy of PARP inhibitors and ADCs in treating TNBC patients with specific genetic alterations is also highlighted.

SUMMARY

The findings discussed in this review have significant implications for clinical practice and research in TNBC. The identification of distinct molecular subtypes through gene expression profiling has enabled a better understanding of TNBC heterogeneity and its clinical implications. This knowledge has the potential to guide treatment decisions, as different subtypes display varying responses to neoadjuvant chemotherapy. Furthermore, the review emphasizes the importance of developing reliable genomic and transcriptomic signatures as biomarkers to refine patient prognostication and optimize therapy selection in TNBC. Integrating these signatures into clinical practice may lead to more personalized treatment approaches, improving outcomes for TNBC patients.

ShallowHRD status acts as an effective prognostic predictor in ovarian cancer patients treated by poly (ADP-ribose) polymerase inhibitors (PARPis)

PURPOSE

Homologous recombination deficiency (HRD) plays a crucial role in ovarian cancer patients who are treated with Poly (ADP-ribose) polymerase inhibitors (PARPis). It could be defined as a prognosis biomarker. However, many high throughput sequencing methods for evaluating HRD, including HRDetect (WGS 10X), SigMA (WGS 40X or panel 1000X), and scarHRD (WGS 30X), are technically complex, time and data-storage consuming, and costly. Herein, we aimed to develop a low-cost method by low sequencing coverage to identify HRD status for precision medication.

METHODS

We utilized ShallowHRD, a software tool to evaluate tumor HRD based on whole genome sequencing (WGS) at low coverage (1X), and established a novel scoring system, ShallowHRD score system.

RESULTS

Compared with negative ShallowHRD status (ShallowHRD score < 15 or BRCAwild), positive ShallowHRD status (ShallowHRD score ≥ 15 or BRCAmut) presented favorable survival after being treated with PARPis.

CONCLUSION

The ShallowHRD status is a good biomarker for predicting prognosis, which could help guide the clinical application of PARPis in ovarian cancer patients by a cost-effective, time and data-storage saving method.

Shallow whole genome sequencing approach to detect Homologous Recombination Deficiency in the PAOLA-1/ENGOT-OV25 phase-III trial

The bevacizumab (bev)/olaparib (ola) maintenance regimen was approved for BRCA1/2-mutated (BRCAmut) and Homologous Recombination Deficient (HRD) high-grade Advanced Ovarian Cancer (AOC) first line setting, based on a significantly improved progression-free survival (PFS) compared to bev alone in the PAOLA-1/ENGOT-ov25 trial (NCT02477644), where HRD was detected by MyChoice CDx PLUS test. The academic shallowHRDv2 test was developed based on shallow whole-genome sequencing as an alternative to MyChoice. Analytical and clinical validities of shallowHRDv2 as compared to MyChoice on 449 PAOLA-1 tumor samples are presented. The overall agreement between shallowHRDv2 and MyChoice was 94% (369/394). Less non-contributive tests were observed with shallowHRDv2 (15/449; 3%) than with MyChoice (51/449; 11%). Patients with HRD tumors according to shallowHRDv2 (including BRCAmut) showed a significantly prolonged PFS with bev+ola versus bev (median PFS: 65.7 versus 20.3 months, hazard ratio (HR): 0.36 [95% CI: 0.24-0.53]). This benefit was significant also for BRCA1/2 wild-type tumors (40.8 versus 19.5 months, HR: 0.45 [95% CI: 0.26-0.76]). ShallowHRDv2 is a performant, clinically validated, and cost-effective test for HRD detection.

Genomic analysis of advanced breast cancer tumors from talazoparib-treated gBRCA1/2mut carriers in the ABRAZO study

These analyses explore the impact of homologous recombination repair gene mutations, including BRCA1/2 mutations and homologous recombination deficiency (HRD), on the efficacy of the poly(ADP-ribose) polymerase (PARP) inhibitor talazoparib in the open-label, two-cohort, Phase 2 ABRAZO trial in germline BRCA1/2-mutation carriers. In the evaluable intent-to-treat population (N = 60), 58 (97%) patients harbor ≥1 BRCA1/2 mutation(s) in tumor sequencing, with 95% (53/56) concordance between germline and tumor mutations, and 85% (40/47) of evaluable patients have BRCA locus loss of heterozygosity indicating HRD. The most prevalent non-BRCA tumor mutations are TP53 in patients with BRCA1 mutations and PIK3CA in patients with BRCA2 mutations. BRCA1- or BRCA2-mutated tumors show comparable clinical benefit within cohorts. While low patient numbers preclude correlations between HRD and efficacy, germline BRCA1/2 mutation detection from tumor-only sequencing shows high sensitivity and non-BRCA genetic/genomic events do not appear to influence talazoparib sensitivity in the ABRAZO trial.ClinicalTrials.gov identifier: NCT02034916.

Management of individuals with germline pathogenic/likely pathogenic variants in CHEK2: A clinical practice resource of the American College of Medical Genetics and Genomics (ACMG)

PURPOSE

Although the role of CHEK2 germline pathogenic variants in cancer predisposition is well known, resources for managing CHEK2 heterozygotes in clinical practice are limited.

METHODS

An international workgroup developed guidance on clinical management of CHEK2 heterozygotes informed by peer-reviewed publications from PubMed.

RESULTS

Although CHEK2 is considered a moderate penetrance gene, cancer risks may be considered as a continuous variable, which are influenced by family history and other modifiers. Consequently, early cancer detection and prevention for CHEK2 heterozygotes should be guided by personalized risk estimates. Such estimates may result in both downgrading lifetime breast cancer risks to those similar to the general population or upgrading lifetime risk to a level at which CHEK2 heterozygotes are offered high-risk breast surveillance according to country-specific guidelines. Risk-reducing mastectomy should be guided by personalized risk estimates and shared decision making. Colorectal and prostate cancer surveillance should be considered based on assessment of family history. For CHEK2 heterozygotes who develop cancer, no specific targeted medical treatment is recommended at this time.

CONCLUSION

Systematic prospective data collection is needed to establish the spectrum of CHEK2-associated cancer risks and to determine yet-unanswered questions, such as the outcomes of surveillance, response to cancer treatment, and survival after cancer diagnosis.

Predicting regional somatic mutation rates using DNA motifs

How the locus-specificity of epigenetic modifications is regulated remains an unanswered question. A contributing mechanism is that epigenetic enzymes are recruited to specific loci by DNA binding factors recognizing particular sequence motifs (referred to as epi-motifs). Using these motifs to predict biological outputs depending on local epigenetic state such as somatic mutation rates would confirm their functionality. Here, we used DNA motifs including known TF motifs and epi-motifs as a surrogate of epigenetic signals to predict somatic mutation rates in 13 cancers at an average 23kbp resolution. We implemented an interpretable neural network model, called contextual regression, to successfully learn the universal relationship between mutations and DNA motifs, and uncovered motifs that are most impactful on the regional mutation rates such as TP53 and epi-motifs associated with H3K9me3. Furthermore, we identified genomic regions with significantly higher mutation rates than the expected values in each individual tumor and demonstrated that such cancer-related regions can accurately predict cancer types. Interestingly, we found that the same mutation signatures often have different contributions to cancer-related and cancer-independent regions, and we also identified the motifs with the most contribution to each mutation signature.

Patient Assessment and Therapy Planning Based on Homologous Recombination Repair Deficiency

Defects in genes involved in the DNA damage response cause homologous recombination repair deficiency (HRD). HRD is found in a subgroup of cancer patients for several tumor types, and it has a clinical relevance to cancer prevention and therapies. Accumulating evidence has identified HRD as a biomarker for assessing the therapeutic response of tumor cells to poly(ADP-ribose) polymerase inhibitors and platinum-based chemotherapies. Nevertheless, the biology of HRD is complex, and its applications and the benefits of different HRD biomarker assays are controversial. This is primarily due to inconsistencies in HRD assessments and definitions (gene-level tests, genomic scars, mutational signatures, or a combination of these methods) and difficulties in assessing the contribution of each genomic event. Therefore, we aim to review the biological rationale and clinical evidence of HRD as a biomarker. This review provides a blueprint for the standardization and harmonization of HRD assessments.

A novel marker integrating multiple genetic alterations better predicts platinum sensitivity in ovarian cancer than HRD score

Introduction: Platinum-based chemotherapy is the first-line treatment strategy for ovarian cancer patients. The dismal prognosis of ovarian cancer was shown to be stringently associated with the heterogeneity of tumor cells in response to this therapy, therefore understanding platinum sensitivity in ovarian cancer would be helpful for improving patients' quality of life and clinical outcomes. HRDetect, utilized to characterize patients' homologous recombination repair deficiency, was used to predict patients' response to platinum-based chemotherapy. However, whether each of the single features contributing to HRD score is associated with platinum sensitivity remains elusive. Methods: We analyzed the whole-exome sequencing data of 196 patients who received platinum-based chemotherapy from the TCGA database. Genetic features were determined individually to see if they could indicate patients' response to platinum-based chemotherapy and prognosis, then integrated into a Pt-score employing LASSO regression model to assess its predictive performance. Results and discussion: Multiple genetic features, including bi-allelic inactivation of BRCA1/2 genes and genes involved in HR pathway, multiple somatic mutations in genes involved in DNA damage repair (DDR), and previously reported HRD-related features, were found to be stringently associated with platinum sensitivity and improved prognosis. Higher contributions of mutational signature SBS39 or ID6 predicted improved overall survival. Besides, arm-level loss of heterozygosity (LOH) of either chr4p or chr5q predicted significantly better disease-free survival. Notably, some of these features were found independent of HRD. And SBS3, an HRD-related feature, was found irrelevant to platinum sensitivity. Integrated all candidate markers using the LASSO model to yield a Pt-score, which showed better predictive ability compared to HRDetect in determining platinum sensitivity and predicting patients' prognosis, and this performance was validated in an independent cohort. The outcomes of our study will be instrumental in devising effective strategies for treating ovarian cancer with platinum-based chemotherapy.

Homologous Recombination Deficiency Landscape of Breast Cancers and Real-World Effectiveness of Poly ADP-Ribose Polymerase Inhibitors in Patients With Somatic BRCA1/2, Germline PALB2, or Homologous Recombination Deficiency Signature

PURPOSE

Poly ADP-ribose polymerase inhibitors (PARPi) are approved for patients with human epidermal growth factor receptor 2-negative metastatic breast cancer (mBC) and germline pathogenic/likely pathogenic variant (hereafter mutation) in the BRCA1/2 genes (gBRCA); however, clinical benefit has also been demonstrated in mBC with somatic BRCA1/2 mutations (sBRCA) or germline PALB2 mutations (gPALB2). This study aims to describe the genomic landscape of homologous recombination repair (HRR) gene alterations in mBC and assess PARPi treatment outcomes for patients with gBRCA compared with other HRR genes and by status of a novel homologous recombination deficiency signature (HRDsig).

METHODS

A real-world (RW) clinico-genomic database (CGDB) of comprehensive genomic profiling (CGP) linked to deidentified, electronic health record-derived clinical data was used. CGP was analyzed for HRR genes and HRDsig. The CGDB enabled cohort characterization and outcomes analyses of 177 patients exposed to PARPi. RW progression-free survival (rwPFS) and RW overall survival (rwOS) were compared.

RESULTS

Of 28,920 patients with mBC, gBRCA was detected in 3.4%, whereas the population with any BRCA alteration or gPALB2 increased to 9.5%. HRDsig+ represented 21% of patients with mBC. BRCA and gPALB2 had higher levels of biallelic loss and HRDsig+ than other HRR alterations. Outcomes on PARPi were assessed for 177 patients, and gBRCA and sBRCA/gPALB2 cohorts were similar: gBRCA versus sBRCA/gPALB2 rwPFS was 6.3 versus 5.4 months (hazard ratio [HR], 1.37 [0.77-2.43]); rwOS was 16.2 versus 21.2 months (HR, 1.45 [0.74-2.86]). Additionally, patients with HRDsig+ versus HRDsig- had longer rwPFS (6.3 v 2.8 months; HR, 0.62 [0.42-0.92]) and numerically longer rwOS (17.8 v 13.0 months; HR, 0.72 [0.46-1.14]).

CONCLUSION

Patients with sBRCA and gPALB2 derive similar benefit from PARPi as those with gBRCA alterations. In combination, HRDsig+, sBRCA, and gPALB2 represent an additional 19% of mBC that can potentially benefit from PARPi. Randomized trials exploring a more inclusive biomarker such as HRDsig are warranted.

Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

Genetic characteristics of platinum-sensitive ovarian clear cell carcinoma

OBJECTIVE

Most ovarian clear cell carcinomas are resistant to platinum-based chemotherapy, while a small subset shows a positive response. The aim of this study was to clarify the clinical, pathological and genetic characteristics of platinum-sensitive ovarian clear cell carcinomas.

METHODS

The study included 53 patients with stage III-IV ovarian clear cell carcinoma who had residual tumours after primary surgery and received platinum-based therapy between 2009 and 2018. A retrospective examination of platinum sensitivity was performed using the criterion of ≥6 months from the last day of first-line platinum therapy until recurrence/progression. Cases determined to be platinum-sensitive were subjected to immunohistochemical staining, genomic analyses using target sequencing (i.e. NCC Oncopanel) and homologous recombination deficiency (myChoice® HRD Plus) assays.

RESULTS

Of the 53 stage III-IV ovarian clear cell carcinoma cases, 11 (21%) were platinum-sensitive. These cases showed better progression-free and overall survival than platinum-resistant cases (hazard ratio = 0.16, P < 0.001). Among the seven sensitive cases whose tumour tissues were available for molecular profiling, five were pure ovarian clear cell carcinoma based on pathological and genetic features, whereas the remaining two cases were re-diagnosed as high-grade serous ovarian carcinoma. The pure ovarian clear cell carcinomas lacked BRCA1 and BRCA2 mutations, consistent with the absence of the homologous recombination deficiency phenotype, whereas two cases (40%) had ATM mutations. By contrast, the two high-grade serous ovarian carcinoma cases had BRCA1 or BRCA2 mutations associated with the homologous recombination deficiency phenotype.

CONCLUSION

The subset of platinum-sensitive ovarian clear cell carcinomas includes a majority with pure ovarian clear cell carcinoma features that lack the homologous recombination deficiency phenotype.

Clinicogenomic characteristics and synthetic lethal implications of germline homologous recombination-deficient hepatocellular carcinoma

BACKGROUNDS AND AIMS

We performed an in-depth examination of pathogenic germline variants (PGVs) and somatic variants in DNA damage response (DDR) genes in hepatocellular carcinoma (HCC) to explore their clinical and genomic impacts.

APPROACH AND RESULTS

We used a merged whole-exome or RNA sequencing data set derived from in-house ( n = 230) and The Cancer Genome Atlas ( n = 362) databases of multiethnic HCC samples. We also evaluated synthetic lethal approaches targeting mutations in homologous recombination (HR) genes using HCC cells selected from five genomic databases of cancer cell lines. A total of 110 PGVs in DDR pathways in 96 patients were selected. Of the PGV carriers, 44 were HR-altered and found to be independently associated with poorer disease-free survival after hepatectomy. The most frequently altered HR gene in both germline and somatic tissues was POLQ , and this variant was detected in 22.7% (10/44) and 23.8% (5/21) of all the corresponding carriers, respectively. PGVs in HR were significantly associated with upregulation of proliferation and replication-related genes and familial risk of HCC. Samples harboring PGVs in HR with loss of heterozygosity were most strongly correlated with the genomic footprints of deficient HR, such as mutation burden and denovoSig2 (analogous to Catalogue of Somatic Mutations in Cancer [COSMIC] 3), and poor outcome. Pharmacologic experiments with HCC cells defective in BRCA2 or POLQ suggested that tumors with this phenotype are synthetic lethal with poly(ADP-ribose) polymerase inhibitors.

CONCLUSIONS

Our findings suggest that germline HR defects in HCC tend to confer a poor prognosis and result in distinctive genomic scarring. Tests of the clinical benefits of HR-directed treatments in the affected patients are needed.

Mutational signatures reveal mutual exclusivity of homologous recombination and mismatch repair deficiencies in colorectal and stomach tumors

Decomposing somatic mutation spectra into mutational signatures and their corresponding etiologies provides a powerful approach for investigating the mechanism of DNA damage and repair. Assessing microsatellite (in)stability (MSI/MSS) status and interpreting their clinical relevance in different malignancies offers significant diagnostic and prognostic value. However, little is known about microsatellite (in)stability and its interactions with other DNA repair mechanisms such as homologous recombination (HR) in different cancer types. Based on whole-genome/exome mutational signature analysis, we showed HR deficiency (HRd) and mismatch repair deficiency (MMRd) occur in a significantly mutually exclusive manner in stomach and colorectal adenocarcinomas. ID11 signature with currently unknown etiology was prevalent in MSS tumors, co-occurred with HRd and was mutually exclusive with MMRd. Apolipoprotein B mRNA editing enzyme, Catalytic polypeptide-like (APOBEC) signature co-occurred with HRd and was mutually exclusive with MMRd in stomach tumors. The HRd signature in MSS tumors and the MMRd signature in MSI tumors were the first or second dominant signatures wherever detected. HRd may drive a distinct subgroup of MSS tumors and lead to poor clinical outcome. These analyses offer insight into mutational signatures in MSI and MMS tumors and reveal opportunities for improved clinical diagnosis and personalized treatment of MSS tumors.

Copy number alteration features in pan-cancer homologous recombination deficiency prediction and biology

Homologous recombination deficiency (HRD) renders cancer cells vulnerable to unrepaired double-strand breaks and is an important therapeutic target as exemplified by the clinical efficacy of poly ADP-ribose polymerase (PARP) inhibitors as well as the platinum chemotherapy drugs applied to HRD patients. However, it remains a challenge to predict HRD status precisely and economically. Copy number alteration (CNA), as a pervasive trait of human cancers, can be extracted from a variety of data sources, including whole genome sequencing (WGS), SNP array, and panel sequencing, and thus can be easily applied clinically. Here we systematically evaluate the predictive performance of various CNA features and signatures in HRD prediction and build a gradient boosting machine model (HRD_CNA) for pan-cancer HRD prediction based on these CNA features. CNA features BP10MB[1] (The number of breakpoints per 10MB of DNA is 1) and SS[ > 7 & <=8] (The log10-based size of segments is greater than 7 and less than or equal to 8) are identified as the most important features in HRD prediction. HRD_CNA suggests the biallelic inactivation of BRCA1, BRCA2, PALB2, RAD51C, RAD51D, and BARD1 as the major genetic basis for human HRD, and may also be applied to effectively validate the pathogenicity of BRCA1/2 variants of uncertain significance (VUS). Together, this study provides a robust tool for cost-effective HRD prediction and also demonstrates the applicability of CNA features and signatures in cancer precision medicine.

FET fusion oncoproteins disrupt physiologic DNA repair networks in cancer

While oncogenes promote cancer cell growth, unrestrained proliferation represents a significant stressor to cellular homeostasis networks such as the DNA damage response (DDR). To enable oncogene tolerance, many cancers disable tumor suppressive DDR signaling through genetic loss of DDR pathways and downstream effectors (e.g., ATM or p53 tumor suppressor mutations). Whether and how oncogenes can help "self-tolerize" by creating analogous functional deficiencies in physiologic DDR networks is not known. Here we focus on Ewing sarcoma, a FET fusion oncoprotein (EWS-FLI1) driven pediatric bone tumor, as a model for the class of FET rearranged cancers. Native FET protein family members are among the earliest factors recruited to DNA double-strand breaks (DSBs) during the DDR, though the function of both native FET proteins and FET fusion oncoproteins in DNA repair remains to be defined. Using preclinical mechanistic studies of the DDR and clinical genomic datasets from patient tumors, we discover that the EWS-FLI1 fusion oncoprotein is recruited to DNA DSBs and interferes with native FET (EWS) protein function in activating the DNA damage sensor ATM. As a consequence of FET fusion-mediated interference with the DDR, we establish functional ATM deficiency as the principal DNA repair defect in Ewing sarcoma and the compensatory ATR signaling axis as a collateral dependency and therapeutic target in multiple FET rearranged cancers. More generally, we find that aberrant recruitment of a fusion oncoprotein to sites of DNA damage can disrupt physiologic DSB repair, revealing a mechanism for how growth-promoting oncogenes can also create a functional deficiency within tumor suppressive DDR networks.

A ribosomal gene panel predicting a novel synthetic lethality in non-BRCAness tumors

Poly (ADP-ribose) polymerase (PARP) inhibitors are one of the most exciting classes of targeted therapy agents for cancers with homologous recombination (HR) deficiency. However, many patients without apparent HR defects also respond well to PARP inhibitors/cisplatin. The biomarker responsible for this mechanism remains unclear. Here, we identified a set of ribosomal genes that predict response to PARP inhibitors/cisplatin in HR-proficient patients. PARP inhibitor/cisplatin selectively eliminates cells with high expression of the eight genes in the identified panel via DNA damage (ATM) signaling-induced pro-apoptotic ribosomal stress, which along with ATM signaling-induced pro-survival HR repair constitutes a new model to balance the cell fate in response to DNA damage. Therefore, the combined examination of the gene panel along with HR status would allow for more precise predictions of clinical response to PARP inhibitor/cisplatin. The gene panel as an independent biomarker was validated by multiple published clinical datasets, as well as by an ovarian cancer organoids library we established. More importantly, its predictive value was further verified in a cohort of PARP inhibitor-treated ovarian cancer patients with both RNA-seq and WGS data. Furthermore, we identified several marketed drugs capable of upregulating the expression of the genes in the panel without causing HR deficiency in PARP inhibitor/cisplatin-resistant cell lines. These drugs enhance PARP inhibitor/cisplatin sensitivity in both intrinsically resistant organoids and cell lines with acquired resistance. Together, our study identifies a marker gene panel for HR-proficient patients and reveals a broader application of PARP inhibitor/cisplatin in cancer therapy.

Breast cancer genomes from CHEK2 c.1100delC mutation carriers lack somatic TP53 mutations and display a unique structural variant size distribution profile

BACKGROUND

CHEK2 c.1100delC was the first moderate-risk breast cancer (BC) susceptibility allele discovered. Despite several genomic, transcriptomic and functional studies, however, it is still unclear how exactly CHEK2 c.1100delC promotes tumorigenesis. Since the mutational landscape of a tumor reflects the processes that have operated on its development, the aim of this study was to uncover the somatic genomic landscape of CHEK2-associated BC.

METHODS

We sequenced primary BC (pBC) and normal genomes of 20 CHEK2 c.1100delC mutation carriers as well as their pBC transcriptomes. Including pre-existing cohorts, we exhaustively compared CHEK2 pBC genomes to those from BRCA1/2 mutation carriers, those that displayed homologous recombination deficiency (HRD) and ER- and ER+ pBCs, totaling to 574 pBC genomes. Findings were validated in 517 metastatic BC genomes subdivided into the same subgroups. Transcriptome data from 168 ER+ pBCs were used to derive a TP53-mutant gene expression signature and perform cluster analysis with CHEK2 BC transcriptomes. Finally, clinical outcome of CHEK2 c.1100delC carriers was compared with BC patients displaying somatic TP53 mutations in two well-described retrospective cohorts totaling to 942 independent pBC cases.

RESULTS

BC genomes from CHEK2 mutation carriers were most similar to ER+ BC genomes and least similar to those of BRCA1/2 mutation carriers in terms of tumor mutational burden as well as mutational signatures. Moreover, CHEK2 BC genomes did not show any evidence of HRD. Somatic TP53 mutation frequency and the size distribution of structural variants (SVs), however, were different compared to ER+ BC. Interestingly, BC genomes with bi-allelic CHEK2 inactivation lacked somatic TP53 mutations and transcriptomic analysis indicated a shared biology with TP53 mutant BC. Moreover, CHEK2 BC genomes had an increased frequency of > 1 Mb deletions, inversions and tandem duplications with peaks at specific sizes. The high chromothripsis frequency among CHEK2 BC genomes appeared, however, not associated with this unique SV size distribution profile.

CONCLUSIONS

CHEK2 BC genomes are most similar to ER+ BC genomes, but display unique features that may further unravel CHEK2-driven tumorigenesis. Increased insight into this mechanism could explain the shorter survival of CHEK2 mutation carriers that is likely driven by intrinsic tumor aggressiveness rather than endocrine resistance.

MND1 and PSMC3IP control PARP inhibitor sensitivity in mitotic cells

The PSMC3IP-MND1 heterodimer promotes meiotic D loop formation before DNA strand exchange. In genome-scale CRISPR-Cas9 mutagenesis and interference screens in mitotic cells, depletion of PSMC3IP or MND1 causes sensitivity to poly (ADP-Ribose) polymerase inhibitors (PARPi) used in cancer treatment. PSMC3IP or MND1 depletion also causes ionizing radiation sensitivity. These effects are independent of PSMC3IP/MND1's role in mitotic alternative lengthening of telomeres. PSMC3IP- or MND1-depleted cells accumulate toxic RAD51 foci in response to DNA damage, show impaired homology-directed DNA repair, and become PARPi sensitive, even in cells lacking both BRCA1 and TP53BP1. Epistasis between PSMC3IP-MND1 and BRCA1/BRCA2 defects suggest that abrogated D loop formation is the cause of PARPi sensitivity. Wild-type PSMC3IP reverses PARPi sensitivity, whereas a PSMC3IP p.Glu201del mutant associated with D loop defects and ovarian dysgenesis does not. These observations suggest that meiotic proteins such as MND1 and PSMC3IP have a greater role in mitotic DNA repair.

Targeting homologous recombination deficiency in uterine leiomyosarcoma

BACKGROUND

Uterine leiomyosarcoma (uLMS) is a rare and aggressive gynaecological malignancy, with individuals with advanced uLMS having a five-year survival of < 10%. Mutations in the homologous recombination (HR) DNA repair pathway have been observed in ~ 10% of uLMS cases, with reports of some individuals benefiting from poly (ADP-ribose) polymerase (PARP) inhibitor (PARPi) therapy, which targets this DNA repair defect. In this report, we screened individuals with uLMS, accrued nationally, for mutations in the HR repair pathway and explored new approaches to therapeutic targeting.

METHODS

A cohort of 58 individuals with uLMS were screened for HR Deficiency (HRD) using whole genome sequencing (WGS), whole exome sequencing (WES) or NGS panel testing. Individuals identified to have HRD uLMS were offered PARPi therapy and clinical outcome details collected. Patient-derived xenografts (PDX) were generated for therapeutic targeting.

RESULTS

All 13 uLMS samples analysed by WGS had a dominant COSMIC mutational signature 3; 11 of these had high genome-wide loss of heterozygosity (LOH) (> 0.2) but only two samples had a CHORD score > 50%, one of which had a homozygous pathogenic alteration in an HR gene (deletion in BRCA2). A further three samples harboured homozygous HRD alterations (all deletions in BRCA2), detected by WES or panel sequencing, with 5/58 (9%) individuals having HRD uLMS. All five individuals gained access to PARPi therapy. Two of three individuals with mature clinical follow up achieved a complete response or durable partial response (PR) with the subsequent addition of platinum to PARPi upon minor progression during initial PR on PARPi. Corresponding PDX responses were most rapid, complete and sustained with the PARP1-specific PARPi, AZD5305, compared with either olaparib alone or olaparib plus cisplatin, even in a paired sample of a BRCA2-deleted PDX, derived following PARPi therapy in the patient, which had developed PARPi-resistance mutations in PRKDC, encoding DNA-PKcs.

CONCLUSIONS

Our work demonstrates the value of identifying HRD for therapeutic targeting by PARPi and platinum in individuals with the aggressive rare malignancy, uLMS and suggests that individuals with HRD uLMS should be included in trials of PARP1-specific PARPi.

Clinical and genomic characterization of mutational signatures across human cancers

Mutational signatures, the generic patterns of mutations, are the footprints of both endogenous and exogenous factors that have influenced cancer development. To date, dozens of mutational signatures have been discerned through computational methods. However, the etiology, mutational properties, clonality, immunology and prognostic value of mutation signatures across cancer types are poorly understood. To address this, we extensively characterized mutational signatures across 8836 cancer samples spanning 42 cancer types. We confirmed and extended clinical and genomic features associated with mutation signatures. Mutation distribution analysis showed that most mutation processes were depleted in exons and APOBEC signatures (SBS2 and SBS13), the Pol-η related signature (SBS9) and SBS40 tended to contribute clustered mutations. We observed that age-related signatures (SBS1 and SBS5) and SBS40 tended to induce mutations affecting cancer genes and subclonal drivers posted by specific signatures (eg, mismatch repair deficiency-related signature SBS44) were unlikely subjected to positive selection. We also revealed early mutation signatures (eg, UV light exposure-related signature SBS7a) and signatures (eg, reactive oxygen species-related signature SBS18) predominated in the late stage of tumorigenesis. Comprehensive association analysis of mutation processes with microenvironment revealed that APOBEC- and mismatch repair deficiency-related signatures were positively associated with immune parameters, while age-related signatures showed negative correlations. In addition, prognostic association analysis showed that many signatures were favorable (eg, SBS9) or adverse factors (eg, SBS18) of patient survival. Our findings enhance appreciation of the role of mutational signatures in tumor evolution and underline their potential in immunotherapy guidance and prognostic prediction.

Pathogenic Variants in Adult-Onset Cancer Predisposition Genes in Pediatric Cancer: Prevalence and Impact on Tumor Molecular Features and Clinical Management

PURPOSE

Clinical genomic sequencing of pediatric tumors is increasingly uncovering pathogenic variants in adult-onset cancer predisposition genes (aoCPG). Nevertheless, it remains poorly understood how often aoCPG variants are of germline origin and whether they influence tumor molecular profiles and/or clinical care. In this study, we examined the prevalence, spectrum, and impacts of aoCPG variants on tumor genomic features and patient management at our institution.

EXPERIMENTAL DESIGN

This is a retrospective study of 1,018 children with cancer who underwent clinical genomic sequencing of their tumors. Tumor genomic data were queried for pathogenic variants affecting 24 preselected aoCPGs. Available tumor whole-genome sequencing (WGS) data were evaluated for second hit mutations, loss of heterozygosity (LOH), DNA mutational signatures, and homologous recombination deficiency (HRD). Patients whose tumors harbored one or more pathogenic aoCPG variants underwent subsequent germline testing based on hereditary cancer evaluation and family or provider preference.

RESULTS

Thirty-three patients (3%) had tumors harboring pathogenic variants affecting one or more aoCPGs. Among 21 tumors with sufficient WGS sequencing data, six (29%) harbored a second hit or LOH affecting the remaining aoCPG allele with four of these six tumors (67%) also exhibiting a DNA mutational signature consistent with the altered aoCPG. Two additional tumors demonstrated HRD, of uncertain relation to the identified aoCPG variant. Twenty-one of 26 patients (81%) completing germline testing were positive for the aoCPG variant in the germline. All germline-positive patients were counseled regarding future cancer risks, surveillance, and risk-reducing measures. No patients had immediate cancer therapy changed due to aoCPG data.

CONCLUSIONS

AoCPG variants are rare in pediatric tumors; however, many originate in the germline. Almost one third of tumor aoCPG variants examined exhibited a second hit and/or conferred an abnormal DNA mutational profile suggesting a role in tumor formation. aoCPG information aids in cancer risk prediction but is not commonly used to alter the treatment of pediatric cancers.

Distinct breast cancer phenotypes in BRCA 1/2 carriers based on ER status

PURPOSE

BRCA1/2 genes are the two main genes associated with hereditary breast cancers (BC). In the present study, we explore clinical and molecular characteristics of BRCA-associated BC in relation to estrogen receptor (ER) status.

METHODS

Three BC databases (DB) were evaluated: (i) Hadassah oncogenetics (n = 4826); (ii) METABRIC (n = 1980), and (iii) Nick-Zainal (n = 560). We evaluated age at diagnosis in BRCA positive (BP) and BRCA negative (BN) patients, and tested for mutational signature differences in cohort iii. mRNA differential expression analysis (DEA) and pathway analysis were performed in cohort ii.

RESULTS

Age at diagnosis was lower in BP vs. BN tumors in all cohorts in the ER- group, and only in cohort i for the ER + group. Signature 3 was universal in BP BC, whereas several signatures were associated with ER status. Pathway analysis was performed between BP&BN, and was significant only in ER- tumors: the major activated pathways involved cancer-related processes and were highly significant. The most significant pathway was estrogen-mediated S-phase entry and the most activated upstream regulator was ERBB2.

CONCLUSION

Signature 3 was universal for all BP BC, while other signatures were associated with ER status. ER + BP& BN show similar genomic characteristics, ER- BP differed markedly from BN. This suggests that the initial carcinogenic process is universal for all BRCA carriers, but further insults lead to the development of two genomically distinct subtypes ER- and ER + . This may shed light on possible mechanisms involved in BP and carry preventive and therapeutic implications.

Prognostic relevance of HRDness gene expression signature in ovarian high-grade serous carcinoma; JGOG3025-TR2 study

BACKGROUND

This study aimed to evaluate the homologous recombination repair pathway deficiency (HRD) in ovarian high-grade serous carcinoma (HGSC).

METHODS

In the ovarian cancer data from The Cancer Genome Atlas, we identified genes differentially expressed between tumours with and without HRD genomic scars and named these genes "HRDness signature". We performed SNP array, RNA sequencing, and methylation array analyses on 274 HGSC tumours for which targeted sequencing of 51 genes and clinical data were available to generate JGOG3025-TR2 dataset. The HRDness signature was tested on external datasets, including the JGOG3025-TR2 cohort, by computational scoring and machine-learning prediction.

RESULTS

High scores and positive predictions of the HRDness signature were significantly associated with BRCA alterations, genomic scar scores, and better survival. On the other hand, among cases with high scores and/or positive predictions, those with BRCA1 methylation showed poorer survival. In the JGOG3025-TR2 cohort, HRD status was significantly associated with the use of olaparib after relapse and progression-free survival after its initiation.

CONCLUSIONS

The HRDness gene expression signature is associated with a good prognosis, while BRCA1 methylation is associated with a poor prognosis. The newly generated JGOG3025-TR2 dataset will be useful in future HGSC studies.

Homologous recombination deficiency prediction using low-pass whole genome sequencing in breast cancer

Homologous recombination repair deficiency (HRD) results in a defect in DNA repair and is a frequent driver of tumorigenesis. Poly(ADP-ribose) polymerase inhibitors (PARPi) or platinum-based therapies have increased theraputic effectiveness when treating HRD positive cancers. For breast cancer and ovairan cancer HRD companion diagnostic tests are commonly used. However, the currently used HRD tests are based on high-depth genome sequencing or hybridization-based capture sequencing, which are technically complex and costly. In this study, we modified an existing method named shallowHRD, which uses low-pass whole genome sequencing (WGS) for HRD detection, and estimated the performance of the modified shallowHRD pipeline. Our shallowHRD pipeline achieved an AUC of 0.997 in simulated low-pass WGS data, with a sensitivity of 0.981 and a specificity of 0.964; and achieved a higher HRD risk score in clinical BRCA-deficient breast cancer samples (p = 5.5 × 10^-5, compared with BRCA-intact breast cancer samples). We also estimated the limit of detection the shallowHRD pipeline could accurately predict HRD status with a minimum sequencing depth of 0.1 ×, a tumor purity of > 20%, and an input DNA amount of 1 ng. Our study demostrates using low-pass sequencing, HRD status can be determined with high accuracy using a simple approach with greatly reduced cost.

Analysis of evolutionary dynamics and clonal architecture in prostate cancer

The extent to which clinical and genomic characteristics associate with prostate cancer clonal architecture, tumor evolution, and therapeutic response remains unclear. Here, we reconstructed the clonal architecture and evolutionary trajectories of 845 prostate cancer tumors with harmonized clinical and molecular data. We observed that tumors from patients who self-reported as Black had more linear and monoclonal architectures, despite these men having higher rates of biochemical recurrence. This finding contrasts with prior observations relating polyclonal architecture to adverse clinical outcomes. Additionally, we utilized a novel approach to mutational signature analysis that leverages clonal architecture to uncover additional cases of homologous recombination and mismatch repair deficiency in primary and metastatic tumors and link the origin of mutational signatures to specific subclones. Broadly, prostate cancer clonal architecture analysis reveals novel biological insights that may be immediately clinically actionable and provide multiple opportunities for subsequent investigation.

Statement of significance

Tumors from patients who self-reported as Black demonstrate linear and monoclonal evolutionary trajectories yet experience higher rates of biochemical recurrence. In addition, analysis of clonal and subclonal mutational signatures identifies additional tumors with potentially actionable alterations such as deficiencies in mismatch repair and homologous recombination.

Homologous recombination deficiency signatures in gastrointestinal and thoracic cancers correlate with platinum therapy duration

There is emerging evidence about the predictive role of homologous recombination deficiency (HRD), but this is less defined in gastrointestinal (GI) and thoracic malignancies. We reviewed whole genome (WGS) and transcriptomic (RNA-Seq) data from advanced GI and thoracic cancers in the Personalized OncoGenomics trial (NCT02155621) to evaluate HRD scores and single base substitution (SBS)3, which is associated with BRCA1/2 mutations and potentially predictive of defective HRD. HRD scores were calculated by sum of loss of heterozygosity, telomeric allelic imbalance, and large-scale state transitions scores. Regression analyses examined the association between HRD and time to progression on platinum (TTPp). We included 223 patients with GI (n = 154) or thoracic (n = 69) malignancies. TTPp was associated with SBS3 (p < 0.01) but not HRD score in patients with GI malignancies, whereas neither was associated with TTPp in thoracic malignancies. Tumors with gBRCA1/2 mutations and a somatic second alteration exhibited high SBS3 and HRD scores, but these signatures were also present in several tumors with germline but no somatic second alterations, suggesting silencing of the wild-type allele or BRCA1/2 haploinsufficiency. Biallelic inactivation of an HR gene, including loss of XRCC2 and BARD1, was identified in BRCA1/2 wild-type HRD tumors and these patients had prolonged response to platinum. Thoracic cases with high HRD score were associated with high RECQL5 expression (p ≤ 0.025), indicating another potential mechanism of HRD. SBS3 was more strongly associated with TTPp in patients with GI malignancies and may be complementary to using HRD and BRCA status in identifying patients who benefit from platinum therapy.

Integrated multi-omics approach to distinct molecular characterization and classification of early-onset colorectal cancer

Incidence of early-onset colorectal cancer (EOCRC), defined by a diagnosed age under 50 years, is increasing, but its heterogeneous etiologies that differ from general CRC remain undetermined. We initially characterize the genome, epigenome, transcriptome, and proteome of tumors from 79 patients in a Chinese CRC cohort. Data for an additional 126 EOCRC subjects are obtained from the International Cancer Genome Consortium Chinese cohort and The Cancer Genome Atlas European cohort. We observe that early-onset tumors have a high tumor mutation burden; increased DNA repair features by mutational signature 3 and multi-layer pathway enrichments; strong perturbations at effects of DNA methylation and somatic copy-number alteration on gene expression; and upregulated immune infiltration as hot tumors underlying immunophenotypes. Notably, LMTK3 exhibits ancestral mutation disparity, potentially being a functional modulator and biomarker that drives molecular alterations in EOCRC development and immunotherapies. This integrative omics study provides valuable knowledge for precision oncology of CRC.