 Methodology article
 Open Access
 Published:
NPCMF: Nearest Profilebased Collaborative Matrix Factorization method for predicting miRNAdisease associations
BMC Bioinformatics volume 20, Article number: 353 (2019)
Abstract
Background
Predicting meaningful miRNAdisease associations (MDAs) is costly. Therefore, an increasing number of researchers are beginning to focus on methods to predict potential MDAs. Thus, prediction methods with improved accuracy are under development. An efficient computational method is proposed to be crucial for predicting novel MDAs. For improved experimental productivity, large biological datasets are used by researchers. Although there are many effective and feasible methods to predict potential MDAs, the possibility remains that these methods are flawed.
Results
A simple and effective method, known as Nearest Profilebased Collaborative Matrix Factorization (NPCMF), is proposed to identify novel MDAs. The nearest profile is introduced to our method to achieve the highest AUC value compared with other advanced methods. For some miRNAs and diseases without any association, we use the nearest neighbour information to complete the prediction.
Conclusions
To evaluate the performance of our method, fivefold crossvalidation is used to calculate the AUC value. At the same time, three disease cases, gastric neoplasms, rectal neoplasms and colonic neoplasms, are used to predict novel MDAs on a goldstandard dataset. We predict the vast majority of known MDAs and some novel MDAs. Finally, the prediction accuracy of our method is determined to be better than that of other existing methods. Thus, the proposed prediction model can obtain reliable experimental results.
Background
MicroRNAs (miRNAs) are small noncoding RNAs whose length is generally 19 to 25 nt [1, 2]. In general, miRNAs regulate the expression of mRNA targets through a series of biological processes. However, the imbalance of miRNAs may have a serious impact on humans. Therefore, identifying novel miRNAdisease associations is important for treating complex genetic diseases [3, 4]. The first miRNA, lin4, was discovered in 1993. It is worth noting that lin4 is not the same as a conventional proteincoding gene; instead, lin4 encodes a 22nt regulatory RNA [5, 6]. In 2000, the second miRNA, let7, was discovered by researchers [7]. Since then, thousands of miRNAs have been discovered by biologists through a variety of biological and medical approaches. More than 2000 human miRNAs have been detected. Moreover, the latest version of the miRNA database miRBase contains 38,589 entries.
Recently, many biologists and medical scientists have found that miRNAs play an important role in different biological processes. In addition, an increasing number of miRNAs have been shown to be associated with cancer and other human diseases. For example, invasion and migration of breast cancer cells are inhibited by mir340 by targeting the oncoprotein cMet [8]. In addition, by targeting Cdc42 and Cdk6, mir137 inhibits the proliferation of lung cancer cells [9]. The progression of head and neck carcinomas is promoted by miR211 through the target TGFβR2 [10]. Moreover, in every paediatric brain tumour type, mir25, mir129, and mir142 are differentially expressed [11]. By identifying unknown potential miRNAdisease associations, the molecular mechanisms and pathogenesis of the disease can be elucidated.
In recent years, many researchers have employed computational methods associated with biomolecules and diseases [12,13,14,15]. In previous studies, an important assumption is that miRNAs with similar functions are more likely to be associated with diseases with similar phenotypes [16]. In other words, miRNAs with similar functions may be associated with the same disease. Increasingly effective methods and models are proposed for identifying novel miRNAdisease associations (MDAs). Chen et al. proposed a computational model named RLSMDA (Regularized Least Squares miRNADisease Association) based on semisupervised learning [17]. In this way, the problem of using negative MDAs is overcome. However, this semisupervised model is not perfect for the optimization of some parameters. Importantly, classifiers from the miRNA space and disease space are difficult to combine to predict novel MDAs. Chen et al. proposed a PathBased MiRNADisease Association (PBMDA) prediction model [15]. Specifically, a depthfirst search algorithm is used to predict novel MDAs on a heterogeneous graph consisting of three interlinked subgraphs. Chen et al. proposed a computational model named BNPMDA (Bipartite Network Projection for MiRNADisease Association) to obtain some valuable and reliable results [18]. The degree of preference between miRNA and disease is first described, then agglomerative hierarchical clustering is used, and finally, the BNPMDA method is implemented to predict potential MDAs. Jiang et al. constructed a model based on hypergeometric distribution through miRNA functional similarity, disease similarity and known MDA networks [19]. Then, these researchers analysed the actual effect in the prediction model. However, the shortcoming of this model is the excessive dependence on neighbouring miRNA data [20]. Chen et al. proposed a computational method to predict novel MDAs by using Laplacian regularized sparse subspace learning, and the accuracy of the prediction is improved [21]. Laplacian regularization is used to preserve the local structures. The strength of dimensionality reduction makes it easy to experiment with higherdimensional datasets. Shi et al. proposed a computational method to predict novel MDAs by performing a random walk algorithm [22]. Proteinprotein interactions (PPIs), miRNAtarget interactions and diseasegene associations were used to discover potential MDAs. This model is reliable, but there are still some shortcomings. The model strongly depended on the miRNAtarget interactions. Therefore, the final experimental results may have a high false positive rate or a high false negative rate [23]. Considering this disadvantage, Chen et al. developed a new method to solve this problem. The Random Walk with Restart for MiRNADisease Association (RWRMDA) model was used to map all miRNAs to a miRNA functional similarity network [24]. Mork et al. considered the protein information and proposed the miRPD method [25]. The method relies on proteindisease associations and proteinmiRNA associations to predict novel miRNAs and diseaserelated proteins. Chen et al. proposed an effective method, Heterogeneous Graph Inference MiRNADisease Association (HGIMDA), to predict novel MDAs [26]. In this method, Gaussian interaction profile (GIP) kernel similarity for diseases and miRNAs are integrated into the computational model. According to the final experimental results, this method improves the prediction accuracy. Chen et al. also proposed an effective method, Matrix Decomposition and Heterogeneous Graph Inference (MDHGI), to predict novel MDAs [14]. Among these approaches, the largest contribution is the combination of matrix decomposition and heterogeneous graph inference to predict new MDAs. In addition, Chen et al. proposed a method called inductive matrix completion [13]. The main measure is to complete the missing miRNAdisease association. Xuan et al. proposed an HDMP method based on weighting knearest neighbours [27]. Moreover, the semantic similarity and phenotypic similarity of the diseases were used to participate in the calculation of the functional similarity matrix of miRNAs. In contrast to previous studies, miRNAs of the same cluster have higher weights; therefore, they have the greatest potential to be associated with similar diseases when calculating the miRNA functional similarity matrix. Based on Xuan et al.’s method, Chen et al. proposed an improved method called RKNNMDA to identify potential MDAs [28]. Later, a valuable model named Matrix Completion for MiRNADisease Association prediction (MCMDA) was proposed by Li et al. [29]. However, this approach has certain limitations for new diseases and new miRNAs. These limitations lead to inaccuracies in the prediction results. Chen et al. developed a computational model named Ensemble Learning and Link Prediction for MiRNADisease Association (ELLPMDA) to identify potential MDAs [30]. Integrated similarity networks and integrated learning were used to predict novel MDAs. At the same time, this method is one of the more advanced methods. Chen et al. compiled the most advanced 20 prediction models to illustrate the importance of MDA prediction. Computational models have become an important means for novel MDA identification. The most important point is that the review can be inspired by more researchers [31].
In this paper, a simple but effective Nearest Profilebased Collaborative Matrix Factorization (NPCMF) method is proposed. This computational method can identify potential MDAs based on known MDAs. More importantly, unlike traditional matrix factorization models, considering that a new miRNA or a new disease is affected by their neighbour information when predicted, the nearest profile (NP) [32] is introduced to the CMF. The benefit of NP is that the nearest neighbour information for miRNA and disease is taken into account. The NP performs prediction through relatively reliable similarity functions. More precisely, the association profile of a new miRNA or disease is predicted using its similarities to other miRNAs or diseases, respectively; a new miRNA is one that has no known diseases, and similarly, a new disease is one that has no known interactions with any miRNAs. Notably, the existence of a large number of missing associations will have a negative impact on the final predictions. Weighted K Nearest Known Neighbours (WKNKN) is used as a preprocessing step to solve this problem [33]. Meanwhile, fivefold crossvalidation is performed to evaluate our experimental results. In addition, a simulation experiment is conducted to predict novel MDAs. Finally, the results demonstrate that our proposed method NPCMF is superior to other advanced methods.
The rest of this paper is organized as follows. Section 2 is first described, including our final experimental results and the goldstandard dataset used in this study. Section 3 contains the corresponding discussion. Section 4 contains conclusions for the full paper. Finally, Section 5 outlines our proposed method, specific solution steps and iterative processes.
Results
MDA dataset
The datasets used in the experiments were obtained from the human miRNAdisease database (HMDD), including 383 diseases, 495 miRNAs and 5430 human miRNAdisease associations [20]. The HMDD, which is a wellknown bioinformatics database, has collected thousands of miRNAdisease association pairs. Table 1 lists the specific information for the dataset.
In addition, the dataset contains three matrices: Y ∈ ℝ^{n × m}, S_{m} ∈ ℝ^{n × n} and S_{d} ∈ ℝ^{m × m}. The matrix Y is an adjacency matrix that is used to describe the associations between miRNAs and diseases. There are n miRNAs as rows and m diseases as columns. If miRNA M(i) is associated with disease d(j), the entity Y(M(i), d(j)) is 1; otherwise, it is 0. Moreover, this dataset is still a goldstandard dataset. The matrix Y is expressed as follows:
Performance evaluation metrics
To evaluate our approach, fivefold crossvalidation is conducted 100 times for each method. The known MDA dataset is randomly divided into 5 subsets, 4 of which are used as training sets, and the remaining subset is used as a testing set. It is worth noting that in our approach, WKNKN is used to eliminate unknown missing values. At the same time, the advantage is that the accuracy of the prediction can be improved to some extent.
In previous studies, the area under the curve (AUC) value is a reliable indicator of the evaluation method. Therefore, the AUC value is also used in this study. The area under the receiver operating characteristic (ROC) curve is considered to be the AUC. In general, the value of this area will not be greater than 1. The AUC values between 0.5 and 1 are reasonable. If the AUC is less than 0.5, the predicted results will be meaningless. In general, the ROC curve can be described in terms of true positive rate (TFR, sensitivity) and false positive rate (FPR, 1specificity). Thus, sensitivity and specificity (SPEC) can be expressed as follows:
where, according to the classification of the classifier, TP is the number of positive samples, FN is the number of false negative samples, and N is the number of negative samples. Similarly, TN is the number of negative samples, and FP is the number of false positive samples.
The MDA pairs are randomly removed in the input matrix Y before performing crossvalidation. This method is called CVp (CrossValidation pairs). Moreover, the purpose is to overcome the difficulty of prediction and accurately evaluate our method.
Comparison with other methods
In this study, the NPCMF method was compared with other advanced methods, CMF [34], HDMP [35], WBSMDA [36], HAMDA [37], and ELLPMDA [30]. Table 2 lists the experimental results with CVp. In Table 2, the final experimental results are expressed as the average of 100 fivefold crossvalidation. It is worth noting that AUC is known to be insensitive to skewed class distributions [38]. Considering that the dataset used in this paper is highly unbalanced, there are more negative factors than positive ones. Thus, AUC is a fair and reasonable evaluation indicator for all methods.
As listed in Table 2, the average AUCs of WBSMDA, HDMP, CMF, HAMDA, ELLPMDA, and NPCMF on the goldstandard dataset are 0.8185 ± 0.0009, 0.8342 ± 0.001, 0.8697 ± 0.0011, 0.8965 ± 0.0012, 0.9193 ± 0.0002 and 0.9429 ± 0.0011, respectively. The best value is in bold. Standard deviations are given in parentheses. From the above statistical results, our method achieved the highest AUC value, which was 12.46, 10.89, 7.34, 4.66, and 2.36% higher than WBSMDA, HDMP, CMF, HAMDA, and ELLPMDA, respectively. Compared with the CMF method, our method NPCMF has the best convergence. Furthermore, as shown in Fig. 1, the convergence analysis of CMF and NPCMF is shown by performing 100 iterations. Therefore, based on the above results, our proposed method is better than other existing advanced methods. Thus, the NPCMF method has proven to be effective and reliable. As shown in Fig. 2, in the fivefold crossvalidation experiment, the performance of each method can be demonstrated using the ROC curve.
Sensitivity analysis from WKNKN
Considering that there are some missing unknown associations in the matrix Y, WKNKN preprocessing is used to minimize the error. K represents the number of nearest known neighbours. p represents a decay term where p ≤ 1. These two parameters will be fixed to the optimal value before performing our method NPCMF. The sensitivities regarding K and p are represented by Figs. 3 and 4, respectively. The AUC tends to be stable when K = 5 and p = 0.7.
Comprehensive prediction for novel MDAs
A simulation experiment is conducted in this subsection. The simulation is conducted to obtain the final prediction score matrix. The specific process is divided into four steps. The first step is to execute our method; then, the two matrices A and B are obtained. The second step is to multiply A and B to obtain a predicted score matrix. The third step is to compare the predicted score matrix with the original MDAs matrix Y and the associations whose predicted score changes are filtered and sorted. The fourth step is to use the existing database to verify that our predicted associations are confirmed. Our method is applied to three disease cases, gastric neoplasms, rectal neoplasms and colonic neoplasms. These three diseases are more common among humans. Many miRNAs are closely related to these three diseases. Therefore, the final prediction results are more universal. In addition, the novel MDAs are validated by two popular miRNA disease databases, dbDEMC and miR2Disease.
The first case is gastric neoplasms. Despite a declining incidence [39], gastric neoplasms are a major cause of cancer death worldwide. Gonzalez et al. observed that gastric neoplasms constitute the second most frequent cancer in the world and the fourth most frequent cancer in Europe [40]. More information about the disease is published in http://www.omim.org/entry/613659. In the dataset used in the experiment, there are five MDAs associated with gastric neoplasms. After the simulation experiment is performed, three known associations are successfully predicted. At the same time, seven novel MDAs are predicted. More importantly, five of the seven novel MDAs have been confirmed by dbDEMC or miR2Disease. It is worth noting that miR214 is confirmed by both databases. For example, in 2011, when Oh et al. identified the biological validity of oncogenic miRNA microarray data for gastric neoplasms, miR214 in GC2 miRNAs was observed to be significantly upregulated [41]. In 2013, Lim et al. also found that miR214 is overexpressed in patients with gastric neoplasms compared with normal subjects [42]. It is worth noting that although both miR30b and miR296 are not confirmed by these two databases, they are still strongly associated with gastric neoplasms. Table 3 lists the detailed experimental results. The known associations are in bold.
The second case is rectal neoplasms. Fourteen known miRNAs were successfully predicted. Because there are more miRNAs associated with rectal neoplasms, we only selected the top 20 miRNAs with the highest correlation with rectal neoplasms. In Table 4, the miRNAs are arranged in descending order of the association score. Among the new miRNAs that are predicted, the fifteenth miRNA, miR196a, has the highest association score. Regarding miR196a, it was confirmed in the previous literature that it is associated with lymphoma [43]. Other researchers have found that miR196a is associated with prostate neoplasms [44]. Although the predicted novel MDAs are not confirmed by dbDEMC or miR2Disease, according to our experimental results, these MDAs are closely related to rectal neoplasms. Table 4 lists the detailed experimental results. The known associations are in bold.
The third case is colonic neoplasms. From the goldstandard dataset used in the experiment, there are more than 50 miRNAs related to colonic neoplasms; therefore, the top 50 are selected as the final prediction results according to the association score. Thirty known miRNAs are successfully predicted, and 20 new miRNAs are predicted. Of the 20 predicted new miRNAs, 12 are confirmed by dbDEMC and 8 are unconfirmed. For example, in 2009, Sarver et al. found that miR520 g was overexpressed in patients with colonic neoplasms compared with normal people according to a reliable biological experiment [43]. These researchers also found that miR204, miR206 and miR215 tend to be negatively expressed in colonic neoplasm patients. In addition, some unconfirmed miRNAs are sorted in descending order of association scores, including miR144, miR515, miR211, miR525, miR219, miR339, miR124 and miR340. Table 5 lists the detailed experimental results. The known associations are in bold.
Discussion
Based on the above experimental results, our proposed model NPCMF is superior to the most advanced methods overall. Moreover, although CMF is not as good as NPCMF, it has also achieved good experimental results. It is worth noting that our greatest contribution is to calculate the NP information for each disease and each miRNA to help predict potential MDAs. The shortcomings of CMF are that for new miRNAs and new diseases, the CMF method is unpredictable. However, NPCMF can achieve the prediction of new miRNAs and new diseases by using each miRNA and the nearest neighbour of the disease. Therefore, it is precisely because of the introduction of NP information that some novel MDAs can be predicted. By using NP information, we can obtain the best AUC value. Of course, this finding does not prove that NPCMF has no defects. One of the most obvious drawbacks for NPCMF is that excessive NP information is introduced, which may add additional noise while reducing prediction accuracy.
Conclusions
In this paper, a novel method based on nearest profile collaborative matrix factorization is developed for predicting novel MDAs. When novel MDAs are predicted, the nearest neighbour information for miRNAs and diseases is fully considered. In addition, incorporating the Gaussian interaction profile kernels of miRNAs and diseases also contributed to the improvement of prediction performance. The AUC value is used as a reliable indicator to evaluate our method. In addition, due to technical limitations, we have not used the latest version of the dataset, such as HMDD V3.0; therefore, we will attempt to use the latest dataset for future experiments.
In the future, more effective methods may be used to predict new MDAs. More differentially expressed miRNAs associated with the disease will be identified. At the same time, increasing numbers of valuable datasets are being published by online bioinformatics databases. Thus, more datasets can be tested by researchers. Importantly, NPCMF may be helpful for novel MDA prediction and relevant miRNA research from computational biology.
Methods
Our goal is to develop a matrix factorization method that can predict novel MDAs based on known MDAs. First, a matrix factorization model is constructed to represent the correlation between miRNAs and diseases. Next, the Gaussian interaction profile kernels of miRNA and disease are expressed as their network information. Then, the nearest profile of miRNAs and diseases are obtained. Finally, a prediction score matrix is obtained by multiplying two low rank matrices.
MiRNA functional similarity
Wang et al. developed a method named MISIM for calculating the similarity scores of miRNA functions [45]. Moreover, the dataset that we used is downloaded from the website http://www.cuilab.cn/files/images/cuilab/misim.zip. Then, matrix S_{m} represents the functional similarity matrix of the miRNAs. Since the selfsimilarity of a miRNA is 1, in the matrix S_{m}, the elements on the diagonal are all 1.
Disease semantic similarity
In previous studies, directed acyclic graphs (DAGs) have been used by many researchers to describe diseases. From the National Library of Medicine (http://www.nlm.nih.Gov/), a variety of disease relationships based on the disease DAG can be obtained from the MeSH descriptor of Category C. DAG(DD) = (d, T(DD), E(DD)) is used to describe disease DD. T(DD) is the node set and E(DD) is the corresponding link set. The DD in DAG(DD) formula is defined as
where Δ represents the semantic contribution factor. In this work, based on previous literature [45], the value of Δ is set to 0.5.
In addition, matrix S_{d} represents the semantic similarity matrix of the disease. Similarly, in the matrix S_{d}, the elements on the diagonal are all 1. It is worth noting that if the two diseases d(i) and d(j) have a larger common part of the DAGs, these two diseases will have higher semantic similarity values. The semantic similarity score between two diseases is defined as follows:
Gaussian interaction profile kernel similarity
The method is based on the following assumption. The topological structure of the known MDA network is represented by Gaussian interaction profile kernel similarity [46]. M(i) and M(j) are two miRNAs, and d(i) and d(j) are two diseases. Therefore, the network similarity calculations can be written as
where γ is expressed as a parameter that adjusts the bandwidth of the kernel. In principle, the setting of γ should be implemented by crossvalidation, but according to a previous study [47], γ is simply set to 1. In addition, the interaction profiles of M_{i} and M_{j} can be represented as Y(M_{i}) and Y(M_{j}), respectively. Similarly, the interaction profiles of d_{i} and d_{j} can be represented as Y(d_{i}) and Y(d_{j}), respectively. Thus, the miRNA network similarity matrix can be combined by S_{m} into K_{m}, and the disease network similarity matrix can be combined by S_{d} into K_{d}. The calculation formulas are as follows:
where α ∈ [0, 1] is an adjustable parameter. We perform a sensitivity analysis on α. When α = 0.5, the highest AUC value can be obtained. Figure 5 shows the sensitivity analysis for α. K_{m} is a miRNA kernel matrix, which represents a linear combination of the miRNA functional similarity matrix S_{m} and the miRNA network similarity matrix GIP_{m}. Similarly, K_{d} is similar to K_{m}. K_{d} is a disease kernel matrix. In each crossvalidation, we recalculate the miRNA Gaussian similarity and disease Gaussian similarity. Specifically, the miRNA Gaussian similarity matrix and the disease Gaussian similarity matrix are obtained from a known MDA matrix. Therefore, we ensure that the Gaussian similarity is recalculated each time the crossvalidation is performed so that the Gaussian similarity correctly reflects the characteristics of the MDA matrix.
NPCMF for MDA prediction
The traditional CMF is a reliable method for predicting novel MDAs [34]. Collaborative filtering is introduced to CMF. The objective function of CMF is defined as
where ‖⋅‖_{F} is the Frobenius norm, and λ_{l}, λ_{d} and λ_{t} are nonnegative parameters. It is worth noting that the three parameters are set on the training set by performing crossvalidation. A grid search is used to obtain the optimal parameters from these values: λ_{l} ∈ {2^{−2}, 2^{−1}, 2^{0}, 2^{1}}, λ_{d}/λ_{l} ∈ {0, 10^{−4}, 10^{−3}, 10^{−2}, 10^{−1}}. The MDA matrix Y is decomposed into two matrices A and B, where AB^{T} ≈ Y. The NPCMF method uses regularization terms to request that the potential feature vectors of similar miRNAs and similar diseases are similar, and the potential feature vectors of dissimilar miRNAs and dissimilar diseases are dissimilar, respectively [33]. In this instance, S_{m} ≈ AA^{T} and S_{d} ≈ BB^{T}.
However, the CMF method ignores the network information of miRNAs and diseases. Therefore, GIP is introduced to the CMF [48]. Therefore, K_{m} and K_{d} are substituted into the objective function and written as
Then, the objective function is further written as
More importantly, when predicting novel MDAs, the nearest neighbour information will affect the final results. Therefore, the nearest profile (NP) is introduced to the CMF. For example, the NP for a new miRNA M(i) is computed as
where M_{nearest} is the miRNA most similar to M_{i}, and Y_{NP}(M_{i}) is the association profile of miRNA M_{i}. The NP for a new disease d_{i} is computed as
where d_{nearest} is the disease most similar to d_{i}, and Y_{NP}(d_{i}) is the association profile of disease d_{i}.
The NP process can be performed in four steps. First, the selfsimilarity of the matrices K_{m} and K_{d} is removed. Next, the nearest neighbour of each miRNA and disease is obtained. Then, all miRNA similarities and disease similarities are reset to 0. Finally, the nearest neighbour matrix N_{m} of the K_{m}based miRNA is obtained. In the previous study [49], the definition of the nearest neighbour matrix is given. According to Eq. (14), we can obtain N_{m} = arg max K_{m}(M_{i}). Simultaneously, the nearest neighbour matrix N_{d} of the K_{d}based disease is also obtained. According to Eq. (15), we can obtain N_{d} = arg max K_{d}(d_{i}). Based on objective function (11), the objective function of NPCMF can be written as follows:
where ‖⋅‖_{F} is the Frobenius norm, and λ_{l}, λ_{d} and λ_{t} are nonnegative parameters. The first term is an approximate model of the matrix Y. In the second term, the Tikhonov regularization is used to minimize the norms of A, B. The last two regularization terms minimize the squared error between N_{m} (N_{d}) and AA^{T} (BB^{T}).
Initialization of A and B
For the input MDAs matrix, A and B are initialized by the singular value decomposition (SVD) method. The initialization formula can be written as follows:
where S_{k} is a diagonal matrix, which contains the k largest singular values.
Optimization
Considering that the least squares method is an effective way to update A and B, in this paper, the least squares method is used to update A and B. A and B are updated until convergence. L is represented as the objection function of the NPCMF method. Then, A and B are respectively subjected to partial derivatives. ∂L/∂A and ∂L/∂B are both set to 0. In addition, λ_{l}, λ_{d} and λ_{t} are automatically determined optimal parameter values by the fivefold crossvalidation. The update rules are as follows:
Therefore, the specific algorithm of NPCMF is as follows:
Availability of data and materials
The datasets that support the findings of this study are available in https://github.com/cuizhensdws/npcmf.
Abbreviations
 CMF:

Collaborative matrix factorization method
 CV:

Crossvalidation
 NPCMF:

Nearest Profilebased Collaborative Matrix Factorization
 SVD:

Singular value decomposition
 WKNKN:

Weighted K Nearest Known Neighbours
References
 1.
Ambros V. microRNAs: tiny regulators with great potential. Cell. 2001;107(7):823–6.
 2.
Ambros V. The functions of animal microRNAs. Nature. 2004;431(7006):350.
 3.
Zheng CH, Huang DS, Zhang L, Kong XZ. Tumor clustering using nonnegative matrix factorization with gene selection. IEEE Trans Inf Technol Biomed. 2009;13(4):599–607.
 4.
Sethupathy P, Collins FS. MicroRNA target site polymorphisms and human disease. Trends Genet. 2008;24(10):489–97.
 5.
Lee RC, Feinbaum RL, Ambros V. The C. elegans heterochronic gene lin4 encodes small RNAs with antisense complementarity to lin14. Cell. 1993;75(5):843.
 6.
Wightman B, Ha I, Ruvkun G. Posttranscriptional regulation of the heterochronic gene lin14 by lin4 mediates temporal pattern formation in C. elegans. Cell. 1993;75(5):855–62.
 7.
Reinhart BJ, Slack FJ, Basson M, Pasquinelli AE, Bettinger JC, Rougvie AE, Horvitz HR, Ruvkun G. The 21nucleotide let7 RNA regulates developmental timing in Caenorhabditis elegans. Nature. 2000;403(6772):901–6.
 8.
Wu ZS, Wu Q, Wang CQ, Wang XN, Huang J, Zhao JJ, Mao SS, Zhang GH, Xu XC, Zhang N. miR340 inhibition of breast cancer cell migration and invasion through targeting of oncoprotein cmet. Cancer. 2011;117(13):2842–52.
 9.
Zhu X, Li Y, Shen H, Li H, Long L, Hui L, Xu W. miR137 inhibits the proliferation of lung cancer cells by targeting Cdc42 and Cdk6. FEBS Lett. 2013;587(1):73–81.
 10.
Chu TH, Yang CC, Liu CJ, Lui MT, Lin SC, Chang KW. miR211 promotes the progression of head and neck carcinomas by targeting TGFβRII. Cancer Lett. 2013;337(1):115–24.
 11.
Patel V, Williams D, Hajarnis S, Hunter R, Pontoglio M, Somlo S, Igarashi P. miR17~92 miRNA cluster promotes kidney cyst growth in polycystic kidney disease. Pnas. 2013;110(26):10765–70.
 12.
Chen X, Yan CC, Zhang X, You ZH. Long noncoding RNAs and complex diseases: from experimental results to computational models. Brief Bioinform. 2016;18(4):558–76.
 13.
Chen X, Wang L, Qu J, Guan NN, Li JQ. Predicting miRNA–disease association based on inductive matrix completion. Bioinformatics. 2018;34(24):4256–65.
 14.
Chen X, Yin J, Qu J, Huang L. MDHGI: matrix decomposition and heterogeneous graph inference for miRNAdisease association prediction. PLoS Comput Biol. 2018;14(8):e1006418.
 15.
You ZH, Huang ZA, Zhu Z, Yan GY, Li ZW, Wen Z, Chen X. PBMDA: a novel and effective pathbased computational model for miRNAdisease association prediction. PLoS Comput Biol. 2017;13(3):e1005455.
 16.
Pasquier C, Gardès J. Prediction of miRNAdisease associations with a vector space model. Sci Rep. 2016;6:27036.
 17.
Chen X, Yan GY. Semisupervised learning for potential human microRNAdisease associations inference. Sci Rep. 2014;4:5501.
 18.
Chen X, Xie D, Wang L, Zhao Q, You ZH, Liu H. BNPMDA: bipartite network projection for MiRNA–disease association prediction. Bioinformatics. 2018;34(18):3178–86.
 19.
Jiang Q, Hao Y, Wang G, Juan L, Zhang T, Teng M, Liu Y, Wang Y. Prioritization of disease microRNAs through a human phenomemicroRNAome network. BMC Syst Biol. 2010;4(Suppl 1):S2.
 20.
Chen X, Gong Y, Zhang DH, You ZH, Li ZW. DRMDA: deep representationsbased miRNA–disease association prediction. J Cell Mol Med. 2018;22(1):472–85.
 21.
Chen X, Huang L. LRSSLMDA: Laplacian regularized sparse subspace learning for MiRNAdisease association prediction. PLoS Comput Biol. 2017;13(12):e1005912.
 22.
Shi H, Xu J, Zhang G, Xu L, Li C, Wang L, Zhao Z, Jiang W, Guo Z, Li X. Walking the interactome to identify human miRNAdisease associations through the functional link between miRNA targets and disease genes. BMC Syst Biol. 2013;7(1):101.
 23.
Chen X, Niu YW, Wang GH, Yan GY. MKRMDA: multiple kernel learningbased Kronecker regularized least squares for MiRNA–disease association prediction. J Transl Med. 2017;15(1):251.
 24.
Chen X, Liu MX, Yan GY. RWRMDA: predicting novel human microRNAdisease associations. Mol BioSyst. 2012;8(10):2792–8.
 25.
Mørk S, Pletscherfrankild S, Palleja CA, Gorodkin J, Jensen LJ. Proteindriven inference of miRNAdisease associations. Bioinformatics. 2014;30(3):392.
 26.
Chen X, Yan CC, Zhang X, You ZH, Huang YA, Yan GY. HGIMDA: heterogeneous graph inference for miRNAdisease association prediction. Oncotarget. 2016;7(40):65257–69.
 27.
Xuan P, Han K, Guo M, Guo Y, Li J, Ding J, Liu Y, Dai Q, Li J, Teng Z. Prediction of microRNAs associated with human diseases based on weighted k Most similar neighbors. PLoS One. 2013;8(9):e70204.
 28.
Chen X, Wu QF, Yan GY. RKNNMDA: rankingbased KNN for MiRNAdisease association prediction. RNA Biol. 2017;14(7):952–62.
 29.
Li JQ, Rong ZH, Chen X, Yan GY, You ZH. MCMDA: matrix completion for MiRNAdisease association prediction. Oncotarget. 2017;8(13):21187–99.
 30.
Chen X, Zhou Z, Zhao Y. ELLPMDA: ensemble learning and link prediction for miRNAdisease association prediction. RNA Biol. 2018;15(6):807–18.
 31.
Chen X, Xie D, Zhao Q, You ZH. MicroRNAs and complex diseases: from experimental results to computational models. Brief Bioinform. 2017;20(2):515–39.
 32.
Yamanishi Y, Araki M, Gutteridge A, Honda W, Kanehisa M. Prediction of drug–target interaction networks from the integration of chemical and genomic spaces. Bioinformatics. 2008;24(13):i232–40.
 33.
Ezzat A, Zhao P, Wu M, Li XL, Kwoh CK. Drugtarget interaction prediction with graph regularized matrix factorization. IEEE/ACM Trans Comput Biol Bioinformatics. 2017;14(3):646–56.
 34.
Shen Z, Zhang YH, Han K, Nandi AK, Honig B, Huang DS. miRNAdisease association prediction with collaborative matrix factorization. Complexity. 2017;2017(9):1–9.
 35.
Lucherini OM, Obici L, Ferracin M, Fulci V, Mcdermott MF, Merlini G, Muscari I, Magnotti F, Dickie LJ, Galeazzi M. Correction: first report of circulating MicroRNAs in tumour necrosis factor receptorassociated periodic syndrome (TRAPS). PLoS One. 2013;8(9):e73443.
 36.
Chen X, Yan CC, Zhang X, You ZH, Deng L, Liu Y, Zhang Y, Dai Q. WBSMDA: within and between score for MiRNAdisease association prediction. Sci Rep. 2016;6:21106.
 37.
Chen X, Niu YW, Wang GH, Yan GY. HAMDA: hybrid approach for MiRNAdisease association prediction. J Biomed Inform. 2017;76:50–8.
 38.
Ezzat A, Wu M, Li XL, Kwoh CK. Drugtarget interaction prediction via class imbalanceaware ensemble learning. Bmc Bioinformatics. 2016;17(19):267–76.
 39.
Howson CP, Hiyama T, Wynder EL. The decline in gastric cancer: epidemiology of an unplanned triumph. Epidemiol Rev. 1986;8(1):1–27.
 40.
González CA, Sala N, Capellá G. Genetic susceptibility and gastric cancer risk. Int J Cancer. 2010;100(3):249–60.
 41.
Oh HK, Tan AL, Das K, Ooi CH, Deng NT, Tan IB, Beillard E, Lee J, Ramnarayanan K, Rha SY. Genomic loss of miR486 regulates tumor progression and the OLFM4 antiapoptotic factor in gastric cancer. Clin Can Res. 2011;17(9):2657–67.
 42.
Lim JY, Yoon SO, Seol SY, Hong SW, Kim JW, Choi SH, Lee JS, Cho JY. Overexpression of miR196b and HOXA10 characterize a poorprognosis gastric cancer subtype. World J Gastroenterol. 2013;19(41):7078–88.
 43.
Sarver AL, French AJ, Borralho PM, Thayanithy V, Oberg AL, Silverstein KA, Morlan BW, Riska SM, Boardman LA, Cunningham JM. Human colon cancer profiles show differential microRNA expression depending on mismatch repair status and are characteristic of undifferentiated proliferative states. BMC Cancer. 2009;9(1):401.
 44.
Taylor BS, Schultz N, Hieronymus H, Gopalan A, Xiao Y, Carver BS, Arora VK, Kaushik P, Cerami E, Reva B. Integrative genomic profiling of human prostate Cancer. Cancer Cell. 2010;18(1):11–22.
 45.
Wang D, Wang J, Lu M, Song F, Cui Q. Inferring the human microRNA functional similarity and functional network based on microRNAassociated diseases. Bioinformatics. 2010;26(13):1644–50.
 46.
Chen X, Huang YA, You ZH, Yan GY, Wang XS. A novel approach based on KATZ measure to predict associations of human microbiota with noninfectious diseases. Bioinformatics. 2016;33(5):733–9.
 47.
van Laarhoven T, Nabuurs SB, Marchiori E. Gaussian interaction profile kernels for predicting drug–target interaction. Bioinformatics. 2011;27(21):3036–43.
 48.
Cui Z, Gao YL, Liu JX, Wang J, Shang J, Dai LY. The computational prediction of drugdisease interactions using the dualnetwork L 2, 1CMF method. BMC bioinformatics. 2019;20(1):5.
 49.
Ding H, Takigawa I, Mamitsuka H, Zhu S. Similaritybased machine learning methods for predicting drug–target interactions: a brief review. Brief Bioinform. 2013;15(5):734–47.
Acknowledgements
Thanks go to the editor and the anonymous reviewers for their comments and suggestions.
Funding
This work was supported in part by the NSFC under Grant Nos. 61872220, 61873001, and 61572284. The funder played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Author information
Affiliations
Contributions
YLG and ZC jointly contributed to the design of the study. YLG designed and implemented the NPCMF method, performed the experiments, and drafted the manuscript. JXL gave statistical and computational advice for the project and participated in designing evaluation criteria. JW and CHZ contributed to the data analysis. All authors read and approved the final manuscript.
Corresponding author
Correspondence to JinXing Liu.
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
About this article
Cite this article
Gao, Y., Cui, Z., Liu, J. et al. NPCMF: Nearest Profilebased Collaborative Matrix Factorization method for predicting miRNAdisease associations. BMC Bioinformatics 20, 353 (2019). https://doi.org/10.1186/s1285901929565
Received:
Accepted:
Published:
Keywords
 MiRNAdisease association prediction
 Nearest profile
 Gaussian interaction profile
 Matrix factorization