Dissecting the biological relationship between TCGA miRNA and mRNA sequencing data using MMiRNA-Viewer

BackgroundMicroRNAs (miRNA) are short nucleotides that interact with their target genes through 3′ untranslated regions (UTRs). The Cancer Genome Atlas (TCGA) harbors an increasing amount of cancer genome data for both tumor and normal samples. However, there are few visualization tools focusing on concurrently displaying important relationships and attributes between miRNAs and mRNAs of both cancer tumor and normal samples. Moreover, a deep investigation of miRNA-mRNA target and biological relationships across multiple cancer types by integrating web-based analysis has not been thoroughly conducted.ResultsWe developed an interactive visualization tool called MMiRNA-Viewer that can concurrently present the co-relationships of expression between miRNA-mRNA pairs of both tumor and normal samples into a single graph. The input file of MMiRNA-Viewer contains the expression information including fold changes between normal and tumor samples for mRNAs and miRNAs, the correlation between mRNA and miRNA, and the predicted target relationship by a number of databases. Users can also load their own input data into MMiRNA-Viewer and visualize and compare detailed information about cancer-related gene expression changes, and also changes in the expression of transcription-regulating miRNAs.To validate the MMiRNA-Viewer, eight types of TCGA cancer datasets with both normal and control samples were selected in this study and three filter steps were applied subsequently. We performed Gene Ontology (GO) analysis for genes available in final selected 238 pairs and also for genes in the top 5 % (95 percentile) for each of eight cancer types to report a significant number of genes involved in various biological functions and pathways. We also calculated various centrality measurement matrices for the largest connected component(s) in each of eight cancers and reported top genes and miRNAs with high centrality measurements.ConclusionsWith its user-friendly interface, dynamic visualization and advanced queries, we also believe MMiRNA-Viewer offers an intuitive approach for visualizing and elucidating co-relationships between miRNAs and mRNAs of both tumor and normal samples. We suggest that miRNA and mRNA pairs with opposite fold changes of their expression and with inverted correlation values between tumor and normal samples might be most relevant for explaining the decoupling of mRNAs and their targeting miRNAs in tumor samples for certain cancer types.

[1]  Yuzhang Wu,et al.  miR-1207-5p and miR-1266 suppress gastric cancer growth and invasion by targeting telomerase reverse transcriptase , 2014, Cell Death and Disease.

[2]  Marian Hajduch,et al.  Diagnostic and prognostic potential of miR-21, miR-29c, miR-148 and miR-203 in adenocarcinoma and squamous cell carcinoma of esophagus , 2015, Diagnostic Pathology.

[3]  Yao Lin,et al.  miR-183 inhibits invasion of gastric cancer by targeting Ezrin. , 2014, International journal of clinical and experimental pathology.

[4]  Rong Luo,et al.  A novel centrality method for weighted networks based on the Kirchhoff polynomial , 2015, Pattern Recognit. Lett..

[5]  Fengming Luo,et al.  MiR-940 inhibits hepatocellular carcinoma growth and correlates with prognosis of hepatocellular carcinoma patients , 2015, Cancer science.

[6]  Mustafa Ozen,et al.  The role of miR-145 in stem cell characteristics of human laryngeal squamous cell carcinoma Hep-2 cells , 2016, Tumor Biology.

[7]  Ying Zhu,et al.  MicroRNA-484 is more highly expressed in serum of early breast cancer patients compared to healthy volunteers , 2014, BMC Cancer.

[8]  Panayiota Poirazi,et al.  A new microRNA target prediction tool identifies a novel interaction of a putative miRNA with CCND2 , 2012, RNA biology.

[9]  Patrick Neven,et al.  Genetic variation at CYP3A is associated with age at menarche and breast cancer risk: a case-control study , 2014, Breast Cancer Research.

[10]  Athanasios Fevgas,et al.  DIANA-TarBase v7.0: indexing more than half a million experimentally supported miRNA:mRNA interactions , 2014, Nucleic Acids Res..

[11]  Hui Zhang,et al.  MiR-183/-96/-182 cluster is up-regulated in most breast cancers and increases cell proliferation and migration , 2014, Breast Cancer Research.

[12]  D. Bartel,et al.  Predicting effective microRNA target sites in mammalian mRNAs , 2015, eLife.

[13]  Minoru Kanehisa,et al.  KEGG as a reference resource for gene and protein annotation , 2015, Nucleic Acids Res..

[14]  R. Harris,et al.  Positively correlated miRNA-mRNA regulatory networks in mouse frontal cortex during early stages of alcohol dependence , 2013, BMC Genomics.

[15]  A. Lund,et al.  MicroRNA and cancer , 2012, Molecular oncology.

[16]  Li Haifeng,et al.  Screening of Target Genes and Regulatory Function of miRNAs as Prognostic Indicators for Prostate Cancer , 2015, Medical science monitor : international medical journal of experimental and clinical research.

[17]  Mustafa Ozen,et al.  Role of miR‐145 in human laryngeal squamous cell carcinoma , 2016, Head & neck.

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

[19]  Dahu Wang,et al.  Hsa-miR-21 and Hsa-miR-29 in Tissue as Potential Diagnostic and Prognostic Biomarkers for Gastric Cancer , 2015, Cellular Physiology and Biochemistry.

[20]  Wei Wang,et al.  Five miRNAs Considered as Molecular Targets for Predicting Esophageal Cancer , 2015, Medical science monitor : international medical journal of experimental and clinical research.

[21]  Pornpimol Charoentong,et al.  ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks , 2009, Bioinform..

[22]  Qiang Li,et al.  hsa-miR-141 downregulates TM4SF1 to inhibit pancreatic cancer cell invasion and migration. , 2014, International journal of oncology.

[23]  Matteo Fassan,et al.  A 4‐MicroRNA signature can discriminate primary lymphomas from anaplastic carcinomas in thyroid cytology smears , 2014, Cancer cytopathology.

[24]  Brad T. Sherman,et al.  Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists , 2008, Nucleic acids research.

[25]  C. Sander,et al.  Analysis of microRNA-target interactions across diverse cancer types , 2013, Nature Structural &Molecular Biology.

[26]  Guangfeng Niu,et al.  miR‐454 is down‐regulated in osteosarcomas and suppresses cell proliferation and invasion by directly targeting c‐Met , 2015, Cell proliferation.

[27]  Alan F. Scott,et al.  Online Mendelian Inheritance in Man (OMIM), a knowledgebase of human genes and genetic disorders , 2002, Nucleic Acids Res..

[28]  K. Tomczak,et al.  The Cancer Genome Atlas (TCGA): an immeasurable source of knowledge , 2015, Contemporary oncology.

[29]  Yadong Wang,et al.  miR2Disease: a manually curated database for microRNA deregulation in human disease , 2008, Nucleic Acids Res..

[30]  Thilo Gambichler,et al.  Microarray analysis of microRNA expression in cutaneous squamous cell carcinoma. , 2012, Journal of dermatological science.

[31]  Dechun Liu,et al.  Lymph node ratio and breast cancer prognosis: a meta-analysis , 2013, Breast Cancer.

[32]  Chuan Li,et al.  A six-microRNA set as prognostic indicators for bile duct cancer. , 2015, International journal of clinical and experimental medicine.

[33]  Yang Liu,et al.  Correlating Bladder Cancer Risk Genes with Their Targeting MicroRNAs Using MMiRNA-Tar , 2015, Genom. Proteom. Bioinform..

[34]  Eleazar Eskin,et al.  Genome reassembly with high-throughput sequencing data , 2013, BMC Genomics.

[35]  Syed Haider,et al.  Ensembl BioMarts: a hub for data retrieval across taxonomic space , 2011, Database J. Biol. Databases Curation.

[36]  George S. Krasnov,et al.  CrossHub: a tool for multi-way analysis of The Cancer Genome Atlas (TCGA) in the context of gene expression regulation mechanisms , 2016, Nucleic acids research.

[37]  Mingming Jia,et al.  COSMIC: exploring the world's knowledge of somatic mutations in human cancer , 2014, Nucleic Acids Res..

[38]  Tao Wang,et al.  Negative correlation of ITCH E3 ubiquitin ligase and miRNA-106b dictates metastatic progression in pancreatic cancer , 2015, Oncotarget.

[39]  Q. Cui,et al.  An Analysis of Human MicroRNA and Disease Associations , 2008, PloS one.

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

[41]  Hiroshi Egawa,et al.  miR-629 Targets TRIM33 to Promote TGFβ/Smad Signaling and Metastatic Phenotypes in ccRCC , 2014, Molecular Cancer Research.

[42]  Adhemar Longatto-Filho,et al.  MicroRNA expression as risk biomarker of breast cancer metastasis: a pilot retrospective case-cohort study , 2013, BMC Cancer.

[43]  John F. McDonald,et al.  Evidence for the Complexity of MicroRNA-Mediated Regulation in Ovarian Cancer: A Systems Approach , 2011, PloS one.

[44]  Saiful Miah,et al.  Distinct microRNA alterations characterize high- and low-grade bladder cancer. , 2009, Cancer research.