Next Article in Journal
ECMpy, a Simplified Workflow for Constructing Enzymatic Constrained Metabolic Network Model
Previous Article in Journal
Cancer-Related Intracellular Signalling Pathways Activated by DOXorubicin/Cyclodextrin-Graphene-Based Nanomaterials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting miRNA-Disease Association Based on Neural Inductive Matrix Completion with Graph Autoencoders and Self-Attention Mechanism

1
College of Computer Science, Nankai University, Tianjin 300350, China
2
College of Artificial Intelligence, Nankai University, Tianjin 300350, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Biomolecules 2022, 12(1), 64; https://doi.org/10.3390/biom12010064
Submission received: 14 December 2021 / Revised: 29 December 2021 / Accepted: 31 December 2021 / Published: 2 January 2022
(This article belongs to the Section Bioinformatics and Systems Biology)

Abstract

:
Many studies have clarified that microRNAs (miRNAs) are associated with many human diseases. Therefore, it is essential to predict potential miRNA-disease associations for disease pathogenesis and treatment. Numerous machine learning and deep learning approaches have been adopted to this problem. In this paper, we propose a Neural Inductive Matrix completion-based method with Graph Autoencoders (GAE) and Self-Attention mechanism for miRNA-disease associations prediction (NIMGSA). Some of the previous works based on matrix completion ignore the importance of label propagation procedure for inferring miRNA-disease associations, while others cannot integrate matrix completion and label propagation effectively. Varying from previous studies, NIMGSA unifies inductive matrix completion and label propagation via neural network architecture, through the collaborative training of two graph autoencoders. This neural inductive matrix completion-based method is also an implementation of self-attention mechanism for miRNA-disease associations prediction. This end-to-end framework can strengthen the robustness and preciseness of both matrix completion and label propagation. Cross validations indicate that NIMGSA outperforms current miRNA-disease prediction methods. Case studies demonstrate that NIMGSA is competent in detecting potential miRNA-disease associations.

1. Introduction

Micro RNAs (miRNAs) are a kind of small non-coding RNAs (about 22 nucleotides) that can regulate target mRNA expression during the post-transcriptional stage, via binding to the 3′-untranslated region of target mRNAs [1,2,3]. Thus, miRNAs can influence a series of biological processes (e.g., epigenetic regulation, cell differentiation, and basal metabolism), playing key roles in plenty of human diseases. For instance, previous research [4] has confirmed that the expression of miRNA hsa-mir-21 can facilitate the proliferation of several kinds of tumor cells, such as breast neoplasms, pancreatic neoplasms, and glioblastoma neoplasms. MiRNA mir-34a can suppress neuroblastoma via promoting tumor cell apoptosis [5]. Therefore, predicting potential miRNA-disease associations is crucial for disease prevention, diagnosis, and treatment.
Potential miRNA-disease associations can be discovered by experimental approaches or computational approaches. As computational approaches, especially machine learning algorithms, are more cost-effective and time-efficient, many machine learning-based methods were proposed to predict potential miRNA-disease associations [3,6,7]. These methods can be categorized into the following types.
1.
Matrix analysis-based methods. Two commonly-used matrix analysis methods for predicting associations among biological entities are manifold regularization [8] and matrix completion [9], which respectively suggest that the association matrix follows manifold constraint or low-rank constraint. Manifold regularization based methods have been widely used for link prediction among biological entities [10,11,12]. Chen et al. [13] proposed a manifold regularized subspace learning method for detecting miRNA-disease associations. Xiao et al. [14] proposed a graph regularized non-negative matrix factorization method to predict microRNA-disease associations. Matrix completion based methods have been commonly used to infer associations among biological entities [15,16,17]. Chen et al. [18] proposed an inductive matrix completion-based method for inferring miRNA-disease associations. Li et al. [19] proposed a matrix completion algorithm for miRNA-disease associations prediction. Yu et al. [20] proposed a matrix completion algorithm for low-rank subspace learning, while incorporating label propagation for miRNA-disease associations prediction. Chen et al. [21] adopted neighborhood constraint matrix completion to predict disease-related miRNAs.
2.
Graph analysis-based methods. Since the dependency among biological entities can be depicted via graphs, methods based on graph algorithms, such as bipartite graph algorithms, neighborhood sampling, and random walk, have been commonly applied in the field of bioinformatics [22,23,24]. Zeng et al. [25] proposed a structural perturbation method-based model for inferring disease-related miRNAs on bipartite miRNA-disease graph. Chen et al. [26] proposed a bipartite network projection-based method for miRNA-disease associations prediction. Xuan et al. [23] adopted a weighted neighborhood sampling algorithm for predicting potential disease-associated miRNAs. Chen et al. [24] proposed a matrix decomposition and heterogeneous graph inference-based model for miRNA-disease association prediction. Since random walk is an efficient way to learn graph representation via topologial relationships of graphs, Chen et al. [22] and Xuan et al. [27] adopted the random walk algorithm to identify potential miRNA-disease associations.
3.
Heterogeneous features fusion methods. Integrating multi-source features is an efficient technique for predicting associations among biological entities [7,16,28]. Peng et al. [29] integrated multiple networks to identify potential miRNA-disease associations. Liu et al. [30] predicted disease-related miRNAs on a heterogeneous network with multiple features. Xiao et al. [31] proposed an adaptive heterogeneous feature inference model for predicting potential disease-associated miRNAs. Ha et al. [32] designed a metric learning model to fuse heterogeneous features for predicting miRNA-disease associations. Yu et al. [33] proposed a multi-layer heterogeneous network embedding model to predict potential miRNA-disease associations.
4.
Deep learning methods. Neural networks have been widely used for detecting potential associations among biological entities [28,33,34]. Zeng et al. [35] adopted a neural network-based model to identify potential miRNA-disease associations. Chen et al. [36] proposed a deep-belief network for inferring disease-related miRNAs. Ji et al. [37] proposed an autoencoder for detecting miRNA-disease associations. Tang et al. [38] proposed a multi-view multi-channel graph attention networks to identify potential miRNA-disease associations. Graph Neural Networks (GNN) [39] have been proposed in deep learning on graphs. Thus, there are some recent studies for predicting associations among biological entities based on GNNs [40,41,42]. Li et al. [43] implemented an inductive matrix completion algorithm based on Graph Convolutional Networks (GCN) for predicting miRNA-disease associations. Li et al. [44] adopted graph autoencoders to identify potential miRNA-disease associations.
In this paper, we propose an inductive matrix completion-based method to predict miRNA-disease associations. Varying from the previous study [18], the inductive matrix completion algorithm in our method is implemented through neural networks, so it is called neural inductive matrix completion. Li et al. [43] implemented a neural inductive matrix completion algorithm based on GCN [45]. In previous works [18,43,44], matrix completion algorithms or graph neural networks were assigned to compute the representations of miRNAs and diseases, then the prediction scores of miRNA-disease pairs were obtained through the dot product of miRNA representations and disease representations. Hence, these works ignored the label propagation procedure for inferring potential miRNA-disease associations. Yu et al. [20] proposed a matrix completion algorithm for low-rank subspace learning and assigned label propagation for miRNA-disease associations prediction. However, in [20], matrix completion and label propagation are separate procedures, leading to decrease their effectiveness. To address this issue, our method propose an end-to-end framework by graph autoencoders to integrate matrix completion and label propagation. Graph autoencoders on miRNA graph and disease graph are competent to reconstruct score matrix through initial association matrix, which is equivalent to propagating labels on graphs. Meanwhile, graph autoencoders on miRNA graph and disease graph are capable of low-rank representation learning from miRNA space and disease space, respectively. The two graph autoencoders are trained collaboratively via neural inductive matrix completion. Since Geng et al. [46] suggested that the attention mechanism is linked to matrix factorization, we illuminate that our inductive matrix completion-based method is an implementation of self-attention mechanism. In a word, our method implement a Neural Inductive Matrix completion-based method with Graph Autoencoders (GAE) and Self-Attention mechanism for miRNA-disease associations prediction. Our method, named NIMGSA, has the following advantages.
1.
NIMGSA implements inductive matrix completion through graph autoencoders, which not only ensures the low-rank property of representations from both miRNA space and disease space, but also depicts label propagation procedure through the reconstruction of association matrix.
2.
NIMGSA integrates inductive matrix completion and label propagation through an end-to-end deep learning framework, which enhances the robustness and preciseness of both integrated procedures.
3.
NIMGSA implements self-attention mechanism through inductive matrix completion on two graph autoencoders, which provides theoretical analysis and biological application to enhance the performance of attention-based neural networks.
4.
The inductive matrix completion procedure is equivalent to training two Graph Autoencoders (i.e., GAE on miRNA graph and GAE on disease graph) collaboratively, which improves the capability for representation learning of these two GAEs.
Experiments demonstrate that NIMGSA is superior to the current state-of-the-art methods. Ablation studies demonstrate the superiority of our proposed architecture of networks. Case studies on several diseases demonstrate the capability of NIMGSA to detect new miRNA-disease associations. The source code of our model is available at https://github.com/zhanglabNKU/NIMGSA (accessed on 13 December 2021).

2. Materials and Methods

2.1. Problem Formulation

Let m and n denote the number of miRNAs and diseases, respectively, and Y R m × n denotes the miRNA-disease association matrix. Y ( i , j ) = 1 if miRNA i has been known to be associated with disease j, otherwise Y ( i , j ) = 0 . An algorithm predicating miRNA-disease associations requires matrix Y, along with miRNA similarity matrix S m R m × m (see Section 2.2) and disease similarity matrix S d R n × n (see Section 2.3), then ensures an optimal score matrix F R m × n , where F ( i , j ) [ 0 , 1 ] denotes the predicted score of the association between miRNA i and disease j. A higher score stands for a higher probability if miRNA i is associated with disease j.
In this paper, the dataset is retrieved from the HMDD v2.0 database (http://www.cuilab.cn/static/hmdd3/data/hmdd2.zip (accessed on 13 December 2021)) [47], including 5430 miRNA-disease associations among 495 miRNAs and 383 diseases.

2.2. MiRNA Similarity Matrix

Wang et al. [48] proposed a method to infer miRNA functional similarities from miRNA-related diseases. The similarity data can be downloaded from MISIM database (http://www.cuilab.cn/files/images/cuilab/misim.zip (accessed on 13 December 2021)). However, there are some miRNAs that are not included in this database. So we compute Gaussian kernel similarity for those not included miRNAs. The gaussian kernel similarity of miRNA pair ( i , j ) is defined as:
G M S ( i , j ) = exp ( Y ( i , : ) Y ( j , : ) 2 / θ m ) ,
where Y ( i , : ) denotes the i-th row of Y, and
θ m = 1 m i = 1 m Y ( i , : ) 2 ,
denotes the kernel bandwidth. In summary, miRNA similarities can be computed through:
S m ( i , j ) = M S ( i , j ) if ( i , j ) has functional similarity G M S ( i , j ) otherwise ,
where M S ( i , j ) denotes the functional similarity of miRNA pair ( i , j ) from the MISIM database.

2.3. Disease Similarity Matrix

Wang et al. [48] and Xuan et al. [23] proposed two different models to compute disease semantic similarities. The Medical Subject Headings (MeSH) database (https://www.ncbi.nlm.nih.gov/ (accessed on 13 December 2021)) is commonly used for describing relationships among human diseases, and both of the models above are based on MeSH descriptors. At the beginning, a hierarchical Directed Acyclic Graph (DAG) is constructed from MeSH, where each node in this DAG denotes a disease, and each directed edge i j denotes a link from disease i to disease j.
Wang et al. [48] suggested that the semantic contribution of disease pair ( i , j ) is computed through:
S C 1 ( i , j ) = 1 j = i max { δ · S C 1 ( i , t ) | t children of j } j i ,
where δ is a hyperparameter and is set as 0.5 in [48]. Suppose N ( i ) denotes a node set including node i itself and its ancestor nodes in disease DAG, disease semantic similarity of disease pair ( i , j ) is computed as:
D S 1 = t N ( i ) N ( j ) ( S C 1 ( i , t ) + S C 1 ( j , t ) ) t N ( i ) S C 1 ( i , t ) + t N ( j ) S C 1 ( j , t ) .
Xuan et al. [23] suggested that the semantic contribution of disease pair ( i , j ) is computed through:
S C 2 ( i , j ) = log the number of DAGs including j the number of diseases .
Then, disease semantic similarity of disease pair ( i , j ) is computed as:
D S 2 = t N ( i ) N ( j ) ( S C 2 ( i , t ) + S C 2 ( j , t ) ) t N ( i ) S C 2 ( i , t ) + t N ( j ) S C 2 ( j , t ) .
However, there are some diseases that are not included in the MeSH database. So we compute Gaussian kernel similarity for those unincluded diseases. The Gaussian kernel similarity of disease pair ( i , j ) is defined as:
G D S ( i , j ) = exp ( Y ( : , i ) Y ( : , j ) 2 / θ d ) ,
where Y ( : , j ) denotes the j-th column of Y, and:
θ d = 1 n j = 1 n Y ( : , j ) 2 ,
denotes the kernel bandwidth. In summary, disease similarities can be computed through:
S d ( i , j ) = D S ( i , j ) if ( i , j ) has semantic similarity G D S ( i , j ) otherwise ,
where
D S ( i , j ) = D S 1 ( i , j ) + D S 2 ( i , j ) 2 ,
denotes the semantic similarity of disease pair ( i , j ) . D S 1 and D S 2 is from Equations (5) and (7).
The procedure of similarity computation is summarized as Figure 1.

2.4. Related Works

2.4.1. Label Propagation

Previous research [49] demonstrates that label propagation is equivalent to solving manifold regularization problem [8] through fixed-point iteration. Manifold regularization assumpts that samples are distributed on a manifold, samples with higher feature similarities are closer on the manifold, and tend to share the same labels. The manifold of data can be depicted by graph structure constructed through feature matrix, which leads to graph semi-supervised learning. This type of method for biological association prediction first computes adjacency matrix from biological features to construct a graph, then propagate labels from labeled biological entities to unlabeled ones on this graph iteratively.
Suppose L denotes normalized Laplacian matrix of the graph, minimizing trace ( F T L F ) can obtain the label matrix F following manifold assumption. Belkin et al. [8] added this manifold constraint to least square problem, then derived Laplacian Regularized Least Square (LRLS) method:
min F F Y F 2 + η trace ( F T L F ) ,
where · F denotes the Frobenius norm of a matrix, and η is a hyperparameter. Equation (12) is a trade-off between the accuracy based on labeled data, and the smoothness of the manifold. This is classified as manifold regularization [8]. Label propagation follows the framework of manifold regularization as Equation (12). Xia et al. [10] derived that association matrix F follows manifold assumption, and can be obtained via solving Equation (12).
Numerous research [50,51,52] demonstrate that Graph Neural Networks (GNN) is closely linked to label propagation algorithm. The outputs of GNN follow the manifold constraint. Hence, a graph autoencoder with Y as input and F as output can obtain the optimal solution of Equation (12). Simulating the label propagation algorithm through the reconstruction procedure of graph autoencoder, has been validated as an efficient way for biological association prediction in previous research [41,42].

2.4.2. Inductive Matrix Completion

Natarajan et al. [15,18] proposed inductive matrix completion to predict associations among biological entities. The matrix completion problem is to approximate the initial miRNA-disease association matrix Y through a low-rank matrix Z R m × n . If rank ( Z ) r min ( m , n ) , Z can be factorized into matrix M R m × r and D R n × r , i.e., Z = M D T . Inductive matrix completion is to optimize:
min M , D 1 2 Y S m M D T S d T F 2 + λ M F 2 + λ D F 2 ,
where S m R m × m and S d R n × n denote similarity matrices of miRNAs and diseases, respectively. Equation (13) can be solved through a non-negative matrix factorization algorithm [53]. After obtaining optimal M and D, score matrix F is computed through:
F = S m M D T S d T .

2.4.3. Attention Mechanism

Attention mechanism [54] in the deep learning model is similar to that in cognitive science, which first calculates a probability distribution over the elements in the inputs, then takes the attention score based on this probability distribution while generating outputs.
The self-attention mechanism [55,56] is a commonly used implementation of attention mechanism. In self-attention mechanism, the input Y is transformed into three matrices, the Query (Q), Key (K), and Value (V), by three different functions. The weight assigned to each value is calculated as the dot-product of the query with the corresponding key:
A t t e n t i o n ( Q , K , V ) = S o f t m a x Q K T d V ,
where d is the dimension of the vector K, and T is the transpose operation. The Q, K, and V are obtained by three linear transformations with the same input separately:
Q = W Q Y , K = W K Y , V = W V Y ,
where W Q , W K , and W V are trainable parameters.
Usually, d is set less than the dimension of Y. It is obvious that the rank of A t t e n t i o n ( Q , K , V ) is no more than d, i.e., the attention score matrix is low-rank. Geng et al. [46] suggested that attention mechanism is linked to matrix factorization for obtaining low-rank outputs.

2.5. NIMGSA

2.5.1. Graph Autoencoder

Suppose Z m = S m M , Z d = S d D , Equation (13) can be rewritten as:
min M , D 1 2 Y Z m Z d T F 2 + λ M F 2 + λ D F 2 ,
where Z m and Z d denote the low-rank representations of miRNAs and diseases, respectively. Previous studies [57,58] have found that autoencoders are competent to obtain low-rank representations. Therefore, we adopt graph autoencoders [59] to obtain low-rank representations in our model, NIMGSA, that Z m and Z d are learned by 2-layer graph convolution [45] encoders, respectively. The encoder is defined as:
Enc ( A , X ) = tanh ( A · ReLU ( A X W ( 0 ) ) W ( 1 ) ) ,
where A, X, and W denote adjacency matrix, inputs, and weights, respectively.
Z m = Enc ( A m , Y ) .
Z d = Enc ( A d , Y T ) .
A m and A d denote the normalized adjacency matrices of miRNA graph and disease graph, respectively.
A m = D m 1 / 2 S m D m 1 / 2 ,
where D m is the degree matrix of S m . D m is a diagonal matrix that is computed via D m ( i , i ) = j S m ( i , j ) . Similarly,
A d = D d 1 / 2 S d D d 1 / 2 .
Then, Equation (17) can be rewritten as:
min W 1 2 Y Z m Z d T F 2 + λ W F 2 .
The encoder-decoder architecture [54] is widely applied for reconstructing outputs from representations. The decoder is defined as:
Dec ( A , X ) = sigmoid ( A · ReLU ( A X W ( 2 ) ) W ( 3 ) ) .
Score matrices F m R m × n and F d R n × m can be decoded through Z m and Z d , respectively:
F m = Dec ( A m , Z m ) .
F d = Dec ( A d , Z d ) .
Following previous research [41], since both F m R m × n and F d R n × m are low-rank provided by autoencoders, and through the rank-sum inequality that:
rank ( α F m + ( 1 α ) F d T ) rank ( F m ) + rank ( F d T ) ,
the final result:
F = α F m + ( 1 α ) F d T
is low-rank, where α ( 0 , 1 ) depicts a balance between miRNA space and disease space.

2.5.2. Self-Attention

In NIMGSA, Z m R m × d and Z d R n × d are equivalent to the Query Q and Key K in self-attention mechanism, which can be obtained by transformations with the same input Y. F can be regarded as Value V of self-attention mechanism. Similar to the definition of attention Equation (15), the attention score of the association matrix can be defined as:
T = S o f t m a x Z m Z d T d F ,
where d is the dimension of hidden vectors, ⊙ denotes element-wise product. Then, Equation (23) can be rewritten as:
min W 1 2 Y T F 2 + λ W F 2 .
Then, following previous research [41,42], we add reconstruction error:
L r = α Y F m F 2 + ( 1 α ) Y F d T F 2 ,
into Equation (30).
min W 1 2 Y T F 2 + β L r + λ W F 2 .
In NIMGSA, we set β = 1 , λ = 10 7 .
The architecture of the NIMGSA model is illustrated as Figure 2. The procedure of NIMGSA is summarized as Figure 3 and Algorithm 1, where GAEm and GAEd represent GAEs on the miRNA graph and disease graph respectively, and NIMC denotes neural inductive matrix completion.
Algorithm 1 NIMGSA Algorithm
Input: 
initial association matrix Y, miRNA similarity matrix S m , disease similarity matrix S d
Output: 
score matrix F
1:
Compute the adjacency matrix of miRNA graph A m and disease graph A d via Equations (21) and (22) respectively
2:
repeat
3:
   Compute Query: Z m = Enc ( A m , Y )    // Equation (19)
4:
   Compute Key: Z d = Enc ( A d , Y T )    // Equation (20)
5:
    F m = Dec ( A m , Z m )    // Equation (25)
6:
    F d = Dec ( A m , Z d )    // Equation (26)
7:
   Compute Value: F = α F m + ( 1 α ) F d T    // Equation (28)
8:
   Compute attention score as Equation (29)
9:
   Train GAEm and GAEd through optimizing Equation (32)    // Update W (i.e., parameters of neural networks) in Equation (32) by Adam optimizer
10:
until Convergence
11:
return F

3. Results

3.1. Comparison with Other Methods

We compare our proposed method, NIMGSA, with other five state-of-the-art methods:
  • IMCMDA: Chen et al. [18] proposed an inductive matrix completion-based method to predict miRNA-disease associations.
  • SPM: Zeng et al. [25] proposed a structural perturbation method- based approach to predict miRNA-disease associations on bipartite miRNA-disease graph.
  • NIMCGCN: Li et al. [43] implemented inductive matrix completion algorithm through graph convolutional networks for miRNA-disease associations prediction.
  • MCLPLDA: Yu et al. [20] adopted matrix completion algorithm for low-rank subspace learning, while integrating label propagation for miRNA-disease associations prediction.
  • GAEMDA: Li et al. [44] adopted graph autoencoders for miRNA-disease associations prediction.
We adopt PyTorch (https://pytorch.org/ (accessed on 13 December 2021)) to construct NIMGSA, and apply an Adam optimizer [60] to train the model. Then, we set the dropout rate [61] of neural networks at 0.5. Our model is trained on a single NVIDIA GeForce GTX 2070 GPU with 8GB of memory.
We adopt five-fold cross validation to evaluate the performance, and the metrics are listed below:
Sensitivity ( SEN ) = TP TP + FN = TPR = Recall
Specificity ( SPEC ) = TN TN + FP = 1 FPR
Accuracy ( ACC ) = TN + TP TN + TP + FN + FP
Precision ( PRE ) = TP TP + FP
F 1 - Score = 2 × Precision × Recall Precision + Recall
MCC = TP × TN FP × FN ( TP + FN ) × ( TP + FP ) × ( TN + FN ) × ( TN + FP )
where TP denotes true positive, FN denotes false negative, TN denotes true negative, FP denotes false negative, TPR denotes true positive rate, FPR denotes false positive rate, and Mcc denotes Matthews correlation coefficient. The Receiver Operating Characteristic (ROC) curve can be plotted by TPR and FPR, while the Area Under ROC curve (AUROC) and the Area under Precision-Recall curve (AUPR) are important metrics to measure the performance of a binary classification model.
We plot the ROC curves and PR curves on Figure 4. The mean values and standard deviations of AUROC and AUPR are listed on Table 1. The results show that VGAELDA outperforms the other state-of-the-art methods in both AUROC and AUPR. In AUROC, NIMGSA achieves an AUROC of 0.9354, which is 0.2% higher than GAEMDA (0.9332), 0.6% higher than MCLPMDA (0.9292), 0.8% higher than NIMCGCN (0.9279), 4.3% higher than SPM (0.8960), and 11.2% higher than IMCMDA (0.8329). In AUPR, NIMGSA achieves an AUPR of 0.4567, which is 4.1% higher than MCLPMDA (0.4387), 10.2% higher than GAEMDA (0.4142), 15.8% higher than NIMCGCN (0.3943), 63.9% higher than IMCMDA (0.2785), and 85.3% higher than SPM (0.2464).
To further evaluate the performance of NIMGSA, we test our model at a high stringency level of specificity according to Equation (34). We fix specificity at 0.99, and then compute sensitivity, accuracy, precision, F1-score, and Mcc. The results are listed on Table 2, which illustrate that NIMGSA outperforms other methods at all five metrics. Matthews correlation coefficient (Mcc) is a comprehensive metric in binary classification on imbalanced data [41]. NIMGSA achieves an Mcc of 0.4273, which is higher than GAEMDA (0.4213), MCLPMDA (0.4138), NIMCGCN (0.3645), IMCMDA (0.3239), and SPM (0.2048).

3.2. Hyperparameter Tuning

In NIMGSA, hyperparameter α ( 0 , 1 ) depicts a balance between miRNA space and disease space. After evaluating our model at each α { 0.1 , 0.3 , 0.5 , 0.7 , 0.9 } , we find that NIMGSA performs best at α = 0.5 . The results are shown on Table 3.
Besides, we evaluate our model at a different learning rate in { 0.001 , 0.01 , 0.05 , 0.1 } , and the results are shown on Table 4. The results show that the best value of learning rate is 0.01.
Moreover, we evaluate our model at different dimension d of hidden vectors, and the results are shown on Table 5. The results depict that the performance of NIMGSA is enhanced with the increase of hidden vector dimension. However, when the dimension is larger than 64, there is little increment and the performance remains stable. Therefore, we set the hidden vector dimension at 64 to save the time and space costs.

3.3. Ablation Studies

To evaluate whether the components in our proposed model are necessary, we conduct ablation studies by removing individual component in our model. As shown in Equation (32), the total loss of NIMGSA consists of two parts: The self-attention loss (i.e., matrix completion loss) Y T F 2 , and the reconstruction loss L r . Hence, we evaluate NIMGSA with the following models.
  • Self-attention: Only use self-attention loss to train the model;
  • Without self-attention: Only use reconstruction loss to train the model.
As seen from Table 6, NIMGSA achieves an AUROC of 0.9354, which is 3.4% higher than the model with a self-attention loss only (0.9332), and 4.9% higher than the model without a self-attention loss (0.8916). NIMGSA achieves an AUPR of 0.4567, which is 21.2% higher than the model with a self-attention loss only (0.3768), and 34.6% higher than the model without self-attention loss (0.3392). In summary, both the self-attention loss (i.e., matrix completion loss) Y T F 2 , and the reconstruction loss L r , are essential for NIMGSA. Therefore, NIMGSA is a powerful model combining neural inductive matrix completion, graph autoencoders, and self-attention mechanism, to enhance the preciseness, robustness, and generalization of miRNA-disease associations prediction.

3.4. Case Studies

Case studies are conducted to validate the capability of NIMGSA to predict unknown miRNA-disease associations. The associations in our benchmark dataset are obtained in HMDD v2.0. We adopt three other databases to confirm the predicted miRNA-disease associations which are not included in HMDD v2.0. These three databases are dbDEMC v2.0 [62] (http://www.picb.ac.cn/dbDEMC (accessed on 13 December 2021)), miR2Disease [63] (http://www.mir2disease.org/ (accessed on 13 December 2021)), and HMDD v3.0 [64] (http://www.cuilab.cn/hmdd (accessed on 13 December 2021)). We listed the predictions of potential lncRNA-disease associations with respect to all diseases in Supplementary Table S1. The unknown disease-related lncRNAs of a disease are ranked by their predicted scores. In this paper, we adopt case studies on miRNAs associated with esophageal neoplasms, breast neoplasms, and lung neoplasms.
Esophageal neoplasms (i.e., esophageal cancer) is a major malignant cancer in digestive system [65]. NIMGSA is applied to predict potential miRNAs related to esophageal neoplasms. Supplementary Table S2 lists the top 50 predicted miRNAs associated with esophageal neoplasms. All top 50 predicted miRNAs associated with esophageal neoplasms have been confirmed. Table 7 lists the top 10 predicted miRNAs associated with esophageal neoplasms. For instance, miRNA hsa-mir-125b prevents the progression of esophageal squamous cell carcinoma through the p38-MAPK signaling pathway [66]. MiRNA hsa-mir-17 and hsa-mir-18a are prognostic indicators in esophageal squamous cell carcinoma [67]. MiRNA hsa-miR-16 induces the suppression of cell apoptosis while promoting proliferation in esophageal squamous cell carcinoma [68].
Breast neoplasms (i.e., breast cancer) is the most commonly diagnosed cancer among females worldwide [65]. NIMGSA is applied to predict potential miRNAs related to breast neoplasms. Supplementary Table S3 lists the top 50 predicted miRNAs associated with Breast Neoplasms. A total of 49 of the top 50 predicted miRNAs associated with breast neoplasms have been confirmed. Table 8 lists the top 10 predicted miRNAs associated with breast neoplasms. For instance, miRNA hsa-mir-15b targets the 3′-untranslated region of MTSS1 (metastasis suppressor protein 1), and the low abundance of MTSS1 correlates with a poor patient prognosis of breast neoplasms [69]. MiRNA hsa-mir-192 causes breast cancer cell growth arrest [70]. MiRNA hsa-miR-106a is significantly over-expressed in the breast tumor specimens compared with those in normal controls [71].
Lung neoplasms (i.e., lung cancer) is a major malignant cancer in the respiratory system [65]. NIMGSA is applied to predict potential miRNAs related to lung neoplasms. Supplementary Table S4 lists the top 50 predicted miRNAs associated with lung neoplasms. All top 50 predicted miRNAs associated with lung neoplasms have been confirmed. Table 9 lists the top 10 predicted miRNAs associated with lung neoplasms. For instance, miRNA hsa-mir-15a inhibits metastasis and lipid metabolism by suppressing histone acetylation in lung neoplasms [72]. MiRNA hsa-mir-106b plays a tumorigenesis role in non-small cell lung cancer progression by down-regulating BTG3 expression, which may lead to a novel insight to the potential biomarker and novel therapeutic strategies for non-small cell lung cancer patients [73]. MiRNA hsa-miR-16 regulates proliferation and invasion of lung cancer cells via the ERK/MAPK signaling pathway by targeted inhibition of MAPK kinase 1 (MEK1) [74].

4. Conclusions

Predicting potential miRNA-disease associations is important for understanding the pathogenesis of human diseases. Thus, it is crucial to infer candidate disease-related miRNAs for the scientific discovery of protecting human health. In this paper, we propose a neural network model, NIMGSA, which incorporates inductive matrix completion and graph autoencoders to detect potential miRNA-disease associations. Label propagation can be simulated through the reconstruction procedure of graph autoencoders. Meanwhile, neural inductive matrix completion algorithm not only adds collaborative training to label propagation, but also learns representations in miRNA space and disease space effeciently. Graph autoencoder is a powerful graph representation learning model that ensures the low-rank property of learned representations. Hence, the optimal score matrix can be obtained simply by the linear combination of reconstructed association matrices through GAE on miRNA graph and GAE on disease graph. NIMGSA implements self-attention mechanism through neural inductive matrix completion on two graph autoencoders, which provides theoretical analysis and biological application to enhance the performance of self-attention mechanism. Experiments demonstrated that NIMGSA is superior to the current miRNA-disease associations prediction methods in a series of statistical metrics, such as AUROC, AUPR, and Matthews correlation coefficient. Ablation studies indicate the superiority of our proposed architecture of networks. Case studies on three diseases (esophageal neoplasms, breast neoplasms, and lung neoplasms) indicate that NIMGSA is able to select candidate disease-related miRNAs.
Compared with existing miRNA-disease associations prediction methods, NIMGSA adopts an end-to-end neural network model to integrate inductive matrix completion, self-attention mechanism, and graph autoencoders. This data-driven end-to-end deep learning model not only improves the robustness and preciseness of predicting potential miRNA-disease associations, but also provides a general way for link prediction tasks of other biological entities.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/biom12010064/s1, Table S1: Predictions of potential miRNA-disease association. Table S2: Top 50 predicted miRNAs associated with esophageal neoplasms. Table S3: Top 50 predicted miRNAs associated with breast neoplasms. Table S4: Top 50 predicted miRNAs associated with lung neoplasms.

Author Contributions

Conceptualization, H.Z., C.J. and Z.S.; methodology, C.J., Z.S. and H.Z; software, C.J., Z.S. and K.L.; writing—original draft preparation, C.J., Z.S. and K.L.; writing—review and editing, H.Z.; supervision, H.Z.; funding acquisition, H.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 through Grants (no.61973174).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The miRNA-disease association dataset is retrieved from the HMDD v2.0 database (http://www.cuilab.cn/static/hmdd3/data/hmdd2.zip (accessed on 13 December 2021)). The similarity data can be downloaded from the MISIM database (http://www.cuilab.cn/files/images/cuilab/misim.zip (accessed on 13 December 2021)). The Medical Subject Headings (MeSH) database (https://www.ncbi.nlm.nih.gov/ (accessed on 13 December 2021)) is commonly used for descripting relationship among human diseases. Case studies are conducted using dbDEMC 2.0 (http://www.picb.ac.cn/dbDEMC (accessed on 13 December 2021)), miR2Disease (http://www.mir2disease.org/ (accessed on 13 December 2021)), and HMDD v3.0 (http://www.cuilab.cn/hmdd (accessed on 13 December 2021)). Source code of our model is available at https://github.com/zhanglabNKU/NIMGSA (accessed on 13 December 2021).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ambros, V. The functions of animal microRNAs. Nature 2004, 431, 350–355. [Google Scholar] [CrossRef]
  2. Bartel, D.P. MicroRNAs: Genomics, Biogenesis, Mechanism, and Function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [Green Version]
  3. Chen, X.; Xie, D.; Zhao, Q.; You, Z.H. MicroRNAs and complex diseases: From experimental results to computational models. Brief. Bioinform. 2017, 20, 515–539. [Google Scholar] [CrossRef] [PubMed]
  4. Wiemer, E.A. The role of microRNAs in cancer: No small matter. Eur. J. Cancer 2007, 43, 1529–1544. [Google Scholar] [CrossRef]
  5. Welch, C.; Chen, Y.; Stallings, R.L. MicroRNA-34a functions as a potential tumor suppressor by inducing apoptosis in neuroblastoma cells. Oncogene 2007, 26, 5017–5022. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Zeng, X.; Zhang, X.; Zou, Q. Integrative approaches for predicting microRNA function and prioritizing disease-related microRNA using biological interaction networks. Brief. Bioinform. 2015, 17, 193–203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Zou, Q.; Li, J.; Song, L.; Zeng, X.; Wang, G. Similarity computation strategies in the microRNA-disease network: A survey. Brief Funct. Genom. 2016, 15, 55–64. [Google Scholar] [CrossRef] [PubMed]
  8. Belkin, M.; Niyogi, P.; Sindhwani, V. Manifold Regularization: A Geometric Framework for Learning from Labeled and Unlabeled Examples. J. Mach. Learn. Res. 2006, 7, 2399–2434. [Google Scholar]
  9. Candès, E.; Recht, B. Exact matrix completion via convex optimization. Found. Comput. Math. 2009, 9, 717. [Google Scholar] [CrossRef] [Green Version]
  10. Xia, Z.; Wu, L.Y.; Zhou, X.; Wong, S.T.C. Semi-supervised drug-protein interaction prediction from heterogeneous biological spaces. BMC Syst. Biol. 2010, 4, S6. [Google Scholar] [CrossRef] [Green Version]
  11. You, Z.H.; Lei, Y.K.; Gui, J.; Huang, D.S.; Zhou, X. Using manifold embedding for assessing and predicting protein interactions from high-throughput experimental data. Bioinformatics 2010, 26, 2744–2751. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Chen, X.; Yan, G.Y. Novel human lncRNA—Disease association inference based on lncRNA expression profiles. Bioinformatics 2013, 29, 2617–2624. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Chen, X.; Huang, L. LRSSLMDA: Laplacian Regularized Sparse Subspace Learning for MiRNA-Disease Association prediction. PLoS Comput. Biol. 2017, 13, e1005912. [Google Scholar] [CrossRef] [PubMed]
  14. Xiao, Q.; Luo, J.; Liang, C.; Cai, J.; Ding, P. A graph regularized non-negative matrix factorization method for identifying microRNA-disease associations. Bioinformatics 2018, 34, 239–248. [Google Scholar] [CrossRef] [PubMed]
  15. Natarajan, N.; Dhillon, I.S. Inductive matrix completion for predicting gene–disease associations. Bioinformatics 2014, 30, i60–i68. [Google Scholar] [CrossRef] [Green Version]
  16. Fu, G.; Wang, J.; Domeniconi, C.; Yu, G. Matrix factorization-based data fusion for the prediction of lncRNA–disease associations. Bioinformatics 2017, 34, 1529–1537. [Google Scholar] [CrossRef]
  17. Lu, C.; Yang, M.; Luo, F.; Wu, F.X.; Li, M.; Pan, Y.; Li, Y.; Wang, J. Prediction of lncRNA–disease associations based on inductive matrix completion. Bioinformatics 2018, 34, 3357–3364. [Google Scholar] [CrossRef] [Green Version]
  18. Chen, X.; Wang, L.; Qu, J.; Guan, N.N.; Li, J.Q. Predicting miRNA-disease association based on inductive matrix completion. Bioinformatics 2018, 34, 4256–4265. [Google Scholar] [CrossRef]
  19. Li, J.Q.; Rong, Z.H.; Chen, X.; Yan, G.Y.; You, Z.H. MCMDA: Matrix completion for miRNA-disease association prediction. Oncotarget 2017, 8, 21187–21199. [Google Scholar] [CrossRef] [Green Version]
  20. Yu, S.P.; Liang, C.; Xiao, Q.; Li, G.H.; Ding, P.J.; Luo, J.W. 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]
  21. Chen, X.; Sun, L.G.; Zhao, Y. NCMCMDA: miRNA-disease association prediction through neighborhood constraint matrix completion. Brief. Bioinform. 2021, 22, 485–496. [Google Scholar] [CrossRef] [PubMed]
  22. Chen, X.; Liu, M.X.; Yan, G.Y. RWRMDA: Predicting novel human microRNA–disease associations. Mol. BioSyst. 2012, 8, 2792–2798. [Google Scholar] [CrossRef] [PubMed]
  23. Xuan, P.; Han, K.; Guo, M.; Guo, Y.; Li, J.; Ding, J.; Liu, Y.; Dai, Q.; Li, J.; Teng, Z.; et al. Prediction of microRNAs Associated with Human Diseases Based on Weighted k Most Similar Neighbors. PLoS ONE 2013, 8, e70204. [Google Scholar] [CrossRef]
  24. Chen, X.; Yin, J.; Qu, J.; Huang, L. MDHGI: Matrix Decomposition and Heterogeneous Graph Inference for miRNA-disease association prediction. PLoS Comput. Biol. 2018, 14, e1006418. [Google Scholar] [CrossRef] [PubMed]
  25. Zeng, X.; Liu, L.; Lu, L.; Zou, Q. Prediction of potential disease-associated microRNAs using structural perturbation method. Bioinformatics 2018, 34, 2425–2432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Chen, X.; Xie, D.; Wang, L.; Zhao, Q.; You, Z.H.; Liu, H. BNPMDA: Bipartite Network Projection for MiRNA–Disease Association prediction. Bioinformatics 2018, 34, 3178–3186. [Google Scholar] [CrossRef] [Green Version]
  27. Xuan, P.; Han, K.; Guo, Y.; Li, J.; Li, X.; Zhong, Y.; Zhang, Z.; Ding, J. Prediction of potential disease-associated microRNAs based on random walk. Bioinformatics 2015, 31, 1805–1815. [Google Scholar] [CrossRef]
  28. Wang, L.; You, Z.H.; Huang, Y.A.; Huang, D.S.; Chan, K.C.C. An efficient approach based on multi-sources information to predict circRNA–disease associations using deep convolutional neural network. Bioinformatics 2019, 36, 4038–4046. [Google Scholar] [CrossRef]
  29. Peng, W.; Lan, W.; Yu, Z.; Wang, J.; Pan, Y. A Framework for Integrating Multiple Biological Networks to Predict MicroRNA-Disease Associations. IEEE Trans. NanoBiosci. 2016, 16, 100–107. [Google Scholar] [CrossRef]
  30. Liu, Y.; Zeng, X.; He, Z.; Zou, Q. Inferring MicroRNA-Disease Associations by Random Walk on a Heterogeneous Network with Multiple Data Sources. IEEE/ACM Trans. Comput. Biol. Bioinform. 2017, 14, 905–915. [Google Scholar] [CrossRef]
  31. Xiao, Q.; Zhang, N.; Luo, J.; Dai, J.; Tang, X. Adaptive multi-source multi-view latent feature learning for inferring potential disease-associated miRNAs. Brief. Bioinform. 2020, 22, 2043–2057. [Google Scholar] [CrossRef]
  32. Ha, J.; Park, C. MLMD: Metric Learning for Predicting MiRNA-Disease Associations. IEEE Access 2021, 9, 78847–78858. [Google Scholar] [CrossRef]
  33. Yu, D.L.; Yu, Z.G.; Han, G.S.; Li, J.; Anh, V. Heterogeneous Types of miRNA-Disease Associations Stratified by Multi-Layer Network Embedding and Prediction. Biomedicines 2021, 9, 1152. [Google Scholar] [CrossRef]
  34. Kang, C.; Zhang, H.; Liu, Z.; Huang, S.; Yin, Y. LR-GNN: A graph neural network based on link representation for predicting molecular associations. Brief. Bioinform. 2021. [Google Scholar] [CrossRef] [PubMed]
  35. Zeng, X.; Wang, W.; Deng, G.; Bing, J.; Zou, Q. Prediction of Potential Disease-Associated MicroRNAs by Using Neural Networks. Mol. Ther. Nucl. Acids 2019, 16, 566–575. [Google Scholar] [CrossRef] [Green Version]
  36. Chen, X.; Li, T.H.; Zhao, Y.; Wang, C.C.; Zhu, C.C. Deep-belief network for predicting potential miRNA-disease associations. Brief. Bioinform. 2021, 22, bbaa186. [Google Scholar] [CrossRef] [PubMed]
  37. Ji, C.; Gao, Z.; Ma, X.; Wu, Q.; Ni, J.; Zheng, C. AEMDA: Inferring miRNA-disease associations based on deep autoencoder. Bioinformatics 2020, 37, 66–72. [Google Scholar] [CrossRef] [PubMed]
  38. Tang, X.; Luo, J.; Shen, C.; Lai, Z. Multi-view Multichannel Attention Graph Convolutional Network for miRNA-disease association prediction. Brief. Bioinform. 2021, 22, 1–16. [Google Scholar] [CrossRef]
  39. Scarselli, F.; Gori, M.; Tsoi, A.C.; Hagenbuchner, M.; Monfardini, G. The Graph Neural Network Model. IEEE Trans. Neural Net. 2009, 20, 61–80. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Xuan, P.; Pan, S.; Zhang, T.; Liu, Y.; Sun, H. Graph Convolutional Network and Convolutional Neural Network Based Method for Predicting lncRNA-Disease Associations. Cells 2019, 8, 1012. [Google Scholar] [CrossRef] [Green Version]
  41. Shi, Z.; Zhang, H.; Jin, C.; Quan, X.; Yin, Y. A representation learning model based on variational inference and graph autoencoder for predicting lncRNA-disease associations. BMC Bioinform. 2021, 22, 136. [Google Scholar] [CrossRef]
  42. Jin, C.; Shi, Z.; Zhang, H.; Yin, Y. Predicting lncRNA–protein interactions based on graph autoencoders and collaborative training. In Proceedings of the IEEE International Conference on Bioinformatics and Biomedicine (BIBM), Houston, TX, USA, 9–12 December 2021. [Google Scholar]
  43. 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] [PubMed]
  44. Li, Z.; Li, J.; Nie, R.; You, Z.H.; Bao, W. A graph auto-encoder model for miRNA-disease associations prediction. Brief. Bioinform. 2020, 22, bbaa240. [Google Scholar] [CrossRef]
  45. Kipf, T.N.; Welling, M. Semi-Supervised Classification with Graph Convolutional Networks. In Proceedings of the International Conference on Learning Representations (ICLR), Toulon, France, 24–26 April 2017. [Google Scholar]
  46. Geng, Z.; Guo, M.H.; Chen, H.; Li, X.; Wei, K.; Lin, Z. Is Attention Better Than Matrix Decomposition? In Proceedings of the International Conference on Learning Representations (ICLR), Vienna, Austria, 4 May 2021. [Google Scholar]
  47. Li, Y.; Qiu, C.; Tu, J.; Geng, B.; Yang, J.; Jiang, T.; Cui, Q. HMDD v2.0: A database for experimentally supported human microRNA and disease associations. Nucleic Acids Res. 2013, 42, D1070–D1074. [Google Scholar] [CrossRef]
  48. Wang, D.; Wang, J.; Lu, M.; Song, F.; Cui, Q. Inferring the human microRNA functional similarity and functional network based on microRNA-associated diseases. Bioinformatics 2010, 26, 1644–1650. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Wang, F.; Zhang, C. Label Propagation through Linear Neighborhoods. IEEE Trans. Knowl. Data Eng. 2008, 20, 55–67. [Google Scholar] [CrossRef]
  50. Li, Q.; Han, Z.; Wu, X.M. Deeper insights into graph convolutional networks for semi-supervised learning. In Proceedings of the Association for the Advancement of Artificial Intelligence (AAAI), New Orleans, LA, USA, 2–7 February 2018. [Google Scholar]
  51. Li, Q.; Wu, X.; Liu, H.; Zhang, X.; Guan, Z. Label Efficient Semi-Supervised Learning via Graph Filtering. In Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR), Los Angeles, CA, USA, 16–20 June 2019; pp. 9574–9583. [Google Scholar]
  52. Qu, M.; Bengio, Y.; Tang, J. GMNN: Graph Markov Neural Networks. In Proceedings of the Machine Learning Research, Los Angeles, CA, USA, 9–15 June 2019; Volume 97, pp. 5241–5250. [Google Scholar]
  53. Lee, D.D.; Seung, H.S. Algorithms for Non-Negative Matrix Factorization. In Proceedings of the 13th International Conference on Neural Information Processing Systems, Denver, CO, USA, 28–30 November 2000; pp. 535–541. [Google Scholar]
  54. Bahdanau, D.; Cho, K.; Bengio, Y. Neural machine translation by jointly learning to align and translate. In Proceedings of the International Conference on Learning Representations, San Diego, CA, USA, 7–9 May 2015. [Google Scholar]
  55. Vaswani, A.; Shazeer, N.; Parmar, N.; Uszkoreit, J.; Jones, L.; Gomez, A.N.; Kaiser, L.u.; Polosukhin, I. Attention is All you Need. In Advances in Neural Information Processing Systems 30; MIT Press: Cambridge, MA, USA, 2017; pp. 5998–6008. [Google Scholar]
  56. Jin, C.; Gao, J.; Shi, Z.; Zhang, H. ATTCry: Attention-based neural network model for protein crystallization prediction. Neurocomputing 2021, 463, 265–274. [Google Scholar] [CrossRef]
  57. Ji, P.; Zhang, T.; Li, H.; Salzmann, M.; Reid, I. Deep Subspace Clustering Network. In Proceedings of the Advances in Neural Information Processing Systems, Los Angeles, CA, USA, 4–9 December 2017; Volume 30, pp. 23–32. [Google Scholar]
  58. Jing, L.; Zbontar, J.; LeCun, Y. Implicit Rank-Minimizing Autoencoder. In Proceedings of the Advances in Neural Information Processing Systems, Virtual-only, 6–12 December 2020; Volume 33, pp. 1–11. [Google Scholar]
  59. Kipf, T.; Welling, M. Variational Graph Auto-Encoders. In Proceedings of the NeurIPS Workshop on Bayesian Deep Learning, Barcelona, Spain, 5–10 December 2016. [Google Scholar]
  60. Kingma, D.P.; Ba, J. Adam: A method for stochastic optimization. In Proceedings of the International Conference on Learning Representations—ICLR, Banff, AB, Canada, 14–16 April 2014. [Google Scholar]
  61. Srivastava, N.; Hinton, G.; Krizhevsky, A.; Sutskever, I.; Salakhutdinov, R. Dropout: A simple way to prevent neural networks from overfitting. J. Mach. Learn. Res. 2014, 15, 1929–1958. [Google Scholar]
  62. Yang, Z.; Wu, L.; Wang, A.; Tang, W.; Zhao, Y.; Zhao, H.; Teschendorff, A.E. dbDEMC 2.0: Updated database of differentially expressed miRNAs in human cancers. Nucleic Acids Res. 2016, 45, D812–D818. [Google Scholar] [CrossRef]
  63. 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. 2008, 37, D98–D104. [Google Scholar] [CrossRef] [Green Version]
  64. 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. 2018, 47, D1013–D1017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.; Torre, L.; Jemal, A. Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries: Global Cancer Statistics 2018. CA 5Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [Green Version]
  66. Cheng, C.; Mao, Q.; Shi, M.; Lu, H.; Liu, Y. miR-125b prevent the progression of esophageal squamous cell carcinoma through the p38-MAPK signaling pathway. J. Gastrointest. Oncol. 2020, 11, 1113–1122. [Google Scholar] [CrossRef]
  67. Xu, X.L.; Jiang, Y.H.; Feng, J.G.; Su, D.; Chen, P.C.; Mao, W.M. MicroRNA-17, MicroRNA-18a, and MicroRNA-19a Are Prognostic Indicators in Esophageal Squamous Cell Carcinoma. Ann. Thorac. Surg. 2014, 97, 1037–1045. [Google Scholar] [CrossRef]
  68. Zhu, Y.; Yang, X.; Niu, H.; Chen, Y. MiR-16 Induced the Suppression of Cell Apoptosis While Promote Proliferation in Esophageal Squamous Cell Carcinoma. Cell. Physiol. Biochem. 2014, 33, 1340–1348. [Google Scholar] [CrossRef]
  69. Kedmi, M.; Ben-Chetrit, N.; Korner, C.; Mancini, M.; Ben-Moshe, N.B.; Lauriola, M.; Lavi, S.; Biagioni, F.; Carvalho, S.; Cohen-Dvashi, H. EGF induces microRNAs that target suppressors of cell migration: miR-15b targets MTSS1 in breast cancer. Sci. Signal. 2015, 8, ra29. [Google Scholar] [CrossRef] [Green Version]
  70. Hu, F.; Meng, X.; Tong, Q.; Liang, L.; Xiang, R.; Zhu, T.; Yang, S. BMP-6 inhibits cell proliferation by targeting microRNA-192 in breast cancer. Biochim. D Acta (BBA) Mol. Basis Dis. 2013, 1832, 2379–2390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Wang, F.; Zheng, Z.; Guo, J.; Ding, X. Correlation and quantitation of microRNA aberrant expression in tissues and sera from patients with breast tumor. Gynecol. Oncol. 2010, 119, 586–593. [Google Scholar] [CrossRef] [PubMed]
  72. Ni, Y.; Yang, Y.; Ran, J.; Zhang, L.; Zhang, L. miR-15a-5p inhibits metastasis and lipid metabolism by suppressing histone acetylation in lung cancer. Free Radic. Biol. Med. 2020, 161, 150–162. [Google Scholar] [CrossRef]
  73. Wei, K.; Pan, C.; Yao, G.; Liu, B.; Ma, T.; Xia, Y.; Jiang, W.; Chen, L.; Chen, Y. MiR-106b-5p Promotes Proliferation and Inhibits Apoptosis by Regulating BTG3 in Non-Small Cell Lung Cancer. Cell. Physiol. Biochem. 2017, 44, 1545–1558. [Google Scholar] [CrossRef]
  74. Chen, T.M.; Xiao, Q.; Wang, X.J.; Wang, Z.Q.; Chen, S.L. miR-16 regulates proliferation and invasion of lung cancer cells via the ERK/MAPK signaling pathway by targeted inhibition of MAPK kinase 1 (MEK1). J. Int. Med. Res. 2019, 47, 5194–5204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Flowchart of similarity computation.
Figure 1. Flowchart of similarity computation.
Biomolecules 12 00064 g001
Figure 2. Illustration of NIMGSA. GAEm and GAEd represent graph autoencoders on the miRNA graph and disease graph respectively.
Figure 2. Illustration of NIMGSA. GAEm and GAEd represent graph autoencoders on the miRNA graph and disease graph respectively.
Biomolecules 12 00064 g002
Figure 3. Computation procedure of NIMGSA. NIMC denotes Neural Inductive Matrix Completion.
Figure 3. Computation procedure of NIMGSA. NIMC denotes Neural Inductive Matrix Completion.
Biomolecules 12 00064 g003
Figure 4. ROC and PR curves of different methods.
Figure 4. ROC and PR curves of different methods.
Biomolecules 12 00064 g004
Table 1. Mean values and standard deviations of AUROC and AUPR, compared with different methods.
Table 1. Mean values and standard deviations of AUROC and AUPR, compared with different methods.
METHODAUROCAUPR
IMCMDA0.8329 ± 0.00110.2785 ± 0.0029
SPM0.8960 ± 0.00700.2464 ± 0.0054
NIMCGCN0.9279 ± 0.00060.3943 ± 0.0054
MCLPMDA0.9292 ± 0.00690.4387 ± 0.0106
GAEMDA0.9332 ± 0.00050.4142 ± 0.0034
NIMGSA0.9354 ± 0.00470.4567 ± 0.0147
Table 2. Binary classification metrics of different methods on Dataset2. Sp denotes specificity. Sn denotes sensitivity. Acc denotes accuracy. Pre denotes precision. F1 denotes F1-score. Mcc denotes Matthews correlation coefficient.
Table 2. Binary classification metrics of different methods on Dataset2. Sp denotes specificity. Sn denotes sensitivity. Acc denotes accuracy. Pre denotes precision. F1 denotes F1-score. Mcc denotes Matthews correlation coefficient.
SPECMETHODSENACCPREF1-ScoreMCC
0.99IMCMDA0.26280.96920.43650.32810.3239
SPM0.15510.96610.31370.20750.2048
NIMCGCN0.30390.97030.47250.36990.3645
MCLPMDA0.35670.97190.51270.42070.4138
GAEMDA0.36500.97210.51860.42840.4213
NIMGSA0.37180.97230.52290.43460.4273
Table 3. AUROC and AUPR at different α .
Table 3. AUROC and AUPR at different α .
α 0.10.30.50.70.9
AUROC0.91190.92890.93540.93380.9312
AUPR0.36480.42550.45670.45560.4509
Table 4. AUROC and AUPR at a different learning rate.
Table 4. AUROC and AUPR at a different learning rate.
lr0.0010.010.050.1
AUROC0.91930.93540.76930.5557
AUPR0.40770.45670.27910.0709
Table 5. AUROC and AUPR at a different dimension of hidden vectors.
Table 5. AUROC and AUPR at a different dimension of hidden vectors.
DIMENSION163264128
AUROC0.90120.92280.93540.9357
AUPR0.36420.41270.45670.4589
Table 6. Ablation studies.
Table 6. Ablation studies.
ModelsAUROCAUPR
Self-attention0.90460.3768
Without self-attention0.89160.3392
NIMGSA0.93540.4567
Table 7. Top 10 predicted miRNAs associated with esophageal neoplasms.
Table 7. Top 10 predicted miRNAs associated with esophageal neoplasms.
MiRNA NAMEEVIDENCE
hsa-mir-125bdbDEMC v2.0; HMDD v3.0
hsa-mir-17dbDEMC v2.0
hsa-mir-16dbDEMC v2.0
hsa-mir-18adbDEMC v2.0
hsa-mir-19bdbDEMC v2.0
hsa-mir-29adbDEMC v2.0
hsa-mir-222dbDEMC v2.0
hsa-mir-1dbDEMC v2.0
hsa-mir-29bdbDEMC v2.0
hsa-mir-200bdbDEMC v2.0
Table 8. Top 10 predicted miRNAs associated with breast neoplasms.
Table 8. Top 10 predicted miRNAs associated with breast neoplasms.
MiRNA NAMEEVIDENCE
hsa-mir-142dbDEMC v2.0; HMDD v3.0
hsa-mir-15bdbDEMC v2.0; HMDD v3.0
hsa-mir-192dbDEMC v2.0; HMDD v3.0
hsa-mir-106adbDEMC v2.0; HMDD v3.0
hsa-mir-150dbDEMC v2.0; HMDD v3.0
hsa-mir-130adbDEMC v2.0; HMDD v3.0
hsa-mir-30edbDEMC v2.0; HMDD v3.0
hsa-mir-92bdbDEMC v2.0; HMDD v3.0
hsa-mir-192bdbDEMC v2.0; miR2Disease; HMDD v3.0
hsa-mir-372dbDEMC v2.0; HMDD v3.0
Table 9. Top 10 predicted miRNAs associated with lung neoplasms.
Table 9. Top 10 predicted miRNAs associated with lung neoplasms.
MiRNA NAMEEVIDENCE
hsa-mir-16dbDEMC v2.0; miR2Disease; HMDD v3.0
hsa-mir-15adbDEMC v2.0; HMDD v3.0
hsa-mir-106bdbDEMC v2.0; miR2Disease; HMDD v3.0
hsa-mir-141dbDEMC v2.0; miR2Disease; HMDD v3.0
hsa-mir-15bdbDEMC v2.0; HMDD v3.0
hsa-mir-122dbDEMC v2.0; HMDD v3.0
hsa-mir-429dbDEMC v2.0; miR2Disease; HMDD v3.0
hsa-mir-20bdbDEMC v2.0; HMDD v3.0
hsa-mir-23bdbDEMC v2.0; HMDD v3.0
hsa-mir-130adbDEMC v2.0; miR2Disease; HMDD v3.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jin, C.; Shi, Z.; Lin, K.; Zhang, H. Predicting miRNA-Disease Association Based on Neural Inductive Matrix Completion with Graph Autoencoders and Self-Attention Mechanism. Biomolecules 2022, 12, 64. https://doi.org/10.3390/biom12010064

AMA Style

Jin C, Shi Z, Lin K, Zhang H. Predicting miRNA-Disease Association Based on Neural Inductive Matrix Completion with Graph Autoencoders and Self-Attention Mechanism. Biomolecules. 2022; 12(1):64. https://doi.org/10.3390/biom12010064

Chicago/Turabian Style

Jin, Chen, Zhuangwei Shi, Ken Lin, and Han Zhang. 2022. "Predicting miRNA-Disease Association Based on Neural Inductive Matrix Completion with Graph Autoencoders and Self-Attention Mechanism" Biomolecules 12, no. 1: 64. https://doi.org/10.3390/biom12010064

APA Style

Jin, C., Shi, Z., Lin, K., & Zhang, H. (2022). Predicting miRNA-Disease Association Based on Neural Inductive Matrix Completion with Graph Autoencoders and Self-Attention Mechanism. Biomolecules, 12(1), 64. https://doi.org/10.3390/biom12010064

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