Nothing Special   »   [go: up one dir, main page]
Skip to main content
Download PDF

CyTRANSFINDER: a Cytoscape 3.3 plugin for three-component (TF, gene, miRNA) signal transduction pathway construction

BMC Bioinformatics volume 17, Article number: 157 (2016) Cite this article

Abstract

Background

Biological research increasingly relies on network models to study complex phenomena. Signal Transduction Pathways are molecular circuits that model how cells receive, process, and respond to information from the environment providing snapshots of the overall cell dynamics. Most of the attempts to reconstruct signal transduction pathways are limited to single regulator networks including only genes/proteins. However, networks involving a single type of regulator and neglecting transcriptional and post-transcriptional regulations mediated by transcription factors and microRNAs, respectively, may not fully reveal the complex regulatory mechanisms of a cell. We observed a lack of computational instruments supporting explorative analysis on this type of three-component signal transduction pathways.

Results

We have developed CyTRANSFINDER, a new Cytoscape plugin able to infer three-component signal transduction pathways based on user defined regulatory patterns and including miRNAs, TFs and genes. Since CyTRANSFINDER has been designed to support exploratory analysis, it does not rely on expression data. To show the potential of the plugin we have applied it in a study of two miRNAs that are particularly relevant in human melanoma progression, miR-146a and miR-214.

Conclusions

CyTRANSFINDER supports the reconstruction of small signal transduction pathways among groups of genes. Results obtained from its use in a real case study have been analyzed and validated through both literature data and preliminary wet-lab experiments, showing the potential of this tool when performing exploratory analysis.

Background

Network representation of intracellular biological systems, considering molecular components within a cell as nodes (e.g., genes, proteins, miRNA, etc.) and their direct or indirect interactions as links, is steadily gaining interest because of its potential to represent, characterize, and model a wide range of intricate natural systems and phenomena.

Among the different types of biological network models proposed in the literature [1], we are interested in Signal Transduction Pathways (STP) [2, 3]. A cell is highly responsive to specific chemicals in its environment. Broadly, signal transduction pathways can be viewed as molecular circuits. They model how cells receive, process, and respond to information from the environment toward a biological identified end result, thus providing snapshots of the (overall) cell dynamics. The number of these processes shows how many ways the organism can react and respond to its environment. Therefore, discovering new STPs is an important task to contribute to the current knowledge of the cell behavior.

The traditional approach to identify molecular components of a signaling network is through gene knockout experiments and epistasis analysis [4]. In such experiments, an organism is engineered to suppress the expression of one or more genes in order to study the resulting perturbation in the cell dynamics. Although these experiments are effective to identify simple direct signaling activities, more complex signaling circuitries are difficult to identify and understand. Moreover this analysis is time-consuming, expensive, and sometimes the results can be misinterpreted [5].

Computational approaches for modeling and reconstruction of STPs are currently a hot research area. STPs have been modeled through modular kinetic simulations of biochemical networks [6], and detailed integration of biochemical properties of the pathways [7]. Bayesian networks applied to multi-variate expression data have also been used to infer signaling pathways [8]. More recently, PPI networks have been largely used to reconstruct signaling transduction pathways [9–13]. In general these methods try to extract STPs from PPI networks, which are known to be affected by a high rate of false-positive and false-negative interactions. The use of expression data is used to mitigate this uncertainty.

Most of the attempts to reconstruct STPs focus on gene/protein based networks. However, networks involving a single type of regulator may not fully reveal the complex regulatory mechanisms of a cell. Complexity strongly increases when STPs include post-transcriptional regulation mediated by microRNAs (miRNAs) interacting with different transcription factors (TFs). It is predicted that miRNAs regulate approximately 30 % of the human protein-coding genome [14], they are therefore highly important in modeling the cell regulation. Only a few attempts to reconstruct STPs including miRNAs, TFs, and mRNAs can be found in the literature [15, 16].

Motivated by this, we have developed CyTRANSFINDER, a new Cytoscape 3.3 [17] plugin able to construct three-component signal transduction pathways with the presence of miRNAs, TFs and genes starting from public available regulatory information. Rather than trying to construct big networks as proposed in other studies, CyTRANSFINDER focuses on reconstruction of small signal transduction pathways based on user defined regulatory patterns. These pathways may be of direct use to drive exploratory analysis enabling to better understand experimental data and to further drive laboratory experiments. Formally the problem addressed by CyTRANSFINDER is the following: “Given two set of genes, to discover a set of STPs connecting each gene of the first set with each gene of the second set according to a signaling pattern set by the user.” Recurring signaling patterns have been widely studied in gene regulatory networks as well as other real-world complex systems scenarios [18], because of their central role in driving regulatory responses by specific functions [2]. This assumption is based on the expectation that designs with higher modularity have higher adaptability and therefore higher survival rates [19], thus suggesting that modularity can spontaneously arise under changing environments [20], which eventually results in extremely complex systems made of simple basic building blocks [19].

Since CyTRANSFINDER has been designed to support exploratory analysis, it does not rely on expression data. It includes a data-fusion engine that scrapes information from seven online repositories and integrates them to infer candidate pathways. Different filters can be applied to restrict or enlarge the set of produced results based on the specific use cases. The integration with Cytoscape 3.3 features an intuitive user interface that automates complex tasks and makes the plugin a potential software instrument for biologists with limited skills in computer programming and network analysis. Moreover, it enables to further process and analyze the identified networks with the huge ecosystems of network analysis plugins and functions already available in Cytoscape 3.3.

To the best of our knowledge no other Cytoscape plugin offers the functionalities provided by CyTRANSFINDER. A Cytoscape 2.6 plugin implementing a front-end to BIANA (Biologic Interactions and Network Analysis) is the only tool that somehow offers functionalities related to CyTRANSFINDER [21]. BIANA is a general Python framework aiming at integrating information from several external data-sets in network representations that can be visualized through the Cytoscape plugin. However, differently from CyTRANSFINDER, most of the effort given in BIANA is put on the possibility of describing external data sources and rules to integrate data from different sources. It is therefore a more generic software that does not specifically focuses on the problem of reconstructing STPs, as done instead by CyTRANSFINDER. BIANA standalone application appears discontinued from 2013, while the latest plugin update is dated 2009 and the plugin is only compatible with Cytoscape 2.6, which is becoming obsolete.

To show the capability of the plugin, we have applied it to a study of two miRNAs that are particularly relevant in human melanoma progression, miR-146a and miR-214. Results obtained from CyTRANSFINDER have been analyzed and validated through both literature data and preliminary wet-lab experiments, showing the capability of this tool when performing exploratory analysis.

Implementation

In its basic setup, CyTRANSFINDER implements STP discovery among two sets of genes into Cytoscape integrating regulatory information on the Homo sapiens (human) species. It is developed to work with Cytoscape 3.3. All examples proposed in this paper have been tested with the latest Cytoscape version (Cytoscape 3.3). Once installed from the Cytoscape App Manager, CyTRANSFINDER is available from the Apps menu of Cytoscape. Figure 1 shows a screenshot of CyTRANSFINDER running on a small example whereas Fig. 2 shows the conceptual architecture of the software highlighting its main data sources and computational modules.

Fig. 1
figure 1

CyTRANSFINDER overview. a The main plugin control panel. It allows the user to set the source and destination list of genes and the specific pattern of regulators to search. b Shows the graphical output of the plugin that consists of a network connecting source genes with destination genes. Nodes of this network represent genes, TFs and miRNAs. c This panel allows to define a set of parameters related to the integration of miRNAs into the generated STPs. They can be used to control the size of the generated networks. d This panel allows the user to export the results in the form of a text file including all identified circuits or to delete the current experiment and start with a new one. e The Cytoscape node and edge tables. They can be used to access detailed information on the nodes and arcs of the identified STPs

Full size image
Fig. 2
figure 2

CyTRANSFINDER software architecture. CyTRANSFINDER processes three inputs: (1) the source list of genes (SRL), (2) the destination list of genes (DGL), and (3) the STP pattern (STPP) to be reconstructed. Its fusion engine connects to several on-line repositories to collect regulatory information used to infer STPs connecting source and destination genes according to the selected STP pattern. The identified STPs are then purged to remove duplicated nodes and arcs and the result is returned to the user as a Cytoscape network or exported in the form of a tab separated text file

Full size image

CyTRANSFINDER processes three main inputs:

  1. 1.

    Source Regulator List (SRL): is a list of regulators working as root nodes of the inferred STPs. Regulators can be either genes or miRNAs (DNM3 – Dynamin 3 – NCBI ID: 26052 in the example of Fig. 1).

  2. 2.

    Destination Gene List (DGL): is a list of genes working as leaves of the inferred STPs (ALCAM – activated leukocyte cell adhesion molecule – NCBI ID: 214 in the example of Fig. 1).

  3. 3.

    STP Pattern (STPP): is a pattern provided through an internal STPP specification language. It describes a general template of regulators to be identified to connect a source node s n∈S R L to a destination gene d n∈D G L.

The user can enter the desired inputs through panel (A) of Fig. 1. Both (SRL) and DGL are provided in a text file formatted as described in Fig. 2. Each gene can be defined by either the gene symbol or the NCBI gene ID, while miRNAs are defined using the miRBase identifier (e.g., hsa-mir-214). The STPP can be chosen from a list of default patterns or customized as described in the following sections.

RUN button in panel (A) of Fig. 1 starts the STP search. CyTRANSFINDER fusion engine connects to several on-line repositories to collect regulatory information used to infer STPs connecting source and destination nodes according to the selected STPP. The identified STPs are then purged to remove duplicated nodes and arcs and displayed in the form of a network (panel (B) – Fig. 1).

STPP specification language

A STPP is a general sequence of regulators connecting two genes (or a miRNA to a gene) that are of interest for a given biological question. CyTRANSFINDER provides the user with a simple STPP Specification Language to specify the desired signaling pattern. Three regulators are available in CyTRANSFINDER: (1) gene, (2) TF, and (3) miRNA. The basic regulators are combined into a set of 6 possible interactions that are reported in Table 1 together with the notation used for their specification. They represent realistic biological interactions among the three regulators. The user is free to combine all the interactions of Table 1 in order to describe the desired signaling pattern, which is provided to CyTRANSFINDER in the form of a single text file loaded when setting STPP Type to “Custom STPP” in panel (A) of Fig. 1.

Table 1 STPP Specification Language
Full size table

Five default STPPs representing common recurring patterns often analyzed in the literature are directly embedded in the plugin (Fig. 3); three STPPs starting from genes and two starting from miRNAs. The Direct miRNA STPP is the simplest pattern. A source gene hosts a miRNA which also targets one of the destination genes. The Indirect miRNA STPP, is similar to the Direct miRNA STPP, but it involves a TF as miRNA mediator for the destination genes regulation. The Double miRNA indirect STPP is the most complex pattern that combines the two previous ones into two levels of indirect regulation: the first one is an Indirect miRNA STPP, which regulates a Direct miRNA STPP that targets the destination gene. Additionally, a version of the Indirect miRNA STPP and of the Double miRNA indirect STPP starting from a miRNA instead of a gene are available and named Indirect s. miRNA and Double s. miRNA indirect, respectively.

Fig. 3
figure 3

CyTRANSFINDER built-in STPPs. The figure presents the five default STPPs embedded in the plugin. i) Direct miRNA STPP is the simplest pattern: a source gene hosts a intragenic mirna miRNA or is located close to the region of an intergenic miRNA, which targets one of the destination genes. ii) Indirect miRNA STPP, is pretty similar to the Direct miRNA STPP, but it involves a TF as miRNA mediator for the regulation of the destination genes. iii) the miRNA sourced version of (ii). iv) Double miRNA indirect STPP is the most complex pattern. It involves two levels of regulation; the first indirect regulation is modeled on top of an Indirect miRNA STPP, which regulates a Direct miRNA STPP that targets the destination genes. v) The miRNA sourced version of (iv)

Full size image

Data fusion engine

The RUN button available in panel (A) of Fig. 1 starts the CyTRANSFINDER data fusion engine. The engine first parses the SRL, the DGL and the STPP provided by the user. It then connects to several external repositories to obtain interaction data to search for the existence of the STPP among the genes contained in SRL and DGL.

Figure 4 provides a high-level pseudo-code of the implemented data-fusion algorithm. The main algorithm is described in the STPFinder procedure (Fig. 4 - lines 1–22). This procedure receives as parameters the source and destination node lists (i.e., SRL and DGL) and the STPP. The STPP is an ordered list of regulators S T P P=(r 1,r 2,·,r n ), with r i ∈{TF,gene,miRNA}. The produced STPs are organized into a set of levels (stplevels in Fig. 4). Each level contains a set of nodes and corresponds to one of the elements of STPP. At the beginning of the search the first level is initialized with the nodes contained in SRL (Fig. 4 - lines 2). Nodes of adjacent levels are connected through a set of interactions (stpinters in Fig. 4).

Fig. 4
figure 4

CyTRANSFINDER data fusion algorithm. A pseudocode description of the main steps carried out by the plugin to integrate different data sources and to construct the final STP network

Full size image

The search procedure is an iterative process that analyzes couples of consecutive STPP elements, i.e., S T P P i and S T P P i+1 with i∈[1,|S T P P|−1] (Fig. 4 - lines 4–21). For each node available at level i (Fig. 4 - lines 7–18) the procedure searches a set of target nodes to add to level i+1 through the search function (Fig. 4 - lines 8). The way this search works (Fig. 4 - lines 23–38) depends on the type of regulators to search at level i and i+1 of STPP, which in turn requires to connect to different repositories to obtain interaction data. Additional details regarding this process will be provided later in this section. Each target node identified with this procedure is added to the set of nodes at level i+1 (Fig. 4 - line 14) and the set of interactions between couples of nodes at the two levels is recorded as well (Fig. 4 - line 15). When reaching the last couple of STPP elements, all identified interactions are finally filtered to limit them to those ending to one of the nodes available in DGL (Fig. 4 - lines 9–11).

At the end of this iterative process, stplevels and stpinters are used to build a Cytoscape network representing the inferred STPs. stplevels contains the set of nodes of the final network organized in levels and stpinters contains the set of arcs connecting the different nodes. Panel (B) of Fig. 1 shows an example of network identified when searching for the Indirect miRNA STPP between DNM3 and ALCAM. This network represents the main output provided by CyTRANSFINDER. Different symbols and colors have been used to make it easy to identify the different types of regulators in the network.

The same information can also be exported into a tab separated plaintext file that enumerates all identified STPs (Export button of panel (D) – Fig. 1). Each row of the file reports a single STP (i.e., a signaling chain from one node in SRL to a node in DGL according to the STPP) and each column represents a regulator in the signaling chain (i.e., a gene, miRNA or TF). This file is obtained by searching all possible paths that connect nodes at the source level with nodes at the destination level. This format is particularly helpful for fast data inspection, especially when the number of discovered STPs is high, which eventually results in a very complex network difficult to visualize.

The remaining of this section focuses on the way interactions are obtained and integrated from public repositories.

Transcription factors interaction data

Transcription Factors (TFs) related to gene entities are extracted from TargetMine [22] using its RESTful interface, which allows to search for TFs given a target gene. This information is required to search for interactions of type 1,2 and 6 of Table 1. It is important to highlight that TargetMine does not provide any information regarding the up- or down- regulatory activity of a TF; users must eventually resort to manual validation in order to understand the exact regulatory effect.

miRNA interaction data

Two different repositories are exploited to retrieve miRNA based interactions.

Interactions of type 3 and 4 from Table 1 are obtained through the miRIAD repository [23]. miRIAD is a web search tool designed to access integrated information concerning intragenic microRNAs and their host genes. The miRIAD database references annotated genes from human genome (hg19) and miRNAs annotated from miRBase (version 19). Given a gene, CyTRANSFINDER uses miRIAD to search for miRNAs “hosted” by the gene. Two types miRNA interactions can be identified with this procedure. The main class is represented by the intragenic miRNAs, which are mapped to intragenic loci of protein coding genes (namely “host genes”). Previous studies have suggested that these miRNAs are transcribed in parallel with their host transcripts [24, 25] therefore creating a direct signaling link. In addition to this, miRIAD enables to search for intergenic miRNAs located in an intergenic region close the the analyzed gene. The relationship between a gene and the intergenic miRNAs is weaker than the one of intragenic miRNAs. Nevertheless, it may represent a valuable information when performing explorative analysis. By acting on the miRNA host type control (panel (C) – Fig. 1) the user is free to work with intragenic only interactions or both intragenic and intergenic interactions.

Interactions of type 5 and 6 from Table 1 are instead obtained from TargetHUB [26]. This web-service provides a programmer friendly interface to access multiple repositories of miRNA target genes with a uniform set of APIs. TargetHUB RESTful interface allows users to interrogate information from four different databases: miRTarBase [27], TargetScan [28], PicTar [29], and miRanda [30]. Using TargetHUB, the list of target genes of a miRNA can be easily retrieved.

CyTRANSFINDER exploits TargetHUB functionalities to allow users to filter miRNA targets (miRNA target settings of panel (C) – Fig. 1). Filtering miRNA targets is a very crucial step during STP discovery. In fact, the lack of miRNA specificity and the large amount of possible miRNA targets (in the order of thousands) may easily increase the complexity of the networks generated by the plugin. A wise usage of filters can dramatically reduce the analysis time. Available filters include the possibility of selecting specific miRNA target databases or to perform majority voting across multiple databases to have a mandatory minimum or exact set of confirmations for the target selection. Moreover, the user can decide to work with regulatory information regarding mature o precursor miRNAs.

Finally, miRNA target settings of panel (C) include an additional control that acts in a opposite way with respect to the other controls, and tries to enlarge the obtained network. This control acts after the full STP search is concluded adding to the network all external miRNAs targeting at least one of the nodes identified in the generated network (i.e., not hosted by one of the network nodes). This option is particularly useful whenever users are focusing on the role of miRNAs in the studied phenomena.

Results and discussion

In this section we show the capability of CyTRANSFINDER by presenting its application in the framework of a research activity on human melanoma performed by the authors of this paper. In previous studies we and others identified that miR-146a and miR-214 are involved in melanoma growth and metastasis formation by modulating several target genes. We are therefore interested in performing discovery analysis searching for STPs involving these two miRNAs. This represents a typical biological question for which CyTRANSFINDER can provide explorative analysis support.

STPs involving human miR-146a analysis

Human miR-146a is located on the positive strand of chromosome 5. Although it is an intergenic miRNA and it does not lie inside a host protein-coding gene, it is overlapped to a manually-annotated long-intergenic-noncoding RNA (lincRNA), CTC-231O11.1 ([31] and http://www.ensembl.org). miR-146a has a crucial role in the immune and inflammatory response, as well as in many human pathologies including muscle disorders, cancer and metastasis [31]. We and others found that miR-146a has a dual role during melanoma development and progression, favoring primary tumor growth while inhibiting metastatic dissemination [32]. We are interested in exploring STPs involving miR-146a to identify new regulatory paths of interest for the melanoma progression. In order to exploit CyTRANSFINDER for this purpose we need to create a SRL and a DGL file.

Our SRL list contains the miR-146a (miRBase identifier hsa-mir-146a) [see Additional file 1 – srl.txt]. We considered instead the set of miR-146a conserved target genes according to TargetScan 5.2 algorithm (224 genes) as DGL [see Additional file 1 – dgl.txt].

We performed an analysis using the Indirect s. miRNA STPP using the default setting of the plugin for all filters, and we have been able to identify a set of 312 STPs [see Additional file 1 – Indirect-miRNA-STPP] starting from miR-146a. Among all, we got particularly interested in TFAP2C (AP-2 γ). In our previous studies, we identified the central role of TFAP2C in melanoma progression, and we are particularly interested in studying the STPs involving this transcription factor ([33, 34]).

All the TFAP2C-mediated STPs are listed in Table 2 and the related network is reported in Fig. 5. Notably, miR-146a relation with TFAP2C was completely unknown.

Fig. 5
figure 5

miR-146a Indirect s. miRNA STPs involving TFAP2C. Subnetwork of the Indirect s. miRNA STPs using human miR-146a as source intergenic miRNA, miR-146a targets according to TargetScan 5.1 as destination genes and involving TFAP2C as a hub transcription factor. Results are computed using miRNA targets confirmed in at least one source database

Full size image
Table 2 List of the Indirect s. miRNA STPs obtained with CyTRANSFINDER using human miR-146a as the source intergenic miRNA, miR-146a targets according to TargetScan 5.1 as destination genes and involving TFAP2C as a hub transcription factor
Full size table

Given the interest of this result, we performed gene expression analysis via quantitative Real Time Polymerase Chain Reaction (qRT-PCR) to further investigate this relation. To obtain transient miR-146a or non-specific control (pre-Cntrl) expression, human melanoma MA-2 cells (cultured as in [33]) were transfected using HiPerFect (Qiagen) reagent, according to the manufacturer’s instructions. Total RNA was isolated 48h later from using TRIzol®;Reagent (Invitrogen Life Technologies). 1 μg of DNAse-treated RNA (RQ1 RNase-Free DNase, Promega) was retrotranscribed with High-Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific,) and qRT-PCRs were carried out using gene-specific primers for TFAP2C mRNA detection (fw: TCCACGACATGCCTCACCA, rv: TCCTTCTGACAGGGGAGGTTCA).

Quantitative normalization was performed on the expression of the GAPDH gene (qRT-PCR QuantiTect Primer assay QT01192646, Qiagen). The relative expression levels between samples were calculated using the comparative delta Ct (threshold cycle number) method (2−ΔΔCt) with a control sample as the reference point [35].

Data are presented as mean ± s.e.m. (standard error of the mean) and Two tailed Student’s t-test was used for comparison, with **P <0.01 considered to be statistically significant. Based on these experiments we were able to experimentally verify that miR-146a is able to downmodulate TFAP2C expression upon transient overexpression in human melanoma cells (see Fig. 6).

Fig. 6
figure 6

miR-146a overexpression leads to reduced TFAP2C mRNA levels. Quantitative-Real Time PCR (qRT-PCR) evaluation of TFAP2C mRNA was performed in melanoma cells upon miR-146a overexpression, compared to controls (pre-146a vs pre-Cntrl). Three independent preparations of melanoma cells RNA were used and results were pooled together. **P <0.01

Full size image

Furthermore, by searching the literature for TF ↦ gene STPs involving TFAP2C and our DGL, we found a paper by Woodfield and colleagues where direct regulation by TFAP2C on GGA2, NOTCH2, NPR3 and TDRKH promoter regions was demonstrated by chromatin immunoprecipitation followed by sequencing (ChIP-Seq) analysis [36], as shown in Table 2.

Next, we also searched for Double Indirect s. miRNA STPP involving miR-146a and the selected DGL. Given the complexity of this pattern, we performed the analysis restricting to miRNA targets confirmed in at least two databases out of the four available in TargetHUB. In this case, we obtained a significantly shorter list of records, that is reported in Table 3 and is visually reproduced in Fig. 7.

Fig. 7
figure 7

miR-146a Double Indirect s. miRNA STPs. Network of the Double Indirect s. miRNA STPs using human miR-146a as source intergenic miRNA, miR-146a targets according to TargetScan 5.1 as destination genes and miRNA targets confirmed in at least two of the source databases

Full size image
Table 3 List of Double Indirect s. miRNA STPs obtained with CyTRANSFINDER using human miR-146a as the source intergenic miRNA, targets according to TargetScan 5.1 as destination genes and miRNA targets confirmed in at least two of the source databases
Full size table

Notably, the STPs identified by CyTRANSFINDER seem of biological relevance, since our paths link miR-146a to TFs SMAD4 and BRCA1 (miRNA ↦ TF regulation), as well-established in literature (Table 3). SMAD4 is a key transcription factor involved in the TGF- β mediated response [37], while BRCA1 is involved in the DNA damage repair and is one of the main mutated genes in familial breast and ovarian cancers [38]. SMAD4- and BRCA1- regulated genes (TF target (gene)) that we obtained, POR and CYP19A1 respectively, are again well-established in literature (Table 3), and their involvement downstream of miR-146a could be very interesting, since they both have a role in hormones production and cancer. POR gene codifies for cytochrome P450 oxidoreductase enzyme, which catalyzes the biosynthesis of steroid hormones and metabolize drugs [39], while CYP19A1 codes for the aromatase enzyme that converts androgens into estrogens, and which dysregulation may affect estrogen production in breast cancer cases [40]. POR and CYP19A1 genes host one microRNA each, miR-4651 and miR-4713, respectively (gene ↦ miRNA regulation; they are indicated as Intragenic microRNA). miR-4651 and miR-4713 target genes predicted by at least 2 algorithms are listed in Table 3 (DG), but, unfortunately, since they are recently discovered miRNAs, none of their target genes has been experimentally validated yet.

STPs involving human miR-214 analysis

Human miR-214 gene is located in the chromosomal region 1q24.3, in intron 14 of the Dynamin-3 gene (DNM3) inside an almost 8 kb-long noncoding RNA, named DNM3os. This transcript contains the sequences for miR-214 and miR-199a-2, two clustered miRs that are approximately 6 kb apart. miR-214 is deregulated in a variety of human tumors including melanoma, breast, ovarian, gastric, and hepatocellular carcinomas as reviewed in [41]. In melanoma, we demonstrated that miR-214 has essential roles in regulating invasiveness, extravasation and metastasis formation [33, 34]. In particular, we identified a signature of 73 genes whose expression was driven by miR-214 [33].

In order to identify new molecular pathways underlying miR-214-mediated regulation of these genes we took advantage of CyTRANSFINDER. Differently from the previous case, to show the use of the software on a STPP starting from a gene, we used DNM3, the host gene of miR-214, as SRL [see Additional file 2 – srl.txt] and the miR-214-dependent signature mentioned above as the DGL [see Additional file 2 – dgl.txt]. We searched for Double miRNA Indirect STPs identifying 312 STPs involving different transcription factors (TFs), miRNA host genes (TF target gene) targeted by these TFs and cognate intragenic miRNAs (ImiRNAs) as nodes [see Additional file 2 – Double-indirect-miRNA-STPP.xlsx].

Interestingly enough, the majority of the STPPs were controlled by two of the most relevant TFs for melanoma biology, the transcription factor AP-2 gamma TFAP2C (AP-2 γ) [33] and the cAMP responsive element binding protein 1, CREB1 [42]. We focused our attention on the STPs driven by these two TFs and we selected a subgroup (101) of STPPs, containing well-described intragenic miRNAs (Intragenic microRNAs) as nodes (Table 4 and Fig. 8).

Fig. 8
figure 8

miR-214 Double Indirect miRNA STPs. Sub-Network of a selection of 101 interesting Double Indirect miRNA STPs using human DNM3 as source gene, a signature of 73 genes published in [33] as destination genes, and involving TFAP2C and CREB1 transcription factors. Results are computed using miRNA targets confirmed in at least two databases

Full size image
Table 4 List of a selection (101) of the Double Indirect miRNA STPPs obtained with CyTRANSFINDER using the host gene of miR-214, DNM3, as the source gene, a signature of 73 genes published in [33] as destination genes, and involving TFAP2C and CREB1 transcription factors
Full size table

Searching the literature for potential validations of these STPPs, we were able to find partial validations. The connection between miR-214 and TFAP2C was clearly demonstrated in our previous work [33], where we showed the direct targeting of miR-214 on TFAP2C 3 ′-UTR; while no data linking miR-214 and CREB1 were found. No connections were observed for either TFAP2C or CREB1 and the host genes of intragenic miRNAs (TF target) present in the STPs, except for the phospholipid-dependent protein-serine/threonine kinase PRKCA gene. PRKCA plays a major role in intracellular signaling pathways associated with transformation and tumor progression and its expression was shown to be under the control of TFAP2 transcription factor family [43].

Looking for potential targeting of the analyzed intragenic miRNAs (Intragenic microRNA) on genes of the DGL we found numerous experimental validations in the literature. In particular, we were interested in STPPs driven by CREB1 since potential miR-214 ↦CREB1 connections could open up new lines of research in understanding miR-214-driven metastatization. Among CREB1-controlled STPs we found SREBP2 (TF target) and miR-33a (Intragenic microRNA) that are known to be co-regulated [44] and we previously demonstrated to be downregulated by miR-214 [45]. Very recently, Zhou and colleagues demonstrated miR-33a tumor suppressive role in melanoma, thus suggesting a potential additional effect of miR-214 in promoting melanoma malignancy via the downregulation of another miRNA, miR-33a [46]. We demonstrated the ability of miR-214 to promote melanoma progression by downregulation of miR-148b at least partially via TFAP2C regulation, thus leading to miR-148b targets derepression, such as ALCAM [34]. Interestingly enough, we were able to find miR-148b and ALCAM in one of the STPs, but surprisingly from CyTRANSFINDER analysis CREB1, and not TFAP2C, resulted to be the master regulator of the pathway. These new data are very interesting for us and we would like to investigate this potential pathway more in detail. In fact, it has been demonstrated that CREB1 is able to regulate TFAP2A expression in melanoma [47], so we could hypothesize a double control of miR-214 on TFAP2C, direct, via targeting, and, indirect, via CREB1, thus leading to a strong promotion of melanoma progression. Finally, another STP interestingly linked miR-214 to miR-15b-5p. In particular, 5 (SEMA3A, APP, PAK2, TFAP2A and ARHGAP12) out of 10 DGL genes resulted to be validated targets of this miRNA and moreover, miR-15b was shown to be involved in tumor cell proliferation and apoptosis in malignant melanoma [48].

Conclusions

Here we presented a new plugin for Cytoscape, CyTRANSFINDER that provides support to discover three-component signal transduction pathways with the presence of miRNAs, TFs and genes. Differently from other tools, the plugin is specifically designed to perform exploratory analysis and to identify new biological circuits to be tested in laboratory. Therefore, it only relies on aggregation of complex repositories without requiring any expression data.

To show the capabilities of this plugin we applied it to a real use case involving the study of two miRNAs that are particularly relevant in human melanoma progression. Taken together, our analyses on the STPs generated by CyTRANSFINDER unravelled many relevant potential pathways regulated by miR-146a and miR-214 in human physiology and pathology; some of these are validated in literature, while others were validated by us.

This should give a clear view of the potential this tool has to support biologists in discovering novel signal transduction pathways regulated by miRNAs and transcription factors.

A detailed tutorial containing a step-by-step guide covering all CyTRANSFINDER features is available at http://apps.cytoscape.org/apps/cytransfinder.

We are currently at the first release of this tool. One of the critical aspects of this implementation is the huge amount of data the plugin needs to retrieve from the Internet and to process internally. This in turns requires processing time that ranges from a few minutes for very simple lists of genes to several hours for complex lists of genes such as the one used in the proposed case studies. Future releases will address this specific aspect by introducing several caching mechanism enabling to reduce the network traffic and to significantly increase the computation speed.

Moreover, we are also interested in inserting new features to make the generated networks as specific as possible. Recently, we discovered mimiRNA, a database of miRNA expression profiles data across different tissues and cell lines [49]. mimiRNA incorporates a sample classification algorithm that groups identical miRNA or mRNA experiments from separate sources and provides reliable expression profiles of miRNA in different tissues and cell lines. We plan to integrate mimiRNA into CyTRANSFINDER as an additional filter enabling to select only intra/intergenic miRNAs that are expressed in a set of tissues or cell lines of interest for the user, thus reducing the size of the generated networks.

Finally, we are evaluating the possibility of defining a scoring approach to help users to select STPs based on different search criteria (e.g., involvement of the identified regulator in a disease based on literature, strength of the single identified interactions based on scores available in other external databases, etc.).

Ethics

Our study does not involve humans, human data or animals. For this reason no ethical approvals are required.

Availability and requirements

References

  1. Yu D, Kim M, Xiao G, Hwang TH. Review of biological network data and its applications. Genomics Inform. 2013; 11(4):200–10. doi:10.5808/GI.2013.11.4.200.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Alon U. Network motifs: theory and experimental approaches. Nat Rev Genet. 2007; 8(6):450–61. doi:10.1038/nrg2102.

    Article  CAS  PubMed  Google Scholar 

  3. Hartwell LH, Hopfield JJ, Leibler S, Murray AW. From molecular to modular cell biology. Nature. 1999; 402(6761 Suppl):47–52. doi:10.1038/35011540.

    Article  Google Scholar 

  4. Yeang CH, Ideker T, Jaakkola T. Physical network models. J Comput Biol. 2004; 11(2-3):243–62. doi:10.1089/1066527041410382.

    Article  CAS  PubMed  Google Scholar 

  5. Forsburg SL. The art and design of genetic screens: yeast. Nat Rev Genet. 2001; 2(9):659–8. doi:10.1038/35088500.

    Article  CAS  PubMed  Google Scholar 

  6. Neves SR, Iyengar R. Modeling of signaling networks. Bioessays. 2002; 24(12):1110–7. doi:10.1002/bies.1154.

    Article  CAS  PubMed  Google Scholar 

  7. Choi C, Crass T, Kel A, Kel-Margoulis O, Krull M, Pistor S, Potapov A, Voss N, Wingender E. Consistent re-modeling of signaling pathways and its implementation in the transpath database. Genome Inform. 2004; 15(2):244–54.

    CAS  PubMed  Google Scholar 

  8. Sachs K, Perez O, Pe’er D, Lauffenburger DA, Nolan GP. Causal protein-signaling networks derived from multiparameter single-cell data. Science. 2005; 308(5721):523–9. doi:10.1126/science.1105809.

    Article  CAS  PubMed  Google Scholar 

  9. Steffen M, Petti A, Aach J, D’haeseleer P, Church G. Automated modelling of signal transduction networks. BMC bioinformatics. 2002; 3(1):34. doi:12413400.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Liu Y, Zhao H. A computational approach for ordering signal transduction pathway components from genomics and proteomics data. BMC Bioinformatics. 2004; 5:158. doi:10.1186/1471-2105-5-158.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Koyutürk M, Grama A, Szpankowski W. An efficient algorithm for detecting frequent subgraphs in biological networks. Bioinformatics. 2004; 20 Suppl 1:200–7. doi:10.1093/bioinformatics/bth919.

    Article  Google Scholar 

  12. Shlomi T, Segal D, Ruppin E, Sharan R. Qpath: a method for querying pathways in a protein-protein interaction network. BMC Bioinformatics. 2006; 7:199. doi:10.1186/1471-2105-7-199.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Bebek G, Yang J. Pathfinder: mining signal transduction pathway segments from protein-protein interaction networks. BMC Bioinformatics. 2007; 8:335. doi:10.1186/1471-2105-8-335.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Filipowicz W, Bhattacharyya SN, Sonenberg N. Mechanisms of post-transcriptional regulation by micrornas: are the answers in sight?Nat Rev Genet. 2008; 9(2):102–14. doi:10.1038/nrg2290.

    Article  CAS  PubMed  Google Scholar 

  15. Le TD, Liu L, Liu B, Tsykin A, Goodall GJ, Satou K, Li J. Inferring microrna and transcription factor regulatory networks in heterogeneous data. BMC Bioinformatics. 2013; 14:92. doi:10.1186/1471-2105-14-92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Le Béchec A, Portales-Casamar E, Vetter G, Moes M, Zindy PJ, Saumet A, Arenillas D, Theillet C, Wasserman WW, Lecellier CH, Friederich E. Mir@nt@n: a framework integrating transcription factors, micrornas and their targets to identify sub-network motifs in a meta-regulation network model. BMC Bioinformatics. 2011; 12:67. doi:10.1186/1471-2105-12-67.

    Article  PubMed  PubMed Central  Google Scholar 

  17. 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; 13(11):2498–504. doi:10.1101/gr.1239303.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Milo R, Shen-Orr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U. Network motifs: simple building blocks of complex networks. Science. 2002; 298(5594):824–7. doi:10.1126/science.298.5594.824.

    Article  CAS  PubMed  Google Scholar 

  19. Lenski RE, Ofria C, Pennock RT, Adami C. The evolutionary origin of complex features. Nature. 2003; 423(6936):139–44. doi:10.1038/nature01568.

    Article  CAS  PubMed  Google Scholar 

  20. Lipson H, Pollack JB, Suh NP. On the origin of modular variation. Evolution. 2002; 56(8):1549–56.

    Article  PubMed  Google Scholar 

  21. Garcia-Garcia J, Guney E, Aragues R, Planas-Iglesias J, Oliva B. Biana: a software framework for compiling biological interactions and analyzing networks. BMC Bioinformatics. 2010; 11:56. doi:10.1186/1471-2105-11-56.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Chen YA, Tripathi LP, Mizuguchi K. Targetmine, an integrated data warehouse for candidate gene prioritisation and target discovery. PLoS ONE. 2011; 6(3):17844. doi:10.1371/journal.pone.0017844.

    Article  Google Scholar 

  23. Hinske LC, França GS, Torres HAM, Ohara DT, Lopes-Ramos CM, Heyn J, Reis LFL, Ohno-Machado L, Kreth S, Galante PAF. miRIAD—integrating microRNA inter- and intragenic data. Database. 2014:bau099. doi:10.1093/database/bau099. Accessed 6 October 2014.

  24. Rodriguez A, Griffiths-Jones S, Ashurst JL, Bradley A. Identification of mammalian microrna host genes and transcription units. Genome Res. 2004; 14(10A):1902–10. doi:10.1101/gr.2722704.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Hinske LCG, Galante PAF, Kuo WP, Ohno-Machado L. A potential role for intragenic mirnas on their hosts’ interactome. BMC Genomics. 2010; 11:533. doi:10.1186/1471-2164-11-533.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Manyam G, Ivan C, Calin GA, Coombes KR. targetHub: a programmable interface for miRNA-gene interactions. Bioinformatics. 2013. doi:10.1093/bioinformatics/btt439 http://bioinformatics.oxfordjournals.org/content/early/2013/09/06/bioinformatics.btt439.full.pdf+html.

  27. Hsu SD, Lin FM, Wu WY, Liang C, Huang WC, Chan WL, Tsai WT, Chen GZ, Lee CJ, Chiu CM, Chien CH, Wu MC, Huang CY, Tsou AP, Huang HD. miRTarBase: a database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. 2011; 39(Database issue):163–9. doi:10.1093/nar/gkq1107.

    Article  Google Scholar 

  28. 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. doi:10.1016/j.cell.2004.12.035035.

    Article  CAS  PubMed  Google Scholar 

  29. Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, Rajewsky N. Combinatorial microRNA target predictions. Nat Genet. 2005; 37(5):495–500.

    Article  CAS  PubMed  Google Scholar 

  30. John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS. Human MicroRNA targets. PLoS Biol. 2004; 2(11):363. doi:10.1371/journal.pbio.0020363.

    Article  Google Scholar 

  31. Li L, Chen XP, Li YJ. Microrna-146a and human disease. Scand J Immunol. 2010; 71(4):227–31. doi:10.1111/j.1365-3083.2010.02383.x.

    Article  CAS  PubMed  Google Scholar 

  32. Forloni M, Dogra SK, Dong Y, Conte Jr D, Ou J, Zhu LJ, Deng A, Mahalingam M, Green MR, Wajapeyee N. mir-146a promotes the initiation and progression of melanoma by activating notch signaling. Elife. 2014; 3:01460. doi:10.7554/eLife.01460.

    Article  Google Scholar 

  33. Penna E, Orso F, Cimino D, Tenaglia E, Lembo A, Quaglino E, Poliseno L, Haimovic A, Osella-Abate S, De Pittà C, Pinatel E, Stadler MB, Provero P, Bernengo MG, Osman I, Taverna D. microrna-214 contributes to melanoma tumour progression through suppression of tfap2c. EMBO J. 2011; 30(10):1990–2007. doi:10.1038/emboj.2011.102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Penna E, Orso F, Cimino D, Vercellino I, Grassi E, Quaglino E, Turco E, Taverna D. mir-214 coordinates melanoma progression by upregulating alcam through tfap2 and mir-148b downmodulation. Cancer Res. 2013; 73(13):4098–111. doi:10.1158/0008-5472.CAN-12-36863686.

    Article  CAS  PubMed  Google Scholar 

  35. Bookout AL, Mangelsdorf DJ. Quantitative real-time pcr protocol for analysis of nuclear receptor signaling pathways. Nucl Recept Signal. 2003; 1:012. doi:10.1621/nrs.01012.

    Article  Google Scholar 

  36. Woodfield GW, Chen Y, Bair TB, Domann FE, Weigel RJ. Identification of primary gene targets of tfap2c in hormone responsive breast carcinoma cells. Genes Chromosomes Cancer. 2010; 49(10):948–62. doi:10.1002/gcc.20807.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. He Y, Huang C, Sun X, Long X-R, Lv X-W, Li J. Microrna-146a modulates tgf-beta1-induced hepatic stellate cell proliferation by targeting smad4. Cell Signal. 2012; 24(10):1923–30. doi:10.1016/j.cellsig.2012.06.003.

    Article  CAS  PubMed  Google Scholar 

  38. Gu Y, Zhang M, Peng F, Fang L, Zhang Y, Liang H, Zhou W, Ao L, Guo Z. The brca1/2-directed mirna signature predicts a good prognosis in ovarian cancer patients with wild-type brca1/2. Oncotarget. 2015; 6(4):2397–406.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Tee MK, Huang N, Damm I, Miller WL. Transcriptional regulation of the human p450 oxidoreductase gene: hormonal regulation and influence of promoter polymorphisms. Mol Endocrinol. 2011; 25(5):715–31. doi:10.1210/me.2010-0236.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Ghosh S, Lu Y, Katz A, Hu Y, Li R. Tumor suppressor brca1 inhibits a breast cancer-associated promoter of the aromatase gene (cyp19) in human adipose stromal cells. Am J Physiol Endocrinol Metab. 2007; 292(1):246–52. doi:10.1152/ajpendo.00242.2006.

    Article  Google Scholar 

  41. Penna E, Orso F, Taverna D. mir-214 as a key hub that controls cancer networks: small player, multiple functions. J Invest Dermatol. 2015; 135(4):960–9. doi:10.1038/jid.2014.479.

    Article  CAS  PubMed  Google Scholar 

  42. Poser I, Bosserhoff AK. Transcription factors involved in development and progression of malignant melanoma. Histol Histopathol. 2004; 19(1):173–88.

    CAS  PubMed  Google Scholar 

  43. Clark JH, Haridasse V, Glazer RI. Modulation of the human protein kinase c alpha gene promoter by activator protein-2. Biochemistry. 2002; 41(39):11847–56.

    Article  CAS  PubMed  Google Scholar 

  44. Bommer GT, MacDougald OA. Regulation of lipid homeostasis by the bifunctional srebf2-mir33a locus. Cell Metab. 2011; 13(3):241–7. doi:10.1016/j.cmet.2011.02.004.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Politano G, Benso A, Di Carlo S, Orso F, Savino A, Taverna D. A computational study to identify tp53 and srebf2 as regulation mediators of mir-214 in melanoma progression. In: Proceedings of the International Conference on Bioinformatics Models, Methods and Algorithms (BIOINFORMATICS). Portugal: INSTICC: 2014. p. 49–56. doi:9789897580123.

    Google Scholar 

  46. Zhou J, Xu D, Xie H, Tang J, Liu R, Li J, Wang S, Chen X, Su J, Zhou X, Xia K, He Q, Chen J, Xiong W, Cao P, Cao K. mir-33a functions as a tumor suppressor in melanoma by targeting hif-1 α. Cancer Biol Ther. 2015; 16(6):846–55. doi:10.1080/15384047.2015.1030545.

    Article  CAS  PubMed  Google Scholar 

  47. Melnikova VO, Dobroff AS, Zigler M, Villares GJ, Braeuer RR, Wang H, Huang L, Bar-Eli M. Creb inhibits ap-2alpha expression to regulate the malignant phenotype of melanoma. PLoS One. 2010; 5(8):12452. doi:10.1371/journal.pone.0012452.

    Article  Google Scholar 

  48. Satzger I, Mattern A, Kuettler U, Weinspach D, Voelker B, Kapp A, Gutzmer R. Microrna-15b represents an independent prognostic parameter and is correlated with tumor cell proliferation and apoptosis in malignant melanoma. Int J Cancer. 2010; 126(11):2553–62. doi:10.1002/ijc.24960.

    CAS  PubMed  Google Scholar 

  49. Ritchie W, Flamant S, Rasko JEJ. mimirna: a microrna expression profiler and classification resource designed to identify functional correlations between micrornas and their targets. Bioinformatics. 2010; 26(2):223–7. doi:10.1093/bioinformatics/btp649.

    Article  CAS  PubMed  Google Scholar 

  50. Li J, Huang J, Dai L, Yu D, Chen Q, Zhang X, Dai K. mir-146a, an il-1 β responsive mirna, induces vascular endothelial growth factor and chondrocyte apoptosis by targeting smad4. Arthritis Res Ther. 2012; 14(2):75. doi:10.1186/ar37983.

    Article  Google Scholar 

  51. Liu Z, Lu CL, Cui LP, Hu YL, Yu Q, Jiang Y, Ma T, Jiao DK, Wang D, Jia CY. Microrna-146a modulates tgf- β1-induced phenotypic differentiation in human dermal fibroblasts by targeting smad4. Arch Dermatol Res. 2012; 304(3):195–202. doi:10.1007/s00403-011-1178-0.

    Article  CAS  PubMed  Google Scholar 

  52. Xiao B, Zhu ED, Li N, Lu DS, Li W, Li BS, Zhao YL, Mao XH, Guo G, Yu PW, Zou QM. Increased mir-146a in gastric cancer directly targets smad4 and is involved in modulating cell proliferation and apoptosis. Oncol Rep. 2012; 27(2):559–66. doi:10.3892/or.2011.1514.

    CAS  PubMed  Google Scholar 

  53. Shen J, Ambrosone CB, DiCioccio RA, Odunsi K, Lele SB, Zhao H. A functional polymorphism in the mir-146a gene and age of familial breast/ovarian cancer diagnosis. Carcinogenesis. 2008; 29(10):1963–6. doi:10.1093/carcin/bgn172.

    Article  CAS  PubMed  Google Scholar 

  54. Hu Y, Ghosh S, Amleh A, Yue W, Lu Y, Katz A, Li R. Modulation of aromatase expression by brca1: a possible link to tissue-specific tumor suppression. Oncogene. 2005; 24(56):8343–8. doi::10.1038/sj.onc.1208985.

    Article  CAS  PubMed  Google Scholar 

  55. Lu Y, Kang T, Hu Y. Brca1/bard1 complex interacts with steroidogenic factor 1–a potential mechanism for regulation of aromatase expression by brca1. J Steroid Biochem Mol Biol. 2011; 123(1-2):71–8. doi:10.1016/j.jsbmb.2010.11.006.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Xue Y, Ouyang K, Huang J, Zhou Y, Ouyang H, Li H, Wang G, Wu Q, Wei C, Bi Y, Jiang L, Cai Z, Sun H, Zhang K, Zhang Y, Chen J, Fu XD. Direct conversion of fibroblasts to neurons by reprogramming ptb-regulated microrna circuits. Cell. 2013; 152(1-2):82–96. doi:10.1016/j.cell.2012.11.045.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Whisnant AW, Bogerd HP, Flores O, Ho P, Powers JG, Sharova N, Stevenson M, Chen CH, Cullen BR. In-depth analysis of the interaction of hiv-1 with cellular microrna biogenesis and effector mechanisms. MBio. 2013; 4(2):000193. doi:10.1128/mBio.00193-13.

    Article  Google Scholar 

  58. Kishore S, Jaskiewicz L, Burger L, Hausser J, Khorshid M, Zavolan M. A quantitative analysis of clip methods for identifying binding sites of rna-binding proteins. Nat Methods. 2011; 8(7):559–64. doi:10.1038/nmeth.1608.

    Article  CAS  PubMed  Google Scholar 

  59. Gottwein E, Corcoran DL, Mukherjee N, Skalsky RL, Hafner M, Nusbaum JD, Shamulailatpam P, Love CL, Dave SS, Tuschl T, Ohler U, Cullen BR. Viral microrna targetome of kshv-infected primary effusion lymphoma cell lines. Cell Host Microbe. 2011; 10(5):515–26. doi:10.1016/j.chom.2011.09.012.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Balakrishnan I, Yang X, Brown J, Ramakrishnan A, Torok-Storb B, Kabos P, Hesselberth JR, Pillai MM. Genome-wide analysis of mirna-mrna interactions in marrow stromal cells. Stem Cells. 2014; 32(3):662–73. doi:10.1002/stem.1531.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Zhu S, Pan W, Song X, Liu Y, Shao X, Tang Y, Liang D, He D, Wang H, Liu W, Shi Y, Harley JB, Shen N, Qian Y. The microrna mir-23b suppresses il-17-associated autoimmune inflammation by targeting tab2, tab3 and ikk- α. Nat Med. 2012; 18(7):1077–86. doi:10.1038/nm.2815.

    Article  CAS  PubMed  Google Scholar 

  62. Martin HC, Wani S, Steptoe AL, Krishnan K, Nones K, Nourbakhsh E, Vlassov A, Grimmond SM, Cloonan N. Imperfect centered mirna binding sites are common and can mediate repression of target mrnas. Genome Biol. 2014; 15(3):51. doi:10.1186/gb-2014-15-3-r51.

    Article  Google Scholar 

  63. Pillai MM, Gillen AE, Yamamoto TM, Kline E, Brown J, Flory K, Hesselberth JR, Kabos P. Hits-clip reveals key regulators of nuclear receptor signaling in breast cancer. Breast Cancer Res Treat. 2014; 146(1):85–97. doi:10.1007/s10549-014-3004-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Lipchina I, Elkabetz Y, Hafner M, Sheridan R, Mihailovic A, Tuschl T, Sander C, Studer L, Betel D. Genome-wide identification of microrna targets in human es cells reveals a role for mir-302 in modulating bmp response. Genes Dev. 2011; 25(20):2173–86. doi:10.1101/gad.17221311.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Hafner M, Landthaler M, Burger L, Khorshid M, Hausser J, Berninger P, Rothballer A, Ascano Jr M, Jungkamp AC, Munschauer M, Ulrich A, Wardle GS, Dewell S, Zavolan M, Tuschl T. Transcriptome-wide identification of rna-binding protein and microrna target sites by par-clip. Cell. 2010; 141(1):129–41. doi:10.1016/j.cell.2010.03.009.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Bischoff A, Bayerlová M, Strotbek M, Schmid S, Beissbarth T, Olayioye MA. A global microrna screen identifies regulators of the erbb receptor signaling network. Cell Commun Signal. 2015; 13:5. doi:10.1186/s12964-015-0084-z.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Karginov FV, Hannon GJ. Remodeling of ago2-mrna interactions upon cellular stress reflects mirna complementarity and correlates with altered translation rates. Genes Dev. 2013; 27(14):1624–32. doi:10.1101/gad.215939.113.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This work was supported by grants from Compagnia di San Paolo, Torino (IT), grants no. 2008.1054(DT), AIRC2010, 2013 (IG2010-10104DT; IG2013-14201DT), Fondazione Cassa di Risparmio Torino CRT, Torino (IT), grant no. 2014.1085DT and Italian Ministry of Education, University and Research (IT), grant no. PRIN 2010, MIND.

Author information

Authors and Affiliations

  1. Department of Control and Computer Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, Torino, 10129, Italy

    Gianfranco Politano, Alfredo Benso, Alessandro Savino & Stefano Di Carlo

  2. Molecular Biotechnology Center (MBC), Via Nizza, 52, Torino, 10126, Italy

    Francesca Orso, Monica Raimo & Daniela Taverna

  3. Dept. Molecular Biotechnology and Health Sciences, University of Torino, Via Nizza, 52, Torino, 10126, Italy

    Francesca Orso, Monica Raimo & Daniela Taverna

  4. Center for Complex Systems in Molecular Biology and Medicine, Via Accademia Albertina, 13, Torino, 10123, Italy

    Francesca Orso & Daniela Taverna

Authors
  1. Gianfranco Politano
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Francesca Orso
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Monica Raimo
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Alfredo Benso
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Alessandro Savino
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Daniela Taverna
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Stefano Di Carlo
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Stefano Di Carlo.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

GP and SDC conceived the idea and supervised the study. GP designed and implemented the plugin supported by AS. FO, DT and MR defined the case study and performed the literature validation and wet-lab experiments. GP, AB and SDC analyzed the results. GP and SDC drafted the manuscript. All authors read and approved the final manuscript.

Additional files

Additional file 1

∙ srl.txt: the file containing the SRL. ∙ dgl.txt: the file containing the DGL. ∙ Indirect-miRNA-STPP.xlsx: an excel file containing the full list of 312 Indirect miRNA STPPs identified by CyTRANSFINDER using human miR-146a as the source intergenic miRNA and miR-146a 223 target genes according to TargetScan 5.2 as DGL. Results are computed using miRNA targets confirmed in at lease one source database. SmiRNA: source intragenic miRNA; TF: Transcription Factor; DG: Destination Gene. (ZIP 201 kb)

Additional file 2

∙ srl.txt: the file containing the SRL. ∙ dgl.txt: the file containing the DGL. ∙ Indirect-miRNA-STPP.xlsx: an excel file containing the full list of 312 Indirect miRNA STPPs identified by CyTRANSFINDER using human miR-146a as the source intergenic miRNA and miR-146a 223 target genes according to TargetScan 5.2 as DGL. Results are computed using miRNA targets confirmed in at lease one source database. SmiRNA: source intragenic miRNA; TF: Transcription Factor; DG: Destination Gene.

an excel file containing the full list of 292 Double Indirect miRNA STPPs identified by TransFINDER using DNM3 as the Source Gene, the cognate human intragenic miR-214 as source intragenic miRNA (SmiRNA) and a previously described signature of 73 genes whose expression was driven by miR-214 [33] as destination genes (DG). TF: Transcription Factor; TF target (Gene): Transcription Factor target gene, which is also the host gene for a miRNA; intragenic miRNA: miRNA located inside the TF target gene; DG: Destination Gene, list of targets of the intragenic miRNAs predicted by at least two algorithms TS: Transcription Factor. (ZIP 192 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Politano, G., Orso, F., Raimo, M. et al. CyTRANSFINDER: a Cytoscape 3.3 plugin for three-component (TF, gene, miRNA) signal transduction pathway construction. BMC Bioinformatics 17, 157 (2016). https://doi.org/10.1186/s12859-016-0964-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12859-016-0964-2

Keywords

BMC Bioinformatics

ISSN: 1471-2105

Contact us