Skip to main content

Inferring miRNA-disease associations using collaborative filtering and resource allocation on a tripartite graph

Abstract

Background

Developing efficient and successful computational methods to infer potential miRNA-disease associations is urgently needed and is attracting many computer scientists in recent years. The reason is that miRNAs are involved in many important biological processes and it is tremendously expensive and time-consuming to do biological experiments to verify miRNA-disease associations.

Methods

In this paper, we proposed a new method to infer miRNA-disease associations using collaborative filtering and resource allocation algorithms on a miRNA-disease-lncRNA tripartite graph. It combined the collaborative filtering algorithm in CFNBC model to solve the problem of imbalanced data and the method for association prediction established multiple types of known associations among multiple objects presented in TPGLDA model.

Results

The experimental results showed that our proposed method achieved a reliable performance with Area Under Roc Curve (AUC) and Area Under Precision-Recall Curve (AUPR) values of 0.9788 and 0.9373, respectively, under fivefold-cross-validation experiments. It outperformed than some other previous methods such as DCSMDA and TPGLDA. Furthermore, it demonstrated the ability to derive new associations between miRNAs and diseases among 8, 19 and 14 new associations out of top 40 predicted associations in case studies of Prostatic Neoplasms, Heart Failure, and Glioma diseases, respectively. All of these new predicted associations have been confirmed by recent literatures. Besides, it could discover new associations for new diseases (or miRNAs) without any known associations as demonstrated in the case study of Open-angle glaucoma disease.

Conclusion

With the reliable performance to infer new associations between miRNAs and diseases as well as to discover new associations for new diseases (or miRNAs) without any known associations, our proposed method can be considered as a powerful tool to infer miRNA-disease associations.

Background

MicroRNA (miRNA) is a small RNA, about 22–26 nucleotides, which belongs to the noncoding RNA class [1]. Recent researches have shown that miRNAs are involved in many crucial biological processes like cell differentiation, proliferation, signal transduction, viral infection, and so on [2]. Identifying miRNA-disease associations could not only help us understand disease mechanism at miRNA level but also facilitate us in detecting disease biomarkers and discovering drugs for disease diagnosis, treatment, prognosis, and prevention. It has been confirmed that the dysregulations of the miRNAs are associated with the development and progression of various complex human diseases [3,4,5,6]. Until now, there are only a few known miRNA-disease associations in comparison with the number of newly discovered miRNAs. It is also tremendously expensive and time-consuming to do biological experiments to verify miRNA-disease associations. Therefore, expanding effective and outstanding computational methods to predict potential miRNA-disease associations is urgently needed and is attracting many computer scientists in recent years [7].

Recently, various computational methods to forecast possible miRNA-disease associations have been developed. For example, Liu et al. [8] proposed PBMDA prediction model which integrated known human miRNA-disease associations, miRNA functional similarity, disease semantic similarity and Gaussian interaction profile kernel similarity for miRNAs and diseases. They constructed a heterogeneous graph and further adopted depth-first search algorithm to figure out probable miRNA-disease associations. Chen et al. [9] presented a model called Graphlet Interaction for miRNA-Disease Association prediction (GIMDA) to predict miRNA-disease associations by measuring the graphlet interaction among miRNAs and among diseases. Graphlet is a type of subgraph with a few connections in a large network. GIMDA achieved a decisive performance but it was significantly time-consuming. Liang et al. [10] proposed a miRNA-disease association prediction method based on adaptive multi-view multi-label learning (AMVML). It learned a new affinity graph for miRNAs and diseases from multiple data sources. However, the integration of unreliable similarity matrices might weaken its overall prediction accuracy. The above mentioned methods for predicting miRNA-disease associations strongly relied on known human miRNA-disease associations. Most of existing methods need to use the similarity matrices such as the disease semantic similarity matrix and miRNA functional similarity matrix but they are not directly related to the miRNA-disease associations [11]. Besides, they have to deal with the problem of sparse similarity matrices which affected the prediction accuracies [12]. One other problem is that the miRNA-target interactions usually have a high rate of false-positive and false-negative [9, 13].

In fact, diseases are caused by the disturbance of a complex of interacting multiple biomolecules rather than the abnormity of a single biomolecule. The functionally dependent molecular components in human cells form a complex biological network, in which lncRNAs and proteins are important parts of human tissues and cells. It is the reason that some computational methods have recently based on multiple types of known associations or interactions among multiple objects to predict potential miRNA-disease associations. For example, Zhao et al. [7] developed a computational method based on a distance correlation set to predict miRNA-disease associations (DCSMDA) by integrating known lncRNA-disease associations, known miRNA-lncRNA associations, disease semantic similarity, and various lncRNA and disease similarity measures. DCSMDA did not require known miRNA-disease associations but it required the calculation of various similarity matrices and its performance depended on the pre-given threshold parameter. Mørk et al. [14] relied on known miRNA–protein associations and known protein–disease associations to infer miRNA–disease associations. Marissa Sumathipala and Weiss [15] integrated miRNA-gene, protein–protein, and gene-disease network information into a multi-level complex network to predict and prioritize biologically relevant miRNAs for diseases. Ji et al. [16] constructed a heterogeneous information network by integrating the known associations among lncRNAs, drugs, proteins, diseases, and miRNAs. They further employed the network embedding method which learned graph representations with global structural information to predict miRNA-disease associations. In general, the computational methods for predicting miRNA-disease associations based on multiple types of known associations among multiple objects are usually helpful for improving prediction accuracy. However, the number of known associations among biological objects is very limited in comparison with the number of objects in each type. Therefore, once again, these models have to be considered with the sparsity data problem.

In recent years, a variety of recommender systems have been developed to increase the association prediction reliability based on collaborative filtering methods. These methods rely on prior actions to predict user-item relationships to solve the problem of scarce known associations among different objects [17, 18]. Up to date, recommender algorithms have been appended into some computational models of prediction to identify different potential disease related biological objects. For example, Yu et al. [19] proposed a collaborative filtering model for lncRNA-disease association prediction based on the Naïve Bayesian classifier. Shen et al. [2] predicted miRNA-disease association with Collaborative Matrix Factorization model which caused bias to miRNAs with more known associated diseases. Li et al. [11] presented a collaborative filtering-based miRNA-disease association prediction model (CFMDA) to predict miRNA-disease association. CFMDA was straight and robust by considering a minimal amount of related information and no tunable parameters were defined. However, CFMDA’s association prediction performance was subjective because it only relies on miRNA-disease associations to execute predictions.

To solve the sparsity data problem and to take advantages of the integration of multiple types of known associations among multiple objects in improving prediction accuracy, in this paper, we proposed a new method to infer miRNA-disease associations using collaborative filtering and resource allocation algorithms on a tripartite graph. Our method is inspired by combining the collaborative filtering algorithm in CFNBC model introduced by Yu et al. [19] to solve the problem of imbalanced data and the method for association prediction established multiple types of known associations among multiple objects presented in TPGLDA model which introduced by Ding et al.[20] and the model in our former study [21]. Firstly, we constructed a tripartite graph which based on the known miRNA-disease associations, the known lncRNA-disease associations, and the known miRNA-lncRNA interactions. Secondly, we used a collaborative filtering algorithm to recommend miRNAs for lncRNAs and diseases, respectively. Next, we employed a resource allocation algorithm to infer miRNA-disease associations. Finally, we ranked all candidate miRNAs for each disease in descending order to suggest associations between miRNAs and diseases for further giving the evidence in the future. Our method achieved a trustworthy prediction performance under fivefold-cross-validation experiments with an Area Under Roc Curve (AUC) averaged value of 0.9788 and an Area Under Precision-Recall Curve (AUPR) averaged value of 0.9373. It is outperformance in comparison to several previous methods such as the DCSMDA [7] and the TPGLDA [20].

Methods

Materials

In this paper, we used datasets which came from the study of Zhao et al. [7]. We downloaded and used the Additional files 1, 2, 3, 4, and 5 from this study. These datasets contain 190 diseases, 111 lncRNAs and 264 miRNAs as described as follows:

Known lncRNA-miRNA associations

The known lncRNA-miRNA associations were collected from the starBasev2.0 [22] in February, 2017 and provided the most comprehensive experimentally confirmed lncRNA-miRNA interactions based on large-scale CLIP-Seq data. After eliminating duplicate values and erroneous data and also removing lncRNAs not included in DS2 dataset, we obtained the DS1 dataset which contains 1880 known lncRNA-miRNA associations.

Known lncRNA-disease associations

The known lncRNA-disease associations were collected from 8842 known disease-lncRNA associations in the MNDR database [23] and 2934 known disease-lncRNA associations in the LncRNADisease database [24]. After eliminating diseases without any MeSH descriptors because the disease names came from two different databases, merging the diseases with the same MeSH descriptors and removing the lncRNAs which were not included in the lncRNA-miRNA dataset (DS1), 936 known associations between diseases and lncRNAs (DS2) remained.

Known disease-miRNA associations

The known human miRNA-disease associations were downloaded from the HMDD V2.0 database [25]. This dataset (DS3) contains 3252 quality miRNA-disease associations after we eliminated the duplicate associations and miRNA-disease associations involving with other diseases or lncRNAs which were not contained in the DS1 or DS2 datasets.

Method overview

In this paper, we proposed a new method to infer miRNA-disease associations. The flowchart of the proposed method is illustrated in Fig. 1. Generally, our proposed method contains four main stages. At the first stage, we constructed a tripartite graph G0 based on known miRNA-disease associations, known lncRNA-disease associations, and known miRNA-lncRNA interactions. The tripartite graph G0 is represented by three adjacency matrices: A0MD, A0ML and A0DL where A0MD is the adjacency matrix between miRNAs and diseases, A0ML is the adjacency matrix between miRNAs and lncRNAs, A0DL is the adjacency matrix between diseases and lncRNAs. During the second stage, to solve the imbalance data problem, we employed a collaborative filtering algorithm on the tripartite graph G0 to obtain a tripartite graph Gu. The tripartite graph Gu is represented by three adjacency matrices: AuMD, AuML and A0DL where AuMD, AuML are the adjacency matrices obtained by updating A0MD and A0ML after using collaborative filtering algorithm. The tripartite graph Gu is used in a resource allocation algorithm at the third stage to calculate final resource score (Rscore_final) of miRNA candidates for each disease. At the final stage, we ranked all miRNA candidates’ Rscore_final for each disease in descending order so that the candidate with greater Rscore_final will have higher possibility to be verified in the future.

Fig. 1
figure 1

The flowchart of the proposed method

Construction of a tripartite graph G0

Inspired by previous studies [19, 20] to infer lncRNA-disease associations by using a tripartite graph, in this paper, we firstly construct a miRNA-disease-lncRNA tripartite graph G0 as follows:

Construction of known miRNA-disease association graph

Let M = {mk; k = 1,…,nm} denotes the set of miRNAs, D = {dj; j = 1,…, nd} denotes the set of diseases where nm, nd represent the number of miRNAs and diseases, respectively. We build a MD0 graph based on the known miRNA-disease associations. The MD0 graph is represented by a matrix A0MD which is the adjacency matrix of known miRNA-disease associations. The entity A0MD(mk, dj) is the element in kth row and jth column of A0MD, and A0MD(mk, dj) = 1 if miRNA mk is associated with disease dj, otherwise, A0MD(mk, dj) = 0.

Construction of known miRNA-lncRNA interaction graph

In the same way, let M = {mk; k = 1,…,nm} denotes the set of miRNAs, L = {li; i = 1,…, nl} denotes the set of lncRNAs where nm, nl represent number of miRNAs and lncRNAs, respectively. We can obtain ML0 graph and A0ML matrix. ML0 graph is built on known miRNA-lncRNA interactions. A0ML is the adjacency matrix of known miRNA-lncRNA interactions. The entity A0ML(mk, li) is the element in kth row and ith column of A0ML, and A0ML(mk, li) = 1 if miRNA mk interacts with lncRNA li, otherwise, A0ML(mk, li) = 0.

Construction of known disease-lncRNA association graph

Similarly, let D = {dj; j = 1,…, nd} denotes the set of diseases, L = {li; i = 1,…,nl} denotes the set of lncRNAs, where nd, nl represent number of diseases and lncRNAs, respectively. We can obtain DL0 graph and A0DL matrix where DL0 graph is built on known disease-lncRNA associations and A0DL is the adjacency matrix of known disease-lncRNA associations. The entity A0DL(dj, li) is the element in jth row and ith column of A0DL, and A0DL(dj, li) = 1 if disease dj is associated with lncRNA li, otherwise, A0DL(dj, li) = 0.

Construction of a tripartite graph G0

From the integration of the three MD0, ML0, DL0 graphs, we obtain a tripartite graph G0. The tripartite graph G0 is represented by three adjacency matrices: A0MD, A0ML and A0DL as mentioned before.

Construction of a tripartite graph G u

In the tripartite graph G0, the number of known associations between miRNAs and diseases as well as between miRNAs and lncRNAs are small. So that, for any given lncRNA node li and disease node dj, it is clear that the number of miRNA nodes which associated with both li and dj will be very small. To improve it, in our method, we use a collaborative filtering algorithm for recommending suitable miRNA nodes to corresponding lncRNA nodes and disease nodes, respectively. By considering that a recommender system may involve various input data including users and items [18], in our proposed method, we take lncRNAs and diseases as users, while miRNAs as items. For the two adjacency matrices A0ML and A0MD obtained above, it is easy for us to construct another adjacency matrix A0MLD = [A0ML, A0MD] by splicing A0ML and A0MD together because the number of rows in both A0ML and A0MD are same. It is clear that the row vector of A0MLD consists of the row vectors in A0ML and A0MD while the column vectors in A0MLD is the same as the column vectors in A0ML or A0MD.

On the basis of A0MLD and tripartite graph G0, we can obtain a co-occurrence matrix Rm x m, in which, the entity R(mk, mr) indicates the element in kth row and rth column of Rm x m where R(mk, mr) = 1 if and only if the miRNA mk and miRNA mr have at least one common neighboring node in G0, otherwise R(mk, mr) = 0. The common neighboring node can be an lncRNA or a disease in G0. So, a similarity matrix Rnor can be calculated by normalizing Rm x m as the following equation:

$${\mathrm{R}}^{nor}\left({m}_{k}, {m}_{r}\right)=\frac{\left|N\left({m}_{k}\right)\bigcap N({m}_{r})\right|}{\sqrt{\left|N({m}_{k})\right|*\left|N({m}_{r})\right|}}$$
(1)

where k, r are the number of miRNAs. \(\left|N\left({m}_{k}\right)\right|\) indicates the number of known lncRNAs and diseases associated to mk in G0, which means the number of elements with value equaling to 1 in kth row of A0MLD. \(\left|N\left({m}_{r}\right)\right|\) indicates the number of known lncRNAs and diseases associated to mr in G0, which means the number of elements with value equaling to 1 in rth row of A0MLD. ∣N(mk) ∩ N(mr)∣ indicates the number of known lncRNAs and diseases associated with both miRNA mk and miRNA mr simultaneously in G0.

Based on the similarity matrix Rnor and the adjacency matrix A0MLD, we calculate a new recommender matrix AuMLD as follows:

$$A^{u}_{MLD} = \, R^{nor} * \, A^{0}_{MLD}$$
(2)

Specifically, for a particular lncRNA li or disease dj in G0, if there is a miRNA mk satifying A0MLD(mk, li) = 1 or A0MLD(mk, dj) = 1 in A0MLD, then we firstly calculate the sum of the values of all elements in the ith or jth column in AuMLD, respectively. Therefore, we will have its averaged value P. Next, if the ith or jth column of AuMLD contains a miRNA \({m}_{\theta }\) which satisfies AuMLD(\({m}_{\theta }\), li) > P or AuMLD(\({m}_{\theta }\), dj) > P then we recommend miRNA \({m}_{\theta }\) for lncRNA li or disease dj, respectively. Also, we will add new edge between \({m}_{\theta }\) and li or \({m}_{\theta }\) and dj into the tripartite graph G0.

Finally, we obtain a tripartite graph Gu. The tripartite graph Gu contains three graphs: MDupdate, MLupdate and DL0 and can be represented by three adjacency matrices: AuMD, AuML and A0DL. MDupdate is the updated graph of MD0 after adding new edge between recommended miRNAs and diseases. MLupdate is the updated graph of ML0 after adding new edge between recommended miRNAs and lncRNAs. AuMD is the adjacency matrix which represents MDupdate graph. It contains 10,310 known and recommended associations and 39,850 unknown remained associations. AuML is the adjacency matrix which represents MLupdate graph.

Employing resource allocation process on the tripartite graph G u to infer miRNA-disease associations

To infer miRNA-disease association, we employ the resource allocation algorithm on the tripartite graph Gu as described in the following steps:

Step 1: Calculating resource allocation between miRNAs and diseases

For a specific miRNA mk, we define the initial resources located on disease dj as:

$$fd\left( {m_{k} } \right) = A^{u}_{MD} \left( {m_{k} , \, d_{j} } \right),\quad \, j = 1,2, \ldots ,n_{d}$$
(3)

where nd is the number of diseases.

Then we calculate the resource moved back from D to M by using a weight matrix W = {wkt}nm x nm to indicate the resource allocation process between miRNAs and diseases as follows:

$$w_{kt} = \frac{1}{{\deg A_{MD}^{u} \left( {m_{k} } \right)}}*\mathop \sum \limits_{j = 1}^{{n_{d} }} \frac{{A_{MD }^{u} \left( {m_{k} , d_{j} } \right) * A_{MD }^{u} \left( {m_{t} , d_{j} } \right)}}{{\deg A_{MD}^{u} \left( {d_{j} } \right)}}$$
(4)

where \({w}_{kt}\) is the contribution resource moved from tth node to kth node in M, and it can be understood as the similarity between miRNA mk and miRNA mt in MDupdate graph. \(\mathit{deg}{A}_{MD}^{u}\left({m}_{k}\right)\) is the degree of miRNA mk in MDupdate graph and it represents the number of associated diseases for miRNA mk. Similarly, \(\mathit{deg}{A}_{MD}^{u}\left({d}_{j}\right)\) is the degree of disease dj in MDupdate graph and it represents the number of associated miRNAs for disease dj.

With respect to previous study [20], we also modify the resource allocation algorithm by considering the level of consistency between the contribution of resource transferred in both directions. It shows the impact of co-selection (mk, mt) between the contribution of resource from mk to mt and the contribution of resource from mt to mk. A consistence-based resource allocation to represent a final miRNA-disease weight matrix W’ = {w’kt} can be defined as in the following equation:

$$W_{kt}^{^{\prime}} = W_{kt} + \frac{{W_{tk} }}{{\mathop \sum \nolimits_{s = 1}^{{n_{m} }} W_{sk} }}$$
(5)

From the combination of the final miRNA-disease weight matrix W’ and the adjacency matrix AuMD, we define a final resource Rscore_ondisease_1 located on D as follows:

$$Rscore\_ondisease\_1 = W^{{\prime }} *A^{u}_{MD}$$
(6)

Step 2: Calculating resource allocation between diseases and lncRNAs

In regard to resource allocation between genes and diseases in TPGLDA [20], the same initial resources located on M nodes are allocated from nodes in M to nodes in D and then moved back, and the final resource matrix Rscore_ondisease_2 located on D nodes are issued by:

$$Rscore\_ondisease\_2 = \mathop \sum \limits_{s = 1}^{{n_{l} }} \frac{{A_{DL }^{0} \left( {d_{j} , l_{s} } \right) }}{{\deg A_{DL}^{0} \left( {l_{i} } \right)}}*\mathop \sum \limits_{k = 1}^{{n_{d} }} \frac{{A_{MD}^{u} \left( { m_{k} , d_{j} } \right)}}{{\deg A_{DL}^{0} \left( {d_{j} } \right)}}$$
(7)

where \(\mathrm{deg}{A}_{DL}^{0}\left({l}_{i}\right)={\sum }_{j=1}^{{n}_{d}}{A}_{DL}^{0}({d}_{j}, {l}_{i})\) is the number of related diseases for lncRNA li or the degree of lncRNA li in DL0 graph. \(\mathrm{deg}{A}_{DL}^{0}\left({d}_{j}\right)\)=\({\sum }_{i=1}^{{n}_{l}}{A}_{DL}^{0}({d}_{j}, {l}_{i})\) is the number of related lncRNAs for disease dj or the degree of disease dj in DL0 graph.

Step 3: Calculating the final resource score Rscore_final to infer the potential disease-related miRNAs

We calculate the final resource score Rscore_final which is used to measure latent disease-related miRNAs as follows:

$$Rscore\_final = \gamma * Rscore\_ondisease\_1 + \, \left( {1 - \gamma } \right) \, *Rscore\_ondisease\_2$$
(8)

where γ is a tunable parameter with value in [0, 1]. Our model achieves the best prediction performance when γ = 0.9.

Ranking all candidate miRNAs’ Rscores for each disease in descending order

Finally, we sort all candidate miRNAs’ Rscore_final for each disease in descending order so that a higher score candidate will have more chances to be verified in the future.

Results

Performance measures

To evaluate our method performance in inferring miRNA-disease associations, we performed the fivefold-cross-validation experiments and evaluated the Area under roc curve (AUC) and the Area under precision-recall curve (AUPR) as described in following sections:

Evaluating the AUC under 5-fold-cross validation

After applying a collaborative filtering algorithm on tripatite graph G0, we obtained a tripartite graph Gu which contained three subgraphs: MDupdated graph, MLupdated graph and DL0 graph. By employing the resource allocation algorithm on the tripartite graph Gu, we predicted potential miRNA-disease associations. To evaluate our model performance in AUC term [26], we compared the inferred miRNA-disease associations resulted in Rscore_final matrix with the adjacency matrix AuMD of MDupdated graph. In MDupdated graph, we considered 10,310 associations of known and recommended associations as positive samples and the 39,850 remained unknown associations as negative samples. Then we randomly divided all positive and negative samples into 5 equal parts to perform fivefold-cross-validation. Next, in each running time, we used 4 parts of positive and negative samples for training and the remain part for testing. Our model is trained to recalculate Rscore_final in each running time. Basically, we computed the false positive rate (FPR) and true positive rate (TPR) with different γ values where FPR indicates the proportion of the real negative samples in predicted positive samples to all negative samples and TPR indicates the proportion of the real positive samples in all predicted positive samples. The FPR and TPR are calculated by the following equations:

$$FPR=\frac{FP}{FP+TN}$$
(9)
$$TPR= \frac{TP}{TP+FN}$$
(10)

where TP (true positive) means that a positive sample is correctly predicted as positive sample; FN (false negative) means that a positive sample is incorrectly predicted as negative sample; FP (false positive) indicates that a negative sample incorrectly predicted as positive sample; TN (true negative) indicates that a negative sample is correctly predicted as negative sample. We use TPR as vertical axis and FPR as horizontal axis to draw the receiver operating characteristic (ROC) curve [32], and the AUC value of our model achieves 0.9788 with γ = 0.9 after we perform the experiment for 10 times under fivefold-cross-validation. Figure 2 illustrates AUC curve with γ = 0.9 in one experimental running time.

Fig. 2
figure 2

AUC curve with γ = 0.9 in one experimental running time

Evaluate AUPR under 5-fold-cross validation

As previously mentioned, the data to evaluate our model performance is not balanced. Therefore, we also draw precision-recall curve and calculate the AUPR curve to evaluate prediction performance [27]. The Precision reflects the percentage of the accurately predicted positive samples in all predicted positive samples, and the Recall reflects the percentage of the accurately predicted positive samples in all real positive samples. We calculate Precision and Recall as follows:

$$Precision=\frac{TP}{TP+FP}$$
(11)
$$Recall=\frac{TP}{TP+FN}$$
(12)

After we perform the experiment under fivefold-cross-validation for 10 times, our model achieves the best AUPR value 0.9373 with γ = 0.9. Figure 3 illustrates AUPR curve with γ = 0.9 one experimental running time.

Fig. 3
figure 3

AUPR curve with γ = 0.9 in one experimental running time

Performance comparison with other related models

To demonstrate the outperformance of our model, we compare our model performance with the performance of DCSMDA method proposed by Zhao et al. [7]. We also implements predicting miRNA-disease associations by applying the resource allocation process introduced in [20] without applying collaborative filtering algorithm. The performances of these methods are shown in the Table 1.

Table 1 Performance comparison with other related models

As can be seen, our proposed method achieves better performance in comparison with DCSMDA and the method of applying TPGLDA in prediction of miRNA-disease associations for both AUC and AUPR values. Because of the sparsity data problem, AUC value usually achieves high score. However, in our proposed method, by using collaborative filtering algorithm to improve the density of miRNA-disease associations so that the updated adjacency matrix AuMD becomes more balanced which implies that the AUPR value (0.9373) could significantly be improved in comparison to AUPR value (0.7421) in case of applying TPGLDA model to predict miRNA-disease associations without using collaborative filtering algorithm. It demonstrates that our model achieves a more reliable performance than other previous methods.

Case studies

In addition to fivefold-cross-validation experiments, we also employed some case studies on our proposed model by sorting all candidate miRNAs for each disease. These predictions are utilized for further validation. In consistence with the previous study [20], all known and recommended miRNA-disease associations are considered as training samples, then the Rscore_final for each potential miRNA-disease association is calculated in sequence. Higher Rscore_final value indicates greater potential miRNA-disease association. In more detail, case studies on Prostatic Neoplasms, Heart Failure, Glioma and Open-angle Glaucoma are constructed to show the ability of our model in order to identify new disease-associated miRNAs.

Prostatic neoplasms, also known as Prostate Cancer, is the second-most prevalent type of cancers and the fifth-leading cause of cancer-related death in men [28]. miRNAs have been shown to play an important role in predicting prognosis of Prostate Cancer. Up to now, a variety of miRNAs have been reported to be associated with Prostatic Neoplasms /Prostate Cancer. For example, a target gene of miR-653-5p represses the proliferation and invasion of prostate cancer cells [29]. The dual action of miR-125b as a Tumor Suppressor and OncomiR-22 promotes Prostate Cancer tumorigenesis [30]. As shown in Table 2, there are 8 new miRNA-disease associations out of top predicted 40 miRNAs by applying our proposed method. All of new 8 miRNA-disease associations were confirmed by recent literatures.

Table 2 Top 40 predicted miRNAs for Prostatic Neoplasms

Heart failure (HF), also known as congestive heart failure (CHF) and congestive cardiac failure (CCF), is when the heart is unable to pump sufficiently to maintain blood flow to meet the body's needs. It is a widely prevalent syndrome imposing a significant burden of morbidity and mortality world-wide [31]. Unravelling the functional relevance of miRNAs within pathogenic pathways is a major challenge in cardiovascular research. Recently, a numerous miRNAs have been reported to be associated with heart failure. For instance, plasma miR-126 levels are up-regulated in HF patients [32]. MicroRNA-34 family members (miR-34a, -34b, and -34c) are up-regulated in the heart in response to stress [33]. Local microRNA-133a downregulation is associated with hypertrophy in the dyssynchronous heart [34]. Table 3 shows top 40 predicted heart failure related miRNAs by applying our proposed method. As can be seen, it contains 19 new miRNAs associated with Heart failure. All of these predicted associations were confirmed by literatures.

Table 3 Top 40 predicted miRNAs for Heart failure

Glioma is the most common central nervous system tumor and associated with poor prognosis. Identifying effective diagnostic biomarkers for glioma is particularly important in order to guide optimizing treatment [35]. Many studies have shown that some miRNAs are correlated with the diagnosis and prognosis of gliomas. For example, MiR-34a acts as tumor-suppressor by targeting many oncogenes related to proliferation, apoptosis, and invasion of gliomas [36]. MicroRNA (miR) 125b regulates cell growth and invasion in pediatric low grade glioma [37]. MicroRNA-21 promotes migration and invasion of glioma cells via activation of Sox2 and β-catenin signaling [38]. Therefore, in this study, we chose glioma as a case study to demonstrate our model’s ability in prediction associations between miRNAs and diseases. Table 4 lists top 40 glioma associated miRNAs inferred by our model. As illustrated, there are 14 new miRNAs associated with glioma, which are uncovered by applying our proposed method and all of them have been validated by literatures.

Table 4 Top 40 predicted miRNAs for Glioma

Glaucoma is the second leading cause of blindness in the United States of America [39]. The most common types of open-angle glaucoma (OAG) are primary open-angle glaucoma (POAG) and exfoliation glaucoma (XFG) [40]. Recent studies have shown that miRNAs may play a role in pathways implicated in glaucoma and act as biomarkers for disease pathogenesis [41]. In this paper, open-angle glaucoma is considered as an isolated disease because it is not associated with any miRNAs in the used datasets. However, our proposed method can be used to discover new associations for new diseases (or miRNAs) without any known associations before. As illustrated in Table 5, by applying our proposed method, 11 out of top 20 predicted open-angle glaucoma-related miRNAs have been confirmed by recent literatures.

Table 5 Top 20 miRNAs for GlaucomaOpen-Angle

Discussions

Although our proposed method achieved a reliable performance, it still exists some limitations which require further research. Firstly, our method still focuses on unweighted tripartite graph, so it may be improved by weighting the known lncRNA-disease associations, known miRNA-disease associations, and verified lncRNA-miRNA interactions. Secondly, enhancing the algorithm of appropriating resources can integrate the updated lncRNA-miRNA interactions into resource allocation process. Finally, the latest useful datasets should be collected to update our dataset library (Additional files 1, 2, 3, 4, 5).

Conclusion

In this paper, we proposed a new method to infer miRNA-disease associations using collaborative filtering and resource allocation on a miRNA-disease-lncRNA tripartite graph. By applying our proposed method, we can improve prediction accuracy, solve the sparsity data problem, and have not to use subjective and not directly related to association prediction information. The experimental results show that our method achieves a reliable performance with AUC and AUPR values 0.9788 and 0.9373, respectively, which is more impressive than several mentioned previously methods. It demonstrates the ability to infer new associations between miRNAs and diseases as indicated in case studies of Prostatic Neoplasms, Heart Failure, and Glioma diseases. Besides, it can discover new associations for new diseases (or miRNAs) without any known associations as indicated in the case study of Open-angle glaucoma disease. It suggests that our method can be considered as a powerful tool to predict miRNA-disease associations.

Availability of data and materials

The datasets used in our research were collected from Zhao et al.’s study https://doi.org/10.1186/s12859-018-2146-x [7].

Abbreviations

AUC:

Area under Roc curve

AUPR:

Area under precision-recall curve

FN:

False negative

FP:

False positive

FPR:

False positive rate

TP:

True positive

TPR:

True positive rate

lncRNA:

Long non-coding RNA

miRNA:

Micro RNA

OAG:

Open-angle glaucoma

POAG:

Primary open-angle glaucoma

References

  1. Ambros V. The functions of animal microRNAs. Nature. 2004;431(7006):350–5.

    Article  CAS  PubMed  Google Scholar 

  2. Shen Z, Zhang YH, Han K, Nandi AK, Honig B, Huang DS. MiRNA-disease association prediction with collaborative matrix factorization. Complexity. 2017;2017.

  3. Chen X, Xie D, Zhao Q, You ZH. MicroRNAs and complex diseases: from experimental results to computational models. Brief Bioinform. 2019;20(2):515–39.

    Article  CAS  PubMed  Google Scholar 

  4. Giannopoulou E, Alves P, Tewari AK, Gerstein MB. Epigenetic repression of miR-31 disrupts androgen receptor homeostasis and contributes to prostate cancer progression. Cancer Res. 2014;73(3):1232–44.

    Google Scholar 

  5. Masson S, Batkai S, Beermann J, Bär C, Pfanne A, Thum S, et al. Circulating microRNA-132 levels improve risk prediction for heart failure hospitalization in patients with chronic heart failure. Eur J Heart Fail. 2018;20(1):78–85.

    Article  CAS  PubMed  Google Scholar 

  6. Shi L, Cheng Z, Zhang J, Li R, Zhao P, Fu Z, et al. Hsa-Mir-181a and Hsa-Mir-181B function as tumor suppressors in human glioma cells. Brain Res. 2008;1236:185–93.

    Article  CAS  PubMed  Google Scholar 

  7. Zhao H, Kuang L, Wang L, Ping P, Xuan Z, Pei T, et al. Prediction of microRNA-disease associations based on distance correlation set. BMC Bioinform. 2018;19(1):1–14.

    Article  CAS  Google Scholar 

  8. Liu Y, Zeng X, He Z, Zou Q. Inferring MicroRNA-disease associations by random walk on a heterogeneous network with multiple data sources. IEEE/ACM Trans Comput Biol Bioinform. 2014;14(4):905–15.

    Article  Google Scholar 

  9. Chen X, Guan NN, Li JQ, Yan GY. GIMDA: graphlet interaction-based MiRNA-disease association prediction. J Cell Mol Med. 2018;22(3):1548–61.

    Article  CAS  PubMed  Google Scholar 

  10. Liang C, Yu S, Luo J. Adaptive multi-view multi-label learning for identifying disease-associated candidate miRNAs. PLoS Comput Biol. 2019;15(4):1–18.

    Article  Google Scholar 

  11. Li ZS, Liu B, Yan C. CFMDA: collaborative filtering-based MiRNA-disease association prediction. Multimed Tools Appl. 2019;78(1):605–18.

    Article  Google Scholar 

  12. Yu SP, Liang C, Xiao Q, Li GH, Ding PJ, Luo JW. MCLPMDA: a novel method for miRNA-disease association prediction based on matrix completion and label propagation. J Cell Mol Med. 2019;23(2):1427–38.

    Article  CAS  PubMed  Google Scholar 

  13. Chen X, Niu YW, Wang GH, Yan GY. MKRMDA: multiple kernel learning-based Kronecker regularized least squares for MiRNA-disease association prediction. J Transl Med. 2017;15(1):1–14. https://doi.org/10.1186/s12967-017-1340-3.

    Article  CAS  Google Scholar 

  14. Mørk S, Pletscher-frankild S, Caro AP, Gorodkin J, Jensen LJ. Protein-driven inference of miRNA-disease associations. Bioinformatics. 2014;30(3):392–7.

    Article  PubMed  Google Scholar 

  15. Sumathipala M, Weiss ST. Predicting miRNA-based disease-disease relationships through network diffusion on multi-omics biological data. Sci Rep. 2020;10(1):1–12. https://doi.org/10.1038/s41598-020-65633-6.

    Article  CAS  Google Scholar 

  16. Ji BY, You ZH, Cheng L, Zhou JR, Alghazzawi D, Li LP. Predicting miRNA-disease association from heterogeneous information network with GraRep embedding model. Sci Rep. 2020;10(1):1–12.

    Article  Google Scholar 

  17. Sarwar B, Karypis G, Konstan J. Item-based collaborative filtering recommendation algorithms. In: WWW ’01 Proceedings of the 10th international conference on world wide web. 2001;285–295.

  18. Liu NN, He L, Zhao M. Social temporal collaborative ranking for context aware movie recommendation. ACM Trans Intell Syst Technol. 2013;4(1).

  19. Yu J, Xuan Z, Feng X, Zou Q, Wang L. A novel collaborative filtering model for LncRNA-disease association prediction based on the Naïve Bayesian classifier. BMC Bioinform. 2019;20(1):1–13.

    Article  Google Scholar 

  20. Ding L, Wang M, Sun D, Li A. TPGLDA: novel prediction of associations between lncRNAs and diseases via lncRNA-disease-gene tripartite graph. Sci Rep. 2018;8(1):1–11.

    Google Scholar 

  21. Nguyen VT, Le TTK, Tran DH. A new method on lncRNA-disease-miRNA tripartite graph to predict lncRNA-disease associations. In: KSE2020 [Internet]. IEEE; 2020. p. 287–93. https://ieeexplore.ieee.org/document/9287563

  22. Li J, Liu S, Zhou H, Qu L, Yang J. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein—RNA interaction networks from large-scale CLIP-Seq data. 2014;42(December 2013):92–7. https://doi.org/10.1093/nar/gkt1248

  23. Cui T, Zhang L, Huang Y, Yi Y, Tan P, Zhao Y, et al. MNDR v2.0: an updated resource of ncRNA-disease associations in mammals. Nucl Acids Res. 2018;46(D1):D371–4. https://doi.org/10.1093/nar/gkx1025

  24. Bao Z, Yang Z, Huang Z, Zhou Y, Cui Q, Dong D. LncRNADisease 2.0: an updated database of long non-coding RNA-associated diseases. Nucl Acids Res. 2019;47(D1):D1034–7. https://doi.org/10.1093/nar/gky905

  25. Li Y, Qiu C, Tu J, Geng B, Yang J, Jiang T. HMDD v2.0: a database for experimentally supported human microRNA and disease associations. 2014;42(November 2013):1070–4. https://doi.org/10.1093/nar/gkt1023

  26. Hajian-Tilaki K et al. Receiver operating characteristic (ROC) curve analysis for medical diagnostic test evaluation. Casp J Intern Med 2013; 2013;4(2):627–35.

  27. Saito T, Rehmsmeier M. The precision-recall plot is more informative than the ROC plot when evaluating binary classifiers on imbalanced datasets. PLoS ONE. 2015;1–21.

  28. McGuire S. World Cancer Report 2014. Geneva, Switzerland: World Health Organization, International Agency for Research on Cancer, WHO Press, 2015. Adv Nutr. 2016;7(2):418–9.

  29. Fu Q, Sun Z, Yang F, Mao T, Gao Y, Wang H. SOX30, a target gene of miR-653-5p, represses the proliferation and invasion of prostate cancer cells through inhibition of Wnt/β-catenin signaling. Cell Mol Biol Lett. 2019;24(1):1–13.

    Article  CAS  Google Scholar 

  30. Budd WT, Seashols-Williams SJ, Clark GC, Weaver D, Calvert V, Petricoin E, et al. Dual action of miR-125b as a tumor suppressor and OncomiR-22 promotes prostate cancer tumorigenesis. PLoS ONE. 2015;10(11):1–21.

    Article  Google Scholar 

  31. Chen Y, Wang J, Sing K, Lee L, Wah O, Mark A. The association of heart failure-related microRNAs with neurohormonal signaling. BBA Mol Basis Dis. 2017;1863(8):2031–40. https://doi.org/10.1016/j.bbadis.2016.12.019.

    Article  CAS  Google Scholar 

  32. Wei XJ, Han M, Yang FY, Wei GC, Liang ZG, Yao H, et al. Biological significance of miR-126 expression in atrial fibrillation and heart failure. Braz J Med Biol Res. 2015;48:983–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Bernardo BC, Gao XM, Winbanks CE, Boey EJH, Tham YK, Kiriazis H, et al. Therapeutic inhibition of the miR-34 family attenuates pathological cardiac remodeling and improves heart function. Proc Natl Acad Sci U S A. 2012;109(43):17615–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. van Middendorp LB, Kuiper M, Munts C, Wouters P, Maessen JG, van Nieuwenhoven FA, et al. Local microRNA-133a downregulation is associated with hypertrophy in the dyssynchronous heart. ESC Hear Fail. 2017;4(3):241–51.

    Article  Google Scholar 

  35. Zhou Q, Liu J, Quan J, Liu W, Tan H, Li W. MicroRNAs as potential biomarkers for the diagnosis of glioma: a systematic review and meta-analysis. Cancer Sci. 2018;109(9):2651–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Vaitkiene P, Pranckeviciene A, Stakaitis R, Steponaitis G, Tamasauskas A, Bunevicius A. Association of miR-34a expression with quality of life of glioblastoma patients: a prospective study. Cancers (Basel). 2019;11(3):1–11.

    Article  Google Scholar 

  37. Yuan M, Da Silva ACAL, Arnold A, Okeke L, Ames H, Correa-Cerro LS, et al. MicroRNA (miR) 125b regulates cell growth and invasion in pediatric low grade glioma. Sci Rep. 2018;8(1):1–14.

    Google Scholar 

  38. Luo G, Luo W, Sun X, Lin J, Wang M, Zhang Y, et al. MicroRNA-21 promotes migration and invasion of glioma cells via activation of Sox2 and β-catenin signaling. Mol Med Rep. 2017;15(1):187–93.

    Article  CAS  PubMed  Google Scholar 

  39. Shaikh Y, Yu F, Coleman AL. Burden of undetected and untreated glaucoma in the United States. Am J Ophthalmol. 2014;158(6):1121–9. https://doi.org/10.1016/j.ajo.2014.08.023.

    Article  PubMed  Google Scholar 

  40. Drewry MD, Challa P, Kuchtey JG, Navarro I, Helwa I, Hu Y, et al. Differentially expressed microRNAs in the aqueous humor of patients with exfoliation glaucoma or primary open-angle glaucoma. Hum Mol Genet. 2018;27(7):1263–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Hindle AG, Thoonen R, Jasien JV, Grange RMH, Amin K, Wise J, et al. Identification of candidate miRNA biomarkers for glaucoma. Investig Ophthalmol Vis Sci. 2019;60(1):134–46.

    Article  CAS  Google Scholar 

  42. Qin W, Xie W, Yang X, Yang K, Zhou Q, Meng C. Down-regulation of miR-34a promotes the cell proliferation and inhibits apoptosis in glaucoma. Int J Clin Exp Pathol. 2016;9(2):1368–75.

    CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

About this supplement

This article has been published as part of BMC Medical Genomics Volume 14 Supplement 3 2021: Selected articles from the 19th Asia Pacific Bioinformatics Conference (APBC 2021): medical genomics The full contents of the supplement are available at https://bmcmedgenomics.biomedcentral.com/articles/supplements/volume-14-supplement-3.

Funding

Publication costs are funded by the Vietnam Ministry of Education and Training, project No. B2021-SPH-01. The funders did not play any role in the design of the study, the collection, analysis, and interpretation of data, or in writing of the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

VTN, TTKL, DHT, TQVN conceived and designed the study; VTN, DHT performed computational analyses; VTN, TTKL collected data and performed experiments. VTN wrote the first draft of the manuscript. All authors contributed to writing the paper, read, and approved the final manuscript.

Corresponding author

Correspondence to Dang Hung Tran.

Ethics declarations

Ethics approval and consent to participate

Not applicable. The study does not involve human subjects, only used public data.

Consent for publication

Not applicable.

Competing interests

The authors declare that there is no competing interest in relation to the publication of this article.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: For known lncRNA-disease associations.

Additional file 2

: For known lncRNA-disease associations.

Additional file 3

: For known lncRNA-disease associations.

Additional file 4

: For known lncRNA-miRNA associations.

Additional file 5

: For known disease-miRNA associations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nguyen, V.T., Le, T.T.K., Nguyen, T.Q.V. et al. Inferring miRNA-disease associations using collaborative filtering and resource allocation on a tripartite graph. BMC Med Genomics 14 (Suppl 3), 225 (2021). https://doi.org/10.1186/s12920-021-01078-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-021-01078-8

Keywords