MicroRNA-Target Network Inference and Local Network Enrichment Analysis Identify Two microRNA Clusters with Distinct Functions in Head and Neck Squamous Cell Carcinoma

MicroRNAs represent ~22 nt long endogenous small RNA molecules that have been experimentally shown to regulate gene expression post-transcriptionally. One main interest in miRNA research is the investigation of their functional roles, which can typically be accomplished by identification of mi-/mRNA interactions and functional annotation of target gene sets. We here present a novel method “miRlastic”, which infers miRNA-target interactions using transcriptomic data as well as prior knowledge and performs functional annotation of target genes by exploiting the local structure of the inferred network. For the network inference, we applied linear regression modeling with elastic net regularization on matched microRNA and messenger RNA expression profiling data to perform feature selection on prior knowledge from sequence-based target prediction resources. The novelty of miRlastic inference originates in predicting data-driven intra-transcriptome regulatory relationships through feature selection. With synthetic data, we showed that miRlastic outperformed commonly used methods and was suitable even for low sample sizes. To gain insight into the functional role of miRNAs and to determine joint functional properties of miRNA clusters, we introduced a local enrichment analysis procedure. The principle of this procedure lies in identifying regions of high functional similarity by evaluating the shortest paths between genes in the network. We can finally assign functional roles to the miRNAs by taking their regulatory relationships into account. We thoroughly evaluated miRlastic on a cohort of head and neck cancer (HNSCC) patients provided by The Cancer Genome Atlas. We inferred an mi-/mRNA regulatory network for human papilloma virus (HPV)-associated miRNAs in HNSCC. The resulting network best enriched for experimentally validated miRNA-target interaction, when compared to common methods. Finally, the local enrichment step identified two functional clusters of miRNAs that were predicted to mediate HPV-associated dysregulation in HNSCC. Our novel approach was able to characterize distinct pathway regulations from matched miRNA and mRNA data. An R package of miRlastic was made available through: http://icb.helmholtz-muenchen.de/mirlastic.

[1]  H. Zitzelsberger,et al.  Integrative analysis of the microRNA-mRNA response to radiochemotherapy in primary head and neck squamous cell carcinoma cells , 2015, BMC Genomics.

[2]  Fabian J Theis,et al.  Integrative Analysis of MicroRNA and mRNA Data Reveals an Orchestrated Function of MicroRNAs in Skeletal Myocyte Differentiation in Response to TNF-α or IGF1 , 2015, PloS one.

[3]  R. Gregory,et al.  MicroRNA biogenesis pathways in cancer , 2015, Nature Reviews Cancer.

[4]  A. Mäkitie,et al.  Expression of toll-like receptors in HPV-positive and HPV-negative oropharyngeal squamous cell carcinoma—an in vivo and in vitro study , 2015, Tumor Biology.

[5]  W. Xiong,et al.  The Oncogenic MicroRNA miR-21 Promotes Regulated Necrosis in Mice , 2015, Nature Communications.

[6]  Debora S. Marks,et al.  MicroRNA control of protein expression noise , 2015, Science.

[7]  Steven J. M. Jones,et al.  Comprehensive genomic characterization of head and neck squamous cell carcinomas , 2015, Nature.

[8]  Sam Griffiths-Jones,et al.  Bias in microRNA functional enrichment analysis , 2015, Bioinform..

[9]  Xiaowei Wang,et al.  miRDB: an online resource for microRNA target prediction and functional annotations , 2014, Nucleic Acids Res..

[10]  Christopher J. Cheng,et al.  MicroRNA silencing for cancer therapy targeted to the tumor microenvironment , 2014, Nature.

[11]  Hushan Yang,et al.  Current management of locally advanced head and neck cancer: the combination of chemotherapy with locoregional treatments. , 2014, Seminars in oncology.

[12]  Fabian J Theis,et al.  Proteome-wide analysis reveals an age-associated cellular phenotype of in situ aged human fibroblasts , 2014, Aging.

[13]  M. Zavolan,et al.  Identification and consequences of miRNA–target interactions — beyond repression of gene expression , 2014, Nature Reviews Genetics.

[14]  N. Dimitrova,et al.  A novel saliva-based microRNA biomarker panel to detect head and neck cancers , 2014, Cellular Oncology.

[15]  Philip A. Ewels,et al.  Understanding functional miRNA–target interactions in vivo by site-specific genome engineering , 2014, Nature Communications.

[16]  Ryan K. Orosco,et al.  Multi-tiered genomic analysis of head and neck cancer ties TP53 mutation to 3p loss , 2014, Nature Genetics.

[17]  Lili Xiong,et al.  Genome-wide survey of tissue-specific microRNA and transcription factor regulatory networks in 12 tissues , 2014, Scientific Reports.

[18]  J. Machiels,et al.  Targeted therapies for squamous cell carcinoma of the head and neck: current knowledge and future directions. , 2014, Cancer treatment reviews.

[19]  J. Grandis,et al.  Integration of molecular targeted therapy with radiation in head and neck cancer. , 2014, Pharmacology & therapeutics.

[20]  M. Rietbergen,et al.  Cancer stem cell enrichment marker CD98: a prognostic factor for survival in patients with human papillomavirus-positive oropharyngeal cancer. , 2014, European journal of cancer.

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

[22]  C. Belka,et al.  Changes in circulating microRNAs after radiochemotherapy in head and neck cancer patients , 2013, Radiation Oncology.

[23]  Wan Li,et al.  Unraveling the characteristics of microRNA regulation in the developmental and aging process of the human brain , 2013, BMC Medical Genomics.

[24]  A. Contreras-Paredes,et al.  Role of Innate Immunity against Human Papillomavirus (HPV) Infections and Effect of Adjuvants in Promoting Specific Immune Response , 2013, Viruses.

[25]  C. Farah,et al.  MicroRNAs in Head and Neck Cancer , 2013, International journal of dentistry.

[26]  Michael Orth,et al.  Current concepts in clinical radiation oncology , 2013, Radiation and Environmental Biophysics.

[27]  Bryan R. Cullen,et al.  MicroRNA Target Site Identification by Integrating Sequence and Binding Information , 2013, Nature Methods.

[28]  Fabian J. Theis,et al.  The human transcriptome is enriched for miRNA-binding sites located in cooperativity-permitting distance , 2013, RNA biology.

[29]  L. Waldron,et al.  Potentially Prognostic miRNAs in HPV-Associated Oropharyngeal Carcinoma , 2013, Clinical Cancer Research.

[30]  M. Zavolan,et al.  A biophysical miRNA-mRNA interaction model infers canonical and noncanonical targets , 2013, Nature Methods.

[31]  A. Hui,et al.  MicroRNA-193b Enhances Tumor Progression via Down Regulation of Neurofibromin 1 , 2013, PloS one.

[32]  Xiaowei Wang,et al.  A microRNA expression signature for the prognosis of oropharyngeal squamous cell carcinoma , 2013, Cancer.

[33]  Steven J. M. Jones,et al.  Comprehensive genomic characterization of squamous cell lung cancers , 2012, Nature.

[34]  S. Syrjänen,et al.  Biology of Human Papillomavirus Infections in Head and Neck Carcinogenesis , 2012, Head and Neck Pathology.

[35]  N. Barkai,et al.  Efficiency and specificity in microRNA biogenesis , 2012, Nature Structural &Molecular Biology.

[36]  F. Nielsen,et al.  The role of miRNAs in human papilloma virus (HPV)-associated cancers: bridging between HPV-related head and neck cancer and cervical cancer , 2012, British Journal of Cancer.

[37]  A. Pasquinelli MicroRNAs and their targets: recognition, regulation and an emerging reciprocal relationship , 2012, Nature Reviews Genetics.

[38]  A. Luttun,et al.  Quantification of miRNA-mRNA Interactions , 2012, PloS one.

[39]  Susumu Goto,et al.  KEGG for integration and interpretation of large-scale molecular data sets , 2011, Nucleic Acids Res..

[40]  Fang Zhang,et al.  Methylation of microRNA-9 is a specific and sensitive biomarker for oral and oropharyngeal squamous cell carcinomas. , 2012, Oral oncology.

[41]  T. Furukawa,et al.  MicroRNAs Associated with Mitogen-Activated Protein Kinase in Human Pancreatic Cancer , 2011, Molecular Cancer Research.

[42]  Minghua Deng,et al.  A Lasso regression model for the construction of microRNA-target regulatory networks , 2011, Bioinform..

[43]  Fabian J. Theis,et al.  MicroRNAs coordinately regulate protein complexes , 2011, BMC Systems Biology.

[44]  A. McKenna,et al.  The Mutational Landscape of Head and Neck Squamous Cell Carcinoma , 2011, Science.

[45]  Fabian J Theis,et al.  miTALOS: analyzing the tissue-specific regulation of signaling pathways by human and mouse microRNAs. , 2011, RNA.

[46]  Elizabeth E. Hoskins,et al.  Alteration of microRNA profiles in squamous cell carcinoma of the head and neck cell lines by human papillomavirus , 2011, Head & neck.

[47]  G. de Rosa,et al.  The role of human papillomavirus in the pathogenesis of head & neck squamous cell carcinoma: an overview , 2011, Infectious Agents and Cancer.

[48]  D. Hanahan,et al.  Hallmarks of Cancer: The Next Generation , 2011, Cell.

[49]  Hui Zhou,et al.  starBase: a database for exploring microRNA–mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data , 2010, Nucleic Acids Res..

[50]  S. Datta,et al.  Modeling microRNA-mRNA Interactions Using PLS Regression in Human Colon Cancer , 2011, BMC Medical Genomics.

[51]  Tuan D Pham,et al.  Integrative analysis of next generation sequencing for small non-coding RNAs and transcriptional regulation in Myelodysplastic Syndromes , 2011, BMC Medical Genomics.

[52]  F. Slack,et al.  OncomiR addiction in an in vivo model of microRNA-21-induced pre-B-cell lymphoma , 2010, Nature.

[53]  Anjali J. Koppal,et al.  Supplementary data: Comprehensive modeling of microRNA targets predicts functional non-conserved and non-canonical sites , 2010 .

[54]  Cary A Moody,et al.  Human papillomavirus oncoproteins: pathways to transformation , 2010, Nature Reviews Cancer.

[55]  Nicholas T. Ingolia,et al.  Mammalian microRNAs predominantly act to decrease target mRNA levels , 2010, Nature.

[56]  Gabriele Sales,et al.  MAGIA, a web-based tool for miRNA and Genes Integrated Analysis , 2010, Nucleic Acids Res..

[57]  Wei Shi,et al.  Comprehensive MicroRNA Profiling for Head and Neck Squamous Cell Carcinomas , 2010, Clinical Cancer Research.

[58]  Trevor Hastie,et al.  Regularization Paths for Generalized Linear Models via Coordinate Descent. , 2010, Journal of statistical software.

[59]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[60]  Y. Tomari,et al.  Making RISC. , 2010, Trends in biochemical sciences.

[61]  Fabian J Theis,et al.  PhenomiR: a knowledgebase for microRNA expression in diseases and biological processes , 2010, Genome Biology.

[62]  S. Cohen,et al.  Use of microRNA sponges to explore tissue-specific microRNA functions in vivo , 2009, Nature Methods.

[63]  Kathryn A. O’Donnell,et al.  Therapeutic microRNA Delivery Suppresses Tumorigenesis in a Murine Liver Cancer Model , 2009, Cell.

[64]  W. Westra,et al.  Molecular pathology of head and neck cancer: implications for diagnosis, prognosis, and treatment. , 2009, Annual review of pathology.

[65]  S. Weed,et al.  Actin cytoskeletal mediators of motility and invasion amplified and overexpressed in head and neck cancer , 2008, Clinical & Experimental Metastasis.

[66]  Sanghyuk Lee,et al.  miRGator: an integrated system for functional annotation of microRNAs , 2007, Nucleic Acids Res..

[67]  D. DiMaio,et al.  Human papillomavirus in cervical and head-and-neck cancer , 2008, Nature Clinical Practice Oncology.

[68]  Carole Fakhry,et al.  Case-control study of human papillomavirus and oropharyngeal cancer. , 2007, The New England journal of medicine.

[69]  T. Gilmore Introduction to NF-κB: players, pathways, perspectives , 2006, Oncogene.

[70]  B. Wollenberg,et al.  Cytokines in head and neck cancer. , 2006, Cytokine & growth factor reviews.

[71]  F. Slack,et al.  Oncomirs — microRNAs with a role in cancer , 2006, Nature Reviews Cancer.

[72]  Gábor Csárdi,et al.  The igraph software package for complex network research , 2006 .

[73]  H. Zou,et al.  Regularization and variable selection via the elastic net , 2005 .

[74]  D. Bartel,et al.  Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes. , 2005, RNA.

[75]  C. Burge,et al.  Conserved Seed Pairing, Often Flanked by Adenosines, Indicates that Thousands of Human Genes are MicroRNA Targets , 2005, Cell.

[76]  Anton J. Enright,et al.  Human MicroRNA Targets , 2004, PLoS biology.

[77]  Salvatore Pece,et al.  Endocytosis and cancer. , 2004, Current opinion in cell biology.

[78]  C. Burge,et al.  Prediction of Mammalian MicroRNA Targets , 2003, Cell.

[79]  C. Van Waes,et al.  Effects of pharmacologic antagonists of epidermal growth factor receptor, PI3K and MEK signal kinases on NF‐κB and AP‐1 activation and IL‐8 and VEGF expression in human head and neck squamous cell carcinoma lines , 2002, International journal of cancer.

[80]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[81]  D. Brachman,et al.  Molecular biology of head and neck cancer. , 1994, Seminars in oncology.

[82]  Edsger W. Dijkstra,et al.  A note on two problems in connexion with graphs , 1959, Numerische Mathematik.