Search

Author(s): Shivan Sivakumar 1, Ines de Santiago 1, Leon Chlon 1,2, Florian Markowetz 1,2,* Introduction Pancreatic ductal adenocarcinoma (PDAC) is the most lethal human malignancy, with a 5-y survival of [...], Background KRAS is the most frequently mutated gene in pancreatic ductal adenocarcinoma (PDAC), but the mechanisms underlying the transcriptional response to oncogenic KRAS are still not fully understood. We aimed to uncover transcription factors that regulate the transcriptional response of oncogenic KRAS in pancreatic cancer and to understand their clinical relevance. Methods and Findings We applied a well-established network biology approach (master regulator analysis) to combine a transcriptional signature for oncogenic KRAS derived from a murine isogenic cell line with a coexpression network derived by integrating 560 human pancreatic cancer cases across seven studies. The datasets included the ICGC cohort (n = 242), the TCGA cohort (n = 178), and five smaller studies (n = 17, 25, 26, 36, and 36). 55 transcription factors were coexpressed with a significant number of genes in the transcriptional signature (gene set enrichment analysis [GSEA] p < 0.01). Community detection in the coexpression network identified 27 of the 55 transcription factors contributing to three major biological processes: Notch pathway, down-regulated Hedgehog/Wnt pathway, and cell cycle. The activities of these processes define three distinct subtypes of PDAC, which demonstrate differences in survival and mutational load as well as stromal and immune cell composition. The Hedgehog subgroup showed worst survival (hazard ratio 1.73, 95% CI 1.1 to 2.72, coxPH test p = 0.018) and the Notch subgroup the best (hazard ratio 0.62, 95% CI 0.42 to 0.93, coxPH test p = 0.019). The cell cycle subtype showed highest mutational burden (ANOVA p < 0.01) and the smallest amount of stromal admixture (ANOVA p < 2.2e-16). This study is limited by the information provided in published datasets, not all of which provide mutational profiles, survival data, or the specifics of treatment history. Conclusions Our results characterize the regulatory mechanisms underlying the transcriptional response to oncogenic KRAS and provide a framework to develop strategies for specific subtypes of this disease using current therapeutics and by identifying targets for new groups.