Detection of Simultaneous Group Effects in MicroRNA Expression and Related Target Gene Sets

Expression levels of mRNAs are among other factors regulated by microRNAs. A particular microRNA can bind specifically to several target mRNAs and lead to their degradation. Expression levels of both, mRNAs and microRNAs, can be obtained by microarray experiments. In order to increase the power of detecting microRNAs that are differentially expressed between two different groups of samples, we incorporate expression levels of their related target gene sets. Group effects are determined individually for each microRNA, and by enrichment tests and global tests for target gene sets. The resulting lists of p-values from individual and set-wise testing are combined by means of meta analysis. We propose a new approach to connect microRNA-wise and gene set-wise information by means of p-value combination as often used in meta-analysis. In this context, we evaluate the usefulness of different approaches of gene set tests. In a simulation study we reveal that our combination approach is more powerful than microRNA-wise testing alone. Furthermore, we show that combining microRNA-wise results with ‘competitive’ gene set tests maintains a pre-specified false discovery rate. In contrast, a combination with ‘self-contained’ gene set tests can harm the false discovery rate, particularly when gene sets are not disjunct.

[1]  Chi-Ying F. Huang,et al.  miRTarBase: a database curates experimentally validated microRNA–target interactions , 2010, Nucleic Acids Res..

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

[3]  Y. Asmann,et al.  miRNA Expression in Colon Polyps Provides Evidence for a Multihit Model of Colon Cancer , 2011, PloS one.

[4]  Stijn van Dongen,et al.  miRBase: tools for microRNA genomics , 2007, Nucleic Acids Res..

[5]  G. Calin,et al.  Complex Patterns of Altered MicroRNA Expression during the Adenoma-Adenocarcinoma Sequence for Microsatellite-Stable Colorectal Cancer , 2011, Clinical Cancer Research.

[6]  S. Dudoit,et al.  Multiple Hypothesis Testing in Microarray Experiments , 2003 .

[7]  J. Barker,et al.  Integrating microRNA and mRNA expression profiles of neuronal progenitors to identify regulatory networks underlying the onset of cortical neurogenesis , 2009, BMC Neuroscience.

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

[9]  Y. Benjamini,et al.  THE CONTROL OF THE FALSE DISCOVERY RATE IN MULTIPLE TESTING UNDER DEPENDENCY , 2001 .

[10]  M. Kendall.  Statistical Methods for Research Workers , 1937, Nature.

[11]  S. Stouffer.  Adjustment during army life , 1977 .

[12]  Anton J. Enright,et al.  Reciprocal regulation of microRNA and mRNA profiles in neuronal development and synapse formation , 2009, BMC Genomics.

[13]  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 .

[14]  Guillemette Marot,et al.  Statistical Applications in Genetics and Molecular Biology Sequential Analysis for Microarray Data Based on Sensitivity and Meta-Analysis , 2011 .

[15]  Sergio Contrino,et al.  ArrayExpress—a public repository for microarray gene expression data at the EBI , 2004, Nucleic Acids Res..

[16]  H. Le,et al.  Comparative Expression Profile of miRNA and mRNA in Primary Peripheral Blood Mononuclear Cells Infected with Human Immunodeficiency Virus (HIV-1) , 2011, PloS one.

[17]  Peter N. Robinson,et al.  GOing Bayesian: model-based gene set analysis of genome-scale data , 2010, Nucleic acids research.

[18]  Gary A. Churchill,et al.  Analysis of Variance for Gene Expression Microarray Data , 2000, J. Comput. Biol..

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

[20]  L. Lim,et al.  MicroRNA targeting specificity in mammals: determinants beyond seed pairing. , 2007, Molecular cell.

[21]  Tsung-Cheng Chang,et al.  microRNAs in vertebrate physiology and human disease. , 2007, Annual review of genomics and human genetics (Print).

[22]  Kim W. Carter,et al.  Integrated Analysis of miRNA and mRNA Expression in Childhood Medulloblastoma Compared with Neural Stem Cells , 2011, PloS one.

[23]  Terence P. Speed,et al.  A comparison of normalization methods for high density oligonucleotide array data based on variance and bias , 2003, Bioinform..

[24]  Di Wu,et al.  ROAST: rotation gene set tests for complex microarray experiments , 2010, Bioinform..

[25]  U. Mansmann,et al.  Testing Differential Gene Expression in Functional Groups , 2005, Methods of Information in Medicine.

[26]  John D. Storey,et al.  Empirical Bayes Analysis of a Microarray Experiment , 2001 .

[27]  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.

[28]  Edgar Brunner,et al.  Comparison of global tests for functional gene sets in two-group designs and selection of potentially effect-causing genes , 2011, Bioinform..

[29]  G. A. Whitmore,et al.  Importance of replication in microarray gene expression studies: statistical methods and evidence from repetitive cDNA hybridizations. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

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

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

[32]  Belinda Phipson,et al.  Opposing roles of polycomb repressive complexes in hematopoietic stem and progenitor cells. , 2010, Blood.

[33]  Peter Bühlmann,et al.  Analyzing gene expression data in terms of gene sets: methodological issues , 2007, Bioinform..

[34]  Jelle J. Goeman,et al.  A global test for groups of genes: testing association with a clinical outcome , 2004, Bioinform..

[35]  Ulrich Mansmann,et al.  GlobalANCOVA: exploration and assessment of gene group effects , 2008, Bioinform..

[36]  Ulrich Mansmann,et al.  Multiple testing on the directed acyclic graph of gene ontology , 2008, Bioinform..

[37]  Øyvind Langsrud,et al.  Rotation tests , 2005, Stat. Comput..