Identification of key regulators of pancreatic cancer progression through multidimensional systems-level analysis

BackgroundPancreatic cancer is an aggressive cancer with dismal prognosis, urgently necessitating better biomarkers to improve therapeutic options and early diagnosis. Traditional approaches of biomarker detection that consider only one aspect of the biological continuum like gene expression alone are limited in their scope and lack robustness in identifying the key regulators of the disease. We have adopted a multidimensional approach involving the cross-talk between the omics spaces to identify key regulators of disease progression.MethodsMultidimensional domain-specific disease signatures were obtained using rank-based meta-analysis of individual omics profiles (mRNA, miRNA, DNA methylation) related to pancreatic ductal adenocarcinoma (PDAC). These domain-specific PDAC signatures were integrated to identify genes that were affected across multiple dimensions of omics space in PDAC (genes under multiple regulatory controls, GMCs). To further pin down the regulators of PDAC pathophysiology, a systems-level network was generated from knowledge-based interaction information applied to the above identified GMCs. Key regulators were identified from the GMC network based on network statistics and their functional importance was validated using gene set enrichment analysis and survival analysis.ResultsRank-based meta-analysis identified 5391 genes, 109 miRNAs and 2081 methylation-sites significantly differentially expressed in PDAC (false discovery rate ≤ 0.05). Bimodal integration of meta-analysis signatures revealed 1150 and 715 genes regulated by miRNAs and methylation, respectively. Further analysis identified 189 altered genes that are commonly regulated by miRNA and methylation, hence considered GMCs. Systems-level analysis of the scale-free GMCs network identified eight potential key regulator hubs, namely E2F3, HMGA2, RASA1, IRS1, NUAK1, ACTN1, SKI and DLL1, associated with important pathways driving cancer progression. Survival analysis on individual key regulators revealed that higher expression of IRS1 and DLL1 and lower expression of HMGA2, ACTN1 and SKI were associated with better survival probabilities.ConclusionsIt is evident from the results that our hierarchical systems-level multidimensional analysis approach has been successful in isolating the converging regulatory modules and associated key regulatory molecules that are potential biomarkers for pancreatic cancer progression.

[1]  Jin‐Young Jang,et al.  Erratum to: Integrative analysis of multi-omics data for identifying multi-markers for diagnosing pancreatic cancer , 2015, BMC Genomics.

[2]  Xia Li,et al.  Integrating analysis reveals microRNA-mediated pathway crosstalk among Crohn's disease, ulcerative colitis and colorectal cancer. , 2014, Molecular bioSystems.

[3]  Jignesh M. Patel,et al.  Michigan molecular interactions r2: from interacting proteins to pathways , 2008, Nucleic Acids Res..

[4]  Ali S. Hadi,et al.  Finding Groups in Data: An Introduction to Chster Analysis , 1991 .

[5]  A. Feinberg,et al.  Hypomethylation distinguishes genes of some human cancers from their normal counterparts , 1983, Nature.

[6]  N. Jamieson,et al.  Tissue Biomarkers for Prognosis in Pancreatic Ductal Adenocarcinoma: A Systematic Review and Meta-analysis , 2011, Clinical Cancer Research.

[7]  A. Barabasi,et al.  Lethality and centrality in protein networks , 2001, Nature.

[8]  Douglas G Altman,et al.  Key Issues in Conducting a Meta-Analysis of Gene Expression Microarray Datasets , 2008, PLoS medicine.

[9]  Rafael A Irizarry,et al.  Exploration, normalization, and summaries of high density oligonucleotide array probe level data. , 2003, Biostatistics.

[10]  Benjamin Haibe-Kains,et al.  A comparative study of survival models for breast cancer prognostication based on microarray data: does a single gene beat them all? , 2008, Bioinform..

[11]  Jianwei Niu,et al.  A Novel Complex Networks Clustering Algorithm Based on the Core Influence of Nodes , 2014, TheScientificWorldJournal.

[12]  Andrea Califano,et al.  Reverse engineering of TLX oncogenic transcriptional networks identifies RUNX1 as tumor suppressor in T-ALL , 2011, Nature Medicine.

[13]  Ka Yee Yeung,et al.  Principal component analysis for clustering gene expression data , 2001, Bioinform..

[14]  Mariano J. Alvarez,et al.  Identification of Causal Genetic Drivers of Human Disease through Systems-Level Analysis of Regulatory Networks , 2014, Cell.

[15]  Mariano J. Alvarez,et al.  A human B-cell interactome identifies MYB and FOXM1 as master regulators of proliferation in germinal centers , 2010, Molecular systems biology.

[16]  Giovanni Scardoni,et al.  Analyzing biological network parameters with CentiScaPe , 2009, Bioinform..

[17]  C. Borror.  Nonparametric Statistical Methods, 2nd, Ed. , 2001 .

[18]  Jianmin Wu,et al.  Genome‐wide DNA methylation patterns in pancreatic ductal adenocarcinoma reveal epigenetic deregulation of SLIT‐ROBO, ITGA2 and MET signaling , 2014, International journal of cancer.

[19]  R. Jove,et al.  miR‐194 is a marker of hepatic epithelial cells and suppresses metastasis of liver cancer cells in mice , 2010, Hepatology.

[20]  Gordon K Smyth,et al.  Statistical Applications in Genetics and Molecular Biology Linear Models and Empirical Bayes Methods for Assessing Differential Expression in Microarray Experiments , 2011 .

[21]  Audrey Kauffmann,et al.  Bioinformatics Applications Note Arrayqualitymetrics—a Bioconductor Package for Quality Assessment of Microarray Data , 2022 .

[22]  A. Jemal,et al.  Cancer statistics, 2014 , 2014, CA: a cancer journal for clinicians.

[23]  Philip H. Ramsey.  Nonparametric Statistical Methods , 1974, Technometrics.

[24]  Rainer Breitling,et al.  Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments , 2004, FEBS letters.

[25]  Alex E. Lash,et al.  Gene Expression Omnibus: NCBI gene expression and hybridization array data repository , 2002, Nucleic Acids Res..

[26]  N. Dubrawsky.  Cancer statistics , 1989, CA: a cancer journal for clinicians.

[27]  K. Sneppen,et al.  Specificity and Stability in Topology of Protein Networks , 2002, Science.

[28]  R. Bourgon,et al.  miR-221/222 Targeting of Trichorhinophalangeal 1 (TRPS1) Promotes Epithelial-to-Mesenchymal Transition in Breast Cancer , 2011, Science Signaling.

[29]  Nini Rao,et al.  [Principal component analysis for exploring gene expression patterns]. , 2007, Sheng wu yi xue gong cheng xue za zhi = Journal of biomedical engineering = Shengwu yixue gongchengxue zazhi.

[30]  C. Croce,et al.  A microRNA expression signature of human solid tumors defines cancer gene targets , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[31]  K. Ohuchida,et al.  MicroRNA-10a is Overexpressed in Human Pancreatic Cancer and Involved in Its Invasiveness Partially via Suppression of the HOXA1 Gene , 2012, Annals of Surgical Oncology.

[32]  A. Jemal,et al.  Cancer Statistics, 2009 , 2009, CA: a cancer journal for clinicians.

[33]  Trey Ideker,et al.  Cytoscape 2.8: new features for data integration and network visualization , 2010, Bioinform..

[34]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[35]  Mariza de Andrade,et al.  Leukocyte DNA Methylation Signature Differentiates Pancreatic Cancer Patients from Healthy Controls , 2011, PloS one.

[36]  E. Giovannetti,et al.  microRNAs with prognostic significance in pancreatic ductal adenocarcinoma: A meta-analysis. , 2015, European journal of cancer.

[37]  R. Tsien,et al.  Specificity and Stability in Topology of Protein Networks , 2022 .

[38]  T. Ideker,et al.  Network-based classification of breast cancer metastasis , 2007, Molecular systems biology.

[39]  Chin-Lee Wu,et al.  Aberrant splicing of cables gene, a CDK regulator, in human cancers , 2005, Cancer biology & therapy.

[40]  Xuerui Yang,et al.  An Extensive MicroRNA-Mediated Network of RNA-RNA Interactions Regulates Established Oncogenic Pathways in Glioblastoma , 2011, Cell.

[41]  M. Korc,et al.  Angiogenic gene signature in human pancreatic cancer correlates with TGF-beta and inflammatory transcriptomes , 2015, Oncotarget.

[42]  Kenneth M. Yamada,et al.  Tumor Suppressor PTEN Inhibits Integrin- and Growth Factor–mediated Mitogen-activated Protein (MAP) Kinase Signaling Pathways , 1998, The Journal of cell biology.

[43]  Chengcheng Hu,et al.  Differentially Expressed MicroRNAs in Postpartum Breast Cancer in Hispanic Women , 2015, PloS one.

[44]  J. Habermann,et al.  Serum biomarkers for improved diagnostic of pancreatic cancer: a current overview , 2011, Journal of Cancer Research and Clinical Oncology.

[45]  Anne E Carpenter,et al.  ZFHX4 interacts with the NuRD core member CHD4 and regulates the glioblastoma tumor-initiating cell state. , 2014, Cell reports.

[46]  H. Xiong,et al.  TLR signaling by tumor and immune cells: a double-edged sword , 2008, Oncogene.

[47]  Wolfgang Schramm,et al.  Team , 2018, Spaces of Intensity.

[48]  Xiangyu Kong,et al.  Aberrant MicroRNAs in Pancreatic Cancer: Researches and Clinical Implications , 2014, Gastroenterology research and practice.

[49]  J. Mesirov,et al.  Interpreting patterns of gene expression with self-organizing maps: methods and application to hematopoietic differentiation. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[50]  Gary D Bader,et al.  GeneMANIA: Fast gene network construction and function prediction for Cytoscape , 2014, F1000Research.

[51]  Rafael A Irizarry,et al.  Processing of Agilent microRNA array data , 2010, BMC Research Notes.

[52]  D. Hwang,et al.  Characterization of gene expression and activated signaling pathways in solid-pseudopapillary neoplasm of pancreas , 2014, Modern Pathology.

[53]  How Early Can Pancreatic Cancer Be Recognized? A Case Report and Review of the Literature , 2015, Case Reports in Oncology.

[54]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[55]  Rainer Breitling,et al.  RankProd: a bioconductor package for detecting differentially expressed genes in meta-analysis , 2006, Bioinform..

[56]  P. López-Romero.  Pre-processing and differential expression analysis of Agilent microRNA arrays using the AgiMicroRna Bioconductor library , 2011, BMC Genomics.

[57]  Philippe Bardou,et al.  jvenn: an interactive Venn diagram viewer , 2014, BMC Bioinformatics.

[58]  H. Horvitz,et al.  MicroRNA expression profiles classify human cancers , 2005, Nature.

[59]  Li Wang,et al.  Identifying core gene modules in glioblastoma based on multilayer factor-mediated dysfunctional regulatory networks through integrating multi-dimensional genomic data , 2015, Nucleic acids research.

[60]  Joshua M. Stuart,et al.  The Cancer Genome Atlas Pan-Cancer analysis project , 2013, Nature Genetics.

[61]  Syed Haider,et al.  International Cancer Genome Consortium Data Portal—a one-stop shop for cancer genomics data , 2011, Database J. Biol. Databases Curation.

[62]  N. Shen,et al.  MicroRNA-195 targets ADP-ribosylation factor-like protein 2 to induce apoptosis in human embryonic stem cell-derived neural progenitor cells , 2013, Cell Death and Disease.

[63]  Mariano J. Alvarez,et al.  Cross-species regulatory network analysis identifies a synergistic interaction between FOXM1 and CENPF that drives prostate cancer malignancy. , 2014, Cancer cell.

[64]  Paul C. Boutros,et al.  NanoStringNorm: an extensible R package for the pre-processing of NanoString mRNA and miRNA data , 2012, Bioinform..

[65]  Andrea Califano,et al.  The TLX1 oncogene drives aneuploidy in T-cell transformation , 2010, Nature Medicine.

[66]  Kyoung Mii Park,et al.  CDA: Combinatorial Drug Discovery Using Transcriptional Response Modules , 2012, PloS one.

[67]  Chittibabu Guda,et al.  A Meta Analysis of Pancreatic Microarray Datasets Yields New Targets as Cancer Genes and Biomarkers , 2014, PloS one.

[68]  Multiple functions of hypoxia-regulated miR-210 in cancer , 2014, Journal of experimental & clinical cancer research : CR.

[69]  L. Tanoue.  Cancer Statistics, 2009 , 2010 .

[70]  R. Dalla‐Favera,et al.  Mutations of multiple genes cause deregulation of NF-κB in diffuse large B-cell lymphoma , 2009, Nature.

[71]  Pan Du,et al.  lumi: a pipeline for processing Illumina microarray , 2008, Bioinform..