Comprehensive genomic analysis of a BRCA2 deficient human pancreatic cancer

Uncovering novel mutational signatures by de novo extraction with SigProfilerExtractor

Mutational signature analysis is commonly performed in cancer genomic studies. Here, we present SigProfilerExtractor, an automated tool for de novo extraction of mutational signatures, and benchmark it against another 13 bioinformatics tools by using 34 scenarios encompassing 2,500 simulated signatures found in 60,000 synthetic genomes and 20,000 synthetic exomes. For simulations with 5% noise, reflecting high-quality datasets, SigProfilerExtractor outperforms other approaches by elucidating between 20% and 50% more true-positive signatures while yielding 5-fold less false-positive signatures. Applying SigProfilerExtractor to 4,643 whole-genome- and 19,184 whole-exome-sequenced cancers reveals four novel signatures. Two of the signatures are confirmed in independent cohorts, and one of these signatures is associated with tobacco smoking. In summary, this report provides a reference tool for analysis of mutational signatures, a comprehensive benchmarking of bioinformatics tools for extracting signatures, and several novel mutational signatures, including one putatively attributed to direct tobacco smoking mutagenesis in bladder tissues.

A first-in-class Polymerase Theta Inhibitor selectively targets Homologous-Recombination-Deficient Tumors

DNA polymerase theta (POLθ) is synthetic lethal with Homologous Recombination (HR) deficiency and thus a candidate target for HR-deficient cancers. Through high-throughput small molecule screens we identified the antibiotic Novobiocin (NVB) as a specific POLθ inhibitor that selectively kills HR-deficient tumor cells in vitro and in vivo. NVB directly binds to the POLθ ATPase domain, inhibits its ATPase activity, and phenocopies POLθ depletion. NVB kills HR-deficient breast and ovarian tumors in GEMM, xenograft and PDX models. Increased POLθ levels predict NVB sensitivity, and BRCA-deficient tumor cells with acquired resistance to PARP inhibitors (PARPi) are sensitive to NVB in vitro and in vivo. Mechanistically, NVB-mediated cell death in PARPi-resistant cells arises from increased double-strand break end resection, leading to accumulation of single-strand DNA intermediates and non-functional RAD51 foci. Our results demonstrate that NVB may be useful alone or in combination with PARPi in treating HR-deficient tumors, including those with acquired PARPi resistance. (151/150).

A Preclinical Trial and Molecularly Annotated Patient Cohort Identify Predictive Biomarkers in Homologous Recombination–deficient Pancreatic Cancer

Purpose:

Pancreatic ductal adenocarcinoma (PDAC) arising in patients with a germline BRCA1 or BRCA2 (gBRCA) mutation may be sensitive to platinum and PARP inhibitors (PARPi). However, treatment stratification based on gBRCA mutational status alone is associated with heterogeneous responses.

Experimental Design:

We performed a seven-arm preclinical trial consisting of 471 mice, representing 12 unique PDAC patient-derived xenografts, of which nine were gBRCA mutated. From 179 patients whose PDAC was whole-genome and transcriptome sequenced, we identified 21 cases with homologous recombination deficiency (HRD), and investigated prognostic biomarkers.

Results:

We found that biallelic inactivation of BRCA1/BRCA2 is associated with genomic hallmarks of HRD and required for cisplatin and talazoparib (PARPi) sensitivity. However, HRD genomic hallmarks persisted in xenografts despite the emergence of therapy resistance, indicating the presence of a genomic scar. We identified tumor polyploidy and a low Ki67 index as predictors of poor cisplatin and talazoparib response. In patients with HRD PDAC, tumor polyploidy and a basal-like transcriptomic subtype were independent predictors of shorter survival. To facilitate clinical assignment of transcriptomic subtype, we developed a novel pragmatic two-marker assay (GATA6:KRT17).

Conclusions:

In summary, we propose a predictive and prognostic model of gBRCA-mutated PDAC on the basis of HRD genomic hallmarks, Ki67 index, tumor ploidy, and transcriptomic subtype.

Strong suppression of gene conversion with increasing DNA double-strand break load delimited by 53BP1 and RAD52

In vertebrates, genomic DNA double-strand breaks (DSBs) are removed by non-homologous end-joining processes: classical non-homologous end-joining (c-NHEJ) and alternative end-joining (alt-EJ); or by homology-dependent processes: gene-conversion (GC) and single-strand annealing (SSA). Surprisingly, these repair pathways are not real alternative options restoring genome integrity with equal efficiency, but show instead striking differences in speed, accuracy and cell-cycle-phase dependence. As a consequence, engagement of one pathway may be associated with processing-risks for the genome absent from another pathway. Characterization of engagement-parameters and their consequences is, therefore, essential for understanding effects on the genome of DSB-inducing agents, such as ionizing-radiation (IR). Here, by addressing pathway selection in G2-phase, we discover regulatory confinements in GC with consequences for SSA- and c-NHEJ-engagement. We show pronounced suppression of GC with increasing DSB-load that is not due to RAD51 availability and which is delimited but not defined by 53BP1 and RAD52. Strikingly, at low DSB-loads, GC repairs ∼50% of DSBs, whereas at high DSB-loads its contribution is undetectable. Notably, with increasing DSB-load and the associated suppression of GC, SSA gains ground, while alt-EJ is suppressed. These observations explain earlier, apparently contradictory results and advance our understanding of logic and mechanisms underpinning the wiring between DSB repair pathways.

The spectrum of genetic variants in hereditary pancreatic cancer includes Fanconi anemia genes

Approximately 5-10% of all pancreatic cancer patients carry a predisposing mutation in a known susceptibility gene. Since >90% of patients present with late stage disease, it is crucial to identify high risk individuals who may be amenable to early detection or other prevention. To explore the spectrum of hereditary pancreatic cancer susceptibility, we evaluated germline DNA from pancreatic cancer participants (n = 53) from a large hereditary cancer registry. For those without a known predisposition mutation gene (n = 49), germline next generation sequencing was completed using targeted capture for 706 candidate genes. We identified 16 of 53 participants (30%) with a pathogenic (P) or likely pathogenic (LP) variant that may be related to their hereditary pancreatic cancer predisposition; seven had mutations in genes associated with well-known cancer syndromes (13%) [ATM (2), BRCA2 (3), MSH2 (1), MSH6 (1)]. Many had mutations in Fanconi anemia complex genes [BRCA2 (3 participants), FANCF, FANCM]. Eight participants had rare protein truncating variants of uncertain significance with no other P or LP variants. Earlier age of pancreatic cancer diagnosis (57.5 vs 64.8 years) was indicative of possessing a P or LP variant, as was cancer family history (p values <0.0001). Our multigene panel approach for identifying known cancer predisposing genetic susceptibility in those at risk for hereditary pancreatic cancer may have direct applicability to clinical practice in cases with mutations in actionable genes. Future pancreatic cancer predisposition studies should include evaluation of the Fanconi anemia genes.

Biomarkers and pathways of chemoresistance and chemosensitivity for personalized treatment of pancreatic adenocarcinoma

Pancreatic carcinoma is usually diagnosed late when treatment options are limited and is considered a chemo-resistant malignancy. However, early stage, good performance status and specific patient subgroup are thought to have a more favorable prognosis. Search for novel molecular biomarkers, which could predict treatment resistance, represents a major opportunity, but also a challenge in further research. This review summarizes most aspects of individualized therapy of pancreatic cancer including promising biomarkers, BRCA-deficient pancreatic cancer and its etiology. It may be estimated that nearly a third of metastatic pancreatic ductal adenocarcinoma patients could benefit from treatment other than gold standard chemotherapy. Thus, other aspects of an individualized approach concerning the main factors for the choice of the best therapy for individual pancreatic cancer patient (surgery and chemotherapy), as well as the future directions (target therapy and immunotherapy), are also addressed.

TGF-βRII Knock-down in Pancreatic Cancer Cells Promotes Tumor Growth and Gemcitabine Resistance. Importance of STAT3 Phosphorylation on S727

Pancreatic adenocarcinoma (PDAC) is one of the most deadly cancers in the Western world because of a lack of early diagnostic markers and efficient therapeutics. At the time of diagnosis, more than 80% of patients have metastasis or locally advanced cancer and are therefore not eligible for surgical resection. Pancreatic cancer cells also harbour a high resistance to chemotherapeutic drugs such as gemcitabine that is one of the main palliative treatments for PDAC. Proteins involved in TGF-β signaling pathway (SMAD4 or TGF-βRII) are frequently mutated in PDAC (50–80%). TGF-β signalling pathway plays antagonistic roles during carcinogenesis by initially inhibiting epithelial growth and later promoting the progression of advanced tumors and thus emerged as both tumor suppressor and oncogenic pathways. In order to decipher the role of TGF-β in pancreatic carcinogenesis and chemoresistance, we generated CAPAN-1 and CAPAN-2 cell lines knocked down for TGF-βRII (first actor of TGF-β signaling). The impact on biological properties of these TGF-βRII-KD cells was studied both in vitro and in vivo. We show that TGF-βRII silencing alters tumor growth and migration as well as resistance to gemcitabine. TGF-βRII silencing also leads to S727 STAT3 and S63 c-Jun phosphorylation, decrease of MRP3 and increase of MRP4 ABC transporter expression and induction of a partial EMT phenotype. These markers associated with TGF-β signaling pathways may thus appear as potent therapeutic tools to better treat/manage pancreatic cancer.

Distinct binding of BRCA2 BRC repeats to RAD51 generates differential DNA damage sensitivity

BRCA2 is a multi-faceted protein critical for the proper regulation of homology-directed repair of DNA double-strand breaks. Elucidating the mechanistic features of BRCA2 is crucial for understanding homologous recombination and how patient-derived mutations impact future cancer risk. Eight centrally located BRC repeats in BRCA2 mediate binding and regulation of RAD51 on resected DNA substrates. Herein, we dissect the biochemical and cellular features of the BRC repeats tethered to the DNA binding domain of BRCA2. To understand how the BRC repeats and isolated domains of BRCA2 contribute to RAD51 binding, we analyzed both the biochemical and cellular properties of these proteins. In contrast to the individual BRC repeat units, we find that the BRC5-8 region potentiates RAD51-mediated DNA strand pairing and provides complementation functions exceeding those of BRC repeats 1-4. Furthermore, BRC5-8 can efficiently repair nuclease-induced DNA double-strand breaks and accelerate the assembly of RAD51 repair complexes upon DNA damage. These findings highlight the importance of the BRC5-8 domain in stabilizing the RAD51 filament and promoting homology-directed repair under conditions of cellular DNA damage.

Whole-exome DNA sequence analysis of Brca2- and Trp53-deficient mouse mammary gland tumours

Germline mutations in the tumour suppressor BRCA2 predispose to breast, ovarian and a number of other human cancers. Brca2-deficient mouse models are used for preclinical studies but the pattern of genomic alterations in these tumours has not yet been described in detail. We have performed whole-exome DNA sequencing analysis of mouse mammary tumours from Blg-Cre Brca2(f/f) Trp53(f/f) animals, a model of BRCA2-deficient human cancer. We also used the sequencing data to estimate DNA copy number alterations in these tumours and identified a recurrent copy number gain in Met, which has been found amplified in other mouse mammary cancer models. Through a comparative genomic analysis, we identified several mouse Blg-Cre Brca2(f/f) Trp53(f/f) mammary tumour somatic mutations in genes that are also mutated in human cancer, but few of these genes have been found frequently mutated in human breast cancer. A more detailed analysis of these somatic mutations revealed a set of genes that are mutated in human BRCA2 mutant breast and ovarian tumours and that are also mutated in mouse Brca2-null, Trp53-null mammary tumours. Finally, a DNA deletion surrounded by microhomology signature found in human BRCA1/2-deficient cancers was not common in the genome of these mouse tumours. Although a useful model, there are some differences in the genomic landscape of tumours arising in Blg-Cre Brca2(f/f) Trp53(f/f) mice compared to human BRCA-mutated breast cancers. Therefore, this needs to be taken into account in the use of this model.

TET proteins and the control of cytosine demethylation in cancer

The discovery that ten-eleven translocation (TET) proteins are α-ketoglutarate-dependent dioxygenases involved in the conversion of 5-methylcytosines (5-mC) to 5-hydroxymethylcytosine (5-hmC), 5-formylcytosine and 5-carboxycytosine has revealed new pathways in the cytosine methylation and demethylation process. The description of inactivating mutations in TET2 suggests that cellular transformation is in part caused by the deregulation of this 5-mC conversion. The direct and indirect deregulation of methylation control through mutations in DNA methyltransferase and isocitrate dehydrogenase (IDH) genes, respectively, along with the importance of cytosine methylation in the control of normal and malignant cellular differentiation have provided a conceptual framework for understanding the early steps in cancer development. Here, we review recent advances in our understanding of the cytosine methylation cycle and its implication in cellular transformation, with an emphasis on TET enzymes and 5-hmC. Ongoing clinical trials targeting the activity of mutated IDH enzymes provide a proof of principle that DNA methylation is targetable, and will trigger further therapeutic applications aimed at controlling both early and late stages of cancer development.

Exome-wide somatic microsatellite variation is altered in cells with DNA repair deficiencies

Microsatellites (MST), tandem repeats of 1–6 nucleotide motifs, are mutational hot-spots with a bias for insertions and deletions (INDELs) rather than single nucleotide polymorphisms (SNPs). The majority of MST instability studies are limited to a small number of loci, the Bethesda markers, which are only informative for a subset of colorectal cancers. In this paper we evaluate non-haplotype alleles present within next-gen sequencing data to evaluate somatic MST variation (SMV) within DNA repair proficient and DNA repair defective cell lines. We confirm that alleles present within next-gen data that do not contribute to the haplotype can be reliably quantified and utilized to evaluate the SMV without requiring comparisons of matched samples. We observed that SMV patterns found in DNA repair proficient cell lines without DNA repair defects, MCF10A, HEK293 and PD20 RV:D2, had consistent patterns among samples. Further, we were able to confirm that changes in SMV patterns in cell lines lacking functional BRCA2, FANCD2 and mismatch repair were consistent with the different pathways perturbed. Using this new exome sequencing analysis approach we show that DNA instability can be identified in a sample and that patterns of instability vary depending on the impaired DNA repair mechanism, and that genes harboring minor alleles are strongly associated with cancer pathways. The MST Minor Allele Caller used for this study is available at https://github.com/zalmanv/MST_minor_allele_caller.

A whole-genome massively parallel sequencing analysis of BRCA1 mutant oestrogen receptor-negative and -positive breast cancers

BRCA1 encodes a tumour suppressor protein that plays pivotal roles in homologous recombination (HR) DNA repair, cell-cycle checkpoints, and transcriptional regulation. BRCA1 germline mutations confer a high risk of early-onset breast and ovarian cancer. In more than 80% of cases, tumours arising in BRCA1 germline mutation carriers are oestrogen receptor (ER)-negative; however, up to 15% are ER-positive. It has been suggested that BRCA1 ER-positive breast cancers constitute sporadic cancers arising in the context of a BRCA1 germline mutation rather than being causally related to BRCA1 loss-of-function. Whole-genome massively parallel sequencing of ER-positive and ER-negative BRCA1 breast cancers, and their respective germline DNAs, was used to characterize the genetic landscape of BRCA1 cancers at base-pair resolution. Only BRCA1 germline mutations, somatic loss of the wild-type allele, and TP53 somatic mutations were recurrently found in the index cases. BRCA1 breast cancers displayed a mutational signature consistent with that caused by lack of HR DNA repair in both ER-positive and ER-negative cases. Sequencing analysis of independent cohorts of hereditary BRCA1 and sporadic non-BRCA1 breast cancers for the presence of recurrent pathogenic mutations and/or homozygous deletions found in the index cases revealed that DAPK3, TMEM135, KIAA1797, PDE4D, and GATA4 are potential additional drivers of breast cancers. This study demonstrates that BRCA1 pathogenic germline mutations coupled with somatic loss of the wild-type allele are not sufficient for hereditary breast cancers to display an ER-negative phenotype, and has led to the identification of three potential novel breast cancer genes (ie DAPK3, TMEM135, and GATA4).

C. A modified method for whole exome resequencing from minimal amounts of starting DNA

Next generation DNA sequencing (NGS) technologies have revolutionized the pace at which whole genome and exome sequences can be generated. However, despite these advances, many of the methods for targeted resequencing, such as the generation of high-depth exome sequences, are somewhat limited by the relatively large amounts of starting DNA that are normally required. In the case of tumour analysis this is particularly pertinent as many tumour biopsies often return submicrogram quantities of DNA, especially when tumours are microdissected prior to analysis. Here, we present a method for exome capture and resequencing using as little as 50 ng of starting DNA. The sequencing libraries generated by this minimal starting amount (MSA-Cap) method generate datasets that are comparable to standard amount (SA) whole exome libraries that use three micrograms of starting DNA. This method, which can be performed in most laboratories using commonly available reagents, has the potential to enhance large scale profiling efforts such as the resequencing of tumour exomes.

The DNA damage response and cancer therapy

Genomic instability is one of the most pervasive characteristics of tumour cells and is probably the combined effect of DNA damage, tumour-specific DNA repair defects, and a failure to stop or stall the cell cycle before the damaged DNA is passed on to daughter cells. Although these processes drive genomic instability and ultimately the disease process, they also provide therapeutic opportunities. A better understanding of the cellular response to DNA damage will not only inform our knowledge of cancer development but also help to refine the classification as well as the treatment of the disease.