MiRComb: An R Package to Analyse miRNA-mRNA Interactions. Examples across Five Digestive Cancers

MicroRNAs (miRNAs) are small RNAs that regulate the expression of target mRNAs by specific binding on the mRNA 3'UTR and promoting mRNA degradation in the majority of cases. It is often of interest to know the specific targets of a miRNA in order to study them in a particular disease context. In that sense, some databases have been designed to predict potential miRNA-mRNA interactions based on hybridization sequences. However, one of the main limitations is that these databases have too many false positives and do not take into account disease-specific interactions. We have developed an R package (miRComb) able to combine miRNA and mRNA expression data with hybridization information, in order to find potential miRNA-mRNA targets that are more reliable to occur in a specific physiological or disease context. This article summarizes the pipeline and the main outputs of this package by using as example TCGA data from five gastrointestinal cancers (colon cancer, rectal cancer, liver cancer, stomach cancer and esophageal cancer). The obtained results can be used to develop a huge number of testable hypotheses by other authors. Globally, we show that the miRComb package is a useful tool to deal with miRNA and mRNA expression data, that helps to filter the high amount of miRNA-mRNA interactions obtained from the pre-existing miRNA target prediction databases and it presents the results in a standardised way (pdf report). Moreover, an integrative analysis of the miRComb miRNA-mRNA interactions from the five digestive cancers is presented. Therefore, miRComb is a very useful tool to start understanding miRNA gene regulation in a specific context. The package can be downloaded in http://mircomb.sourceforge.net.

[1]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[2]  Elias Campo Guerri,et al.  International network of cancer genome projects , 2010 .

[3]  B. Frey,et al.  Using expression profiling data to identify human microRNA targets , 2007, Nature Methods.

[4]  Tim Beißbarth,et al.  Graph based fusion of miRNA and mRNA expression data improves clinical outcome prediction in prostate cancer , 2011, BMC Bioinformatics.

[5]  Peter Dalgaard,et al.  R Development Core Team (2010): R: A language and environment for statistical computing , 2010 .

[6]  Stijn van Dongen,et al.  miRBase: microRNA sequences, targets and gene nomenclature , 2005, Nucleic Acids Res..

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

[8]  Xinxia Peng,et al.  Computational identification of hepatitis C virus associated microRNA-mRNA regulatory modules in human livers , 2009, BMC Genomics.

[9]  Jeffrey A. Thompson,et al.  Common features of microRNA target prediction tools , 2014, Front. Genet..

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

[11]  S. Pfeffer,et al.  miR-122--a key factor and therapeutic target in liver disease. , 2015, Journal of hepatology.

[12]  A. Roses,et al.  O3-02-07: Identification of miRNA changes in Alzheimer's disease brain and CSF yields putative biomarkers and insights into disease pathways , 2008, Alzheimer's & Dementia.

[13]  Cheng Li,et al.  Adjusting batch effects in microarray expression data using empirical Bayes methods. , 2007, Biostatistics.

[14]  Charity W. Law,et al.  voom: precision weights unlock linear model analysis tools for RNA-seq read counts , 2014, Genome Biology.

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

[16]  Aedín C. Culhane,et al.  RamiGO: an R/Bioconductor package providing an AmiGO Visualize interface , 2013, Bioinform..

[17]  Robert Gentleman,et al.  Using GOstats to test gene lists for GO term association , 2007, Bioinform..

[18]  Gary D Bader,et al.  International network of cancer genome projects , 2010, Nature.

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

[20]  Walter Fontana,et al.  Fast folding and comparison of RNA secondary structures , 1994 .

[21]  Roland Eils,et al.  circlize implements and enhances circular visualization in R , 2014, Bioinform..

[22]  R. Giegerich,et al.  Fast and effective prediction of microRNA/target duplexes. , 2004, RNA.

[23]  A. Roses,et al.  Identification of miRNA Changes in Alzheimer's Disease Brain and CSF Yields Putative Biomarkers and Insights into Disease Pathways , 2008 .

[24]  Russell Bowler,et al.  The multiMiR R package and database: integration of microRNA–target interactions along with their disease and drug associations , 2014, Nucleic acids research.

[25]  Raphael Gottardo,et al.  Orchestrating high-throughput genomic analysis with Bioconductor , 2015, Nature Methods.