 Research
 Open Access
 Published:
MiRNAdisease association prediction via hypergraph learning based on highdimensionality features
BMC Medical Informatics and Decision Making volume 21, Article number: 133 (2021)
Abstract
Background
MicroRNAs (miRNAs) have been confirmed to have close relationship with various human complex diseases. The identification of diseaserelated miRNAs provides great insights into the underlying pathogenesis of diseases. However, it is still a big challenge to identify which miRNAs are related to diseases. As experimental methods are in general expensive and time‐consuming, it is important to develop efficient computational models to discover potential miRNAdisease associations.
Methods
This study presents a novel prediction method called HFHLMDA, which is based on highdimensionality features and hypergraph learning, to reveal the association between diseases and miRNAs. Firstly, the miRNA functional similarity and the disease semantic similarity are integrated to form an informative highdimensionality feature vector. Then, a hypergraph is constructed by the KNearestNeighbor (KNN) method, in which each miRNAdisease pair and its k most relevant neighbors are linked as one hyperedge to represent the complex relationships among miRNAdisease pairs. Finally, the hypergraph learning model is designed to learn the projection matrix which is used to calculate uncertain miRNAdisease association score.
Result
Compared with four stateoftheart computational models, HFHLMDA achieved best results of 92.09% and 91.87% in leaveoneout cross validation and fivefold cross validation, respectively. Moreover, in case studies on Esophageal neoplasms, Hepatocellular Carcinoma, Breast Neoplasms, 90%, 98%, and 96% of the top 50 predictions have been manually confirmed by previous experimental studies.
Conclusion
MiRNAs have complex connections with many human diseases. In this study, we proposed a novel computational model to predict the underlying miRNAdisease associations. All results show that the proposed method is effective for miRNA–disease association predication.
Background
MicroRNAs (miRNAs) are endogenous noncoding singlestranded RNA molecules that play important roles in eukaryotic gene expression through posttranscriptional regulation [1,2,3]. Functional studies indicate that miRNA plays a significant role in manifold biological processes, such as cell proliferation, stem cell maintenance, immune responses and so on [4,5,6]. Dysregulation of miRNA expression and function is reported in various diseases including cancer, metabolic disorders as well as neurological disorders [7]. Therefore, identifying diseaserelated miRNAs is important to treat, diagnose, and prevent human complex diseases [8, 9].
Generally, researchers use biological experimental methods such as quantitative reverse transcription, microarray analysis, or deep sequencing of small RNAs to explore miRNAs that are differentially expressed in a disease state. For example, Pan et al. used microarray analysis and found that miR130a3p, miR4245p, miR5745p, and miR146a presented significant difference between tuberculous meningitis and healthy controls [10]. However, experimental identification of diseaserelated miRNAs by existing techniques is expensive and timeconsuming. So, based on vast amount of biological data about miRNAs, researchers have developed computational methods for predicting miRNAdisease associations [11,12,13,14,15,16,17,18,19,20,21], which can select most promising miRNAs for further analysis and hence decrease the number of the experiments.
For predicting diseaserelated miRNAs, many methods are based on a credible assumption that functionally similar miRNAs tend to have associations with phenotypically similar diseases and vice versa. Xiao et al. proposed a method called GRNMF, which based on graph regularized nonnegative matrix factorization from the similarity and association perspective of miRNAs and diseases to discover potential associations [22]. Liu et al. proposed the method for predicting miRNA–disease associations by performing random walks on heterogeneous omics data [23]. You et al. presented the prediction model of PBMDA by constructing a heterogeneous graph consisting of three interlinked subgraphs, and performing a depthfirst search algorithm on the heterogeneous network to infer diseaserelated miRNAs [24]. PBMDA integrated different types of heterogeneous biological datasets, so it can be applied to the new diseases/miRNAs without known associated miRNAs/diseases. Subsequently, Chen et al. proposed a novel method based on Hybrid Approach for MiRNADisease Association prediction (HAMDA) [25]. They considered network structure, information propagation, and node attribution, and used the hybrid graphbased recommendation algorithm to uncover diseaserelated miRNAs. In addition, Chen et al. devised a computational approach by Graphlet Interaction to predict diseaserelated miRNAs (GIMDA) [26]. In this method, graphlet interaction was utilized to analyze the complex relationships between two nodes in a graph. However, HAMDA and GIMDA are not applicable to predicting a new association between a new miRNA and a new disease. Furthermore, Chen et al. developed a method of Graph Regression for MiRNADisease Association prediction (GRMDA) [27]. The graph regression was synchronously performed in three latent spaces, by using Singular Value Decomposition (SVD) and Partial LeastSquares (PLS) to extract important related attributes and filter the noise. But it is difficulties to choosing parameters in SVD and PLS. Lately, Jiang et al. implemented a improved collaborative filteringbased method to infer miRNAdisease associations (ICFMDA) [28]. They improved collaborative filtering algorithm by combining the similarity matrices, and defined significance SIG between pairs of diseases or miRNAs to predict diseaserelated miRNAs even new diseases without known association.
In addition, several computational models used machine learning to uncover the association between miRNAs and diseases. Xu et al. introduced an approach based on the miRNA target–dysregulated network (MTDN) to prioritize novel disease miRNAs [29]. They applied Support vector machine classifier to miRNAs in the MTDN. However, negative samples required by the classifier are difficult to obtain. To overcome this limitation, Chen et al. introduced a semisupervised method named RLSMDA [30]. It is developed under the framework of regularized least squares and can predict new miRNAs for diseases which do not have any known related miRNAs. Similarly, Luo et al. developed another semisupervised method named KRLSM based on Kronecker regularized least squares [31]. KRLSM integrated different omics data, combined the disease and miRNA space, and used the semisupervised classifier of regularized least squares to predict diseaserelated miRNAs. However, this approach involves multiple parameters and establishing the optimal parameter values remains a challenging problem. Chen et al. designed a method based on restricted Boltzmann machine for predicting miRNAdisease associations [32]. This approach can also predict association types of miRNAdisease pairs, but can not applicable to a new disease with no known associated miRNAs. Furthermore, Chen et al. developed an effective method called HGIMDA [33]. HGIMDA calculated the diseasemiRNA association possibility by investigating all the 3length paths in the constructed heterogeneous graph. Recently, Chen et al. utilized Extreme Gradient Boosting Machine to uncover diseaserelated miRNAs and named EGBMMDA [34]. In this method, based on statistical measures, graph theoretical, and matrix factorization, they constructed an informative feature vector for each miRNAdisease pair and used a decision tree model to predict diseaserelated miRNAs.
Although existing methods have made great contributions to uncover diseaserelated miRNAs, there are still some limitations that could be improved. For example, many methods are difficult to extract the deep feature representation of the multiple kinds of data. In this study, we propose a novel prediction method via hypergraph learning based on highdimensionality features and refer to it as HFHLMDA. Hypergraph learning, which can capture the highorder relationships of samples, has been widely used in clustering, classification and information retrieval tasks. In a hypergraph, an edge connects more than two vertices, thus it can well encode the relationship among more than two vertices. We construct highdimensionality feature vectors for all the miRNAdisease pairs, and utilize KNearestNeighbor (KNN) method to form a hypergraph to predict potential miRNAdisease association. To demonstrate the effectiveness of our method, we apply Leaveoneout cross validation (LOOCV) and fivefold cross validation to measure the prediction performance. We compare our method with four state‐of‐the‐art methods and the results indicate that our method can achieve better performance. In addition, case studies of three common diseases are implemented to further verify the reliability and robustness of HFHLMDA.
Methods
Human MiRNAdisease associations network
The human miRNAdisease associations used in this work come from the HMDDv2.0 [35], which contains 5430 experimentally associations between 495 miRNAs and 383 diseases. Technically, we use an adjacency matrix A with 495 (nm) rows and 383 (nd) columns to clearly describe the relation of each miRNAdisease pairs. The element A(m(i), d(j)) is equal to 1 if miRNA m(i) is verified to be associated with disease d(j), and 0 otherwise. Finally, 5430 entries of matrix A are assigned 1, the rest ones are assigned 0. Our goal is to confirm the uncertain associations between miRNAs and diseases.
MiRNA similarity matrix
Wang et al. developed a method named MISIM for calculating the function similarity scores of miRNA [36]. Here, we directly downloaded the miRNA functional similarity scores from http://www.cuilab.cn/files/images/cuilab/misim.zip. Then, an adjacency matrix SM with 495 rows and 495 columns is built to denote the similarity of miRNAs, in which the larger the SM(m(i), m(j)) is, the more similar m(i) and m(j) are.
However, SM has the problem of sparsity. Sparse matrix is difficult to provide more effective information, which will seriously affect the prediction performance of the computational model. So we calculate the Gaussian interaction profile kernel similarity of miRNAs [37]. Specifically, a binary vector BV(m(i)), i.e. the ith row of matrix A, is recorded as the interaction profiles of miRNA m(i) for representing the associations between m(i) itself and each disease. All known miRNAdisease associations in matrix A will be used to calculate similarity, two miRNAs would likely have greater similarities if they share more disease associations. Thus, the Gaussian interaction profile kernel similarity GKM(m(i), m(j)) of miRNA m(i) and miRNA m(j) is defined as
where γ_{m} is a parameter used to control the kernel bandwidth, which is set as
By integrating SM and GKM, a new complete miRNA similarity matrix SM can be obtained as
Disease similarity matrix
The association between different diseases can be represented by a directed acyclic graph (DAG), which consists of some nodes and links. Each node represents a disease while a link represents the association of two diseases. For a given disease D, DAG = (D, T_{D}, E_{D}), where T_{D} represents its ancestor nodes and itself while E_{D} is the set of corresponding edges. The contribution values of disease d(t) to the semantic value of disease d(i) can be calculated as follows:
where D(d(i)) is the node set in DAG(d(i)) including node d(i) itself. Therefore, the semantic similarity between disease d(i) and d(j) can be defined as follows:
Similarly, we also calculate the Gaussian interaction profile kernel similarity GKD for diseases by the follow formulas
where BV(d(i)) and BV(d(j)) denote the ith column and the jth column of A. At last, the disease similarity matrix SD is obtained by
HFHLMDA
The HFHLMDA model can be separated into three steps (see Fig. 1). First, feature factor construction, in which a feature factor x for each miRNAdisease pair consisting of corresponding rows of SM and SD. Second, hypergraph construction, where a hypergraph G is constructed to formulate the relationship between these feature vectors. Third, hypergraph learning, to learn the projection matrix P, which map the original feature x to the relevance score S = x^{.}P, and thus it can be used to predict the association for the unknown miRNAdisease pair x^{unk}.
Feature factor construction
According to the biological observation that miRNAs with more functional similarity tend to be more associated with similar diseases and vice versa, so the topologic information of miRNA/disease similarity network can be used to construct feature factor directly.
For each miRNA, there are 495 similarity scores. We use similarity scores as features to represent each miRNA by a 495dimensional feature vector. For example, we represent miRNA m(i) by a feature vector, SM(m(i)) = (m_{1}, m_{2}, …, m_{495}), where SM(m(i)) is the ith row vector of SM and represents the similarities between m(i) and all the miRNAs.
For each disease, we can obtain a 383dimensonal feature vector in a similar way to miRNA, SD(d(j)) = (d_{1}, d_{2}, …, d_{383}), where SD(d(j)) is the jth row of matrix SD. Therefore, each miRNAdisease pair can be described by an 878dimensional vector x = (SM(m(i)), SD(d(j))). Furthermore, we consider (SM(m(i)), SD(d(j))) as a positive sample if miRNA m(i) is associated with disease d(j), otherwise as a negative sample. To construct the balanced dataset, the training set have 5,430 positive samples, and an equal number of samples were randomly selected as negative training examples from the pool of unknown associations. It is possible to use unconfirmed miRNAdisease pairs with association as negative samples, from the perspective of probability, because the miRNAdisease pairs we selected as negative samples account for only 5430 ÷ (495 × 383) ≈ 2.86% of all miRNAdisease pairs, which is negligible [38].
Hypergraph construction
Firstly, we briefly introduce the hypergraph learning theory. As a generalization of graph, hypergraph represents the structure of data via measuring the similarity between groups of points. Different from a simple graph, an edge in a hypergraph can connect three or more vertices, it can model highorder relations between their vertices by hyperedges, whose influence can be assessed by properly estimating their weights. Obviously, modeling the highorder relationship among objects can improve the predicting performance significantly. Moreover, the quality of the hypergraph structure plays an important role for data modeling. A well constructed hypergraph structure can represent the data correlation accurately, and leading to better performance.
A hypergraph is defined as G = (V, E, w), where V is a set of vertices, E is a set of hyperedges and each hyperedge e is given a positive weight w(e). The hypergraph G can be denoted by a V ×E incidence matrix H, in which each entry is defined by
The degree of vertex v ∈ V and hyperedge e ∈ E can be respectively represented as:
Accordingly, denote Dv and De as two diagonal matrices of the vertex degrees and the hyperedge degrees, respectively.
Zhou et al. proposed a regularization framework on hypergraph [39], which is defined as
where f is the tobelearned function, Ω(f) is a regularizer on the hypergraph, R_{emp}(f) is an empirical loss, and λ > 0 is the tradeoff parameter. Usually, the empirical loss R_{emp}(f) is defined as
where Y is the label matrix of samples. The regularizer on the hypergraph is defined by
Let Θ = D_{v}^{−(1/2)}HWD_{e}^{−1}H^{T}D_{v}^{−(1/2)}, the normalized cost function can be written as
where Δ = I – Θ, which is a positive semidefinite matrix.
In this study, given a set of training samples {x_{i} i = 1,…, n} ∈ \({\mathbb{R}}\)^{878}, the data matrix X = [x_{1},..., x_{i},..., x_{n}]^{T} ∈ \({\mathbb{R}}\)^{n×878} contains n samples in its rows, the corresponding labels matrix Y = [y_{1},..., y_{2},..., y_{l}] ∈ \({\mathbb{R}}\)^{n×l}, y_{i} is the label vector of the ith class. A miRNAdisease pairs hypergraph \({\mathcal{G}} = \left( {{\mathcal{V}},{ \mathcal{E}},{ \mathcal{W}}} \right)\) is constructed, and its hyperedge is generated based on the KNN algorithm. Concretely, for each vertex v, we search its corresponding k nearest neighbors, and use these nearest neighbors to form a hyperedge e(v). We initialize k as 15 here empirically. An illustration on the hyperedge generation process is shown in Fig. 2. Moreover, the diagonal matrix \({\mathcal{W}}\) denote the weights of the hyperedges. All the hyperedges are initialized with an equal weight, e.g., w(e) = 1/n_{e}, where n_{e} is the number of hyperedges.
Hypergraph learning
The hypergraph learning targets on learning a regularized projection to discriminate different categories. According to Zhang et al. introduction [40], the cost function F for learning the projection matrix P can be formulated as:
where λ and μ are positive parameters, and we empirically set them as 10^{1},10^{0} respectively, which can achieve the best performance. Specifically, hypergraph Laplacian regularizer Ω(P) is calculated as
where function tr(·) returns the trace of matrix. The empirical loss term R_{emp} (P) is defined as
Φ(P) is a l_{2} norm regularizer to avoid overfitting for P, which is defined as:
Consequently, Eq. (17) can be reformed as:
Such problem is a typical Least Square problem which can be efficiently solved, its solution is as follows:
where I is an identity matrix. Based on the learned P, the relevance score of the unknown miRNAdisease pair x^{unk} can be obtained by
Results
Effect of parameters on the performance of HFHLMDA
In this work, we used KNN algorithm to generate hyperedge, one parameters k was included, which represent the number of nearest neighbors of miRNA or disease. In the hypergraph learning section of the Methods, we defined two parameters, namely, λ and μ to balance the items in Eq. (17), the values of λ and μ ranged from 10^{–2}, 10^{–1}, 10^{0}, 10^{1} to 10^{2}. We conducted a series of experiments on the above parameters to acquire the effects of these parameters. The experimental results are shown in Figs. 3 and 4. In Fig. 3, we can see that regardless of how k change, the AUC of fivefold cross validation keep around 0.9187. Thus, for efficiency, we set k = 15. Furthermore, Fig. 4 describes the prediction performances of HFHLMDA with different values of λ and μ. We can see that HFHLMDA obtains the best prediction performance when λ is set to be 10^{1} and μ is set to be 10^{0}.
Performance evaluation
Based on the known miRNA–disease associations in HMDDv2.0 database, two validation schemas were used to evaluate the performance of HFHLMDA: LOOCV and fivefold cross validation. We selected four classical computational methods: EGBMMDA [34], ICFMDA [28], RLSMDA [30], and SACMDA [41] to compete with HFHLMDA in cross validation. Specifically, LOOCV selected a known miRNAdisease association in turn as a test sample, and the rest of the associations were considered as training samples. All unknown associations were used as candidate samples. Considering that the Gaussian interaction profile kernel similarity depend on known miRNAdisease associations, the corresponding value of a test sample in matrix A should be set to 0. The predicted score for the test sample was ranked relative to the scores for candidate samples and, each ranking will take turns as a threshold in each fold, if test ranking was above a given threshold, we obtained a successful prediction made by the model. By changing the threshold, we could calculate the corresponding true positive rate (TPR) and false positive rate (FPR). Furthermore, receiveroperating characteristics (ROC) curve could be drawn according to TPR against FPR. The areas under the ROC curve (AUC) was used to evaluate the whole prediction performance. Figure 5 shows the global LOOCV ROC curves for HFHLMDA and other methods. HFHLMDA, EGBMMDA, ICFMDA, RLSMDA and SACMDA obtained AUCs of 0.9209, 0.9123, 0.9067, 0.8426 and 0.8770, respectively. HFHLMDA achieved the better prediction performance.
As for fivefold cross validation, in order to make the validation more accurate, we repeated fivefold cross validation procedure 100 times. The average AUC values of the five methods (HFHLMDA, EGBMMDA, ICFMDA, RLSMDA, SACMDA) were 0.9187(± 0.0009), 0.9048(± 0.0012), 0.9045(± 0.0008), 0.8569(± 0.0020) and 0.8767(± 0.0011), respectively (see Fig. 6). In summary, under the same dataset, our model outperformed other competitive methods.
Case studies
Case studies were conducted to further verify the capability of HFHLMDA to predict miRNAdisease associations. We implemented three different kinds of case studies in this study. In the first case study, we conducted HFHLMDA to predict potential diseasemiRNA associations taking advantages of known diseasesmiRNAs associations included in HMDD v2.0 database. Subsequently, top 50 miRNAs for the investigated disease ranked according to their predicted scores were verified using another two wellknown miRNAdisease association databases of dbDEMC [42] and miR2Disease [43]. In the second case study, we simulated the situation where HFHLMDA was conducted for disease without known miRNA associations. More concretely, we removed the known miRNA associations of the disease of interest, after which HFHLMDA was implemented according newly obtained association records. The prediction results were also verified by other databases. The final case study investigated the robustness of HFHLMDA prediction performance. We evaluated the model with a smaller and earlier version HMDDv1.0 database [44].
Esophageal cancer (EC) is one of the most common cancers worldwide, and its 5year survival rate is about 20% [45]. Study indicate that miR130b plays an oncogenic role in esophageal squamous cell carcinoma cells by repressing phosphatase and tensin homolog expression and Akt phosphorylation [46]. Therefore, specific and sensitive biomarkers for diagnosis and targeted therapy of EC are urgently needed. As the first type of case study, 10 out of top 10, 28 out of top 30, 45 out of top 50 predicted esophageal neoplasms related miRNAs were confirmed by dbDEMC (See Table 1).
Hepatocellular carcinoma (HC) is a complex polygenetic disease ascribed to the interactions between genetic predisposition and environmental factors [47]. The discovery of vital target for genetic therapy are of great clinical significance to the improvement of the comprehensive effect of HC. For example, miR122, let7 family, and miR101 are downregulated in HC, suggesting that it is a potential tumor suppressor of HC. miR221 and miR222 are upregulated in HCC and may act as oncogenic miRNAs in hepatocarcinogenesis [48]. We took hepatocellular carcinoma as the second kind of case study. Finally, 49 out of top 50 miRNAs were experimentally confirmed by HMDD v2.0, dbDEMC and miR2Disease (See Table 2).
Breast Neoplasms is the most common malignancy in women, accounting more than 40,000 deaths each year [49]. Data have shown that the number of affected people is climbing, and a forecast deemed that there will be nearly 3.2 million new patients per year by 2050 [50]. In breast cancer, approximately onefifth of metastatic patients survive 5 years [51]. Researchers have found that many miRNAs are associated with breast neoplasms by clinical experiments, such as mir‐155 and mir‐21, both of which can lead to Breast Neoplasms tumorigenesis or metastasis [52]. We took breast neoplasms as the last kind of case study, in which we got the prediction with HFHLMDA using HMDDv1.0 database. Then, we verified the predicted potential breast neoplasms related miRNAs in other databases. At last, 48 out of top 50 miRNAs were experimentally confirmed by HMDD v2.0, dbDEMC and miR2Disease (See Table 3).
The aforementioned case studies indicate that HFHLMDA has good prediction performance. HFHLMDA can efficiently predict diseaserelated miRNAs based on known miRNAdisease associations, disease semantic similarity and miRNA functional similarity, and a disease without known associations also can be predicted.
Discussion
In this work, we developed a new computational model based on hypergraph learning to predict potential miRNA‐disease associations. Several important factors contribute to the excellent performance of our model. First, highdimensionality features. Based on a credible assumption that functionally similar miRNAs tend to have associations with phenotypically similar diseases. We use the miRNAs or diseases similarity scores directly as a feature factor, with a dimension of up to 878, which contains all similar information about miRNAs or diseases. Second, hypergraph is suitable to represent local group information and the highorder relationship of data, and can completely represent the complex relationships among miRNAdisease pairs. Different from the simplegraph learning methods consider only the pairwise relationship between two samples, and they ignore the relationship in a higherorder, hypergraph learning aims to get the relationship between several samples in a higher order. Hypergraph learning is a kind of graph clustering algorithm, the process of graph clustering is actually the optimization of graph partition. The purpose of optimization is to reduce the similarity between subgraphs and increase the similarity within subgraphs. Hypergraphbased models have proven to be beneficial for a variety of classification/clustering tasks, and we think it can also be applied to different fields of bioinformatics, such as drugdisease associations [53], miRNA–drug interactions [54].
Despite the practicability and efficiency of HFHLMDA, there still has some limitations. Since our method is based on machine learning techniques, negative samples are required during the training process. However, experimentally confirmed negative samples are difficult to obtain. To resolve this issue, we have randomly selected a subset of unknown miRNA–disease associations as negative instances. In addition, in our method, after the hypergraph has been constructed, it never changes during the learning process, leading to a static hypergraph structure learning mechanism. However, it is uneasy to guarantee that the generated hypergraph structure is optimal and suitable for all applications. In future work, it is necessary to investigate the hypergraph structure optimization, leading to a dynamic hypergraph structure learning scheme.
Conclusion
Increasing evidence indicates that aberrant expression of miRNAs is closely related to the occurrence and development of human complex diseases. Understanding the underlying mechanisms of miRNAs in diseases is becoming an urgent problem worldwide. Compared with traditional methods, the computational model developed for processing heterogeneous biological big data is more efficient and convenient. To predict potentially diseaserelated miRNAs, we proposed a hypergraph learning method called HFHLMDA. Both crossvalidation and case studies had proved the effectiveness of HFHLMDA in predicting potential miRNAdisease associations.
Availability of data and materials
The datasets used during this study is provided by Li et al. [35]. Please download the data from http://www.cuilab.cn/hmdd/ or contact the authors for data requests.
Abbreviations
 KNN:

KNearestNeighbor
 HMDD:

Human microRNA disease database
 dbDEMC:

Database of differentially expressed miRNAs in human cancers
 ROC:

Receiver operating characteristics
 AUC:

The area under the ROC curve
 LOOCV:

Leaveoneout cross validation
References
 1.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.
 2.
Kye MJ, Gonçalves ICG. The role of miRNA in motor neuron disease. Front Cell Neurosci. 2014;8:15.
 3.
Adams BD, Kasinski AL, Slack FJ. Aberrant regulation and function of microRNAs in cancer. Curr Biol. 2014;24(16):R762–76.
 4.
Cheng AM, Byrom MW, Shelton J, Ford LP. Antisense inhibition of human miRNAs and indications for an involvement of miRNA in cell growth and apoptosis. Nucleic Acids Res. 2005;33(4):1290–7.
 5.
Karp X, Ambros V. Encountering microRNAs in cell fate signaling. Science. 2005;310(5752):1288–9.
 6.
Shivdasani RA. MicroRNAs: regulators of gene expression and cell differentiation. Blood. 2006;108(12):3646–53.
 7.
Sayed D, Abdellatif M. MicroRNAs in development and disease. Physiol Rev. 2011;91(3):827–87.
 8.
Tricoli JV, Jacobson JW. MicroRNA: potential for cancer detection, diagnosis, and prognosis. Cancer Res. 2007;67(10):4553–5.
 9.
Cho WCS. MicroRNAs: potential biomarkers for cancer diagnosis, prognosis and targets for therapy. Int J Biochem Cell Biol. 2010;42(8):1273–81.
 10.
Pan LP, Liu F, Zhang JL, et al. Genomewide miRNA analysis identifies potential biomarkers in distinguishing tuberculous and viral meningitis. Front Cell Infect Microbiol. 2019;9:323.
 11.
Jiang QH, Hao YY, Wang GH, Juan LR, Zhang TJ, et al. Prioritization of disease microRNAs through a human phenomemicroRNAome network. BMC Syst Biol. 2010;4(Suppl 1):S2.
 12.
Chen X, Liu MX, Yan GY. RWRMDA: predicting novel human microRNAdisease associations. Mol BioSyst. 2012;8(10):2792–8.
 13.
Shi HB, Xu J, Zhang GD, Xu LD, Li CQ, Wang L, et al. Walking the interactome to identify human miRNAdisease associations through the functional link between miRNA targets and disease genes. BMC Syst Biol. 2013;7(1):101.
 14.
Zhao XM, Liu KQ, Zhu G, et al. Identifying cancerrelated microRNAs based on gene expression data. Bioinformatics. 2015;31(8):1226–34.
 15.
Qin GM, Li RY, Zhao XM. Identifying disease associated miRNAs based on protein domains. IEEE/ACM Trans Comput Biol Bioinform. 2016;13(6):1027–35.
 16.
Chen X, Huang L. LRSSLMDA: Laplacian regularized sparse subspace learning for MiRNAdisease association prediction. PLoS Comput Biol. 2017;13(12):e1005912.
 17.
Chen X, Xie D, Wang L, Zhao Q, You ZH, Liu H. BNPMDA: bipartite network projection for MiRNAdisease association prediction. Bioinformatics. 2018;34(18):3178–86.
 18.
Chen X, Wang L, Qu J, Guan NN, Li JQ. Predicting miRNAdisease association based on inductive matrix completion. Bioinformatics. 2018;34(24):4256–65.
 19.
Chen X, Yin J, Qu J, Huang L. MDHGI: matrix decomposition and heterogeneous graph inference for miRNAdisease association prediction. PLoS Comput Biol. 2018;14(8):e1006418.
 20.
Chen X, Zhu CC, Yin J. Ensemble of decision tree reveals potential miRNAdisease associations. PLoS Comput Biol. 2019;15(7):e1007209.
 21.
Chen X, Xie D, Zhao Q, You ZH. MicroRNAs and complex diseases: from experimental results to computational models. Brief Bioinform. 2019;20(2):515–39.
 22.
Xiao Q, Luo JW, Liang C, Cai J, Ding PJ. A graph regularized nonnegative matrix factorization method for identifying microRNAdisease associations. Bioinformatics. 2018;34(2):239–48.
 23.
Liu Y, Zeng X, He Z, Zou Q. Inferring microRNAdisease associations by random walk on a heterogeneous network with multiple data sources. IEEE/ACM Trans Comput Boil Bioinform. 2017;14(4):905–15.
 24.
You ZH, Huang ZA, Zhu Z, et al. PBMDA: a novel and effective pathbased computational model for miRNAdisease association prediction. PLoS Comput Biol. 2017;13(3):e1005455.
 25.
Chen X, Niu YW, Wang GH, et al. HAMDA: hybrid approach for MiRNAdisease association prediction. J Biomed Inform. 2017;76:50–8.
 26.
Chen X, Guan NN, Li JQ, et al. GIMDA: graphlet interactionbased MiRNAdisease association prediction. J Cell Mol Med. 2018;22(3):1548–61.
 27.
Chen X, Yang JR, Guan NN, Li JQ. GRMDA: graph regression for MiRNAdisease association prediction. Front Physiol. 2018;9:92.
 28.
Jiang YD, Liu BT, Yu LH, et al. Predict MiRNAdisease association with collaborative filtering. Neuroinformatics. 2018;16:363–72.
 29.
Xu J, Li CX, Lv JY, et al. Prioritizing candidate disease miRNAs by topological features in the miRNA targetdysregulated network: case study of prostate cancer. Mol Cancer Ther. 2011;10(10):1857–66.
 30.
Chen X, Yan GY. Semisupervised learning for potential human microRNAdisease associations inference. Sci Rep. 2014;4:5501.
 31.
Luo JW, Xiao Q, Liang C, Ding PJ. Predicting microRNAdisease associations using Kronecker regularized least squares based on heterogeneous omics data. IEEE Access. 2017;5:2503–13.
 32.
Chen X, Yan CC, Zhang X, Li Z, Deng L, et al. RBMMMDA: predicting multiple types of diseasemicroRNA associations. Sci Rep. 2015;5:13877.
 33.
Chen X, Yan CC, Zhang X, You ZH, et al. HGIMDA: heterogeneous graph inference for miRNAdisease association prediction. Oncotarget. 2016;7(40):65257–69.
 34.
Chen X, Huang L, Xie D, Zhao Q. EGBMMDA: extreme gradient boosting machine for MiRNAdisease association prediction. Cell Death Dis. 2018;9(1):3.
 35.
Li Y, Qiu CX, Tu J, Geng B, Yang JC, et al. HMDD v2.0: a database for experimentally supported human microRNA and disease associations. Nucleic Acids Res. 2014;42:D1070–4.
 36.
Wang D, Wang J, Lu M, Song F, Cui QH. Inferring the human microRNA functional similarity and functional network based on microRNAassociated diseases. Bioinformatics. 2010;26(13):1644–50.
 37.
Van Laarhoven T, Nabuurs SB, Marchiori E. Gaussian interaction profile kernels for predicting drugtarget interaction. Bioinformatics. 2011;27(21):3036–43.
 38.
Wang L, You ZH, Huang YA, Huang DS, Chan KCC. An efficient approach based on multisources information to predict CircRNAdisease associations using deep convoltional neural network. Bioinformatics. 2020;36(13):4038–46.
 39.
Zhou DY, Huang JY, Schlkopf B. Learning with hypergraphs: clustering, classification, and embedding. Adv Neural Inf Process Syst. 2006;19:1601–8.
 40.
Zhang ZZ, Liu HJ, Zhao XB, Ji RR, Gao Y. Inductive multihypergraph learning and its application on viewbased 3D object classification. IEEE Trans Image Process. 2018;27(12):5957–68.
 41.
Shao BY, Liu BT, Yan CG. SACMDA: MiRNAdisease association prediction with short acyclic connections in heterogeneous graph. Neuroinformatics. 2018;16:373–82.
 42.
Yang Z, Ren F, Liu CN, He SM, Sun G, Gao Q, et al. dbDEMC: a database of differentially expressed miRNAs in human cancers. BMC Genomics. 2010;11(Suppl 4):S5.
 43.
Jiang QH, Wang YD, Hao YY, Juan LR, Teng MX, et al. miR2Disease: a manually curated database for microRNA deregulation in human disease. Nucleic Acids Res. 2009;37(1):D98–104.
 44.
Lu M, Zhang QP, Deng M, Miao J, Guo YH, et al. An analysis of human microRNA and disease associations. PLoS ONE. 2008;3(10):e3420.
 45.
Torre LA, Bray F, Siegel RL, Ferlay J, LortetTieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65:87–108.
 46.
Yu T, Cao R, Li S, et al. MiR130b plays an oncogenic role by repressing PTEN expression in esophageal squamous cell carcinoma cells. BMC Cancer. 2015;15:29.
 47.
Nie X, Liu Y, Chen WD, Wang YD. Interplay of miRNAs and canonical Wnt signaling pathway in hepatocellular carcinoma. Front Pharmacol. 2018;9:657.
 48.
Saito Y, Suzuki H, Matsuura M, Sato A, Kasai Y, et al. MicroRNAs in hepatobiliary and pancreatic cancers. Front Gene. 2011;2:66.
 49.
Desantis CE, Fedewa SA, et al. Breast cancer statistics, 2015: convergence of incidence ratesbetween black and white women. CA Cancer J Clin. 2016;66(1):31–42.
 50.
Gomella LG. Prostate cancer statistics: anything you want them to be. Can J Urol. 2017;24(1):8603–4.
 51.
Lee JH, Zhao XM, Yoon I, et al. Integrative analysis of mutational and transcriptional profiles reveals driver mutations of metastatic breast cancers. Cell Discov. 2016;2:16025.
 52.
Feber A, Xi L, Luketich JD, et al. MicroRNA expression profiles of esophageal cancer. J Thorac Cardiovasc Surg. 2008;135(2):255–60.
 53.
Yang K, Zhao X, Waxman D, Zhao XM. Predicting drugdisease associations with heterogeneous network embedding. Chaos. 2019;29(12):123109.
 54.
Xie WB, Yan H, Zhao XM. EmDL: extracting miRNA–drug interactions from literature. IEEE/ACM Trans Comput Biol Bioinform. 2019;16(5):1722–8.
Acknowledgments
We would like to thank the editor and referees for the thoughtful and insightful comments.
About this supplement
This article has been published as part of BMC Medical Informatics and Decision Making Volume 21 Supplement 1, 2021: Proceedings of the 2019 International Conference on Intelligent Computing (ICIC 2019): medical informatics and decision making. The full contents of the supplement are available online at https://bmcmedinformdecismak.biomedcentral.com/articles/supplements/volume21supplement1.
Funding
Publication costs are funded by the National Natural Science Foundation of China (Nos. 61873001, U19A2064, 61872220, 61672037, 61861146002, 11701318 and 61732012), the Key Project of Anhui Provincial Education Department (No. KJ2017ZD01), and the Xinjiang Autonomous Region University Research Program (XJEDU2019Y002). The funding bodies did not play any role in the design of the study or collection, analysis and interpretation of data or in writing the manuscript.
Author information
Affiliations
Contributions
JCN and CHZ supervised the entire project. YTW and QWW conceptualized and designed the study. YTW and ZG undertook data collection. YTW, QWW and ZG performed the data analysis. QWW drafted the initial version. JCN and CHZ revised the manuscript iteratively for important intellectual content. All authors edited the paper and gave final approval for the version to be published. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wang, YT., Wu, QW., Gao, Z. et al. MiRNAdisease association prediction via hypergraph learning based on highdimensionality features. BMC Med Inform Decis Mak 21, 133 (2021). https://doi.org/10.1186/s1291102001320w
Received:
Accepted:
Published:
Keywords
 MicroRNA
 Disease
 MiRNAdisease association
 Knearestneighbor
 Hypergraph learning