Increased frequency of CHD1 deletions in prostate cancers of African American men is associated with rapid disease progression without inducing homologous recombination deficiency.

We analyzed genomic data derived from the prostate cancer of African and European American men in order to identify differences that may contribute to racial disparity of outcome and that could also define novel therapeutic strategies. In addition to analyzing patient derived next generation sequencing data, we performed FISH based confirmatory studies of Chromodomain helicase DNA-binding protein 1 (CHD1) loss on prostate cancer tissue microarrays. We created CRISPR edited, CHD1 deficient prostate cancer cell lines for genomic, drug sensitivity and functional homologous recombination (HR) activity analysis. We found that subclonal deletion of CHD1 is nearly three times as frequent in prostate tumors of African American men than in men of European ancestry and it associates with rapid disease progression. We further showed that CHD1 deletion is not associated with homologous recombination deficiency associated mutational signatures in prostate cancer. In prostate cancer cell line models CHD1 deletion did not induce HR deficiency as detected by RAD51 foci formation assay or mutational signatures, which was consistent with the moderate increase of olaparib sensitivity. CHD1 deficient prostate cancer cells, however, showed higher sensitivity to talazoparib. CHD1 loss may contribute to worse outcome of prostate cancer in African American men. A deeper understanding of the interaction between CHD1 loss and PARP inhibitor sensitivity will be needed to determine the optimal use of targeted agents such as talazoparib in the context of castration resistant prostate cancer.

Mutational signature-based identification of DNA repair deficient gastroesophageal adenocarcinomas for therapeutic targeting.

Homologous recombination (HR) and nucleotide excision repair (NER) are the two most frequently disabled DNA repair pathways in cancer. HR-deficient breast, ovarian, pancreatic and prostate cancers respond well to platinum chemotherapy and PARP inhibitors. However, the frequency of HR deficiency in gastric and esophageal adenocarcinoma (GEA) still lacks diagnostic and functional validation. Using whole exome and genome sequencing data, we found that a significant subset of GEA, but very few colorectal adenocarcinomas, show evidence of HR deficiency by mutational signature analysis (HRD score). High HRD gastric cancer cell lines demonstrated functional HR deficiency by RAD51 foci assay and increased sensitivity to platinum chemotherapy and PARP inhibitors. Of clinical relevance, analysis of three different GEA patient cohorts demonstrated that platinum treated HR deficient cancers had better outcomes. A gastric cancer cell line with strong sensitivity to cisplatin showed HR proficiency but exhibited NER deficiency by two photoproduct repair assays. Single-cell RNA-sequencing revealed that, in addition to inducing apoptosis, cisplatin treatment triggered ferroptosis in a NER-deficient gastric cancer, validated by intracellular GSH assay. Overall, our study provides preclinical evidence that a subset of GEAs harbor genomic features of HR and NER deficiency and may therefore benefit from platinum chemotherapy and PARP inhibitors.

Biologically informed deep learning for explainable epigenetic clocks.

Ageing is often characterised by progressive accumulation of damage, and it is one of the most important risk factors for chronic disease development. Epigenetic mechanisms including DNA methylation could functionally contribute to organismal aging, however the key functions and biological processes may govern ageing are still not understood. Although age predictors called epigenetic clocks can accurately estimate the biological age of an individual based on cellular DNA methylation, their models have limited ability to explain the prediction algorithm behind and underlying key biological processes controlling ageing. Here we present XAI-AGE, a biologically informed, explainable deep neural network model for accurate biological age prediction across multiple tissue types. We show that XAI-AGE outperforms the first-generation age predictors and achieves similar results to deep learning-based models, while opening up the possibility to infer biologically meaningful insights of the activity of pathways and other abstract biological processes directly from the model.

ATM deficiency confers specific therapeutic vulnerabilities in bladder cancer.

Ataxia-telangiectasia mutated (ATM) plays a central role in the cellular response to DNA damage and ATM alterations are common in several tumor types including bladder cancer. However, the specific impact of ATM alterations on therapy response in bladder cancer is uncertain. Here, we combine preclinical modeling and clinical analyses to comprehensively define the impact of ATM alterations on bladder cancer. We show that ATM loss is sufficient to increase sensitivity to DNA-damaging agents including cisplatin and radiation. Furthermore, ATM loss drives sensitivity to DNA repair-targeted agents including poly(ADP-ribose) polymerase (PARP) and Ataxia telangiectasia and Rad3 related (ATR) inhibitors. ATM loss alters the immune microenvironment and improves anti-PD1 response in preclinical bladder models but is not associated with improved anti-PD1/PD-L1 response in clinical cohorts. Last, we show that ATM expression by immunohistochemistry is strongly correlated with response to chemoradiotherapy. Together, these data define a potential role for ATM as a predictive biomarker in bladder cancer.

Nucleotide excision repair deficiency is a targetable therapeutic vulnerability in clear cell renal cell carcinoma.

Due to a demonstrated lack of DNA repair deficiencies, clear cell renal cell carcinoma (ccRCC) has not benefitted from targeted synthetic lethality-based therapies. We investigated whether nucleotide excision repair (NER) deficiency is present in an identifiable subset of ccRCC cases that would render those tumors sensitive to therapy targeting this specific DNA repair pathway aberration. We used functional assays that detect UV-induced 6-4 pyrimidine-pyrimidone photoproducts to quantify NER deficiency in ccRCC cell lines. We also measured sensitivity to irofulven, an experimental cancer therapeutic agent that specifically targets cells with inactivated transcription-coupled nucleotide excision repair (TC-NER). In order to detect NER deficiency in clinical biopsies, we assessed whole exome sequencing data for the presence of an NER deficiency associated mutational signature previously identified in ERCC2 mutant bladder cancer. Functional assays showed NER deficiency in ccRCC cells. Some cell lines showed irofulven sensitivity at a concentration that is well tolerated by patients. Prostaglandin reductase 1 (PTGR1), which activates irofulven, was also associated with this sensitivity. Next generation sequencing data of the cell lines showed NER deficiency-associated mutational signatures. A significant subset of ccRCC patients had the same signature and high PTGR1 expression. ccRCC cell line-based analysis showed that NER deficiency is likely present in this cancer type. Approximately 10% of ccRCC patients in the TCGA cohort showed mutational signatures consistent with ERCC2 inactivation associated NER deficiency and also substantial levels of PTGR1 expression. These patients may be responsive to irofulven, a previously abandoned anticancer agent that has minimal activity in NER-proficient cells.

Nucleotide excision repair deficiency is a targetable therapeutic vulnerability in clear cell renal cell carcinoma.

Purpose: Due to a demonstrated lack of DNA repair deficiencies, clear cell renal cell carcinoma (ccRCC) has not benefitted from targeted synthetic lethality-based therapies. We investigated whether nucleotide excision repair (NER) deficiency is present in an identifiable subset of ccRCC cases that would render those tumors sensitive to therapy targeting this specific DNA repair pathway aberration. Experimental Design: We used functional assays that detect UV-induced 6-4 pyrimidine-pyrimidone photoproducts to quantify NER deficiency in ccRCC cell lines. We also measured sensitivity to irofulven, an experimental cancer therapeutic agent that specifically targets cells with inactivated transcription-coupled nucleotide excision repair (TC-NER). In order to detect NER deficiency in clinical biopsies, we assessed whole exome sequencing data for the presence of an NER deficiency associated mutational signature previously identified in ERCC2 mutant bladder cancer. Results: Functional assays showed NER deficiency in ccRCC cells. Irofulven sensitivity increased in some cell lines. Prostaglandin reductase 1 (PTGR1), which activates irofulven, was also associated with this sensitivity. Next generation sequencing data of the cell lines showed NER deficiency-associated mutational signatures. A significant subset of ccRCC patients had the same signature and high PTGR1 expression. Conclusions: ccRCC cell line based analysis showed that NER deficiency is likely present in this cancer type. Approximately 10% of ccRCC patients in the TCGA cohort showed mutational signatures consistent with ERCC2 inactivation associated NER deficiency and also substantial levels of PTGR1 expression. These patients may be responsive to irofulven, a previously abandoned anticancer agent that has minimal activity in NER-proficient cells.

Hierarchy and control of ageing-related methylation networks.

DNA methylation provides one of the most widely studied biomarkers of ageing. Since the methylation of CpG dinucleotides function as switches in cellular mechanisms, it is plausible to assume that by proper adjustment of these switches age may be tuned. Though, adjusting hundreds of CpG methylation levels coherently may never be feasible and changing just a few positions may lead to biologically unstable state. A prominent example of methylation-based age estimators is provided by Horvath's clock, based on 353 CpG dinucleotides, showing a high correlation (not necessarily causation) with chronological age across multiple tissue types. On this small subset of CpG dinucleotides we demonstrate how the adjustment of one methylation level leads to a cascade of changes at other sites. Among the studied subset, we locate the most important CpGs (and related genes) that may have a large influence on the rest of the sub-system. According to our analysis, the structure of this network is way more hierarchical compared to what one would expect based on ensembles of uncorrelated connections. Therefore, only a handful of CpGs is enough to modify the system towards a desired state. When propagation of the change over the network is taken into account, the resulting modification in the predicted age can be significantly larger compared to the effect of isolated CpG perturbations. By adjusting the most influential single CpG site and following the propagation of methylation level changes we can reach up to 5.74 years in virtual age reduction, significantly larger than without taking into account of the network control. Extending our approach to the whole methylation network may identify key nodes that have controller role in the ageing process.

RAF1 amplification drives a subset of bladder tumors and confers sensitivity to MAPK-directed therapeutics.

Bladder cancer is a genetically heterogeneous disease, and novel therapeutic strategies are needed to expand treatment options and improve clinical outcomes. Here, we identified a unique subset of urothelial tumors with focal amplification of the RAF1 (CRAF) kinase gene. RAF1-amplified tumors had activation of the RAF/MEK/ERK signaling pathway and exhibited a luminal gene expression pattern. Genetic studies demonstrated that RAF1-amplified tumors were dependent upon RAF1 activity for survival, and RAF1-activated cell lines and patient-derived models were sensitive to available and emerging RAF inhibitors as well as combined RAF plus MEK inhibition. Furthermore, we found that bladder tumors with HRAS- or NRAS-activating mutations were dependent on RAF1-mediated signaling and were sensitive to RAF1-targeted therapy. Together, these data identified RAF1 activation as a dependency in a subset making up nearly 20% of urothelial tumors and suggested that targeting RAF1-mediated signaling represents a rational therapeutic strategy.

Comparative Assessment of Diagnostic Homologous Recombination Deficiency-Associated Mutational Signatures in Ovarian Cancer.

PURPOSE: Homologous recombination (HR) deficiency (HRD) is one of the key determinants of PARP inhibitor response in ovarian cancer, and its accurate detection in tumor biopsies is expected to improve the efficacy of this therapy. Because HRD induces a wide array of genomic aberrations, mutational signatures may serve as a companion diagnostic to identify PARP inhibitor-responsive cases. EXPERIMENTAL DESIGN: From the The Cancer Genome Atlas (TCGA) whole-exome sequencing (WES) data, we extracted different types of mutational signature-based HRD measures, such as the HRD score, genome-wide LOH, and HRDetect trained on ovarian and breast cancer-specific sequencing data. We compared their performance to identify BRCA1/2-deficient cases in the TCGA ovarian cancer cohort and predict survival benefit in platinum-treated, BRCA1/2 wild-type ovarian cancer. RESULTS: We found that the HRD score, which is based on large chromosomal alterations alone, performed similarly well to an ovarian cancer-specific HRDetect, which incorporates mutations on a finer scale as well (AUC = 0.823 vs. AUC = 0.837). In an independent cohort these two methods were equally accurate predicting long-term survival after platinum treatment (AUC = 0.787 vs. AUC = 0.823). We also found that HRDetect trained on ovarian cancer was more accurate than HRDetect trained on breast cancer data (AUC = 0.837 vs. AUC = 0.795; P = 0.0072). CONCLUSIONS: When WES data are available, methods that quantify only large chromosomal alterations such as the HRD score and HRDetect that captures a wider array of HRD-induced genomic aberrations are equally efficient identifying HRD ovarian cancer cases.

Genetic Determinants of Somatic Selection of Mutational Processes in 3,566 Human Cancers.

The somatic landscape of the cancer genome results from different mutational processes represented by distinct "mutational signatures." Although several mutagenic mechanisms are known to cause specific mutational signatures in cell lines, the variation of somatic mutational activities in patients, which is mostly attributed to somatic selection, is still poorly explained. Here, we introduce a quantitative trait, mutational propensity (MP), and describe an integrated method to infer genetic determinants of variations in the mutational processes in 3,566 cancers with specific underlying mechanisms. As a result, we report 2,314 candidate determinants with both significant germline and somatic effects on somatic selection of mutational processes, of which, 485 act via cancer gene expression and 1,427 act through the tumor-immune microenvironment. These data demonstrate that the genetic determinants of MPs provide complementary information to known cancer driver genes, clonal evolution, and clinical biomarkers. SIGNIFICANCE: The genetic determinants of the somatic mutational processes in cancer elucidate the biology underlying somatic selection and evolution of cancers and demonstrate complementary predictive power across cancer types.

A subset of lung cancer cases shows robust signs of homologous recombination deficiency associated genomic mutational signatures.

PARP inhibitors are approved for the treatment of solid tumor types that frequently harbor alterations in the key homologous recombination (HR) genes, BRCA1/2. Other tumor types, such as lung cancer, may also be HR deficient, but the frequency of such cases is less well characterized. Specific DNA aberration profiles (mutational signatures) are induced by homologous recombination deficiency (HRD) and their presence can be used to assess the presence or absence of HR deficiency in a given tumor biopsy even in the absence of an observed alteration of an HR gene. We derived various HRD-associated mutational signatures from whole-genome and whole-exome sequencing data in the lung adenocarcinoma and lung squamous carcinoma cases from TCGA, and in a patient of ours with stage IVA lung cancer with exceptionally good response to platinum-based therapy, and in lung cancer cell lines. We found that a subset of the investigated cases, both with and without biallelic loss of BRCA1 or BRCA2, showed robust signs of HR deficiency. The extreme platinum responder case also showed a robust HRD-associated genomic mutational profile. HRD-associated mutational signatures were also associated with PARP inhibitor sensitivity in lung cancer cell lines. Consequently, lung cancer cases with HRD, as identified by diagnostic mutational signatures, may benefit from PARP inhibitor therapy.

Detection of Molecular Signatures of Homologous Recombination Deficiency in Bladder Cancer.

PURPOSE: Poly (ADP ribose)-polymerase (PARP) inhibitors are approved for use in breast, ovarian, prostate, and pancreatic cancers, which are the solid tumor types that most frequently have alterations in key homologous recombination (HR) genes, such as BRCA1/2. However, the frequency of HR deficiency (HRD) in other solid tumor types, including bladder cancer, is less well characterized. EXPERIMENTAL DESIGN: Specific DNA aberration profiles (mutational signatures) are induced by HRD, and the presence of these "genomic scars" can be used to assess the presence or absence of HRD in a given tumor biopsy even in the absence of an observed alteration of an HR gene. Using whole-exome and whole-genome data, we measured various HRD-associated mutational signatures in bladder cancer. RESULTS: We found that a subset of bladder tumors have evidence of HRD. In addition to a small number of tumors with biallelic BRCA1/2 events, approximately 10% of bladder tumors had significant evidence of HRD-associated mutational signatures. Increased levels of HRD signatures were associated with promoter methylation of RBBP8, which encodes CtIP, a key protein involved in HR. CONCLUSIONS: A subset of bladder tumors have genomic features suggestive of HRD and therefore may be more likely to benefit from therapies such as platinum agents and PARP inhibitors that target tumor HRD.

Strand Orientation Bias Detector to determine the probability of FFPE sequencing artifacts.

Formalin-fixed paraffin-embedded tissue, the most common tissue specimen stored in clinical practice, presents challenges in the analysis due to formalin-induced artifacts. Here, we present Strand Orientation Bias Detector (SOBDetector), a flexible computational platform compatible with all the common somatic SNV-calling pipelines, designed to assess the probability whether a given detected mutation is an artifact. The underlying predictor mechanism is based on the posterior distribution of a Bayesian logistic regression model trained on The Cancer Genome Atlas whole exomes. SOBDetector is a freely available cross-platform program, implemented in Java 1.8.

Identification of a Synthetic Lethal Relationship between Nucleotide Excision Repair Deficiency and Irofulven Sensitivity in Urothelial Cancer.

PURPOSE: Cisplatin-based chemotherapy is a first-line treatment for muscle-invasive and metastatic urothelial cancer. Approximately 10% of bladder urothelial tumors have a somatic missense mutation in the nucleotide excision repair (NER) gene, ERCC2, which confers increased sensitivity to cisplatin-based chemotherapy. However, a significant subset of patients is ineligible to receive cisplatin-based therapy due to medical contraindications, and no NER-targeted approaches are available for platinum-ineligible or platinum-refractory ERCC2-mutant cases. EXPERIMENTAL DESIGN: We used a series of NER-proficient and NER-deficient preclinical tumor models to test sensitivity to irofulven, an abandoned anticancer agent. In addition, we used available clinical and sequencing data from multiple urothelial tumor cohorts to develop and validate a composite mutational signature of ERCC2 deficiency and cisplatin sensitivity. RESULTS: We identified a novel synthetic lethal relationship between tumor NER deficiency and sensitivity to irofulven. Irofulven specifically targets cells with inactivation of the transcription-coupled NER (TC-NER) pathway and leads to robust responses in vitro and in vivo, including in models with acquired cisplatin resistance, while having minimal effect on cells with intact NER. We also found that a composite mutational signature of ERCC2 deficiency was strongly associated with cisplatin response in patients and was also associated with cisplatin and irofulven sensitivity in preclinical models. CONCLUSIONS: Tumor NER deficiency confers sensitivity to irofulven, a previously abandoned anticancer agent, with minimal activity in NER-proficient cells. A composite mutational signature of NER deficiency may be useful in identifying patients likely to respond to NER-targeting agents, including cisplatin and irofulven.See related commentary by Jiang and Greenberg, p. 1833.

Detection of Molecular Signatures of Homologous Recombination Deficiency in Prostate Cancer with or without BRCA1/2 Mutations.

PURPOSE: Prostate cancers with mutations in genes involved in homologous recombination (HR), most commonly BRCA2, respond favorably to PARP inhibition and platinum-based chemotherapy. We investigated whether other prostate tumors that do not harbor deleterious mutations in these particular genes can similarly be deficient in HR, likely rendering those sensitive to HR-directed therapies. EXPERIMENTAL DESIGN: Homologous recombination deficiency (HRD) levels can be estimated using various mutational signatures derived from next-generation sequencing data. We used this approach on whole-genome sequencing (WGS; n = 311) and whole-exome sequencing (WES) data (n = 498) of both primary and metastatic prostate adenocarcinomas to determine whether prostate cancer cases display clear signs of HRD in somatic tumor biopsies. RESULTS: Known BRCA-deficient samples showed all previously described HRD-associated mutational signatures in the WGS data. HRD-associated mutational signatures were also detected in a subset of patients who did not harbor germline or somatic mutations in BRCA1/2 or other HR-related genes. Similar results, albeit with lower sensitivity and accuracy, were also obtained from WES data. CONCLUSIONS: These findings may expand the number of cases likely to respond to PARP inhibitor treatment. On the basis of the HR-associated mutational signatures, 5% to 8% of localized prostate cancer cases may be good candidates for PARP-inhibitor treatment (including those with BRCA1/2 mutations).