Skip to main content

Identifying potential association on gene-disease network via dual hypergraph regularized least squares

Abstract

Background

Identifying potential associations between genes and diseases via biomedical experiments must be the time-consuming and expensive research works. The computational technologies based on machine learning models have been widely utilized to explore genetic information related to complex diseases. Importantly, the gene-disease association detection can be defined as the link prediction problem in bipartite network. However, many existing methods do not utilize multiple sources of biological information; Additionally, they do not extract higher-order relationships among genes and diseases.

Results

In this study, we propose a novel method called Dual Hypergraph Regularized Least Squares (DHRLS) with Centered Kernel Alignment-based Multiple Kernel Learning (CKA-MKL), in order to detect all potential gene-disease associations. First, we construct multiple kernels based on various biological data sources in gene and disease spaces respectively. After that, we use CAK-MKL to obtain the optimal kernels in the two spaces respectively. To specific, hypergraph can be employed to establish higher-order relationships. Finally, our DHRLS model is solved by the Alternating Least squares algorithm (ALSA), for predicting gene-disease associations.

Conclusion

Comparing with many outstanding prediction tools, DHRLS achieves best performance on gene-disease associations network under two types of cross validation. To verify robustness, our proposed approach has excellent prediction performance on six real-world networks. Our research work can effectively discover potential disease-associated genes and provide guidance for the follow-up verification methods of complex diseases.

Background

Identification of the association between disease and human gene has attracted more attention in the field of biomedicine, and has become an important research topic. A great deal of evidence shows that understanding genes related to diseases is of great help to prevent and treat diseases. However, identifying the relationship between disease and gene by biological experiments has to spend a long time and cost. Many computational models have been proposed to solve some similar biologically related problems. For example, in the fields of biology [13], pharmacy [4], and medicine [5, 6], machine learning methods help solve many analytical tasks.

In order to explore the relationship between gene and disease, a variety of algorithms have been proposed for association prediction. The typical machine learning methods [710] is to extract relevant features of known genetic data of each disease and train the model to determine which disease is related to those genes, so these algorithms are usually single-task algorithms for each disease. This model needs to be trained separately. Therefore, for a new disease or an existing disease with few known genes, due to the lack of known association data or the relevant information between various diseases, it is difficult to train the learning model. As a machine learning method, the matrix completion methods [1113] can solve the above problem by calculating the similarity information and predicting the association between disease and gene, but the matrix completion method usually takes a long time to converge the local optimal solution. The other type is network-based model [1417]. Li et al. [17] predicted the association by systematically embedding a heterogeneous network of genes and diseases into Graph Convolutional Network. This model usually divides genes and diseases into two heterogeneous networks. The edges in network represent the similarity between nodes. The model is based on the assumption that genes with high similarity are easily related to similar diseases. However, they are biased by the network topology, and it is necessary to rely on effective similarity information. It is not easy for these methods to integrate related sources of multiple genes and diseases.

Multiple Kernel Learning (MKL) is an important machine learning method, which can effectively combine multi-source information to improve the model effect, and is applied to many biological problems. For instance, Yu et al. [8] implemented one-class of Support Vector Machine while optimizing the linear combination of the gene nucleus and the MKL method. Ding et al. [1821] proposed multiple information fusion models to identify drug-target and drug-side effect associations. Wang et al. [22] proposed a novel Multiple Kernel Support Vector Machine (MKSVM) classifier based on Hilbert Schmidt Independence Criterion to identify membrane proteins. Shen [23] and Ding et al. [24] proposed a MKSVM model to identify multi-label protein subcellular localization. Ding et al also employ fuzzy-besd model to predict DNA-binding proteins [25] and protein crystallization [26]. Zhang et al. [27] developed an ensemble predictive model of classifier chain to identify anti-inflammatory peptides.

LapRLS framework [28] is often used in various fields based on machine learning model, such as the prediction of Human Microbe-Disease Association [29] and the detection of human microRNA-disease association [30]. At the same time, Hypergraph learning [3133] is becoming popular. Hypergraphs can represent more complex relationships among various objects. Bai et al. [34] introduced two end-to-end trainable operators to the family of graph neural networks, i.e., hypergraph convolution and hypergraph attention. Whilst hypergraph convolution defines the basic formulation of performing convolution on a hypergraph, hypergraph attention further enhances the capacity of representation learning by leveraging an attention module. Zhang et al. [35] developed a new self-attention based graph neural network called Hyper-SAGNN applicable to homogeneous and heterogeneous hypergraphs with variable hyperedge sizes. Ding et al. [36] predicted miRNAs-disease associations by a hypergraph regularized bipartite local model, which is based on hypergraph embedded Laplacian support vector machine.

Inspired by what is mentioned above, we propose a novel prediction method named Dual Graph Hypergraph Least Squares model (DHRLS) to predict gene-disease associations. Some computational models based on graph learning can effectively solve various network problems. In this paper, the gene-disease association detection can be defined as the link prediction problem in bipartite network [3739]. Furthermore, two feature spaces are described by similarity information of multiple genes and diseases. Multiple kernel learning is also used to combine multiple informations linearly. Here, we use the Centered Kernel Alignment-based Multiple Kernel Learning (CKA-MKL) [40] to obtain weights of multiple kernels and then combine these kernels via optimal weights in two spaces, respectively. In addition, we also embed hypergraphs in graph regular terms to preserve high-order information of genes and diseases, using more complex information to improve prediction performance. To prove the effectiveness of our proposed method, six types of real networks and one gene-disease associations network are employed to test our predictive model. On the gene-disease associations dataset, our method has been compared with some methods under two types of cross-validation (CV). Comparing DHRLS with other state-of-the-art methods on predicting gene-disease associations, including CMF, GRMF and Spa-LapRLS, our model achieves the highest AUC and AUPR in 10-fold cross validation under CV1, but our model achieves lower AUC under CV2 compared with Spa-LapRLS. At the same time, DHRLS has excellent prediction performance on six benchmark datasets.

Results

In order to better test the performance of our method, our proposed approach is verified on real gene-disease associations dataset under two types of cross validation. We also test the capability of DHRLS in predicting novel disease after confirming the excellent performance of our method based on cross validation. Furthermore, we employ benchmark datasets to evaluate our approach and compare it with other existing methods.

Dataset

We download the dataset of gene-disease associations from [41] (http://cssb2.biology.gatech.edu/knowgene). Since the number of genes is too large and the information is insufficient, we remove some redundant gene data. Finally, our dataset contains 31,519 associations of 960 diseases and 3,118 genes, where 279 genes are associated with only one disease.

Evaluation measurements

The 10-fold Cross Validation (CV) is usually used to verify the bipartite network detection. In order to compare the prediction performance with other methods under the same evaluation measurement, we will also use 10-fold CV for verification. At the same time, Area under the receiver operating characteristic curve (AUC) and Area Under the Precision-Recall curve (AUPR) as the major evaluation indicator, will also be applied to evaluate methods. There are two CV settings as follows:

CV1: Pair prediction. All gene-disease associations are randomly divided into test set and training set, and the associations in the test set are removed.

CV2: Disease prediction. All diseases are randomly divided into test set and training set, and all associations of diseases in the test set are removed.

Parameter settings

In our study, DHRLS has some parameters λd,λg,β, k and number of iterations. In the parameter selection, we consider all combinations of following values: number of k-Nearest Neighbor is from 10 to 100 (with step 10); number of iterations is {1,2,...,15}; {2−5,...,20,...,25} for λd and λg; β=1.

Figure 1 shows the results of our model obtained under different iteration times and k values. For the number of k-Nearest Neighbor, we select the optimal k under the highest AUPR value and can clearly find that AURP reaches its peak when k=50. For the number of iterations, it basically converges at the four times. In order to train the model more fully, we finally choose the number of iterations to be 10.

Fig. 1
figure 1

Different AUPR values under different k-Nearest Neighbor and iterations

Figure 2 shows the results of AUC and AUPR in grid search for parameters λd and λg. The optimal λd and λs are also selected under highest AUPR value. In this study, the optimal parameters of Hypergraph Laplace regular terms are obtained on λd=1 and λg=0.25. Under this parameter selection, the AUC value is relatively high.

Fig. 2
figure 2

The AUC (a) and AUPR (b) of models with different λd and λg under CV1. λd (horizontal axis) and λg (vertical axis) are set from 2−5 to 25 with step 21. The yellow color is the higher value, and blue color is the lower value

Evaluation on gene-disease association data

Performance analysis

We evaluate the different performance of CKA-MKL, mean weighted-based MKL and single kernel \(\left (\mathbf {K}^{d}_{SEM} \& \mathbf {K}^{g}_{GO} \text {and} \mathbf {K}^{d}_{GIP} \& \mathbf {K}^{g}_{GIP}\right)\). The testing results are shown in Table 1 and Fig. 3.

Fig. 3
figure 3

The performance of different models under CV1

Table 1 The performance of different models under CV1

Obviously, the model of CKA-MKL on DHRLS obtains the best performance with AUC of 0.9742 and AUPR of 0.8092. Comparing with mean weighted on DHRLS, AUPR and AUC are increased by 0.0086 and 0.0039. This means that CKA combines multi-kernel information more effectively than simple average combination. What’s more, DHRLS with single kernel \(\left (\mathbf {K}^{d}_{SEM} \& \mathbf {K}^{g}_{GO}\right)\) obtains lower performance than the model with GIP kernel. Therefore, GIP is an effective method to calculate the kernel matrix. By comparing the results of single kernel and multi-kernel models, combining multiple information is an effective method to improve the prediction effect of the model.

Furthermore, Fig. 4 shows the weights of each kernel matrix in the gene space and disease space. The weight of the kernel indicates the degree of contribution of the corresponding kernel matrix. Comparing the weights in the gene and disease spaces, the GIP kernel has a higher weight in both spaces, which is consistent with the results in Table 1. In gene space, except for GIP kernel, the kernel weight of \(\mathbf {K}^{g}_{GO}\) is higher than \(\mathbf {K}^{g}_{PPI}\) and \(\mathbf {K}^{g}_{SW}\). This means that \(\mathbf {K}^{g}_{GO}\)’s contribution to the overall is better than the other two kernel matrices.

Fig. 4
figure 4

The kernel weights of drug and gene space respectively

Comparison to existing predictors

Many excellent methods have been proposed to predict the bipartite network link, including Spa-LapRLS [30], GRMF [42] and CMF [43]. Our method is compared to the existing methods and DGRLS under CV1 and CV2, respectively. Under CV1, the results are shown in Table 2 and Fig. 5. Our method achieves the best AUC (0.9742) and AUPR (0.8092). For AUC, DHRLS is not much different from DGRLS and Spa-LapRLS, which is about 0.01 higher than GRMF and CMF. As for AUPR, DHRLS achieves better performance than other methods. Comparing the results of DHRLS and DGRLS, it can be seen that the hypergraph-based model is better than the normal graph model, which shows that the high-level graph information constructed by the hypergraph is helpful for the predict performance. This is related to the ability of hypergraph to effectively find similar information between nodes. At the same time, the methods based on LapRLS (DHRLS, DGRLS and Spa-LapRLS) are higher than those based on matrix factorization (GRMF and CMF), indicating that the model framework of LapRLS has more advantages in the prediction of gene-disease associations.

Fig. 5
figure 5

The performance of different methods under CV1

Table 2 The performance of different methods under CV1

In order to test the performance of our method detecting new diseases, the associations for new diseases (CV2) are not observed in the training set. Table 3 and Fig. 6 show the results of CV2. Under CV2, our method obtains best AUPR (0.1413). However, the performance of our model on AUC (0.8987) is secondary best, which is about 0.02 lower than that of Spa-LapRLS.

Fig. 6
figure 6

The performance of different methods under CV2

Table 3 The performance of different methods under CV2

Comparing the results of DGRLS and DHRLS under CV1 and CV2, we clearly find that utilizing hypergraph to establish higher-order relationships greatly improves the predictive ability of the model.

Case study

Our model can predict genes associated with new diseases. Here, we use DHRLS to rank the predicted values of genes related to new diseases in descending order. The higher the ranking, the more likely it is to interact. We set the value of a disease in the correlation matrix to 0 as a new disease. One example is Lung Diseases. We intercepted the top 50 predicted genes and 40 (80%) known related genes in the predicted results. All predicted ranking results are shown in Table 4.

Table 4 Predicted top 50 genes for Lung Diseases by our method

Evaluation on six benchmark datasets

To test the performance of our proposed method, we consider six real-world networks: (i) G-protein coupled receptors (GPC Receptors): the biological network of drugs binding GPC receptors; (ii)Ion channels: the biological network of drugs binding ion channel proteins; (iii) Enzymes: the biological network of drugs binding enzyme proteins; (iv) Southern Women (referred here as “SW”): the social relations network of women and events; (v) Drug-target: the chemical network of drug-target interaction; (vi) Country-organization (referred here as “CO”): the network of organization most related to the country. Detailed information about six datasets is described in Table 5.

Table 5 Statistics of six real-world networks

Since there is only the data of interaction matrix of binary network, in order not to introduce additional data, we directly use the GIP kernel extracted from the interaction matrix as the kernel matrix for each real-world network. The kernel is defined as follows:

$$\begin{array}{*{20}l} \mathbf{K}^{*}_{1}(i,j) &= \exp \left(-\gamma_{1}||\mathbf{Y}_{i}-\mathbf{Y}_{j}||^{2}\right) \end{array} $$
(1a)
$$\begin{array}{*{20}l} \mathbf{K}^{*}_{2}(m,n) &= \exp \left(-\gamma_{2}||\mathbf{Y}^{T}_{m}-\mathbf{Y}^{T}_{n}||^{2}\right) \end{array} $$
(1b)

where Y is the train set of binary network, and Yi is the vector of associations.

We test our method on above six datasets and compare results with other methods [44]. Wang et al. [44] proposed a framework, called Similarity Regularized Nonnegative Matrix Factorization (SRNMF), for link prediction in bipartite networks by combining the similarity based structure and the latent feature model from a new perspective. Tables 6 and 7 show the comparison of precision and AUC for six real-world networks. DHRLS performs better than other methods on Enzymes and Ionchannel networks, and values of our precision and AUC are higher than others. For GPC and Drug-target networks, the precision is same, but AUC is slightly higher. This directly indicates the clear performance advantage of our approach in real-world binary networks.

Table 6 Precision by different methods on six real networks
Table 7 AUC by different methods on six real networks

Discussion

We developed the model DHRLS for the gene-disease association prediction. In order to evaluate our model, we test not only on real gene-disease associations dataset, but also on some benchmark datasets. By comparing the results of single-kernel model and multi-kernel model, MKL can effectively combine multi-kernel information to improve the predictive ability of the model. By adjusting different kernel weights, different kernel matrices can express different levels of information. However, MKL needs to be applied to samples with multiple feature information, and the application effect is not obvious for problems with fewer features. The comparison of DHRLS and DGRLS can illustrate the effectiveness of hypergraph. After adding the hypergraph, the result of the model is obviously improved, which is caused by the characteristics of the hypergraph. Hypergraph uses high-order information between nodes, that is, a hyperedge can connect more than two nodes, which can better indicate the degree of similarity between nodes. Comparing DHRLS with other state-of-the-art methods on predicting gene-disease associations, including CMF, GRMF and Spa-LapRLS, our model achieves the highest AUC and AUPR in 10-fold cross validation under CV1, but our model achieves lower AUC under CV2 compared with Spa-LapRLS. At the same time, DHRLS has excellent prediction performance on six benchmark datasets.

Nevertheless, our model still has some flaws.First of all, the model contains a large number of matrix operations and optimization problems, and lacks a certain degree of simplicity. Secondly, we need to calculate the multi-kernel information of the sample. Therefore, we cannot achieve predictions for samples without features. At present, most of the computational methods are developed to predict the associations of gene-disease, and there is still a great possibility to improve the prediction performance. For example, hypergraph can be considered in the graph based method. In the future, for optimizing the model and improving the prediction performance, we can add some data preprocessing and calculate simplification on the basis of DHRLS, as well as better method to build hypergraph.

Conclusion

In summary, we propose a Dual Hypergraph Regularized Least Squares (DHRLS) based on CKA-MKL algorithm, for the gene-disease association prediction. We use multiple kernels to describe gene and disease spaces. The weights of these kernels are obtained by CKA-MKL and used to combine kernels. We use hypergraph to describe more complex information to improve our prediction. Our purpose is to establish an accurate and effective prediction model of gene-disease association based on the existing data of gene-disease associations, and provide guidance for the follow-up verification methods of complex diseases.

Methods

In this study, we first use two disease kernels and four gene kernels to reveal potential associations of genes and diseases. Then, the MKL method CKA is used to combine above kernels into one disease kernel and one gene kernel. Finally, we use Dual Hypergraph Regularized Least Squares to identify gene-disease associations. Figure 7 show the flowchart of our method DHRLS.

Fig. 7
figure 7

The overview of our proposed method

Problem definition

The prediction of gene-disease associations can be regarded as a recommendation system. Given n diseases D={d1,d2,,...,dn}, m genes S={g1,g2,,...,gm} and gene-disease associations. The association between gene and disease items can be expressed as an adjacent matrix YRn×m. The element of adjacent matrix Y is the relationship between genes and diseases. If disease dj(1≤jm) is associated with gene gi(1≤in), the value of Yi,j is set as 1, otherwise it is 0. Genes, diseases, and their associations are formulated as a bipartite network.

Related work

LapRLS framework [28] is often used in various fields based on machine learning model, such as the prediction of Human Microbe-Disease Association [29] and the detection of human microRNA-disease association [30]. At the same time, Hypergraph learning [3133] is becoming popular. Hypergraphs can represent more complex relationships among various objects.

The Laplacian Regularized Least Squares (LapRLS) model [45] based on graph regularization is employed to predict potential associations in a bipartite network. The functions of model can be defined as follows:

$$ \begin{aligned} &{}\mathbf{F}^{*}_{a}=\arg\min{J}(\mathbf{F}_{a}) = ||\mathbf{Y}_{train}-\mathbf{F}_{a}||^{2}_{F}+\lambda_{a}tr\left(\mathbf{F}^{T}_{a}\mathbf{L}_{a}\mathbf{F}_{a}\right)\\ &{}\mathbf{F}^{*}_{b}=\arg \min {J}(\mathbf{F}_{b}) = ||\mathbf{Y}_{train}-\mathbf{F}_{b}||^{2}_{F}+\lambda_{b}tr\left(\mathbf{F}^{T}_{b}\mathbf{L}_{b}\mathbf{F}_{b}\right) \end{aligned} $$
(2)

where \(\mathbf {F}^{*}_{a}=\mathbf {K}_{a}\pmb {\alpha }^{*}_{a}, \mathbf {F}^{*}_{b}=\mathbf {K}_{b}\pmb {\alpha }^{*}_{b}\), and \(\mathbf {F}_{a}, \pmb {\alpha }^{*}_{a}, \mathbf {F}^{T}_{b}, {\pmb {\alpha }^{*}_{b}}^{T}, \mathbf {Y}_{train} \in \mathbf {R}^{n\times m}\). KaRn×n and KbRm×m are kernels in two feature space, separately.

LaRn×n and LbRm×m are the normalized Laplacian matrices as follows:

$$ \begin{aligned} \mathbf{L}_{a}=\mathbf{D}^{-1/2}_{a}\Delta_{a}\mathbf{D}^{1/2}_{a},&\Delta_{a}=\mathbf{D}_{a}-\mathbf{K}_{a}\\ \mathbf{L}_{b}=\mathbf{D}^{-1/2}_{b}\Delta_{b}\mathbf{D}^{1/2}_{b},&\Delta_{b}=\mathbf{D}_{b}-\mathbf{K}_{b} \end{aligned} $$
(3)

where Da and Db are diagonal matrices, \(\mathbf {D}_{a}(k,k)=\sum _{l=1}^{n}\mathbf {K}_{a}(k,l), \mathbf {D}_{b}(k,k)=\sum _{l=1}^{m}\mathbf {K}_{b}(k,l)\)

The variables αa and \(\pmb {\alpha }^{*}_{b}\) of LapRLS can be solved as follows:

$$ \begin{aligned} \pmb{\alpha}^{*}_{a}&=\left(\mathbf{K}_{a}+\lambda_{a}\mathbf{L}_{a}\mathbf{K}_{a}\right)^{-1}\mathbf{Y}_{train}\\ \pmb{\alpha}^{*}_{b}&=(\mathbf{K}_{b}+\lambda_{b}\mathbf{L}_{b}\mathbf{K}_{b})^{-1}(\mathbf{Y}_{train})^{T} \end{aligned} $$
(4)

And \(\mathbf {F}^{*}_{a}\) and \(\mathbf {F}^{*}_{b}\) can be calculated as follows:

$$ \begin{aligned} \mathbf{F}^{*}_{a}&=\mathbf{K}_{a}(\mathbf{K}_{a}+\lambda_{a}\mathbf{L}_{a}\mathbf{K}_{a})^{-1}\mathbf{Y}_{train}\\ \mathbf{F}^{*}_{b}&=\mathbf{K}_{b}\left(\mathbf{K}_{b}+\lambda_{b}\mathbf{L}_{b}\mathbf{K}_{b}\right)^{-1}\left(\mathbf{Y}_{train}\right)^{T} \end{aligned} $$
(5)

The predictions from two feature spaces are combined into:

$$ \mathbf{F}^{*}=\frac{\mathbf{F}^{*}_{a}+(\mathbf{F}^{*}_{b})^{T}}{2} $$
(6)

Feature extraction

To improve effectiveness of detecting gene-disease associations, We use two and four types of similarity for disease and gene separately. In our work, we constructed the multiple kernels of diseases and genes to represent the feature sets. Table 8 summarizes whole kernels, including two feature spaces.

Table 8 Summary of kernels in two feature spaces

Disease space

We calculate two classes of disease kernels, including semantic similarity kernel and Gaussian Interaction Profile (GIP) kernel (for disease).

a) Semantic similarity The disease semantic similarity kernel is calculated by the relative positions in the MeSH [46] disease. Directed Acyclic Graph (DAG) [47] can describe disease di as a node. A disease di can be described as a node in DAG and denoted as \({DAG}_{d_{i}} = (d_{i}, T_{d_{i}}, E_{d_{i}}) \), where \(T_{d_{i}}\) is the set of all ancestor nodes of di including node di itself and \(E_{d_{i}}\) is the set of corresponding links. A semantic score of each disease \( t \in T_{d_{i}}\) can be calculated as follows:

$$ D_{d_{i}}(t) \left\{\begin{array}{ll} 1& \text{if~} t =d_{i} \\ \max \{\Delta \ast D_{d_{i}}(t^{'})|t^{'} \in children \ of \ t \}& \text{if~} t \not= d_{i} \end{array}\right. $$
(7)

where Δ is the semantic contribution factor, which is set to 0.5 in this paper.

Then, the semantic score of disease di can be calculated as follows:

$$ DV(d_{i})={\sum\nolimits}_{t \in T_{d_{i}}}D_{d_{i}}(t) $$
(8)

So, the disease semantic similarity kernel \(K^{d}_{SEM} \in \mathbf {R}^{n\times n}\) is calculated as follows:

$$ K^{d}_{SEM}\left(d_{i},d_{j}\right)= \frac{{\sum\nolimits}_{t \in T_{d_{i}} \bigcap T_{d_{j}}}\left(D_{d_{i}}(t)+D_{d_{j}}(t)\right)} {DV(d_{i})+DV(d_{j})} $$
(9)

b) GIP kernel similarity The similarity between diseases can also be calculated by GIP. Given two diseases di and dj(i,j=1,2,...,n), the GIP kernel can be calculated as follows:

$$ K^{d}_{GIP}\left(d_{i},d_{j}\right)= \exp \left(-\gamma_{d}||\mathbf{Y}_{d_{i}}-\mathbf{Y}_{d_{j}}||^{2}\right) $$
(10)

where \(\mathbf {Y}_{d_{i}}\) and \(\mathbf {Y}_{d_{j}}\) are the information of associations for vector disease di and dj. γd (set as 0.5) is the bandwidth of GIP kernel.

Gene space

Four types of gene kernels, including Gene Ontology (GO) [48] similarity, Protein-protein interactions (PPIs) network similarity, sequence similarity kernel and GIP kernel (for gene) are utilized to represent the relationship between genes.

a) GO similarity The information of GO is obtained through DAVID [49]. GO similarity \(\left (K^{g}_{GO} \in \mathbf {R}^{m \times m}\right)\) is the overlap of GO annotations on two genes, and we simply use GOSemSim [50] to get it. We consider one option of GO: cellular component (CC) to represent gene functional annotation.

b) PPIs similarity We download the protein-protein interactions network from previous research [41] and select the sub-networks related to our genes. Give the topological feature vectors pi and pj of two genes in the PPIs network. The Cosine similarity of PPIs network can be calculated as follows:

$$ K^{g}_{PPI}\left(p_{i},p_{j}\right)= \frac{p_{i} \cdot p_{j}}{||p_{i}||||p_{j}||} $$
(11)

c) Sequence similarity We use the normalized Smith Waterman (SW) score [51] to measure the sequence similarity between the two gene sequences, which is calculated as follows:

$$ K^{g}_{SW}\left(g_{i},g_{j}\right)= \frac{SW\left(S_{g_{i}},S_{g_{j}}\right)}{\sqrt{SW\left(S_{g_{i}},S_{g_{i}}\right)} \sqrt{SW\left(S_{g_{j}},S_{g_{j}}\right)}} $$
(12)

where SW(.,.) is Smith Waterman score. \(S_{g_{i}}\) is the information of sequence for gene gi.

d) GIP kernel similarity GIP is also employed to build gene GIP kernel \(\left (K^{g}_{GIP}\right)\). Given two genes gi and gj(i,j=1,2,...,m), the GIP kernel can be calculated as follows:

$$ K^{g}_{GIP}\left(g_{i},g_{j}\right)= \exp \left(-\gamma_{g}||\mathbf{Y}_{g_{i}}-\mathbf{Y}_{g_{j}}||^{2}\right) $$
(13)

where \(\mathbf {Y}_{g_{i}}\) and \(\mathbf {Y}_{g_{j}}\) are the information of associations for vector gene gi and gj. γg (set as 0.5) is the bandwidth of GIP kernel.

Multiple kernel learning

In our work, two kernels in the disease space including \(K^{d}_{SEM}\) and \(K^{d}_{GIP}\), and four kernels of gene space including \(K^{g}_{SW}, K^{g}_{GO}, K^{g}_{PPI}\) and \(K^{d}_{GIP}\). We then need to combine these kernels by means of linear combination in order to achieve the optimal ones.

$$ \begin{aligned} \mathbf{K^{*}}=\sum_{i=1}^{k}\omega_{i}\mathbf{K}_{i}\\ \mathbf{K}_{i}\in \mathbf{R}^{N\times N}\\ \sum_{i=1}^{k}\omega_{i} = 1 \end{aligned} $$
(14)

where k is the number of kernels and ωi is the weight of the kernel Ki. N is the number of samples in kernel Ki.

The method CKA-MKL is utilized to combine gene kernels and disease kernels, respectively. The cosine similarity between K1 and K2 is defined as follows:

$$ CA\left(\mathbf{K}_{1},\mathbf{K}_{2}\right)= \frac{<\mathbf{K}_{1},\mathbf{K}_{2}>_{F}}{||\mathbf{K}_{1}||_{F}||\mathbf{K}_{2}||_{F}} $$
(15)

where K1,K2Rn×n,<K1,K2>F=Trace(K1TK2) is the Frobenius inner product and \(||\mathbf {K}_{1}||_{F}=\sqrt {<\mathbf {K}_{1},\mathbf {K}_{1}>_{F}}\) is Frobenius norm.

The higher the cosine value, the greater the similarity between the kernels. CKA is based on the assumption that the combined kernel (feature space) should be similar to the ideal kernel (label space). Therefore, the alignment score between the combined kernel and the ideal kernel should be maximized. The objective function of centered kernel alignment is as follows:

$$ \begin{aligned} \max\limits_{\boldsymbol{\omega} \ge 0} CA\left(\mathbf{K}_{\boldsymbol{\omega}},\mathbf{K}_{ideal}\right) \\ =\max\limits_{\boldsymbol{\omega}\ge 0}\frac{<\mathbf{K}^{c}_{\boldsymbol{\omega}},\mathbf{K}_{ideal}>_{F}}{||\mathbf{K}^{c}_{\boldsymbol{\omega}}||_{F}||\mathbf{K}_{ideal}||_{F}}\\ \text{subject to} \mathbf{K}_{\boldsymbol{\omega}} = \sum_{i=1}^{k}\omega_{i}\mathbf{K}_{i}\\ \mathbf{K}^{c}_{\boldsymbol{\omega}} = \mathbf{U}_{N}\mathbf{K}_{\omega}\mathbf{U}_{N} \end{aligned} $$
(16)

where \(\sum _{i=1}^{k}\omega _{i}=1, \omega _{i} \ge 0, i=1,2,...,k\). \(\mathbf {U}_{N} = \mathbf {I}_{N}-(1/N)\mathbf {l}_{N}\mathbf {l}^{T}_{N}\) denotes a centering matrix, and INRN×N is the N-order identity matrix, lN is the N-order vector with all entries equal to one. \(\mathbf {K}^{c}_{\boldsymbol {\omega }}\) is the centered kernel matrix associated with Kω. Equation 16 can be written as follow:

$$ \min\limits_{\boldsymbol{\omega} \ge 0} \boldsymbol{\omega}^{T}\mathbf{M}\boldsymbol{\omega} - 2\boldsymbol{\omega}^{T}\mathbf{a} $$
(17)

where \(\mathbf {a}\,=\,\!\left (\!<\!\!\mathbf {K}^{c}_{1},\!\mathbf {K}_{ideal}\!\!>_{F}\!,\!<\!\mathbf {K}^{c}_{2},\!\mathbf {K}_{ideal}\!\!>\!\!_{F},...\!<\!\!\mathbf {K}^{c}_{k},\!\mathbf {K}_{ideal}\!\!>\!\!_{F}\!\right)\!^{T} \!\!\in \mathbf {R^{k\times 1}}\) and M denotes the matrix defined by \(\mathbf {M}_{ij}=<\mathbf {K}^{c}_{i},\mathbf {K}^{c}_{j}>_{F}\), for i,j=1,...,k. We can obtain the weight (ω) by solving this simple quadratic programming problem.

CKA-MKL estimates the weights of \(w_{d} \in \mathbf {R}^{k_{d}\times 1},w_{g} \in \mathbf {R}^{k_{g}\times 1}\), to combine disease \(\left (\mathbf {K}^{d}_{SEM},\mathbf {K}^{d}_{GIP} \in \mathbf {R}^{n\times n}\right)\) and gene \(\left (\mathbf {K}^{g}_{SW},\mathbf {K}^{g}_{GO},\mathbf {K}^{g}_{PPI},\mathbf {K}^{g}_{GIP} \in \mathbf {R}^{m\times m}\right)\) kernels, separately. kd and kg are the number of kernels in disease space and gene space. In order to obtain the optimal kernel matrix \(\mathbf {K}^{*}_{d}\) and \(\mathbf {K}^{*}_{g}\) in the two spaces, first calculate the weights of kernel matrices in each space by Eq. 17, and then combine them by Eq. 14. Here, \(\mathbf {K}^{d}_{ideal} = \mathbf {Y}_{train}\mathbf {Y}^{T}_{train} \in \mathbf {R}^{n\times n}\) in the disease space; and \(\mathbf {K}^{g}_{ideal} = \mathbf {Y}^{T}_{train}\mathbf {Y}_{train} \in \mathbf {R}^{m\times m}\) in the gene space.

Hypergraph learning

In graph theory, a graph represents the pairwise relationship between a group of objects. In traditional graph structures, vertices represent objects, and edges represent relationships between two objects. However, traditional graph structures cannot express more complex relationships. For example, they cannot express more than three relationships in pairs. Hypergraph [31] solves this problem well. In hypergraph theory, this kind of multi-object relationship is represented by using a subset of vertex sets as super edges. In this study, we use hypergraph to establish this higher-order relationship. In Fig. 8 (left), {v1,v2,...,,v7} represents the vertex set, and {v2,v4,v6} are contained in hyperedge e1. Each hyperedge may comprise two or more vertices. The hyperedge will degenerate into a normal edge, when there are only two vertices in the hyperedge.

Fig. 8
figure 8

Hypergraph model in our method

The construction of hypergraph is similar to that of ordinary graph. Hypergraph also needs a vertex set V, an hyperedge set E and the weight of hyperedge \(\phantom {\dot {i}\!}\pmb {w} \in \mathbf {R}^{N_{e}\times 1}\). Here, each hyperedge ei(i=1,2,...,Ne) is given a weight w(ei). The difference is that the hyperedge set of a hypergraph is actually a set of vertices. Therefore, a hypergraph can be represented by G=(V,E,w).

For the hypergraph G, the incidence matrix H conveys the affinity between vertices and hyperedges. And, each element of H can be given by the following formula:

$$ H(v,e)= \left\{\begin{array}{ll} 1& \text{if} v \in e \\ 0& \text{if} v \not\in e \end{array}\right. $$
(18)

The matrix H describes the relationship between vertices and is shown in Fig. 2 (right). Specifically, Hi,j=1 means the vertex vi is included in the hyperedge ej. On the contrary, Hi,j=0 means that the vertex vi is not in the hyperedge ej

In a hypergraph G. The degree of each vertex and hyperedge and the weight of hpyperdege are expressed as follows:

$$ \begin{aligned} d(v) &= \sum\limits_{e \in \mathbf{E}} H(v,e)\\ \delta(e) &=\sum\limits_{v \in \mathbf{V}} H(v,e)\\ w(e_{j}) &= \sum_{i=1}^{k}K^{*}(v_{i},v_{j}) \end{aligned} $$
(19)

where K is the combined kernel.

The hypergraph is constructed using the k Nearest Neighbor (kNN) algorithm. Specifically, each vertex as the center point, and find the k vertices with the largest similarity according to the kernel matrix to form a hyperedge. Assuming that there are N samples, we can construct N hyperedges. In this study, we define the weight of each hyperedge is the sum of kernel values of the k vertices closest to center point, and finally the weight is normalized.

Then, we compute three matrices Dv,De and Dw, where Dv and De are the diagonal matrices of d(v) and d(e). Dw is the matrix of hyperedge weights.

$$ \begin{aligned} &\mathbf{D}_{v} = diag(d)\\ &\mathbf{D}_{e} = diag(\delta)\\ &\mathbf{D}_{w} = diag(w) \end{aligned} $$
(20)

The hypergraph Laplacian matrix Lh [31] is defined as follows:

$$ \begin{aligned} \mathbf{L}^{h} &= \mathbf{I}-\pmb{\Theta}\\ \pmb{\Theta} &= \mathbf{D}^{-1/2}_{v}\mathbf{H}\mathbf{D}_{w}\mathbf{D}^{-1}_{e}\mathbf{H}^{T}\mathbf{D}^{-1/2}_{v} \end{aligned} $$
(21)

where I is the identity matrix.

Consequently, we can obtain the hypergraph Laplacian matrix \(\mathbf {L}^{h}_{d}\) and \(\mathbf {L}^{h}_{g}\) about the disease and gene spaces, respectively.

Dual hypergraph regularized least squares

Baesd on LapRLS method, we propose a novel model to predict the associations of genes and diseases, named Dual Hypergraph Regularized Least Squares (DHRLS), through incorporation of the multiple informations of gene and disease feature spaces into the dual hypergraph regularized least squares framework. The objective function can be written as follow:

$$ \min E(\mathbf{F}^{*}) = ||\mathbf{K}^{*}_{d}\pmb{\alpha}_{d}+\left(\mathbf{K}^{*}_{g}\pmb{\alpha}_{g}\right)^{T}-2\mathbf{Y}_{train}||^{2}_{F} $$
(22)

where \(\mathbf {F}^{*}_{d}=\mathbf {K}^{*}_{d}\pmb {\alpha }_{d}\) and \(\mathbf {F}^{*}_{g}=\mathbf {K}^{*}_{g}\pmb {\alpha }_{g}\). The F could be calculated by \(\mathbf {F}^{*}=\left (\mathbf {F}^{*}_{d}+(\mathbf {F}^{*}_{g})^{T}\right)/2\). F is an average combination of gene and disease space evaluation as the final prediction result.

Then to avoid overfitting of αd and αg to training data, we apply L2 (Tikhonov) regularization to Eq. 22 by adding two terms regarding αd and αg.

$$ \begin{aligned} \min E(\mathbf{F}^{*}) &= ||\mathbf{K}^{*}_{d}\pmb{\alpha}_{d}+\left(\mathbf{K}^{*}_{g}\pmb{\alpha}_{g}\right)^{T}-2\mathbf{Y}_{train}||^{2}_{F}\\ &+\beta\left(||\pmb{\alpha}_{d}||^{2}_{F}+||\pmb{\alpha}_{g}||^{2}_{F}\right) \end{aligned} $$
(23)

where β is a regularization coefficient.

Since previous studies [52] have shown that graph regularization terms are beneficial to improve the prediction effect of the model, graph regularization terms related to genes and diseases are added to the model. According to the local invariance assumption [53], if two data points are close in the intrinsic geometry of the data distribution, then the representations of these two points with respect to the new basis, are also close to each other. This assumption plays an essential role in the development of various kinds of algorithms. In our model, we minimize the distance between the potential feature vectors of two adjacent diseases and genes respectively

$$ \begin{aligned} \min\limits_{\pmb{\alpha}_{d}}\sum_{i,r}^{n}\mathbf{K}^{*}_{d}(i,r)||{\mathbf{F}^{*}_{d}}^{i}-{\mathbf{F}^{*}_{d}}^{r}||^{2}\\ =tr\left(\pmb{\alpha}^{T}_{d}\mathbf{K}^{*}_{d}\mathbf{L}_{d}\mathbf{K}^{*}_{d}\pmb{\alpha}_{d}\right)\\ \min\limits_{\pmb{\alpha}_{g}}\sum_{j,q}^{m}\mathbf{K}^{*}_{g}(j,q)||{\mathbf{F}^{*}_{g}}^{j}-{\mathbf{F}^{*}_{g}}^{q}||^{2}\\ =tr\left(\pmb{\alpha}^{T}_{g}\mathbf{K}^{*}_{g}\mathbf{L}_{g}\mathbf{K}^{*}_{g}\pmb{\alpha}_{g}\right) \end{aligned} $$
(24)

where \({\mathbf {F}^{*}_{d}}^{i}\) is the i-th row vector of \(\mathbf {F}^{*}_{d}=\mathbf {K}^{*}_{d}\pmb {\alpha }_{d}\in \mathbf {R}^{n\times m}, i,r=1,2,...,n\). Similarly, \({\mathbf {F}^{*}_{g}}^{j}\) is the j-th row vector of \(\mathbf {F}^{*}_{g}=\mathbf {K}^{*}_{g}\pmb {\alpha }_{g}\in \mathbf {R}^{m\times n},j,q=1,2,...,m\). \({\mathbf {F}^{*}_{d}}^{i}\) and \({\mathbf {F}^{*}_{g}}^{j}\) mean the representations of the new base. \(\mathbf {K}^{*}_{d}(i,r)\) and \(\mathbf {K}^{*}_{g}(j,q)\) are the weights of two points in two spaces respectively. After adding the graph regular term, the objective function is redefined as follows:

$$ \begin{aligned} \min E(\mathbf{F}^{*}) &= ||\mathbf{K}^{*}_{d}\pmb{\alpha}_{d}+\left(\mathbf{K}^{*}_{g}\pmb{\alpha}_{g}\right)^{T}-2\mathbf{Y}_{train}||^{2}_{F}\\ &+\lambda_{d}tr\left(\pmb{\alpha}^{T}_{d}\mathbf{K}^{*}_{d}\mathbf{L}_{d}\mathbf{K}^{*}_{d}\pmb{\alpha}_{d}\right) \\ &+\lambda_{g}tr\left(\pmb{\alpha}^{T}_{g}\mathbf{K}^{*}_{g}\mathbf{L}_{g}\mathbf{K}^{*}_{g}\pmb{\alpha}_{g}\right)\\ &+\beta\left(||\pmb{\alpha}_{d}||^{2}_{F}+||\pmb{\alpha}_{g}||^{2}_{F}\right) \end{aligned} $$
(25)

where λd and λg are the coefficients of graph regular terms.

We take formula 25 as a model, called Dual Graph Regularized Least Squares (DGRLS). In order to be able to express the high-order relationship between nodes, while improving the prediction effect, Hypergraph Laplacian matrix is applied to our final model DHRLS. Thus, the final objective function can be described as follows:

$$ \begin{aligned} \min E(\mathbf{F}^{*}) &= ||\mathbf{K}^{*}_{d}\pmb{\alpha}_{d}+\left(\mathbf{K}^{*}_{g}\pmb{\alpha}_{g}\right)^{T}-2\mathbf{Y}_{train}||^{2}_{F}\\ &+\lambda_{d}tr\left(\pmb{\alpha}^{T}_{d}\mathbf{K}^{*}_{d}\mathbf{L}^{h}_{d}\mathbf{K}^{*}_{d}\pmb{\alpha}_{d}\right)\\ &+\lambda_{g}tr\left(\pmb{\alpha}^{T}_{g}\mathbf{K}^{*}_{g}\mathbf{L}^{h}_{g}\mathbf{K}^{*}_{g}\pmb{\alpha}_{g}\right)\\ &+\beta\left(||\pmb{\alpha}_{d}||^{2}_{F}+||\pmb{\alpha}_{g}||^{2}_{F}\right) \end{aligned} $$
(26)

where Lh is the hypergraph laplacian matrix, it can be calculated by Eq. 21.

Objective function optimization for DHRLS

We select alternating least squares to estimate αd and αg, and then run alternatingly until convergence.

Optimizingαd Suppose αg are known, to obtain the optimal αd by setting E(F)/αd=0:

$$ {} \begin{aligned} \frac{\partial E(\mathbf{F}^{*})}{\pmb{\alpha}_{d}}&=0\\ \mathbf{K}^{*}_{d}\left(\mathbf{K}^{*}_{d}\pmb{\alpha}_{d}+\pmb{\alpha}^{T}_{g}(\mathbf{K}^{*}_{g})^{T}-2\mathbf{Y}_{train}\right)&+\beta\pmb{\alpha}_{d}+\lambda_{d}\mathbf{K}^{*}_{d}\mathbf{L}^{h}_{d}\mathbf{K}^{*}_{d}\pmb{\alpha}_{d}=0\\ \left(\mathbf{K}^{*}_{d}\mathbf{K}^{*}_{d}+\beta\mathbf{I}+\lambda_{d}\mathbf{K}^{*}_{d}\mathbf{L}^{h}_{d}\mathbf{K}^{*}_{d}\right)\pmb{\alpha}_{d} &= 2\mathbf{K}^{*}_{d}\mathbf{Y}_{train}-\mathbf{K}^{*}_{d}\pmb{\alpha}^{T}_{g}\mathbf{K}^{*}_{g}\\ \pmb{\alpha}_{d} = \left(\mathbf{K}^{*}_{d}\mathbf{K}^{*}_{d}+\beta\mathbf{I}+\lambda_{d}\mathbf{K}^{*}_{d}\mathbf{L}^{h}_{d}\mathbf{K}^{*}_{d}\right)&^{-1}\left(2\mathbf{K}^{*}_{d}\mathbf{Y}_{train}-\mathbf{K}^{*}_{d}\pmb{\alpha}^{T}_{g}\mathbf{K}^{*}_{g}\right) \end{aligned} $$
(27)

Optimizingαg Similarly, suppose αd are known, to obtain the optimal αg by setting E(F)/αg=0:

$$ {} \begin{aligned} \frac{\partial E(\mathbf{F}^{*})}{\pmb{\alpha}_{g}}&=0\\ \mathbf{K}^{*}_{g}\left(\mathbf{K}^{*}_{g}\pmb{\alpha}_{g}+\pmb{\alpha}^{T}_{d}(\mathbf{K}^{*}_{d})^{T}-2\mathbf{Y}^{T}_{train}\right)&+\beta\pmb{\alpha}_{g}+\lambda_{g}\mathbf{K}^{*}_{g}\mathbf{L}^{h}_{g}\mathbf{K}^{*}_{g}\pmb{\alpha}_{g}=0\\ \left(\mathbf{K}^{*}_{g}\mathbf{K}^{*}_{g}+\beta\mathbf{I}+\lambda_{g}\mathbf{K}^{*}_{g}\mathbf{L}^{h}_{g}\mathbf{K}^{*}_{g}\right)\pmb{\alpha}_{g} &= 2\mathbf{K}^{*}_{g}\mathbf{Y}^{T}_{train}-\mathbf{K}^{*}_{g}\pmb{\alpha}^{T}_{d}\mathbf{K}^{*}_{d}\\ \pmb{\alpha}_{g} = \left(\mathbf{K}^{*}_{g}\mathbf{K}^{*}_{g}+\beta\mathbf{I}+\lambda_{g}\mathbf{K}^{*}_{g}\mathbf{L}^{h}_{g}\mathbf{K}^{*}_{g}\right)&^{-1}\left(2\mathbf{K}^{*}_{g}\mathbf{Y}^{T}_{train}-\mathbf{K}^{*}_{g}\pmb{\alpha}^{T}_{d}\mathbf{K}^{*}_{d}\right) \end{aligned} $$
(28)

The final prediction result is by combining the matrices in the two spaces:

$$ \mathbf{F}^{*}=\frac{\mathbf{K}^{*}_{d}\pmb{\alpha}_{d}+(\mathbf{K}^{*}_{g}\pmb{\alpha}_{g})^{T}}{2} $$
(29)

The overview of our method is shown in Table 9.

Table 9 The algorithm of our proposed method

Availability of data and materials

The datasets, codes and corresponding results of our model are available at https://github.com/guofei-tju/DHRLS. Associations for all genes associated with the 960 diseases are available at http://cssb2.biology.gatech.edu/knowgene. Six types of real networks are published at [44].

Abbreviations

DHRLS:

Dual Hypergraph Regularized Least Squares

DGRLS:

Dual Graph Regularized Least Squares

CKA:

Centered Kernel Alignment

MKL:

Multiple Kernel Learning

ALSA:

Alternating Least squares algorithm

MKSVM:

Multiple Kernel Support Vector Machine

CV:

cross-validation

AUPR:

Area Under the Precision-Recall curve

AUC:

Area under the receiver operating characteristic curve

GIP:

Gaussian Interaction Profile

kNN:

k Nearest Neighbor

LapRLS:

Laplacian Regularized Least Squares

PPIs:

Protein-protein interactions

GO:

Gene Ontology

CC:

cellular component

SW:

Smith Waterman

References

  1. Wei L, Liao M, Gao Y, Ji R, He Z, Zou Q. Improved and promising identification of human micrornas by incorporating a high-quality negative set. IEEE/ACM Trans Compu Biol Bioinforma. 2013; 11(1):192–201.

    Article  Google Scholar 

  2. Wei L, Su R, Wang B, Li X, Zou Q, Gao X. Integration of deep feature representations and handcrafted features to improve the prediction of n6-methyladenosine sites. Neurocomputing. 2019; 324:3–9.

    Article  Google Scholar 

  3. Liu H, Ren G, Chen H, Liu Q, Yang Y, Zhao Q. Predicting lncrna–mirna interactions based on logistic matrix factorization with neighborhood regularized. Knowl-Based Syst. 2020; 191:105261.

    Article  Google Scholar 

  4. Wang J, Wang H, Wang X, Chang H. Predicting drug-target interactions via fm-dnn learning. Curr Bioinforma. 2020; 15(1):68–76.

    CAS  Article  Google Scholar 

  5. Huang Y, Yuan K, Tang M, Yue J, Bao L, Wu S, Zhang Y, Li Y, Wang Y, Ou X, et al. Melatonin inhibiting the survival of human gastric cancer cells under er stress involving autophagy and ras-raf-mapk signalling. J Cell Mol Med. 2021; 25(3):1480–92.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  6. Xiao Y, Wu J, Lin Z, Zhao X. A deep learning-based multi-model ensemble method for cancer prediction. Comput Methods Prog Biomed. 2018; 153:1–9.

    Article  Google Scholar 

  7. Mordelet F, Vert J-P. Prodige: Prioritization of disease genes with multitask machine learning from positive and unlabeled examples. BMC Bioinforma. 2011; 12(1):389.

    Article  Google Scholar 

  8. Yu S, Falck T, Daemen A, Tranchevent L-C, Suykens J, De Moor B, Moreau Y. L 2-norm multiple kernel learning and its application to biomedical data fusion. BMC Bioinforma. 2010; 11(1):309.

    Article  CAS  Google Scholar 

  9. Deo R, Musso G, Tasan M, Tang P, Poon A, Yuan C, Felix J, Vasan R, Beroukhim R, De Marco T, et al.Prioritizing causal disease genes using unbiased genomic features. Genome Biol. 2014; 15(12):534.

    PubMed  Article  PubMed Central  Google Scholar 

  10. Yang P, Li X-L, Mei J-P, Kwoh C-K, Ng S-K. Positive-unlabeled learning for disease gene identification. Bioinforma. 2012; 28(20):2640–7.

    CAS  Article  Google Scholar 

  11. Natarajan N, Dhillon I. Inductive matrix completion for predicting gene–disease associations. Bioinforma. 2014; 30(12):60–8.

    Article  CAS  Google Scholar 

  12. Zakeri P, Simm J, Arany A, ElShal S, Moreau Y. Gene prioritization using bayesian matrix factorization with genomic and phenotypic side information. Bioinforma. 2018; 34(13):447–56.

    Article  CAS  Google Scholar 

  13. Zeng X, Ding N, Rodríguez-Patón A, Zou Q. Probability-based collaborative filtering model for predicting gene–disease associations. BMC Med Genet. 2017; 10(5):76.

    Google Scholar 

  14. Singh-Blom U, Natarajan N, Tewari A, Woods J, Dhillon I, Marcotte E. Prediction and validation of gene-disease associations using methods inspired by social network analyses. PloS one. 2013; 8(5):e58977.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. Luo P, Ding Y, Lei X, Wu F-X. deepdriver: predicting cancer driver genes based on somatic mutations using deep convolutional neural networks. Front Genet. 2019; 10:13.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  16. Rao A, Saipradeep V, Joseph T, Kotte S, Sivadasan N, Srinivasan R. Phenotype-driven gene prioritization for rare diseases using graph convolution on heterogeneous networks. BMC Med Genet. 2018; 11(1):57.

    Google Scholar 

  17. Li Y, Kuwahara H, Yang P, Song L, Gao X. Pgcn: Disease gene prioritization by disease and gene embedding through graph convolutional neural networks. bioRxiv. 2019; 532226.

  18. Ding Y, Tang J, Guo F. Identification of drug–target interactions via fuzzy bipartite local model. Neural Comput Applic. 2019:1–17.

  19. Ding Y, Tang J, Guo F. Identification of drug–target interactions via dual laplacian regularized least squares with multiple kernel fusion. Knowl-Based Syst. 2020; 204:106254.

    Article  Google Scholar 

  20. Ding Y, Tang J, Guo F. Identification of drug-side effect association via semisupervised model and multiple kernel learning. IEEE J Biomed Health Inform. 2018; 23(6):2619–32.

    PubMed  Article  PubMed Central  Google Scholar 

  21. Ding Y, Tang J, Guo F. Identification of drug-side effect association via multiple information integration with centered kernel alignment. Neurocomputing. 2019; 325:211–24.

    Article  Google Scholar 

  22. Wang H, Ding Y, Tang J, Guo F. Identification of membrane protein types via multivariate information fusion with hilbert–schmidt independence criterion. Neurocomputing. 2020; 383:257–69.

    Article  Google Scholar 

  23. Shen Y, Tang J, Guo F. Identification of protein subcellular localization via integrating evolutionary and physicochemical information into chou’s general pseaac. J Theor Biol. 2019; 462:230–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  24. Ding Y, Tang J, Guo F. Human protein subcellular localization identification via fuzzy model on kernelized neighborhood representation. Appl Soft Comput. 2020:106596.

  25. Yi Zou, XGLPYDJT HongjieWu, Guo F. Mk-fsvm-svdd: A multiple kernel-based fuzzy svm model for predicting dna-binding proteins via support vector data description. Curr Bioinforma. 2020; 16:274–83.

    Article  CAS  Google Scholar 

  26. Ding Y, Tang J, Guo F. Protein crystallization identification via fuzzy model on linear neighborhood representation. IEEE/ACM Trans Comput Biol Bioinforma. 2019::1.

  27. Zhang J, Zhang Z, Pu L, Tang J, Guo F. Aiepred: an ensemble predictive model of classifier chain to identify anti-inflammatory peptides. IEEE/ACM Trans Comput Biol Bioinforma. 2020::1.

  28. 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–434.

    Google Scholar 

  29. Wang F, Huang Z-A, Chen X, Zhu Z, Wen Z, Zhao J, Yan G-Y. Lrlshmda: Laplacian regularized least squares for human microbe–disease association prediction. Sci Rep. 2017; 7(1):1–11.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  30. Jiang L, Xiao Y, Ding Y, Tang J, Guo F. Fkl-spa-laprls: an accurate method for identifying human microrna-disease association. BMC Genomics. 2018; 19(10):11–25.

    Google Scholar 

  31. Zhou D, Huang J, Schölkopf B. Learning with hypergraphs: Clustering, classification, and embedding. In: Advances in Neural Information Processing Systems. Cambridge: MIT Press: 2007. p. 1601–1608.

    Google Scholar 

  32. Wu W, Kwong S, Zhou Y, Jia Y, Gao W. Nonnegative matrix factorization with mixed hypergraph regularization for community detection. Inf Sci. 2018; 435:263–81.

    Article  Google Scholar 

  33. Xu X-X, Dai L-Y, Kong X-Z, Liu J-X. A low-rank representation method regularized by dual-hypergraph laplacian for selecting differentially expressed genes. Hum Hered. 2019; 84(1):1–13.

    Article  CAS  Google Scholar 

  34. Bai S, Zhang F, Torr P. Hypergraph convolution and hypergraph attention. Pattern Recog. 2021; 110:107637.

    Article  Google Scholar 

  35. Zhang R, Zou Y, Ma J. Hyper-sagnn: a self-attention based graph neural network for hypergraphs. arXiv preprint arXiv:1911.02613. 2019.

  36. Ding Y, Jiang L, Tang J, Guo F. Identification of human microrna-disease association via hypergraph embedded bipartite local model. Comput Biol Chem. 2020; 89:107369.

    CAS  PubMed  Article  Google Scholar 

  37. Lü L, Zhou T. Link prediction in complex networks: A survey. Physica A Stat Mech its Appl. 2011; 390(6):1150–70.

    Article  Google Scholar 

  38. Holme P, Liljeros F, Edling C, Kim B. Network bipartivity. Phys Rev E. 2003; 68(5):056107.

    Article  CAS  Google Scholar 

  39. Kunegis J, De Luca E, Albayrak S. The link prediction problem in bipartite networks. In: International Conference on Information Processing and Management of Uncertainty in Knowledge-based Systems. Berlin Heidelberg: Springer: 2010. p. 380–9.

    Google Scholar 

  40. Lu Y, Wang L, Lu J, Yang J, Shen C. Multiple kernel clustering based on centered kernel alignment. Pattern Recog. 2014; 47(11):3656–64.

    Article  Google Scholar 

  41. Zhou H, Skolnick J. A knowledge-based approach for predicting gene–disease associations. Bioinforma. 2016; 32(18):2831–8.

    CAS  Article  Google Scholar 

  42. Ezzat A, Zhao P, Wu M, Li X-L, Kwoh C-K. Drug-target interaction prediction with graph regularized matrix factorization. IEEE/ACM Trans Comput Biol Bioinforma. 2016; 14(3):646–56.

    Article  Google Scholar 

  43. Zheng X, Ding H, Mamitsuka H, Zhu S. Collaborative matrix factorization with multiple similarities for predicting drug-target interactions. In: Proceedings of the 19th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. ACM: 2013. p. 1025–33.

  44. Wang W, Chen X, Jiao P, Jin D. Similarity-based regularized latent feature model for link prediction in bipartite networks. Sci Rep. 2017; 7(1):1–12.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  45. Xia Z, Wu L-Y, Zhou X, Wong S. Semi-supervised drug-protein interaction prediction from heterogeneous biological spaces. In: BMC Systems Biology. BioMed Central: 2010. p. 6.

  46. Lowe H, Barnett G. Understanding and using the medical subject headings (mesh) vocabulary to perform literature searches. Jama. 1994; 271(14):1103–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  47. Wang D, Wang J, Lu M, Song F, Cui Q. Inferring the human microrna functional similarity and functional network based on microrna-associated diseases. Bioinforma. 2010; 26(13):1644–50.

    CAS  Article  Google Scholar 

  48. Ashburner M, Ball C, Blake J, Botstein D, Butler H, Cherry J, Davis A, Dolinski K, Dwight S, Eppig J, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000; 25(1):25–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  49. Sherman B, Lempicki R, et al. Systematic and integrative analysis of large gene lists using david bioinformatics resources. Nat Protocol. 2009; 4(1):44.

    Article  CAS  Google Scholar 

  50. Yu G, Li F, Qin Y, Bo X, Wu Y, Wang S. Gosemsim: an r package for measuring semantic similarity among go terms and gene products. Bioinforma. 2010; 26(7):976–8.

    CAS  Article  Google Scholar 

  51. Yamanishi Y, Araki M, Gutteridge A, Honda W, Kanehisa M. Prediction of drug–target interaction networks from the integration of chemical and genomic spaces. Bioinforma. 2008; 24(13):232–40.

    Article  CAS  Google Scholar 

  52. Gu Q, Zhou J, Ding C. Collaborative filtering: Weighted nonnegative matrix factorization incorporating user and item graphs. In: Proceedings of the 2010 SIAM International Conference on Data Mining. Columbus: SIAM: 2010. p. 199–210.

    Google Scholar 

  53. Cai D, He X, Han J, Huang T. Graph regularized nonnegative matrix factorization for data representation. IEEE Trans Pattern Anal Mach Intell. 2010; 33(8):1548–60.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

No application.

Funding

This work is supported by a grant from the National Natural Science Foundation of China (NSFC 61772362, 61902271 and 61972280), and National Key R&D Program of China (2020YFA0908400), and the Natural Science Research of Jiangsu Higher Education Institutions of China (19KJB520014).

Author information

Authors and Affiliations

Authors

Contributions

HY conceived and designed the experiments; YD performed the experiments and analyzed the data; FG wrote the paper; JT reviewed the manuscript. All authors have read and approved the whole manuscript.

Corresponding authors

Correspondence to Yijie Ding or Fei Guo.

Ethics declarations

Ethics approval and consent to participate

No application.

Consent for publication

No application.

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.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yang, H., Ding, Y., Tang, J. et al. Identifying potential association on gene-disease network via dual hypergraph regularized least squares. BMC Genomics 22, 605 (2021). https://doi.org/10.1186/s12864-021-07864-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-021-07864-z

Keywords

  • Gene-disease association network
  • Hypergraph learning
  • Dual Laplacian regularized least squares
  • Bipartite network
  • Multiple kernel learning