Prioritizing potential circRNA biomarkers for bladder cancer and bladder urothelial cancer based on an ensemble model

Abstract

Bladder cancer is the most common cancer of the urinary system. Bladder urothelial cancer accounts for 90% of bladder cancer. These two cancers have high morbidity and mortality rates worldwide. The identification of biomarkers for bladder cancer and bladder urothelial cancer helps in their diagnosis and treatment. circRNAs are considered oncogenes or tumor suppressors in cancers, and they play important roles in the occurrence and development of cancers. In this manuscript, we developed an Ensemble model, CDA-EnRWLRLS, to predict circRNA-Disease Associations (CDA) combining Random Walk with restart and Laplacian Regularized Least Squares, and further screen potential biomarkers for bladder cancer and bladder urothelial cancer. First, we compute disease similarity by combining the semantic similarity and association profile similarity of diseases and circRNA similarity by combining the functional similarity and association profile similarity of circRNAs. Second, we score each circRNA-disease pair by random walk with restart and Laplacian regularized least squares, respectively. Third, circRNA-disease association scores from these models are integrated to obtain the final CDAs by the soft voting approach. Finally, we use CDA-EnRWLRLS to screen potential circRNA biomarkers for bladder cancer and bladder urothelial cancer. CDA-EnRWLRLS is compared to three classical CDA prediction methods (CD-LNLP, DWNN-RLS, and KATZHCDA) and two individual models (CDA-RWR and CDA-LRLS), and obtains better AUC of 0.8654. We predict that circHIPK3 has the highest association with bladder cancer and may be its potential biomarker. In addition, circSMARCA5 has the highest association with bladder urothelial cancer and may be its possible biomarker.

1 Introduction

Bladder cancer is considered to be the most common cancer in the urinary system (Kamat et al., 2016). It is the fourth most common malignant tumor in men and the eighth most common in women in the Western world. In the United States and Europe, it accounts for 5%–10% among all malignancies in men. The risk with the bladder cancer infection at less than 75 years is 2%–4% in men and 0.5%–1% for women (Kirkali et al., 2005). The incidence of bladder cancer has been increasing (Kamat et al., 2016). The majority of patients with bladder cancer suffer from the less aggressive non-muscle-invasive disease, while 30% of patients suffer from muscle-invasive disease (Lopez-Beltran et al., 2021; Tran et al., 2021; Yang et al., 2021).

Bladder cancer has a metastatic potential, and thus presents a worse prognosis. It is usually grouped into three pathological categories: bladder urothelial carcinoma, bladder squamous cell carcinoma, and bladder adenocarcinoma (Black and Black, 2020; Lopez-Beltran et al., 2021). Bladder urothelial carcinoma accounts for over 90% among all cases of bladder cancer. Furthermore, bladder urothelial carcinoma can be categorized into muscle-invasive bladder cancer, which accounts for about 75% of all cases, and non-muscle-invasive bladder cancer (Kirkali et al., 2005). The all-stage five-year survival rate of bladder urothelial cancer remains approximately 80% (Lopez-Beltran et al., 2021).

Recently, the treatment of bladder cancer has obtained great progresses worldwide. Besides traditional surgical resection, radiotherapy, and chemotherapy, immunotherapy is also a promising avenue for bladder cancer treatment (Gao et al., 2021; Mancini et al., 2021). However, postoperative recurrence and distant metastasis cause five-year survival rates to still be very low for advanced bladder cancer (Fabiano et al., 2021; Roviello et al., 2021). Advanced disease or relapse of radical cystectomy is closely associated with the poor outcomes (Nouhaud et al., 2021). The first-line therapy of metastatic bladder urothelial cancer usually adopts cisplatin-based combinations, and has been unaltered over the last decades (Powles et al., 2021; Renner et al., 2021; Walia et al., 2021). Unfortunately, almost all patients with bladder urothelial cancer will finally progress and die from bladder cancer, despite their initial response to cisplatin-based combinations (Bin Riaz et al., 2021; Lopez-Beltran et al., 2021). Consequently, inferring potential biomarkers for bladder cancer is a good way to diagnose and treat it (Peng et al., 2017; Peng et al., 2018).

With the advance of sequencing technology, there are now massive amounts of RNA data (Ozsolak and Milos, 2011; Peng et al., 2020; Yang et al., 2020; Peng et al., 2022a), which help the prognosis and treatment of various diseases (Xu et al., 2020; Li et al., 2021). Circular RNAs (circRNAs) are a class of single-stranded noncoding RNA molecules that are lack of terminal 5′ caps and 3′ poly(A) tails (He et al., 2017). circRNAs are widely distributed in various organisms. They have circular features, and thus demonstrate more resistance to degradation by exonucleases and stronger stability than linear RNAs (Xia et al., 2018; Li G. et al., 2020). The estimated total number of circRNAs is approximately 1% of one of poly (A) molecules. In addition, the expression levels of the majority of circRNAs are estimated to be 5%–10% of the corresponding linear RNAs (Jeck and Sharpless, 2014; Zhang J. et al., 2021).

Although circRNAs were found in 1976, they were originally considered to be functionless by-products from aberrant RNA splicing and thus did not obtain enough attention over the past 3 decades. However, with the rapid advance of high-throughput sequencing technologies, massive differentially expressed circRNAs have been increasing discovered in human normal and malignant cells (Zhang et al., 2018; Li G. et al., 2020; Yang et al., 2021). circRNAs exist widely in various tissues, serum, and urine. The expression profiles of circRNAs demonstrate strong specificity in cell types, tissues, and developmental stages (Yang et al., 2021). Furthermore, circRNAs can regulate transcription or splicing, translate proteins, interact with RNA-binding proteins, and act as miRNA sponges (Sheng et al., 2018). A large body of evidence shows that circRNAs have dense associations with various diseases, including neurological dysfunction, cardiovascular diseases, and cancer. Here, circRNAs, as miRNA sponges, can inhibit the regulation from downstream cancer target genes. For instance, circCDR1as and circMTO1 can control gene regulation and further indirectly stimulate or inhibit tumors by binding to miR-7 and miR-9 (Vromman et al., 2021).

circRNAs have abundant associations with cancers and thus can be used as candidate cancer biomarkers (Zhang et al., 2018). An increasing amount of evidence has reported that circRNAs present in human biofluids and exosomes, and are a class of potential biomarkers of noninvasive liquid biopsies. For instance, circ-ZEB1.33 is overexpressed in hepatocellular cancer and has close links with the survival of hepatocellular cancer patients (Gong et al., 2018). In particular, substantial studies have demonstrated that circRNAs play key roles in the carcinogenesis and progression of bladder cancer. For example, circRNAs Cdr1as performs anti-oncogenic functions in bladder cancer through microRNA 135a (Li et al., 2018), BCRC-3 suppresses bladder cancer proliferation via sponging miR-182-5p/p27 (Xie et al., 2018), MYLK and circPDSS1 promote bladder cancer progression separately by modulating VEGFA/VEGFR2 signaling pathway and down-regulating miR-16 (Zhong et al., 2017; Yu et al., 2020), PRMT5 supports metastasis of bladder urothelial cancer through Sponging miR-30c (Chen et al., 2018), circSLC8A1 suppresses bladder cancer progression through regulating PTEN (Lu et al., 2019), and circMTO1 inhibits bladder cancer metastasis through sponging miR-221 (Li G. et al., 2019).

Many computational methods have been proposed to identify possible CDAs and further discovered possible circRNA biomarkers for various complex diseases including cancers by case studies (Wang CC. et al., 2021). For example, Lei et al. (Lei et al., 2018) designed a path weighted-based CDA prediction approach (PWCDA). Li et al. (Li Y. et al., 2019; Li J. et al., 2020) explored two CDA identification models (NCPCDA and DWNCPCDA) based on network consistency projection. Zhang et al. (Zhang et al., 2019) developed a linear neighborhood label propagation algorithm for CDA identification. Deepthi et al. (Deepthi and Jereesh (2020) used autoencoder and deep neural network and explored an ensemble model to predict CDAs. Lu et al. (Lu et al. (2021) improved CDA prediction using convolutional and recurrent neural networks. Wang et al. (Wang et al., 2020; Wang et al. 2021b; Wang et al., 2021c) proposed three CDA identification methods (GCNCDA, MGRCDA, and SGANRDA) based on graph convolutional network, metagraph recommendation, and semi-supervised generative adversarial network, respectively. These methods efficiently predicted possible CDAs.

In this study, inspired by computational CDA prediction methods, we develop an ensemble model, CDA-RWLRLS, to find potential circRNA biomarkers for bladder cancer and bladder urothelial cancer based on known CDAs. CDA-EnRWLRLS first computes circRNA similarity by integrating their functional similarity and association profile similarity, and it computes disease similarity by integrating their semantic similarity and association profile similarity. Second, CDA-EnRWLRLS computes the association probability for each circRNA-disease pair based on random walk with restart and Laplacian regularized least squares. Third, the prediction results obtained by these two models are integrated by the soft voting method. We finally use the proposed CDA-EnRWLRLS model to identify possible circRNAs associated with bladder cancer and bladder urothelial cancer.

2 Materials and methods

2.1 Materials

2.1.1 Human circRNA-disease associations

circRNA-disease association data can be downloaded from the circR2Disease database (Fan et al., 2018a). This database provides 739 experimentally confirmed CDAs from 661 circRNAs and 100 diseases. We remove redundant elements related to mice and rats and achieve a human circRNA-disease association dataset containing 650 associations between 585 circRNAs and 88 diseases. In particular, suppose that and separately denote the sets of circRNAs and diseases, then we construct a binary matrix to depict circRNA-disease associations by Eq. 1:

2.1.2 Disease semantic similarity

Many studies have computed disease semantic similarity to screen credible noncoding RNAs for a query disease. Inspired by these methods, we investigate disease similarity to improve the prediction performance. Disease semantic similarity can be computed based on corresponding disease ontology. The disease ontology is often represented using a directed acyclic graph and can be downloaded from http://disease-ontology.org/. For two query diseases and corresponding ontology term sets from the two diseases and , their semantic similarity can be scored by the “doSim” function in the DOSE software package, which can be downloaded from http://www.bioconductor.org/packages/release/bioc/html/DOSE.html (Yu et al., 2015). Finally, we compute the semantic similarity matrix among diseases.

2.1.3 circRNA functional similarity

To measure the functional similarity between two circRNAs, we utilize the semantic similarity of two diseases linking to the two circRNAs. In particular, suppose that and denote the disease groups linking to circRNAs and , the functional similarity between and can be computed by Eq. 2:

andwhere denotes the similarity between disease linking to circRNA and disease set linking to circRNA .

2.2 Methods

In this manuscript, we develop circRNA-Disease Association prioritization method (CDA-EnRWLRLS) by an Ensemble of Random Walk with restart and Laplacian Regularization Least Squares. First, CDA-EnRWLRLS measures circRNA functional similarity and disease semantic similarity. Second, it computes association profile similarity of circRNAs and diseases, respectively. Third, functional similarity and association profile similarity of circRNAs are combined to obtain the final circRNA similarity. Similarly, disease similarity is fused. Fourth, random walk with restart and Laplacian regularization least squares are used to score each circRNA-disease pair. Fifth, the final association score matrix is obtained by integrating the results from random walk with restart and Laplacian regularization least squares based on the soft voting strategy. Finally, CDA-EnRWLRLS is applied to find possible circRNA biomarkers for bladder cancer and bladder urothelial cancer. The flowchart of CDA-EnRWLRLS is shown in Figure 1.

FIGURE 1

2.2.1 Association profile similarity of circRNAs and diseases

For two diseases with known ontology terms, we can compute their semantic similarity based on their ontology terms. However, semantic similarity computation may fail for two diseases without ontology terms. Thus, we introduce association profile similarity to further complement similarity measurement of circRNAs and diseases.

Suppose that the association profile of a circRNA is represented as the th row of a CDA matrix . describes information from all diseases associated with . Association profile similarity between two circRNAs (i.e., ) can be computed by Eq. 4:where is bandwidth parameter and set as the default value of 1.

Similarly, association profile similarity between two diseases (i.e., ) can be computed by Eq. 5:where indicates bandwidth parameter and set as the default value of 1.

2.2.2 Similarity fusion

circRNA functional similarity , disease semantic similarity , and association profile similarity of circRNAs and diseases ( and ) are fused to obtain the final circRNA similarity matrix and disease similarity by Eqs 6, 7:

The parameter is used to balance the importance between functional similarity and association profile similarity of circRNAs in Eq. 6 and is used to balance the important between semantic similarity and association profile similarity of diseases in Eq. 7.

2.2.3 Random walk with restart for CDA prediction

Random walk algorithm has been widely used and obtained better performance in various association prediction fields (Peng et al., 2021a). In this study, we utilize Random Walk with Restart for CDA prediction on the heterogeneous circRNA-disease network (CDA-RWR). We first train the random walk with restart model on the CDA dataset and screen possible CDAs with the highest association probability from unknown circRNA-disease pairs on the dataset.

First, circRNA similarity network , disease similarity network , and CDA network are used to build a heterogeneous circRNA-disease network., , and correspond to adjacency matrices of the three networks, respectively. Consequently, the heterogeneous circRNA-disease network can be represented as: , where is the transpose of .

Second, we compute the transition probability of random walk on the heterogeneous circRNA-disease network. Suppose that denote the transition matrix, where and separately indicate the walk within the circRNA network and the disease network, and separately represent the jump from the circRNA network to the disease network and the disease network to the circRNA network. For a known jumping probability from the circRNA network to the disease network or from the disease network to the circRNA network, the transition probability from circRNAs to can be calculated by Eq. 8:

The transition probability from circRNA to disease can be calculated by Eq. 9:

The transition probability from diseases to can be calculated by Eq. 10:

The transition probability from disease to circRNA can be calculated by Eq. 11:

For a query circRNA/disease, it can either stay in the current network with a restart probability or jump to another network graph. Consequently, we can compute association probability for each circRNA-disease pair at the -th step by Eq. 12:where denotes the association probability matrix at the -th step, denotes the initial probability and , where and indicate the initial probability on the circRNA and disease network, respectively. When we want to discover possible circRNAs associated with a query disease , it is regarded as a seed in the disease network. Consequently, is assigned as 1 and other disease nodes are 0, thereby building the initial probability matrix of the disease network . All nodes in the circRNA network are assigned as an equal probability whose sum is 1. The parameter is used to balance the importance of the circRNA network and the disease network.

2.3 Laplacian regularized least squares for CDA prediction

We can calculate association probability for each circRNA-disease pair based on random walk with restart. However, for random walk with restart, the jump probability is measured by known CDAs and the circRNA and disease similarity matrices. For a circRNA in a CDA network, if two other circRNAs and have the equal similarity with , and will contribute to the jump between nodes at an equal probability. However, the circRNA that exhibits lower similarities with other circRNAs should have more contribution to the jump. Thus, we further use Laplacian regularized least squares (Shen et al., 2022) to compute association probability for each circRNA-disease pair.

First, we compute the circRNA Laplacian matrix and the disease Laplacian matrix by Eqs 13, 14:where / indicates the diagonal matrix of circRNA/disease similarity matrix and / is the summation of the -th/-th row of /.

Second, we define the loss functions of Laplacian regularization least squares in the circRNA and disease spaces based on the Laplacian matrices and by Eqs 15, 16, respectively:where , , and separately indicate the transposes of , , and , indicates the Frobenius norm, and and indicate trade-off parameters. The Laplacian regularized least square models (15) and (16) can be solved by Eqs 17, 18:

Finally, the association probability for each circRNA-disease pair by Laplacian regularized least squares can be computed by Eq. 19:

2.4 Ensemble learning for CDA prediction

Ensemble learning integrates multiple results from individual models and demonstrates better performance compared to individual models (Zhou et al., 2021a; Peng et al., 2022b). Therefore, in this study, we develop an ensemble learning model by combining random walk with restart and Laplacian regularized least squares to improve the CDA’s prediction performance by Eq. 20:where denotes the predicted final CDA score matrix, and denote the computed CDA probability matrices based on random walk with restart and Laplacian regularized least squares, respectively. is used to weigh the importance of results computed by the above two models.

3 Experiments

3.1 Experimental settings

For similarity computation, the weights between biological feature similarity and association profile similarity and are set as 0.5. For random walk with restart, the restart probability is set as 0.2, and and are set as 0.1 and 0.6, respectively. For Laplacian regularized least squares, both and are set as 0.95 and 0.2, respectively. For ensemble learning model, is set as 0.3. The parameters in other three comparative methods are set as defaults provided by the corresponding methods. We conduct 5-fold cross validation for 10 times. The final prediction performance is from the average value of the 10 experiments. AUC (area under the receiver operating characteristic curve) has been widely used to evaluate the performance of CDA prediction methods. Larger AUC denotes better performance. Thus, we use AUC to measure the performance of our proposed method.

3.2 Performance comparison with five CDA prediction methods

Several comparative experiments are conducted to measure the performance of our proposed CDA-EnRWLRLS model. CD-LNLP (Zhang et al., 2019), DWNN-RLS (Yan et al., 2018), KATZHCDA (Fan et al., 2018b), and CDA-EnRWLRLS are conducted on the preprocessed CDA dataset. CD-LNLP (Zhang et al., 2019) is a linear neighborhood label propagation-based algorithm for CDA prediction. DWNN-RLS (Yan et al., 2018) used regularized least squares to predict possible CDAs. KATZHCDA (Fan et al., 2018b) discovered CDA candidates based on the KATZ measurement (Zhou et al., 2020). Figure 2 shows the AUC values computed by these four CDA prediction methods.

FIGURE 2

From Figure 2, we can find that CDA-EnRWLRLS is significantly better than CD-LNLP (Zhang et al., 2019), DWNN-RLS (Yan et al., 2018), and KATZHCDA (Fan et al., 2018b) based on the AUC value. Compared to the three models, CDA-EnRWLRLS obtains the highest AUC of 0.8654, outperforming 7.60%, 24.86%, and 0.25%, respectively. In particular, DWNN-RLS used regularized least squares with Kronecker product kernel for CDA prediction. Disease similarity was computed by their semantic similarity and Gaussian association profile similarity. Meanwhile, circRNA similarity was computed by their Gaussian association profiles. CDA-EnRWLRLS uses an ensemble model to identify possible CDAs. Similar to DWNN-RLS, CDA-EnRWLRLS computes disease similarity. However, CDA-EnRWLRLS computes circRNA similarity by their functional similarity and Gaussian association profile similarity. Furthermore, CDA-EnRWLRLS still computes association score between each circRNA-disease pair using random walk with restart except Laplacian regularized least squares and integrates the results from the two models by the soft voting technique. Therefore, CDA-EnRWLRLS outperforms DWNN-RLS, which demonstrates its powerful CDA prediction ability.

3.3 Performance evaluation of ensemble learning model with individual models

Our proposed CDA-EnRWLRLS model is an ensemble of two state-of-the-art models (i.e., random walk with restart and Laplacian regularized least squares). To evaluate the performance of ensemble learning model and individual models, we conducted 5-fold cross validation experiment for CDA-EnRWLRLS and random walk with restart (CDA-RWR) and Laplacian regularized least squares (CDA-LRLS) on the CDA dataset. Figure 3 shows the AUC values computed by CDA-EnRWLRLS, CDA-RWR, and CDA-LRLS. From Figure 3, we can find that CDA-EnRWLRLS obtains better AUC than two individual models, CDA-RWR and CDA-LRLS, which shows that the proposed ensemble learning-based model can outperforms individual models.

FIGURE 3

3.4 Evaluation of parameter sensitivity

In this study, we ensemble two individual models, random walk with restart and Laplacian regularized least squares. However, the two models may have different effects on the CDA prediction performance. To evaluate their effect on the performance, we consider in the range of [0.1, 0.9] with stride of 0.1. The results are shown in Figure 4.

FIGURE 4

From Figure 4, we can find that AUC computed by CDA-EnRWLRLS gradually increases when the parameter is from 0.1 to 0.3. Its computed AUCs gradually decrease when the parameter is from 0.3 to 0.9. In other words, CDA-EnRWLRLS obtains the best AUC when the parameter is 0.3. Thus, the parameter is finally set as 0.3.

3.5 Case study

We consequently compute the association score for each circRNA-disease pair. In particular, we compute association abilities between all circRNAs and bladder cancer and bladder urothelial cancer to analyze any possible associations between these circRNAs and the two cancers, and to further screen for potential circRNA biomarkers for them.

3.5.1 circRNA biomarker analysis for bladder cancer

Bladder cancer is a heterogeneous disease with high morbidity and mortality rates (Kamat et al., 2016). It has been estimated that about 73,510 new cases of bladder cancer were diagnosed in the United States in 2012. During the same period, about 14,880 patients died from bladder cancer (Clark et al., 2013). To analyze circRNA biomarkers for bladder cancer, we compute association between all circRNAs and bladder cancer after training CDA-EnRWLRLS. Table 1 gives the top 20 circRNAs that are predicted to have the highest association scores with bladder cancer.

In the CDA dataset, 15 circRNAs are known to associate with bladder cancer among 585 circRNAs. From Table 1, we can find that the 15 circRNAs are predicted to have the highest association scores with bladder cancer and are ranked as the top 15. Furthermore, we predict that circHIPK3 may associate with bladder cancer with the ranking of 16. Furthermore, circHIPK3 is a promising cancer-related circRNA (Zhang et al., 2020). It can regulate cell growth through sponging multiple miRNAs (Zheng et al., 2016). For instance, circHIPK3 can regulate cell proliferation and migration in hepatocellular cancer by sponging miR-124 (Chen X. et al., 2018), modulate autophagy in STK11 mutant lung cancer (Chen et al., 2020), and promote glioma progression as a prognostic marker (Jin et al., 2018). The overexpression of circHIPK3 can accelerate the proliferation and invasion of prostate cancer cells (Cai et al., 2019). Its inhibition can block angiotensin II-induced cardiac fibrosis (Ni et al., 2019). In this study, we infer that circHIPK3 may be a biomarker of bladder cancer and need experimental validation. Figure 5 shows the association information between the top 20 circRNAs with bladder cancer.

FIGURE 5

3.5.2 circRNA biomarker analysis for bladder urothelial cancer

Over 90% bladder cancer is bladder urothelial cancer. Bladder urothelial cancer is a common malignancy with high morbidity and mortality worldwide (Cancer Genome Atlas Research Network, 2014). In the United Sates, bladder urothelial cancer is one of the main histologic subtypes (Clark et al., 2013). However, no molecularly targeted agent has been applied to the treatment, until now. To infer potential circRNA biomarkers for bladder urothelial cancer, we compute association scores between all circRNAs and bladder urothelial cancer using CDA-EnRWLRLS. Table 2 gives the top 20 circRNAs that are predicted to have the highest association scores with bladder urothelial cancer.

In the CDA dataset, only one circRNA, circHIPK3, associates with bladder urothelial cancer among all potential 585 circRNAs. We predict that SMARCA5 may associate with bladder urothelial cancer with the ranking of 2. SMARCA5 is a member of the ISWI family that is involved in chromatin remodeling. It can regulate chromosome remodeling through diverse mechanisms, hinder cell proliferation, and assist apoptosis by sponging miRNAs. Its expression may boost the susceptibility of cells to chemotherapy, boost the sensitivity of cancer detection, promote early diagnosis, and help the treatment of chemotherapy-resistant cancers (Qin and Wan, 2022). Its expression level has a certain association with clinical features of many cancers. For instance, SMARCA5 can promote cell proliferation in bladder cancer and prostate cancer (Tan et al., 2019), suppress colorectal cancer progression (Miao et al., 2020), inhibit tumor metastasis in cervical cancer (Zhang X. et al., 2021) and inhibit cell proliferation, migration, and invasion in non-small cell lung cancer (Wang et al., 2019), and boost cell migration and invasion as well as inhibit cell apoptosis in bladder cancer (Kong et al., 2017; Tan et al., 2019). Many studies have reported that circSMARCA5 plays a key role in the occurrence and development of cancer. Moreover, it also serves as a reliable indicator of tumor screening or cancer prognosis evaluation (Qin and Wan, 2022). Therefore, SMARCA5 is a diagnostic and prognostic biomarker of cancer and has obtained wide attention. In this study, we predict that SMARCA5 may be potential biomarker of bladder urothelial cancer; however, this needs validation. Figure 6 shows the association information between the top 20 circRNAs with bladder urothelial cancer.

FIGURE 6

4 Discussion and conclusion

Bladder cancer, including bladder urothelial cancer, is a common and complex disease. These cancers have caused high morbidity and mortality. The identification of biomarkers for bladder cancer and bladder urothelial cancer can help in their prognosis and treatment. In this manuscript, we developed an ensemble learning model, CDA-EnRWLRLS, to discover potential circRNA biomarkers for the two cancers based on CDA association prediction.

CDA-EnRWLRLS first computes circRNA similarity and disease similarity by fusing semantic similarity and association profile similarity of diseases and functional similarity and association profile similarity of circRNAs. Second, it scores each circRNA-disease pair by random walk with restart and Laplacian regularized least squares, respectively. Third, the results computed by random walk with restart and Laplacian regularized least squares are integrated by the soft voting approach based on ensemble learning. Finally, it is applied to discover potential circRNA biomarkers for bladder cancer and bladder urothelial cancer.

CDA-EnRWLRLS is compared to three classical CDA prediction methods (CD-LNLP, DWNN-RLS, and KATZHCDA) and two individual models (CDA-RWR and CDA-LRLS). The results show that CDA-EnRWLRLS computes relatively better AUC, which demonstrates its relatively powerful CDA prediction ability. We predict that circHIPK3 and SMARCA5 may be potential biomarkers of bladder cancer and bladder urothelial cancer, respectively.

CDA-EnRWLRLS has two advantages: on the one hand, it better fuses biological features and association features of diseases and circRNAs; while on the other hand, it combines two individual classical association prediction models to obtain the powerful association prediction performance from different bioinformatics tools. Although CDA-EnRWLRLS computed better CDA inference ability, the circRNA functional similarity was calculated indirectly by disease semantic similarity. Moreover, its prediction performance needs further improvement. In the future, we will consider biological features of circRNAs and develop more efficient machine learning, especially ensemble learning models (Zhou et al., 2021a; Peng et al., 2022a) and deep learning models (Peng et al., 2021b; Zhou et al., 2021b; Sun et al., 2022; Yang et al., 2022) to discover potential biomarkers for bladder cancer and bladder urothelial cancer.

References

• 1

  Bin RiazI.KhanA. M.CattoJ. W. F.HussainS. A. (2021). Bladder cancer: Shedding light on the most promising investigational drugs in clinical trials.. Expert Opin. Investig. Drugs30 (8), 837–855. 10.1080/13543784.2021.1948999

  • CrossRef
  • Google Scholar

• 2

  BlackA. J.BlackP. C. (2020). Variant histology in bladder cancer: Diagnostic and clinical implications. Transl. Cancer Res.9 (10), 6565–6575. 10.21037/tcr-20-2169

  • CrossRef
  • Google Scholar

• 3

  CaiC.ZhiY.WangK.ZhangP.JiZ.XieC.et al (2019). CircHIPK3 overexpression accelerates the proliferation and invasion of prostate cancer cells through regulating miRNA-338-3p.. Onco. Targets. Ther.12, 3363–3372. 10.2147/OTT.S196931

  • CrossRef
  • Google Scholar

• 4

  Cancer Genome Atlas Research Network (2014). Comprehensive molecular characterization of urothelial bladder carcinoma[J]. Nature507 (7492), 315. 10.1038/nature12965

  • CrossRef
  • Google Scholar

• 5

  ChenG.ShiY.LiuM.SunJ. (2018). circHIPK3 regulates cell proliferation and migration by sponging miR-124 and regulating AQP3 expression in hepatocellular carcinoma.. Cell. Death Dis.9 (2), 1–13. 10.1038/s41419-017-0204-3

  • CrossRef
  • Google Scholar

• 6

  ChenX.ChenR. X.WeiW. S.LiY. H.FengZ. H.TanL.et al (2018). PRMT5 circular RNA promotes metastasis of urothelial carcinoma of the bladder through sponging miR-30c to induce epithelial-mesenchymal transition.. Clin. Cancer Res.24 (24), 6319–6330. 10.1158/1078-0432.CCR-18-1270

  • CrossRef
  • Google Scholar

• 7

  ChenX.MaoR.SuW.YangX.GengQ.GuoC.et al (2020). Circular RNA circHIPK3 modulates autophagy via MIR124-3p-STAT3-PRKAA/AMPKα signaling in STK11 mutant lung cancer.. Autophagy16 (4), 659–671. 10.1080/15548627.2019.1634945

  • CrossRef
  • Google Scholar

• 8

  ClarkP. E.AgarwalN.BiagioliM. C.EisenbergerM. A.GreenbergR. E.HerrH. W.et al (2013). Bladder cancer.. J. Natl. Compr. Canc. Netw.11 (4), 446–475. 10.6004/jnccn.2013.0059

  • CrossRef
  • Google Scholar

• 9

  DeepthiK.JereeshA. S. (2020). An ensemble approach for CircRNA-disease association prediction based on autoencoder and deep neural network.. Gene762, 145040. 10.1016/j.gene.2020.145040

  • CrossRef
  • Google Scholar

• 10

  FabianoE.DurduxC.DufourB.MejeanA.ThiounnN.ChretienY.et al (2021). Long-term outcomes after bladder-preserving tri-modality therapy for patients with muscle-invasive bladder cancer.. Acta Oncol.60 (6), 794–802. 10.1080/0284186X.2021.1915498

  • CrossRef
  • Google Scholar

• 11

  FanC.LeiX.FangZ.JiangQ.WuF-X. (2018a). CircR2Disease: A manually curated database for experimentally supported circular RNAs associated with various diseases[J]. Database2018, bay044. 10.1093/database/bay044

  • CrossRef
  • Google Scholar

• 12

  FanC.LeiX.WuF. X. (2018b). Prediction of CircRNA-disease associations using KATZ model based on heterogeneous networks.. Int. J. Biol. Sci.14 (14), 1950–1959. 10.7150/ijbs.28260

  • CrossRef
  • Google Scholar

• 13

  GaoS.YangX.XuJ.QiuN.ZhaiG. (2021). Nanotechnology for boosting cancer immunotherapy and remodeling tumor microenvironment: The horizons in cancer treatment. ACS Nano15 (8), 12567–12603. 10.1021/acsnano.1c02103

  • CrossRef
  • Google Scholar

• 14

  GongY.MaoJ.WuD. I.WangX.LiL.ZhuL.et al (2018). Circ-ZEB1.33 promotes the proliferation of human HCC by sponging miR-200a-3p and upregulating CDK6.. Cancer Cell. Int.18 (1), 116–119. 10.1186/s12935-018-0602-3

  • CrossRef
  • Google Scholar

• 15

  HeJ.XieQ.XuH.LiJ.LiY. (2017). Circular RNAs and cancer.. Cancer Lett.396, 138–144. 10.1016/j.canlet.2017.03.027

  • CrossRef
  • Google Scholar

• 16

  JeckW. R.SharplessN. E. (2014). Detecting and characterizing circular RNAs.. Nat. Biotechnol.32 (5), 453–461. 10.1038/nbt.2890

  • CrossRef
  • Google Scholar

• 17

  JinP.HuangY.ZhuP.ZouY.ShaoT.WangO. (2018). CircRNA circHIPK3 serves as a prognostic marker to promote glioma progression by regulating miR-654/IGF2BP3 signaling.. Biochem. Biophys. Res. Commun.503 (3), 1570–1574. 10.1016/j.bbrc.2018.07.081

  • CrossRef
  • Google Scholar

• 18

  KamatA. M.HahnN. M.EfstathiouJ. A.LernerS. P.MalmstromP. U.ChoiW.et al (2016). Bladder cancer.. Lancet388 (10061), 2796–2810. 10.1016/S0140-6736(16)30512-8

  • CrossRef
  • Google Scholar

• 19

  KirkaliZ.ChanT.ManoharanM.AlgabaF.BuschC.ChengL.et al (2005). Bladder cancer: Epidemiology, staging and grading, and diagnosis.. Urology66 (6), 4–34. 10.1016/j.urology.2005.07.062

  • CrossRef
  • Google Scholar

• 20

  KongZ.WanX.ZhangY.ZhangP.ZhangY.ZhangX.et al (2017). Androgen-responsive circular RNA circSMARCA5 is up-regulated and promotes cell proliferation in prostate cancer.. Biochem. Biophys. Res. Commun.493 (3), 1217–1223. 10.1016/j.bbrc.2017.07.162

  • CrossRef
  • Google Scholar

• 21

  LiG.LuoJ.WangD.LiangC.XiaoQ.DingP.et al (2020). Potential circRNA-disease association prediction using DeepWalk and network consistency projection.. J. Biomed. Inf.112, 103624. 10.1016/j.jbi.2020.103624

  • CrossRef
  • Google Scholar

• 22

  LiG.YueY.LiangC.XiaoQ.DingP.LuoJ. (2019a). Ncpcda: Network consistency projection for circRNA-disease association prediction.. RSC Adv.9 (57), 33222–33228. 10.1039/c9ra06133a

  • CrossRef
  • Google Scholar

• 23

  LiJ.SunD.PuW.WangJ.PengY. (2020). Circular RNAs in cancer: Biogenesis, function, and clinical significance.. Trends Cancer6 (4), 319–336. 10.1016/j.trecan.2020.01.012

  • CrossRef
  • Google Scholar

• 24

  LiP.YangX.YuanW.YangC.ZhangX.HanJ.et al (2018). CircRNA-Cdr1as exerts anti-oncogenic functions in bladder cancer by sponging MicroRNA-135a.. Cell.. Physiol. biochem.46 (4), 1606–1616. 10.1159/000489208

  • CrossRef
  • Google Scholar

• 25

  LiT.HuangT.GuoC.WangA.ShiX.MoX.et al (2021). Genomic variation, origin tracing, and vaccine development of SARS-CoV-2: A systematic review.. Innovation.2 (2), 100116. 10.1016/j.xinn.2021.100116

  • CrossRef
  • Google Scholar

• 26

  LiY.WanB.LiuL.ZhouL.ZengQ. (2019b). Circular RNA circMTO1 suppresses bladder cancer metastasis by sponging miR-221 and inhibiting epithelial-to-mesenchymal transition.. Biochem. Biophys. Res. Commun.508 (4), 991–996. 10.1016/j.bbrc.2018.12.046

  • CrossRef
  • Google Scholar

• 27

  Lopez-BeltranA.CimadamoreA.BlancaA.MassariF.VauN.ScarpelliM.et al (2021). Immune checkpoint inhibitors for the treatment of bladder cancer. Cancers13 (1), 131. 10.3390/cancers13010131

  • CrossRef
  • Google Scholar

• 28

  LuC.ZengM.WuF. X.LiM.WangJ. (2021). Improving circRNA–disease association prediction by sequence and ontology representations with convolutional and recurrent neural networks. Bioinformatics36 (24), 5656–5664. 10.1093/bioinformatics/btaa1077

  • CrossRef
  • Google Scholar

• 29

  LuQ.LiuT.FengH.YangR.ZhaoX.ChenW.et al (2019). Circular RNA circSLC8A1 acts as a sponge of miR-130b/miR-494 in suppressing bladder cancer progression via regulating PTEN.. Mol. Cancer18 (1), 111–113. 10.1186/s12943-019-1040-0

  • CrossRef
  • Google Scholar

• 30

  ManciniM.RighettoM.NoessnerE. (2021). Checkpoint inhibition in bladder cancer: Clinical expectations, current evidence, and proposal of future strategies based on a tumor-specific immunobiological approach.. Cancers13 (23), 6016. 10.3390/cancers13236016

  • CrossRef
  • Google Scholar

• 31

  MiaoX.XiZ.ZhangY.LiZ.HuangL.XinT.et al (2020). Circ-SMARCA5 suppresses colorectal cancer progression via downregulating miR-39-3p and upregulating ARID4B.. Dig. Liver Dis.52 (12), 1494–1502. 10.1016/j.dld.2020.07.019

  • CrossRef
  • Google Scholar

• 32

  NiH.LiW.ZhugeY.XuS.WangY.ChenY.et al (2019). Inhibition of circHIPK3 prevents angiotensin II-induced cardiac fibrosis by sponging miR-29b-3p.. Int. J. Cardiol.292, 188–196. 10.1016/j.ijcard.2019.04.006

  • CrossRef
  • Google Scholar

• 33

  NouhaudF. X.ChakrounM.LenormandC.ouzaidI.PeyronnetB.GrynA.et al (2021). Comparison of the prognosis of primary vs. progressive muscle invasive bladder cancer after radical cystectomy: Results from a large multicenter study. Urologic Oncol. Seminars Orig. Investigations39 (3), 195.e1–195.e6. 10.1016/j.urolonc.2020.09.006

  • CrossRef
  • Google Scholar

• 34

  OzsolakF.MilosP. M. (2011). RNA sequencing: Advances, challenges and opportunities.. Nat. Rev. Genet.12 (2), 87–98. 10.1038/nrg2934

  • CrossRef
  • Google Scholar

• 35

  PengL.ShenL.XuJ.TianX.LiuF.WangJ.et al (2021a). Prioritizing antiviral drugs against SARS-CoV-2 by integrating viral complete genome sequences and drug chemical structures[J]. Sci. Rep.11 (1), 1–11. 10.1038/s41598-021-83737-5

  • CrossRef
  • Google Scholar

• 36

  PengL.TanJ.TianX.ZhouL. (2022b). EnANNDeep: An ensemble-based lncRNA–protein interaction prediction framework with adaptive k-nearest neighbor classifier and deep models[J]. Interdiscip. Sci. Comput. Life Sci., 1–24. 10.1007/s12539-021-00483-y

  • CrossRef
  • Google Scholar

• 37

  PengL.TianX.TianG.XuJ.HuangX.WengY.et al (2020). Single-cell RNA-seq clustering: Datasets, models, and algorithms.. RNA Biol.17 (6), 765–783. 10.1080/15476286.2020.1728961

  • CrossRef
  • Google Scholar

• 38

  PengL.WangF.WangZ.TanJ.HuangL.TianX.et al (2022a). Cell-cell communication inference and analysis in the tumour microenvironments from single-cell transcriptomics: Data resources and computational strategies.. Brief. Bioinform.23 (4), bbac234. 10.1093/bib/bbac234

  • CrossRef
  • Google Scholar

• 39

  PengL. H.ChenY. Q.MaN.ChenX. (2017). Narrmda: Negative-aware and rating-based recommendation algorithm for miRNA-disease association prediction.. Mol. Biosyst.13 (12), 2650–2659. 10.1039/c7mb00499k

  • CrossRef
  • Google Scholar

• 40

  PengL. H.SunC. N.GuanN. N.QiangJ.ChenX. (2018). Hnmda: Heterogeneous network-based miRNA-disease association prediction.. Mol. Genet. Genomics293 (4), 983–995. 10.1007/s00438-018-1438-1

  • CrossRef
  • Google Scholar

• 41

  PengL.WangC.TianX.ZhouL.LiK. (2021b). Finding lncRNA-protein interactions based on deep learning with dual-net neural architecture[J]. IEEE/ACM Trans. Comput. Biol. Bioinforma.2021, 29. 10.1109/TCBB.2021.3116232

  • CrossRef
  • Google Scholar

• 42

  PowlesT.CsősziT.ÖzgüroğluM.MatsubaraN.GecziL.ChengS. Y. S.et al (2021). Pembrolizumab alone or combined with chemotherapy versus chemotherapy as first-line therapy for advanced urothelial carcinoma (KEYNOTE-361): A randomised, open-label, phase 3 trial.. Lancet. Oncol.22 (7), 931–945. 10.1016/S1470-2045(21)00152-2

  • CrossRef
  • Google Scholar

• 43

  QinH.WanR. (2022). The emerging roles of circSMARCA5 in cancer[J]. J. Oncol., 2022.

  • Google Scholar

• 44

  RennerA.BurottoM.ValdesJ. M.RomanJ. C.Walton-DiazA. (2021). Neoadjuvant immunotherapy for muscle invasive urothelial bladder carcinoma: Will it change current standards?Ther. Adv. Urol.13, 17562872211029779. 10.1177/17562872211029779

  • CrossRef
  • Google Scholar

• 45

  RovielloG.CatalanoM.SantiR.PalmieriV. E.VanniniG.GalliI. C.et al (2021). Immune checkpoint inhibitors in urothelial bladder cancer: State of the art and future perspectives.. Cancers13 (17), 4411. 10.3390/cancers13174411

  • CrossRef
  • Google Scholar

• 46

  ShenL.LiuF. X.HuangL.LiuG. Y.ZhouL. Q.PengL. H. (2022). VDA-RWLRLS: An anti-SARS-CoV-2 drug prioritizing framework combining an unbalanced bi-random walk and Laplacian regularized least squares. Comput. Biol. Med.140, 105119. 10.1016/j.compbiomed.2021.105119

  • CrossRef
  • Google Scholar

• 47

  ShengJ. Q.LiuL.WangM. R.LiP. Y. (2018). Circular RNAs in digestive system cancer: Potential biomarkers and therapeutic targets.. Am. J. Cancer Res.8 (7), 1142–1156.

  • Google Scholar

• 48

  SunF.SunJ.ZhaoQ. (2022). A deep learning method for predicting metabolite-disease associations via graph neural network.. Brief. Bioinform.23 (4), bbac266. 10.1093/bib/bbac266

  • CrossRef
  • Google Scholar

• 49

  TanY.ZhangT.LiangC. (2019). Circular RNA SMARCA5 is overexpressed and promotes cell proliferation, migration as well as invasion while inhibits cell apoptosis in bladder cancer.. Transl. Cancer Res.8 (5), 1663–1671. 10.21037/tcr.2019.08.08

  • CrossRef
  • Google Scholar

• 50

  TranL.XiaoJ. F.AgarwalN.DuexJ. E.TheodorescuD. (2021). Advances in bladder cancer biology and therapy.. Nat. Rev. Cancer21 (2), 104–121. 10.1038/s41568-020-00313-1

  • CrossRef
  • Google Scholar

• 51

  VrommanM.VandesompeleJ.VoldersP. J. (2021). Closing the circle: Current state and perspectives of circular RNA databases.. Brief. Bioinform.22 (1), 288–297. 10.1093/bib/bbz175

  • CrossRef
  • Google Scholar

• 52

  WaliaA. S.SweisR. F.AgarwalP. K.KaderA. K.ModiP. K. (2021). Cost-effectiveness of immune checkpoint inhibitors in urothelial carcinoma-A review.. Cancers14 (1), 73. 10.3390/cancers14010073

  • CrossRef
  • Google Scholar

• 53

  WangC. C.HanC. D.ZhaoQ.ChenX. (2021a). Circular RNAs and complex diseases: From experimental results to computational models.. Brief. Bioinform.22 (6), bbab286. 10.1093/bib/bbab286

  • CrossRef
  • Google Scholar

• 54

  WangL.YanX.YouZ. H.ZhouX.LiH. Y.HuangY. A. (2021c). Sganrda: Semi-supervised generative adversarial networks for predicting circRNA-disease associations.. Brief. Bioinform.22 (5), bbab028. 10.1093/bib/bbab028

  • CrossRef
  • Google Scholar

• 55

  WangL.YouZ. H.HuangD. S.LiJ-Q. (2021b). Mgrcda: Metagraph recommendation method for predicting CircRNA-disease association[J]. IEEE Trans. Cybern.2021. 10.1109/TCYB.2021.3090756

  • CrossRef
  • Google Scholar

• 56

  WangL.YouZ. H.LiY. M.ZhengK.HuangY. A. (2020). Gcncda: A new method for predicting circRNA-disease associations based on graph convolutional network algorithm.. PLoS Comput. Biol.16 (5), e1007568. 10.1371/journal.pcbi.1007568

  • CrossRef
  • Google Scholar

• 57

  WangY.LiH.LuH.QinY. (2019). Circular RNA SMARCA5 inhibits the proliferation, migration, and invasion of non-small cell lung cancer by miR-19b-3p/HOXA9 axis.. Onco. Targets. Ther.12, 7055–7065. 10.2147/OTT.S216320

  • CrossRef
  • Google Scholar

• 58

  XiaL.SongM.SunM.WangF.YangC. (2018).. Circular RNAs, 171–187. 10.1007/978-981-13-1426-1_14Circular RNAs as biomarkers for cancer[J]

  • CrossRef
  • Google Scholar

• 59

  XieF.LiY.WangM.HuangC.TaoD.ZhengF.et al (2018). Circular RNA BCRC-3 suppresses bladder cancer proliferation through miR-182-5p/p27 axis.. Mol. Cancer144, 1–12. 10.1186/s12943-018-0892-z

  • CrossRef
  • Google Scholar

• 60

  XuJ.CaiL.LiaoB.ZhuW.YangJ. (2020). CMF-impute: An accurate imputation tool for single-cell RNA-seq data.. Bioinformatics36 (10), 3139–3147. 10.1093/bioinformatics/btaa109

  • CrossRef
  • Google Scholar

• 61

  YanC.WangJ.WuF. X. D. W. N. N-R. L. S. (2018). Regularized least squares method for predicting circRNA-disease associations[J]. BMC Bioinforma.19 (19), 73–81.

  • Google Scholar

• 62

  YangJ.JuJ.GuoL.JiB.ShiS.YangZ.et al (2022). Prediction of HER2-positive breast cancer recurrence and metastasis risk from histopathological images and clinical information via multimodal deep learning.. Comput. Struct. Biotechnol. J.20, 333–342. 10.1016/j.csbj.2021.12.028

  • CrossRef
  • Google Scholar

• 63

  YangJ.PengS.ZhangB.HoutenS.SchadtE.ZhuJ.et al (2020). Human geroprotector discovery by targeting the converging subnetworks of aging and age-related diseases.. Geroscience42 (1), 353–372. 10.1007/s11357-019-00106-x

  • CrossRef
  • Google Scholar

• 64

  YangX.YeT.LiuH.LvP.DuanC.WuX.et al (2021). Expression profiles, biological functions and clinical significance of circRNAs in bladder cancer[J]. Mol. cancer20 (1), 1–25. 10.1186/s12943-020-01300-8

  • CrossRef
  • Google Scholar

• 65

  YuG.WangL. G.YanG. R.HeQ. Y. (2015). Dose: An R/bioconductor package for disease ontology semantic and enrichment analysis.. Bioinformatics31 (4), 608–609. 10.1093/bioinformatics/btu684

  • CrossRef
  • Google Scholar

• 66

  YuQ.LiuP.HanG.XueX.MaD. (2020). CircRNA circPDSS1 promotes bladder cancer by down-regulating miR-16.. Biosci. Rep.40 (1), BSR20191961. 10.1042/BSR20191961

  • CrossRef
  • Google Scholar

• 67

  ZhangJ.HouL.ZuoZ.JiP.ZhangX.XueY.et al (2021a). Comprehensive profiling of circular RNAs with nanopore sequencing and CIRI-long.. Nat. Biotechnol.39 (7), 836–845. 10.1038/s41587-021-00842-6

  • CrossRef
  • Google Scholar

• 68

  ZhangW.YuC.WangX.LiuF. (2019). Predicting CircRNA-disease associations through linear neighborhood label propagation method. Ieee Access7, 83474–83483. 10.1109/access.2019.2920942

  • CrossRef
  • Google Scholar

• 69

  ZhangX.ZhangQ.ZhangK.WangF.QiaoX.CuiJ. (2021b). Circ SMARCA5 inhibited tumor metastasis by interacting with SND1 and downregulating the YWHAB gene in cervical cancer. Cell. Transpl.30, 096368972098378. 10.1177/0963689720983786

  • CrossRef
  • Google Scholar

• 70

  ZhangY.LiuQ.LiaoQ. (2020). CircHIPK3: A promising cancer-related circular RNA.. Am. J. Transl. Res.12 (10), 6694–6704.

  • Google Scholar

• 71

  ZhangZ.YangT.XiaoJ. (2018). Circular RNAs: Promising biomarkers for human diseases.. EBioMedicine34, 267–274. 10.1016/j.ebiom.2018.07.036

  • CrossRef
  • Google Scholar

• 72

  ZhaoQ.YangY.RenG.GeE.FanC. (2019). Integrating bipartite network projection and KATZ measure to identify novel CircRNA-disease associations.. IEEE Trans. Nanobioscience18 (4), 578–584. 10.1109/TNB.2019.2922214

  • CrossRef
  • Google Scholar

• 73

  ZhengQ.BaoC.GuoW.LiS.ChenJ.ChenB.et al (2016). Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs.. Nat. Commun.7 (1), 1–13. 10.1038/ncomms11215

  • CrossRef
  • Google Scholar

• 74

  ZhongZ.HuangM.LvM.HeY.DuanC.ZhangL.et al (2017). Circular RNA MYLK as a competing endogenous RNA promotes bladder cancer progression through modulating VEGFA/VEGFR2 signaling pathway.. Cancer Lett.403, 305–317. 10.1016/j.canlet.2017.06.027

  • CrossRef
  • Google Scholar

• 75

  ZhouL.WangJ.LiuG.LuQ.DongR.TianG.et al (2020). Probing antiviral drugs against SARS-CoV-2 through virus-drug association prediction based on the KATZ method.. Genomics112 (6), 4427–4434. 10.1016/j.ygeno.2020.07.044

  • CrossRef
  • Google Scholar

• 76

  ZhouL.DuanQ.TianX.TangJ.PengL. H. (2021a). LPI-HyADBS: A hybrid framework for lncRNA-protein interaction prediction integrating feature selection and classification[J]. BMC Bioinforma.22 (1), 1–31.

  • Google Scholar

• 77

  ZhouL.WangZ.TianX.PengL. (2021b). LPI-deepGBDT: A multiple-layer deep framework based on gradient boosting decision trees for lncRNA-protein interaction identification.. BMC Bioinforma.22, 479. 10.1186/s12859-021-04399-8

  • CrossRef
  • Google Scholar