Skip to main content

miRTrail - a comprehensive webserver for analyzing gene and miRNA patterns to enhance the understanding of regulatory mechanisms in diseases

Abstract

Background

Expression profiling provides new insights into regulatory and metabolic processes and in particular into pathogenic mechanisms associated with diseases. Besides genes, non-coding transcripts as microRNAs (miRNAs) gained increasing relevance in the last decade. To understand the regulatory processes of miRNAs on genes, integrative computer-aided approaches are essential, especially in the light of complex human diseases as cancer.

Results

Here, we present miRTrail, an integrative tool that allows for performing comprehensive analyses of interactions of genes and miRNAs based on expression profiles. The integrated analysis of mRNA and miRNA data should generate more robust and reliable results on deregulated pathogenic processes and may also offer novel insights into the regulatory interactions between miRNAs and genes. Our web-server excels in carrying out gene sets analysis, analysis of miRNA sets as well as the combination of both in a systems biology approach. To this end, miRTrail integrates information on 20.000 genes, almost 1.000 miRNAs, and roughly 280.000 putative interactions, for Homo sapiens and accordingly for Mus musculus and Danio rerio. The well-established, classical Chi-squared test is one of the central techniques of our tool for the joint consideration of miRNAs and their targets. For interactively visualizing obtained results, it relies on the network analyzers and viewers BiNA or Cytoscape-web, also enabling direct access to relevant literature. We demonstrated the potential of miRTrail by applying our tool to mRNA and miRNA data of malignant melanoma. MiRTrail identified several deregulated miRNAs that target deregulated mRNAs including miRNAs hsa-miR-23b and hsa-miR-223, which target the highest numbers of deregulated mRNAs and regulate the pathway "basal cell carcinoma". In addition, both miRNAs target genes like PTCH1 and RASA1 that are involved in many oncogenic processes.

Conclusions

The application on melanoma samples demonstrates that the miRTrail platform may open avenues for investigating the regulatory interactions between genes and miRNAs for a wide range of human diseases. Moreover, miRTrail cannot only be applied to microarray based expression profiles, but also to NGS-based transcriptomic data. The program is freely available as web-server at mirtrail.bioinf.uni-sb.de.

Background

Gene expression profiles have gained increasing relevance over the last three decades and have become essential in modern biomedical sciences. About two decades ago, a further class of RNAs has been discovered: these non-coding oligonucleotides are indeed transcribed from the human genome, but no proteins are assembled according to their blueprints. MicroRNAs are a subgroup of these non-coding RNAs, currently attracting more and more attention. They have first been reported in a work by Ruvkun [1] and their first appearance in experiments has been associated with Lee et al [2].

MicroRNAs usually consist of 17 to 23 nucleotides and are detectable in the majority of human tissues and almost all bodily fluids [3–5]. It is known today that microRNAs influence the expression of target genes by binding to the corresponding mRNA, leading to its inactivation. Over 50% of all human coding genes seem to be targets of these short non-coding RNAs. MicroRNAs hereby help to control and fine-tune physiological cellular processes like differentiation, proliferation, or apoptosis. Nowadays, it also became apparent that microRNAs have a strong impact on pathological processes as well: Various microRNAs show altered expression patterns in human disorders including malignant [6–10], neurological [11], cardiovascular [12, 13], or rheumatic diseases [14, 15]. In order to get new insights into the molecular mechanisms leading to a specific disease, increasing attention is paid to the interaction of microRNAs and mRNAs of target genes.

The technologies that are most commonly applied to measure miRNA expression profiles are closely related to the methods for measuring gene expression profiles, namely quantitative real-time polymerase chain reaction (qRT-PCR) [16, 17], oligonucleotide microarrays [18, 19], and high-throughput sequencing [20, 21]. These three technologies allow measuring the expression of sets of miRNA very efficiently. While qRT-PCR is mostly applied to rather small subsets of miRNAs, microarrays enable to profile the whole human miRNome and high-throughput sequencing is additionally applied to detect novel mature forms of miRNAs. Remarkably, with the still growing number of miRNAs, and the likewise growing number of biological experiments carried out with the above-mentioned high-throughput methods, and the manifold of possible interactions between miRNAs and mRNAs, computer aided analyses are essential to grasp the information hidden in the large data sets. Therefore, much ongoing work focuses on the combined analysis of miRNAs and their targets. Two classes of bioinformatics approaches related to this topic are 1) tools that aim at discovering the targets of miRNAs and 2) tools that aim at an integrative analysis of miRNA and mRNA sets. Algorithms belonging to the first class usually rely on sequence-complementarity and often also include thermodynamical aspects [22], machine learning [23–25], or experimental validation steps [26]. An overview of respective programs, including a comparison, can be found in [27]. Additionally, approaches primarily based on experiments are becoming prominent in recent years [28, 29]. Naturally, these approaches are more likely to reveal significant miRNA - mRNA interaction pairs than computational approaches. However, they usually require not unimportant amounts of time and resources and, e.g. by design, might also miss relevant interactions. While not strictly being a tool for the discovery of targets or for an integrated analysis of miRNA and mRNA sets, TAM [30] offers enrichment analyses on miRNA sets, thus potentially paving the way to link common functions with related miRNAs. Tools for the second purpose, an integrative analysis of genes and their miRNA regulators, include MMIA [31], DIANA-mirExTra [32], or miRGator [33]. MMIA, allows to combine expression profiles of miRNA and mRNA experiments and then performs a pathway analysis on the intersection of the predicted target mRNAs and the according inversely correlated mRNAs. Additional analyses include Transcription Factor Binding Sites enrichment and diseases that are found to be associated with the inversely deregulated miRNAs. DIANA-mirExTra web-server integrates the potentially novel prediction of miRNAs having one or more of the submitted genes as their targets. This, likewise, allows shedding light on the function of the miRNAs. In detail, the algorithm investigates the 3' UTR sequences of deregulated genes and searches for over-represented six nucleotide long motifs, thus, enabling the identification of matching miRNAs. Finally, miRGator uses public expression data to analyze expression correlation between miRNA and target mRNA/proteins. The miRNA - target interactions are based on miRanda [22], PicTar [34], and TargetScanS [35] and the function of miRNAs is inferred from the related target mRNAs. To this end, a statistical enrichment analysis is performed for the established GO-terms, pathways, and also disease associations. Moreover, it integrates a first approach towards a manual inspection of the underlying network, offering vertex- or edge-filtering but, to-date, no ways to further cope with this information. Here, we present miRTrail (freely available to non-commercial users at mirtrail.bioinf.uni-sb.de), a knowledge-based tool for integrative network analysis that allows for studying the interactions between microRNAs and their target genes, and especially in the case of diseases, the implications of expression changes on pathogenic processes. Our tool excels by its broad functionality, as (1) it can be applied to a single disease or a group of diseases, (2) it covers a wide variety of biochemical categories, and it can be used to evaluate (3) qRT-PCR, microarray, as well as NGS-based transcriptome data. In its current stage, the organisms of Homo sapiens, Mus musculus, and Danio rerio are supported and further extension is continuing. While many solutions exist that provide either analyses of miRNAs, or mRNAs, or a combination of both, miRTrail allows for the simultaneous, combined statistical analysis of all of these three. A schematic description of its workflow is presented in Figure 1, depicting the integration of the provided data about miRNA and mRNA deregulation and the offered statistical analyses intended to facilitate the work with such complex information, especially when used in combination. As such, our tool is able to not only give initial but also thorough insights, even for a very detailed inspection of the given input based on the network analysis.

Figure 1
figure 1

Workflow. Workflow of miRTrail. User submits two RNA sets (one is the set of deregulated miRNAs, the other is the set of deregulated mRNAs, both for the same disease). Orange color represents information flow of miRNA-related information: For each provided miRNA, the target mRNAs are determined (based on microCosm predictions or custom, uploaded interactions). This information then is used by miRTrail, indicated by the red arrow. In general, red color represents flow of mRNA-related information: The uploaded mRNA set as well as the miRNA targets are used in GeneTrail to perform ORAs as described in the "Methods"-section. Blue color represents information flow of results-related information, e.g. for the overlap of pathway sets. Finally in the results, the network analysis allows for targeted inspection of the provided information, based e.g. on miRNA families-related subnetworks. The modular design of miRTrail becomes visible here, also allowing for convenient extension of future analyses and usage for a diversity of different organisms.

One of the original goals of our research was to improve the understanding on the molecular level of melanoma. Thus, as a first application, we investigated miRNA and mRNA expression profiles of this cancer entity, integrating information from 1) the gene expression omnibus GEO [36], 2) the PhenomiR 2.0 human miRNA and diseases database [37], 3) target prediction algorithms [38], and 4) biochemical pathway information of different resources integrated via the BNDB and GeneTrail [39]. By applying miRTrail to these data, we found highly significant coherences between dysregulated miRNAs and matching dysregulated targets of these miRNAs. An additional network analysis highlighted the potential implications of eight miRNAs via their target mRNAs on pathogenic processes in melanomas.

Implementation

In this section, we start by describing the general idea behind miRTrail, followed by the data, techniques and tools that are used to provide the rich functionality, as well as information on the exemplary input data. Here, the input is originating from publicly available services like NCBI GEO and PhenomiR. Our tool is not restricted to these services, as they are intended for demonstration purposes. Especially, all of miRTrail's functionality is available for the organisms of Homo sapiens, Mus musculus, and Danio rerio, and can easily be extended to support other organisms in the future.

Methodoloy - Multipartite graph

Our webservice miRTrail allows for the joint/integrated analysis of miRNA and mRNA entities, in respect to given diseases - since the latter protein-coding RNAs are targets of the former non-coding RNAs. We decided to realize the integration of data by constructing a graph or network. To be exact, an r-partite or multipartite graph G:

G = V , E

with:

  • V = VmiR ∪ VmR ∪ Vdi

  • V X ∩ V Y = ∅ for X ≠ Y and X, Y ∈ {miR, mR, di}

  • E = EmR-di ∪ EmiR-di ∪ EmiR-mR

  • E X − Y = {(u, v)| u ∈ V X , v ∈ V Y , X ≠ Y} for X, Y ∈ {miR, mR, di}

where the vertex-set V is the union of disjoint vertex-subsets, the edge-set E is the union of disjoint edge-subsets, and an edge only connects vertices of different vertex-subsets.

For each user and the according uploads, individual networks are created. An efficient and open-source interface for the creation of graphs is offered by the C++ Boost Graph Library and the associated adjacency_list-construct.

MicroRNA - Target mRNA Interactions

The actually known 20,000 genes and 1,000 miRNAs allow for 20 million possible interaction pairs, where a miRNA may regulate a gene. To find the most reliable candidates, prediction algorithms have been developed. One of the most prominent algorithms for miRNA - target mRNA interactions is the miRanda algorithm and the respective web-resource microCosm [38]. The miRanda algorithm is sequence-complementarity based and includes a thermodynamic analysis of the miRNA - target mRNA complex. The results are then post-processed by a filtering on conservation of the target site. MicroCosm offers miRNA - target mRNA interactions in combination with a p-value threshold. In the beginning, we decided to perform analyses for three thresholds (0.01, 0.001, 0.0001), and extracted all interactions having a value smaller than the respective alpha level, yielding 279,225, 85,050, and 26,984 interactions, respectively. Because of the heterogenity of the expression data, we finally chose to use a threshold of 0.01, in turn leading to approximately 400 target mRNAs per miRNA in human. The appropriate predictions for the other supported organisms are automatically selected by miRTrail according to the organism in the identifiers of the uploaded miRNA deregulation information.

Alternatively, custom pairwise miRNA - target mRNA interactions can be uploaded in a tab-delimited format, thus allowing e.g. for the use of experimentally validated interactions. Details on the exact format for this input can be found on the homepage of miRTrail, especially regarding gene and miRNA identifiers.

Analysis of independance

Based on the pairwise miRNA - target mRNA interactions (for a custom prediction-threshold (default of 0.01) or from a custom list provided as upload), the miRNA and (target) mRNA of each pair is compared to the input in order to see if it is up-, down- or not deregulated. This information is tabulated in a contingency table to provide an overview to the user. An according p-value is calculated, based on a Chisq-distribution with 4 (6 - 2) degrees of freedom. Given that the uploaded information about dysregulated genes/mRNAs only contains entries with the same direction of deregulation, the computation of an according p-value is not allowed by miRTrail and no p-value will be displayed, but instead a note for the user. However, the table will nevertheless be displayed as an overview.

To help the researcher get an impression about the influence of the used miRNA - target mRNA interactions in this step, especially when using data from prediction algorithms, we offer an option to randomize upon the provided data of miRNA and mRNA deregulation. The deregulation pattern (genes/miRs being up- or downregulated) is kept as-is while the identifiers are sampled at random. This functionality is available via the "Randomize"-button.

MicroRNAs from PhenomiR

For integrating dysregulated miRNAs, we used PhenomiR 2.0 (last update: 2011-02-15). This service offers manually curated data about differential regulation for a variety of diseases.

Specifically, we used the data for entry/ID: 639, concluding the results of a published melanoma study based on microRNA low density arrays, including 666 microRNAs. Selection of the statistically significantly dysregulated miRNAs in the miRNA extracts of adult melanoma patients and benign nevi controls was done with univariate Two-sample T-test and a significance level of 0.05. The size of the patient samples is 10 and 4 for the control, respectively. An overview of the numbers of up- and downregulated miRNAs can be found in Table 1.

Table 1 Summary of deregulated genes and miRNAs

MessengerRNAs from NCBI's GEO

The NCBI Gene Expression Omnibus (GEO) is a public repository providing data from microarray experiments, next-generation sequencing, and other high-throughput functional genomic data. The microarray experiments, in particular, must comply with MIAME guidelines in order to be accepted by NCBI's GEO.

We extracted the microarray expression profiles from data set GDS1375 (series published: 2005-08-25), including 63 arrays for 45 melanoma and 18 benign nevi samples. Due to possible variations between the experiments, we carried out a quantile-normalization of the expression values for all genes present on the respective data set. Selection of differentially expressed genes was performed on the normalized data using the univariate Two-sample T-test and a significance level of 0.05. An overview of the numbers of up- and downregulated genes can be found in Table 1.

GeneTrail

The gene set analysis tool GeneTrail has been developed to help in the analysis of readily available or newly created high-throughput data. It allows for a comprehensive and efficient statistical evaluation of large genomic or proteomic datasets and covers a plethora of biological categories and pathways, e.g. KEGG, TRANSPATH, TRANSFAC, and GO. Analyses can be either performed via an 'Over-Representation Analysis' (ORA) comparing a reference set of genes to a test set or a 'Gene Set Enrichment Analysis' (GSEA) based on a sorted list of genes. While the calculation of ORA p-values relies on Hypergeometric distribution, many existing tools offer the calculation of GSEA p-values based on permutation tests, usually limited to a fixed number of permutations for performance reasons. GeneTrail integrates an exact calculation [40] corresponding to a commonly used non-parametric unweighted permutation test. This calculation is based on dynamic programming and thus allows, especially for large sets, a higher accuracy than by using a fixed number of permutations.

Recently, GeneTrail has been extended to directly allow the analysis of expression data originating from the NCBI GEO, resulting in GeneTrail Express [41]. This integration greatly facilitates the selection of differentially regulated genes and allows for a fast evaluation of the expression profiles in respect to biological categories and pathways.

Visualization: BiNA and Cytoscape-web

While computational approaches are very important in contemporary research, manual inspection is often desireable to support the automatic analyses or to identify new aspects. To this end, we decided to include the visualization of the resulting interaction network of miRNAs and their (putative) targets. Due to the large amount of integrated data, efficient means for focusing are crucial. Therefore, we provide respective subnetworks, depending either on the choice of individual miRNAs or on members of miRNA-families contained in the input. Furthermore, only deregulated miRNAs are respected that are connected to deregulated target mRNAs, either by the prediction algorithm or the provided custom interactions, as we envision these entities and relations as the most relevant. For the actual visualisation, two selections are available for the user: BiNA and Cytoscape-web.

BiNA is a visualization and analysis tool for various biological networks. We developed a plug-in for the Java Webstart version of BiNA, which takes the miRTrail results and uses the visualization capabilities of BiNA for presenting the network. The user can choose between different graph layouts (organic, hierarchic, and orthogonal) and can modify the visualization in many ways. By default, the target-mRNA nodes are sized according to their degree for easier retrieval of high-degree nodes. It is also possible to save the network in different file formats for reusing the data in other tools or BiNA again. For larger graphs, this visualization-option is probably beneficial.

Cytoscape-web is modeled after the Cytoscape Java network visualization and analysis software [42]. Its JavaScript API allows for an integration into HTML-pages and convenient display of networks. We offer the user a choice of three different graph layouts (Circular, Radial, Tree) and the possibility to select the first neighbors of a selected node. Zoom and pan functionality is available and target-nodes are also sized according to their degree. As the graph is directly displayed in the browser-window, this visualization is especially suitable for a quick inspection of the network. Finally, we implemented context-menu items that greatly facilitate the search for related publications by performing NCBI PubMed queries ("inclusive" or "exclusive") for a custom selection of miRNA and mRNA nodes, given a disease was specified in the input.

Results and Discussion

In the following, we will describe the range of different functions offered by miRTrail. Subsequently, an analysis of cutaneous malignant melanoma versus benign nevi is performed to illustrate the potential of our tool.

Functionality of miRTrail

The miRTrail webserver recieves two dysregulation sets in separate text-files as input, one is the set of dysregulated miRNAs and the other is the set of dysregulated genes. For each uploaded identifier (for miRNAs, the standard annotation of miRBase is used, for genes the HGNC GeneSymbol annotation, respectively), the information whether the respective gene/miRNA is upregulated ('1') or down-regulated ('−1') has to be provided in the files by the user. Optionally, the disease of interest can be specified to allow for convenient NCBI PubMed queries for related information.

As the next step, the user can either choose a target-prediction threshold for microCosm targets predictions or can provide a list of custom pairwise miRNA - target mRNA interactions, potentially originating from proprietary experiments or other prediction algorithms. The default threshold for microCosm targets predictions is 0.01, amounting to around 280,000 miRNA - target mRNA interactions. Here, the user can also opt-in for a thorough GeneTrail analysis. Based on this information, the analyses are then carried out and, finally, the user is directed to the results. These will be stored uniquely for each analysis performed and can be shared with others by simply providing them with the link of the results page. The results presented herein can be reproduced using the example files provided by miRTrail.

The first provided analysis computes a contingency table relating the dysregulation of miRNAs and the dysregulation of target mRNAs and calculates the according p-value, based on a χ2 distribution. This analysis allows for estimating whether there is an independance in the deregulation of the miRNAs and the target mRNAs.

Second, a Venn diagram is computed, providing the dysregulated genes that are targets of dysregulated miRNAs (overlap of the diagram), the not-dysregulated targets of dysregulated miRNAs (left part of the diagram) and the dysregulated genes that are not targets of the dysregulated miRNAs (right part of the diagram). For this Venn diagram, a p-value using the Hypergeometric distribution is calculated to show whether there exists a significant overlap between dysregulated genes and targets of dysregulated miRNAs. Third, gene set enrichment analyses for three gene sets are carried out using the comprehensive functionality of GeneTrail. Independently of each other, a so-called Over-Representation Analysis (ORA) - based on the Hypergeometric distribution - is carried out for the dysregulated genes, targets of dysregulated miRNAs and dysregulated targets of dysregulated miRNAs. In all cases, the gene sets are tested for significant enrichments/depletions in KEGG pathways. If the user previously decided to perform all GeneTrail analyses, the results will also include information about GO terms, TransPath pathways, transcription factors from Transfac, SNPs, and chromosomal location, among many others. By clicking on the 'details' button, the complete list of results is provided. Moreover, an overview showing the biological categories being significant in at least two of the three analyses is created. The "code" represents in which of the pairwise overlaps the respective category was found, similar to the file-permission scheme in Linux. So, a code of "2" e.g. shows that a category was found in the enrichment analysis of the dysregulated targets of dysregulated miRNAs as well as in the results of the dysregulated genes. A code of "3" would hence mean that this category was additionally found in the results of the targets of dysregulated miRNAs. Accordingly in Table 2, e.g. the "DNA replication" pathway was found to be enriched for dysregulated targets of dysregulated miRNAs as well as for the dysregulated genes/mRNAs.

Table 2 Overlapping pathways

Finally, we carry out an integrative network-analysis approach on the comprehensive network containing dysregulated genes, dysregulated miRNAs, and the target interactions between them. A subset of interesting miRNAs and their according targets is selected as well as a custom degree constraint. The subset can be constructed either by selecting individual miRNAs, s. Figure 2, or miRNA families based on an Over-Representation Analysis of miRNA-family data from miRBase (miFam.dat) [43], s. Figure 3. The custom degree constraint allows the selection of the genes being the target of at least as many miRNAs as specified by the parameter. Based on this selection, using the Java Webstart-based viewer BiNA [44] or the web-based viewer Cytoscape-web [45], we show the resulting network, allowing for a manual inspection of the inherent interactions. In the network visualizations, nodes with rectangular shapes belong to miRNAs, nodes with round shapes to genes, red color means up-regulation, green color means down-regulation, and genes and miRNAs are connected by edges if a putative miRNA - target interaction exists. Additionally, a more fine-grained ORA is available, being performed only on the genes that are contained in the custom selection, which is also separately available as a list. These genes are assumed to be the most disease-relevant as they are found to be deregulated and simultaneously putative targets of deregulated miRNAs while, at the same time, being central to the network, according to their degree. This list is analysed for enrichments/depletions in KEGG pathways, Gene Ontology terms, OMIM disease relations, and NIA human disease gene sets. A thorough GeneTrail analysis can also be chosen here.

Figure 2
figure 2

MiRNA selection (individual). Demonstrates the selection of miRNAs of interest. Link next to the each miRNA shows the respective dysregulated target mRNAs.

Figure 3
figure 3

MiRNA selection (families). Demonstrates the selection of miRNAs based on enriched miRNA families. All significant families (p(adj) < 0.05) are preselected.

Melanoma case study

We compared the expression profiles of cutaneous malignant melanoma to those of benign skin nevi samples from adult patients. While the proportion of melanoma cases among skin cancer patients is rather low (4%), it accounts for almost 75% of all skin cancer-related deaths. Even more, the prognosis for advanced melanoma is very poor (5-year survival-rate is only 5%) [46]. Hence, we decided to validate our tool based on melanoma data and to identify new aspects of this disease, potentially helping in the creation of promising new therapies for advanced melanoma patients.

An illustration of the the results page for the miRTrail analysis on the melanoma miRNA and mRNA samples, as mentioned in the previous section, can be seen in Figure 4.

Figure 4
figure 4

Results Page. This illustrates the results page of a miRTrail analysis. The network visualization and subsequent analyses are available in the subpanel at the bottom.

Analysis of independance

Using our tool, for the melanoma samples and a prediction threshold of 0.01 for the human miRNA - target mRNA interactions from MicroCosm targets, we were able to find statistical evidence about the dependance of deregulation of miRNAs and target mRNAs. The contingency table yielded a p-value of 0.025 (α = 0.05). Interestingly, independant of the miRNA being up- or downregulated, similar amounts of interactions were found for targets, then, being up (450 and 480), down (394 and 355), or not deregulated (3193 and 3385), respectively. In turn, more targets were found to be up-regulated for dysregulated miRNAs as well as for not dysregulated miRNAs. Surprisingly, 805 (450 + 355) interactions were found with both miRNA and predicted target, being deregulated in the same direction. Finally, as expected, a large number of interactions was found were both, the miRNA and the predicted target, were not deregulated (212782).

ORA of the three mRNA sets

Inspecting the ORA results, showed one KEGG pathway (GPI-anchor biosynthesis) to be significantly enriched for the set of targets of dysregulated miRNAs. In turn, the over-representation analysis of dysregulated targets of dysregulated miRNAs revealed nine significant pathways, including cancer related categories, like Non-small cell lung cancer, Prostate cancer, Small cell lung cancer, Endometrial cancer, and Glioma as well as enrichments in the Lysosome pathway and DNA replication. The analysis of the dysregulated mRNAs revealed the highest number of statistically signifi cant KEGG categories with a total of 24. Here again, several cancer-related pathways were found to be enriched as well as pathways like Cell cycle, Focal adhesion, or even signaling pathways (e.g. TGF-beta signaling pathway).

The result of the pairwise overlaps of the resulting pathway sets is described in Table 2. Among those, the DNA replication pathway [47] as well as the Lysosome pathway [48, 49] have already been attributed to melanoma.

Network analysis

From the 33 input miRNAs, 21 were found to have dysregulated targets for a prediction threshold of 0.01. The miRNA with the least dysregulated targets was hsa-miR-211 (52 targets) while hsa-miR-23b was the miRNA with the most dysregulated targets (115). For this analysis, we decided to use the eight miRNAs having more than 80 dysregulated targets (miR-23b [50], miR-223, miR-193b [51], miR-424, miR-20a [52], miR-98, miR-891a, and miR-566), see Figure 2. We left the custom degree constraint at the default of 1 for the subsequent ORA. The resulting mRNA set comprises the dysregulated mRNAs that were predicted targets of at least one of the eight earlier miRNAs. Specifying a higher constraint would lead to a smaller network with only the mRNA nodes being targets of at least as many miRNAs as specified by this parameter and the according miRNA nodes, respectively.

ORA of subnetwork

The analysis of KEGG pathways showed significant enrichments for the three cancer-related categories: Prostate cancer, Non-small cell lung cancer, and Endometrial cancer. These categories were found to be enriched for genes that were deregulated while being targets of deregulated miRNAs, hence, the genes that are assumed to be the most disease-relevant due to their joint deregulation.

A total of 274 GO terms were found to be enriched or depleted for all of the three GO-trees, with enrichments in anti-apoptosis, cell proliferation, cell cycle, transcript initiation, RNA elongation, and regulation of translational initiation among others in the biological subtree.

Furthermore, an enrichment (RASA1 and PTCH1) for "Susceptibility to basal cell carcinoma" was found in the OMIM categories.

Visualization

For this step, we decided to focus on smaller miRNA and mRNA sets to increase the visibility. However, also large selections can be efficiently handled and used for detailed manual inspections. Exemplary visualization can be found in Figure 5 and 6, for BiNA and for Cytoscape-web, respectively.

Figure 5
figure 5

Visualization: BiNA. Subnetwork of top-8 miRNAs and a degree constraint for the target mRNAs of 3, thus, only seven miRNA nodes are displayed. Round shape represents mRNAs, rectangular miRNAs, respectively. Red color indicates upregulation, green color downregulation, respectively. Size of the mRNAs is according to their degree.

Figure 6
figure 6

Visualization: Cytoscape-web. Subnetwork of top-8 miRNAs and a degree constraint for the target mRNAs of 3, thus, only seven miRNA nodes are displayed. Round shape represents mRNAs, rectangular miRNAs, respectively. Red color indicates upregulation, green color downregulation, respectively. Size of the mRNAs is according to their degree.

Conclusions

The constantly increasing availability of data from different origins and of different nature allows for more complex and comprehensive analyses. To this end, we developed miRTrail to integrate information about RNA deregulation in diseases and putative interactions of miRNAs and mRNAs. Our tool provides a large collection of different analyses and performs in a transparent way, requiring only minor activity by the user, while offering versatile results. This greatly facilitates the adaption of this tool as it does not require complicated initial learning. Via the visualisation component, miRTrail enables the user to easily inspect the interactions and, thus, also to further process upon the selection.

MiRTrail - results will also be of great help in any scheme that aims in experimental confirmation of miRNA-targets. The final proof, here, requires extended experiments including the identification of the specific targeted region of a gene by in vitro binding and the analysis of in vivo effects by altered miRNA expression. The melanoma case study shows that we were able to detect highly significant results, despite the fact that we did not use autologous samples. This sets the ground for specific experimental assays that focus on significant miRNA - mRNA interactions in this tumor type. Hence, miRTrail is of great interest for the life sciences community as it can use data from next-generation sequencing, qRT-PCR, or microarray experiments.

Availability and requirements

Project name: miRTrail

Project home page: http://mirtrail.bioinf.uni-sb.de

Operating system(s): Platform independent

Programming languages: C++, php

Other requirements: JavaWS version 1.6 or higher

References

  1. Ruvkun G: Molecular biology. Glimpses of a tiny RNA world. Science 2001, 294(5543):797–9. 10.1126/science.1066315

    Article  CAS  PubMed  Google Scholar 

  2. Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell 1993, 75(5):843–54. 10.1016/0092-8674(93)90529-Y

    Article  CAS  PubMed  Google Scholar 

  3. Zen K, Zhang CY: Circulating MicroRNAs: a novel class of biomarkers to diagnose and monitor human cancers. Med Res Rev 2010.

    Google Scholar 

  4. Iguchi H, Kosaka N, Ochiya T: Secretory microRNAs as a versatile communication tool. Commun Integr Biol 2010, 3(5):478–81. 10.4161/cib.3.5.12693

    Article  PubMed Central  PubMed  Google Scholar 

  5. Weber JA, Baxter DH, Zhang S, Huang DY, Huang KH, Lee MJ, Galas DJ, Wang K: The microRNA spectrum in 12 body fluids. Clin Chem 2010, 56(11):1733–41. 10.1373/clinchem.2010.147405

    Article  CAS  PubMed  Google Scholar 

  6. Medina PP, Slack FJ: MicroRNAs and cancer: an overview. Cell Cycle 2008, 7: 2485–2492. 10.4161/cc.7.16.6453

    Article  CAS  PubMed  Google Scholar 

  7. Zhang B, Pan X, Cobb GP, Anderson TA: MicroRNAs as oncogenes and tumor suppressors. Dev Biol 2007, 302: 1–12. 10.1016/j.ydbio.2006.08.028

    Article  CAS  PubMed  Google Scholar 

  8. Roth P, Wischhusen J, Happold C, Chandran A, Hofer S, Eisele G, Weller M, Keller A: A specific miRNA signature in the peripheral blood of glioblastoma patients. J Neurochem 2011.

    Google Scholar 

  9. Keller A, Leidinger P, Borries A, Wendschlag A, Wucherpfennig F, Scheffler M, Huwer H, Lenhof HP, Meese E: miRNAs in lung cancer: Studying complex fingngerprints in patient's blood cells by microarray experiments. BMC Cancer 2009.

    Google Scholar 

  10. Häusler S, Keller A, Chandran A, Ziegler K, Zipp K, Heuer S, Krockenberger M, Engel J, Hönig A, Scheffler M, Dietl J, Wischhusen J: Whole blood-derived miRNA profiles as potential new tools for ovarian cancer screening. Br J Cancer 2010.

    Google Scholar 

  11. Saugstad JA: MicroRNAs as effectors of brain function with roles in ischemia and injury, neuroprotection, and neurodegeneration. J Cereb Blood Flow Metab 2010, 30(9):1564–76. 10.1038/jcbfm.2010.101

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Frost RJ, van Rooij E: miRNAs as therapeutic targets in ischemic heart disease. J Cardiovasc Transl Res 2010, 3(3):280–9. 10.1007/s12265-010-9173-y

    Article  PubMed  Google Scholar 

  13. Meder B, Keller A, Vogel B, Sedaghat F, Kayvanpour E, Haas J, Just S, Borries A, Rudloff J, Leidinger P, Meese E, Katus H, Rottbauer W: MicroRNA Signatures as Novel Biomarkers for Acute Myocardial Infarction. Basic Res Cardiol 2010.

    Google Scholar 

  14. Alevizos I, Illei GG: MicroRNAs as biomarkers in rheumatic diseases. Nat Rev Rheumatol 2010, 6(7):391–8. 10.1038/nrrheum.2010.81

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Keller A, Leidinger P, Lange A, Borries J, Schroers H, Scheffler M, Lenhof HP, Ruprecht K, Meese E: Multiple sclerosis: MicroRNA expression profiles accurately differentiate patients with relapsing-remitting disease from healthy controls. PLoS One 2009.

    Google Scholar 

  16. Fiedler SD, Carletti MZ, Christenson LK: Quantitative RT-PCR Methods for Mature microRNA Expression Analysis. Methods Mol Biol 2010, 630: 49–64. 10.1007/978-1-60761-629-0_4

    Article  CAS  PubMed  Google Scholar 

  17. Chen C, Ridzon DA, Broomer A, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, Lao KQ, Livak KJ, Guegler KJ: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res 2005, 33(20):e179. 10.1093/nar/gni178

    Article  PubMed Central  PubMed  Google Scholar 

  18. Thomson JM, Parker J, Perou CM, Hammond SM: A custom microarray platform for analysis of microRNA gene expression. Nat Methods 2004, 1: 47–53. 10.1038/nmeth704

    Article  CAS  PubMed  Google Scholar 

  19. Miska1 EA, Alvarez-Saavedra E, Townsend M, Yoshii A, Å estan N, Rakic P, Constantine-Paton M, Horvitz HR: Microarray analysis of microRNA expression in the developing mammalian brain. Genome Biol 2004, 5: R68. 10.1186/gb-2004-5-9-r68

    Article  PubMed Central  PubMed  Google Scholar 

  20. Motameny S, Wolters S, NÄurnberg P, Schumacher B: Next Generation Sequencing of miRNAs - Strategies, Resources and Methods. Genes 2010, 1(1):70–84. 10.3390/genes1010070

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Morin RD, O'Connor MD, Griffith M, Kuchenbauer F, Delaney A, Prabhu AL, Zhao Y, McDonald H, Zeng T, Hirst M, Eaves CJ, Marra MA: Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Res 2008, 18: 610–621. 10.1101/gr.7179508

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human MicroRNA targets. PLoS Biol 2004, 2: e363. 10.1371/journal.pbio.0020363

    Article  PubMed Central  PubMed  Google Scholar 

  23. Wang X, El Naqa IM: Prediction of both conserved and nonconserved microRNA targets in animals. Bioinformatics 2008, 24(3):325–332. 10.1093/bioinformatics/btm595

    Article  PubMed  Google Scholar 

  24. Wang X: miRDB: a microRNA target prediction and functional annotation database with a wiki interface. RNA 2008, 14(6):1012–1017. 10.1261/rna.965408

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Huang JC, Babak T, Corson TW, Chua G, Khan S, Gallie BL, Hughes TR, Blencowe BJ, Frey BJ, Morris QD: Using expression profiling data to identify human microRNA targets. Nat Methods 2007, 4: 1045–1049. 10.1038/nmeth1130

    Article  CAS  PubMed  Google Scholar 

  26. Papadopoulos GL, Reczko M, Simossis VA, Sethupathy P, Hatzigeorgiou AG: The database of experimentally supported targets: a functional update of TarBase. Nucleic Acids Res 2009, (37 Database):D155–8.

  27. Thomas M, Lieberman J, Lal A: Desperately seeking microRNA targets. Nat Struct Mol Biol 2010, 1169–74.

    Google Scholar 

  28. Su WL, Kleinhanz RR, Schadt EE: Characterizing the role of miRNAs within gene regulatory networks using integrative genomics techniques. Mol Syst Biol 2011, 7: 490.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Jayaswal V, Lutherborrow M, Ma DD, Yang YH: Identification of microRNA-mRNA modules using microarray data. BMC Genomics 2011, 12: 138. 10.1186/1471-2164-12-138

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Lu M, Shi B, Wang J, Cao Q, Cui Q: TAM: A method for enrichment and depletion analysis of a microRNA category in a list of microRNAs. BMC Bioinformatics 2010, 11: 419. 10.1186/1471-2105-11-419

    Article  PubMed Central  PubMed  Google Scholar 

  31. Nam S, Li M, Choi K, Balch C, Kim S, Nephew KP: MicroRNA and mRNA integrated analysis (MMIA): a web tool for examining biological functions of microRNA expression. Nucleic Acids Res 2009, (37 Web Server):W356-W362.

    Google Scholar 

  32. Alexiou P, Maragkakis M, Papadopoulos GL, Simmosis VA, Zhang L, Hatzigeorgiou AG: The DIANA-mirExTra web server: from gene expression data to microRNA function. PLoS ONE 2010., 5(2):

    Google Scholar 

  33. Nam S, Kim B, Shin S, Lee S: miRGator: an integrated system for functional annotation of microRNAs. Nucleic Acids Res 2008, (36 Database):D159–64.

    Google Scholar 

  34. Krek A, Grün D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade1 I, Gunsalus KC, Stoffel M, Rajewsky N: Combinatorial microRNA target predictions. Nat Genet 2005, 37: 495–500. 10.1038/ng1536

    Article  CAS  PubMed  Google Scholar 

  35. Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 2005, 120(1):15–20. 10.1016/j.cell.2004.12.035

    Article  CAS  PubMed  Google Scholar 

  36. Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Muertter RN, Edgar R: NCBI GEO: archive for high-throughput functional genomic data. Nucleic Acids Res 2009, (37 Database):D5–15.

    Google Scholar 

  37. Ruepp A, Kowarsch A, Schmidl D, Buggenthin F, Brauner B, Dunger I, Fobo G, Frishman G, Montrone C, Theis FJ: PhenomiR: a knowledgebase for microRNA expression in diseases and biological processes. Genome Biol 2010, 11: R6. 10.1186/gb-2010-11-1-r6

    Article  PubMed Central  PubMed  Google Scholar 

  38. Computational Prediction Protocol of EMBL-EBI microCosm targets[http://www.ebi.ac.uk/enright-srv/microcosm/htdocs/targets/v5/info.html]

  39. Backes C, Keller A, Kuentzer J, Kneissl B, Comtesse N, Elnakady YA, Müller R, Meese E, Lenhof HP: GeneTrail - advanced gene set enrichment analysis. Nucleic Acids Res 2007., (Web Server):

    Google Scholar 

  40. Keller A, Backes C, Lenhof HP: Computation of significance scores of unweighted Gene Set Enrichment Analyses. BMC Bioinformatics 2007., 8:

    Google Scholar 

  41. Keller A, Backes C, Al-Awadhi M, Gerasch A, Kuentzer J, Kohlbacher O, Kaufmann M, Lenhof HP: GeneTrailExpress: a web-based pipeline for the statistical evaluation of microarray experiments. BMC Bioinformatics 2008.

    Google Scholar 

  42. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003, 2498–504.

    Google Scholar 

  43. Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res 2008, (36 Database):D154-D158.

    Google Scholar 

  44. Kuentzer J, Blum T, Gerasch A, Backes C, Hildebrandt A, Kaufmann M, Kohlbacher O, Lenhof HP: BN++ - A Biological Information System. J Integr Bioinform 2006, 3(2):34.

    Google Scholar 

  45. Lopes CT, Franz M, Kazi F, Donaldson SL, Morris Q, Bader GD: Cytoscape Web: an interactive web-based network browser. Bioinformatics 2010.

    Google Scholar 

  46. Mueller DW, Bosserhoff AK: Role of miRNAs in the progression of malignant melanoma. Br J Cancer 2009, 101(4):551–6. 10.1038/sj.bjc.6605204

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Kauffmann A, Rosselli F, Lazar V, Winnepenninckx V, Mansuet-Lupo A, Dessen P, van den Oord JJ, Spatz A, Sarasin A: High expression of DNA repair pathways is associated with metastasis in melanoma patients. Oncogene 2008, 27(5):565–73. 10.1038/sj.onc.1210700

    Article  CAS  PubMed  Google Scholar 

  48. Ma XH, Piao S, Wang D, McAfee QW, Nathanson KL, Lum J, Li LZ, Amaravadi RK: Measurements of tumor cell autophagy predict invasiveness, resistance to chemotherapy, and survival in melanoma. Clin Cancer Res 2011, 17(10):3478–89. 10.1158/1078-0432.CCR-10-2372

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Tormo D, Checińska A, Alonso-Curbelo D, Pérez-Guijarro E, Cañón E, Riveiro-Falkenbach E, Calvo TG, Larribere L, Megías D, Mulero F, Piris MA, Dash R, Barral PM, Rodríguez-Peralto JL, Ortiz-Romero P, Tüting T, Fisher PB, Soengas MS: Targeted activation of innate immunity for therapeutic induction of autophagy and apoptosis in melanoma cells. Cancer Cell 2009, 16(2):103–14. 10.1016/j.ccr.2009.07.004

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Philippidou D, Schmitt M, Moser D, Margue C, Nazarov PV, Muller A, Vallar L, Nashan D, Behrmann I, Kreis S: Signatures of microRNAs and selected microRNA target genes in human melanoma. Cancer Res 2010, 70(10):4163–73. 10.1158/0008-5472.CAN-09-4512

    Article  CAS  PubMed  Google Scholar 

  51. Caramuta S, Egyházi S, Rodolfo M, Witten D, Hansson J, Larsson C, Lui WO: MicroRNA expression profiles associated with mutational status and survival in malignant melanoma. J Invest Dermatol 2010, 130(8):2062–70. 10.1038/jid.2010.63

    Article  CAS  PubMed  Google Scholar 

  52. Zhang L, Huang J, Yang N, Greshock J, Megraw MS, Giannakakis A, Liang S, Naylor TL, Barchetti A, Ward MR, Yao G, Medina A, O'brien-Jenkins A, Katsaros D, Hatzigeorgiou A, Gimotty PA, Weber BL, Coukos G: MicroRNAs exhibit high frequency genomic alterations in human cancer. Proc Natl Acad Sci USA 2006, 103(24):9136–41. 10.1073/pnas.0508889103

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work has been funded by DFG Priority Program SPP 1335: LE 952/3- 1, DFG LE 952/5-1, and by DFG Me 917/20-1.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Andreas Keller.

Additional information

Authors' contributions

CL implemented the underlying framework as well as the web-interface. PL and NL performed the melanoma data analysis. AG and MK developed the BiNA visualization component. JH, BV, and BM tested and benchmarked miRTrail and took part in the writing of this manuscript. HAK collaborated in the project and the overall study and contributed to this manuscript. Support in the implementation of GeneTrail-related components came from CB. CS contributed the conceptual design of this study. EM, HPL, and AK developed the overall study and contributed to this manuscript. AK further supported the analysis and implementation of the herein presented solution. EM, HPL, and AK are the senior authors. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Laczny, C., Leidinger, P., Haas, J. et al. miRTrail - a comprehensive webserver for analyzing gene and miRNA patterns to enhance the understanding of regulatory mechanisms in diseases. BMC Bioinformatics 13, 36 (2012). https://doi.org/10.1186/1471-2105-13-36

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-13-36

Keywords