Next Article in Journal
Near-Hexaploid and Near-Tetraploid Aneuploid Progenies Derived from Backcrossing Tetraploid Parents Hibiscus syriacus × (H. syriacus × H. paramutabilis)
Previous Article in Journal
H1299R Variant in Factor V and Recurrent Pregnancy Loss: A Systematic Review and Meta-Analysis Protocol
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MDSCMF: Matrix Decomposition and Similarity-Constrained Matrix Factorization for miRNA–Disease Association Prediction

1
Network Information Center, Qufu Normal University, Qufu 273165, China
2
School of Cyber Science and Engineering, Qufu Normal University, Qufu 273165, China
3
School of Artifial Intelligence, Anhui University, Hefei 230601, China
*
Authors to whom correspondence should be addressed.
Genes 2022, 13(6), 1021; https://doi.org/10.3390/genes13061021
Submission received: 1 May 2022 / Revised: 1 June 2022 / Accepted: 2 June 2022 / Published: 6 June 2022

Abstract

:
MicroRNAs (miRNAs) are small non-coding RNAs that are related to a number of complicated biological processes, and numerous studies have demonstrated that miRNAs are closely associated with many human diseases. In this study, we present a matrix decomposition and similarity-constrained matrix factorization (MDSCMF) to predict potential miRNA–disease associations. First of all, we utilized a matrix decomposition (MD) algorithm to get rid of outliers from the miRNA–disease association matrix. Then, miRNA similarity was determined by utilizing similarity kernel fusion (SKF) to integrate miRNA function similarity and Gaussian interaction profile (GIP) kernel similarity, and disease similarity was determined by utilizing SKF to integrate disease semantic similarity and GIP kernel similarity. Furthermore, we added L2 regularization terms and similarity constraint terms to non-negative matrix factorization to form a similarity-constrained matrix factorization (SCMF) algorithm, which was applied to make prediction. MDSCMF achieved AUC values of 0.9488, 0.9540, and 0.8672 based on fivefold cross-validation (5-CV), global leave-one-out cross-validation (global LOOCV), and local leave-one-out cross-validation (local LOOCV), respectively. Case studies on three common human diseases were also implemented to demonstrate the prediction ability of MDSCMF. All experimental results confirmed that MDSCMF was effective in predicting underlying associations between miRNAs and diseases.

1. Introduction

MiRNAs are 17–24 nt non-coding RNAs that play a pivotal role in controlling the expression of genes through RNA cleavage or translation repression [1,2,3]. Lin-4 was the first miRNA inspected experimentally, by Lee et al. [4] in 1993. Since that time, a large number of miRNAs have been discovered experimentally by researchers [4,5,6]. Researchers have found that various miRNAs are bound up with several crucial biological processes, such as cell development, cell differentiation, cell proliferation, etc. [7,8,9,10]. Developmental defects can be the result of the dysregulation of miRNAs that are associated with the progression of diseases [11]. In the meantime, numerous studies have indicated that miRNAs are connected with a serious of human neoplasms, including breast neoplasms, lung neoplasms, prostate neoplasms, etc. [12,13,14]. Hence, distinguishing miRNAs associated with diseases can deepen the understanding of the genetic causes of complex diseases. Strong connections between miRNAs and diseases have been found by a variety of traditional experiments in the past few years [15,16]. Traditional manual models can infer the associations between miRNAs and diseases, but these are time-consuming, laborious, and have a high failure rate. Therefore, showing the potential relationships between miRNAs and diseases requires effective and stable computational methods, which can obtain increasingly reliable miRNA–disease associations.
In the past, a great number of heterogeneous-network-based algorithms and methods have been applied to predict potential miRNA–disease relationships [17,18,19]. Under the assumption that miRNAs with similar functions have a high probability of being related to diseases with similar phenotypes, and vice versa [20], Jiang et al. [21] established a new calculation-based model that identified potential miRNA–disease connections by applying hypergeometric distribution. However, the similarity information utilized in this model excluded similarity scores. Li et al. [22] constructed a new model that could be used to prioritize human cancer miRNAs by measuring the associations between cancer and miRNAs based on the functional consistency scores of the miRNA target genes and the cancer-related genes. To infer the miRNA–protein connections and disease–protein connections, Mørk et al. [23] built the miPRD model. This model used selective connections to predict the relationships between miRNAs and diseases. Chen et al. [24] utilized the within and between scores of each miRNA–disease combination in the WBSMDA model to predict underlying miRNAs related to diseases. The WBSMDA model also predicted the possible relationships between new diseases and new miRNAs. Yu et al. [25] proposed an identifiable model to infer potential miRNA–disease relationships. This model combined miRNA functional similarity, disease semantic similarity and disease phenotypic similarity to create a modified information flow method. In a phenome–microRNAome network, possible connections and validated relationships between miRNAs and diseases were adopted. Chen et al. [26] introduced the Jaccard similarity between miRNAs and diseases into the BLHARMDA model to investigate prospective miRNA–disease relationships. For improving the prediction efficiency, BLHARMDA used a bipartite local model with a KNN architecture. Ha et al. [27] proposed a computational framework of metric learning named MLMD for predicting potential miRNA–disease associations. MLMD exploited distance metric learning on a miRNA–disease bipartite graph to infer unconfirmed miRNA–disease associations. The excellent performance of MLMD could be attributed to two factors: On the one hand, the implementation of metric learning overcame the violation of triangle inequality. On the other hand, the miRNA expression data were adequately trained in metric learning. Li et al. [28] proposed a similarity-constrained matrix factorization method to infer unconfirmed disease-associated miRNAs. To construct an information-rich similarity matrix, they utilized similarity network fusion to integrate various kinds of similarities. Then, similarity-based regularization terms were added to common non-negative matrix factorization to form a similarity-constrained matrix factorization algorithm, which was applied to make accuracy predictions. The above methods are mainly based on the construction of heterogeneous networks to identify and speculate on the potential disease-related miRNAs, and after cross-validation and case analysis experiments, it was proven that they can be used to observe the potential association between miRNA and disease, but their prediction performance still needs to be improved.
Recently, methods based on the random walk method have gradually been proposed, and more accurate prediction results have been obtained. Shi et al. [29] utilized the function links between human disease genes and miRNA targets to devise a novel model. A random walk algorithm and global network distance measurement were applied to search for feasible miRNA–disease relationships. Chen et al. [30] utilized a random walk with restart algorithm to construct the RWRMDA model. Because the prediction performance calculated by global network similarity was better than the of the local network [31], RWRMDA employed global network similarity to determine the feasible interactions between microRNAs and diseases. Unfortunately, RWRMDA was inappropriate for the diseases without known associated miRNAs. Liu et al. [32] also implemented a random walk with restart algorithm in their model to make prediction results to a higher degree. They employed the random walk with restart algorithm on a heterogeneous graph established by utilizing disease similarity and miRNA similarity. Luo et al. [33] employed an imbalanced bi-random walk method on a heterogeneous network with information on miRNAs and diseases to identify feasible miRNA–disease interactions. When the random walk algorithm is used for association prediction, the initial state of disease nodes and miRNA nodes in the network is very important. Researchers have proposed many design methods for the initial state of nodes in recent years, but the prediction performance has not been greatly improved.
As artificial intelligence technology has developed, machine-learning-based models have increasingly been employed for the accurate prediction of miRNA–disease relationships. To obtain accurate results in matrix completion for miRNA–disease association prediction, Li et al. [34] avoided using negative samples in MCMDA. To infer unknown miRNA–disease interactions, the probabilistic matrix factorization (PMF) algorithm was applied [35] to make predictions. The PMF algorithm is a machine learning technique commonly employed in recommender systems, and can effectively utilize all available data to recommend miRNAs linked to the disease in question. Ha et al. [36] utilized a matrix completion with network regularization method to recognize potential disease-related miRNAs. They considered an miRNA network as additional implicit feedback, and made predictions for disease associations with a given miRNA relying on its direct neighbors. Guo et al. [37] introduced MLPMDA—a novel model for predicting miRNA–disease associations using multilayer linear projection. They processed miRNA–disease interaction information by processing the top nearest neighbors of entities, and then used the updated miRNA–disease interactions and disease similarity to construct a heterogeneous matrix. In this heterogeneous matrix, the multilayer projection and layer-stacking strategy were employed to make predictions. However, in order to obtain dependable and steady performance, MLPMDA requires high-quality biological data. Ding et al. [38] presented a novel computational model named VGAMF for predicting miRNA–disease associations. VGAMF first integrated several different types of information about miRNAs and diseases into comprehensive similarity networks of miRNAs and diseases, which were used to extract the nonlinear representations of miRNAs and diseases based on the variational graph autoencoders. Then, VGAMF obtained the linear representations of miRNAs and diseases by implementing non-negative matrix factorization to process the miRNA–disease association matrix. Finally, a fully connected neural network combined linear representations with nonlinear representations to generate the predicted miRNA–disease association scores. Wang et al. [39] presented a novel method called NMCMDA to observe unknown disease-related miRNAs. The encoder and decoder were the two essential components in NMCMDA. The encoder was developed using a graph neural network to extract latent miRNA and disease characteristics from a heterogeneous miRNA–disease network. These latent features were used by the decoder to generate miRNA–disease association scores. For NMCMDA, a variety of encoders and decoders have been proposed. Finally, in NMCMDA, the combination of a relational graph convolutional network encoder and a neural multirelational decoder achieved the best prediction results. In summary, machine-learning-based models can produce more accurate prediction results, but most of them have difficulties in adjusting the optimal parameters and selecting negative samples, which seriously affect the training efficiency of the model.
Despite their outstandingly good performance, the abovementioned prediction models have several limitations, such as inadequate measurement of similarity, excessive noise in experimental data, and inaccurate prediction results. To overcome these limitations, we present a novel model called MDSCMF, which combines matrix decomposition with similarity-constrained matrix factorization to predict unobserved miRNA–disease associations. To construct information-rich miRNA similarity and disease similarity, we applied SKF to integrate various kinds of miRNA similarity data and disease similarity data. In addition, because the unknown miRNA–disease associations were much more numerous than the known associations, an MD algorithm was used to get rid of outliers from the miRNA–disease association matrix. Furthermore, we added L 2 regularization terms and similarity constraint terms to non-negative matrix factorization to form an SCMF algorithm, which was implemented to obtain the final association scores of each miRNA–disease pair. To evaluate the effectiveness of MDSCMF, 5-CV, global LOOCV, and local LOOCV were carried out on the known miRNA–disease association data downloaded from HMDD v3.2 [40]. Furthermore, we performed case studies on colon neoplasms, breast neoplasms, and lung neoplasms for prediction. As a result, 29, 29, and 28 out of the top 30 miRNAs potentially connected to these high-risk human diseases, respectively, were confirmed by miR2Disease [41] and dbDEMC v2.0 [42]. Experimental results showed that MDSCMF was effective for inferring possible relationships between miRNAs and diseases.

2. Results

2.1. Performance Evaluation

In this section, based on the verified associations between miRNAs and diseases in the HMDD v3.2 database, 5-CV, global LOOCV, and local LOOCV were implemented to evaluate the prediction performance of MDSCMF.
In the framework of 5-CV, we compared MDSCMF with other previous computational methods, including GCAEMDA [43], MSCHLMDA [44], NIMCGCN [45], and HFHLMDA [46]. The full set of verified miRNA–disease associations were divided into five parts in a random manner, where the test set was held by each part in turn, while the training set consisted of the other four parts. The full set of unknown miRNA–disease associations were considered as candidate samples. We applied our method to determine the ranking of the test set relative to candidate samples. Furthermore, for the purpose of reducing potential deviations resulting in random sample segmentations, we applied 100 repeated segmentations to verify the miRNA–disease associations. When the ranking of all test samples was higher than a certain threshold, MDSCMF was regarded as a valid method. Then we could utilize the receiver operating characteristic (ROC) curve that was obtained by plotting the true positive rate (TPR) against the false positive rate (FPR) to effectively evaluate the performance of MDSCMF. We could calculate the areas under the ROC curve (AUCs) of these methods, whose values were between 0 and 1. Figure 1 indicates that MDSCMF, GCAEMDA, MSCHLMDA, NIMCGCN, and HFHLMDA had AUC values of 0.9488, 0.9415, 0.9324, 0.9378, and 0.9301, respectively. The AUC value of MDSCMF was clearly higher than that of the other methods.
In the framework of global LOOCV, MDSCMF was also compared with GCAEMDA, MSCHLMDA, NIMCGCN, and HFHLMDA. The test set was held by each verified miRNA–disease association in turn, while the training set was composed of the other verified associations. The full set of unknown miRNA–disease associations were considered as candidate samples. In addition, we applied MDSCMF to obtain all predicted association scores so that the ranking of the test set relative to the candidate samples could be determined. Similar to 5-CV, we also calculated the AUCs of these methods so as to effectively evaluate their performance. From Figure 2, we can see that MDSCMF, GCAEMDA, MSCHLMDA, NIMCGCN, and HFHLMDA had AUC values of 0.9540, 0.9505, 0.9378, 0.9410, and 0.9321, respectively. Hence, the AUC value of MDSCMF was also higher than that of the other methods.
In the framework of local LOOCV, we also compared MDSCMF with other previous models (i.e., RFMDA [47], BNPMDA [48], ABMDA [49] and VGAMF [38]) to objectively evaluate its performance. In this way, we could determine the ability of MDSCMF to predict the associations between miRNAs and diseases without any verified related miRNAs. For random diseases in the HMDD v3.2 database, the confirmed associations between each disease and all miRNAs were considered as the test set, and remaining associations were regarded as the training set. Similar to the previous two cross-validation methods, the AUC value in local LOOCV still served as the evaluation criterion to reflect the ability of these models. The specific results are shown in Figure 3, which shows that the prediction performance of MDSCMF was better than that of the other models.

2.2. Parameter Analysis

In this section, the parameters ϑ and σ were quantitatively analyzed to research their effects on the prediction performance. ϑ and σ were set as the regularization parameters, which were applied to control the overfitting degree and the smoothness of similarity consistency, respectively. We utilized all combinations of two values ϑ { 2 3 , 2 2 , , 2 3 } and σ { 2 3 , 2 2 , , 2 3 } to conduct MDSCMF. The AUC values of 5-CV were applied to evaluate the performance of the model under different combinations of parameters. After various tests were conducted, we concluded that the model obtained the best performance when ϑ = 2 2 and σ = 2 0 , as shown in Figure 4.

2.3. Effects of Matrix Decomposition Analysis

In this section, we evaluated the effect of the pre-processing MD step for known miRNA–disease association matrix A on the model’s performance. The AUC values of 5-CV were considered as indicators, and the corresponding ROC curves are shown in Figure 5. In MDSCMF, MD considers the sparsity of the miRNA–disease association matrix, thereby improving the prediction ability of the model. Conversely, MDSCMF without MD disregards the sparsity of the original association matrix; thus, the noise data in the matrix may reduce the accuracy of the prediction. As shown in Figure 5, the AUC value of MDSCMF under the 5-CV framework was 0.9488. In contrast, the AUC value of MDSCMF without MD under the 5-CV framework was 0.9291. The results of the comparison distinctly show that MDSCMF with MD has a higher AUC value compared to that without MD.

2.4. Case Studies

For the purpose of demonstrating the effectiveness and accuracy of MDSCMF, we applied an evaluation experiment in this study. We implemented several types of human diseases—i.e., colon neoplasms, breast neoplasms, and lung neoplasms—as case studies to validate the performance of our method. Colon neoplasms are malignancies in the field of medicine that have been confirmed to be associated with several miRNAs [50,51]. Breast neoplasms, which have been observed to be associated with several miRNAs in clinical experiments, have a high incidence rate among women [52]. Lung neoplasms are among the most dangerous malignancies, with the fastest increases in morbidity and mortality [53]. A growing body of evidence indicates that these diseases have close relationships with several miRNAs. The miRNAs associated with these diseases were ranked in line with the prediction scores. Moreover, we utilized two databases—miR2Disease [41] and dbDEMC v2.0 [42]—to check miRNAs that had been ranked.
As a result, 29, 29, and 28 of the top 30 miRNAs inferred by our model were individually confirmed to be associated with colon neoplasms, breast neoplasms, and lung neoplasms, respectively, according to the miR2Disease [41] and dbDEMC v2.0 [42] databases. Table 1, Table 2 and Table 3 show the corresponding prediction results.

3. Materials and Methods

In this paper, we utilized the biological information of miRNAs and diseases to propose a novel method called MDSCMF, which fully extends the advantages of matrix decomposition and similarity-constrained matrix factorization to predict possible miRNA–disease associations. The flowchart of MDSCMF is clearly shown in Figure 6.

3.1. Human miRNA–Disease Associations

In this study, we took advantage of miRNA–disease association data from the HMDD v3.2 database [40], which contained 12,446 verified associations between 853 miRNAs and 591 diseases. To make calculation more convenient, we constructed an adjacency matrix A R n m × n d to indicate the miRNA–disease relationships. We set n d and n m to stand for the numbers of diseases and miRNAs, respectively. Specifically, the element A ( i , j ) is equal to 1 when miRNA m i is proved to be connected with disease d j , and otherwise it is equal to 0. Therefore, the matrix A contains 12,446 entries that are equal to 1.

3.2. MiRNA Functional Similarity

The miRNAs with similar functions have a high probability of being related to diseases that are similar, and vice versa [20]. Therefore, we downloaded the miRNA functional similarity data from http://www.cuilab.cn/files/images/cuilab/misim.zip, accessed on 1 June 2022. For ease of calculation, we constructed the matrix S M 1 to store the data. The element S M 1 ( m i , m j ) represents the value of similarity between miRNA m i and miRNA m j .

3.3. Disease Semantic Similarity

The directed acyclic graph (DAG) based on the MeSH descriptor [54] can be utilized to describe diseases. D A G ( D ) = ( D , T ( D ) , E ( D ) ) represents the DAG of disease D. T ( D ) denotes the nodes in the DAG that include D itself and its ancestor nodes. E ( D ) denotes the edges in the DAG that connect child nodes with parent nodes directly. The formula to calculate the semantic score of disease D is defined as follows:
D V 1 ( D ) = d ϵ T ( D ) D D 1 ( d ) ,
where the formula to calculate the contribution value D D 1 ( d ) of disease d is as follows:
D D 1 ( d ) = {                                   1                               i f   d = D max { Δ D D 1 ( d ) | d ϵ c h i l d r e n   o f   d }   i f   d D ,
where Δ is the semantic contribution factor, which was equal to 0.5 in our paper, based on previous literature [55].
The formula to obtain the semantic similarity score between disease d i and disease d j is defined as follows:
S S 1 ( d i , d j ) = t T ( d i ) T ( d j ) ( D d i 1 ( t ) + D d j 1 ( t ) ) D V 1 ( d i ) + D V 1 ( d j ) .
Furthermore, for the two diseases of the same layer in a DAG, assuming they have different occurrences in DAGs, it does not make sense to define the semantic contributions of the two diseases for this DAG to be consistent. Objectively speaking, the semantic contribution of high-incidence diseases should be less than that of low-incidence diseases. Consequently, to further optimize the similarity information between diseases, another strategy was introduced to calculate disease semantic similarity following this method [56]. Specifically, the formulae to calculate the semantic score of disease D and the contribution values of disease d are as follows:
D V 2 ( D ) = d ϵ T ( D ) D D 2 ( d ) ,
D D 2 ( d ) = log ( t h e   n u m b e r   o f   D A G s   i n c l u d i n g   d   t h e   n u m b e r   o f   d i s e a s e s ) .
Then, the formula to obtain the semantic similarity score between d i and disease d j is as follows:
S S 2 ( d i , d j ) = t T ( d i ) T ( d j ) ( D d i 2 ( t ) + D d j 2 ( t ) ) D V 2 ( d i ) + D V 2 ( d j ) .
For the purpose of making the results more accurate, we set two kinds of semantic similarity that were equally important. Therefore, if disease d i and d j had semantic similarity, we calculated the average S D 1 ( d i , d j ) of S S 1 ( d i , d j ) and S S 2 ( d i , d j ) by the following formula:
S D 1 ( d i , d j ) = S S 1 ( d i , d j ) + S S 2 ( d i , d j ) 2 .

3.4. Gaussian Interaction Profile Kernel Similarity

The miRNAs with similar functions have a high probability of being related to similar diseases, and vice versa [20]. Therefore, the Gaussian interaction profile kernel similarity was applied to determine the miRNA similarity and disease similarity [57,58]. We made vector K ( d i ) to represent the interaction profile of disease d i in accordance with whether or not d i had a verified association with each miRNA. Similarly, we made vector K ( m i ) to represent the interaction profile m i in accordance with whether or not m i had a verified association with each disease. The equation to calculate the GIP kernel similarity of diseases is defined as follows:
S D 2 ( d i , d j ) = exp ( ρ d K ( d i ) K ( d j ) 2 ) ,
where ρ d is applied to control the kernel bandwidth. The ρ d   is obtained by normalizing the original bandwidth ρ d by the average number of verified associations with miRNAs per disease, as follows:
ρ d = ρ d / ( 1 n d i = 1 n d K ( d i ) 2 ) .
Similarly, we used the following equations to calculate the GIP kernel similarity of miRNAs:
S M 2 ( m i , m j ) = exp ( ρ m K ( m i ) K ( m j ) 2 ) ,
ρ m = ρ m / ( 1 n m i = 1 n m K ( m i ) 2 ) .

3.5. Integrating Similarity for miRNAs and Diseases

In this section, the similarity kernel fusion [59] was implemented to integrate miRNA functional similarity and GIP kernel similarity into ultimate miRNA similarity. The concrete integration process of miRNA similarity matrices can be divided into the following major steps:
In the first step, two different miRNA similarities are treated as original miRNA similarity kernels, which are defined as S M n , n = 1 ,   2 in the above sections. Each miRNA similarity is normalized by the following equation:
F n ( m i , m j ) = S M n ( m i , m j ) m k M S M n ( m k , m j ) ,
where F n ( m i , m j ) denotes the normalized kernel that satisfies m k M F n ( m i , m j ) = 1 , and M = { m i } i = 1 n m indicates the set of miRNAs.
In the second step, the neighbor-constraint kernel for each miRNA original kernel can be constructed as follows:
S n ( m i , m j ) = { S M n ( m i , m j ) m k N i S M n ( m i , m k )   i f   m j N i                 0                   i f   m j N i ,
where S n ( m i , m j ) denotes a neighbor-constraint kernel that obeys m k M S n ( m i , m j ) = 1 , and N i denotes the collection of all neighbors of miRNA m i , including itself.
In the third step, the normalized kernels and neighbor-constraint kernels are integrated as follows:
F n l + 1 = τ ( S n × t n F t l 2 × S n T ) + ( 1 τ ) t n F t 0 2 ,
where F n l + 1 represents the value of n-th kernel after l + 1 iterations, P t 0 represents the initial value of F t , and the weight parameter τ ( 0 , 1 ) is used to balance the rate. After F n l + 1 , n = 1 ,   2 is obtained, the overall kernel S M * can be calculated by the following formula:
S M * = 1 2 n = 1 2 F n l + 1 .
In the fourth step, a weighted matrix W is applied to further eliminate noises in the overall kernel S M * . The construction process of W is as follows:
W ( m i , m j ) = { 1         i f   m i N j m j N i 0         i f   m i N j m j N i 0.5             o t h e r w i s e             .
In the last step, the ultimate miRNA similarity kernel S M R n m × n m can be calculated by the following formula:
S M = W × S M * .
In the same light, we could obtain the ultimate disease similarity kernel as S D R n d × n d .

3.6. Matrix Decomposition

From the published literature [60], we found that the data used in experiments were far from perfect. Several real data of miRNA–disease associations were redundant and/or missing. Therefore, we decomposed the adjacency matrix A into two sections: The linear combination of the adjacency matrix A and low-rank matrix Y was the first section. The second section was the sparse matrix X, which included a large number of zero values. Clearly, the data of the sparse matrix X can be regarded as outliers. The matrix decomposition method was applied to acquire the lowest-rank matrix, which was employed to reconstruct a novel adjacency matrix. The formula to decompose the adjacency matrix A is defined as follows:
A = A Y + X .
For the purpose of making the Y become low-rank, we could enforce nuclear norm on Y . In addition, the L 2 , 1 norm was enforced on the X so that X became sparse. The specific process can be represented by the following formula:
min Y , X Y * + φ X 2 , 1 s . t . A = A Y + X ,
where Y * = i β i ( i . e . ,   β i   i s   t h e   s i g u l a r   v a l u e s   o f   Y ) represents the nuclear norm of Y , X 2 , 1 = j = 1 n i = 1 n ( X i j ) 2 represents the L 2 , 1 norm of X , and the positive free parameter φ is applied to balance the weights of Y * and X 2 , 1 . Furthermore, minimizing the nuclear norm of Y and the L 2 , 1 norm of X contributes to convenient calculation.
If the matrix A combined with Y is treated as an identity matrix, the algorithm will become the robust PCA. Therefore, Equation (19) can be treated as a robust PCA generalization [61] and changed into a comparable problem, as follows:
min Y , X , J J * + φ X 2 , 1 s . t . A = A Y + X ,   Y = J .
Equation (20) is a constraint and convex optimization problem. Therefore, both the first-order information and the special properties of these convex optimization problems can be employed to solve the issue of scalability. The inexact augmented Lagrange multipliers (IALM) algorithm [62] can be utilized to convert Equation (20) to an unconstraint problem. Then, the augmented Lagrange function is adopted to minimize this problem, as follows:
L = J * + φ X 2 , 1 + t r ( F 1 T ( A A Y X ) ) + t r ( F 2 T ( Y J ) ) + α 2 ( A A Y X F 2 + Y J F 2 ) ,
where the penalty parameter α 0 . After minimization with respect to J , Y , and X , the above problem can be settled effectively. In addition, Equation (21) can be solved when the other variables are fixed and the Lagrange multipliers F 1 and F 2 are updated. The specific steps for solving Equation (21) are displayed in Algorithm 1.
Algorithm 1: Solving Equation (21) by IALM
Input: Given an incomplete matrix A and parameters φ { 0 , 1 }
Output: Y * and X *
Initialize: Y = 0 ,   X = 0 ,   F 1 = 0 ,   F 2 = 0 ,   α = 10 4 , m a x α = 10 10 , γ = 1.1 , ε = 10 8
while not converged do
1: Fix the other and update J by J = a r g   m i n 1 α J * + 1 2 J ( Y + F 2 α ) F 2
2: Fix the other and update Y by Y = ( I + A T A ) 1 ( A T A A T X + J + ( A T F 1 F 2 ) / α )
3: Fix the other and update X by X = arg m i n φ α | | X | | 2 , 1 + 1 2 X ( A Y + F 1 / α ) F 2
4: Update the multiplier F 1 = F 1 + α ( A A Y X ) ;   F 2 = F 2 + α ( Y J )
5: Update parameter α by α = min ( γ α , m a x α )
6: Check the convergence condition | | A A Y X | | < ε   a n d   | | Y J | | < ε  
end while
We defined the solution of Equation (21) as Y * and X * . The A ( i , j ) was used to represent the association between miRNA m i and disease d j , so Y * R n d × n d could be applied to represent the similarity between diseases. When Y * was obtained, the adjacency matrix A * denoted new associations between miRNAs and diseases that could be calculated by the following equation:
A * = A Y * .

3.7. Similarity-Constrained Matrix Factorization

In this section, the L 2 regularization terms and similarity constraint terms were added to a traditional non-negative matrix factorization algorithm to form similarity-constrained matrix factorization, which was applied to observe more unknown miRNA–disease interactions. The matrix A * R n m × n d can be factorized into U R n m × γ and V R n d × γ , where γ represents the dimensions of miRNA features and disease features. Concretely, the miRNA–disease association can be regarded as the inner product between the miRNA feature vector and the disease feature vector: a * i j u i v j T , where a * i j indicates the ( i , j ) th element of matrix A * , while u i and v j indicate the i th row of U and the j th row of V , respectively. The corresponding objective function is defined as follows:
m i n 1 2 i j ( a * i j u i v j T ) 2 .
In what follows, the L 2 regularization terms of u i and v j are added to above function for preventing overfitting in the model:
m i n 1 2 i j ( a * i j u i v j T ) 2 + ϑ 2 i u i 2 + ϑ 2 j v j 2 ,
where ϑ denotes the regularization parameter for controlling the balance.
When data points are mapped from high-rank space into low-rank space, the geometric properties of the data points will most likely stay the same [63,64]. Owing to the miRNA similarity S M and disease similarity S D being able to represent the geometric structure of the data points, the similarity constraint terms S U and S V are proposed as follows:
S U = 1 2 i j u i u j 2 S M i j ,
S V = 1 2 i j v i v j 2 S D i j ,
where S M i j represents the similarity between miRNAs m i and m j , while S D i j denotes the similarity between diseases d i and d j . Because the degree of similarity between two random data points is determined by the distance between them, S U will incur a heavy penalty if the distance between m i and m j is close in the miRNA feature space. Thus, we minimized the S U to keep the geometric structure of the miRNA data points, which would give rise to m i and m j being mapped closely in low-dimensional space. The same is true for disease data nodes, so we also minimized the S V . Based on the above analysis, the objective function of SCMF can be defined by adding S U and S V to Equation (24) as follows:
min U , V L = 1 2 i j ( a * i j u i v j T ) 2 + ϑ 2 i u i 2 + ϑ 2 j v j 2 + σ 2 i j u i u j 2 S M i j + σ 2 i j v i v j 2 S D i j ,
where σ denotes the hyperparameter to control the smoothness degree of similarity consistency. Subsequently, an efficacious optimization algorithm is proposed to calculate the objective function of SCMF.
First, the partial derivatives of L with respect to u i and v j can be calculated by the following formulae:
                    u i L = j ( u i v j T a * i j ) v j + ϑ u i + σ ( j ( u i u j ) S M i j j ( u j u i ) S M j i )     = u i ( V T V + ϑ I + σ ( j S M i j + j S M j i ) I ) A * ( i , : ) V σ j ( S M i j + S M j i ) u j ,
                    v j L = j ( v j u i T a * i j ) u i + ϑ v j + σ ( i ( v j v i ) S D j i i ( v i v j ) S D i j )     = v j ( U T U + ϑ I + σ ( i S D i j + i S D j i ) I ) A * ( : , j ) T U σ i ( S D i j + S D j i ) v i ,
where A * ( i , : ) and A * ( : , j ) indicate the   i th row and j th column of matrix A * , respectively.
Next, the calculation of the second derivatives of L with respect to u i and v j is presented as follows:
u i 2 L = V T V + ϑ I + σ ( j S M i j + j S M j i ) I ,
v j 2 L = U T U + ϑ I + σ ( i S D i j + i S D j i ) I .
Then, u i and v j can be iteratively updated according to Newton’s method, as follows:
u i u i u i L ( u i 2 L ) 1 ,
v j v j v j L ( v j 2 L ) 1 .
More specifically, the update of u i and v j can be performed by the below formulas:
u i ( A * ( i , : ) V + σ j ( S M i j + S M j i ) u j ) ( V T V + ϑ I + σ ( j S M i j + j S M j i ) I ) 1 ,
v j ( A * ( : , j ) T U + σ i ( S D i j + S D j i ) v i ) ( U T U + ϑ I + σ ( i S D i j + i S D j i ) I ) 1 .
The update of u i and v j will stop when the convergence condition is satisfied. After that, the predicted association matrix can be calculated by the following formula:
A = U V T .
The value of A i j denotes the predicted association score between miRNA m i and disease d j . The higher the prediction score, the greater the association probability.

4. Discussion

To solve the problems of inadequate measurement of similarity, excessive noise in experimental data, and inaccurate prediction results existing in previous prediction models, we developed a computational model for predicting miRNA–disease associations based on matrix decomposition and similarity-constrained matrix factorization (MDSCMF). Because the miRNA–disease association matrix was a sparse matrix, we applied the MD algorithm to complete it. Our results demonstrated that the MD algorithm could improve the prediction performance to some extent. In addition, we applied SKF to integrate various types of similarities for constructing information-rich miRNA similarity and disease similarity. Furthermore, L 2 regularization terms and similarity constraint terms were added to non-negative matrix factorization to form the SCMF algorithm, which was utilized to generate association scores of each miRNA–disease pair. In the frameworks of 5-CV, global LOOCV, and local LOOCV, the AUCs of MDSCMF achieved 0.9488, 0.9540, and 0.8672, respectively, indicating that the performance of our method had a significant improvement relative to previous methods. Furthermore, the predicted miRNAs related to colon neoplasms, prostate neoplasms, and lung neoplasms were confirmed by the experimental literature, so the prediction results generated by our method were proven to be reliable.
It should be noted that the following factors may contribute to the reliable performance of MDSCMF: First of all, the MD algorithm, which greatly alleviated the influence of the inherent noise existing in the current dataset, was utilized to refine the miRNA–disease association matrix. In addition, when we used SCMF to make predictions, the L 2 regularization terms and similarity constraint terms could avoid overfitting problems and generate robustness of the data richness, respectively.
However, several limitations may influence the performance of MDSCMF. First of all, although the amount of data had increased, we still ought to spare no effort to expand the experimental data. Furthermore, the data we utilized included miRNA function similarity data and disease semantic similarity data, which may contain noise and outliers. Therefore, we should continuously optimize our model to improve its performance in the future.

Author Contributions

Conceptualization, J.N. and L.L.; methodology, J.N.; validation, C.J., Y.W. and C.Z.; investigation, J.N.; data curation, L.L.; writing—original draft preparation, J.N.; writing—review and editing, L.L. and C.Z.; supervision, L.L.; funding acquisition, J.N., Y.W. and C.Z.; All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (grant numbers 61873001, U19A2064, 11701318), the Natural Science Foundation of Shandong Province (grant number ZR2020KC022), and the Open Project of Anhui Provincial Key Laboratory of Multimodal Cognitive Computation, Anhui University (grant number MMC202006).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are present within the manuscript or available by request to corresponding author, Lei Li ([email protected]).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vasques, L.R.; Pereira, L.V.; Izzotti, A.; Schoof, C.R.G.; Ribeiro, A.O. MicroRNAs: Modulators of cell identity, and their applications in tissue engineering. MicroRNA 2014, 3, 45–53. [Google Scholar]
  2. Bartel, D.P. MicroRNAs: Genomics, Biogenesis, Mechanism, and Function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [Green Version]
  3. He, L.; Hannon, G.J. MicroRNAs: Small RNAs with a big role in gene regulation. Nat. Rev. Genet. 2004, 5, 522–531. [Google Scholar] [CrossRef] [PubMed]
  4. Lee, R.C.; Feinbaum, R.L.; Ambros, V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell 1993, 75, 843–854. [Google Scholar] [CrossRef]
  5. Ana, K.; Sam, G.J. miRBase: Annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014, 42, 68–73. [Google Scholar]
  6. Jopling, C.L.; Yi, M.; Lancaster, A.M.; Lemon, S.M.; Sarnow, P. Modulation of hepatitis C virus RNA abundance by a liver-specific MicroRNA. Science 2005, 309, 1577–1581. [Google Scholar] [CrossRef] [Green Version]
  7. Cheng, A.M.; Byrom, M.W.; Shelton, J.; Ford, L.P. Antisense inhibition of human miRNAs and indications for an involvement of miRNA in cell growth and apoptosis. Nucleic Acids Res. 2005, 33, 1290–1297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Karp, X.; Ambros, V. Developmental biology. Encountering MicroRNAs in Cell Fate Signaling. Science 2005, 310, 1288–1289. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Miska, E.A. How microRNAs control cell division, differentiation and death. Curr. Opin. Genet. Dev. 2005, 15, 563–568. [Google Scholar] [CrossRef]
  10. Xu, P.; Guo, M.; Hay, B.A. MicroRNAs and the regulation of cell death. Trends Genet. 2004, 20, 617–624. [Google Scholar] [CrossRef]
  11. Lynam-Lennon, N.; Maher, S.; Reynolds, J.V. The roles of microRNA in cancer and apoptosis. Biol. Rev. Camb. Philos. Soc. 2009, 84, 55–71. [Google Scholar] [CrossRef]
  12. Meola, N.; Gennarino, V.A.; Banfi, V. microRNAs and genetic diseases. PathoGenetics 2009, 2, 7. [Google Scholar] [CrossRef] [Green Version]
  13. Yanaihara, N.; Caplen, N.; Bowman, E.; Seike, M.; Kumamoto, K. Unique microRNA molecular profiles in lung cancer diagnosis and prognosis. Cancer Cell 2006, 9, 189–198. [Google Scholar] [CrossRef] [Green Version]
  14. Yanaihara, N.; Caplen, N.; Bowman, E.; Seike, M.; Kumamoto, K.; Yi, M.; Stephens, R.M.; Okamoto, A.; Yokota, J.; Tanaka, T.; et al. Circulating microRNAs as potential new biomarkers for prostate cancer. Cancer Cell 2013, 108, 1925–1930. [Google Scholar]
  15. Thomson, J.M.; Parker, J.S.; Hammond, S.M. Microarray Analysis of miRNA Gene Expression. Methods Enzymol. 2007, 427, 107–122. [Google Scholar] [CrossRef]
  16. Han, K.; Xuan, P.; Ding, J.; Zhao, Z.J.; Hui, L.; Zhong, Y.L. Prediction of disease-related microRNAs by incorporating functional similarity and common association information. Genet. Mol. Res. 2014, 13, 2009–2019. [Google Scholar] [CrossRef]
  17. You, Z.; Huang, Z.; Zhu, Z.; Yan, G.; Li, Z.; Wen, Z.; Chen, X. PBMDA: A novel and effective path-based computational model for miRNA-disease association prediction. PLoS Comput. Biol. 2017, 13, e1005455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Chen, X.; Gong, Y.; Zhang, D.-H.; You, Z.; Li, Z.-W. DRMDA: Deep representations-based miRNA-disease association prediction. J. Cell. Mol. Med. 2018, 22, 472–485. [Google Scholar] [CrossRef]
  19. Zeng, X.; Xuan, X.; Quan, B. Integrative approaches for predicting microRNA function and prioritizing disease-related mi-croRNA using biological interaction networks. Brief Bioinform. 2016, 17, 193–203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Lu, M.; Zhang, Q.; Min, D.; Jing, M.; Guo, Y.; Guo, W.; Cui, Q. An Analysis of Human MicroRNA and Disease Associations. PLoS ONE 2008, 3, e3420. [Google Scholar] [CrossRef] [Green Version]
  21. Jiang, Q.; Hao, Y.; Wang, G.; Juan, L.; Wang, Y. Prioritization of disease microRNAs through a human phenome-microRNAome network. BMC Syst. Biol. 2010, 4, S2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Li, X.; Wang, Q.; Zheng, Y.; Lv, S.; Ning, S.; Sun, J.; Huang, T.; Zheng, Q.; Ren, H.; Xu, J.; et al. Prioritizing human cancer microRNAs based on genes’ functional consistency between microRNA and cancer. Nucleic Acids Res. 2011, 39, e153. [Google Scholar] [CrossRef] [PubMed]
  23. Mørk, S.; Sune, P.F.; Albert, P.C.; Jan, G.; Jensen, L.J. Protein-driven inference of miRNA–disease associations. Bioinformatics 2014, 30, 392–397. [Google Scholar] [CrossRef] [Green Version]
  24. Chen, X.; Yan, C.; Zhang, X.; You, Z.; Deng, L.; Liu, Y.; Zhang, Y.; Dai, Q. WBSMDA: Within and Between Score for MiRNA-Disease Association prediction. Sci. Rep. 2016, 6, 21106. [Google Scholar] [CrossRef] [PubMed]
  25. Yu, H.; Chen, X.; Lu, L. Large-scale prediction of microRNA-disease associations by combinatorial prioritization algorithm. Sci. Rep. 2017, 7, 43792. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Chen, X.; Cheng, J.; Yin, J. Predicting microRNA-disease associations using bipartite local models and hubness-aware regression. RNA Biol. 2018, 15, 1192–1205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Ha, J.; Park, C. MLMD: Metric Learning for predicting MiRNA-Disease associations. IEEE Access. 2021, 9, 78847–78858. [Google Scholar] [CrossRef]
  28. Li, L.; Gao, Z.; Wang, Y.; Zhang, M.; Ni, J.; Zheng, C.; Su, Y. SCMFMDA: Predicting microRNA-disease associations based on similarity constrained matrix factorization. PLoS Comput. Biol. 2021, 17, e1009165. [Google Scholar] [CrossRef]
  29. 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 miRNA-disease associations through the functional link between miRNA targets and disease genes. BMC Syst. Biol. 2013, 7, 101. [Google Scholar] [CrossRef] [Green Version]
  30. Chen, X.; Liu, M.; Yan, G. RWRMDA: Predicting novel human microRNA-disease associations. Mol. Biosyst. 2012, 8, 2792–2798. [Google Scholar] [CrossRef]
  31. Köhler, S.; Bauer, S.; Horn, D.; Robinson, P.N. Walking the Interactome for Prioritization of Candidate Disease Genes. Am. J. Hum. Genet. 2008, 82, 949–958. [Google Scholar] [CrossRef] [Green Version]
  32. Liu, Y.; Zeng, X.; He, Z.; Zou, Q. Inferring microRNA-disease associations by random walk on a heterogeneous network with multiple data sources. ACM Trans. Comput. Biol. Bioinform. 2017, 14, 905–915. [Google Scholar] [CrossRef] [PubMed]
  33. Luo, J.; Xiao, Q. A novel approach for predicting microRNA-disease associations by unbalanced bi-random walk on heterogeneous network. J. Biomed. Inform. 2017, 66, 194–203. [Google Scholar] [CrossRef]
  34. Li, J.; Rong, Z.; Chen, X.; Yan, G.; You, Z. MCMDA: Matrix completion for MiRNA-disease association prediction. Oncotarget 2017, 8, 21187–21199. [Google Scholar] [CrossRef] [Green Version]
  35. Xu, J.; Cai, L.; Liao, B.; Zhu, W.; Yang, J. Identifying Potential miRNAs-Disease Associations with Probability Matrix Factorization. Front. Genet. 2019, 10, 1234. [Google Scholar] [CrossRef] [PubMed]
  36. Hua, J.; Park, C.; Park, C.; Park, S. Improved Prediction of miRNA-Disease Associations Based on Matrix Completion with Network Regularization. Cells 2020, 9, 881. [Google Scholar] [CrossRef] [Green Version]
  37. Guo, L.; Shi, K.; Lin, W. MLPMDA: Multi-layer linear projection for predicting miRNA-disease association. Knowl. Based Syst. 2021, 214, 106718. [Google Scholar] [CrossRef]
  38. Ding, Y.; Lei, X.; Liao, B.; Wu, F. Predicting miRNA-Disease Associations Based on Multi-View Variational Graph Auto-Encoder with Matrix Factorization. IEEE J. Biomed. Health Inform. 2021, 26, 446–457. [Google Scholar] [CrossRef] [PubMed]
  39. Wang, J.; Li, J.; Yue, K.; Wang, L.; Ma, Y.; Li, Q. NMCMDA: Neural multicategory MiRNA-disease association prediction. Brief Bioinform. 2021, 22, bbab074. [Google Scholar] [CrossRef] [PubMed]
  40. Huang, Z.; Shi, J.; Gao, Y.; Cui, C.; Zhang, S.; Li, J.; Zhou, Y.; Cui, Q. HMDD v3. 0: A database for experimentally supported human microRNA-disease associations. Nucleic Acids Res. 2019, 47, D1013–D1017. [Google Scholar] [CrossRef] [Green Version]
  41. Jiang, Q.; Wang, Y.; Hao, Y.; Juan, L.; Teng, M.; Zhang, X.; Li, M.; Wang, G.; Liu, Y. miR2Disease: A manually curated database for microRNA deregulation in human disease. Nucleic Acids Res. 2009, 37, D98–D104. [Google Scholar] [CrossRef] [Green Version]
  42. Yang, Z.; Wu, L.; Wang, A.; Tang, W.; Zhao, Y.; Zhao, H.; Teschendorff, A.E. dbDEMC 2.0: Updated database of differ-entially expressed miRNAs in human cancers. Nucleic Acids Res. 2017, 45, D812–D818. [Google Scholar] [CrossRef] [PubMed]
  43. Li, L.; Wang, Y.; Ji, C.; Zheng, C.; Ni, J.; Su, Y. GCAEMDA: Predicting miRNA-disease associations via graph convolutional. PLoS Comput. Biol. 2021, 17, e1009655. [Google Scholar] [CrossRef]
  44. Wu, Q.; Wang, Y.; Gao, Z.; Ni, J.; Zheng, C. MSCHLMDA: Multi-Similarity Based Combinative Hypergraph Learning for Predicting MiRNA-Disease Association. Front. Genet. 2020, 11, 354. [Google Scholar] [CrossRef]
  45. Li, J.; Zhang, S.; Liu, T.; Ning, C.; Zhang, Z.; Zhou, W. Neural Inductive Matrix Completion with Graph Convolutional Networks for miRNA-disease Association Prediction. Bioinformatics 2020, 36, 2538–2546. [Google Scholar] [CrossRef]
  46. Wang, Y.; Wu, Q.; Gao, Z.; Ni, J.; Zheng, C. MiRNA-disease association prediction via hypergraph learning based on high-dimensionality features. BMC Med. Inform. Decis. Mak. 2021, 21, 133. [Google Scholar] [CrossRef] [PubMed]
  47. Chen, X.; Wang, C.; Yin, J.; You, Z. Novel Human miRNA-Disease Association Inference Based on Random Forest. Mol. Ther. -Nucleic Acids 2018, 13, 568–579. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Chen, X.; Xie, D.; Wang, L.; Zhao, Q.; You, Z.; Liu, H. BNPMDA: Bipartite Network Projection for MiRNA–Disease Association prediction. Bioinformatics 2018, 34, 3178–3186. [Google Scholar] [CrossRef] [Green Version]
  49. Zhao, Y.; Chen, X.; Yin, J. Adaptive boosting-based computational potential miRNA-disease associations. Bioinformatics 2020, 36, 330. [Google Scholar] [CrossRef]
  50. DeSantis, C.E.; Miller, K.D.; Sauer, A.G.; Jemal, A.; Siegel, R.L. Cancer statistics for African Americans, 2019. CA Cancer J. Clin. 2019, 69, S7. [Google Scholar] [CrossRef] [Green Version]
  51. Thackeray, E.W.; Charatcharoenwitthaya, P.; Elfaki, D.; Sinakos, E.; Lindor, K.D. Colon Neoplasms Develop Early in the Course of Inflammatory Bowel Disease and Primary Sclerosing Cholangitis. Clin. Gastroenterol. Hepatol. 2011, 9, 52–56. [Google Scholar] [CrossRef] [PubMed]
  52. Fu, S.W.; Chen, L.; Man, Y.-G. miRNA Biomarkers in Breast Cancer Detection and Management. J. Cancer 2011, 2, 116–122. [Google Scholar] [CrossRef] [PubMed]
  53. Wang, X.; Yin, H.; Zhang, L.; Zhang, D.; Zhu, J. The construction and analysis of the aberrant lncRNA-miRNA-mRNA network in non-small cell lung cancer. J. Thorac. Dis. 2019, 11, 1772–1778. [Google Scholar] [CrossRef]
  54. Lipscomb, C.E. Medical Subject Headings (MeSH). Bull. Med. Libr. Assoc. 2000, 88, 265–266. [Google Scholar]
  55. Chen, X.; Yin, J.; Qu, J.; Huang, L.; Wang, E. MDHGI: Matrix Decomposition and Heterogeneous Graph Inference for miRNA-disease association prediction. PLoS Comput. Biol. 2018, 14, e1006418. [Google Scholar] [CrossRef]
  56. Xuan, P.; Han, K.; Guo, M.; Guo, Y.; Huang, Y. Prediction of microRNAs Associated with Human Diseases Based on Weighted k Most Similar Neighbors. PLoS ONE 2013, 8, e70204. [Google Scholar] [CrossRef]
  57. Van, L.T.; Nabuurs, S.B.; Marchiori, E. Gaussian inter-action profile kernels for predicting drug-target interaction. Bioinformatics 2011, 27, 3036–3043. [Google Scholar]
  58. Chen, X.; Huang, Y.-A.; You, Z.-H.; Yan, G.-Y.; Wang, X. A novel approach based on KATZ measure to predict associations of human microbiota with non-infectious diseases. Bioinformatics 2017, 33, 733–739. [Google Scholar] [CrossRef] [PubMed]
  59. Wang, B.; Mezlini, A.M.; Demir, F.; Fiume, M.; Tu, Z.W.; Brudno, M.; Haibe-Kains, B.; Goldenberg, A. Similarity network fusion for aggregating data types on a genomic scale. Nat. Methods 2014, 11, 333–337. [Google Scholar] [CrossRef]
  60. Yu, S.; Liang, C.; Xiao, Q.; Li, G.; Ding, P.; Luo, J. MCLPMDA: A novel method for miRNA-disease association prediction based on matrix completion and label propagation. J. Cell. Mol. Med. 2019, 23, 1427–1438. [Google Scholar] [CrossRef] [PubMed]
  61. Chen, J.; Yang, J. Robust Subspace Segmentation Via Low-Rank Representation. IEEE Trans. Cybern. 2014, 44, 1432–1445. [Google Scholar] [CrossRef]
  62. Meng, F.; Yang, X.; Zhou, C. The Augmented Lagrange Multipliers Method for Matrix Completion from Corrupted Samplings with Application to Mixed Gaussian-Impulse Noise Removal. PLoS ONE 2014, 9, e108125. [Google Scholar] [CrossRef] [PubMed]
  63. Zhang, W.; Liu, X.R.; Chen, Y.L.; Wu, W.J.; Wang, W.; Li, X.H. Feature-derived graph regularized matrix factorization for predicting drug side effects-Science Direct. Neurocomputing 2018, 287, 154–162. [Google Scholar] [CrossRef]
  64. Rana, B.; Juneja, A.; Saxena, M.; Gudwani, S.; Kumaran, S.S.; Behari, M.; Agrawal, R.K. Graph-Theory-based Spectral Feature Selection for Computer Aided Diagnosis of Parkinson’s Disease Using T1-weighted MRI. Int. J. Imaging Syst. Technol. 2015, 25, 245–255. [Google Scholar] [CrossRef]
Figure 1. AUC of 5-CV compared with those of GCAEMDA, MSCHLMDA, NIMCGCN, and HFHLMDA.
Figure 1. AUC of 5-CV compared with those of GCAEMDA, MSCHLMDA, NIMCGCN, and HFHLMDA.
Genes 13 01021 g001
Figure 2. AUC of global LOOCV compared with those of GCAEMDA, MSCHLMDA, NIMCGCN, and HFHLMDA.
Figure 2. AUC of global LOOCV compared with those of GCAEMDA, MSCHLMDA, NIMCGCN, and HFHLMDA.
Genes 13 01021 g002
Figure 3. Comparisons between MDSCMF and other computational models by local LOOCV.
Figure 3. Comparisons between MDSCMF and other computational models by local LOOCV.
Genes 13 01021 g003
Figure 4. AUCs at different values of ϑ and σ .
Figure 4. AUCs at different values of ϑ and σ .
Genes 13 01021 g004
Figure 5. The ROC curves of MDSCMF and MDSCMF without MD.
Figure 5. The ROC curves of MDSCMF and MDSCMF without MD.
Genes 13 01021 g005
Figure 6. Flowchart of MDSCMF.
Figure 6. Flowchart of MDSCMF.
Genes 13 01021 g006
Table 1. The top 30 potential miRNAs associated with colon neoplasms.
Table 1. The top 30 potential miRNAs associated with colon neoplasms.
miRNAEvidencemiRNAEvidence
hsa-mir-630dhsa-mir-29bm; d
hsa-mir-20am; dhsa-mir-141m; d
hsa-mir-143m; dhsa-mir-132m; d
hsa-mir-584dhsa-mir-19bm; d
hsa-mir-506dhsa-mir-29am; d
hsa-mir-552dhsa-mir-223d
hsa-mir-128unconfirmedhsa-let-125bd
hsa-mir-7im; dhsa-mir-622d
hsa-mir-127m; dhsa-mir-18ad
hsa-mir-1290dhsa-mir-143d
hsa-mir-493dhsa-mir-125am; d
hsa-mir-498dhsa-mir-21m; d
hsa-mir-107m; dhsa-mir-137m; d
hsa-mir-191m; dhsa-mir-424d
hsa-mir-32m; dhsa-mir-200bd
m: miR2Disease database; d: dbDEMC v2.0 database.
Table 2. The top 30 potential miRNAs associated with breast neoplasms.
Table 2. The top 30 potential miRNAs associated with breast neoplasms.
miRNAEvidencemiRNAEvidence
hsa-mir-99am; dhsa-mir-663m
hsa-mir-542dhsa-mir-520hd
hsa-mir-96dhsa-mir-519dd
hsa-mir-98m; dhsa-mir-186d
hsa-mir-185dhsa-mir-381d
hsa-mir-130adhsa-mir-32d
hsa-mir-708dhsa-mir-590unconfirmed
hsa-mir-150dhsa-mir-330d
hsa-mir-192dhsa-mir-433d
hsa-mir-196bdhsa-mir-942d
hsa-mir-888dhsa-mir-661m; d
hsa-mir-9m; dhsa-mir-337d
hsa-mir-130bdhsa-mir-494d
hsa-mir-592dhsa-mir-212d
hsa-mir-99bdhsa-mir-618d
m: miR2Disease database; d: dbDEMC v2.0 database.
Table 3. The top 30 potential miRNAs associated with lung neoplasms.
Table 3. The top 30 potential miRNAs associated with lung neoplasms.
miRNAEvidencemiRNAEvidence
hsa-mir-96dhsa-mir-937unconfirmed
hsa-mir-145m; dhsa-mir-30em
hsa-mir-99am; dhsa-mir-151d
hsa-mir-9m; dhsa-mir-614d
hsa-mir-185dhsa-mir-1323d
hsa-mir-130adhsa-mir-32d
hsa-mir-7m; dhsa-mir-1298d
hsa-mir-150m; dhsa-mir-330d
hsa-mir-192m; dhsa-mir-433d
hsa-mir-769unconfirmedhsa-mir-522d
hsa-mir-939dhsa-mir-449ad
hsa-mir-98m; dhsa-mir-143m; d
hsa-mir-130bm; dhsa-mir-564d
hsa-mir-638dhsa-mir-212m; d
hsa-mir-99bdhsa-mir-615unconfirmed
m: miR2Disease; d: dbDEMC v2.0 database.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ni, J.; Li, L.; Wang, Y.; Ji, C.; Zheng, C. MDSCMF: Matrix Decomposition and Similarity-Constrained Matrix Factorization for miRNA–Disease Association Prediction. Genes 2022, 13, 1021. https://doi.org/10.3390/genes13061021

AMA Style

Ni J, Li L, Wang Y, Ji C, Zheng C. MDSCMF: Matrix Decomposition and Similarity-Constrained Matrix Factorization for miRNA–Disease Association Prediction. Genes. 2022; 13(6):1021. https://doi.org/10.3390/genes13061021

Chicago/Turabian Style

Ni, Jiancheng, Lei Li, Yutian Wang, Cunmei Ji, and Chunhou Zheng. 2022. "MDSCMF: Matrix Decomposition and Similarity-Constrained Matrix Factorization for miRNA–Disease Association Prediction" Genes 13, no. 6: 1021. https://doi.org/10.3390/genes13061021

APA Style

Ni, J., Li, L., Wang, Y., Ji, C., & Zheng, C. (2022). MDSCMF: Matrix Decomposition and Similarity-Constrained Matrix Factorization for miRNA–Disease Association Prediction. Genes, 13(6), 1021. https://doi.org/10.3390/genes13061021

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop