Skip to main content

Advertisement

Case-oriented pathways analysis in pancreatic adenocarcinoma using data from a sleeping beauty transposon mutagenesis screen

Article metrics

Abstract

Background

Mutation studies of pancreatic ductal adenocarcinoma (PDA) have revealed complicated heterogeneous genomic landscapes of the disease. These studies cataloged a number of genes mutated at high frequencies, but also report a very large number of genes mutated in lower percentages of tumors. Taking advantage of a well-established forward genetic screening technique, with the Sleeping Beauty (SB) transposon, several studies produced PDA and discovered a number of common insertion sites (CIS) and associated genes that are recurrently mutated at high frequencies. As with human mutation studies, a very large number of genes were found to be altered by transposon insertion at low frequencies. These low frequency CIS associated genes may be very valuable to consider for their roles in cancer, since collectively they might emerge from a core group of genetic pathways.

Result

In this paper, we determined whether the genetic mutations in SB-accelerated PDA occur within a collated group of biological processes defined as gene sets. The approach considered both genes mutated in high and lower frequencies. We implemented a case-oriented, gene set enrichment analysis (CO-GSEA) on SB altered genes in PDA. Compared to traditional GSEA, CO-GSEA enables us to consider individual characteristics of mutation profiles of each PDA tumor. We identified genetic pathways with higher numbers of genetic mutations than expected by chance. We also present the correlations between these significant enriched genetic pathways, and their associations with CIS genes.

Conclusion

These data suggest that certain pathway alterations cooperate in PDA development.

Background

The molecular analysis of human cancer cells has revealed a startling amount of genetic and epigenetic heterogeneity. In recent years, forward genetic screens have taken place in mice using DNA transposons, primarily Sleeping Beauty (SB) [1]. The SB-based approach has been successfully employed to induce many different forms of cancer such as brain tumors, sarcomas, hematopoietic malignancies, and carcinomas [2, 3] via insertional mutagenesis. A large number of of loci recurrently mutated by insertion of SB transposons called common insertion sites (CIS) have been identified [4]. The general impression from these studies is one of tremendous genetic complexity.

Recent large-scale analyses of human cancer genomes mirrors these results in general. Most types of human cancer harbor a small number of genes that are altered in a high percentage of cases, so called “mountains”, and a large number of genes altered in a low percentage of cases, so called “hills”. In addition, two patients diagnosed with the same type of cancer often show distinct genetic alternations, however, the disrupted pathways tend to be similar among patients [5].

Conventional pathway analysis approaches usually obtain gene-based scores by summarizing data across tumor cases, then calculating pathway statistics using the scores of the genes in the pathway. However, these approaches could potentially lose information regarding whether multiple mutations in a pathway are from a single patient or multiple cases with a single mutation at various genes in the pathway. In contrast, case-oriented gene set analysis (CO-GSEA) can consider the two situations differently and hence can incorporate heterogeneity of each tumor case into the analysis. This approach provides a case-based score for each pathway and further enhances the study of correlation of mutation events between pathways, as well as between genes and pathways. It has recently been applied to the analysis of human tumors [6].

Pancreatic ductal adenocarcinoma (PDA) is the fourth leading cause of death due to cancer, with over a 98 % case-fatality rate. The crucial molecular events, required for progression from a pre-invasive and non-life threatening state to an invasive and metastatic lethal condition, are not well-understood. We previously reported the results of a SB transposon-based forward genetic screen for drivers of PDA in mice expressing the KrasG12D oncogene in epithelial cells of the pancreas [7]. 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 chromosome-linked deubiquitinase Usp9x, which was inactivated in over 50 % of the tumors. In addition, several hundred candidate PDA genes were identified as CIS in this screen.

In this paper, we report analyses intended to determine whether a core group biological processes or pathways are populated by genes from CIS. We applied a less stringent criterion to consider CIS associated genes that mutated both at high frequencies (mountains) and at lower frequencies (hills). Secondly, we determined whether non-random associations between alteration of genes in certain pathways or biological pathways exist by analysis of CIS from individual tumors.

Results and discussion

Certain pathways are enriched in CIS-associated genes

We collected insertional mutatgenesis data of tumor samples from 146 KrasLSL-G12D; Pdx1-cre; T2/Onc; Rosa26-LSL-SB13 mice. To determine whether a core group of pathways were enriched with CIS-associated genes than reported previously [7], we analyszed 968 CIS with uncorrected p value <10−4 from TapDance. Among these, 239 genes were mapped an grouped into 281 KEGG curated pathways categories. After excluding pathways with less than 6 genes, 272 KEGG pathways remain in the following analysis.

Using the CO-GSEA described in the Method Section, we found 95 KEGG pathways that are enriched with CIS-associated genes with permutated p value <10−7 listed in Table 1 (more details about the disrupted genes in each pathway can be found in Additional file 1). In Table 1, “# of genes” records the number of genes defined in the pathway from KEGG; “# of CIS” (third column) reports the number of CIS genes in the pathway; and “# of mutated cases” (fourth column) records the number of cases that the pathway was disrupted. A histogram of the sizes of each of the KEGG pathways is shown in Additional file 2. In Figs. 1 and 2, we plotted the KEGG diagrams of two pathways that are enriched CIS-associated genes.

Fig. 1
figure1

Frequently mutated genes in ubiquitin mediated proteolysis pathway. Darker red color indicates higher mutation frequencies in mice

Fig. 2
figure2

Frequently mutated genes in ErbB signaling pathway. Darker color indicates higher mutation frequencies in mice

Table 1 Pathways that are enriched with CIS-associated genes (permuted p value <10−7)

The genetic screen was designed to discover genes that when altered would cause acceleration of PDA in pancreatic ductal epithelial cells expressing an activated form of the Kras oncogene, KrasG12D. As such it was not surprising that KEGG pathways with the strongest statistical support for CIS associated gene enrichment were many cancer associated pathways. As expected, we found some of the same pathways previously reported and which were expected [7, 8]. An informal prior analysis [7] suggested that TGF β signaling was enriched in CIS-associated genes and indeed we found that this KEGG pathway is enriched. Similarly, Rb1/p16Inka4a pathway was suggested to be recurrently altered by CIS-associated genes [7]. Indeed, we found that the KEGG pathway CELL CYCLE was enriched in CIS-associated genes. Many other cancer-associated pathways were enriched in CIS-associated genes including the RAS, PI3K-AKT, HIPPO, VEGF, HEDGEHOG, MAPK, FOXO1, and MTOR pathways. Moreover, the human disease KEGG pathway PANCREATIC CANCER and several other human cancer pathways were enriched in CIS-associated genes.

In addition to these expected KEGG pathways, many involving metabolism have not been strongly linked to pancreatic cancer development or cancer development in general. However, recent studies revealed evidence of metabolic reprogramming to sustain tumor survival in KRAS-mutated PDA tumors [9]. For example, KRAS-dependent tumor cells compensated the energy loss through increasing glycolysis, amino acid and lipid biosynthesis [10]. In particular, TERPENOID BIOSYNTHESIS, LYSINE DEGRADATION and the SULFUR RELAY SYSTEM are significantly altered in the SB-accelerated tumor models. To date KRAS remains a poorly druggable target, hence, targeting the downstream metabolic regulation could be effective alternatives in inhibiting tumor growth.

Several organismal systems KEGG pathways were also enriched in CIS-associated genes despite not being strongly linked to pancreatic cancer development. These include OXYTOCIN SIGNALING, CHOLINERGIC SYNAPSE, and MELANOGENESIS. Our recent work helped show that the AXON GUIDANCE pathway is enriched for CIS-associated genes, a result which led to the discovery that these genes and the pathways they participate in are altered in human PDA [8]. This result was reproduced in this current analysis. Thus, it is clear that the broadened definition of CIS allows for the identification of many known and novel candidate cancer pathways. These data suggest many new hypotheses to be tested in PDA development.

Analysis of individual tumor reveals significant co-altered pathways

We and others have published results of SB screens in which we found that individual CIS tended to be co-mutated by transposon insertion more than expected by chance (e.g. [11]). We wondered whether an analysis of individual tumors would reveal that specific pathways would be co-altered in this same manner. Figure 3 shows a heat map of adjusted correlation between pair of pathways, which are co-altered by transposon insertions within/near genes in those pathways. We observed that there are two major clusters of strongly co-altered pathways. Within these clusters certain specific pathways show strong associations, being altered by transposon insertion in the same tumors more often than would be expected by chance. These data provide the basis for developing specific hypotheses about pathways that interact to cause cancer. Thus, alterations of one pathway may allow the other pathway to exert its full oncogenic effects.

Fig. 3
figure3

Heat map of correlation between pair of pathways. Legend indicates strength of correlation coefficient (red: high correlation; black: weak correlation). a Heat map of correlation between all pairs of pathways. b Zoom in of block 1 shown in panel a. c Zoom in of block 2 shown in panel a. Pathway names from left to right (bottom to top): 1: Ubiquitin mediated proteolysis, 2: One carbon pool by folate, 3: Wnt signaling pathway, 4: Cell cycle, 5: Hippo signaling pathway, 6: Protein processing in endoplasmic reticulum, 7: MAPK signaling pathway, 8: Lysine degradation, 9: RNA polymerase, 10: N-Glycan biosynthesis, 11: mRNA surveillance pathway, 12: Hedgehog signaling pathway, 13: Dopaminergic synapse, 14: GPI-anchor, 15: RNA degradation, 16: Terpenoid backbone biosynthesis, 17: Mucin type O-Glycan biosynthesis, 18: Rap1 signaling pathway, 19: Adherens junction, 20: Leukocyte transendothelial migration, 21: TGF-beta signaling pathway, 22: Axon guidance, 23: MicroRNAs in cancer, 24: Tight junction, 25: PI3K-Akt signaling pathway, 26: Ras signaling pathway, 27: Chemokine signaling pathway, 28: Serotonergic synapse, 29: Glutamatergic synapse, 30: Cholinergic synapse, 31: Retrograde endocannabinoid signaling, 32: Circadian entrainment, 33: Sulfur relay system, 34: Notch signaling pathway, 35: Alanine, aspartate and glutamate metabolism, 36: Aldosterone-regulated sodium reabsorption, 37: Natural killer cell mediated cytotoxicity, 38: VEGF signaling pathway, 39: Fc epsilon RI signaling pathway, 40: Gap junction, 41: Melanogenesis, 42: FoxO signaling pathway, 43: HIF-1 signaling pathway, 44: Fc gamma R-mediated phagocytosis, 45: Estrogen signaling pathway, 46: Platelet activation, 47: Oxytocin signaling pathway, 48: Vascular smooth muscle contraction, 49: PIP, 50: mTOR signaling pathway, 51: Focal adhesion, 52: Regulation of actin cytoskeleton, 53: Insulin signaling pathway, 54: Thyroid hormone signaling pathway, 55: ErbB signaling pathway, 56: T cell receptor signaling pathway, 57: Neurotrophin signaling pathway, 58: B cell receptor signaling pathway, 59: Prolactin signaling pathway

A careful analysis of some of the associations reveals pairs of pathways that might be predicted to interact based on what is known about their functions and regulation already. For example, block 1, labeled in Fig. 3, contains strong associations between the ubiquitin processing pathway and several pathways including ErbB, Insulin and mTOR signaling. It is known that cell signaling pathways that transmit signals from the extracellular space into the cell cytoplasm and nucleus are regulated by the abundance and stability of certain proteins. In many cases, the stability of these proteins is regulated by ubiquitination and degradation by the proteosome. Well known examples, include NF κB and Wnt/ β-catenin signaling pathways. Work shows that members of the ErbB family of receptors are downregulated by ubiquitination involving the E3 ubiquitin ligase Cbl [12]. Ubiquitination also regulates Akt-mTOR signaling in multiple myeloma [13] and Akt-mTOR is activated by insulinsignaling [14].

Block 2, labeled in Fig. 3 contains several other intriguing pathway-pathway associations. For example, we see a strong association between cell cycle control and miRNAs known to be involved in cancer. Indeed, there are several well studied examples of miRNAs that regulate the mRNA transcripts of cell cycle regulators such as MYC [15], RB1 [16] and CCND1 [17]. Also in block 2, we see evidence for TGF β pathway and MAPK pathway co-dysregulation. Abundant evidence for crosstalk between these pathways has been published [18, 19]. Thus, it is entirely plausible that co-alteration between these pathways is specifically selected for during PDA progression. Specific hypotheses can, or have been, tested in the laboratory. For example, MAPK activation, via expression of the KrasG12D oncogene, cooperates strongly with Smad4 inactivation, which alters/inactivates TGF β signaling, in a mouse model of PDA [20]. This functionally confirms the observation from the analyses done here. We can thus predict, that many other pathway-pathway associations observed in Fig. 3 can be functionally validated. More speculative, but of tremendous therapeutic significance, is the idea that targeting one pathway of a pathway-pathway pair observed in Fig. 3 would alter the ability of the second pathway to exert its oncogenic effects. Indeed, co-targeting both of such pairs of altered pathways may be the most effective way to treat individual cases of PDA. These ideas wait functional testing in the laboratory using model systems.

Association of CIS-associated genes and enriched pathways

Several of the most commonly altered genes in the PDA screen (i.e. the top ranked CIS-associated genes) have little published functional data. We speculate that by finding which pathways they most often interact with, something could be learned about their function in general and in PDA development. The associations between the top ranked CIS and enriched pathways are shown in Fig. 4. In Fig. 4, several CIS-associated genes such as Stag2, Arhgap5, Usp9x, Magi1, Arid1a have few connections to enriched pathways then other CIS-associated genes. In Additional file 3, we listed these connections and corresponding estimates from regression model, p values and FDR. For example, in Additional file 1, Usp9x is associated with PI3K-AKT signaling pathway, DOPAMINERGICSYNAPSE, HIPPO signaling, and TIGHT JUNCTION pathways. Thus, it seems likely that Usp9x mutation or down regulation has to cooperate with alterations in these other pathways in order for PDA to develop. The CIS-associated genes that also demonstrated association with the these Usp9x-associated pathways are Gsk3b, Ctnna1, Mll5, Pten, Arfip1, Magil.

Fig. 4
figure4

Associations between CIS and enriched pathways. Red nodes represent CIS associated genes, and blue nodes indicate pathways. Abbreviations: miRNAs: MicroRNAs in cancer, GPI: Glycosylphosphatidylinositol-anchor biosynthesis, T.J.: Tight junction, Hedgehog: Hedgehog signaling pathway, PIP: Phosphatidylinositol signaling system, B cell: B cell receptor signaling pathway, Leuk. M.: Leukocyte transendothelial migration, DA: Dopaminergic synapse, Axon: Axon guidance, Ubiquitin: Ubiquitin mediated proteolysis, F.A: Focal adhesion, M.G.: Melanogenesis, Lysine D.: Lysine degradation, Protein: Protein processing in endoplasmic reticulum, T cell: T cell receptor signaling pathway, T.H.: Thyroid hormone signaling pathway, N.T: Neurotrophin signaling pathway, A.J.: Adherens junction, Actin: Regulation of actin cytoskeleton, Ck: Chemokine signaling pathway

Conclusion

In this work, we demonstrate the non-random enrichment of CIS-associated genes from a transposon-based screen for PDA into certain KEGG signaling pathways, disease states and biological processes.

Methods

To assess whether a pathway harbors more CIS-associated genes than expected by chance, we use CO-GSEA approach. For each tumor sample, we considered a pathway is altered (coded: 1) if at least 1 gene in the pathway was mutated; coded zero if it’s not. A score for each pathway was calculated to be the number of tumors in which the pathway is altered. We assessed whether the score of a pathway was statistically significant through random permutation. For example, if a mouse tumor contains 100 mutations, we randomly assigned the 100 mutations to 100 different genes. A score for each pathway can be obtained by counting the number of altered tumor samples after the permutation. We repeated the permutation 107 times to obtain the distribution of score under the null for each pathway and calculated a p value based on the permuted null distribution. A similar approach was also applied in mutation analysis of human tumor samples in [6].

The ability to detect a significant pathway using the CO-GSEA approach depends on the background mutation rate and the size of the pathway under consideration. The relationship between the number of total cases, and the expected score of a given pathway under random permutation can be described as: \(N- \Sigma _{i=1}^{N}\tiny {\frac {\left (\! \begin {array}{c} G-n_{i} \\ P_{s} \end {array} \!\right)} {\left (\! \begin {array}{c} G\\ P_{s} \end {array} \!\right)}}\), where N is the total number of cases; G is the number of genes considered in the pathway analysis; n i is the number of events in sample i and P s is the number of genes in the pathway [6].

Analysis of co-altered pathways

To investigate whether a pair of pathways was co-altered in a significant manner, we remove the CIS-associated genes that are present in both pathways, and for each sample, we calculated the mutation frequency in each pathway using the remaining non-overlapping CIS as: \(\frac {\text {\# of mutations in sample } i\text { in the pathway}}{\text {\# of non-overlapping CIS in the pathway}}\). For each pair of pathways, Pearson correlations were calculated to present the correlation between pathways characterized by non-overlapping CIS using the mutation counts.

Association between top CIS-associated genes and enriched pathways

Among the top 20 CIS-associated genes previously reported [7], 12 of them listed do not map to any KEGG pathways. We conducted association analysis between the top 20 CIS-associated genes and the enriched KEGG pathway using quasi-Poisson regression models with over-dispersion. For each CIS-associated gene, we examined whether mutation status of the CIS-associated gene (code 1 if mutated; 0 otherwise) is associated with higher mutation counts (the number of altered CIS-associated genes) for a pathway under consideration. We reported the CIS and pathways associations with FDR < 0.001.

References

  1. 1

    Mann KM, Ward JM, Yew CCK, Kovochich A, Dawson DW, Black MA, Brett BT, Sheetz TE, Dupuy AJ, Australian Pancreatic Cancer Genome Initiative, Chang DK, Biankin AV, Waddell N, Kassahn KS, Grimmond SM, Rust AG, Adams DJ, Jenkins NA, Copeland NG. Sleeping beauty mutagenesis reveals cooperating mutations and pathways in pancreatic adenocarcinoma. Proc Natl Acad Sci U S A. 2012; 109(16):5934–41. doi:10.1073/pnas.1202490109.

  2. 2

    Copeland NG, Jenkins NA. Harnessing transposons for cancer gene discovery. Nat Rev Cancer. 2010; 10(10):696–706. doi:10.1038/nrc2916.

  3. 3

    O’Donnell KA, Keng VW, York B, Reineke EL, Seo D, Fan D, Silverstein KAT, Schrum CT, Xie WR, Mularoni L, Wheelan SJ, Torbenson MS, O’Malley BW, Largaespada DA, Boeke JD. A sleeping beauty mutagenesis screen reveals a tumor suppressor role for ncoa2/src-2 in liver cancer. Proc Natl Acad Sci U S A. 2012; 109(21):1377–86. doi:10.1073/pnas.1115433109.

  4. 4

    Abbott KL, Nyre ET, Abrahante J, Ho YY, Isaksson Vogel R, Starr TK. The candidate cancer gene database: a database of cancer driver genes from forward genetic screens in mice. Nucleic Acids Res. 2015; 43(Database issue):844–8. doi:10.1093/nar/gku770.

  5. 5

    Vogelstein B, Papadopoulos N, Velculescu VE, Zhou S, Diaz LA Jr, Kinzler KW. Cancer genome landscapes. Science. 2013; 339(6127):1546–58. doi:10.1126/science.1235122.

  6. 6

    Boca SM, Kinzler KW, Velculescu VE, Vogelstein B, Parmigiani G. Patient-oriented gene set analysis for cancer mutation data. Genome Biol. 2010; 11(11):112. doi:10.1186/gb-2010-11-11-r112.

  7. 7

    Pérez-Mancera PA, Rust AG, van der Weyden L, Kristiansen G, Li A, Sarver AL, Silverstein KAT, Grützmann R, Aust D, Rümmele P, Knösel T, Herd C, Stemple DL, Kettleborough R, Brosnan JA, Li A, Morgan R, Knight S, Yu J, Stegeman S, Collier LS, ten Hoeve JJ, de Ridder J, Klein AP, Goggins M, Hruban RH, Chang DK, Biankin AV, Grimmond SM, Australian Pancreatic Cancer Genome Initiative, Wessels LFA, Wood SA, Iacobuzio-Donahue CA, Pilarsky C, Largaespada DA, Adams DJ, Tuveson DA. The deubiquitinase usp9x suppresses pancreatic ductal adenocarcinoma. Nature. 2012; 486(7402):266–70. doi:10.1038/nature11114.

  8. 8

    Biankin AV, Waddell N, Kassahn KS, Gingras MC, Muthuswamy LB, Johns AL, Miller DK, Wilson PJ, Patch AM, Wu J, Chang DK, Cowley MJ, Gardiner BB, Song S, Harliwong I, Idrisoglu S, Nourse C, Nourbakhsh E, Manning S, Wani S, Gongora M, Pajic M, Scarlett CJ, Gill AJ, Pinho AV, Rooman I, Anderson M, Holmes O, Leonard C, Taylor D, Wood S, Xu Q, Nones K, Fink JL, Christ A, Bruxner T, Cloonan N, Kolle G, Newell F, Pinese M, Mead RS, Humphris JL, Kaplan W, Jones MD, Colvin EK, Nagrial AM, Humphrey ES, Chou A, Chin VT, Chantrill LA, Mawson A, Samra JS, Kench JG, Lovell JA, Daly RJ, Merrett ND, Toon C, Epari K, Nguyen NQ, Barbour A, Zeps N, Australian Pancreatic Cancer Genome Initiative, Kakkar N, Zhao F, Wu YQ, Wang M, Muzny DM, Fisher WE, Brunicardi FC, Hodges SE, Reid JG, Drummond J, Chang K, Han Y, Lewis LR, Dinh H, Buhay CJ, Beck T, Timms L, Sam M, Begley K, Brown A, Pai D, Panchal A, Buchner N, De Borja R, Denroche RE, Yung CK, Serra S, Onetto N, Mukhopadhyay D, Tsao MS, Shaw PA, Petersen GM, Gallinger S, Hruban RH, Maitra A, Iacobuzio-Donahue CA, Schulick RD, Wolfgang CL, Morgan RA, Lawlor RT, Capelli P, Corbo V, Scardoni M, Tortora G, Tempero MA, Mann KM, Jenkins NA, Perez-Mancera PA, Adams DJ, Largaespada DA, Wessels LFA, Rust AG, Stein LD, Tuveson DA, Copeland NG, Musgrove EA, Scarpa A, Eshleman JR, Hudson TJ, Sutherland RL, Wheeler DA, Pearson JV, McPherson JD, Gibbs RA, Grimmond SM. Pancreatic cancer genomes reveal aberrations in axon guidance pathway genes. Nature. 2012; 491(7424):399–405. doi:10.1038/nature11547.

  9. 9

    Viale A, Pettazzoni P, Lyssiotis CA, Ying H, Sánchez N, Marchesini M, Carugo A, Green T, Seth S, Giuliani V, Kost-Alimova M, Muller F, Colla S, Nezi L, Genovese G, Deem AK, Kapoor A, Yao W, Brunetto E, Kang Y, Yuan M, Asara JM, Wang YA, Heffernan TP, Kimmelman AC, Wang H, Fleming JB, Cantley LC, DePinho RA, Draetta GF. Oncogene ablation-resistant pancreatic cancer cells depend on mitochondrial function. Nature. 2014; 514(7524):628–32. doi:10.1038/nature13611.

  10. 10

    Cohen R, Neuzillet C, Tijeras-Raballand A, Faivre S, de Gramont A, Raymond E. Targeting cancer cell metabolism in pancreatic adenocarcinoma. Oncotarget. 2015; 6(19):16832–47. doi:10.18632/oncotarget.4160.

  11. 11

    Rahrmann EP, Watson AL, Keng VW, Choi K, Moriarity BS, Beckmann DA, Wolf NK, Sarver A, Collins MH, Moertel CL, Wallace MR, Gel B, Serra E, Ratner N, Largaespada DA. Forward genetic screen for malignant peripheral nerve sheath tumor formation identifies new genes and pathways driving tumorigenesis. Nat Genet. 2013; 45(7):756–66. doi:10.1038/ng.2641.

  12. 12

    Cao Z, Wu X, Yen L, Sweeney C, Carraway KL 3rd. Neuregulin-induced erbb3 downregulation is mediated by a protein stability cascade involving the e3 ubiquitin ligase nrdp1. Mol Cell Biol. 2007; 27(6):2180–188. doi:10.1128/MCB.01245-06.

  13. 13

    Gu Y, Kaufman JL, Bernal L, Torre C, Matulis SM, Harvey RD, Chen J, Sun SY, Boise LH, Lonial S. Mln4924, an nae inhibitor, suppresses akt and mtor signaling via upregulation of redd1 in human myeloma cells. Blood. 2014; 123(21):3269–76. doi:10.1182/blood-2013-08-521914.

  14. 14

    Asnaghi L, Bruno P, Priulla M, Nicolin A. mtor: a protein kinase switching between life and death. Pharmacol Res. 2004; 50(6):545–9. doi:10.1016/j.phrs.2004.03.007.

  15. 15

    Sampson VB, Rong NH, Han J, Yang Q, Aris V, Soteropoulos P, Petrelli NJ, Dunn SP, Krueger LJ. Microrna let-7a down-regulates myc and reverts myc-induced growth in burkitt lymphoma cells. Cancer Res. 2007; 67(20):9762–70. doi:10.1158/0008-5472.CAN-07-2462.

  16. 16

    Scarola M, Schoeftner S, Schneider C, Benetti R. mir-335 directly targets rb1 (prb/p105) in a proximal connection to p53-dependent stress response. Cancer Res. 2010; 70(17):6925–33. doi:10.1158/0008-5472.CAN-10-0141.

  17. 17

    Sun F, Fu H, Liu Q, Tie Y, Zhu J, Xing R, Sun Z, Zheng X. Downregulation of ccnd1 and cdk6 by mir-34a induces cell cycle arrest. FEBS Lett. 2008; 582(10):1564–68. doi:10.1016/j.febslet.2008.03.057.

  18. 18

    Dziembowska M, Danilkiewicz M, Wesolowska A, Zupanska A, Chouaib S, Kaminska B. Cross-talk between smad and p38 mapk signalling in transforming growth factor beta signal transduction in human glioblastoma cells. Biochem Biophys Res Commun. 2007; 354(4):1101–06. doi:10.1016/j.bbrc.2007.01.113.

  19. 19

    Wenner CE, Yan S. Biphasic role of tgf-beta1 in signal transduction and crosstalk. J Cell Physiol. 2003; 196(1):42–50. doi:10.1002/jcp.10243.

  20. 20

    Bardeesy N, Cheng KH, Berger JH, Chu GC, Pahler J, Olson P, Hezel AF, Horner J, Lauwers GY, Hanahan D, DePinho RA. Smad4 is dispensable for normal pancreas development yet critical in progression and tumor biology of pancreas cancer. Genes Dev. 2006; 20(22):3130–46. doi:10.1101/gad.1478706.

Download references

Acknowledgements

Y-YH is partially supported by NIH/NCI grants: P30-CA077598, P30-DK050456, P50-CA101955, U19-CA157345, R01 CA179246 and American Cancer Society grant 125627-RSG-14-074-01-TBG. TKS is supported by funding from the NIH NCI (5R00CA151672-04) and the Masonic Cancer Center NIH support grant (P30-CA77598). DL is partially supported by P50-CA101955.

Author information

Correspondence to Yen-Yi Ho.

Additional information

Competing interests

DL is share holder and consultant for NeoClone Biotechnology, Inc and Discovery Genomics, Inc. and received grant support from Genentech, Inc. Y-Y H, TKS, and RL declare that they have no competing interests.

Authors’ contributions

Y-YH and DL developed analysis approaches. Y-YH wrote R scripts and performed the analysis of the data. RL, TKS, DL and Y-YH interpreted the analysis results and wrote the manuscript. All authors read and approved the final manuscript.

Additional files

Additional file 1

Significant CIS enriched pathways and disrupted genes. Number in the parenthesis represents # of cases with the gene mutated. (CSV 13.4 kb)

Additional file 2

Histogram of gene set sizes. (PNG 25.3 kb)

Additional file 3

CIS and pathways association table. (TXT 3.87 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ho, Y., Starr, T.K., LaRue, R.S. et al. Case-oriented pathways analysis in pancreatic adenocarcinoma using data from a sleeping beauty transposon mutagenesis screen. BMC Med Genomics 9, 16 (2016) doi:10.1186/s12920-016-0176-7

Download citation

Keywords

  • Forward genetic screen
  • Sleeping Beauty transposon
  • Case-oriented gene set analysis
  • Pathways correlations
  • CIS
  • Common insertion sites