# A Novel Model for Predicting Associations between Diseases and LncRNA-miRNA Pairs Based on a Newly Constructed Bipartite Network.

1. IntroductionMicroRNAs (miRNAs) are endogenous small and nonencoding RNA molecules, which can regulate gene expression at the posttranscriptional level by combining the 3' untranslated regions (UTRs) of target mRNAs (UTR) and lead the translation inhibited cleavage of the target mRNAs [1]. Moreover, long-noncoding RNAs (lncRNAs), as the biggest class of noncoding RNAs with length greater than 200 nt, can also regulate gene expression at different levels including transcriptional, posttranscriptional, and epigenetic regulation. Recently, increasing studies demonstrate that lncRNAs and miRNAs play a signification role in the cell proliferation and cell differentiation [2-5] and that the interactions between lncRNAs and microRNAs may have consequences for diseases, explain disease processes, and present opportunities for new therapies [6]. For example, Dey et al. proved that lncRNA H19 would give rise to microRNAs miR-675-3p and miR-675-5p to promote skeletal muscle differentiation and regeneration [7]. Yao et al. discovered that knockdown of lncRNA XIST could exert tumor-suppressive functions in human glioblastoma stem cells by upregulating miR-152 [8]. Wang et al. demonstrated that silencing of lncRNA MALAT1 by miR-101 and miR-217 would inhibit proliferation, migration, and invasion of esophageal squamous cell carcinoma cells [9]. Zhang et al. presented that lncRNA ANRIL indicated a poor prognosis of gastric cancer and promoted tumor growth by epigenetically silencing ofmiR-99a/miR-449a [10]. You et al. found that miR-449a inhibited cell growth in lung cancer and regulated lncRNA NEAT1 [11]. Emmrich et al. discovered that lncRNAs MONC and MIR100HG would act as oncogenes in AMKL blasts [12]. Leung et al. found that miR-222 and miR-221 upregulated by Ang II were transcribed from a large transcript and knockdown of Lnc-Ang362 would decrease expression of miR-221 and miR-222 and reduce cell proliferation [13]. Zhu et al. discovered that lncRNA H19 and H19-derived miRNA-675 were significantly downregulated in the metastatic prostate cancer cell line M12 compared with the non-meta-static prostate epithelial cell line [14]. Hirata et al. found that lncRNA MALAT1 was associated with miR-205 and promoted aggressive renal cell carcinoma [15]. Zhao and Ren demonstrated that TUG1 knockdown was significantly associated with decreased cell proliferation and promoted apoptosis of breast cancer cells through the regulation of miR 9 [16].

More and more researches have indicated that lncRNA-miRNA interactions are associated with the development of complex diseases, but until now, as far as we know, no prediction models have been proposed for large-scale forecasting of the associations between diseases and LMPairs. However, some prediction models have been reported to infer the associations between diseases and miRNA-miRNA pairs [17-21]. Moreover, there are researches showing that miRNA-miRNA pairs can work cooperatively to regulate an individual gene or cohort of genes that participate in similar processes [18, 22]. Inspired by these existing state-of-the-art methods and ideas for large-scale prediction of the associations between diseases and miRNA-miRNA pairs and based on the reasonable assumption that functionally similar LMPairs tend to be associated with similar diseases, in this paper, a new model named PADLMP is proposed to predict potential associations between diseases and LMPairs. To date, it is the first computational model used to predict disease-LMPairs associations. PADLMP can predict novel disease-LMPairs associations in a large scale by combining the known lncRNA-disease, miRNA-disease, and lncRNA-miRNA associations. To evaluate the prediction performance of the proposed model, evaluation frameworks of leave-one-out cross validation (LOOCV), 2-fold, and 5-fold cross validation were adopted based on the known disease-LMPairs. A series of comparison experiments were also implemented to evaluate the influence of the number of walks on prediction performance. As a result, PADLMP achieved its best performance when the number of walks was set as 2. Specifically, PADLMP achieved value of AUCs of 0.9318, 0.9090 [+ or -] 0.0264, and 0.8950 [+ or -] 0.0027 in the LOOCV, 2-fold, and 5-fold cross validation framework, respectively. The results of the prediction show that the PADLMP model is feasible and effective in predicting broad-scale disease-LMPairs associations by considering the topology information of the known disease-LMPairs dichotomous network.

2. Materials

2.1. LncRNA-Disease Associations. Known lncRNA-disease associations were downloaded from different databases such as the lncRNA-disease database lncRNADisease [23], MNDR [24], and Lnc2Cancer [25], respectively, and then, after preprocessing (getting rid of duplicate associations), 2048 distinct experimentally confirmed lncRNA-disease associations that including 1126 lncRNAs and 356 diseases were finally obtained (see Supplementary Table 1). Then we further constructed an adjacency matrix A1 of size 1126 x 356 as the information source.

2.2. miRNA-Disease Associations. We also downloaded known disease-miRNA associations from three different databases such as the miR2Disease [26], HMDD [27], and miRCancer [28], respectively. And then, after preprocessing (getting rid of duplicate associations) and mapping these newly obtained miRNAs and diseases to databases of miRBasev21 [29] and Disease Ontology (DO) [30] separately, we finally obtained 4041 disease-miRNA associations including 438 miRNAs and 263 diseases from HMDD, 1839 disease-miRNA associations including 83 cancers and 327 miRNAs from miRCancer, and 1487 disease-miRNA associations including 107 diseases and 276 miRNAs from miR2Disease (see Supplementary Table 2).

2.3. LncRNA-miRNA Associations. In this section, we downloaded two versions (2015 Version and 2017 Version) of lncRNA-miRNA association datasets from the starBasev2.0 database [31], which provided the most comprehensive experimentally confirmed lncRNA-miRNA interactions based on large-scale CLIP-Seq data. And then, after preprocessing (including elimination of duplicate values, erroneous data, and disorganized data), 20324 lncRNA-miRNA interactions including 494 miRNAs and 1127 lncRNAs were obtained finally (see Supplementary Table 3).

3. Methods

3.1. Methods Overview. In order to predict potential novel associations between diseases and LMPairs, a new model named PADLMP is proposed, which consists of three steps (Figure 1). First, the construction of association network and data integrate. Second, the similarities for lncRNAs, diseases, miRNAs, and lncRNA-miRNA pairs are calculated based on the association network. Finally, potential associations between disease and LMPairs are inferred.

3.2. Construct the Associated Network

3.2.1. LncRNA-Disease Network, Disease-miRNA Network, and LncRNA-miRNA Network. Based on these newly obtained known lncRNA-disease associations, we constructed the lncRNA-disease bipartite network [G.sub.1] = ([V.sub.1], [E.sub.1]) according to the following steps.

Step 1. Let [V.sub.l1] be the set of newly collected 1126 lncRNAs, let [V.sub.d1] be the set of newly collected 356 diseases, and [V.sub.1] = [V.sub.l1] [union] [V.sub.d1], then we can obtain the vertex set [V.sub.1] of [G.sub.1].

Step 2. [for all][l.sub.i] [member of] [V.sub.l1], if there is [d.sub.j] [member of] [V.sub.d1] satisfying the fact that the association between [l.sub.i] and [d.sub.j] belongs to the set of newly collected 2048 lncRNA-disease associations, then we define that there is an edge between [l.sub.i] and [d.sub.j] in [G.sub.1], and by this way, we can obtain the edge set [E.sub.1] of [G.sub.1]. Obviously, [E.sub.1] is composed of these newly collected 2048 lncRNA-disease associations.

Similar to [G.sub.1], we constructed the disease-miRNA bipartite network [G.sub.2] = ([V.sub.2], [E.sub.2]) according to the following steps.

Step 1. Let [V.sub.m1] be the set of all these newly collected miRNAs, let [V.sub.d2] be the set of all these newly collected diseases, and [V.sub.2] = [V.sub.m1] [union] [V.sub.d2], then we can obtain the vertex set [V.sub.2] of [G.sub.2].

Step 2. [for all][m.sub.i] [member of] [V.sub.m1], if there is [d.sub.j] [member of] [V.sub.d2] satisfying the fact that the association between [m.sub.i] and [d.sub.j] belongs to the set of all these newly collected disease-miRNA associations, then we define that there is an edge between mi and [d.sub.j] in [G.sub.2], and by this way, we can obtain the edge set [E.sub.2] of [G.sub.2]. Obviously, [E.sub.2] is composed of all these newly collected disease-miRNA associations.

We also constructed the lncRNA-miRNA bipartite network [G.sub.3] = ([V.sub.3], [E.sub.3]) according to the following steps.

Step 1. Let [V.sub.l2] be the set of newly collected 1127 lncRNAs, let [V.sub.m2] be the set of newly collected 494 miRNAs, and [V.sub.3] = [V.sub.m2] [union] [V.sub.l2], then we can obtain the vertex set [V.sub.3] of [G.sub.3].

Step 2. [for all][l.sub.i] [member of] [V.sub.l2], if there [m.sub.j] [member of] [V.sub.m2] satisfying the fact that the association between [l.sub.i] and [m.sub.j] belongs to the set of newly collected 18286 lncRNA-miRNA associations, then we define that there is an edge between [l.sub.i] and [m.sub.j] in [G.sub.3], and by this way, we can obtain the edge set [E.sub.3] of [G.sub.3]. Obviously, [E.sub.3] is composed of these newly collected 20324 lncRNA-miRNA associations.

3.2.2. Disease-LncRNA-miRNA Network. Based on above newly constructed bipartite networks such as [G.sub.1], [G.sub.2], and [G.sub.3], we constructed a new tripartite network [G.sub.4] = ([V.sub.4], [E.sub.4]) according to the following steps.

Step 1 Let [V.sub.l'] = [V.sub.l1] [intersection] [V.sub.l2] [V.sub.m'] = [V.sub.m1] [intersection] [V.sub.m2], and [V.sub.d3] = [V.sub.d1] [intersection] [V.sub.d2]. [for all][d.sub.i] [member of] [V.sub.d3], if there are [l.sub.j] [member of] [V.sub.l'] and [m.sub.k] [member of] [V.sub.m'] satisfying the fact that the association between [d.sub.i] and [l.sub.j] belongs to [E.sub.1], the association between [d.sub.i] and [m.sub.k] belongs to [E.sub.2], and the association between [l.sub.j] and [m.sub.k] belongs to [E.sub.3] simultaneously. Then we define that there are an edge between [d.sub.i] and [l.sub.j], an edge between [d.sub.i] and [m.sub.k], and an edge between [l.sub.j] and [m.sub.k] in [G.sub.4] separately, and by this way, we can obtain the edge set [E.sub.4] of [G.sub.4].

Step 2. Let [V.sub.l] [subset or equal to] [V.sub.l'] satisfying the fact that [for all][l.sub.i] [member of] [V.sub.l] there is [d.sub.j] [member of] [V.sub.d3] satisfying the fact that the association between [d.sub.j] and [l.sub.i] belongs to [E.sub.4]. Let [V.sub.m] [subset or equal to] [V.sub.m'] satisfying the fact that [for all][m.sub.i] [member of] [V.sub.m] there is [d.sub.j] [member of] [V.sub.d3] satisfying that the association between [d.sub.j] and [m.sub.i] belongs to [E.sub.4]. Let [V.sub.4] = [V.sub.4] [union] [V.sub.m] [union] [V.sub.d3], then we can obtain the vertex set [V.sub.4] of [G.sub.4].

3.2.3. Disease-LMPairs Network. Based on above newly obtained tripartite Disease-LncRNA-miRNA network [G.sub.4], we constructed a new bipartite disease-LMPairs network G = (V, E) according to the following steps.

Step 1. [for all][l.sub.i] [member of] [V.sub.l] and [m.sub.j] [member of] [V.sub.m], let [p.sub.ij] = ([l.sub.i], [m.sub.j]) and [V.sub.p] = {[p.sub.ij]} where i [member of] [1, [absolute value of ([V.sub.l])]] and j [member of] [1, [absolute value of ([V.sub.m])]], then we define V = [V.sub.d3] [union] [V.sub.p], and by this way, we can obtain the vertex set V of G.

Step 2. [V.sub.dk] [member of] [V.sub.d3], there is [p.sub.ij] = ([l.sub.i], [m.sub.j]) [member of] [V.sub.p] satisfying

the fact that the association between [d.sub.k] and [l.sub.i] belongs to [E.sub.4], the association between [d.sub.k] and [m.sub.j] belongs to [E.sub.2], and the association between [l.sub.i] and [m.sub.j] belongs to [E.sub.3] simultaneously. Then we define that there is an edge between [d.sub.k] and [p.sub.ij] in G, and by this way, we can obtain the edge set E of G.

To make it easier to understand the construction of the network, we list in "The Meaning of Vertex and Edges in the Networks" each of the vertices, edges, and their meanings that appear in Sections 3.2.1, 3.2.2, and 3.2.3.

3.3. Calculation the Similarity of Disease

3.3.1. Calculation of the Disease Semantic Similarity (DisSemSim). Firstly, we downloaded MeSH descriptors from the National Library of Medicine and curated the names of diseases using the standard MeSH disease terms. Next, we represented the relationship of different diseases by a structure of directed acyclic graph (DAG) such as DAG(D) = (T(D), E(D)). Here, T(D) represented the node set including node D and its ancestor nodes, and E(D) denoted the edge set of corresponding direct links from a parent node to a child node, which represented the relationship between different diseases [32]. Then, based on the disease DAG, the contribution of an ancestor node d to the semantic value of disease D and the contribution of the semantic value of disease D itself can be calculated by the following two equations, respectively:

[D.sub.D](d) = l if d = D

[D.sub.D](d) = max {[DELTA] * [D.sub.D](d') | d' [member of] children of d} (1)

if d [not equal to] D,

DV(D)= [summation over (d[member of]T(D))] [D.sub.D](d), (2)

where [D.sub.D](d) represents the contribution of an ancestor node d to the semantic value of disease D, DV(D) represents the contribution of the semantic value of disease D itself, and [DELTA] is the semantic contribution decay factor with value between 0 and 1. The function of parameter [DELTA] is to guarantee that, as the distances between disease D and its ancestor disease d increase, the contribution of d to D will progressively decrease. Moreover, from the above formula (1), it is easy to see that it is also reasonable to define the contribution of D to itself as 1. In addition, according to the experimental results of some previous state-of-the-art methods [33,34], we will set the value of [DELTA] as 0.5 in this paper.

In order to measure disease semantic similarity that two diseases with more common ancestor nodes in the DAG shall have higher semantic similarity, based on the assumption, we can define the semantic similarity between two diseases [d.sub.i] and [d.sub.j] as follows:

[mathematical expression not reproducible] (3)

where T([d.sub.i]) and T([d.sub.j]) represented the node sets of the DAG of [d.sub.i] and [d.sub.j], respectively

3.3.2. Calculation of the Gaussian Interaction Profile Kernel Similarity for Diseases (GIPSim). According to the assumption that functionally similar genes tend to be associated with similar diseases, we can integrate the topologic information of known miRNA-disease association network and lncRNA-disease association network to measure the disease similarity. Moreover, in this section, we will adopt Gaussian Interaction Profile Kernel to calculate the similarity of diseases. Firstly, based on the networks such as [G.sub.1] and [G.sub.2] constructed above, we can obtain two adjacency matrices such as [Y.sub.1] (or [Y.sub.2]) as follows. For any given lncRNA [l.sub.i] (or miRNA [m.sub.i]) and disease [d.sub.j], while k takes 1 or 2, we define that

[mathematical expression not reproducible] (4)

Hence, let I[P.sub.k]([d.sub.i]) denote the ith column of matrix [Y.sub.k], then we can calculate the Gaussian Kernel Similarity between the diseases [d.sub.i] and [d.sub.j] based on their interaction profiles as follows:

[mathematical expression not reproducible] (5)

where the parameter [n.sub.k] denotes the number of diseases in [G.sub.k](k = 1, 2).

Based on formula (5), we can adopt squared root approach to calculate the Gaussian Interaction Profile Kernel Similarity for diseases as follows:

GIPSim ([d.sub.i], [d.sub.j]) = [([GIP.sub.1]([d.sub.i], [d.sub.j]) x [GIP.sub.2]([d.sub.i], [d.sub.j])).sup.1/2]. (6)

3.3.3. Calculation of the Integrated Similarity between Disease. Based on these formulas presented above, we can finally define the similarity measurement between diseases [d.sub.i] and [d.sub.j] as follows:

[mathematical expression not reproducible] (7)

3.4. Calculation of the Similarity between LncRNAs (miRNAs)

3.4.1. Calculation of the LncRNA (miRNA) Functional Similarity. For any given two lncRNAs (miRNAs) such as [l.sub.i]([m.sub.i]) and [l.sub.j]([m.sub.j]), let D[T.sub.1] = [d[t.sub.11], d[t.sub.12], ..., d[t.sub.1m]} be all diseases related to [l.sub.i]([m.sub.i]) in [G.sub.1]([G.sub.2]) and let D[T.sub.2] = [d[t.sub.21], d[t.sub.22], ..., d[t.sub.2n]} be all diseases related to [l.sub.j]([m.sub.j]) in [G.sub.1]([G.sub.2]), then we can define the functional similarity between [l.sub.i]([m.sub.i]) and [l.sub.j]([m.sub.j]) as follows (k = 1, v = l or k = 2, v = m):

[FunSim.sub.k]([v.sub.i], [v.sub.j])

= [[summation].sub.1[less than or equal to]p[less than or equal to]m] SemSims (d[t.sub.1p], D[T.sub.2]) + [[summation].sub.1[less than or equal to]p[less than or equal to]n] SemSims (d[t.sub.2p], D[T.sub.1])/m + n, (8)

where

[mathematical expression not reproducible] (9)

3.4.2. Calculation of the Gaussian Interaction Profile Kernel Similarity for IncRNAs (miRNA). For any given two lncRNAs (miRNAs) such as [l.sub.i]([m.sub.i]) and [l.sub.j]([m.sub.j]), in a similar way to the calculation of [GIP.sub.1], [GIP.sub.2] can be obtained as follows (k = 1, v = l or k = 2, v = m):

[mathematical expression not reproducible] (10)

where [IP.sub.k]([v.sub.i]) and [IP.sub.k]([v.sub.j]) are the ith row and the jth row in matrix [Y.sub.k], respectively, and [n.sub.k] is the number of lncRNAs (miRNA) in [G.sub.k].

3.4.3. Calculation of the Integrated Similarity between IncRNAs (miRNAs). Based on these formulas presented above, we can finally define the similarity measurement between lncRNAs [l.sub.i] and [l.sub.j] as follows:

[mathematical expression not reproducible] (11)

3.5. Similarity for LncRNA-miRNA Pairs (LMPairSim). Based on the bipartite disease-LMPairs network G constructed above, for any given two lncRNA-miRNA pairs [p.sub.ij] = ([l.sub.i], [m.sub.j]) and [p.sub.ab] = ([l.sub.a], [m.sub.b]), we can calculate the similarity between them according to the following three different ways:

(1) Average Approach

LMPairSim ([P.sub.ij], [P.sub.ab])

= (lncSim ([l.sub.i], [l.sub.a]) + miRSim ([m.sub.j], [m.sub.b]))/2. (12)

(2) Squared Root Approach

LMPairSim ([P.sub.ij], [P.sub.ab])

= [(lncSim ([l.sub.i], [l.sub.a]) x miRSim ([m.sub.j], [m.sub.b])).sup.1/2]. (13)

(3) Centre Distance Approach

LMPairSim ([d.sub.i], [d.sub.j]) = [square root of ([(lncSim ([l.sub.i], [l.sub.a]) - AvglncSim).sup.2] + [(miRSim ([m.sub.j], [m.sub.b]) - AvgmiRSim).sup.2], (14)

where

[mathematical expression not reproducible] (15)

3.6. Prediction of Potential Associations between Diseases and LMPairs. Inspired by the KATZ method in social networks [35], disease-gene correlation prediction [36], and lncRNA-association prediction of disease [37], we explored the PADLMP measure by developing a new computational model for predicting disease-LMPairs associations (see Figure 1). Obviously, based on the formulas (12), (13), (14), and (15), let [N.sub.d] denote the number of diseases in G, [N.sub.p] denote the number of LMPairs in G, [N.sub.l] denote the number of lncRNAs in G, and [N.sub.m] denote the number of miRNAs in G, respectively, thenwe canobtaina [N.sub.d] x [N.sub.d] dimensional matrix DisSim and [N.sub.p] x [N.sub.p] dimensional matrix PairSim. Moreover, we can construct [N.sub.p] x [N.sub.p] dimensional adjacency matrices DP as follows:

[mathematical expression not reproducible] (16)

where [d.sub.i] denotes the ith disease in G and [p.sub.j] denotes the jth LMPair in G

Hence, inspired by the approach based on KATZHMDA [38] and KATZ [35], we can construct an integrated matrix [DP.sup.*] for further predicting the potential associations between diseases and LMPairs as follows:

[mathematical expression not reproducible] (17)

Based on the integrated matrix [DP.sup.*] constructed above and letting [V.sub.p] = {[P.sub.1], [P.sub.2], ..., [P.sub.Np]}, lncRNA-miRNA pair [p.sub.i] [member of] [V.sub.p] and diseases node [d.sub.j] [member of] [V.sub.d], the probability of potential association between [p.sub.i] and [d.sub.k] can be obtained as follows:

S(i,j) = [K.summation over (l=1)][[gamma].sup.l] x [DP.sup.*l] (i,j), (18)

where the parameter K is an integer bigger than 1 and the parameter [gamma] satisfies 0 < [gamma] <1.

Additionally, according to the above formula (18), it is obvious that the ([N.sub.p] + [N.sub.d]) x ([N.sub.p] + [N.sub.d]) dimensional matrix S depicts the possibilities of all associations between diseases and LMPairs in G, and it can be further modified into the following form:

[mathematical expression not reproducible] (19)

where [S.sub.11] is [N.sub.p] x [N.sub.p] dimensional matrix, [S.sub.12] is [N.sub.p] x [N.sub.d] dimensional matrix, [S.sub.21] is [N.sub.d] x [N.sub.p] dimensional matrix, and [S.sub.22] is [N.sub.d] x [N.sub.d] dimensional matrix.

From formula (19), it is easily to know that [S.sub.12] is exactly the final prediction result matrix, which includes all of the potential associations between diseases and LMPairs in G. In addition, considering that a long walker in a sparse network maybe less meaningful, it will disrupt association prediction, so we set K to 2, 3, and 4 here. Then, final prediction result matrix could be represented by matrix DP, PairSim, and DisSim based on aforementioned equation (19).

While K = 2, there is

[S.sub.122] = [gamma] x DP + [[gamma].sup.2]

x (PairSim x DP + DP x DisSim). (20)

While K = 3, there is

[S.sub.123] = [S.sub.122] + [[gamma].sup.3] x (DP x [DP.sup.T] x DP + [PairSim.sup.2] x DP + PairSim x DP x DisSim + DP x [DisSim.sup.2]. (21)

While K = 4, there is

[mathematical expression not reproducible] (22)

4. Results

In order to estimate the prediction performance of our newly proposed model PADLMP, the leave-one-out cross validation (LOOCV) procedure was adopted based on the positive samples of disease-LMPair associations. In the LOOCV validation framework, each known disease-LMPair association is used as a test sample, and the remaining disease-LMPairs association is used as a training sample for model learning. In particular, all the disease-LMPairs without known relevance proofs will be considered as candidate samples. In the LOOCV, we can obtain the rank of each left-out testing sample relative to candidate samples, and if the test samples are with a prediction level higher than a given threshold, then it will be considered to be successfully predicted. The corresponding true positive rates (TPR, sensitivity) and false positive rates (FPR, 1 - specificity) could be obtained by setting different thresholds. Here, sensitivity measures the percentage of test samples which are predicted with a higher rank than given threshold, specificity is calculated as the percentage of negative samples ranked below a given threshold. The receiver operating characteristics (ROC) curves can be drawn by plotting TPR versus FPR by different thresholds. In order to evaluate the predictive performance of PADLMP, the areas under the ROC curve (AUC) were further calculated. 1 of the AUC value showed a perfect prediction, while 0.5 of the AUC value represented purely random performance.

From the above, we can find that there are some parameters such as K, y adopted in our prediction model PADLMP. It is obvious that these parameters are critical to the prediction performance of our model. Moreover, in Section 3.5, three different ways have been proposed to calculate the similarity for lncRNA-miRNA pairs (LMPairSim), then we need to further evaluate the performances of these three different ways also. Hence, in this section, based on the validation framework of LOOCV, we implemented a series of comparison experiments to evaluate the influence of these parameters, and the simulation results were shown in Figure 2. As a result, from Figure 2, it is easy to see that PADLMP can achieve the best prediction performance while K was set to 2. Additionally, as for other parameters [gamma], during simulations, we will set [gamma] as 0.01 based on the empirical values given by previous state-of-the-art works [37, 39-41]. Moreover, in the LOOCV, for the similarity calculation of LMPairSim, we use formulas (12), (13), and (14) in order and then select the formula that obtains the maximum AUC value. As a result, the AUC value of 0.9318, 0.9262, and 0.9247 were obtained when selecting formulas (12), (14), and (13), respectively.

Furthermore, we also compared the performance of our prediction model PADLMP with that of the RLSMDA [42], WBSMDA [39], and LRLSLDA [41] in LOOCV, since negative samples were not required in PADLMP, RLSMDA, WBSMDA, and LRLSLDA. The simulation results were shown in Figure 3. It is easy to see that PADLMP can achieve a reliable AUC of 0.9318, which is much higher than the AUC of 0.8104 and 0.9281 achieved by RLSMDA, WBSMDA, LRLSLDA, respectively, In addition, we can clearly see that the AUC value of the model LRSLDA is less than 0.5, which is obviously unreasonable. So based on prior knowledge [43], we subtract this value less than 0.5 from 1 and then we get the AUC value of LRSLDA being 0.5254.

Moreover, in order to further evaluate the prediction performance of PADLMP, the k-fold cross validation was also implemented, in which all the known disease-LMPair association samples were randomly equally divided into k parts, and k -1 parts were then used as training samples for model learning while the rest part was used as testing samples for model evaluation. Specifically, in this section, considering time complexity and costs, we would only implement 2-fold and 5-fold cross validation to evaluate the prediction performance of PADLMP. In a similar way to that of LOOCV, all the disease-LMPairs without known relevance evidences would be considered as candidate samples in the k-fold cross validation. Next, in case of the prediction performance bias caused by random division of the testing samples, we would repeat the random division of the testing samples and our simulations for 100 times, and then, the corresponding ROC curves and AUCs would be obtained in a similar way to that of LOOCV. Simulation results were shown in Table 1, and as a result, from the Table 1, it is easy to see that PADLMP can achieve the best prediction performance with average AUCs of 0.9090 and 0.8950 with Standard Deviation (STD) of 0.0264 and 0.0027 in the 2-fold and 5-fold cross validation, respectively, while setting K = 2.

From the above descriptions, it is obvious that the newly proposed model PADLMP can achieve a reliable and effective prediction performance in both LOOCV and k-fold cross validation. Therefore, we released the potential disease-LMPair associations with higher predicted relevance scores publicly (see Supplementary Table 4) and anticipated that these disease-LMPair associations may offer valuable information and clues for corresponding biological experiments and would be confirmed by experimental observations in the future.

5. Case Studies

Colon cancer is a malignant tumor that is usually found at the borders of rectum and sigmoid colon [44]. This is the third most common cancer and the third leading cause of cancer death in men and women in the United States [45]. However, patients with early colon tumors only suffer from subtle symptoms [46], which make the disease difficult to be detected. In addition, worse, it is reported that its incidence has an upward trend in recent years [47]. Therefore, there is an urgent need to predict potential miRNAs and lncRNAs associated with colon tumors. With the help of modern medicine, many miRNAs have been shown to be associated with colon tumors. For example, miRNA-145 targets the insulin receptor substrate-1 and thus inhibits the growth of colon cancer cells [48].

Moreover, as the second largest cause of cancer deaths in women, breast cancer accounts for the total number of cancers in women 22% [49, 50]. Breast cancer is caused by a variety of molecular changes, traditionally diagnosed by histopathological features such as tumor size, grading, and lymph node status [49]. Studies have shown that lncRNAs and miRNAs play important role in many biological processes and are closely related to the formation of various cancers, including breast cancer [51, 52]. In order to better diagnose and treat breast cancer, it is necessary to predict breast cancer-related lncRNA or miRNAs and to identify lncRNA and miRNA biomarkers [52].

In addition, prostate cancer is a malignant tumor derived from prostate epithelial cells [53]. There are many factors, including age, family history of disease, and race, which may increase the risk of prostate neoplasms [54]. So far, many miRNAs and lncRNAs, such as miRNA has-let-7a-5p and lncRNA XIST in the prostate, have been found to be associated with prostate tumors.

As described previously, PADLMP has been demonstrated that it can achieve a reliable and effective prediction performance. Hence, in this section, case studies about above three kinds of important cancers based on top 5% of predicted results will be implemented to show the prediction performance of PADLMP. As illustrated in Table 2, the prediction results have been verified based on the recent updates in the databases such as lncRNADisease, MNDR v2.0, starBase v2.0, HMDD, miR2Disease, and miRCancer.

In Table 2, "#" and stand for databases of lncRNAdisease and MNDR v2.0, respectively, which consist of known disease-lncRNA associations. "$" stands for starBase v2.0 database, which consists of known lncRNA-miRNA associations. "!," "&," and "+" stand for databases of HMDD, miR2Disease, and miRCancer, respectively, which consist of known disease-miRNA associations.

6. Discussion and Conclusion

Accumulating evidences show that the interaction of lncRNA-miRNAs is involved in the formation of many complex human diseases, such as breast cancer [16]; however, to our knowledge, there are no prediction models proposed for large scale forecasting the associations between diseases and LMPairs. Hence, based on the existing miRNA-disease associations, lncRNA-disease associations, lncRNA-miRNA interactions, and the assumption that genes with similar functions are often associated with similar diseases, we proposed a novel prediction model PADLMP to infer potential associations between diseases and LMPairs.

In this paper, we achieved the following contributions mainly: (1) we proposed the first computational model PADLMP for large-scale prediction of disease-LMPair associations, which can predict potential associations between diseases and lncRNA-miRNA pairs effectively. (2) We transformed the tripartite Disease-LncRNA-miRNA network into a bipartite disease-LMPair network, which greatly reduced the complexity of our prediction model. (3) Negative samples were not required in our prediction model.

However, although PADLMP is a powerful tool to infer novel associations between diseases and lncRNA-miRNA pairs, there are some limitations still existing in our method. For example, firstly, although we introduced semantic similarity for diseases and LMPairs, but the calculation of Gaussian Interaction Profile Kernel Similarity greatly relied on known disease-lncRNA associations, disease-miRNA associations, and disease-LMPairs associations. Therefore, it would cause inevitable bias towards those well-investigated diseases and LMPairs. Secondly, PADLMP could not be applied to unknown diseases and LMPairs, which were poorly investigated and had not any known associations. In the future, we will try to design new methods that do not rely on the topological information of disease-LMPair association network to solve these limitations.

The Meaning of Vertex and Edges in the Networks

[G.sub.1]: LncRNA-disease bipartite network [G.sub.2]: Disease-miRNA bipartite network [G.sub.3]: LncRNA-miRNA bipartite network [G.sub.4]: Disease-LncRNA-miRNA tripartite network G: Disease-LMPairs bipartite network [V.sub.11]: LncRNA in lncRNA-disease association [V.sub.12]: LncRNA in lncRNA-miRNA association [V.sub.d1]: Disease in lncRNA-disease association [V.sub.d2]: Disease in miRNA-disease association [V.sub.m1]: miRNA in miRNA-disease association [V.sub.m2]: Disease in lncRNA-miRNA association [V.sub.l']: [V.sub.l'] = [V.sub.l1] [intersection] [V.sub.l2] [V.sub.m']: [V.sub.m'] = [V.sub.m1] [intersection] [V.sub.m2] [V.sub.d3]: [V.sub.d3] = [V.sub.d1] [intersection] [V.sub.d2] [V.sub.l]: [V.sub.l] [subset or equal to] [V.sub.l'], [V.sub.m]: [V.sub.m] [subset or equal to] [V.sub.m'] [p.sub.ij]: LncRNA-miRNA pair ([l.sub.i], [m.sub.j]) [E.sub.1]: Edge of [G.sub.1], lncRNA [l.sub.i] associated with disease [d.sub.j] if edge <[l.sub.i], [d.sub.j]> [member of] [E.sub.1] [E.sub.2]: Edge of [G.sub.2], lncRNA [m.sub.i] associated with disease [d.sub.j] if edge <[m.sub.i], [d.sub.j]> [member of] [E.sub.2] [E.sub.3]: Edge of [G.sub.3], lncRNA [l.sub.i] associated with miRNA [m.sub.j] if edge <[l.sub.i], [m.sub.j]> [member of] [E.sub.3] [E.sub.4]: Edge of [G.sub.4], lncRNA [l.sub.i] associated with disease [d.sub.j] if edge <[l.sub.i], [d.sub.j]> [member of] [E.sub.4] or lncRNA [l.sub.i] associated with miRNA [m.sub.k] if edge <[l.sub.i], [m.sub.k]> [member of] [E.sub.4] or miRNA [m.sub.k] associated with disease [d.sub.j] if edge <[m.sub.k], [d.sub.j]> [member of] [E.sub.4] [E.sub.5]: Edge of [G.sub.4], disease [d.sub.i] associated with LMPair [p.sub.jk] if edge <[d.sub.i], [p.sub.jk]> [member of] [E.sub.5].

https://doi.org/10.1155/2018/6789089

Conflicts of Interest

The authors declare no conflicts of interest in this work.

Acknowledgments

The project is partly sponsored by the Natural Science Foundation of Hunan Province (no. 2018JJ4058 and no. 2017JJ5036), the National Natural Science Foundation of China (no. 61640210 and no. 61672447), and the CERNET Next Generation Internet Technology Innovation Project (no. NGII20160305).

Supplementary Materials

There are four supplementary tables in this manuscript. And among them, Supplementary Table 1 is a description of the lncRNA-disease associations, in which there are 2048 lncRNA-disease associations included, and the 1st column represents the lncRNAs and the 2nd column represents the diseases. Supplementary Table 2 is a description of the miRNA-disease associations, which contains the ID of diseases, the ID of miRNAs, and the associations between diseases and miRNAs. Supplementary Table 3 is a description of the lncRNA-miRNA associations, which contains 20343 lncRNA-miRNA associations and in which the first column represents the lncRNAs and the second column represents the miRNAs. Supplementary Table 4 is a description of the predictive results of associations between diseases and lncRNA-miRNA pairs while adopting PADLMP to execute prediction. (Supplementary Materials)

References

[1] E. Berezikov and E. Al, "Approaches to microrna discovery," Nature Genetics, vol. 38, p. 52, 2006.

[2] P. J. Batista and H. Y. Chang, "Long noncoding RNAs: cellular address codes in development and disease," Cell, vol. 152, no. 6, pp. 1298-1307, 2013.

[3] K. C. Wang and H. Y. Chang, "Molecular mechanisms of long noncoding RNAs," Molecular Cell, vol. 43, no. 6, pp. 904-914, 2011.

[4] C. P. Ponting, P. L. Oliver, and W. Reik, "Evolution and functions oflong noncoding RNAs," Cell, vol. 136, no. 4, pp. 629-641, 2009.

[5] J. Zheng, H. Peng, and L. Wang, "Similarities/dissimilarities analysis of protein sequences based on recurrence quantification analysis," Current Bioinformatics, vol. 10, no. 1, pp. 112-119, 2015.

[6] L. Salmena, L. Poliseno, Y. Tay, L. Kats, and P. P. Pandolfi, "A ceRNA hypothesis: the rosetta stone of a hidden RNA language?" Cell, vol. 146, no. 3, pp. 353-358, 2011.

[7] B. K. Dey, K. Pfeifer, and A. Dutta, "The H19 long noncoding RNA gives rise to microRNAs miR-675-3p and miR-675-5p to promote skeletal muscle differentiation and regeneration," Genes & Development, vol. 28, no. 5, pp. 491-501, 2014.

[8] Y. Yao, J. Ma, Y. Xue et al., "Knockdown of long non-coding RNA XIST exerts tumor-suppressive functions in human glioblastoma stem cells by up-regulating miR-152," Cancer Letters, vol. 359, no. 1, pp. 75-86, 2015.

[9] X. Wang, M. Li, Z. Wang et al., "Silencing of long noncoding rna malat1 by mir-101 and mir-217 inhibits proliferation, migration, and invasion of esophageal squamous cell carcinoma cells," The Journal of Biological Chemistry, vol. 290, no. 7, pp. 3925-3935, 2015.

[10] E.-B. Zhang, R. Kong, D.-D. Yin et al., "Long noncoding RNA ANRIL indicates a poor prognosis of gastric cancer and promotes tumor growth by epigenetically silencing of miR-99a/miR-449a," Oncotarget, vol. 5, no. 8, pp. 2276-2292, 2014.

[11] J. You, Y. Zhang, B. Liu, Y. Li, N. Fang, L. Zu et al., "Microrna-449a inhibits cell growth in lung cancer and regulates long noncoding rna nuclear enriched abundant transcript 1," Indian Journal of Cancer, vol. 51, no. 7, p. e77, 2014.

[12] S. Emmrich, A. Streltsov, F. Schmidt, V. R. Thangapandi, D. Reinhardt, and J.-H. Klusmann, "LincRNAs MONC and MIR100HG act as oncogenes in acute megakaryoblastic leukemia," Molecular Cancer, vol. 13, no. 1, article 171, 2014.

[13] A. Leung, C. Trac, W. Jin et al., "Novel long noncoding RNAs are regulated by angiotensin II in vascular smooth muscle cells," Circulation Research, vol. 113, no. 3, pp. 266-278, 2013.

[14] M. Zhu, Q. Chen, X. Liu et al., "LncRNA H19/miR-675 axis represses prostate cancer metastasis by targeting TGFBI," FEBS Journal, vol. 281, no. 16, pp. 3766-3775, 2014.

[15] H. Hirata, Y. Hinoda, V. Shahryari et al., "Long noncoding RNA MALAT1 promotes aggressive renal cell carcinoma through Ezh2 and interacts with miR-205," Cancer Research, vol. 75, no. 7, pp. 1322-1331, 2015.

[16] X. B. Zhao and G. S. Ren, "Lncrna tug1 promotes breast cancer cell proliferation via inhibiting mir-9," Cancer Biomarkers, pp. 1-8, 2016.

[17] H. Peng, C. Lan, Y. Zheng, G. Hutvagner, D. Tao, and J. Li, "Cross disease analysis of co-functional microRNA pairs on a reconstructed network of disease-gene-microRNA tripartite," BMC Bioinformatics, vol. 18, no. 1, article 193, 2017

[18] X. Li, J. Xu, Y. Li et al., "Dissection of the potential characteristic of miRNA-miRNA functional synergistic regulations," Molecular BioSystems, vol. 9, no. 2, pp. 217-224, 2013.

[19] X. Yun, C. Xu, J. Guan, Y. Ping, H. Fan, Y. Li et al., "Discovering dysfunction of multiple micrornas cooperation in disease by a conserved microrna co-expression network," PLoS ONE, vol. 7, no. 2, p. e32201, 2012.

[20] S. Yoon and G. D. Micheli, Prediction of Regulatory Modules Comprising microRNAs and Target Genes, Oxford University Press, Oxford, UK, 2005.

[21] B. Wu, C. Li, P. Zhang, Q. Yao, J. Wu, J. Han et al., "Dissection of mirna-mirna interaction in esophageal squamous cell carcinoma," PLoS ONE, vol. 8, no. 9, p. e73191, 2013.

[22] X. Lai, U. Schmitz, S. K. Gupta et al., "Computational analysis of target hub gene repression regulated by multiple and cooperative miRNAs," Nucleic Acids Research, vol. 40, no. 18, pp. 8818-8834, 2012.

[23] G. Chen, Z. Wang, D. Wang et al., "LncRNADisease: a database for long-non-coding RNA-associated diseases," Nucleic Acids Research, vol. 41, no. 1, pp. D983-D986, 2013.

[24] Y. Wang, L. Chen, B. Chen, X. Li, J. Kang, K. Fan et al., "Mammalian ncrna-disease repository: a global view of ncrna-mediated disease network," Cell Death and Disease, vol. 4, no. 8, p. e765, 2013.

[25] S. Ning, J. Zhang, P. Wang et al., "Lnc2Cancer: a manually curated database of experimentally supported lncRNAs associated with various human cancers," Nucleic Acids Research, vol. 44, no. 1, pp. D980-D985, 2016.

[26] Q. Jiang, Y. Wang, Y. Hao et al., "miR2Disease: a manually curated database for microRNA deregulation in human disease," Nucleic Acids Research, vol. 37, no. 1, pp. D98-D104, 2009.

[27] Y. Li, C. Qiu, J. Tu et al., "HMDD v2.0: a database for experimentally supported human microRNA and disease associations," Nucleic Acids Research, vol. 42, pp. D1070-D1074, 2014.

[28] B. Xie, Q. Ding, H. Han, and D. Wu, "miRCancer: a micro RNA-cancer association database constructed by text mining on literature," Bioinformatics, vol. 29, no. 5, pp. 638-644, 2013.

[29] A. Kozomara and S. Griffiths-Jones, "miRBase: annotating high confidence microRNAs using deep sequencing data," Nucleic Acids Research, vol. 42, pp. D68-D73, 2014.

[30] L. M. Schriml, C. Arze, S. Nadendla et al., "Disease ontology: a backbone for disease semantic integration," Nucleic Acids Research, vol. 40, no. 1, pp. D940-D946, 2012.

[31] J. H. Li, S. Liu, H. Zhou, L. H. Qu, and J. H. Yang, "Starbase v2.0: decoding mirna-cerna, mirna-ncrna and protein-rna interaction networks from large-scale clip-seq data," Nucleic Acids Research, vol. 42, p. D92, 2014.

[32] D. Wang, J. Wang, M. Lu, F. Song, and Q. Cui, "Inferring the human microRNA functional similarity and functional network based on microRNA-associated diseases," Bioinformatics, vol. 26, no. 13, pp. 1644-1650, 2010.

[33] X. Chen, Z. You, G. Yan, and D. Gong, "IRWRLDA: improved random walk with restart for lncRNA-disease association prediction," Oncotarget, vol. 7, no. 36, pp. 57919-57931, 2016.

[34] Y. A. Huang, X. Chen, Z. H. You, D. S. Huang, and K. C. Chan, "ILNCSIM: improved lncRNA functional similarity calculation model," Oncotarget, vol. 7, no. 18, pp. 25902-25914, 2016.

[35] L. Katz, "A new status index derived from sociometric analysis," Psychometrika, vol. 18, no. 1, pp. 39-43, 1953.

[36] X. Yang, L. Gao, X. Guo, X. Shi, H. Wu, F. Song et al., "A network based method for analysis of lncrna-disease associations and prediction of lncrnas implicated in diseases," PLoS ONE, vol. 9, no. 1, p. e87797, 2014.

[37] X. Chen, "KATZLDA: KATZ measure for the lncRNA-disease association prediction," Scientific Reports, vol. 5, Article ID 16840, 2014.

[38] X. Chen, Y. A. Huang, Z. H. You, G. Y. Yan, and X. S. Wang, "A novel approach based on katz measure to predict associations of human microbiota with non-infectious diseases," Bioinformatics, vol. 33, no. 5, pp. 733-739, 2016.

[39] X. Chen, C. C. Yan, X. Zhang et al., "WBSMDA: within and between score for MiRNA-disease association prediction," Scientific Reports, vol. 6, Article ID 21106, 2016.

[40] X. Chen, Y. C. Clarence, X. Zhang, Z. H. You, Y. A. Huang, and G. Y. Yan, "HGIMDA: Heterogeneous graph inference for miRNA-disease association prediction," Oncotarget, vol. 7, no. 40, pp. 65257-65269, 2016.

[41] X. Chen and G.-Y. Yan, "Novel human lncRNA-disease association inference based on lncRNA expression profiles," Bioinformatics, vol. 29, no. 20, pp. 2617-2624, 2013.

[42] X. Chen and G. Y. Yan, "Semi-supervised learning for potential human microrna-disease associations inference," Scientific Reports, vol. 4, p. 5501, 2014.

[43] T. Fawcett, "An introduction to ROC analysis," Pattern Recognition Letters, vol. 27, no. 8, pp. 861-874, 2006.

[44] A. I. Phipps, N. M. Lindor, M. A. Jenkins et al., "Colon and rectal cancer survival by tumor location and microsatellite instability: The colon cancer family registry," Diseases of the Colon & Rectum, vol. 56, no. 8, pp. 937-944, 2013.

[45] F. Liu, D. Yuan, Y. Wei, W. Wang, L. Yan, T. Wen et al., "Systematic review and meta-analysis of the relationship between ephxl polymorphisms and colorectal cancer risk," PLoS ONE, vol. 7, no. 8, p. e43821, 2012.

[46] S. Pita-Fernandez, S. Pertega-DIaz, B. Lopez-Calvmo et al., "Diagnostic and treatment delay, quality of life and satisfaction with care in colorectal cancer patients: a study protocol," Health and Quality of Life Outcomes, vol. 11, no. 1, article 117, 2013.

[47] V. H. Chong, M. S. Abdullah, P. U. Telisinghe, and A. Jalihal, "Colorectal cancer: incidence and trend in Brunei Darussalam," Singapore Medical Journal, vol. 50, no. 11, pp. 1085-1089, 2009.

[48] B. Shi, L. Sepp-Lorenzino, M. Prisco, P. Linsley, T. Deangelis, and R. Baserga, "Micro RNA 145 targets the insulin receptor substrate-1 and inhibits the growth of colon cancer cells," The Journal of Biological Chemistry, vol. 282, no. 45, pp. 32582-32590, 2007.

[49] H. J. Donahue and D. C. Genetos, "Genomic approaches in breast cancer research," Briefings in Functional Genomics, vol. 12, no. 5, pp. 391-396, 2013.

[50] K. Karagoz, R. Sinha, and K. Y. Arga, "Triple negative breast cancer: a multi-omics network discovery strategy for candidate targets and driving pathways," OMICS: A Journal of Integrative Biology, vol. 19, no. 2, pp. 115-130, 2015.

[51] J. Meng, P. Li, Q. Zhang, Z. Yang, and S. Fu, "A four-long non-coding rna signature in predicting breast cancer survival," Journal of Experimental and Clinical Cancer Research, vol. 33, no. 1, article 84, 2014.

[52] N. Xu, F. Wang, M. Lv, and L. Cheng, "Microarray expression profile analysis of long non-coding rnas in human breast cancer: a study of chinese women," Biomedicine Pharmacotherapy, vol. 69, no. 3, pp. 221-227, 2015.

[53] G. A. Gmyrek, M. Walburg, C. P. Webb et al., "Normal and malignant prostate epithelial cells differ in their response to hepatocyte growth factor/scatter factor," The American Journal of Pathology, vol. 159, no. 2, pp. 579-590, 2001.

[54] P. C. Walsh and A. W. Partin, "Family history facilitates the early diagnosis of prostate carcinoma," Cancer, vol. 80, no. 9, pp. 1871-1874, 1997

Shunxian Zhou, (1,2) Zhanwei Xuan, (2) Lei Wang (iD), (2) Pengyao Ping, (2) and Tingrui Pei (2)

(1) College of Software and Communication Engineering, Xiangnan University, Chenzhou 423000, China

(2) College of Information Engineering, Xiangtan University, Xiangtan 411105, China

Correspondence should be addressed to Lei Wang; wanglei@xtu.edu.cn

Received 21 December 2017; Revised 16 March 2018; Accepted 26 March 2018; Published 6 May 2018

Academic Editor: Yu Xue

Caption: Figure 1: Flowchart of PADLMP based on known miRNA-disease, lncRNA-disease, and lncRNA-miRNA association network.

Caption: Figure 2: Prediction performance of PADLMP while K takes different values in LOOCV.

Caption: Figure 3: Comparison between PADLMP and RLSMDA, WBSMDA, and LRLSLDA in LOOCV.

Table 1: Prediction performance of PADLMP while K was set to different values in the 2-fold and 5-fold cross validation, respectively. 5-fold K = 2 K = 3 K = 4 AUC 0.8950 0.8367 0.7724 STD 0.0027 0.0050 0.0109 2-fold K = 2 K = 3 K = 4 AUC 0.9090 0.8709 0.8518 STD 0.0264 0.0361 0.0441 Table 2: PADLMP was applied to three kinds of important cancer. Disease LncRNA miRNA Evidence Colon cancer MALAT1 hsa-miR-145-5p #, $,! Colon cancer MALAT1 hsa-miR-181a-5p #, $, + Colon cancer MALAT1 hsa-miR-155-5p #, $,! Colon cancer MALAT1 hsa-miR-101-3p #, $,! Colon cancer MALAT1 hsa-miR-25-3p #, $, + Colon cancer MALAT1 hsa-miR-143-3p #, $,! Colon cancer MALAT1 hsa-miR-200c-3p #, $,! Colon cancer MALAT1 hsa-miR-429 #, $, + Colon cancer MALAT1 hsa-miR-22-3p #, $,! Colon cancer MALAT1 hsa-miR-320a #, $, + Breast cancer XIST hsa-let-7b-5p #,$,! Breast cancer XIST hsa-let-7a-5p #,$,! Breast cancer XIST hsa-miR-146a-5p #,$,! Breast cancer XIST hsa-miR-27a-3p #,$,! Breast cancer XIST hsa-let-7c-5p #,$,! Breast cancer XIST hsa-miR-181b-5p #,$,! Breast cancer XIST hsa-miR-181a-5p #,$,! Breast cancer XIST hsa-miR-34a-5p #,$,! Breast cancer XIST hsa-miR-25-3p #,$,! Breast cancer XIST hsa-miR-30a-5p #,$,! Prostate cancer XIST hsa-let-7b-5p #, $, & Prostate cancer XIST hsa-miR-146a-5p *, $, & Prostate cancer XIST hsa-miR-27a-3p *, $, & Prostate cancer XIST hsa-miR-7a-5p *, $, & Prostate cancer XIST hsa-miR-30a-5p *, $, & Prostate cancer XIST hsa-miR-34a-5p *, $, & Prostate cancer XIST hsa-miR-155-5p *,$, + Prostate cancer XIST hsa-miR-124-3p *,$, + Prostate cancer XIST hsa-miR-181b-5p *, $, & Prostate cancer XIST hsa-miR-25-3p *, $, &

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Zhou, Shunxian; Xuan, Zhanwei; Wang, Lei; Ping, Pengyao; Pei, Tingrui |

Publication: | Computational and Mathematical Methods in Medicine |

Date: | Jan 1, 2018 |

Words: | 7963 |

Previous Article: | Potential Genes and Pathways of Neonatal Sepsis Based on Functional Gene Set Enrichment Analyses. |

Next Article: | Novel Signal Noise Reduction Method through Cluster Analysis, Applied to Photoplethysmography. |

Topics: |