Large-scale inference of competing endogenous RNA networks with sparse partial correlation

Abstract Motivation MicroRNAs (miRNAs) are important non-coding post-transcriptional regulators that are involved in many biological processes and human diseases. Individual miRNAs may regulate hundreds of genes, giving rise to a complex gene regulatory network in which transcripts carrying miRNA binding sites act as competing endogenous RNAs (ceRNAs). Several methods for the analysis of ceRNA interactions exist, but these do often not adjust for statistical confounders or address the problem that more than one miRNA interacts with a target transcript. Results We present SPONGE, a method for the fast construction of ceRNA networks. SPONGE uses ’multiple sensitivity correlation’, a newly defined measure for which we can estimate a distribution under a null hypothesis. SPONGE can accurately quantify the contribution of multiple miRNAs to a ceRNA interaction with a probabilistic model that addresses previously neglected confounding factors and allows fast P-value calculation, thus outperforming existing approaches. We applied SPONGE to paired miRNA and gene expression data from The Cancer Genome Atlas for studying global effects of miRNA-mediated cross-talk. Our results highlight already established and novel protein-coding and non-coding ceRNAs which could serve as biomarkers in cancer. Availability and implementation SPONGE is available as an R/Bioconductor package (doi: 10.18129/B9.bioc.SPONGE). Supplementary information Supplementary data are available at Bioinformatics online.

[1]  Alexander van Oudenaarden,et al.  Genome-wide dissection of microRNA functions and cotargeting networks using gene set signatures. , 2010, Molecular cell.

[2]  Yuming Luo,et al.  Linc00511 acts as a competing endogenous RNA to regulate VEGFA expression through sponging hsa‐miR‐29b‐3p in pancreatic ductal adenocarcinoma , 2017, Journal of cellular and molecular medicine.

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

[4]  Seongho Kim.  ppcor: An R Package for a Fast Calculation to Semi-partial Correlation Coefficients. , 2015, Communications for statistical applications and methods.

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

[6]  P. Pandolfi,et al.  A ceRNA Hypothesis: The Rosetta Stone of a Hidden RNA Language? , 2011, Cell.

[7]  C. Ponting,et al.  Single-Cell Multiomics: Multiple Measurements from Single Cells , 2017, Trends in genetics : TIG.

[8]  Beth Israel,et al.  Decision letter: Replication Study: A coding-independent function of gene and pseudogene mRNAs regulates tumour biology , 2010 .

[9]  Mary Goldman,et al.  The UCSC Xena platform for public and private cancer genomics data visualization and interpretation , 2018, bioRxiv.

[10]  Junpeng Zhang,et al.  Inferring miRNA sponge co-regulation of protein-protein interactions in human breast cancer , 2017, BMC Bioinformatics.

[11]  Xia Li,et al.  Identification of lncRNA-associated competing triplets reveals global patterns and prognostic markers for cancer , 2015, Nucleic acids research.

[12]  Yvonne Tay,et al.  Competing endogenous RNA networks: tying the essential knots for cancer biology and therapeutics , 2015, Journal of Hematology & Oncology.

[13]  J. Rinn,et al.  Integrative analyses reveal a long noncoding RNA-mediated sponge regulatory network in prostate cancer , 2016, Nature Communications.

[14]  Mihaela Zavolan,et al.  Argonaute CLIP--a method to identify in vivo targets of miRNAs. , 2012, Methods.

[15]  Prahlad T. Ram,et al.  Cupid: simultaneous reconstruction of microRNA-target and ceRNA networks , 2015, Genome research.

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

[17]  Tao Huang,et al.  Cancer-Related Triplets of mRNA-lncRNA-miRNA Revealed by Integrative Network in Uterine Corpus Endometrial Carcinoma , 2017, BioMed research international.

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

[19]  R. Fisher.  FREQUENCY DISTRIBUTION OF THE VALUES OF THE CORRELATION COEFFIENTS IN SAMPLES FROM AN INDEFINITELY LARGE POPU;ATION , 1915 .

[20]  Tingting Shao,et al.  The mRNA related ceRNA–ceRNA landscape and significance across 20 major cancer types , 2015, Nucleic acids research.

[21]  Weining Yang,et al.  Versican 3′‐untranslated region (3′‐UTR) functions as a ceRNA in inducing the development of hepatocellular carcinoma by regulating miRNA activity , 2013, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[22]  Angel Rubio,et al.  Joint analysis of miRNA and mRNA expression data , 2013, Briefings Bioinform..

[23]  Rui-Xi Hua,et al.  Long Intergenic Noncoding RNA 00511 Acts as an Oncogene in Non–small-cell Lung Cancer by Binding to EZH2 and Suppressing p57 , 2016, Molecular therapy. Nucleic acids.

[24]  C. Burge,et al.  Most mammalian mRNAs are conserved targets of microRNAs. , 2008, Genome research.

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

[26]  Mihaela Zavolan,et al.  Single‐cell mRNA profiling reveals the hierarchical response of miRNA targets to miRNA induction , 2018, Molecular systems biology.

[27]  Guanming Lu,et al.  Long noncoding RNA LINC00511 contributes to breast cancer tumourigenesis and stemness by inducing the miR-185-3p/E2F1/Nanog axis , 2018, Journal of Experimental & Clinical Cancer Research.

[28]  H. Seitz,et al.  microRNA target prediction programs predict many false positives , 2017, Genome research.

[29]  Mary Goldman,et al.  Toil enables reproducible, open source, big biomedical data analyses , 2017, Nature Biotechnology.

[30]  Ling Fang,et al.  Expression of CD44 3′-untranslated region regulates endogenous microRNA functions in tumorigenesis and angiogenesis , 2010, Nucleic acids research.

[31]  Ellen T. Gelfand,et al.  The Genotype-Tissue Expression (GTEx) project , 2013, Nature Genetics.

[32]  Marcel H. Schulz,et al.  JAMI: fast computation of conditional mutual information for ceRNA network analysis , 2018, Bioinform..

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

[34]  Panayiotis Tsanakas,et al.  DIANA-LncBase v2: indexing microRNA targets on non-coding transcripts , 2015, Nucleic Acids Res..

[35]  C. Sander,et al.  Target mRNA abundance dilutes microRNA and siRNA activity , 2010, Molecular systems biology.

[36]  Hsien-Da Huang,et al.  miRTarBase 2016: updates to the experimentally validated miRNA-target interactions database , 2015, Nucleic Acids Res..

[37]  R. Fisher.  Statistical methods for research workers , 1927, Protoplasma.

[38]  Guido Caldarelli,et al.  Web Search Queries Can Predict Stock Market Volumes , 2011, PloS one.

[39]  Xia Li,et al.  Comprehensive characterization of lncRNA-mRNA related ceRNA network across 12 major cancers , 2016, Oncotarget.

[40]  Thomas M. Keane,et al.  The BRAF Pseudogene Functions as a Competitive Endogenous RNA and Induces Lymphoma In Vivo , 2015, Cell.

[41]  Jiayi Wang,et al.  CREB up-regulates long non-coding RNA, HULC expression through interaction with microRNA-372 in liver cancer , 2010, Nucleic acids research.

[42]  John T. Powers,et al.  Multiple mechanisms disrupt the let-7 microRNA family in neuroblastoma , 2016, Nature.

[43]  Junpeng Zhang,et al.  Computational methods for identifying miRNA sponge interactions , 2016, Briefings Bioinform..

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

[45]  Debora S. Marks,et al.  miRcode: a map of putative microRNA target sites in the long non-coding transcriptome , 2012, Bioinform..

[46]  Hui Zhou,et al.  starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein–RNA interaction networks from large-scale CLIP-Seq data , 2013, Nucleic Acids Res..

[47]  Ziv Bar-Joseph,et al.  Reconstructing dynamic microRNA-regulated interaction networks , 2013, Proceedings of the National Academy of Sciences.

[48]  P. Pandolfi,et al.  The multilayered complexity of ceRNA crosstalk and competition , 2014, Nature.

[49]  D. Bartel.  MicroRNAs: Target Recognition and Regulatory Functions , 2009, Cell.

[50]  A. Dreher.  Modeling Survival Data Extending The Cox Model , 2016 .

[51]  Lorenzo Farina,et al.  Computational analysis identifies a sponge interaction network between long non-coding RNAs and messenger RNAs in human breast cancer , 2014, BMC Systems Biology.

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