Multi-omics analysis reveals distinct non-reversion mechanisms of PARPi resistance in BRCA1- versus BRCA2-deficient mammary tumors.

BRCA1 and BRCA2 both function in DNA double-strand break repair by homologous recombination (HR). Due to their HR defect, BRCA1/2-deficient cancers are sensitive to poly(ADP-ribose) polymerase inhibitors (PARPis), but they eventually acquire resistance. Preclinical studies yielded several PARPi resistance mechanisms that do not involve BRCA1/2 reactivation, but their relevance in the clinic remains elusive. To investigate which BRCA1/2-independent mechanisms drive spontaneous resistance in vivo, we combine molecular profiling with functional analysis of HR of matched PARPi-naive and PARPi-resistant mouse mammary tumors harboring large intragenic deletions that prevent reactivation of BRCA1/2. We observe restoration of HR in 62% of PARPi-resistant BRCA1-deficient tumors but none in the PARPi-resistant BRCA2-deficient tumors. Moreover, we find that 53BP1 loss is the prevalent resistance mechanism in HR-proficient BRCA1-deficient tumors, whereas resistance in BRCA2-deficient tumors is mainly induced by PARG loss. Furthermore, combined multi-omics analysis identifies additional genes and pathways potentially involved in modulating PARPi response.

Truncated FGFR2 is a clinically actionable oncogene in multiple cancers.

Somatic hotspot mutations and structural amplifications and fusions that affect fibroblast growth factor receptor 2 (encoded by FGFR2) occur in multiple types of cancer1. However, clinical responses to FGFR inhibitors have remained variable1-9, emphasizing the need to better understand which FGFR2 alterations are oncogenic and therapeutically targetable. Here we apply transposon-based screening10,11 and tumour modelling in mice12,13, and find that the truncation of exon 18 (E18) of Fgfr2 is a potent driver mutation. Human oncogenomic datasets revealed a diverse set of FGFR2 alterations, including rearrangements, E1-E17 partial amplifications, and E18 nonsense and frameshift mutations, each causing the transcription of E18-truncated FGFR2 (FGFR2ΔE18). Functional in vitro and in vivo examination of a compendium of FGFR2ΔE18 and full-length variants pinpointed FGFR2-E18 truncation as single-driver alteration in cancer. By contrast, the oncogenic competence of FGFR2 full-length amplifications depended on a distinct landscape of cooperating driver genes. This suggests that genomic alterations that generate stable FGFR2ΔE18 variants are actionable therapeutic targets, which we confirmed in preclinical mouse and human tumour models, and in a clinical trial. We propose that cancers containing any FGFR2 variant with a truncated E18 should be considered for FGFR-targeted therapies.

Comparative oncogenomics identifies combinations of driver genes and drug targets in BRCA1-mutated breast cancer.

BRCA1-mutated breast cancer is primarily driven by DNA copy-number alterations (CNAs) containing large numbers of candidate driver genes. Validation of these candidates requires novel approaches for high-throughput in vivo perturbation of gene function. Here we develop genetically engineered mouse models (GEMMs) of BRCA1-deficient breast cancer that permit rapid introduction of putative drivers by either retargeting of GEMM-derived embryonic stem cells, lentivirus-mediated somatic overexpression or in situ CRISPR/Cas9-mediated gene disruption. We use these approaches to validate Myc, Met, Pten and Rb1 as bona fide drivers in BRCA1-associated mammary tumorigenesis. Iterative mouse modeling and comparative oncogenomics analysis show that MYC-overexpression strongly reshapes the CNA landscape of BRCA1-deficient mammary tumors and identify MCL1 as a collaborating driver in these tumors. Moreover, MCL1 inhibition potentiates the in vivo efficacy of PARP inhibition (PARPi), underscoring the therapeutic potential of this combination for treatment of BRCA1-mutated cancer patients with poor response to PARPi monotherapy.

XenofilteR: computational deconvolution of mouse and human reads in tumor xenograft sequence data.

BACKGROUND: Mouse xenografts from (patient-derived) tumors (PDX) or tumor cell lines are widely used as models to study various biological and preclinical aspects of cancer. However, analyses of their RNA and DNA profiles are challenging, because they comprise reads not only from the grafted human cancer but also from the murine host. The reads of murine origin result in false positives in mutation analysis of DNA samples and obscure gene expression levels when sequencing RNA. However, currently available algorithms are limited and improvements in accuracy and ease of use are necessary. RESULTS: We developed the R-package XenofilteR, which separates mouse from human sequence reads based on the edit-distance between a sequence read and reference genome. To assess the accuracy of XenofilteR, we generated sequence data by in silico mixing of mouse and human DNA sequence data. These analyses revealed that XenofilteR removes > 99.9% of sequence reads of mouse origin while retaining human sequences. This allowed for mutation analysis of xenograft samples with accurate variant allele frequencies, and retrieved all non-synonymous somatic tumor mutations. CONCLUSIONS: XenofilteR accurately dissects RNA and DNA sequences from mouse and human origin, thereby outperforming currently available tools. XenofilteR is open source and available at https://github.com/PeeperLab/XenofilteR .

Transcriptomics and Transposon Mutagenesis Identify Multiple Mechanisms of Resistance to the FGFR Inhibitor AZD4547.

In human cancers, FGFR signaling is frequently hyperactivated by deregulation of FGF ligands or by activating mutations in the FGFR receptors such as gene amplifications, point mutations, and gene fusions. As such, FGFR inhibitors are considered an attractive therapeutic strategy for patients with mutations in FGFR family members. We previously identified Fgfr2 as a key driver of invasive lobular carcinoma (ILC) in an in vivo insertional mutagenesis screen using the Sleeping Beauty transposon system. Here we explore whether these FGFR-driven ILCs are sensitive to the FGFR inhibitor AZD4547 and use transposon mutagenesis in these tumors to identify potential mechanisms of resistance to therapy. Combined with RNA sequencing-based analyses of AZD4547-resistant tumors, our in vivo approach identified several known and novel potential resistance mechanisms to FGFR inhibition, most of which converged on reactivation of the canonical MAPK-ERK signaling cascade. Observed resistance mechanisms included mutations in the tyrosine kinase domain of FGFR2, overexpression of MET, inactivation of RASA1, and activation of the drug-efflux transporter ABCG2. ABCG2 and RASA1 were identified only from de novo transposon insertions acquired during AZD4547 treatment, demonstrating that insertional mutagenesis in mice is an effective tool for identifying potential mechanisms of resistance to targeted cancer therapies.Significance: These findings demonstrate that a combined approach of transcriptomics and insertional mutagenesis in vivo is an effective method for identifying potential targets to overcome resistance to therapy in the clinic. Cancer Res; 78(19); 5668-79. ©2018 AACR.

Selective Loss of PARG Restores PARylation and Counteracts PARP Inhibitor-Mediated Synthetic Lethality.

Inhibitors of poly(ADP-ribose) (PAR) polymerase (PARPi) have recently entered the clinic for the treatment of homologous recombination (HR)-deficient cancers. Despite the success of this approach, drug resistance is a clinical hurdle, and we poorly understand how cancer cells escape the deadly effects of PARPi without restoring the HR pathway. By combining genetic screens with multi-omics analysis of matched PARPi-sensitive and -resistant Brca2-mutated mouse mammary tumors, we identified loss of PAR glycohydrolase (PARG) as a major resistance mechanism. We also found the presence of PARG-negative clones in a subset of human serous ovarian and triple-negative breast cancers. PARG depletion restores PAR formation and partially rescues PARP1 signaling. Importantly, PARG inactivation exposes vulnerabilities that can be exploited therapeutically.

Identifying transposon insertions and their effects from RNA-sequencing data.

Insertional mutagenesis using engineered transposons is a potent forward genetic screening technique used to identify cancer genes in mouse model systems. In the analysis of these screens, transposon insertion sites are typically identified by targeted DNA-sequencing and subsequently assigned to predicted target genes using heuristics. As such, these approaches provide no direct evidence that insertions actually affect their predicted targets or how transcripts of these genes are affected. To address this, we developed IM-Fusion, an approach that identifies insertion sites from gene-transposon fusions in standard single- and paired-end RNA-sequencing data. We demonstrate IM-Fusion on two separate transposon screens of 123 mammary tumors and 20 B-cell acute lymphoblastic leukemias, respectively. We show that IM-Fusion accurately identifies transposon insertions and their true target genes. Furthermore, by combining the identified insertion sites with expression quantification, we show that we can determine the effect of a transposon insertion on its target gene(s) and prioritize insertions that have a significant effect on expression. We expect that IM-Fusion will significantly enhance the accuracy of cancer gene discovery in forward genetic screens and provide initial insight into the biological effects of insertions on candidate cancer genes.

Insertional mutagenesis identifies drivers of a novel oncogenic pathway in invasive lobular breast carcinoma.

Invasive lobular carcinoma (ILC) is the second most common breast cancer subtype and accounts for 8-14% of all cases. Although the majority of human ILCs are characterized by the functional loss of E-cadherin (encoded by CDH1), inactivation of Cdh1 does not predispose mice to develop mammary tumors, implying that mutations in additional genes are required for ILC formation in mice. To identify these genes, we performed an insertional mutagenesis screen using the Sleeping Beauty transposon system in mice with mammary-specific inactivation of Cdh1. These mice developed multiple independent mammary tumors of which the majority resembled human ILC in terms of morphology and gene expression. Recurrent and mutually exclusive transposon insertions were identified in Myh9, Ppp1r12a, Ppp1r12b and Trp53bp2, whose products have been implicated in the regulation of the actin cytoskeleton. Notably, MYH9, PPP1R12B and TP53BP2 were also frequently aberrated in human ILC, highlighting these genes as drivers of a novel oncogenic pathway underlying ILC development.

PTEN Loss in E-Cadherin-Deficient Mouse Mammary Epithelial Cells Rescues Apoptosis and Results in Development of Classical Invasive Lobular Carcinoma.

Invasive lobular carcinoma (ILC) is an aggressive breast cancer subtype with poor response to chemotherapy. Besides loss of E-cadherin, a hallmark of ILC, genetic inactivation of PTEN is frequently observed in patients. Through concomitant Cre-mediated inactivation of E-cadherin and PTEN in mammary epithelium, we generated a mouse model of classical ILC (CLC), the main histological ILC subtype. While loss of E-cadherin induced cell dissemination and apoptosis, additional PTEN inactivation promoted cell survival and rapid formation of invasive mammary tumors that recapitulate the histological and molecular features, estrogen receptor (ER) status, growth kinetics, metastatic behavior, and tumor microenvironment of human CLC. Combined inactivation of E-cadherin and PTEN is sufficient to cause CLC development. These CLCs showed significant tumor regression upon BEZ235-mediated inhibition of PI3K signaling. In summary, this mouse model provides important insights into CLC development and suggests inhibition of phosphatidylinositol 3-kinase (PI3K) signaling as a potential therapeutic strategy for targeting CLC.

Replication fork stability confers chemoresistance in BRCA-deficient cells.

Cells deficient in the Brca1 and Brca2 genes have reduced capacity to repair DNA double-strand breaks by homologous recombination and consequently are hypersensitive to DNA-damaging agents, including cisplatin and poly(ADP-ribose) polymerase (PARP) inhibitors. Here we show that loss of the MLL3/4 complex protein, PTIP, protects Brca1/2-deficient cells from DNA damage and rescues the lethality of Brca2-deficient embryonic stem cells. However, PTIP deficiency does not restore homologous recombination activity at double-strand breaks. Instead, its absence inhibits the recruitment of the MRE11 nuclease to stalled replication forks, which in turn protects nascent DNA strands from extensive degradation. More generally, acquisition of PARP inhibitors and cisplatin resistance is associated with replication fork protection in Brca2-deficient tumour cells that do not develop Brca2 reversion mutations. Disruption of multiple proteins, including PARP1 and CHD4, leads to the same end point of replication fork protection, highlighting the complexities by which tumour cells evade chemotherapeutic interventions and acquire drug resistance.