In vivo identification of tumor- suppressive PTEN ceRNAs in an oncogenic BRAF-induced mouse model of melanoma.

We recently proposed that competitive endogenous RNAs (ceRNAs) sequester microRNAs to regulate mRNA transcripts containing common microRNA recognition elements (MREs). However, the functional role of ceRNAs in cancer remains unknown. Loss of PTEN, a tumor suppressor regulated by ceRNA activity, frequently occurs in melanoma. Here, we report the discovery of significant enrichment of putative PTEN ceRNAs among genes whose loss accelerates tumorigenesis following Sleeping Beauty insertional mutagenesis in a mouse model of melanoma. We validated several putative PTEN ceRNAs and further characterized one, the ZEB2 transcript. We show that ZEB2 modulates PTEN protein levels in a microRNA-dependent, protein coding-independent manner. Attenuation of ZEB2 expression activates the PI3K/AKT pathway, enhances cell transformation, and commonly occurs in human melanomas and other cancers expressing low PTEN levels. Our study genetically identifies multiple putative microRNA decoys for PTEN, validates ZEB2 mRNA as a bona fide PTEN ceRNA, and demonstrates that abrogated ZEB2 expression cooperates with BRAF(V600E) to promote melanomagenesis.

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.

Function of C/EBPdelta in a regulatory circuit that discriminates between transient and persistent TLR4-induced signals.

The innate immune system is like a double-edged sword: it is absolutely required for host defense against infection, but when uncontrolled, it can trigger a plethora of inflammatory diseases. Here we use systems-biology approaches to predict and confirm the existence of a gene-regulatory network involving dynamic interaction among the transcription factors NF-kappaB, C/EBPdelta and ATF3 that controls inflammatory responses. We mathematically modeled transcriptional regulation of the genes encoding interleukin 6 and C/EBPdelta and experimentally confirmed the prediction that the combination of an initiator (NF-kappaB), an amplifier (C/EBPdelta) and an attenuator (ATF3) forms a regulatory circuit that discriminates between transient and persistent Toll-like receptor 4-induced signals. Our results suggest a mechanism that enables the innate immune system to detect the duration of infection and to respond appropriately.

The mutational landscapes of genetic and chemical models of Kras-driven lung cancer.

Next-generation sequencing of human tumours has refined our understanding of the mutational processes operative in cancer initiation and progression, yet major questions remain regarding the factors that induce driver mutations and the processes that shape mutation selection during tumorigenesis. Here we performed whole-exome sequencing on adenomas from three mouse models of non-small-cell lung cancer, which were induced either by exposure to carcinogens (methyl-nitrosourea (MNU) and urethane) or by genetic activation of Kras (Kras(LA2)). Although the MNU-induced tumours carried exactly the same initiating mutation in Kras as seen in the Kras(LA2) model (G12D), MNU tumours had an average of 192 non-synonymous, somatic single-nucleotide variants, compared with only six in tumours from the Kras(LA2) model. By contrast, the Kras(LA2) tumours exhibited a significantly higher level of aneuploidy and copy number alterations compared with the carcinogen-induced tumours, suggesting that carcinogen-induced and genetically engineered models lead to tumour development through different routes. The wild-type allele of Kras has been shown to act as a tumour suppressor in mouse models of non-small-cell lung cancer. We demonstrate that urethane-induced tumours from wild-type mice carry mostly (94%) Kras Q61R mutations, whereas those from Kras heterozygous animals carry mostly (92%) Kras Q61L mutations, indicating a major role for germline Kras status in mutation selection during initiation. The exome-wide mutation spectra in carcinogen-induced tumours overwhelmingly display signatures of the initiating carcinogen, while adenocarcinomas acquire additional C > T mutations at CpG sites. These data provide a basis for understanding results from human tumour genome sequencing, which has identified two broad categories of tumours based on the relative frequency of single-nucleotide variations and copy number alterations, and underline the importance of carcinogen models for understanding the complex mutation spectra seen in human cancers.

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.

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.

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 FOXO3-IRF7 gene regulatory circuit limits inflammatory sequelae of antiviral responses.

Antiviral responses must be tightly regulated to defend rapidly against infection while minimizing inflammatory damage. Type 1 interferons (IFN-I) are crucial mediators of antiviral responses and their transcription is regulated by a variety of transcription factors; principal among these is the family of interferon regulatory factors (IRFs). The IRF gene regulatory networks are complex and contain multiple feedback loops. The tools of systems biology are well suited to elucidate the complex interactions that give rise to precise coordination of the interferon response. Here we have used an unbiased systems approach to predict that a member of the forkhead family of transcription factors, FOXO3, is a negative regulator of a subset of antiviral genes. This prediction was validated using macrophages isolated from Foxo3-null mice. Genome-wide location analysis combined with gene deletion studies identified the Irf7 gene as a critical target of FOXO3. FOXO3 was identified as a negative regulator of Irf7 transcription and we have further demonstrated that FOXO3, IRF7 and IFN-I form a coherent feed-forward regulatory circuit. Our data suggest that the FOXO3-IRF7 regulatory circuit represents a novel mechanism for establishing the requisite set points in the interferon pathway that balances the beneficial effects and deleterious sequelae of the antiviral response.

The deubiquitinase USP9X suppresses pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDA) remains a lethal malignancy despite much progress concerning its molecular characterization. PDA tumours harbour four signature somatic mutations in addition to numerous lower frequency genetic events of uncertain significance. Here we use Sleeping Beauty (SB) transposon-mediated insertional mutagenesis in a mouse model of pancreatic ductal preneoplasia to identify genes that cooperate with oncogenic Kras(G12D) to accelerate tumorigenesis and promote progression. Our screen revealed new candidate genes for PDA and confirmed the importance of many genes and pathways previously implicated in human PDA. The most commonly mutated gene was the X-linked deubiquitinase Usp9x, which was inactivated in over 50% of the tumours. Although previous work had attributed a pro-survival role to USP9X in human neoplasia, we found instead that loss of Usp9x enhances transformation and protects pancreatic cancer cells from anoikis. Clinically, low USP9X protein and messenger RNA expression in PDA correlates with poor survival after surgery, and USP9X levels are inversely associated with metastatic burden in advanced disease. Furthermore, chromatin modulation with trichostatin A or 5-aza-2'-deoxycytidine elevates USP9X expression in human PDA cell lines, indicating a clinical approach for certain patients. The conditional deletion of Usp9x cooperated with Kras(G12D) to accelerate pancreatic tumorigenesis in mice, validating their genetic interaction. We propose that USP9X is a major tumour suppressor gene with prognostic and therapeutic relevance in PDA.

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.

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.

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.

Chromatin landscapes of retroviral and transposon integration profiles.

The ability of retroviruses and transposons to insert their genetic material into host DNA makes them widely used tools in molecular biology, cancer research and gene therapy. However, these systems have biases that may strongly affect research outcomes. To address this issue, we generated very large datasets consisting of ~ 120,000 to ~ 180,000 unselected integrations in the mouse genome for the Sleeping Beauty (SB) and piggyBac (PB) transposons, and the Mouse Mammary Tumor Virus (MMTV). We analyzed ~ 80 (epi)genomic features to generate bias maps at both local and genome-wide scales. MMTV showed a remarkably uniform distribution of integrations across the genome. More distinct preferences were observed for the two transposons, with PB showing remarkable resemblance to bias profiles of the Murine Leukemia Virus. Furthermore, we present a model where target site selection is directed at multiple scales. At a large scale, target site selection is similar across systems, and defined by domain-oriented features, namely expression of proximal genes, proximity to CpG islands and to genic features, chromatin compaction and replication timing. Notable differences between the systems are mainly observed at smaller scales, and are directed by a diverse range of features. To study the effect of these biases on integration sites occupied under selective pressure, we turned to insertional mutagenesis (IM) screens. In IM screens, putative cancer genes are identified by finding frequently targeted genomic regions, or Common Integration Sites (CISs). Within three recently completed IM screens, we identified 7%-33% putative false positive CISs, which are likely not the result of the oncogenic selection process. Moreover, results indicate that PB, compared to SB, is more suited to tag oncogenes.

Insertional mutagenesis and deep profiling reveals gene hierarchies and a Myc/p53-dependent bottleneck in lymphomagenesis.

Retroviral insertional mutagenesis (RIM) is a powerful tool for cancer genomics that was combined in this study with deep sequencing (RIM/DS) to facilitate a comprehensive analysis of lymphoma progression. Transgenic mice expressing two potent collaborating oncogenes in the germ line (CD2-MYC, -Runx2) develop rapid onset tumours that can be accelerated and rendered polyclonal by neonatal Moloney murine leukaemia virus (MoMLV) infection. RIM/DS analysis of 28 polyclonal lymphomas identified 771 common insertion sites (CISs) defining a 'progression network' that encompassed a remarkably large fraction of known MoMLV target genes, with further strong indications of oncogenic selection above the background of MoMLV integration preference. Progression driven by RIM was characterised as a Darwinian process of clonal competition engaging proliferation control networks downstream of cytokine and T-cell receptor signalling. Enhancer mode activation accounted for the most efficiently selected CIS target genes, including Ccr7 as the most prominent of a set of chemokine receptors driving paracrine growth stimulation and lymphoma dissemination. Another large target gene subset including candidate tumour suppressors was disrupted by intragenic insertions. A second RIM/DS screen comparing lymphomas of wild-type and parental transgenics showed that CD2-MYC tumours are virtually dependent on activation of Runx family genes in strong preference to other potent Myc collaborating genes (Gfi1, Notch1). Ikzf1 was identified as a novel collaborating gene for Runx2 and illustrated the interface between integration preference and oncogenic selection. Lymphoma target genes for MoMLV can be classified into (a) a small set of master regulators that confer self-renewal; overcoming p53 and other failsafe pathways and (b) a large group of progression genes that control autonomous proliferation in transformed cells. These findings provide insights into retroviral biology, human cancer genetics and the safety of vector-mediated gene therapy.

Nuclear receptor binding protein 1 regulates intestinal progenitor cell homeostasis and tumour formation.

Genetic screens in simple model organisms have identified many of the key components of the conserved signal transduction pathways that are oncogenic when misregulated. Here, we identify H37N21.1 as a gene that regulates vulval induction in let-60(n1046gf), a strain with a gain-of-function mutation in the Caenorhabditis elegans Ras orthologue, and show that somatic deletion of Nrbp1, the mouse orthologue of this gene, results in an intestinal progenitor cell phenotype that leads to profound changes in the proliferation and differentiation of all intestinal cell lineages. We show that Nrbp1 interacts with key components of the ubiquitination machinery and that loss of Nrbp1 in the intestine results in the accumulation of Sall4, a key mediator of stem cell fate, and of Tsc22d2. We also reveal that somatic loss of Nrbp1 results in tumourigenesis, with haematological and intestinal tumours predominating, and that nuclear receptor binding protein 1 (NRBP1) is downregulated in a range of human tumours, where low expression correlates with a poor prognosis. Thus NRBP1 is a conserved regulator of cell fate, that plays an important role in tumour suppression.

Sleeping Beauty mutagenesis in a mouse medulloblastoma model defines networks that discriminate between human molecular subgroups.

The Sleeping Beauty (SB) transposon mutagenesis screen is a powerful tool to facilitate the discovery of cancer genes that drive tumorigenesis in mouse models. In this study, we sought to identify genes that functionally cooperate with sonic hedgehog signaling to initiate medulloblastoma (MB), a tumor of the cerebellum. By combining SB mutagenesis with Patched1 heterozygous mice (Ptch1(lacZ/+)), we observed an increased frequency of MB and decreased tumor-free survival compared with Ptch1(lacZ/+) controls. From an analysis of 85 tumors, we identified 77 common insertion sites that map to 56 genes potentially driving increased tumorigenesis. The common insertion site genes identified in the mutagenesis screen were mapped to human orthologs, which were used to select probes and corresponding expression data from an independent set of previously described human MB samples, and surprisingly were capable of accurately clustering known molecular subgroups of MB, thereby defining common regulatory networks underlying all forms of MB irrespective of subgroup. We performed a network analysis to discover the likely mechanisms of action of subnetworks and used an in vivo model to confirm a role for a highly ranked candidate gene, Nfia, in promoting MB formation. Our analysis implicates candidate cancer genes in the deregulation of apoptosis and translational elongation, and reveals a strong signature of transcriptional regulation that will have broad impact on expression programs in MB. These networks provide functional insights into the complex biology of human MB and identify potential avenues for intervention common to all clinical subgroups.

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(+/-).

Transposon mutagenesis identifies genes that transform neural stem cells into glioma-initiating cells.

Neural stem cells (NSCs) are considered to be the cell of origin of glioblastoma multiforme (GBM). However, the genetic alterations that transform NSCs into glioma-initiating cells remain elusive. Using a unique transposon mutagenesis strategy that mutagenizes NSCs in culture, followed by additional rounds of mutagenesis to generate tumors in vivo, we have identified genes and signaling pathways that can transform NSCs into glioma-initiating cells. Mobilization of Sleeping Beauty transposons in NSCs induced the immortalization of astroglial-like cells, which were then able to generate tumors with characteristics of the mesenchymal subtype of GBM on transplantation, consistent with a potential astroglial origin for mesenchymal GBM. Sequence analysis of transposon insertion sites from tumors and immortalized cells identified more than 200 frequently mutated genes, including human GBM-associated genes, such as Met and Nf1, and made it possible to discriminate between genes that function during astroglial immortalization vs. later stages of tumor development. We also functionally validated five GBM candidate genes using a previously undescribed high-throughput method. Finally, we show that even clonally related tumors derived from the same immortalized line have acquired distinct combinations of genetic alterations during tumor development, suggesting that tumor formation in this model system involves competition among genetically variant cells, which is similar to the Darwinian evolutionary processes now thought to generate many human cancers. This mutagenesis strategy is faster and simpler than conventional transposon screens and can potentially be applied to any tissue stem/progenitor cells that can be grown and differentiated in vitro.

Sleeping Beauty mutagenesis reveals cooperating mutations and pathways in pancreatic adenocarcinoma.

Pancreatic cancer is one of the most deadly cancers affecting the Western world. Because the disease is highly metastatic and difficult to diagnosis until late stages, the 5-y survival rate is around 5%. The identification of molecular cancer drivers is critical for furthering our understanding of the disease and development of improved diagnostic tools and therapeutics. We have conducted a mutagenic screen using Sleeping Beauty (SB) in mice to identify new candidate cancer genes in pancreatic cancer. By combining SB with an oncogenic Kras allele, we observed highly metastatic pancreatic adenocarcinomas. Using two independent statistical methods to identify loci commonly mutated by SB in these tumors, we identified 681 loci that comprise 543 candidate cancer genes (CCGs); 75 of these CCGs, including Mll3 and Ptk2, have known mutations in human pancreatic cancer. We identified point mutations in human pancreatic patient samples for another 11 CCGs, including Acvr2a and Map2k4. Importantly, 10% of the CCGs are involved in chromatin remodeling, including Arid4b, Kdm6a, and Nsd3, and all SB tumors have at least one mutated gene involved in this process; 20 CCGs, including Ctnnd1, Fbxo11, and Vgll4, are also significantly associated with poor patient survival. SB mutagenesis provides a rich resource of mutations in potential cancer drivers for cross-comparative analyses with ongoing sequencing efforts in human pancreatic adenocarcinoma.

Cake: a bioinformatics pipeline for the integrated analysis of somatic variants in cancer genomes.

We have developed Cake, a bioinformatics software pipeline that integrates four publicly available somatic variant-calling algorithms to identify single nucleotide variants with higher sensitivity and accuracy than any one algorithm alone. Cake can be run on a high-performance computer cluster or used as a stand-alone application. Availabilty: Cake is open-source and is available from http://cakesomatic.sourceforge.net/