- Open Access
miRTissue ce: extending miRTissue web service with the analysis of ceRNA-ceRNA interactions
BMC Bioinformatics volume 21, Article number: 199 (2020)
Non-coding RNAs include different classes of molecules with regulatory functions. The most studied are microRNAs (miRNAs) that act directly inhibiting mRNA expression or protein translation through the interaction with a miRNAs-response element. Other RNA molecules participate in the complex network of gene regulation. They behave as competitive endogenous RNA (ceRNA), acting as natural miRNA sponges to inhibit miRNA functions and modulate the expression of RNA messenger (mRNA). It became evident that understanding the ceRNA–miRNA–mRNA crosstalk would increase the functional information across the transcriptome, contributing to identify new potential biomarkers for translational medicine.
We present miRTissue ce, an improvement of our original miRTissue web service. By introducing a novel computational pipeline, miRTissue ce provides an easy way to search for ceRNA interactions in several cancer tissue types. Moreover it extends the functionalities of previous miRTissue release about miRNA-target interaction in order to provide a complete insight about miRNA mediated regulation processes. miRTissue ce is freely available at http://tblab.pa.icar.cnr.it/mirtissue.html.
The study of ceRNA networks and its dynamics in cancer tissue could be applied in many fields of translational biology, as the investigation of new cancer biomarker, both diagnostic and prognostic, and also in the investigation of new therapeutic strategies of intervention. In this scenario, miRTissue ce can offer a powerful instrument for the analysis and characterization of ceRNA-ceRNA interactions in different tissue types, representing a fundamental step in order to understand more complex regulation mechanisms.
The cross-talk between coding and non-coding RNA molecules represents a novel layer of gene regulation. MicroRNAs (miRNA) are defined as "master regulators" of gene expression, as they can influence gene expression of RNA messenger (mRNA) through the direct binding on their MRE (miRNA Responsive Element) sites . Apart from mRNA, other non-coding RNA (ncRNA) molecules can bind miRNAs through MRE sites. They are called competitive endogenous RNA (ceRNA) or miRNA sponge because they compete with mRNA target for miRNA binding [2, 3]. To date, the following ncRNA molecules have been experimentally validated having a ceRNA behaviour: long non-coding RNAs (lncRNAs), pseudogenes, and circular RNAs (circRNA). Some other classes of ncRNA molecules that work as ceRNAs, such as the tRNA-derived small RNAs (tsRNA) that have been recently associated with gene expression regulation in cancer progression [4–6], have a few pieces of evidence not yet experimentally validated.
ncRNAs are defined as regulative molecules because they can interfere with miRNAs in gene regulation process. As a consequence, an ncRNA action can alter a specific miRNA-target interaction through miRNA sequestration. Indeed they seem to have a relevant function in different physiological and pathological conditions as the development of human cancers [7, 8].
Salmena et al.  presented a unifying "ceRNA-hypothesis", in which different ncRNAs classes would actively talk to each other, through a "ceRNA language", composing a large-scale interaction network, and guiding their respective expression levels . This molecular network is dependent on different features as concentration and subcellular distribution of the different RNA molecules considered, and also on the cellular type. Figure 1 shows a schematic overview of that process. All those features can have a different impact on general developmental and metabolic processes in multiple tissues, especially on transcription-related functions . Applying this RNA regulative network in the context of cancer, considering a gene with tumor-suppressor function, a miRNA molecule targeting this mRNA could block its expression, acting as an oncogenic miRNA; in this case, a ceRNA molecule could compete with the tumor-suppressor for the same miRNA and restore the gene function. In the next part of this Section we discuss in detail all the ncRNA molecules with ceRNA function.
lncRNA are RNA molecules, more than 200 nucleotides long, involved in regulation of gene expression. In cancer context, several scientific works demonstrated that they are differentially regulated in tumor tissues compared with normal ones. One example is lncRNA Highly Up-regulated in Liver Cancer (HULC), highly over-expressed in hepatocellular carcinoma (HCC) . MEG3 is another lncRNA which is involved in modulation of apoptosis and autophagy. It acts on p53 gene, enhancing its transcriptional activation, and on MDM2 gene, causing its downregulation. The final result of this molecular cascade is the block of cell-cycle and regulation of autophagy (Fig. 2).
Pseudogenes are segments of DNA highly homologous to their gene counterpart but "apparently" with no functional activity. Poliseno et al. hypothesized for the first time a ceRNA role for the phosphatase and tensin homolog (PTEN) pseudogene 1 (PTENP1) : its sponge effect on PTEN gene leads to a decreased gene expression of PTEN and of other tumor suppressor genes. Moreover authors proved that cells presenting a PTENP1 overexpression show a growth inhibition, deducing a potential role of this pseudogene as tumor suppressor  (Fig. 2).
circular RNA (circRNA) is another class of non coding RNA molecules, called this way because of its covalently closed ring structure . For instance, circTP63 is a circular RNA identified in differ cancer types .
More recently, there are evidences that mRNAs were associated with ceRNA functions. For example Lee et al.  showed that Versican (VCAN) 3’ UTR modulates retinoblatoma 1 (RB1) expression through specific miRNA interaction. It binds miR-199a and miR-144, upregulating the expression of this tumor suppressor gene (Fig. 2).
Considering the biological relevance of ceRNA molecules, in this paper we present an upgrade of the original miRTissue web service , called miRTissue ce, that stands for miRTissue "ceRNA edition", integrating novel knowledge about ceRNA interactions. Built on the same platform and service architecture of its first release, miRTissue ce now provides the possibility to search for ceRNA interaction networks in different tissues. Together with the basic functionalities of original miRTissue, that is the characterization of miRNA target interactions in different tissues and in different conditions (normal/tumour), miRTissue ce offers a novel and complete insight about miRNA mediated regulation processes. Moreover, the study of intra-tissue (normal versus tumour) and inter-tissues (tumor A versus tumor B) correlations allows the user to better address its research of biological cancer markers in the field of precision medicine. Indeed, each tissue type has a proper biologic profile, with a proper expression profile for the same molecules under study, thus, a molecular marker identified for a specific tissue, could not be a marker for another tissue type.
miRTissue ce was developed in order to improve the previous version of the tool, providing information on potential molecular markers in tumor biology, supporting the search for ceRNA interaction networks in many tissues. Its main impact is the research of potential endogenous molecules to use in precision medicine, in the field of prognosis, prediction, and cancer treatments.
The rest of the paper is organized as follows. In the next section, a set of similar bioinformatics services is reviewed. In “Materials and methods” section, we present the data sources, software packages and computational pipeline employed in our service. In “Results” section we describe miRTissue ce’s new features and new functionalities and then we introduce some case studies. In “Discussion” section, we point out the potential impact of our system and provide a comparison with other similar services. Finally “Conclusion” are drawn.
In recent years, some bioinformatics services focused on the analysis of miRNA-mediated interaction networks. One of the first database collecting these interaction networks is miRTarBase . miRTarBase collects only experimentally validated interactions obtained from published papers. Curators also collect information from other databases such as Human MicroRNA Disease Database (HMDD) , National Center for Biotechnology Information (NCBI) Entrez Gene  and NCBI Reference Sequence (RefSeq) , The Cancer Genome Atlas (TCGA) , Gene Expression Omnibus (GEO) , Kyoto Encyclopedia of Genes and Genomes (KEGG)  and Database of Annotation, Visualization and Integrated Discovery (DAVID) . HMDD information are aimed to investigate the relationship between miRNA and human diseases. Among them NCBI RefSeq and Entrez Gene are used to obtain information about the target gene, while information about miRNA expression profiles are obtained from TCGA and Gene Expression Omnibus; DAVID and KEGG are used to obtain functional annotation of miRNAs.
ceRDB  is another database based on similar features. ceRDB provides putative ceRNA interactions starting from putative miRNA-mRNA interactions predicted by TargetScan . The potential for miRNA-mRNA competition for each mRNA is ranked using the count of number of miRNA binding sites shared between the two RNA molecules (miRNA and mRNA). ceRDB considers as putative ceRNA only mRNA molecules.
lnCeDB  collects from TargetScan the putative mRNA targets of human miRNAs, and from starBase  the targets predicted from available Argonaute (AGO) PAR-CLIP datasets, that is data obtained from Argonaute photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation (PAR-CLIP) sequencing. Furthermore, the lncRNA targets of human miRNAs (up to GENCODE 11) are downloaded from miRcode database . To find seed-matched target sites, miRNA targets on the rest of the GENCODE 19 lncRNAs are predicted by the developed prediction algorithm. These putative miRNA-lncRNA interactions are mapped to the AGO protein interacting regions within lncRNAs. The likelihood of a lncRNA-mRNA pair for actually being ceRNA is obtained by a two steps method. First, a ceRNA score is calculated from the ratio of the number of shared MREs between the pair with the total number of MREs of the individual candidate gene. Second, the hypergeometric test, using the number of shared miRNAs between the ceRNA pair against the number of miRNAs interacting with the individual RNAs, is used to calculate the p-value for each ceRNA pair. With the discovery of new ncRNA classes other databases collected and integrated information also on these new ncRNAs.
miRSponge  is a manually curated database providing an experimentally supported resource for miRNA sponges. The data are manually curated, and collected from PubMed literature, using about 1.200 published articles, and the database contains data on 599 miRNA-sponge interactions and 463 ceRNA relationships from 11 species. Database classes include endogenously generated molecules as coding genes, pseudogenes, long noncoding RNAs and circular RNAs, along with exogenously introduced molecules, including viral RNAs and artificial engineered sponges. Considering that lncRNA are not yet fully understood a wide spectrum of information is collected for each entry, such as pathway interactions from KEGG  and BioCarta , also integrated by information from TarBase , miR2Disease  and miRTarBase . The functional information is obtained from Gene Expression Omnibus (GEO)  and supported by a "guilty-by-association" strategy.
LncACTdb 2.0  is another database curating many ceRNA types such as circular RNAs and pseudogenes. It provides a comprehensive information on ceRNAs interaction networks in 23 species and 213 diseases/phenotypes, including 2663 experimentally supported, and manually curated, ceRNA interactions from more than 5000 published works of literature. LncACTdb also identifies and scores candidate lncRNA-associated ceRNA interactions across 33 cancer types from TCGA data, and providing illustration of survival, network and cancer hallmark information for ceRNAs. The system recalls the ceRNA interactions related to a lncRNA given as query, and then the user can access a set of tools to obtain more information. In particular new ceRNA interactions can be identified by integrating the expression profiles associated to a disease or a phenotype, and lncRNa functions can be studied through the pathways downloaded from repositories like KEGG, Biocarta and Reactome. Moreover new biomarkers can be evidenced by searching the survival time related to 33 kinds of cancer disease in TCGA database.
miRcode  is a searchable map of putative target sites across the whole set of sequences stored in GENCODE database of annotated transcriptome. The putative sites are obtained scanning the sequences for seed complementary fragments and considering the evolutionary conservation of the fragment. Evolutionary conservation is assessed by using multiple alignment of vertebrate sequences. Conservation level is evaluated checking the presence of the sequence in primates, non-primates mammals, and non-mammals vertebrates.
starBase v2.0  is one of the first database collecting ceRNA interaction networks and it is also one of the most complete. starBase data sets are generated by 37 independent studies and the system can systematically identify the RNA–RNA, and protein–RNA interaction networks from 108 CLIP-Seq (PAR-CLIP, HITS-CLIP, iCLIP, CLASH). starBase v2.0 includes two web servers, miRFunction and ceRNAFunction, to predict the function of miRNAs and other ncRNAs from the miRNA-mediated regulatory networks. The database content includes, among others, the annotation and identification of miRNA-mRNA and miRNA-ncRNA interactions, and the annotation and identification of miRNA-mediated ceRNA regulatory networks. Regarding this content, the pipeline proposed in starBase combines CLIP-supported miRNA-mRNA, miRNA-lncRNA, miRNA-circRNA and miRNA-pseudogene interactions, then the hypergeometric test is implemented to predict ceRNA pairs among mRNAs, lncRNAs, circRNAs and pseudogenes. Recently, starBase v2.0 has been evolved and integrated into ENCORI: The Encyclopedia of RNA Interactomes (http://starbase.sysu.edu.cn). ENCORI is an open-source platform for studying the miRNA-ncRNA, miRNA-mRNA, ncRNA-mRNA, mRNA-mRNA, RBP-ncRNA, and RBP-mRNA interactions from CLIP-seq, degradome-seq and RNA-RNA interactome data.
Materials and methods
In this Section, we describe data sources, computational tools and pipeline used in order to develop the new release of miRTissue web service. Moreover we provide a brief summary of the main functionalities of the earlier miRTissue release.
One of the data source of the presented pipeline are a set of RNA-target interactions, summarized in Table 1. To be more precise, the interactions are composed of a set of validated miRNA-target interactions (miRTarBase); a set of experimental validated ceRNA interactions (miRSponge and LncACTdb 2.0); a set of predicted miRNA-lncRNa and miRNA-pseudogenes interactions (miRcode); a set of experimentally validated miRNA-target interactions (miRTarBase).
RNA expression profiles
The other data sources of the computational pipeline are the expression profiles of several RNA molecules. Those data are extracted by the Cancer Genome Atlas project (TCGA, https://cancergenome.nih.gov/). TCGA repository was created in 2005 to improve the understanding of genetic bases of cancer disease, building a catalog of all genetic mutation responsible for cancer. In the following years the project demonstrated that a mutation atlas could be created for specific cancer types and that the public availability of the material enable researcher to make and validate research discoveries. Today TCGA stores data of more than 30 cancer types, including expression profiles of both coding and non-coding genes of tumor and healthy tissues In this work we consider expression profiles of the following bio types: protein coding, miRNA, pseudogenes, lncRNA. Among miRNA sponges, we will not discuss circRNA and tsRNA. The former because their expression values are not available in TCGA, the latter because there are only few pieces of evidence validating them as ceRNAs .
miRTissue original release
miRTissue is a web service that allows to search for a tissue-specific characterization of miRNA-target interactions in human. Given a set of validated miRNA-target interactions provided by miRTarBase, the interaction type is computed according to a statistical correlation measure, using the global test , among the expression profiles of miRNAs, their mRNA targets and their corresponding proteins. Resulting interactions can be sorted by p-value or interaction types, and it is possible to organize the results for one or more tissue types, for example breast or colon. Given a single tissue, it is possible to check if the interaction type changes according to the tissue condition, that is normal or tumour. More details about the first release of miRTissue are available in .
SPONGE software package
In the last years, many software packages have been released in order to compute in silico ceRNA interactions. For a review we recommend these papers [36–38]. For instance, authors in  proposed a new computational method, called ceRNA predIction Algorithm (CERNIA). CERNIA can be used to study the ceRNA competition among different tissue types, and different classes of genes. The ceRNAs prediction method is based on the DT-Hybrid recommendation algorithm [39, 40].
One of the most recent and most performing ceRNA predictor, according to its authors, is SPONGE (Sparse Partial correlation ON Gene Expression) . SPONGE is an R software library that allows for large-scale inference of ceRNA interactions using a statistical correlation measure called multiple sensibility correlation (mscor). mscor needs an in-deep analysis because it is one of the main steps of miRTissue ce processing pipeline. In particular, mscor is an extension of the basic sensitivity correlation, introduced by , that takes into account the effect of multiple miRNAs for the regulation of ceRNA interactions. It is defined as follows:
where ga and gb are the two genes involved in the ceRNA interaction, M is the set of shared miRNAs between ga and gb,cor() is the Pearson correlation between expression profiles, pcor() is the partial correlation that estimates how two variables are correlated when they control additional variables. In this context, mscor gives an indication about the direct interaction, or less, between the two competing genes (ga and gb) . Indeed, a value of mscor close to zero means a direct interaction between the two observed variables, that is they are low sensitive to the presence of miRNAs. On the other hand, a value of mscor close to cor(ga,gb) means that there is a great contribution of the explanatory variables, that is the miRNAs, leading to an indirect correlation between ga and gb, or, in other words, their interaction is miRNA-mediated. As explained in [34, 36], computational methods for inferring ceRNA interactions, based on statistical correlation, consider the following events: positive correlation between expression profiles of the two involved ceRNAs; negative correlation between the expression profiles of shared miRNA and the ceRNAs. Moreover, SPONGE method defines an mscor null distribution that allows to estimate an empirical p-value for mscor.
miRTissue ce computational pipeline
The new functionalities of miRTissue ce service are provided by means of a three-step novel computational pipeline (Fig. 3) with regards to the one presented in . The data sources are a set of RNA-interactions (see Table 1) and the expression profiles of considered RNA molecules taken from TCGA repository, as previously explained.
In order to compute a set of predicted ceRNA interactions, the first step (A) of our pipeline is to filter a statistical meaningful set of triplets in the form of ceRNA-A (coding gene) - ceRNA-B (RNA acting as sponge, including coding genes, pseudogenes and lncRNA) - list of putative miRNAs interacting with both ceRNA-A and ceRNA-B. As suggested in , this first step is done using the hypergeometric test. Given a couple of RNA, ceRNA-A and ceRNA-B, the hypergeometric test checks for the significance of the set of miRNAs interacting with both ceRNAs. The associated p-value to the test can be computed using the following formula:
where N is the total number of available miRNAs, K is the number of miRNAs that interact with ceRNA-A, O is the number of miRNAs interacting with ceRNA-B, and x is the number of miRNAs shared by both ceRNA-A and ceRNA-B. The hypergeometric test is computed using the GDCRNATools R library . The next filtering step (B) considers only the miRNA-target interactions that show a negative correlation between their expression values. That because, as explained in [16, 44], a miRNA and a target gene, in case they actually interacts, should have anti-correlated expression values. This is computed following the same approach as in miRTissue first release, that is using the globaltest statistical test [35, 44]. More details about the global test and its applications in bioinformatics can be found in [16, 35, 44]. The last step (C) of our computational pipeline is the application of the SPONGE method to the candidate ceRNA interactions obtained in the previous steps. The result of the pipeline is a table of ceRNA interactions in the form of ceRNA-A - ceRNA-B - number of shared miRNAs - list of shared miRNAs - mscor - p-value (Table 2).
In this section, we first present the new functionalities of miRTissue ce, then we discuss the "bioinformatics scenarios", previously presented in the earlier release of MirTissue web service, that now integrates the ceRNA interaction analysis: (1) ceRNA therapeutics analysis in cancer, (2) biomarker discovery in cancer, (3) ceRNA interaction network analysis in cancer.
miRTissue ce new functionalities
All the functionalities of miRTissue ce, including those of the previous release, are explained in detail in Additional file 1. In this Section, we briefly describe the new functionalities. miRTissue ce offers two new use cases that can be selected from the home page.
In the first use case, called ceRNA interaction analysis-compare among different tissues, the user can input from a self completing list one or more mRNA name (ceRNA-A) and/or one or more ceRNA molecule name (ceRNA-B) and can select a list of different tissue. Then, for each tissue, it is possible to visualize the desired ceRNA pairs, if available, and the corresponding p-value. Moreover, by clicking on the p-value referred to a specific tissue, it is possible to switch to the visualization of the involved miRNA list in that tissue (Figure A3 of Additional file 1).
In the second use case, called ceRNA interaction analysis-details for a specific tumor type, the user can input from a self completing list one or more mRNA name (ceRNA-A) and/or one or more ceRNA molecule name (ceRNA-B). Then, after selecting one tissue type it is possible to obtain the related list of ceRNA interactions. Each pair is scored by a p-value and there is also a list of miRNAs which ceRNA-A and ceRNA-B compete for. By clicking on a miRNA in the list, it is possible to visualize all its mRNA targets for that tissue. The last column shows the number of other available tissues that exhibit that ceRNA pair (Figure A4 of Additional file 1).
Finally, the original miRNA-target use case, called details for a specific tumor type has been improved. Now, in fact, it is possible to realize if a miRNA-gene pair is involved in a ceRNA-ceRNA interaction (Figure A2 of Additional file 1).
Bioinformatics scenario 1: ceRNA therapeutics analysis in cancer
In last decade many studies focused their attention on the clinical use of synthetic ncRNA molecules for cancer treatments, as miRNA mimics or antagomirs, that mimic the functions of endogenous miRNAs, acting respectively on oncogenes or tumor-suppressor genes . The identification of ncRNAs that are involved in cellular processes, contribute to oncogenesis, or tumor suppression, or act on miRNAs (well known regulators of gene expression), provides opportunities to develop novel therapeutics for cancer based on targeting lncRNA [46–48]. Several in silico molecules have been developed for this purpose, as small interfering RNAs (siRNAs), acting on specific ncRNAs , Natural Antisense Transcripts (NATs) , and AntiSense Oligonucleotide (ASO) . Indeed, the development of this in silico molecules, can be an RNA-based drugs strategy for precision medicine, not only for treating cancer but also for modulating cancer treatment sensitivity. These molecules offer many advantages compared to other targeted therapies, such as fewer side effects and high specificity compared to other chemical drugs . Moreover some cancer sub-types, as triple negative breast cancer (TNBC), do not have adequate therapies, because of the molecular characteristics of this tumor type, and ceRNA molecules could offer a great opportunity to individuate and intervene on this BC sub-type . Taking into account this scenario, we designed miRTissue ce to allow the researcher to investigate on ceRNA therapeutics analysis in cancer (Fig. 4). To address this issue, few steps are required: first, selection of genes of a specific pathway of interest; second, identification of RNA expression and ceRNA network interaction starting from the selected genes in a specific tumor tissue; third, identification of potential ceRNA molecules (mRNA, lncRNA or pseudogene) acting on oncogenes or tumor-suppressors. The second step can be achieved through correlation and statistical analysis offered by miRTissue ce (Fig. 4). This scenario could be useful to investigate new synthetic molecules acting as ceRNA antagonists or agonists to use in cancer therapeutics.
In order to better illustrate the use of miRTissue ce in the analysis of ceRNA therapeutics application in cancer we report two specific examples:
A) filtering for Breast Cancer (BRCA) tissue, ceRNA network analysis evidenced an interaction between LIMA1 gene, has-miR-20b and has-miR-17 cluster and ncRNA RP11-1000B6.3 (Ensembl ID: ENSG00000261064). See Fig. 5 for interaction details. Both has-miR-20b and has-miR-17 are OncomiRs , acting on LIMA1 tumor suppressor gene . Their interaction leads to LIMA1 degradation, thus decreasing its expression in breast cancer tissue (p-value <0.05). miRTissue ce evidenced also RP11-1000B6.3 lncRNA, competing with LIMA1 for has-miR-20b and has-miR-17 cluster (p-value <0.05). A synthetic ncRNA molecule mimic RP11-1000B6.3 function could be used as therapeutic strategy to sequester those miRNAs and avoiding LIMA1 degradation.
B) Another example, of ceRNA interaction evidenced by miRTissue ce is in colorectal adenocarcinoma (COAD). See Fig. 6 for interaction details. MEIS2 oncogene interacts with has-miR-192, thus decreasing its expression. Experimental evidences support the role of MEIS2 as oncogene  in colon cancer tissue and the role of tumor suppressor of has-miR-192 . They interact with lncRNA RP11-553L6.5. Synthetic ncRNA acting as antagonist could be used to inhibit RP11-553L6.5 function and restore the tumor suppressive miRNA function.
Bioinformatics scenario 2: biomarker discovery in cancer
Cancer biomarkers are biological molecules widely used in different cancer stages, from diagnosis, disease staging or prediction, to clinical prognosis. Indeed they can indicate cancer evolution. There are biomarkers clearly related to disease initiation, or progression, and treatment effects . ncRNAs could be considered good biomarkers for many reasons : for example, their biological role is associated to gene expression regulation, and they are present in different body fluids, including whole blood, plasma, urine, saliva, and gastric juice [58, 59], making them easily detectable .
To address the identification of new molecular biomarkers in cancer, few steps are required: first, identification of expression profiles of genes, miRNAs and other ncRNA molecules; second, reconstruction of ceRNA interaction network (miRNA-target, miRNA-ncRNA, ncRNA-target); third, functional characterization and pathway analysis of cancer specific ceRNA network. In this context, the use of miRTissue ce would easily allow the researchers to get information on ceRNA interactions in the specific tissue under investigation, both considering expression values for all type of RNA molecules for each tissue, and evaluating the correlation between each couple of interacting molecules of the ceRNA network (Fig. 7). Indeed, this allows the user to know what are the exacts molecules expressed in a given tissue, as they could potentially interact with each other but not be expressed equally in all different tissues. Finally, functional characterization and pathway analysis would gain insight into the underlying biology of the ceRNA network and would determine the potential roles of lncRNAs that are aberrantly expressed in cancer. As result of this 3 step procedure, new potential ncRNA biomarkers could be evidenced. In order to better illustrate the use of miRTissue ce in the analysis of biomarker discovery in cancer, we report a specific example: First of all, filtering for Breast Cancer (BRCA) tissue and for PTEN tumor suppressor gene, ceRNA network analysis evidenced an interaction between PTEN, many miRNAs, and ceRNA interactors (mRNAs, lncRNAs and pseudogenes). See Figs. 8 and 9 for interaction details. Figure 9 shows just the top 10 enriched categories, ranked for Fold Enrichment values. The full list of enriched genes is available as supplementary material (Additional file 2). We downloaded all ceRNAs interacting with PTEN gene and then we analyzed them through gene enrichment and pathway analysis. Gene enrichment evidenced an over-representation of the different Biological Processes (BP), Molecular Functions (MF) and cellular Component (CC) linked to regulation of the cellular bio-synthetic process, regulation of signal transduction, apoptotic process, Ras protein signal transduction, and cell death. The pathway analysis also showed enriched terms as p53 protein pathway, RAS protein pathway and apoptosis signalling. Indeed PTEN is a tumor suppressor gene which, in breast cancer cells, is involved in suppression of cell growth by phosphatase activity-dependent G1 arrest followed by cell death . Moreover, a loss of PTEN expression is associated with poor outcome in this tumor type [62, 63]. Some of ceRNAs, as TMBIM6, are linked to pathways previously cited, and they interact with oncogenic miRNAs (hsa-miR-17/has-miR-20b) that regulate these biological processes. This analysis can, thus, evidence ceRNA biomarkers able to bind and block the action of oncogenic miRNA on tumor suppressor genes as PTEN in specific molecular networks (Figs. 8 and 9).
Bioinformatics scenario 3: ceRNA interaction network analysis in cancer
According to recent discoveries, mRNA can be considered a ceRNA molecule . A mRNA-mRNA pair can, indeed, compete for the same miRNA or for the same miRNA cluster. Considering two different tissue conditions, the alteration of mRNA-mRNA interaction, translating it into protein-protein interaction (PPI) network, would potentially lead to modification of biological networks and cellular pathway dynamics. The impact of ceRNA network analysis in clinical research, and consequently the use of miRTissue ce is linked to the deregulation and complexity of molecular interactions rather than a single molecule in complex diseases as is cancer. Through miRTissue ce it is possible to evidence and study changes in the dynamic of the network in two different tissue types, by using these few steps: first, selection of differentially expressed (DE) genes of a specific pathway in two different tissue types (for instance normal and cancer tissue for the same tissue); second, reconstruction of ceRNA interaction network in the two conditions; and third, comparison of network modules and investigation of pathway dynamics (Fig. 10). miRTissue ce allows to evidence and compare RNA interactions and consequently PPI network alterations.
Tumor biology is a complex scenario to investigate, due to its inter and intra-individual variability related to tumor tissue composition, in terms of both the wide numbers of molecules belonging to each tissue, and of quantitative expression. Those features can significantly modify pathway networks and the biological behaviour of a given tumor. In the context of precision medicine, in fact, each cancer patient has a unique patho-physiologic profile, strictly dependent on its molecular and genetic profiles. The knowledge of a cancer patient’s transcriptome can offer a comprehensive view of molecular patterns linked to cancer and can allow identifying biological molecules used as risk, diagnostic or prognostic markers as well as therapeutic targets.
miRTissue ce was developed in order to provide information on ceRNA interaction networks in many tissues, and thus to address clinical research studies focusing on single patient molecular profile. It has features and functionalities that make it different with regards to other bioinformatics databases that provides similar information. Table 3 reports a comparison of miRTissue ce features with other ceRNA databases. In particular, we identified the following seven main features:
ceRNA interaction data source: the data repositories used to find ceRNA interactions;
ceRNA classes: the considered bio-types of the molecules involved in ceRNA interactions;
ceRNA network prediction algorithm: the algorithm, or computational pipelines, used to infer ceRNA networks;
TCGA expression profiles for ceRNA interactions: whether or not TCGA expression profiles of ceRNA molecules are used;
ceRNA interaction score p-value: whether or not a p-value over the score of a ceRNA interaction is used;
multiple ceRNA analysis: whether or not it is possible to carry on analysis considering more than one ceRNA at a time;
multiple tissue selection: whether or not it is possible to carry on analysis considering more than a tissue at a time.
miRTissue ce, through the SPONGE algorithm, is the only system that provides a p-value related to a score, the mscor in this case, associated to a ceRNA interaction. Other systems that implement hypergeometric test, like for instance starBase v2.0, just provides a simple p-value associated to a candidate triplet ceRNA-A/ceRNA-B/miRNA, but there is not a score associated to that triplet, nor a corresponding p-value. One of the main characteristics of miRTissue ce tool is the possibility to analyse many ceRNA triplets (ceRNA-A/miRNA/ceRNA-B) at the same time. Moreover the tool allows to analyse many tissues at the same time, giving the possibility to compare different tissues for the same ceRNA triplet under investigation. Other databases, such as ceRDB or miRSponge, give the user the possibility to analyse just one ceRNA at time, showing the miRNA cluster that interacts with the RNA target, and a list of putative ceRNAs. Moreover, miRTissue ce is one of the few services that allows the user to rank ceRNA interactions by using a p-value score.
ceRNA interactions are all extracted from TCGA repository, this implies that ceRNA interaction analysis is homogeneous. Finally miRTissue ce, with respect to the most complete databases such as LncACTdb and ENCORI, implements in its pipeline the SPONGE algorithm that, at the best of our knowledge, is the most performing tools in order to compute ceRNA interactions.
We presented miRTissue ce, a web service that represents an improvement of our original miRTissue release, extending its focus on the characterization of miRNA-target interactions. miRTissue ce, in fact, integrates data about RNA-target interactions and their expression profiles in order to provide, for different tissue types, a set of ceRNA-ceRNA interactions. RNA molecules acting as ceRNA, indeed, are fundamental elements in the cross-talk mechanism that involves gene expression regulation. The analysis of ceRNA networks in cancer tissues enables to widen clinical research in field such as biomarker discoveries and therapeutic strategies, as explained in the bioinformatics scenarios. With regards to existing similar bioinformatics services, miRTissue ce implements a computation pipeline based on the state-of-the-art algorithm for inferring ceRNA-ceRNA interactions and it provides an easy way to search for ceRNA interactions in several cancer tissue types. Together with the functionalities of the previous release, miRTissue ce offers a complete insight about miRNA mediated regulation processes.
Availability of data and materials
miRTissue ce is freely available at http://tblab.pa.icar.cnr.it/mirtissue.html.
Vompeting endogenous RNA database
Competitive endogenous RNA
ceRNA predIction Algorithm
Cross-linking, ligation and sequencing of hybrids
Database of annotation, visualization and integrated discovery
The Encyclopedia of RNA interactomes
Gene expression omnibus
Human MicroRNA disease database
Highly up-regulated in liver cancer
Kyoto encyclopedia of genes and genomes
Long non-coding RNA
miRNA responsive element
Multiple sensibility correlation
Natural antisense transcript
National center for biotechnology information
Photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation
- PTEN Phosphatase and tensin; PTENP1:
Phosphatase and tensin homolog pseudogene 1
NCBI reference sequence
Small interfering RNA
Sparse partial correlation on gene expression
The cancer genome atlas
Triple negative breast cancer
tRNA-derived small RNA
Alshalalfa M. Microrna response elements-mediated mirna-mirna interactions in prostate cancer. Adv Bioinforma. 2012; 2012. https://doi.org/10.1155/2012/839837.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A cerna hypothesis: the rosetta stone of a hidden rna language?Cell. 2011; 146(3):353–8.
Zhang J, Le TD, Liu L, Li J. Inferring mirna sponge co-regulation of protein-protein interactions in human breast cancer. BMC Bioinformatics. 2017; 18(1):243.
Zhang Y, Deng Q, Tu L, Lv D, Liu D. trna-derived small rnas: A novel class of small rnas in human hypertrophic scar fibroblasts. Int J Mol Med. 2020; 45(1):115–30.
Balatti V, Pekarsky Y, Croce CM. Role of the trna-derived small rnas in cancer: new potential biomarkers and target for therapy. In: Advances in Cancer Research, vol. 135. Elsevier: 2017. p. 173–87. https://doi.org/10.1016/bs.acr.2017.06.007.
La Ferlita A, Alaimo S, Veneziano D, Nigita G, Balatti V, Croce CM, Ferro A, Pulvirenti A. Identification of trna-derived ncrnas in tcga and nci-60 panel cell lines and development of the public database trfexplorer. Database. 2019; 2019. https://doi.org/10.1093/database/baz115.
Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of cerna crosstalk and competition. Nature. 2014; 505(7483):344.
Le TD, Zhang J, Liu L, Li J. Computational methods for identifying mirna sponge interactions. Brief Bioinforma. 2016; 18(4):577–90.
Poliseno L, Salmena L, Zhang J, Carver B, Haveman WJ, Pandolfi PP. A coding-independent function of gene and pseudogene mrnas regulates tumour biology. Nature. 2010; 465(7301):1033.
Xu J, Feng L, Han Z, Li Y, Wu A, Shao T, Ding N, Li L, Deng W, Di X, et al.Extensive cerna–cerna interaction networks mediated by mirnas regulate development in multiple rhesus tissues. Nucleic Acids Res. 2016; 44(19):9438–51.
Wang J, Liu X, Wu H, Ni P, Gu Z, Qiao Y, Chen N, Sun F, Fan Q. Creb up-regulates long non-coding rna, hulc expression through interaction with microrna-372 in liver cancer. Nucleic Acids Res. 2010; 38(16):5366–83.
Zhou Y, Zhong Y, Wang Y, Zhang X, Batista DL, Gejman R, Ansell PJ, Zhao J, Weng C, Klibanski A. Activation of p53 by meg3 non-coding rna. J Biol Chem. 2007; 282(34):24731–42.
Kristensen LS, Andersen MS, Stagsted LV, Ebbesen KK, Hansen TB, Kjems J. The biogenesis, biology and characterization of circular rnas. Nat Rev Genet. 2019; 20(11):675–91.
Cheng Z, Yu C, Cui S, Wang H, Jin H, Wang C, Li B, Qin M, Yang C, He J, et al.c irctp63 functions as a cerna to promote lung squamous cell carcinoma progression by upregulating foxm1. Nat Commun. 2019; 10(1):1–13.
Lee DY, Jeyapalan Z, Fang L, Yang J, Zhang Y, Yee AY, Li M, Du WW, Shatseva T, Yang BB. Expression of versican 3’-untranslated region modulates endogenous microrna functions. PloS ONE. 2010; 5(10):13599.
Fiannaca A, La Rosa M, La Paglia L, Urso A. miRTissue: a web application for the analysis of miRNA-target interactions in human tissues. BMC Bioinformatics. 2018; 19(S15):434. https://doi.org/10.1186/s12859-018-2418-5.
Huang H-Y, Lin Y-C-D, Li J, Huang K-Y, Shrestha S, Hong H-C, Tang Y, Chen Y-G, Jin C-N, Yu Y. mirtarbase 2020: updates to the experimentally validated microrna–target interaction database. Nucleic Acids Res. 2019. https://doi.org/10.1093/nar/gkz896.
Huang Z, Shi J, Gao Y, Cui C, Zhang S, Li J, Zhou Y, Cui Q. Hmdd v3. 0: a database for experimentally supported human microrna–disease associations. Nucleic Acids Res. 2018; 47(D1):1013–7.
Maglott D, Ostell J, Pruitt KD, Tatusova T. Entrez gene: gene-centered information at ncbi. Nucleic Acids Res. 2005; 33(suppl_1):54–8.
Pruitt KD, Tatusova T, Brown GR, Maglott DR. Ncbi reference sequences (refseq): current status, new features and genome annotation policy. Nucleic Acids Res. 2011; 40(D1):130–5.
Tomczak K, Czerwińska P, Wiznerowicz M. The cancer genome atlas (tcga): an immeasurable source of knowledge. Contemp Oncol. 2015; 19(1A):68.
Clough E, Barrett T. The gene expression omnibus database. In: Statistical Genomics. Springer: 2016. p. 93–110. https://doi.org/10.1007/978-1-4939-3578-9_5.
Kanehisa M, et al.The kegg database. In: Novartis Foundation Symposium. Wiley Online Library: 2002. p. 91–100. https://doi.org/10.1002/0470857897.ch8.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. David: database for annotation, visualization, and integrated discovery. Genome Biol. 2003; 4(9):60.
Sarver AL, Subramanian S. Competing endogenous rna database. Bioinformation. 2012; 8(15):731.
Grimson A, Farh KK-H, Johnston WK, Garrett-Engele P, Lim LP, Bartel DP. Microrna targeting specificity in mammals: determinants beyond seed pairing. Mol Cell. 2007; 27(1):91–105.
Das S, Ghosal S, Sen R, Chakrabarti J. lncedb: Database of human long noncoding rna acting as competing endogenous rna. PLoS ONE. 2014; 9(6):1–7. https://doi.org/10.1371/journal.pone.0098965.
Li J-H, Liu S, Zhou H, Qu L-H, Yang J-H. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein–RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2013; 42(D1):92–7. https://doi.org/10.1093/nar/gkt1248. http://oup.prod.sis.lan/nar/article-pdf/42/D1/D92/25891302/gkt1248.pdf.
Jeggari A, Marks DS, Larsson E. miRcode: a map of putative microRNA target sites in the long non-coding transcriptome. Bioinformatics. 2012; 28(15):2062–3. https://doi.org/10.1093/bioinformatics/bts344.
Wang P, Zhi H, Zhang Y, Liu Y, Zhang J, Gao Y, Guo M, Ning S, Li X. miRSponge: a manually curated database for experimentally supported miRNA sponges and ceRNAs. Database J Biol Databases Curation. 2015; 2015. https://doi.org/10.1093/database/bav098.
Nishimura D. Biocarta. Biotech Softw Internet Rep Comput Softw J Scient. 2001; 2(3):117–20.
Sethupathy P, Corda B, Hatzigeorgiou AG. Tarbase: A comprehensive database of experimentally supported animal microrna targets. Rna. 2006; 12(2):192–7.
Jiang Q, Wang Y, Hao Y, Juan L, Teng M, Zhang X, Li M, Wang G, Liu Y. mir2disease: a manually curated database for microrna deregulation in human disease. Nucleic Acids Res. 2008; 37(suppl_1):98–104.
Wang P, Li X, Gao Y, Guo Q, Wang Y, Fang Y, Ma X, Zhi H, Zhou D, Shen W, Liu W, Wang L, Zhang Y, Ning S, Li X. LncACTdb 2.0: an updated database of experimentally supported ceRNA interactions curated from low- and high-throughput experiments. Nucleic Acids Res. 2019; 47(D1):121–7. https://doi.org/10.1093/nar/gky1144.
Goeman JJ, van de Geer SA, de Kort F, van Houwelingen HC. A global test for groups of genes: testing association with a clinical outcome. Bioinformatics. 2004; 20(1):93–9.
Le TD, Zhang J, Liu L, Li J. Computational methods for identifying miRNA sponge interactions. Brief Bioinforma. 2016:042. https://doi.org/10.1093/bib/bbw042.
Veneziano D, Marceca GP, Di Bella S, Nigita G, Distefano R, Croce CM. Investigating miRNA–lncRNA Interactions: Computational Tools and Resources. In: MicroRNA Target Identification. New York, New York, USA: Humana Press: 2019. p. 251–77. https://doi.org/10.1007/978-1-4939-9207-2_14.
Sardina DS, Alaimo S, Ferro A, Pulvirenti A, Giugno R. A novel computational method for inferring competing endogenous interactions. Brief Bioinforma. 2017; 18(6):1071–81. https://doi.org/10.1093/bib/bbw084.
Alaimo S, Bonnici V, Cancemi D, Ferro A, Giugno R, Pulvirenti A. DT-Web: a web-based application for drug-target interaction and drug combination prediction through domain-tuned network-based inference. BMC Syst Biol. 2015; 9(3):4. https://doi.org/10.1186/1752-0509-9-S3-S4.
Alaimo S, Pulvirenti A, Giugno R, Ferro A. Drug–target interaction prediction through domain-tuned network-based inference. Bioinformatics. 2013; 29(16):2004–8. https://doi.org/10.1093/bioinformatics/btt307. https://academic.oup.com/bioinformatics/article-pdf/29/16/2004/16913926/btt307.pdf. .
List M, Dehghani Amirabad A, Kostka D, Schulz MH. Large-scale inference of competing endogenous RNA networks with sparse partial correlation. Bioinformatics. 2019; 35(14):596–604. https://doi.org/10.1093/bioinformatics/btz314.
Paci P, Colombo T, Farina L. Computational analysis identifies a sponge interaction network between long non-coding RNAs and messenger RNAs in human breast cancer. BMC Syst Biol. 2014; 8(1):83. https://doi.org/10.1186/1752-0509-8-83.
Li R, Qu H, Wang S, Wei J, Zhang L, Ma R, Lu J, Zhu J, Zhong W-D, Jia Z. GDCRNATools: an R/Bioconductor package for integrative analysis of lncRNA, miRNA and mRNA data in GDC. Bioinformatics. 2018; 34(14):2515–7. https://doi.org/10.1093/bioinformatics/bty124.
van Iterson M, Bervoets S, de Meijer EJ, Buermans HP, ’t Hoen PAC, Menezes RX, Boer JM. Integrated analysis of microRNA and mRNA expression: adding biological significance to microRNA target predictions. Nucleic Acids Res. 2013; 41(15):146. https://doi.org/10.1093/nar/gkt525.
Nguyen D-D, Chang S. Development of novel therapeutic agents by inhibition of oncogenic micrornas. Int J Mol Sci. 2018; 19(1):65.
Kimura T, Aikata H, Takahashi S, Takahashi I, Nishibuchi I, Doi Y, Kenjo M, Murakami Y, Honda Y, Kakizawa H, et al.Stereotactic body radiotherapy for patients with small hepatocellular carcinoma ineligible for resection or ablation therapies. Hep Res. 2015; 45(4):378–86.
Sasaki YT, Ideue T, Sano M, Mituyama T, Hirose T. Men ε/ β noncoding rnas are essential for structural integrity of nuclear paraspeckles. Proc Natl Acad Sci. 2009; 106(8):2525–30.
Parasramka MA, Maji S, Matsuda A, Yan IK, Patel T. Long non-coding rnas as novel targets for therapy in hepatocellular carcinoma. Pharmacol Ther. 2016; 161:67–78.
Modarresi F, Faghihi MA, Lopez-Toledano MA, Fatemi RP, Magistri M, Brothers SP, Van Der Brug MP, Wahlestedt C. Inhibition of natural antisense transcripts in vivo results in gene-specific transcriptional upregulation. Nat Biotechnol. 2012; 30(5):453.
Gutschner T, Hämmerle M, Eißmann M, Hsu J, Kim Y, Hung G, Revenko A, Arun G, Stentrup M, Groß M, et al.The noncoding rna malat1 is a critical regulator of the metastasis phenotype of lung cancer cells. Cancer Res. 2013; 73(3):1180–9.
Ahmadzada T, Reid G, McKenzie DR. Fundamentals of sirna and mirna therapeutics and a review of targeted nanoparticle delivery systems in breast cancer. Biophys Rev. 2018; 10(1):69–86.
Bobbili MR, Mader RM, Grillari J, Dellago H. Oncomir-17-5p: alarm signal in cancer?Oncotarget. 2017; 8(41):71206.
Collins RJ, Jiang WG, Hargest R, Mason MD, Sanders AJ. Eplin: a fundamental actin regulator in cancer metastasis?Cancer Metastasi Rev. 2015; 34(4):753–64.
Wan Z, Chai R, Yuan H, Chen B, Dong Q, Zheng B, Mou X, Pan W, Tu Y, Yang Q, et al.Meis2 promotes cell migration and invasion in colorectal cancer. Oncol Rep. 2019; 42(1):213–23.
Zheng X-F, Liu K-X, Wang X-M, Zhang R, Li X. Microrna-192 acts as a tumor suppressor in colon cancer and simvastatin activates mir-192 to inhibit cancer cell growth. Mol Med Rep. 2019; 19(3):1753–60.
Group BDW, Atkinson Jr AJ, Colburn WA, DeGruttola VG, DeMets DL, Downing GJ, Hoth DF, Oates JA, Peck CC, Schooley RT, et al.Biomarkers and surrogate endpoints: preferred definitions and conceptual framework. Clin Pharmacol Ther. 2001; 69(3):89–95.
Bolha L, Ravnik-Glavač M, Glavač D. Long noncoding rnas as biomarkers in cancer. Dis Markers. 2017; 2017. https://doi.org/10.1155/2017/7243968.
Sartori DA, Chan DW. Biomarkers in prostate cancer: what’s new?Curr Opin Oncol. 2014; 26(3):259.
Shao Y, Ye M, Jiang X, Sun W, Ding X, Liu Z, Ye G, Zhang X, Xiao B, Guo J. Gastric juice long noncoding rna used as a tumor marker for screening gastric cancer. Cancer. 2014; 120(21):3320–8.
Silva A, Bullock M, Calin G. The clinical relevance of long non-coding rnas in cancer. Cancers. 2015; 7(4):2169–82.
Weng L-P, Smith WM, Dahia PL, Ziebold U, Gil E, Lees JA, Eng C. Pten suppresses breast cancer cell growth by phosphatase activity-dependent g1 arrest followed by cell death. Cancer Res. 1999; 59(22):5808–14.
Depowski PL, Rosenthal SI, Ross JS. Loss of expression of the pten gene protein product is associated with poor outcome in breast cancer. Mod Pathol. 2001; 14(7):672.
Bose S, Crane A, Hibshoosh H, Mansukhani M, Sandweis L, Parsons R. Reduced expression of pten correlates with breast cancer progression. Hum Pathol. 2002; 33(4):405–9.
Chou C-H, Shrestha S, Yang C-D, Chang N-W, Lin Y-L, Liao K-W, Huang W-C, Sun T-H, Tu S-J, Lee W-H, Others. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2017; 46(D1):296–302.
Tomczak K, Czerwińska P, Wiznerowicz M. The Cancer Genome Atlas (TCGA): an immeasurable source of knowledge. Contemp Oncol. 2015; 1A:68–77. https://doi.org/10.5114/wo.2014.47136.
The results shown here are in whole or part based upon data generated by the TCGA Research Network: https://www.cancer.gov/tcga.
About this supplement
This article has been published as part of Volume 21, Supplement 8 2020: Italian Society of Bioinformatics (BITS): Annual Meeting 2019. The full contents of the supplement are available at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume-21-supplement-8.
The publication cost of this paper was funded by the Bioinformatics and Computational Biology for Precision Medicine group (BCB4PM) belonging to the Institute of High Performance Computing and Networking of National Research Council of Italy (ICAR-CNR).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have not competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1
Help pages for miRTissue ce main functionalities.
Additional file 2
Full list of enriched genes according to BP, MF, CC categories, and pathway analysis.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Fiannaca, A., Paglia, L.L., Rosa, M.L. et al. miRTissue ce: extending miRTissue web service with the analysis of ceRNA-ceRNA interactions. BMC Bioinformatics 21 (Suppl 8), 199 (2020). https://doi.org/10.1186/s12859-020-3520-z
- ceRNA interaction
- miRNA-target interaction
- miRNA sponge