Next Article in Journal
Molecular Mechanism of Blood Pressure Regulation through the Atrial Natriuretic Peptide
Next Article in Special Issue
MFIDMA: A Multiple Information Integration Model for the Prediction of Drug–miRNA Associations
Previous Article in Journal
Muting Bacterial Communication: Evaluation of Prazosin Anti-Quorum Sensing Activities against Gram-Negative Bacteria Pseudomonas aeruginosa, Proteus mirabilis, and Serratia marcescens
Previous Article in Special Issue
Intelligent Epileptic Seizure Detection and Classification Model Using Optimal Deep Canonical Sparse Autoencoder
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

SGCNCMI: A New Model Combining Multi-Modal Information to Predict circRNA-Related miRNAs, Diseases and Genes

1
School of Information Engineering, Xijing University, Xi’an 710123, China
2
College of Grassland and Environment Sciences, Xinjiang Agricultural University, Urumqi 830052, China
3
School of Computer Science, Northwestern Polytechnical University, Xi’an 710129, China
*
Author to whom correspondence should be addressed.
Biology 2022, 11(9), 1350; https://doi.org/10.3390/biology11091350
Submission received: 11 July 2022 / Revised: 21 August 2022 / Accepted: 8 September 2022 / Published: 13 September 2022
(This article belongs to the Special Issue Advanced Computational Models for Clinical Decision Support)

Abstract

:

Simple Summary

With the development of circRNA–miRNA-mediated models, circRNAs have been shown to play a prominent role in the development and treatment of diseases such as cancer, and unearthing potential miRNA-associated circRNAs may provide new insights and ideas for the diagnosis and treatment of complex diseases such as cancer. Large-scale prediction using computer technology can provide an a priori guide to biological experiments and save costs. This paper presents the third computational method in this field with the highest accuracy to date, and we also collected and integrated high-quality datasets from the current database, which we believe will allow future computational innovations to develop.

Abstract

Computational prediction of miRNAs, diseases, and genes associated with circRNAs has important implications for circRNA research, as well as provides a reference for wet experiments to save costs and time. In this study, SGCNCMI, a computational model combining multimodal information and graph convolutional neural networks, combines node similarity to form node information and then predicts associated nodes using GCN with a distributive contribution mechanism. The model can be used not only to predict the molecular level of circRNA–miRNA interactions but also to predict circRNA–cancer and circRNA–gene associations. The AUCs of circRNA—miRNA, circRNA–disease, and circRNA–gene associations in the five-fold cross-validation experiment of SGCNCMI is 89.42%, 84.18%, and 82.44%, respectively. SGCNCMI is one of the few models in this field and achieved the best results. In addition, in our case study, six of the top ten relationship pairs with the highest prediction scores were verified in PubMed.

Graphical Abstract

1. Introduction

Circular RNA (circRNA) is a special kind of single-stranded circular endogenous non-coding RNA (ncRNA). Recent research shows that endogenous circRNAs are widely distributed in mammalian cells and involved in transcriptional and posttranscriptional gene expression regulation [1]. CircRNA was first discovered in RNA viruses as early as 1976 [2], and in 1979, Hsu et al. provided electron microscopic evidence for the circular form of RNA [3]. Over the following three decades, only a handful of circRNAs were discovered by chance [4,5,6], and due to their low levels of expression, circRNAs were typically considered to be products of “noise” of an abnormal RNA splicing process, which resulted in circRNAs not receiving corresponding attention.
However, since 2010, with the development of RNA-seq technologies and specialized computational pipelines, many circRNAs have been widely recognized and discovered in eukaryotes, such as mice [7], archaea [8], and humans [9]. With the progress in circRNA research, many circRNAs have been proven to present tissue-specific expression patterns and have specific biological functions [10]. Emerging experimental results show that endogenous circRNAs widely exist in mammals and can work as miRNA sponges, which means that circRNAs reverse the inhibitory effect of the miRNA on its target gene and consequently repress their function [11]. At present, many types of research have indicated the association between miRNA sponges (circRNAs) and human diseases.
CircRNAs have a prominent role in cancer diagnosis and treatment [12]. For example, in bladder cancer studies, circ-ITCH acted as a miRNA sponge to inhibit bladder cancer progression by directly regulating p21 and PTEN in combination with miR-17 and miR-224. Circ-ITCH expression was also lower than normal in bladder cancer tissues [13]. In addition, high expression of another circRNA, circ-TFRC, was detected in bladder cancer patients, which means that circ-TFRC promotes bladder cancer progression by binding to miR-107 [14]. CircCCDC9 expression was significantly lower than normal in gastric cancer tissue samples, and the study confirmed that circCCDC9 inhibits the progression of gastric cancer by regulating CAV1 in combination with miR-6729-3p [15]. CircRNA also plays an important role in the development of renal cell carcinoma (RCC) by regulating CDKN3/E2F1 in combination with miR-127-3p [16]. These results suggest that investigating circRNA–miRNA interactions could be key to diagnosing and addressing complex diseases such as cancer.
Compared with traditional biological experimental methods, which are limited to small scales and require lots of labor and time, using computational models to predict the association between molecules can provide the basis for biological experiments at a low cost. At present, many computational methods have been proposed and applied to predict the correlation between different molecules. For example, Wang et al. proposed a model, SAEMDA, through a new unsupervised training method named Stacked Autoencoder to predict miRNA–disease associations [17]. Ren et al. developed a model named BioChemDDI, which combines a Natural Language Processing algorithm and Hierarchical Representation Learning to effectively extract information, employed Similarity Network Fusion to fuse multiple features, and finally applied a deep neural network to obtain the predicted results [18]. Wang et al. proposed a new method that extracts deep features of molecular similarity through a deep convolutional neural network and sends them to an extreme learning machine classifier to identify potential circRNA–miRNA associations [19]. Such computational methods have achieved gratifying results and provided an experimental basis for further wet experiments.
Compared with related fields, where new circRNA molecules are constantly discovered and the nomenclature is not fully standardized, there are a few computational methods that predict associations between circRNAs and miRNAs. However, with the rapid development of high-throughput sequencing technology, a large number of databases have been developed to store circRNA-related information, such as circR2Disease [20], circRNAdisease [21], circbank [22], and circBase [23]. The circR2Disease database is a high-quality database containing detailed information about circRNA. The latest version contains about 750 circRNA–disease associations between more than 600 circRNAs and 100 diseases. The circRNAdisease database manually collects verified circRNA–disease pairs from the PubMed database by retrieving circRNA and disease keywords. Circbank is a comprehensive database that contains multiple characteristics of circRNA; more than 140,000 circRNAs from different sources can be retrieved by users from the circbank database. The circBase database is one of the early databases to collect circRNA information, including circRNA data, evidence to support circRNA expression, and scripts for identifying known and new circRNAs in sequencing data. The establishment of these databases has provided the materials for predicting associations between circRNAs and miRNAs by using computational methods.
At present, only few models have been proposed, and the predicted results were confirmed in PubMed. Compared with other fields, there are few computational prediction models for circRNA–miRNA interaction prediction. Therefore, it is urgent to develop new and effective prediction methods for circRNA–miRNA association prediction.
According to our understanding, there are some obstacles to using computational methods to predict circRNA–miRNA interactions: (i) The length of circRNA and miRNA sequences varies greatly, resulting in redundancy or sparsity of biological information collected. (ii) A network composed of confirmed circRNA–miRNA associations is difficult to connect, which means it is difficult to extract effective features from relatively isolated nodes. (iii) The data on circRNA–miRNA interactions are scattered among different databases, so it is difficult to collect comprehensive and reliable data. To solve these problems, we developed a model, SGCNCMI, to predict circRNA–miRNA interactions based on multi-source feature extraction and graph representation learning with a layer contribution mechanism. Specifically, we first adopt a K-mer algorithm to extract the internal attribute features in the sequence by taking the most appropriate K value for different RNAs, and to make full use of RNA molecular biological information, two kinds of kernel functions are added to enrich semantic descriptors. Secondly, we introduce the Sparse Autoencoder (SAE) with a sparsity penalty term to process semantic descriptors to obtain the most valuable molecular biological attribute information. Next, we apply a multilayer graph convolutional neural network (GCN) to project the circRNA–miRNA interaction network into a new space to capture non-linear interactions and hidden associations. Meanwhile, we include a layer contribution mechanism in the graph convolutional layer to ensure the maximum contribution of GCN in each layer. Finally, the predicted score of each pair of circRNA–miRNA is obtained from the inner product of the corresponding potential vectors.
Notably, our model supports training and prediction using two types of training data, one based on circRNA–miRNA molecular sequences and known association data and the other based on circRNA as a cancer marker. This means that our model can be trained and predicted from the perspective of both potential molecular relationships and data on associations between clinical disease and markers.
As a result, in a five-fold cross-validation experiment to measure the ability of the model, 89.42% AUC and 88.87% AUPR were obtained by SGCNCMI, and in the circRNA–miRNA interaction dataset test, the performance of SGCNCMI exceeded that of the only other model at present. In addition, 84.18% AUC and 84.83% AUPR were obtained by SGCNCMI in the circRNA–cancer dataset test, and 82.44% AUC and 85.55% AUPR were obtained in the circRNA–gene dataset. Meanwhile, 7 of the 10 pairs with the top predicted scores of the circRNA–miRNA interaction dataset test was verified in PubMed. Obviously, our model, SGCNCMI, is one of the few accurate and reliable prediction models in the field of circRNA–miRNA interaction prediction and is expected to become a powerful candidate model for biological experiments.

2. Materials and Methods

2.1. Dataset

As research progresses, a number of positive circRNA–miRNA correlations have been identified, and various databases have been established. The CircR2Cancer database [24] is an online database that gathers experimentally validated circRNA–cancer and circRNA–miRNA associations reported in published papers. After rigorous screening, we obtained 318 circRNA–miRNA relational pairs between 238 circRNAs and 230 miRNAs.
At present, the techniques for predicting target gene binding sites are well developed, allowing the selection of candidates that closely match the binding sites, with high accuracy for most binding sites, and the vast majority of these predictions were eventually validated in subsequent experiments. Predicting target gene binding sites is already widely used in a variety of methods and tools; for example, CircInteractome [25] uses a well-established TargetScan Perl script to analyze miRNAs that may be associated with circRNA. These data are extremely valuable. The circBank database [22] performs binding site predictions for 140,790 human circRNAs and 1917 miRNAs using Miranda [26] and TargetScan [27] techniques, resulting in 42,917 relationships with more than five binding sites and 3545 relationships with more than one binding site. We selected the top 9589 pairs of circRNA–miRNA relationship pairs with the highest scores. These data were used in the first computational method [28] in the field and partially validated in PubMed.
After combining data from both databases, we eventually obtained 9905 pairs of high-quality relationships for training in our methods, and for ease of description, we identified this dataset as CMI-9905.
To test SGCNCMI’s ability to predict the association between markers and underlying diseases, we downloaded 1049 experimentally supported circRNA–cancer relationship pairs from the Lnc2Cancer database [29] of 743 circRNAs and 70 cancers.
In addition, we downloaded circRNA–gene-associated data from the TransCirc [30] database and selected the top 2000 pairs with the highest confidence scores as training data.

2.2. CircRNA and miRNA Sequence Similarity Based on K-mer

Counting RNA sequences’ K-mers (substrings of length k) is not only an important and common step in bioinformatics analysis but also widely used in computational methods [31,32]. Related studies have indicated that RNA sequences contain abundant biological information. Converting sequence information into a digital vector is an important method to obtain molecular biological information in order to fully explore hidden features in RNA sequences. The K-mer sparse matrix is used to represent RNAs’ attribute features in our model.
For a circRNA sequence, we apply the best 5-mers as the window to scan the sequence, moving one nucleotide at a time. Due to there being four different nucleotides in circRNA, the window of 5-mers will produce 45 vector representations for each circRNA molecule. Therefore, the K-mer matrix of circRNA can be represented as follows:
K M circ R N A = 2346 × 4 5
For a miRNA sequence, with an average length of 21 nucleotides, the scan window we use is 2-mers to obtain the best vector representations, and the K-mer matrix of miRNA is defined as:
K M m i R N A = 962 × 4 2
The details of the K-mer algorithm are shown in Figure 1.

2.3. Similarity for CircRNA and miRNA

RNAs that can bind to the same molecule often have the same binding sites, which means that a potential unknown association can be inferred by analyzing RNA molecules with the same function. In order to fully express the biological characteristics of RNA molecules, we introduce two kinds of similarities (RNA Gaussian interaction profile kernel similarity and RNA sigmoid kernel similarity) as RNA semantic descriptors.
Firstly, we construct a bipartite graph BC×M to represent the 9905 associations between circRNA and miRNA interaction pairs for 2346 circRNAs and 962 miRNAs. In the matrix BC×M, C and M represent the number of circRNAs and miRNAs. When circRNA i is related to miRNA j, the value of Bi×j is equal to 1 and otherwise equal to 0. Each row and column represent circRNA and miRNA interaction profiles, respectively; the interaction profile binary vector LP(Ci) of circRNA Ci is the row corresponding to the circRNA in the adjacent matrix BC×M, and the GIP kernel of each circRNA can be calculated as:
G c i r c R N A ( C i , C j ) = exp ( α c | | L P ( C i ) L P ( C j ) | | 2 )
where Ci and Cj denote circRNA i and circRNA j, GcircRNA (Ci, Cj) is the GIP kernel similarity between circRNA i and circRNA j, and αc is a variable parameter that controls the bandwidth of the GIP kernel, which is defined as follows:
α c = α c ´ / ( 1 n c i = 1 n c | | L P ( C i ) | | 2 )
In this experiment, αc’ is defined as equal to 0.5.
Similarly, the GIP kernel similarity between miRNA mi and miRNA mj is calculated as
G m i R N A ( M i , M j ) = exp ( α m | | L P ( M i ) L P ( M j ) | | 2 )
α m = α m ´ / ( 1 n m i = 1 n m | | L P ( M i ) | | 2 )
The sigmoid kernel of each circRNA is defined as follows:
K c i r c R N A ( C i , C j ) = tanh { β [ ρ ( C i ) ] × κ [ ρ ( C j ) ] }
where β = 1/V, and V is the dimension of original input data.
In the same way, the sigmoid kernel of each miRNA is defined by the formula below:
K m i R N A ( M i , M j ) = tanh { β [ ρ ( M i ) ] × κ [ ρ ( M j ) ] }

2.4. Integrating Attributes and Similarity for circRNA and miRNA

Feature fusion can incorporate more meaningful information from different aspects, which can comprehensively reflect the characteristics of the circRNA and miRNA. In this section, we construct the characteristic fusion matrices of circRNA and miRNA. First, the different types of circRNA similarity (GIPKS and sigmoid kernel) matrixes are combined into one matrix called FC(ci,cj) by the following formula:
F C ( C i , C j ) = G c i r c R N A ( c i , c j ) + K c i r c R N A ( c i , c j ) 2 G c i r c R N A ( c i , c j ) , K c i r c R N A ( c i , c j ) 0 G c i r c R N A ( c i , c j ) + K c i r c R N A ( c i , c j ) o t h e r w i s e
In the same way, the miRNA similarity matrix is defined as
F M ( M i , M j ) = G m i R N A ( m i , m j ) + K m i R N A ( m i , m j ) 2 G m i R N A ( m i , m j ) , K m i R N A ( m i , m j ) 0 G m i R N A ( m i , m j ) + K m i R N A ( m i , m j ) o t h e r w i s e
We integrate the attribute feature matrix and similarity feature matrix to obtain the heterogeneous network HC×M as follows:
H C × M = K M c i r c R N A F C K M m i R N A F M

2.5. Node Feature Extraction Based on Sparse Autoencoder (SAE)

The features extracted from sequence and similarity often have information redundancy or “noise”. In this section, the Sparse Autoencoder (SAE) [33] is used to reconstruct the eigenmatrix. As an unsupervised autoencoder, SAE can effectively learn the hidden features of input vectors, while the introduction of a sparsity penalty term can also learn relatively sparse features well.
SAE is an unsupervised encoder including an input layer hidden layer and output layer. The input layer maps the input data X to the hidden layer Lh for encoding, where layer Lh is defined as follows:
L h = σ ( W L i X ( i ) + b L i )
where X(i) is the original input data, WL is a connection parameter between the input and hidden layers, and bLi represents an offset of function.
SAE defines σ() as the activate function, which can be represented as:
σ ( X ) = 1 ( 1 + e X )
The average activation of the activated hidden units can be calculated as:
ρ ^ h = 1 n i = 1 n [ a h ( X ( i ) ) ]
where αh() denotes the activation amount of the hidden units.
The sparsity penalty term Ps is added to the target function to keep the hidden layer at low average activation values, which are shown as:
P s = i = 1 L n K L ( ρ | | ρ ^ h )
where Ps is the sum of the degrees of penalization, ρ ^ h deviates from ρ, and Ln represents the number of units in the hidden layer. KL divergence (Kullback–Leibler) represents the sparsity penalty term of SAE and is defined as follows:
K L ( ρ | | ρ ^ h ) = ρ log ρ ρ ^ + ( 1 ρ ) log 1 ρ 1 ρ ^ h
where ρ is the sparsity parameter of KL, which is close to 0; when ρ ^ h is closer to ρ, the value of KL is smaller, and when ρ ^ h is equal to ρ, KL is equal to 0; otherwise, it increases monotonically.
With the sparsity penalty term added, the cost function is defined as:
F cos t ( W , B ) = C L ( w , b ) + δ i = 1 L n K L ( ρ | | ρ ^ )
where δ is the weight of the sparsity penalty term, and CL(w, b) is the cost function of each layer, which is calculated by the backpropagation algorithm:
w ( L ) = w ( L ) ϑ w ( L ) C cos t ( W , B )
b ( L ) = b ( L ) ϑ b ( L ) C cos t ( W , B )
where ϑ denotes the learning rate of the neural networks.
In this work, the heterogeneous network H is processed by SAE as the input data, and the final characteristic matrix DC×M is generated, where each row of DC×M represents the attribute characteristics of the corresponding node.

2.6. SGCNCMI

According to the effective application of graph neural networks in the prediction field, we propose a novel prediction model (SGCNCMI) based on a graph convolutional neural network. SGCNCMI can be described in the following six steps: (1) construct a circRNA–miRNA adjacency matrix, (2) use the RNA sequence and functional similarity to generate the node attribute feature representation, (3) use the Sparse Autoencoder (SAE) to further extract features and generate the final node feature representation, (4) apply GCN to map the relationship network diagram to a new space so as to aggregate the features of potentially associated nodes, (5) apply the weighted cross-entropy loss function to train the whole model in an end-to-end manner, and (6) apply an inner product decoder to score each pair of relationships. Next, the implementation details for each step are shown.
In step 1, we integrated known circRNA–miRNA interactions into an adjacency matrix, which contained 9905 processed high-quality interaction pairs between 2346 circRNAs and 962 miRNAs. We treated all of these 9905 interaction pairs as positive edges between circRNA nodes and miRNA nodes, and we also randomly constructed 9905 negative samples to balance the training set to better train the model. Then, all of the positive edges were labeled 1, and all of the negative samples were labeled 0.
In step 2, in order to fully express the attributes of nodes, we tried to combine multi-source information to extract node features and convert the features into digital vectors. First, related studies have confirmed that RNA molecular sequences contain abundant biological attribute information, and we applied the K-mer algorithm to process sequences to obtain the underlying feature representation. Due to the difference in the length of RNA sequences, we used 5-mers for circRNA and 2-mers for miRNA, and finally, we obtained a 128-dimension circRNA sequence vector and a 16-dimension miRNA sequence vector. Next, based on the assumption that circRNAs with similar functions are likely to be related to miRNAs with similar phenotypes, we increased two kinds of similarity (RNA Gaussian interaction profile kernel similarity and RNA sigmoid kernel similarity) to construct the comprehensive similarity matrix.
In step 3, we used SAE to further process the preliminary multidimensional features. SAE is an unsupervised autoencoder with sparsity penalty terms that can effectively extract potential features from a matrix with redundant information, while the introduction of the sparsity penalty term can obtain more valuable information from the sparse matrix. Finally, we obtained the comprehensive characteristics vectors of each node as below:
V c = ( c 1 , c 2 , , c 2346 )
V m = ( m 1 , m 2 , , m 962 )
where ci represents the features of circRNA i, and mj represents the features of miRNA j.
In step 4, we transformed the prediction of circRNA–miRNA association into a link prediction problem on a heterogeneous bipartite graph, and GCN was used to effectively learn latent graph structure information and the representations of node attributes from an end-to-end model structure. First, for an undirected heterogeneous bipartite graph A, self-connections were added to ensure nodes’ characteristic contributions:
A ^ = A + I
where A is the bipartite graph, and I is the identity matrix. In order to promote the contribution of the association relation in the propagation process of the graph convolutional network, we normalized matrix A ˜ as follows:
A ˜ = D ˜ 1 2 A ^ D ˜ 1 2
where D ˜ is calculated as:
D ˜ i i = d i a g ( j A ˜ i j )
Then, we utilized GCN containing three layers of graph convolutional networks to aggregate node features and generate a corresponding lower-dimensional feature matrix. The specific process is shown in the following formulas:
H ( l + 1 ) = σ ( A ˜ H ( l ) W ( l ) )
where H(l) represents the node feature vector of the lth layer, and H(0) is the comprehensive characteristics vector of each node that is extracted by SAE. W(l) is the lth layer trainable weight matrix, and σ() denotes the ReLU activation function. Meanwhile, to solve the problem that the contributions of different layers’ embeddings are unequal, we introduce the attention mechanism, which is defined as follows:
M c m = l n l H ( l )
where nl is the weight parameter, which is auto-learned by the graph convolutional network, and Mcm is the final embedding representation obtained by GCN. The GCN extraction process is shown in Figure 2.
In step 5, we applied weighted cross-entropy as a loss function to train the model. The loss function is defined as follows:
F l = [ b × log ( s i g m o i d ( b * ) ) × ω + ( 1 b ) × log ( 1 s i o m o i d ( b * ) ) ]
where ω represents a weight parameter, which is equal to the ratio of negative samples to positive samples. This function is used to calculate the weighted cross-entropy between the true value of the label b and the target b* obtained by the model’s internal product algorithm. Figure 2 shows the GCN processing flow.
In step 6, the inner product algorithm based on the principle of matrix factorization (MF) was used to obtain the final score of each pair, and the reconstructed score matrix can be calculated as follows:
S = s i g m o i d ( M c m M c m )
The detailed process of SGCNCMI is shown in Figure 3.
In addition, our model directly uses the similarity of a marker to disease as an attribute feature when predicting the relationship between markers and underlying diseases due to the absence of molecular sequences. This makes our model more functional and robust as a predictor of both states.

3. Results

3.1. Evaluation Criteria

Cross-validation is an important evaluation method in the field of machine learning. This section describes the performance of the model as evaluated by five-fold cross-validation experiments. In the five-fold cross-validation, we first randomly divided the samples into five subsets; in each round of the cross-validation experiment, four subsets were used to train the model, and the last subsets were treated as the test set. Meanwhile, in order to ensure the comprehensiveness and fairness of the results and verify the stability and robustness of the model, we used frequently utilized metrics to fully validate our model, which are Acc. (Accuracy), Precision, and Recall. The calculation formula is defined as:
A c c . = ( T N + T P ) ( T N + T P + F N + F P )
Pr e c i s i o n = T P ( T P + F P )
R e c a l l = T P ( T P + F N )
where TP (true positive) is the count of true samples predicted to have interacting circRNA–miRNA pairs; TN (true negative) is the number of true samples predicted to have non-interacting circRNA–miRNA pairs; FN (false negative) is the count of interacting circRNA–miRNA pairs that are predicted to have no interaction; and FP (false positives) refers to the number of non-interacting circRNA–miRNA pairs that are predicted to interact. In addition, AUC (the area under the ROC curve) and AUPR (the area under PR) were constructed to evaluate our model, and the mean value of five-fold cross-validation was used as the final score of the model.

3.2. Model Performance Evaluation

In this study, SGCNCMI was validated on the CMI-9905 dataset to evaluate the ability to predict potential circRNA–miRNA interactions. The results of the five-fold CV are recorded in Table 1. It can be seen in Table 1 that SGCNCMI achieved a mean AUC of 89.42% and a mean AUPR of 88.87%, where the AUCs of five-fold experiments were 88.41%, 89.10%, 89.57%, 89.86%, and 90.39%, and AUPR of each experiment was 87.44%, 88.27%, 89.37%, 89.71%, and 89.58%, respectively. The ROC curve and PR curve are plotted in Figure 4, which were generated by SGCNCMI using a five-fold CV.
SGCNCMI’s biomarker–disease prediction results based on the circRNA–cancer dataset is presented in Table 2, and the AUC and ACPR curves are shown in Figure 5.
SGCNCMI’s circRNA–gene prediction results based on the TransCirc dataset are presented in Table 3, and the AUC and ACPR curves are shown in Figure 6.

3.3. Discussion on the Effectiveness of GCN

The graph convolutional neural network (GCN) has been proven to be powerful for its ability to learn hidden features from an end-to-end model structure. In this work, we built a deep learning prediction model called SGCNCMI and introduced GCN into the model to aggregate the features of the relevant nodes in the network to mine hidden information for inferring circRNA–miRNA interactions.
In order to express the effectiveness of GCN concretely, in this part, we evaluated the effectiveness of GCN concerning its ability to integrate the features of associated nodes. Specifically, we compared the feature extraction based on GCN with the case in which GCN is removed. To this aim, we removed the fourth step in SGCNCMI, and after the features were extracted by SAE, we directly carried out the sixth step to obtain the final prediction score of each circRNA–miRNA interaction pair. Using the inner product based on the matrix decomposition principle, we obtained the model results without GCN aggregation characteristics, which are shown in Table 4 and Figure 7.
Figure 7 shows that the model accuracy of GCN feature extraction has been greatly improved, which proves the effectiveness of GCN as a feature extraction link in the model. In addition, it is worth noting that our model still has good predictive performance without using GCN, which indicates that our model is scientific and efficient in extracting attribute node features. Meanwhile, the strategy for removing redundancy and extracting valid information from original features through SAE has been proven in previous studies [34].

3.4. Effect of the Number of GCN Layers

GCN is a graph neural network with a certain number of layers, and the network layer of the graph neural network projects the association graph into the spectral domain to aggregate the node information in the space. The number of convolutional layers plays a crucial role in aggregating node features and extracting potential information.
As described in this section, we established GCN models with different layers, namely, one layer, two layers, three layers, four layers, or five layers, for comparative observation and recording so as to explore the influence of different layers on feature aggregation.
Table 5 and Figure 8 show the AUC and AUPR of the model with different GCN layers. From the table, it is not difficult to find that GCN with one layer achieved great performance, which demonstrates the effectiveness of GCN. Next, GCN with two layers achieved the best performance, which indicates that the first two layers of GCN can effectively extract the hidden feature information of nodes. As the number of layers increased to three or more, the performance of the model began to deteriorate significantly, which may be due to the over-smoothing of GCN; at the same time, too many GCN layers may also lead to feature redundancy and “noise”.
Therefore, although GCN can effectively extract and aggregate node information, too few or too many layers will result in less-than-optimal results. Our experiments and records also provide a reference for the use of GCN by recording parameters of different layers.

3.5. Layer Attention Mechanism Analysis

Layer attention plays an important role in controlling and quantifying the contributions of different convolutional layers. Introducing a reasonable layer attention mechanism can maximize the contribution of each layer so as to obtain the best prediction effect.
By building GCN models with different graph convolutional layers, we confirmed that each layer will have different effects on the model. In our model, the two-layer GCN model achieved the best results, which indicates that the first and second layers can effectively aggregate information. When the number of layers exceeded two, the performance of the model began to decline, and more layers often mean more redundant information, but this does not mean that these convolutional layers are not contributing. Therefore, assigning different attention weights to the convolutional layer is conducive to improving the contribution of the layers. Table 6 objectively lists the AUC of SGCNCMI with different parameters. To visually display the data, we projected the table into three-dimensional space, which is shown in Figure 9. Through the net pattern parameter, we assigned 0.7, 0.2, and 0.1 attention weights to three layers, respectively, and the model achieved the best performance.

3.6. Comparison of SGCNCMI with Other Related Models

Furthermore, in order to comprehensively prove the superiority of our model in the prediction of circRNA–miRNA interactions, we compared our model with existing models; specifically, we experimented with four models using a five-fold cross-validation method and the same dataset, and our model achieved the best effect. CMIVGSD [28] is the first calculation framework to predict circRNA–miRNA interactions, which obtains the score by using graph variational autoencoders and singular value decomposition. At present, there are few computational models about circRNA–miRNA interactions, so we also compared SGCNCMI with models in other highly relevant fields. The compared methods include DMFMDA [35], NTSHMDA [36], and AE_RF [37]. DMFMDA obtains a low-dimensional dense vector of microbes and diseases through a neural network and uses a neural network with an embedding layer for matrix factorization, and Bayesian Personalized Ranking is used to obtain the optimal model parameters. AE_RF integrates circRNA and disease similarities as features and extracts hidden biological patterns with a deep autoencoder, and the random forest classifier is used to predict the association. NTSHMDA obtains the heterogeneous network from a known microbe–disease association network by connecting the disease and microbe similarity network and uses random walk to predict human microbe–disease associations.
The specific comparison data are shown in Table 7. As shown in the table, our model results are 2% higher than those of the best model in the field of circRNA–miRNA interaction. Meanwhile, compared with models in other highly relevant fields, our model still has strong competitiveness. Without a doubt, SGCNCMI is one of the few powerful methods for predicting circRNA–miRNA interactions.

4. Case Studies

To verify the predictive ability of SGCNCMI under real conditions, we conducted a case study using 9905 circRNA–miRNA interaction pairs as a benchmark dataset. First, we used known circRNA–miRNA interaction pairs to build feature vectors and train the model. Next, the trained model was used to predict unknown interaction pairs. Finally, we ranked the final predicted scores from large to small. The top ten predicted scores are shown in Table 8. It can be seen in Table 8 that six of the top ten circRNA–miRNA interactions were confirmed in PubMed. The four unconfirmed pairs of interactions have not been confirmed by biological experiments, but the possibility of interaction between them is not ruled out.

5. Conclusions

Recently, accumulating experiments have shown that endogenous circRNAs can work as miRNA sponges, which means that circRNAs bind to miRNAs and repress their functions [38]. Predicting circRNA–miRNA interactions reveal a new mechanism for regulating miRNA activity, which will benefit the diagnosis and treatment of diseases. Predicting circRNA–miRNA interactions by the computational method can not only reduce experimental risk and cost but also provide specific ideas for biological experiments. In this work, we developed a computational model named SGCNCMI to predict potential associations based on known associations. In the model, we construct molecular signatures from a variety of angles and use SAE to extract and fuse the features. Then, based on the known association diagram, the association information of surrounding nodes is fully aggregated by a graph convolutional neural network. Finally, the predicted score is obtained through the inner product decoder. We used a variety of evaluation indicators to evaluate the predictive performance of the model, which proved that our model can effectively predict potential circRNA–miRNA interactions. At the same time, our model achieved the best results in the field of predicting circRNA–miRNA interactions, and the performance was better than the only known model at present. Our model shows promising results in predicting both circRNA–cancer and circRNA–gene associations, meaning that our model can be used not only at the molecular level but also for the diagnosis of clinical diseases and the discovery of potentially associated genes, demonstrating the power of our model.
Limited by the number and availability of datasets, the application of computational methods in the field of circRNA–miRNA interaction prediction is in its infancy, and our model is the second known calculation method. In this work, we not only carried out experiments on the data of previously published methods but also improved and added some new reliable data. In the future, we will continue to collect more comprehensive and reliable data and propose new effective computational methods with higher performance. With circRNA becoming a new hotspot in RNA research, new methods will be constantly proposed, and our model will certainly provide a reference for more reliable methods in the future.

Author Contributions

C.-Q.Y., X.-F.W., L.-P.L., Z.-H.Y. and W.-Z.H.: conceptualization, methodology, software, validation, resources and data curation. Z.-H.R., Y.-J.G. and Y.-C.L.: writing—original draft preparation. All authors contributed to manuscript revision. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Science and Technology Innovation 2030—New Generation Artificial Intelligence Major Project (No. 2018AAA0100103) and in part by the NSFC Program under Grants 61873212, 62072378, and 62002297 and the National Natural Science Foundation of China (No. 62172338).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets described in this paper can be found in Circbank http://www.circbank.cn/ (accessed on 11 July 2022) and CircR2Cancer http://www.biobdlab.cn:8000/ (accessed on 11 July 2022). The data and source code can be found at https://github.com/1axin/SGCNCMI (accessed on 11 July 2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Memczak, S.; Jens, M.; Elefsinioti, A.; Torti, F.; Krueger, J.; Rybak, A.; Maier, L.; Mackowiak, S.D.; Gregersen, L.H.; Munschauer, M.; et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 2013, 495, 333–338. [Google Scholar] [CrossRef] [PubMed]
  2. Sanger, H.L.; Klotz, G.; Riesner, D.; Gross, H.J.; Kleinschmidt, A.K. Viroids are single-stranded covalently closed circular RNA molecules existing as highly base-paired rod-like structures. Proc. Natl. Acad. Sci. USA 1976, 73, 3852–3856. [Google Scholar] [CrossRef] [PubMed]
  3. Hsu, M.T.; Coca-Prados, M. Electron microscopic evidence for the circular form of RNA in the cytoplasm of eukaryotic cells. Nature 1979, 280, 339–340. [Google Scholar] [CrossRef] [PubMed]
  4. Arnberg, A.C.; Van Ommen, G.J.; Grivell, L.A.; Van Bruggen, E.F.J.; Borst, P. Borst Some yeast mitochondrial RNAs are circular. Cell 1980, 19, 313–319. [Google Scholar] [CrossRef]
  5. Nigro, J.M.; Cho, K.R.; Fearon, E.R.; Kern, S.E.; Ruppert, J.M.; Oliner, J.D.; Kinzler, K.W.; Vogelstein, B. Scrambled exons. Cell 1991, 64, 607–613. [Google Scholar] [CrossRef]
  6. Zaphiropoulos, P.G. Circular RNAs from transcripts of the rat cytochrome P450 2C24 gene: Correlation with exon skipping. Proc. Natl. Acad. Sci. USA 1996, 93, 6536–6541. [Google Scholar] [CrossRef]
  7. Capel, B.; Swain, A.; Nicolis, S.; Hacker, A.; Walter, M.; Koopman, P.; Goodfellow, P.; Lovell-Badge, R. Circular transcripts of the testis-determining gene Sry in adult mouse testis. Cell 1993, 73, 1019–1030. [Google Scholar] [CrossRef]
  8. Danan, M.; Schwartz, S.; Edelheit, S.; Sorek, R. Transcriptome-wide discovery of circular RNAs in Archaea. Nucleic Acids Res. 2012, 40, 3131–3142. [Google Scholar] [CrossRef]
  9. Burd, C.E.; Jeck, W.R.; Liu, Y.; Sanoff, H.K.; Wang, Z.; Sharpless, N.E. Expression of linear and novel circular forms of an INK4/ARF-associated non-coding RNA correlates with atherosclerosis risk. PLoS Genet. 2010, 6, e1001233. [Google Scholar] [CrossRef]
  10. Hansen, T.B.; Jensen, T.I.; Clausen, B.H.; Bramsen, J.B.; Finsen, B.; Damgaard, C.K.; Kjems, J. Natural RNA circles function as efficient microRNA sponges. Nature 2013, 495, 384–388. [Google Scholar] [CrossRef]
  11. Rong, D.; Sun, H.; Li, Z.; Liu, S.; Dong, C.; Fu, K.; Tang, W.; Cao, H. An emerging function of circRNA-miRNAs-mRNA axis in human diseases. Oncotarget 2017, 8, 73271. [Google Scholar] [CrossRef] [PubMed]
  12. Chen, L.; Shan, G. CircRNA in cancer: Fundamental mechanism and clinical potential. Cancer Lett. 2021, 505, 49–57. [Google Scholar] [CrossRef] [PubMed]
  13. Yang, C.; Yuan, W.; Yang, X.; Li, P.; Wang, J.; Han, J.; Tao, J.; Li, P.; Yang, H.; Lv, Q.; et al. Circular RNA circ-ITCH inhibits bladder cancer progression by sponging miR-17/miR-224 and regulating p21, PTEN expression. Mol. Cancer 2018, 17, 1–12. [Google Scholar] [CrossRef] [PubMed]
  14. Su, H.; Tao, T.; Yang, Z.; Kang, X.; Zhang, X.; Kang, D.; Wu, S.; Li, C. Circular RNA cTFRC acts as the sponge of MicroRNA-107 to promote bladder carcinoma progression. Mol. Cancer 2019, 18, 1–15. [Google Scholar] [CrossRef]
  15. Luo, Z.; Rong, Z.; Zhang, J.; Zhu, Z.; Yu, Z.; Li, T.; Fu, Z.; Qiu, Z.; Huang, C. Circular RNA circCCDC9 acts as a miR-6792-3p sponge to suppress the progression of gastric cancer through regulating CAV1 expression. Mol. Cancer 2020, 19, 1–21. [Google Scholar] [CrossRef]
  16. Cen, J.; Liang, Y.; Huang, Y.; Pan, Y.; Shu, G.; Zheng, Z.; Liao, X.; Zhou, M.; Chen, D.; Fang, Y.; et al. Circular RNA circSDHC serves as a sponge for miR-127-3p to promote the proliferation and metastasis of renal cell carcinoma via the CDKN3/E2F1 axis. Mol. Cancer 2021, 20, 1–14. [Google Scholar] [CrossRef]
  17. Wang, C.C.; Li, T.H.; Huang, L.; Chen, X. Prediction of potential miRNA–disease associations based on stacked autoencoder. Brief. Bioinform. 2022, 23, bbac021. [Google Scholar] [CrossRef]
  18. Ren, Z.H.; Yu, C.Q.; Li, L.P.; You, Z.H.; Pan, J.; Guan, Y.J.; Guo, L.X. BioChemDDI: Predicting Drug–Drug Interactions by Fusing Biochemical and Structural Information through a Self-Attention Mechanism. Biology 2022, 11, 758. [Google Scholar] [CrossRef]
  19. Wang, L.; You, Z.H.; Huang, Y.A.; Huang, D.S.; Chan, K.C. An efficient approach based on multi-sources information to predict circRNA–disease associations using deep convolutional neural network. Bioinformatics 2020, 36, 4038–4046. [Google Scholar] [CrossRef]
  20. Fan, C.; Lei, X.; Fang, Z.; Jiang, Q.; Wu, F.X. CircR2Disease: A manually curated database for experimentally supported circular RNAs associated with various diseases. Database 2018, 2018, bay044. [Google Scholar] [CrossRef] [Green Version]
  21. Zhao, Z.; Wang, K.; Wu, F.; Wang, W.; Zhang, K.; Hu, H.; Liu, Y.; Jiang, T. circRNA disease: A manually curated database of experimentally supported circRNA-disease associations. Cell Death Dis. 2018, 9, 1–2. [Google Scholar] [CrossRef] [PubMed]
  22. Liu, M.; Wang, Q.; Shen, J.; Yang, B.B.; Ding, X. Circbank: A comprehensive database for circRNA with standard nomenclature. RNA Biol. 2019, 16, 899–905. [Google Scholar] [CrossRef] [PubMed]
  23. Glažar, P.; Papavasileiou, P.; Rajewsky, N. circBase: A database for circular RNAs. RNA 2014, 20, 1666–1670. [Google Scholar] [CrossRef] [PubMed]
  24. Lan, W.; Zhu, M.; Chen, Q.; Chen, B.; Liu, J.; Li, M.; Chen, Y.P.P. CircR2Cancer: A manually curated database of associations between circRNAs and cancers. Database 2020, 2020, baaa085. [Google Scholar] [CrossRef] [PubMed]
  25. Dudekula, D.B.; Panda, A.C.; Grammatikakis, I.; De, S.; Abdelmohsen, K.; Gorospe, M. CircInteractome: A web tool for exploring circular RNAs and their interacting proteins and microRNAs. RNA Biol. 2016, 13, 34–42. [Google Scholar] [CrossRef] [PubMed]
  26. Enright, A.; John, B.; Gaul, U.; Tuschl, T.; Sander, C.; Marks, D. MicroRNA targets in Drosophila. Genome Biol. 2003, 4, 1–27. [Google Scholar] [CrossRef]
  27. Lewis, B.P.; Shih, I.H.; Jones-Rhoades, M.W.; Bartel, D.P.; Burge, C.B. Prediction of mammalian microRNA targets. Cell 2003, 115, 787–798. [Google Scholar] [CrossRef]
  28. Qian, Y.; Zheng, J.; Zhang, Z.; Jiang, Y.; Zhang, J.; Deng, L. CMIVGSD: circRNA-miRNA Interaction Prediction Based on Variational Graph Auto-Encoder and Singular Value Decomposition. In Proceedings of the 2021 IEEE International Conference on Bioinformatics and Biomedicine (BIBM), Houston, TX, USA, 9–12 December 2021; IEEE: Piscataway, NJ, USA, 2022; pp. 205–210. [Google Scholar]
  29. Gao, Y.; Shang, S.; Guo, S.; Li, X.; Zhou, H.; Liu, H.; Sun, Y.; Wang, J.; Wang, P.; Zhi, H.; et al. Lnc2Cancer 3.0: An updated resource for experimentally supported lncRNA/circRNA cancer associations and web tools based on RNA-seq and scRNA-seq data. Nucleic Acids Res. 2021, 49, D1251–D1258. [Google Scholar] [CrossRef]
  30. Huang, W.; Ling, Y.; Zhang, S.; Xia, Q.; Cao, R.; Fan, X.; Fang, Z.; Wang, Z.; Zhang, G. TransCirc: An interactive database for translatable circular RNAs based on multi-omics evidence. Nucleic Acids Res. 2021, 49, D236–D242. [Google Scholar] [CrossRef]
  31. Yi, H.C.; You, Z.H.; Wang, M.N.; Guo, Z.H.; Wang, Y.B.; Zhou, J.R. RPI-SE: A stacking ensemble learning framework for ncRNA-protein interactions prediction using sequence information. BMC Bioinform. 2020, 21, 60. [Google Scholar] [CrossRef] [Green Version]
  32. Wang, L.; You, Z.H.; Chen, X.; Li, Y.M.; Dong, Y.N.; Li, L.P.; Zheng, K. LMTRDA: Using logistic model tree to predict MiRNA-disease associations by fusing multi-source information of sequences and similarities. PLoS Comput. Biol. 2019, 15, e1006865. [Google Scholar] [CrossRef] [PubMed]
  33. Andrew, N.G. Sparse autoencoder. Lect. Notes 2011, 72, 1–19. [Google Scholar]
  34. Jiang, H.J.; Huang, Y.A.; You, Z.H. SAEROF: An ensemble approach for large-scale drug-disease association prediction by incorporating rotation forest and sparse autoencoder deep neural network. Sci. Rep. 2020, 10, 4972. [Google Scholar] [CrossRef] [PubMed]
  35. Liu, Y.; Wang, S.L.; Zhang, J.F.; Zhang, W.; Zhou, S.; Li, W. Dmfmda: Prediction of microbe-disease associations based on deep matrix factorization using bayesian personalized ranking. IEEE/ACM Trans. Comput. Biol. Bioinform. 2020, 18, 1763–1772. [Google Scholar] [CrossRef]
  36. Luo, J.; Long, Y. NTSHMDA: Prediction of human microbe-disease association based on random walk by integrating network topological similarity. IEEE/ACM Trans. Comput. Biol. Bioinform. 2018, 17, 1341–1351. [Google Scholar] [CrossRef]
  37. Deepthi, K.; Jereesh, A.S. Inferring potential CircRNA–disease associations via deep autoencoder-based classification. Mol. Diagn. Ther. 2021, 25, 87–97. [Google Scholar] [CrossRef]
  38. Kulcheski, F.R.; Christoff, A.P.; Margis, R. Circular RNAs are miRNA sponges and can be used as a new class of biomarker. J. Biotechnol. 2016, 238, 42–51. [Google Scholar] [CrossRef]
Figure 1. The K-mer algorithm for sequence feature extraction.
Figure 1. The K-mer algorithm for sequence feature extraction.
Biology 11 01350 g001
Figure 2. GCN extracted features.
Figure 2. GCN extracted features.
Biology 11 01350 g002
Figure 3. The detailed process of SGCNCMI.
Figure 3. The detailed process of SGCNCMI.
Biology 11 01350 g003
Figure 4. (a,b) are the ROC and PR curves generated by the SGCNCMI based on CMI-9905 dataset, respectively.
Figure 4. (a,b) are the ROC and PR curves generated by the SGCNCMI based on CMI-9905 dataset, respectively.
Biology 11 01350 g004
Figure 5. (a,b) are the ROC and PR curves generated by the SGCNCMI based on the circRNA–cancer dataset, respectively.
Figure 5. (a,b) are the ROC and PR curves generated by the SGCNCMI based on the circRNA–cancer dataset, respectively.
Biology 11 01350 g005
Figure 6. (a,b) are the ROC and PR curves generated by the SGCNCMI based on the circRNA–gene dataset, respectively.
Figure 6. (a,b) are the ROC and PR curves generated by the SGCNCMI based on the circRNA–gene dataset, respectively.
Biology 11 01350 g006
Figure 7. The ROC curves generated by SGCNCMI without GCN.
Figure 7. The ROC curves generated by SGCNCMI without GCN.
Biology 11 01350 g007
Figure 8. AUC and AUPR of SGCNCMI with different GCN layers.
Figure 8. AUC and AUPR of SGCNCMI with different GCN layers.
Biology 11 01350 g008
Figure 9. AUC of SGCNCMI with different parameters.
Figure 9. AUC of SGCNCMI with different parameters.
Biology 11 01350 g009
Table 1. Five-fold cross-validation results based on CMI-9905 performed by SGCNCMI.
Table 1. Five-fold cross-validation results based on CMI-9905 performed by SGCNCMI.
SGCNCMIOne-FoldTwo-FoldThree-FoldFour-FoldFive-FoldMean
AUC0.88410.89100.89570.89860.90390.8942
AUPR0.87440.88270.89370.89710.89580.8887
Table 2. Five-fold cross-validation results of circRNA–cancer performed by SGCNCMI.
Table 2. Five-fold cross-validation results of circRNA–cancer performed by SGCNCMI.
SGCNCMIOne-FoldTwo-FoldThree-FoldFour-FoldFive-FoldMean
AUC0.84600.85310.84040.82870.84130.8418
AUPR0.84280.85680.84320.84780.85100.8483
Table 3. Five-fold cross-validation results of circRNA–gene performed by SGCNCMI.
Table 3. Five-fold cross-validation results of circRNA–gene performed by SGCNCMI.
SGCNCMIOne-FoldTwo-FoldThree-FoldFour-FoldFive-FoldMean
AUC0.83750.82520.81130.80460.84390.8244
AUPR0.86470.87240.83270.84360.86420.8555
Table 4. Five-fold cross-validation results obtained by SGCNCMI without GCN.
Table 4. Five-fold cross-validation results obtained by SGCNCMI without GCN.
Non-GCNOne-FoldTwo-FoldThree-FoldFour-FoldFive-FoldMean
AUC0.81400.80470.80840.80170.80530.8060
AUPR0.84540.84410.84080.82580.84120.8393
Table 5. AUC and AUPR of SGCNCMI with different GCN layers.
Table 5. AUC and AUPR of SGCNCMI with different GCN layers.
LayersOne LayerTwo LayersThree LayersFour LayersFive Layers
AUC0.86670.87330.85260.82040.8225
AUPR0.86040.87770.84720.82290.8109
Table 6. AUC of SGCNCMI under different parameters.
Table 6. AUC of SGCNCMI under different parameters.
AUCOne-FoldTwo-FoldThree-FoldFour-FoldFive-FoldMean
Weights
{0.5,0.3,0.2}89.4484.4889.2587.5486.6287.43
{0.5,0.4,0.1}89.0788.7788.6588.8887.3288.49
{0.6,0.3,0.1}88.1989.2987.8888.7789.0188.56
{0.6,0.2,0.2}89.0289.8188.7688.1588.1388.71
{0.7,0.2,0.1}88.4189.1089.5789.8690.3989.42
{0.8,0.1,0.1}87.2890.2088.4986.6186.6388.34
Table 7. Results of comparison with highly relevant models.
Table 7. Results of comparison with highly relevant models.
MethodsAE_RFDMFMDANTSHMDACMIVGSDSGCNCMI
AUC0.76620.79220.85260.88040.9015
AUPR0.82390.82300.87720.86290.9011
Table 8. The top ten prediction results in SGCNCMI based on the dataset.
Table 8. The top ten prediction results in SGCNCMI based on the dataset.
NumcircRNAmiRNAEvidence
1hsa_circ_0003998hsa-miR-326PMID:30764896
2hsa_circ_0000523hsa-miR-31PMID:30403259
3hsa_circ_0044553hsa-miR-4726-5pUnconfirmed
4hsa_circ_0000554hsa-miR-339-5pPMID:27465405
5hsa_circ_0089776hsa-miR-6752-5pUnconfirmed
6hsa_circ_0061537hsa-miR-3913-3pUnconfirmed
7hsa_circ_0010596hsa-miR-660-3pPMID:32584784
8hsa_circ_0000799hsa-miR-31-5pPMID:30103209
9hsa_circ_0068761hsa-miR-4487PMID:33534927
10hsa_circ_0079155hsa-miR-6802-3pUnconfirmed
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yu, C.-Q.; Wang, X.-F.; Li, L.-P.; You, Z.-H.; Huang, W.-Z.; Li, Y.-C.; Ren, Z.-H.; Guan, Y.-J. SGCNCMI: A New Model Combining Multi-Modal Information to Predict circRNA-Related miRNAs, Diseases and Genes. Biology 2022, 11, 1350. https://doi.org/10.3390/biology11091350

AMA Style

Yu C-Q, Wang X-F, Li L-P, You Z-H, Huang W-Z, Li Y-C, Ren Z-H, Guan Y-J. SGCNCMI: A New Model Combining Multi-Modal Information to Predict circRNA-Related miRNAs, Diseases and Genes. Biology. 2022; 11(9):1350. https://doi.org/10.3390/biology11091350

Chicago/Turabian Style

Yu, Chang-Qing, Xin-Fei Wang, Li-Ping Li, Zhu-Hong You, Wen-Zhun Huang, Yue-Chao Li, Zhong-Hao Ren, and Yong-Jian Guan. 2022. "SGCNCMI: A New Model Combining Multi-Modal Information to Predict circRNA-Related miRNAs, Diseases and Genes" Biology 11, no. 9: 1350. https://doi.org/10.3390/biology11091350

APA Style

Yu, C. -Q., Wang, X. -F., Li, L. -P., You, Z. -H., Huang, W. -Z., Li, Y. -C., Ren, Z. -H., & Guan, Y. -J. (2022). SGCNCMI: A New Model Combining Multi-Modal Information to Predict circRNA-Related miRNAs, Diseases and Genes. Biology, 11(9), 1350. https://doi.org/10.3390/biology11091350

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

Article Metrics

Back to TopTop