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.

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.

Identifying regulators of aberrant stem cell and differentiation activity in colorectal cancer using a dual endogenous reporter system.

Aberrant stem cell-like activity and impaired differentiation are central to the development of colorectal cancer (CRC). To identify functional mediators of these key cellular programs, we engineer a dual endogenous reporter system by genome-editing the SOX9 and KRT20 loci of human CRC cell lines to express fluorescent reporters, broadcasting aberrant stem cell-like and differentiation activity, respectively. By applying a CRISPR screen targeting 78 epigenetic regulators with 542 sgRNAs to this platform, we identify factors that contribute to stem cell-like activity and differentiation in CRC. Perturbation single cell RNA sequencing (Perturb-seq) of validated hits nominate SMARCB1 of the BAF complex (also known as SWI/SNF) as a negative regulator of differentiation across an array of neoplastic colon models. SMARCB1 is a dependency and required for in vivo growth of human CRC models. These studies highlight the utility of biologically designed endogenous reporter platforms to uncover regulators with therapeutic potential.

Utilizing a dual endogenous reporter system to identify functional regulators of aberrant stem cell and differentiation activity in colorectal cancer.

Aberrant stem cell-like activity and impaired differentiation are central to the development of colorectal cancer (CRC). To identify functional mediators that regulate these key cellular programs in CRC, we developed an endogenous reporter system by genome-editing human CRC cell lines with knock-in fluorescent reporters at the SOX9 and KRT20 locus to report aberrant stem cell-like activity and differentiation, respectively, and then performed pooled genetic perturbation screens. Constructing a dual reporter system that simultaneously monitored aberrant stem cell-like and differentiation activity in the same CRC cell line improved our signal to noise discrimination. Using a focused-library CRISPR screen targeting 78 epigenetic regulators with 542 sgRNAs, we identified factors that contribute to stem cell-like activity and differentiation in CRC. Perturbation single cell RNA sequencing (Perturb-seq) of validated hits nominated SMARCB1 of the BAF complex (also known as SWI/SNF) as a negative regulator of differentiation across an array of neoplastic colon models. SMARCB1 is a dependency in CRC and required for in vivo growth of human CRC models. These studies highlight the utility of a biologically designed endogenous reporter system to uncover novel therapeutic targets for drug development.

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.

Replicative stress in gastroesophageal cancer is associated with chromosomal instability and sensitivity to DNA damage response inhibitors.

Gastroesophageal adenocarcinoma (GEA) is an aggressive malignancy with chromosomal instability (CIN). To understand adaptive responses enabling DNA damage response (DDR) and CIN, we analyzed matched normal, premalignant, and malignant gastric lesions from human specimens and a carcinogen-induced mouse model, observing activation of replication stress, DDR, and p21 in neoplastic progression. In GEA cell lines, expression of DDR markers correlated with ploidy abnormalities, such as number of high-level focal amplifications and whole-genome duplication (WGD). Integrating TP53 status, ploidy abnormalities, and DDR markers into a compositive score helped predict GEA cell lines with enhanced sensitivity to Chk1/2 and Wee1 inhibition, either alone or combined with irinotecan (SN38). We demonstrate that Chk1/2 or Wee1 inhibition combined with SN38/irinotecan shows greater anti-tumor activity in human gastric cancer organoids and an in vivo xenograft mouse model. These findings indicate that specific DDR biomarkers and ploidy abnormalities may predict premalignant progression and response to DDR pathway inhibitors.

A biallelic multiple nucleotide length polymorphism explains functional causality at 5p15.33 prostate cancer risk locus.

To date, single-nucleotide polymorphisms (SNPs) have been the most intensively investigated class of polymorphisms in genome wide associations studies (GWAS), however, other classes such as insertion-deletion or multiple nucleotide length polymorphism (MNLPs) may also confer disease risk. Multiple reports have shown that the 5p15.33 prostate cancer risk region is a particularly strong expression quantitative trait locus (eQTL) for Iroquois Homeobox 4 (IRX4) transcripts. Here, we demonstrate using epigenome and genome editing that a biallelic (21 and 47 base pairs (bp)) MNLP is the causal variant regulating IRX4 transcript levels. In LNCaP prostate cancer cells (homozygous for the 21 bp short allele), a single copy knock-in of the 47 bp long allele potently alters the chromatin state, enabling de novo functional binding of the androgen receptor (AR) associated with increased chromatin accessibility, Histone 3 lysine 27 acetylation (H3K27ac), and ~3-fold upregulation of IRX4 expression. We further show that an MNLP is amongst the strongest candidate susceptibility variants at two additional prostate cancer risk loci. We estimated that at least 5% of prostate cancer risk loci could be explained by functional non-SNP causal variants, which may have broader implications for other cancers GWAS. More generally, our results underscore the importance of investigating other classes of inherited variation as causal mediators of human traits.

Replicative stress in gastroesophageal adenocarcinoma is associated with chromosomal instability and sensitivity to DNA damage response inhibitors.

Gastroesophageal adenocarcinoma (GEA) is an aggressive, often lethal, malignancy that displays marked chromosomal instability (CIN). To understand adaptive responses that enable CIN, we analyzed paired normal, premalignant, and malignant gastric lesions from human specimens and a carcinogen-induced mouse model, observing activation of replication stress, DNA damage response (DDR), and cell cycle regulator p21 in neoplastic progression. In GEA cell lines, expression of DDR markers correlated with ploidy abnormalities, including high-level focal amplifications and whole-genome duplication (WGD). Moreover, high expression of DNA damage marker H2AX correlated with CIN, WGD, and inferior patient survival. By developing and implementing a composite diagnostic score that incorporates TP53 mutation status, ploidy abnormalities, and H2AX expression, among other genomic information, we can identify GEA cell lines with enhanced sensitivity to DDR pathway inhibitors targeting Chk1/2 and Wee1. Anti-tumor properties were further augmented in combination with irinotecan (SN38) but not gemcitabine chemotherapy. These results implicate specific DDR biomarkers and ploidy abnormalities as diagnostic proxy that may predict premalignant progression and response to DDR pathway inhibitors.

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.

Epigenomic charting and functional annotation of risk loci in renal cell carcinoma.

While the mutational and transcriptional landscapes of renal cell carcinoma (RCC) are well-known, the epigenome is poorly understood. We characterize the epigenome of clear cell (ccRCC), papillary (pRCC), and chromophobe RCC (chRCC) by using ChIP-seq, ATAC-Seq, RNA-seq, and SNP arrays. We integrate 153 individual data sets from 42 patients and nominate 50 histology-specific master transcription factors (MTF) to define RCC histologic subtypes, including EPAS1 and ETS-1 in ccRCC, HNF1B in pRCC, and FOXI1 in chRCC. We confirm histology-specific MTFs via immunohistochemistry including a ccRCC-specific TF, BHLHE41. FOXI1 overexpression with knock-down of EPAS1 in the 786-O ccRCC cell line induces transcriptional upregulation of chRCC-specific genes, TFCP2L1, ATP6V0D2, KIT, and INSRR, implicating FOXI1 as a MTF for chRCC. Integrating RCC GWAS risk SNPs with H3K27ac ChIP-seq and ATAC-seq data reveals that risk-variants are significantly enriched in allelically-imbalanced peaks. This epigenomic atlas in primary human samples provides a resource for future investigation.

Circulating Cell-Free DNA in Renal Cell Carcinoma: The New Era of Precision Medicine.

Tumor biopsy is still the gold standard for diagnosing and prognosis renal cell carcinoma (RCC). However, its invasiveness, costs, and inability to accurately picture tumor heterogeneity represent major limitations to this procedure. Analysis of circulating cell-free DNA (cfDNA) is a non-invasive cost-effective technique that has the potential to ease cancer detection and prognosis. In particular, a growing body of evidence suggests that cfDNA could be a complementary tool to identify and prognosticate RCC while providing contemporary mutational profiling of the tumor. Further, recent research highlighted the role of cfDNA methylation profiling as a novel method for cancer detection and tissue-origin identification. This review synthesizes current knowledge on the diagnostic, prognostic, and predictive applications of cfDNA in RCC, with a specific focus on the potential role of cell-free methylated DNA (cfMeDNA).

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.

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).

A Somatically Acquired Enhancer of the Androgen Receptor Is a Noncoding Driver in Advanced Prostate Cancer.

Increased androgen receptor (AR) activity drives therapeutic resistance in advanced prostate cancer. The most common resistance mechanism is amplification of this locus presumably targeting the AR gene. Here, we identify and characterize a somatically acquired AR enhancer located 650 kb centromeric to the AR. Systematic perturbation of this enhancer using genome editing decreased proliferation by suppressing AR levels. Insertion of an additional copy of this region sufficed to increase proliferation under low androgen conditions and to decrease sensitivity to enzalutamide. Epigenetic data generated in localized prostate tumors and benign specimens support the notion that this region is a developmental enhancer. Collectively, these observations underscore the importance of epigenomic profiling in primary specimens and the value of deploying genome editing to functionally characterize noncoding elements. More broadly, this work identifies a therapeutic vulnerability for targeting the AR and emphasizes the importance of regulatory elements as highly recurrent oncogenic drivers.

Distinct Effects of p19 RNA Silencing Suppressor on Small RNA Mediated Pathways in Plants.

RNA silencing is one of the main defense mechanisms employed by plants to fight viruses. In change, viruses have evolved silencing suppressor proteins to neutralize antiviral silencing. Since the endogenous and antiviral functions of RNA silencing pathway rely on common components, it was suggested that viral suppressors interfere with endogenous silencing pathway contributing to viral symptom development. In this work, we aimed to understand the effects of the tombusviral p19 suppressor on endogenous and antiviral silencing during genuine virus infection. We showed that ectopically expressed p19 sequesters endogenous small RNAs (sRNAs) in the absence, but not in the presence of virus infection. Our presented data question the generalized model in which the sequestration of endogenous sRNAs by the viral suppressor contributes to the viral symptom development. We further showed that p19 preferentially binds the perfectly paired ds-viral small interfering RNAs (vsiRNAs) but does not select based on their sequence or the type of the 5' nucleotide. Finally, co-immunoprecipitation of sRNAs with AGO1 or AGO2 from virus-infected plants revealed that p19 specifically impairs vsiRNA loading into AGO1 but not AGO2. Our findings, coupled with the fact that p19-expressing wild type Cymbidium ringspot virus (CymRSV) overcomes the Nicotiana benthamiana silencing based defense killing the host, suggest that AGO1 is the main effector of antiviral silencing in this host-virus combination.

CAUSEL: an epigenome- and genome-editing pipeline for establishing function of noncoding GWAS variants.

The vast majority of disease-associated single-nucleotide polymorphisms (SNPs) mapped by genome-wide association studies (GWASs) are located in the non-protein-coding genome, but establishing the functional and mechanistic roles of these sequence variants has proven challenging. Here we describe a general pipeline in which candidate functional SNPs are first evaluated by fine mapping, epigenomic profiling, and epigenome editing, and then interrogated for causal function by using genome editing to create isogenic cell lines followed by phenotypic characterization. To validate this approach, we analyzed the 6q22.1 prostate cancer risk locus and identified rs339331 as the top-scoring SNP. Epigenome editing confirmed that the rs339331 region possessed regulatory potential. By using transcription activator-like effector nuclease (TALEN)-mediated genome editing, we created a panel of isogenic 22Rv1 prostate cancer cell lines representing all three genotypes (TT, TC, CC) at rs339331. Introduction of the 'T' risk allele increased transcription of the regulatory factor 6 (RFX6) gene, increased homeobox B13 (HOXB13) binding at the rs339331 region, and increased deposition of the enhancer-associated H3K4me2 histone mark at the rs339331 region compared to lines homozygous for the 'C' protective allele. The cell lines also differed in cellular morphology and adhesion, and pathway analysis of differentially expressed genes suggested an influence of androgens. In summary, we have developed and validated a widely accessible approach that can be used to establish functional causality for noncoding sequence variants identified by GWASs.