Identification of genes that promote PI3K pathway activation and prostate tumour formation.

We have performed a functional in vivo mutagenesis screen to identify genes that, when altered, cooperate with a heterozygous Pten mutation to promote prostate tumour formation. Two genes, Bzw2 and Eif5a2, which have been implicated in the process of protein translation, were selected for further validation. Using prostate organoid models, we show that either Bzw2 downregulation or EIF5A2 overexpression leads to increased organoid size and in vivo prostate growth. We show that both genes impact the PI3K pathway and drive a sustained increase in phospho-AKT expression, with PTEN protein levels reduced in both models. Mechanistic studies reveal that EIF5A2 is directly implicated in PTEN protein translation. Analysis of patient datasets identified EIF5A2 amplifications in many types of human cancer, including the prostate. Human prostate cancer samples in two independent cohorts showed a correlation between increased levels of EIF5A2 and upregulation of a PI3K pathway gene signature. Consistent with this, organoids with high levels of EIF5A2 were sensitive to AKT inhibitors. Our study identified novel genes that promote prostate cancer formation through upregulation of the PI3K pathway, predicting a strategy to treat patients with genetic aberrations in these genes particularly relevant for EIF5A2 amplified tumours.

Transposon Mutagenesis Reveals RBMS3 Silencing as a Promoter of Malignant Progression of BRAFV600E-Driven Lung Tumorigenesis.

Mutationally activated BRAF is detected in approximately 7% of human lung adenocarcinomas, with BRAFT1799A serving as a predictive biomarker for treatment of patients with FDA-approved inhibitors of BRAFV600E oncoprotein signaling. In genetically engineered mouse (GEM) models, expression of BRAFV600E in the lung epithelium initiates growth of benign lung tumors that, without additional genetic alterations, rarely progress to malignant lung adenocarcinoma. To identify genes that cooperate with BRAFV600E for malignant progression, we used Sleeping Beauty-mediated transposon mutagenesis, which dramatically accelerated the emergence of lethal lung cancers. Among the genes identified was Rbms3, which encodes an RNA-binding protein previously implicated as a putative tumor suppressor. Silencing of RBMS3 via CRISPR/Cas9 gene editing promoted growth of BRAFV600E lung organoids and promoted development of malignant lung cancers with a distinct micropapillary architecture in BRAFV600E and EGFRL858R GEM models. BRAFV600E/RBMS3Null lung tumors displayed elevated expression of Ctnnb1, Ccnd1, Axin2, Lgr5, and c-Myc mRNAs, suggesting that RBMS3 silencing elevates signaling through the WNT/ß-catenin signaling axis. Although RBMS3 silencing rendered BRAFV600E-driven lung tumors resistant to the effects of dabrafenib plus trametinib, the tumors were sensitive to inhibition of porcupine, an acyltransferase of WNT ligands necessary for their secretion. Analysis of The Cancer Genome Atlas patient samples revealed that chromosome 3p24, which encompasses RBMS3, is frequently lost in non-small cell lung cancer and correlates with poor prognosis. Collectively, these data reveal the role of RBMS3 as a lung cancer suppressor and suggest that RBMS3 silencing may contribute to malignant NSCLC progression. SIGNIFICANCE: Loss of RBMS3 cooperates with BRAFV600E to induce lung tumorigenesis, providing a deeper understanding of the molecular mechanisms underlying mutant BRAF-driven lung cancer and potential strategies to more effectively target this disease.

Engineering Cancer Antigen-Specific T Cells to Overcome the Immunosuppressive Effects of TGF-ß.

Adoptive T cell therapy with T cells expressing affinity-enhanced TCRs has shown promising results in phase 1/2 clinical trials for solid and hematological tumors. However, depth and durability of responses to adoptive T cell therapy can suffer from an inhibitory tumor microenvironment. A common immune-suppressive agent is TGF-ß, which is secreted by tumor cells and cells recruited to the tumor. We investigated whether human T cells could be engineered to be resistant to inhibition by TGF-ß. Truncating the intracellular signaling domain from TGF-ß receptor (TGFßR) II produces a dominant-negative receptor (dnTGFßRII) that dimerizes with endogenous TGFßRI to form a receptor that can bind TGF-ß but cannot signal. We previously generated specific peptide enhanced affinity receptor TCRs recognizing the HLA-A*02-restricted peptides New York esophageal squamous cell carcinoma 1 (NY-ESO-1)157-165/l-Ag family member-1A (TCR: GSK3377794, formerly NY-ESO-1c259) and melanoma Ag gene A10254-262 (TCR: ADP-A2M10, formerly melanoma Ag gene A10c796). In this article, we show that exogenous TGF-ß inhibited in vitro proliferation and effector functions of human T cells expressing these first-generation high-affinity TCRs, whereas inhibition was reduced or abolished in the case of second-generation TCRs coexpressed with dnTGFßRII (e.g., GSK3845097). TGF-ß isoforms and a panel of TGF-ß-associated genes are overexpressed in a range of cancer indications in which NY-ESO-1 is commonly expressed, particularly in synovial sarcoma. As an example, immunohistochemistry/RNAscope identified TGF-ß-positive cells close to T cells in tumor nests and stroma, which had low frequencies of cells expressing IFN-γ in a non-small cell lung cancer setting. Coexpression of dnTGFßRII may therefore improve the efficacy of TCR-transduced T cells.

Correction: Suppression of interferon gene expression overcomes resistance to MEK inhibition in KRAS-mutant colorectal cancer.

A correction to this paper has been published and can be accessed via a link at the top of the paper.

Suppression of interferon gene expression overcomes resistance to MEK inhibition in KRAS-mutant colorectal cancer.

Despite showing clinical activity in BRAF-mutant melanoma, the MEK inhibitor (MEKi) trametinib has failed to show clinical benefit in KRAS-mutant colorectal cancer. To identify mechanisms of resistance to MEKi, we employed a pharmacogenomic analysis of MEKi-sensitive versus MEKi-resistant colorectal cancer cell lines. Strikingly, interferon- and inflammatory-related gene sets were enriched in cell lines exhibiting intrinsic and acquired resistance to MEK inhibition. The bromodomain inhibitor JQ1 suppressed interferon-stimulated gene (ISG) expression and in combination with MEK inhibitors displayed synergistic effects and induced apoptosis in MEKi-resistant colorectal cancer cell lines. ISG expression was confirmed in patient-derived organoid models, which displayed resistance to trametinib and were resensitized by JQ1 co-treatment. In in vivo models of colorectal cancer, combination treatment significantly suppressed tumor growth. Our findings provide a novel explanation for the limited response to MEK inhibitors in KRAS-mutant colorectal cancer, known for its inflammatory nature. Moreover, the high expression of ISGs was associated with significantly reduced survival of colorectal cancer patients. Excitingly, we have identified novel therapeutic opportunities to overcome intrinsic and acquired resistance to MEK inhibition in colorectal cancer.

Optimised ARID1A immunohistochemistry is an accurate predictor of ARID1A mutational status in gynaecological cancers.

ARID1A is a tumour suppressor gene that is frequently mutated in clear cell and endometrioid carcinomas of the ovary and endometrium and is an important clinical biomarker for novel treatment approaches for patients with ARID1A defects. However, the accuracy of ARID1A immunohistochemistry (IHC) as a surrogate for mutation status has not fully been established for patient stratification in clinical trials. Here we tested whether ARID1A IHC could reliably predict ARID1A mutations identified by next-generation sequencing. Three commercially available antibodies - EPR13501 (Abcam), D2A8U (Cell Signaling), and HPA005456 (Sigma) - were optimised for IHC using cell line models and human tissue, and screened across a cohort of 45 gynaecological tumours. IHC was scored independently by three pathologists using an immunoreactive score. ARID1A mutation status was assessed using two independent sequencing platforms and the concordance between ARID1A mutation and protein expression was evaluated using Receiver Operating Characteristic statistics. Overall, 21 ARID1A mutations were identified in 14/43 assessable tumours (33%), the majority of which were predicted to be deleterious. Mutations were identified in 6/17 (35%) ovarian clear cell carcinomas, 5/8 (63%) ovarian endometrioid carcinomas, 2/5 (40%) endometrial carcinomas, and 1/7 (14%) carcinosarcomas. ROC analysis identified greater than 95% concordance between mutation status and IHC using a modified immunoreactive score for all three antibodies allowing a definitive cut-point for ARID1A mutant status to be calculated. Comprehensive assessment of concordance of ARID1A IHC and mutation status identified EPR13501 as an optimal antibody, with 100% concordance between ARID1A mutation status and protein expression, across different gynaecological histological subtypes. It delivered the best inter-rater agreement between all pathologists, as well as a clear cost-benefit advantage. This could allow patients to be accurately stratified based on their ARID1A IHC status into early phase clinical trials.

Capture Hi-C identifies putative target genes at 33 breast cancer risk loci.

Genome-wide association studies (GWAS) have identified approximately 100 breast cancer risk loci. Translating these findings into a greater understanding of the mechanisms that influence disease risk requires identification of the genes or non-coding RNAs that mediate these associations. Here, we use Capture Hi-C (CHi-C) to annotate 63 loci; we identify 110 putative target genes at 33 loci. To assess the support for these target genes in other data sources we test for associations between levels of expression and SNP genotype (eQTLs), disease-specific survival (DSS), and compare them with somatically mutated cancer genes. 22 putative target genes are eQTLs, 32 are associated with DSS and 14 are somatically mutated in breast, or other, cancers. Identifying the target genes at GWAS risk loci will lead to a greater understanding of the mechanisms that influence breast cancer risk and prognosis.

Transposon mutagenesis identifies genes that cooperate with mutant Pten in breast cancer progression.

Triple-negative breast cancer (TNBC) has the worst prognosis of any breast cancer subtype. To better understand the genetic forces driving TNBC, we performed a transposon mutagenesis screen in a phosphatase and tensin homolog (Pten) mutant mice and identified 12 candidate trunk drivers and a much larger number of progression genes. Validation studies identified eight TNBC tumor suppressor genes, including the GATA-like transcriptional repressor TRPS1 Down-regulation of TRPS1 in TNBC cells promoted epithelial-to-mesenchymal transition (EMT) by deregulating multiple EMT pathway genes, in addition to increasing the expression of SERPINE1 and SERPINB2 and the subsequent migration, invasion, and metastasis of tumor cells. Transposon mutagenesis has thus provided a better understanding of the genetic forces driving TNBC and discovered genes with potential clinical importance in TNBC.

Analyzing tumor heterogeneity and driver genes in single myeloid leukemia cells with SBCapSeq.

A central challenge in oncology is how to kill tumors containing heterogeneous cell populations defined by different combinations of mutated genes. Identifying these mutated genes and understanding how they cooperate requires single-cell analysis, but current single-cell analytic methods, such as PCR-based strategies or whole-exome sequencing, are biased, lack sequencing depth or are cost prohibitive. Transposon-based mutagenesis allows the identification of early cancer drivers, but current sequencing methods have limitations that prevent single-cell analysis. We report a liquid-phase, capture-based sequencing and bioinformatics pipeline, Sleeping Beauty (SB) capture hybridization sequencing (SBCapSeq), that facilitates sequencing of transposon insertion sites from single tumor cells in a SB mouse model of myeloid leukemia (ML). SBCapSeq analysis of just 26 cells from one tumor revealed the tumor's major clonal subpopulations, enabled detection of clonal insertion events not detected by other sequencing methods and led to the identification of dominant subclones, each containing a unique pair of interacting gene drivers along with three to six cooperating cancer genes with SB-driven expression changes.

Sleeping Beauty screen reveals Pparg activation in metastatic prostate cancer.

Prostate cancer (CaP) is the most common adult male cancer in the developed world. The paucity of biomarkers to predict prostate tumor biology makes it important to identify key pathways that confer poor prognosis and guide potential targeted therapy. Using a murine forward mutagenesis screen in a Pten-null background, we identified peroxisome proliferator-activated receptor gamma (Pparg), encoding a ligand-activated transcription factor, as a promoter of metastatic CaP through activation of lipid signaling pathways, including up-regulation of lipid synthesis enzymes [fatty acid synthase (FASN), acetyl-CoA carboxylase (ACC), ATP citrate lyase (ACLY)]. Importantly, inhibition of PPARG suppressed tumor growth in vivo, with down-regulation of the lipid synthesis program. We show that elevated levels of PPARG strongly correlate with elevation of FASN in human CaP and that high levels of PPARG/FASN and PI3K/pAKT pathway activation confer a poor prognosis. These data suggest that CaP patients could be stratified in terms of PPARG/FASN and PTEN levels to identify patients with aggressive CaP who may respond favorably to PPARG/FASN inhibition.

A Genome-Wide Association Study for Regulators of Micronucleus Formation in Mice.

In mammals the regulation of genomic instability plays a key role in tumor suppression and also controls genome plasticity, which is important for recombination during the processes of immunity and meiosis. Most studies to identify regulators of genomic instability have been performed in cells in culture or in systems that report on gross rearrangements of the genome, yet subtle differences in the level of genomic instability can contribute to whole organism phenotypes such as tumor predisposition. Here we performed a genome-wide association study in a population of 1379 outbred Crl:CFW(SW)-US_P08 mice to dissect the genetic landscape of micronucleus formation, a biomarker of chromosomal breaks, whole chromosome loss, and extranuclear DNA. Variation in micronucleus levels is a complex trait with a genome-wide heritability of 53.1%. We identify seven loci influencing micronucleus formation (false discovery rate <5%), and define candidate genes at each locus. Intriguingly at several loci we find evidence for sexual dimorphism in micronucleus formation, with a locus on chromosome 11 being specific to males.

Sleeping Beauty transposon mutagenesis identifies genes that cooperate with mutant Smad4 in gastric cancer development.

Mutations in SMAD4 predispose to the development of gastrointestinal cancer, which is the third leading cause of cancer-related deaths. To identify genes driving gastric cancer (GC) development, we performed a Sleeping Beauty (SB) transposon mutagenesis screen in the stomach of Smad4(+/-) mutant mice. This screen identified 59 candidate GC trunk drivers and a much larger number of candidate GC progression genes. Strikingly, 22 SB-identified trunk drivers are known or candidate cancer genes, whereas four SB-identified trunk drivers, including PTEN, SMAD4, RNF43, and NF1, are known human GC trunk drivers. Similar to human GC, pathway analyses identified WNT, TGF-ß, and PI3K-PTEN signaling, ubiquitin-mediated proteolysis, adherens junctions, and RNA degradation in addition to genes involved in chromatin modification and organization as highly deregulated pathways in GC. Comparative oncogenomic filtering of the complete list of SB-identified genes showed that they are highly enriched for genes mutated in human GC and identified many candidate human GC genes. Finally, by comparing our complete list of SB-identified genes against the list of mutated genes identified in five large-scale human GC sequencing studies, we identified LDL receptor-related protein 1B (LRP1B) as a previously unidentified human candidate GC tumor suppressor gene. In LRP1B, 129 mutations were found in 462 human GC samples sequenced, and LRP1B is one of the top 10 most deleted genes identified in a panel of 3,312 human cancers. SB mutagenesis has, thus, helped to catalog the cooperative molecular mechanisms driving SMAD4-induced GC growth and discover genes with potential clinical importance in human GC.

Adenoma development in familial adenomatous polyposis and MUTYH-associated polyposis: somatic landscape and driver genes.

Familial adenomatous polyposis (FAP) and MUTYH-associated polyposis (MAP) are inherited disorders associated with multiple colorectal adenomas that lead to a very high risk of colorectal cancer. The somatic mutations that drive adenoma development in these conditions have not been investigated comprehensively. In this study we performed analysis of paired colorectal adenoma and normal tissue DNA from individuals with FAP or MAP, sequencing 14 adenoma whole exomes (eight MAP, six FAP), 55 adenoma targeted exomes (33 MAP, 22 FAP) and germline DNA from each patient, and a further 63 adenomas by capillary sequencing (41 FAP, 22 MAP). With these data we examined the profile of mutated genes, the mutational signatures and the somatic mutation rates, observing significant diversity in the constellations of mutated driver genes in different adenomas, and loss-of-function mutations in WTX (9%; p < 9.99e-06), a gene implicated in regulation of the WNT pathway and p53 acetylation. These data extend our understanding of the early events in colorectal tumourigenesis in the polyposis syndromes.

Somatic drivers of B-ALL in a model of ETV6-RUNX1; Pax5(+/-) leukemia.

BACKGROUND: B-cell precursor acute lymphoblastic leukemia (B-ALL) is amongst the leading causes of childhood cancer-related mortality. Its most common chromosomal aberration is the ETV6-RUNX1 fusion gene, with ~25% of ETV6-RUNX1 patients also carrying PAX5 alterations. METHODS: We have recreated this mutation background by inter-crossing Etv6-RUNX1 (Etv6 (RUNX1-SB)) and Pax5(+/-) mice and performed an in vivo analysis to find driver genes using Sleeping Beauty transposon-mediated mutagenesis and also exome sequencing. RESULTS: Combination of Etv6-RUNX1 and Pax5(+/-) alleles generated a transplantable B220 + CD19+ B-ALL with a significant disease incidence. RNA-seq analysis showed a gene expression pattern consistent with arrest at the pre-B stage. Analysis of the transposon common insertion sites identified genes involved in B-cell development (Zfp423) and the JAK/STAT signaling pathway (Jak1, Stat5 and Il2rb), while exome sequencing revealed somatic hotspot mutations in Jak1 and Jak3 at residues analogous to those mutated in human leukemias, and also mutation of Trp53. CONCLUSIONS: Powerful synergies exists in our model suggesting STAT pathway activation and mutation of Trp53 are potent drivers of B-ALL in the context of Etv6-RUNX1;Pax5(+/-).

A novel mouse model identifies cooperating mutations and therapeutic targets critical for chronic myeloid leukemia progression.

The introduction of highly selective ABL-tyrosine kinase inhibitors (TKIs) has revolutionized therapy for chronic myeloid leukemia (CML). However, TKIs are only efficacious in the chronic phase of the disease and effective therapies for TKI-refractory CML, or after progression to blast crisis (BC), are lacking. Whereas the chronic phase of CML is dependent on BCR-ABL, additional mutations are required for progression to BC. However, the identity of these mutations and the pathways they affect are poorly understood, hampering our ability to identify therapeutic targets and improve outcomes. Here, we describe a novel mouse model that allows identification of mechanisms of BC progression in an unbiased and tractable manner, using transposon-based insertional mutagenesis on the background of chronic phase CML. Our BC model is the first to faithfully recapitulate the phenotype, cellular and molecular biology of human CML progression. We report a heterogeneous and unique pattern of insertions identifying known and novel candidate genes and demonstrate that these pathways drive disease progression and provide potential targets for novel therapeutic strategies. Our model greatly informs the biology of CML progression and provides a potent resource for the development of candidate therapies to improve the dismal outcomes in this highly aggressive disease.

PAX5 is a tumor suppressor in mouse mutagenesis models of acute lymphoblastic leukemia.

Alterations of genes encoding transcriptional regulators of lymphoid development are a hallmark of B-progenitor acute lymphoblastic leukemia (B-ALL) and most commonly involve PAX5, encoding the DNA-binding transcription factor paired-box 5. The majority of PAX5 alterations in ALL are heterozygous, and key PAX5 target genes are expressed in leukemic cells, suggesting that PAX5 may be a haploinsufficient tumor suppressor. To examine the role of PAX5 alterations in leukemogenesis, we performed mutagenesis screens of mice heterozygous for a loss-of-function Pax5 allele. Both chemical and retroviral mutagenesis resulted in a significantly increased penetrance and reduced latency of leukemia, with a shift to B-lymphoid lineage. Genomic profiling identified a high frequency of secondary genomic mutations, deletions, and retroviral insertions targeting B-lymphoid development, including Pax5, and additional genes and pathways mutated in ALL, including tumor suppressors, Ras, and Janus kinase-signal transducer and activator of transcription signaling. These results show that in contrast to simple Pax5 haploinsufficiency, multiple sequential alterations targeting lymphoid development are central to leukemogenesis and contribute to the arrest in lymphoid maturation characteristic of ALL. This cross-species analysis also validates the importance of concomitant alterations of multiple cellular growth, signaling, and tumor suppression pathways in the pathogenesis of B-ALL.

Transposon mutagenesis identifies genetic drivers of Braf(V600E) melanoma.

Although nearly half of human melanomas harbor oncogenic BRAF(V600E) mutations, the genetic events that cooperate with these mutations to drive melanogenesis are still largely unknown. Here we show that Sleeping Beauty (SB) transposon-mediated mutagenesis drives melanoma progression in Braf(V600E) mutant mice and identify 1,232 recurrently mutated candidate cancer genes (CCGs) from 70 SB-driven melanomas. CCGs are enriched in Wnt, PI3K, MAPK and netrin signaling pathway components and are more highly connected to one another than predicted by chance, indicating that SB targets cooperative genetic networks in melanoma. Human orthologs of >500 CCGs are enriched for mutations in human melanoma or showed statistically significant clinical associations between RNA abundance and survival of patients with metastatic melanoma. We also functionally validate CEP350 as a new tumor-suppressor gene in human melanoma. SB mutagenesis has thus helped to catalog the cooperative molecular mechanisms driving BRAF(V600E) melanoma and discover new genes with potential clinical importance in human melanoma.

Transposon mutagenesis identifies genes and evolutionary forces driving gastrointestinal tract tumor progression.

To provide a more comprehensive understanding of the genes and evolutionary forces driving colorectal cancer (CRC) progression, we performed Sleeping Beauty (SB) transposon mutagenesis screens in mice carrying sensitizing mutations in genes that act at different stages of tumor progression. This approach allowed us to identify a set of genes that appear to be highly relevant for CRC and to provide a better understanding of the evolutionary forces and systems properties of CRC. We also identified six genes driving malignant tumor progression and a new human CRC tumor-suppressor gene, ZNF292, that might also function in other types of cancer. Our comprehensive CRC data set provides a resource with which to develop new therapies for treating CRC.

BRAF inhibitor resistance mediated by the AKT pathway in an oncogenic BRAF mouse melanoma model.

BRAF (v-raf murine sarcoma viral oncogene homolog B) inhibitors elicit a transient anti-tumor response in ∼ 80% of BRAF(V600)-mutant melanoma patients that almost uniformly precedes the emergence of resistance. Here we used a mouse model of melanoma in which melanocyte-specific expression of Braf(V618E) (analogous to the human BRAF(V600E) mutation) led to the development of skin hyperpigmentation and nevi, as well as melanoma formation with incomplete penetrance. Sleeping Beauty insertional mutagenesis in this model led to accelerated and fully penetrant melanomagenesis and synchronous tumor formation. Treatment of Braf(V618E) transposon mice with the BRAF inhibitor PLX4720 resulted in tumor regression followed by relapse. Analysis of transposon insertions identified eight genes including Braf, Mitf, and ERas (ES-cell expressed Ras) as candidate resistance genes. Expression of ERAS in human melanoma cell lines conferred resistance to PLX4720 and induced hyperphosphorylation of AKT (v-akt murine thymoma viral oncogene homolog 1), a phenotype reverted by combinatorial treatment with PLX4720 and the AKT inhibitor MK2206. We show that ERAS expression elicits a prosurvival signal associated with phosphorylation/inactivation of BAD, and that the resistance of hepatocyte growth factor-treated human melanoma cells to PLX4720 can be reverted by treatment with the BAD-like BH3 mimetic ABT-737. Thus, we define a role for the AKT/BAD pathway in resistance to BRAF inhibition and illustrate an in vivo approach for finding drug resistance genes.