Research Paper Advance Articles
RWR-algorithm-based dissection of microRNA-506-3p and microRNA-140-5p as radiosensitive biomarkers in colorectal cancer
- 1 Department of Gastroenterology, Renmin Hospital of Wuhan University, Wuhan 430060, China
- 2 Department of Oncology, Renmin Hospital of Wuhan University, Wuhan 430060, China
- 3 Department of Gastroenterology, The Affiliated Hospital of Guizhou Medical University, Guiyang 550004, China
Received: May 2, 2020 Accepted: July 21, 2020 Published: October 8, 2020https://doi.org/10.18632/aging.103907
How to Cite
Copyright: © 2020 Liao et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Radiotherapy resistance is one of the main causes for treatment failure in colorectal cancer (CRC), and it is associated with the deregulation of certain microRNAs. In this study, we constructed the microRNA-mRNA network consisting of 2275 microRNAs and 7045 target genes, collected the known microRNAs related to CRC-radiosensitivity (CRCR) (n=18) as the seed nodes, and applied the algorithm of random walk with restart (RWR) to the network to identify novel CRCR-related microRNAs (n=263). In functional analysis, 263 novel microRNAs shared a high proportion of the same biological processes and pathways with the known microRNAs. In topological analysis of the sub-network of the 263 microRNAs and their targets, hsa-mir-506-3p and hsa-mir-140-5p were identified as network hub nodes. In plasma, radiosensitive patients had a higher expression level of hsa-mir-506-3p and hsa-mir-140-5p than radioresistant patients. In experimental validation, both hsa-mir-506-3p and hsa-mir-140-5p over-expression could obviously decrease the cell proliferation, survival rate and colonality in CRC cells after radiation. In conclusion, this study combined the novel network-based method with experimental validation, and identified two novel radiosensitive biomarkers of hsa-mir-506-3p and hsa-mir-140-5p in CRC.
Colorectal cancer (CRC) is one of the most common cancers worldwide, with an estimated 1.4 million new cases and 693,900 deaths per year . Currently, surgical resection remains the only curative treatment for CRC. However, about 20% to 40% CRC patients are initially diagnosed at a locally advanced, unresectable, nonmetastatic stage termed “locally advanced CRC” . Neoadjuvant chemoradiation is the standard therapeutic strategy for these patients. Unfortunately, many CRC patients are radioresistant and have a poor prognosis . Therefore, it is necessary to elucidate the underlying radioresistant mechanisms in CRC, which may ultimately improve the therapeutic outcomes.
MicroRNAs are a class of small non-coding RNAs containing about 18~24 nucleotides. MicoRNAs function as post-transcriptional regulators, and usually cause translational suppression or RNA degradation via specific binding to 3′-untranslated region (3′-UTR) of their target genes . Accumulating evidence has suggested a strong association between microRNAs deregulation and tumor radioresistance. For example, microRNA-21 expression could promote the radioresistance in a variety of cancer cell lines, including breast, lung, glioblastoma and nasopharyngeal cancers [3, 5–8]. Several microRNAs were also found dysregulated in radioresistant CRC cell lines, and upregulation or downregulation of these microRNAs could protect the cells from radiation effects. However, more novel microRNAs biomarkers are needed to illuminate the underlying mechanism of CRC radioresistance.
In recent years, studies have investigated the mechanism of several diseases with the help of advanced computational methods, such as the network-based method. Some studies adopted the gene-gene interaction network, which allowed us to identify novel disease genes based on known disease genes . In this study, we adopted a classic network-based algorithm of random walk with restart (RWR), and developed a three-stage filtration strategy to identify novel CRC-radiosensitivity (CRCR)-related microRNAs based on a microRNA-mRNA network and known CRCR-related microRNAs. For the network hub microRNAs, we also conducted an experimental validation.
CRCR-related microRNAs identified by the RWR algorithm
Eighteen microRNAs were collected by literature review, which has been experimentally validated in association with CRCR (Table 1). By integrating three microRNA databases, a microRNA-mRNA network was constructed consisting of 2275 microRNAs and 7045 target genes. Then, we selecting those 18 microRNAs as seed nodes, and conducted the RWR algorithm analysis on the network (Figure 1). After permutation test, 263 novel CRCR-related microRNAs were identified (Supplementary Table 1). There were 1386 common targets between the 18 microRNAs and 263 microRNAs (Figure 2).
Figure 1. Flowchart of the network-based method to identify radiosensitivity-related microRNAs in colorectal cancer.
Figure 2. Distribution of the targets of 18 experimentally validated microRNAs and 263 novel microRNAs predicted by the random walk with restart (RWR) algorithm. Red dot in the inner ring, experimentally validated microRNAs; green dot in the inner ring, novel microRNAs; red dot in the outer ring, common targets of the experimentally validated microRNAs and novel microRNAs; green dot in the outer ring, unique targets of the experimentally validated microRNAs or novel microRNAs; blue line, the interaction between microRNAs and the targets.
Table 1. Characteristics of radiosensitivity-related microRNAs collected from literature in colorectal cancer.
|Reported ID||Official ID||Year||Target||Radiosensitivity||PMID|
Functional enrichment analyses
Functional enrichment analyses were conducted on the targets of the 18 microRNAs and 263 novel microRNAs respectively. In biological processes, there was a high overlapped rate of 70% between the targets of the two groups (Figure 3). Most of these biological processes focused on the biosynthesis and of RNA and protein, and cell proliferation. In pathways, there was a relatively high overlapped rate of 44% between the targets of the two groups. Most of these pathways focused on the tumor related signaling pathways.
Figure 3. Functional enrichment analyses of the targets of 18 experimentally validated microRNAs and 263 novel microRNAs. Multiple biological processes and pathways were enriched in the targets of 18 experimentally validated microRNAs and 263 novel microRNAs, and there existed a high proportion of overlap.
Topological analysis of the sub-network of 263 microRNAs and their targets
The topological features of the sub-network of 263 microRNAs and their targets were evaluated by the degree, betweenness and closeness centrality. We listed the overlap of top 10 microRNAs with topological features in each dimension (Figure 4). Two microRNAs (hsa-mir-506-3p and hsa-mir-140-5p) were identified as hub nodes in the network (Table 2).
Figure 4. Distribution of the sub-network topological parameters. Topological analysis was conducted in the sub-network of 263 novel microRNAs and their targets, and the Venn diagram showed the overlap in top 10 of each dimension.
Table 2. Topological features of the sub-network of 263 novel microRNAs and the targets. Two microRNAs overlapped in top 10 of each dimension.
Hsa-mir-506-3p and hsa-mir-140-5p over-expression and CRCR
In the quantitative RT-PCR of 18 plasma samples, radiosensitive patients had a significant higher expression level of hsa-mir-506-3p (P=0.0004) and hsa-mir-140-5p (P=0.0009) than radioresistant patients (Figure 5A). Moreover, microRNA-506-3p and microRNA-140-5p showed a good diagnostic performance in CRCR (area under ROC curve (AUC) = 0.925; 0.900).
Figure 5. Validation of the association between hsa-mir-506-3p and hsa-mir-140-5p and the radiosensitivity in colorectal cancer. (A) Plasma expression levels of hsa-mir-506-3p and hsa-mir-140-5p between radiosensitive and radioresistant patients. (B) Transfection efficiency of hsa-mir-506-3p mimic and hsa-mir-140-5p mimic. (C) Survival curves of transfected cells with hsa-mir-506-3p mimic and hsa-mir-140-5p mimic, negative control and untransfected cells at 0, 2, 4 and 6Gy radiation. (D) Cell proliferation assay analysis of transfected cells with hsa-mir-506-3p mimic and hsa-mir-140-5p mimic, negative control and untransfected cells after 4Gy radiation. (E) Clone forming assay analysis of transfected cells with hsa-mir-506-3p mimic and hsa-mir-140-5p mimic, negative control and untransfected cells after 4Gy radiation. Compared with the control or negative control, *P<0.05, **P<0.01, ***P<0.001.
After microRNA mimic transfection into HT29 cells for 48h, the expression levels of hsa-mir-506-3p (P<0.001) and hsa-mir-140-5p (P<0.001) increased significantly (Figure 5B). After 0, 2, 4 and 6Gy radiation, the microRNA mimic group had a lower survival rate than the control group and the negative control group for hsa-mir-506-3p (P<0.01) and hsa-mir-140-5p (P<0.001) (Figure 5C). After 4Gy radiation, the microRNA mimic group had a lower growth rate than the control group and the negative control group for hsa-mir-506-3p (P<0.01) and hsa-mir-140-5p (P<0.001) (Figure 5D). After 4Gy radiation, the microRNA mimic group had a less number of colonies than the control group and the negative control group for hsa-mir-506-3p (P=0.0105) and hsa-mir-140-5p (P=0.0018) (Figure 5E).
In this study, we built a microRNA-mRNA network based on the experimental validation computational prediction databases. Eighteen known CRCR-related microRNAs were obtained from the literature and chose as seed nodes, and then we applied the network-based RWR algorithm to identify novel CRCR-related microRNAs. Finally, 263 novel microRNAs were identified. Several microRNAs have been validated in association with the radiosensitivity in multiple cancers. For example, has-mir-339-5p could increase the radiosensitivity of lung cancer cells by targeting phosphatases of regenerating liver-1 (PRL-1) . Hypoxia-responsive has-mir-301a-5p and has-mir-301b-5p promoted radioresistance of prostate cancer cells via down-regulating NDRG2 . Hsa-mir-99a-5p enhanced the radiation sensitivity of non-small cell lung cancer by targeting mTOR, and inhibition of the glucocorticoid receptor resulted in an enhanced hsa-mir-99a-5p-mediated radiation response in stem-like cells from human prostate cancers [12, 13].
To validate the robustness of the prediction, we adopted two different methods. First, functional analysis was conducted respectively on the targets of the 18 microRNAs and 263 microRNAs. We found a significant overlap between the two groups, especially like the biological processes of the regulation of transcription, biosynthesis, cell proliferation and programmed cell death, as well as the pathways of MAPK signaling, PI3K-AKT signaling and TGF-β receptor signaling. These results indicated a similar role of the 263 microRNAs in the radioresistance of CRC, just like the 18 microRNAs.
To further validate our results, we also extracted the sub-network of the 263 microRNAs and their targets. After topological analysis of the sub-network, we found two hub nodes of hsa-mir-506-3p and hsa-mir-140-5p. Plasma samples of 18 CRC patients were collected to investigate the association between expression levels of the two microRNAs and subsequent radiosensitivity. Finally, we found radiosensitive patients had a significantly higher expression of both hsa-mir-506-3p and hsa-mir-140-5p than radioresistant patients. Both hsa-mir-506-3p and hsa-mir-140-5p over-expression could decrease the cell proliferation, survival rate and clonality in CRC cells after radiation. Hsa-mir-506-3p has been found in association with the progression and prognosis of non-small lung cancer and esophageal squamous cell cancer [14, 15]. Hsa-mir-140-5p could attenuate chemotherapeutic drug-induced cell death by regulating autophagy through IP3K2 in human osteosarcoma cells, and hsa-mir-140-5p could also play a therapeutic role for the treatment of non-small cell lung cancer [16, 17]. These indicated the accuracy and efficacy of our prediction. However, well-designed experiments were also needed to validate our results.
As a classical method, RWR algorithm was not as complicated as the methods in some studies mainly focusing on bioinformatics analyses. Nevertheless, we made further innovations in this study. First, we constructed the microRNA network through the bridge of microRNA-mRNA interaction data. Most of previous studies focused on the application of RWR algorithm to the mRNA network, and identified mRNA or protein biomarkers. Second, our study focused on not only the identification of novel CRC-radiosensitive microRNAs, but also the experimental validation. This combined method could enhance the reliability of bioinformatics analysis, especially in the subsequent clinical practice.
In conclusion, this study combined the novel network-based method with experimental validation, and identified two novel radiosensitive biomarkers of hsa-mir-506-3p and hsa-mir-140-5p in CRC.
Materials and Methods
Random walk with restart (RWR) algorithm analysis
As a classic ranking algorithm, the RWR algorithm always simulates a random walker starting from a seed node or several seed nodes and walking on a certain network . Finally, possible novel nodes are identified and ranked from high to low probabilities. The algorithm has been adopted to search novel disease genes or other related problems [19–21].
After the seed nodes are selected, the initial probability P0 for each seed node is set as 1/N (where N was the number of seed nodes), while zero for non-seed nodes. Then, the RWR algorithm simulates a random walker that moved on certain network starting from seed nodes. For formulation, let Pi denote a vector representing the probability of each node after the ith moving procedure is complete. After each moving procedure, Pi is updated as follows:
where M is the column-wise normalized adjacency matrix of the network; r is the restart probability of returning to the seed nodes at every step (r was set to 0.8 in this study to indicate the importance of the known microRNAs). When the L1-norm of the difference between two successive vectors is less than 1×106, the vector becomes stable, the RWR algorithm stops and output Pi+1 as the final vector. Each component in this vector indicates the probability of a node being trait-related. To remove the false positives, we also conducted a permutation test (1000 iterations). A false discovery rate (FDR) of less than 0.05 is considered significant.
Identification of candidate microRNAs
In this study, we adopted the RWR algorithm in the microRNA-mRNA network to identify candidate CRCR-related microRNAs. First, we reviewed the databases of PubMed and Embase, and chose certain microRNAs as the seed nodes. The microRNAs were included if meeting the following criteria: (i) reported in association with CRCR; (ii) conducted an experimental validation; (iii) availability of the experiment details. We excluded the microRNAs without experimental validation, or only detected by bioinformatics analysis.
The microRNA-mRNA network was constructed according to the interaction data from two computational prediction databases of miRDB (http://mirdb.org/) and TargetScan (http://www.targetscan.org/), and the experimental validation database of miRTarBase (http://http://mirtarbase.cuhk.edu.cn/php/index.php). To obtain more reliable results, we only included the interactions overlapped across the databases.
Then, the sub-network of candidate microRNAs and their targets were extracted from the microRNA-mRNA network. The network degree, betweenness and closeness centrality were calculated respectively to analyze topological features of the sub-network . The overlapped microRNAs with top 10 topological features in each dimension were regarded as novel CRCR biomarkers for experimental validation.
Functional enrichment analyses
To investigate the potential function of novel microRNAs, gene ontology (GO) analyses of biological processes and pathways are performed using the online tool ToppGene (https://toppgene.cchmc.org/). A FDR of less than 0.05 is chosen as the cut-off criteria.
Plasma samples and cell culture
The plasma samples were obtained from 18 CRC patients who subsequently underwent radiotherapy at Renmin Hospital of Wuhan University, from September to October in 2018. All patients received no other therapies before radiotherapy, and were provided written informed consent before sample collection. The present study protocol was approved by the ethics committee of Renmin Hospital of Wuhan University. The radiotherapy efficacy was classified into sensitive (n=10) and resistant (n=8) groups according to response evaluation criteria in solid tumors (RECIST).
Human CRC cell line HT29 was cultured in Dulbecco’s Modified Eagle medium (DMEM, Thermo Fisher, USA) supplemented with 10% fetal bovine serum (FBS, Thermo Fisher, USA) in a humidified incubator with 5% CO2 at 37°C.
RNA extraction and quantitative RT-PCR
Total RNA was extracted using TRIzol reagent (Invitrogen, USA). The RNA concentration and purity were measured by the NanoDrop spectrophotometer (Thermo Fisher, USA). Total RNA was synthesized into first-strand cRNA using a synthesis kit (Thermo Fisher, USA). Gene expression levels were subsequently detected by the ABI Sequence Detection System (ABI, USA) using the SYBR Green Premix Ex Taq II (Takara, Japan). The relative expression levels were analyzed using the 2−ΔΔCt method. Each sample was performed in triplicate. The specific primers were as follows: hsa-mir-506-3p: 5’-ACACTCATAAGGCACCCTTC-3’ (forward): 5’-TCTACTCAGAAGGGGAGTAC-3’ (reverse); hsa-mir-140-5p: 5’-CAGTGGTTTTACCCTATGGTAG-3’ (forward): 5’-ACCATAGGGTAAAACCACTGTT-3’ (reverse); GAPDH: 5’-TATAAATTGAGCCCGCAGCC-3’ (forward): 5’-TACGACCAAATCCGTTGACTC-3’ (reverse).
MicroRNA mimic transfection
Hsa-mir-506-3p mimic, hsa-mir-140-5p mimic and their negative controls were purchased from RiboBio (Guangzhou, China). Cells were transfected with Lipofectamine 2000 (Invitrogen, USA) according to the manufacturer's protocol for 48h, and then harvested for further experiments.
Colony forming assay
After transfection, cells were irradiated at the doses of 2, 4, and 6Gy using 250kV X-rays (Philips, the Netherlands). Then, 1000 cells were incubated in 6-well plate for 2 weeks to form visualized colonies of survived cells. Finally, cells were fixed with formaldehyde, stained with 0.1% crystal violet and air-dried. The number of clones was quantified. Plating efficiency was calculated as the ratio of the visible colonies from all seeded cells. The survival fraction was determined by the percentage of irradiated cells compared with the control.
Cell proliferation assay
After transfection and irradiation, cells were seeded in 96-well plate with 1×105/well and cultured for 24, 48 and 72h. Then, 10μl CCK8 solution (Beyotime, China) was added into each well, and the absorbance at 450nm was measured after incubating for 2h.
All statistical analyses are performed using SPSS 20.0 (IBM, USA) and GraphPad Prism 7.0 (GraphPad, UAS). Data are presented as mean ± SD, and t test and one-way ANOVA are used to compare two or multiple groups. Receiver operating characteristic (ROC) curve analysis is conducted to evaluate the diagnostic ability of microRNAs. P<0.05 is considered statistically significant.
FL and WD conceived and designed the bioinformatics analysis. FL and XC conducted the experimental validation. All authors contributed to the discussion. FL, XC, PP and WD contributed to the drafting of the article. All authors have revised the manuscript critically in terms of its important intellectual content and have approved the final version to be published. FL, XC and WD are responsible for the integrity of the work as a whole.
Conflicts of Interest
The authors have declared that no conflicts of interests exist.
This work was supported by National Natural Science Foundation of China (81602117), Natural Science Foundation of Hubei Province (2018CFB119), and Guidance Foundation of Renmin Hospital of Wuhan University (RMYD2018M47).
- 1. Torre LA, Siegel RL, Ward EM, Jemal A. Global cancer incidence and mortality rates and trends—an update. Cancer Epidemiol Biomarkers Prev. 2016; 25:16–27. https://doi.org/10.1158/1055-9965.EPI-15-0578 [PubMed]
- 2. Hu JL, He GY, Lan XL, Zeng ZC, Guan J, Ding Y, Qian XL, Liao WT, Ding YQ, Liang L. Inhibition of ATG12-mediated autophagy by miR-214 enhances radiosensitivity in colorectal cancer. Oncogenesis. 2018; 7:16. https://doi.org/10.1038/s41389-018-0028-8 [PubMed]
- 3. Afshar S, Najafi R, Sedighi Pashaki A, Sharifi M, Nikzad S, Gholami MH, Khoshghadam A, Amini R, Karimi J, Saidijam M. MiR-185 enhances radiosensitivity of colorectal cancer cells by targeting IGF1R and IGF2. Biomed Pharmacother. 2018; 106:763–69. https://doi.org/10.1016/j.biopha.2018.07.002 [PubMed]
- 4. Chen K, Rajewsky N. The evolution of gene regulation by transcription factors and microRNAs. Nat Rev Genet. 2007; 8:93–103. https://doi.org/10.1038/nrg1990 [PubMed]
- 5. Liu ZL, Wang H, Liu J, Wang ZX. MicroRNA-21 (miR-21) expression promotes growth, metastasis, and chemo- or radioresistance in non-small cell lung cancer cells by targeting PTEN. Mol Cell Biochem. 2013; 372:35–45. https://doi.org/10.1007/s11010-012-1443-3 [PubMed]
- 6. Zhu H, Zhu X, Cheng G, Zhou M, Lou W. Downregulation of microRNA-21 enhances radiosensitivity in nasopharyngeal carcinoma. Exp Ther Med. 2015; 9:2185–89. https://doi.org/10.3892/etm.2015.2403 [PubMed]
- 7. Gwak HS, Kim TH, Jo GH, Kim YJ, Kwak HJ, Kim JH, Yin J, Yoo H, Lee SH, Park JB. Silencing of microRNA-21 confers radio-sensitivity through inhibition of the PI3K/AKT pathway and enhancing autophagy in Malignant glioma cell lines. PLoS One. 2012; 7:e47449. https://doi.org/10.1371/journal.pone.0047449 [PubMed]
- 8. Zhang J, Zhang C, Hu L, He Y, Shi Z, Tang S, Chen Y. Abnormal expression of miR-21 and miR-95 in cancer stem-like cells is associated with radioresistance of lung cancer. Cancer Invest. 2015; 33:165–71. https://doi.org/10.3109/07357907.2015.1019676 [PubMed]
- 9. Cai YD, Zhang Q, Zhang YH, Chen L, Huang T. Identification of genes associated with breast cancer metastasis to bone on a protein-protein interaction network with a shortest path algorithm. J Proteome Res. 2017; 16:1027–38. https://doi.org/10.1021/acs.jproteome.6b00950 [PubMed]
- 10. Wang J, Jiang M, Xia S. miR-339-5p increases radiosensitivity of lung cancer cells by targeting phosphatases of regenerating liver-1 (PRL-1). Med Sci Monit. 2018; 24:8408–16. https://doi.org/10.12659/MSM.910808 [PubMed]
- 11. Wang W, Liu M, Guan Y, Wu Q. Hypoxia-responsive mir-301a and mir-301b promote radioresistance of prostate cancer cells via downregulating NDRG2. Med Sci Monit. 2016; 22:2126–32. https://doi.org/10.12659/msm.896832 [PubMed]
- 12. Rane JK, Erb HH, Nappo G, Mann VM, Simms MS, Collins AT, Visakorpi T, Maitland NJ. Inhibition of the glucocorticoid receptor results in an enhanced miR-99a/100-mediated radiation response in stem-like cells from human prostate cancers. Oncotarget. 2016; 7:51965–80. https://doi.org/10.18632/oncotarget.10207 [PubMed]
- 13. Yin H, Ma J, Chen L, Piao S, Zhang Y, Zhang S, Ma H, Li Y, Qu Y, Wang X, Xu Q. MiR-99a enhances the radiation sensitivity of non-small cell lung cancer by targeting mTOR. Cell Physiol Biochem. 2018; 46:471–81. https://doi.org/10.1159/000488615 [PubMed]
- 14. Li SP, Su HX, Zhao D, Guan QL. Plasma miRNA-506 as a prognostic biomarker for esophageal squamous cell carcinoma. Med Sci Monit. 2016; 22:2195–201. https://doi.org/10.12659/msm.899377 [PubMed]
- 15. Guo S, Yang P, Jiang X, Li X, Wang Y, Zhang X, Sun B, Zhang Y, Jia Y. Genetic and epigenetic silencing of mircoRNA-506-3p enhances COTL1 oncogene expression to foster non-small lung cancer progression. Oncotarget. 2017; 8:644–57. https://doi.org/10.18632/oncotarget.13501 [PubMed]
- 16. Yang Y, Liu D, Xi Y, Li J, Liu B, Li J. Upregulation of miRNA-140-5p inhibits inflammatory cytokines in acute lung injury through the MyD88/NF-κB signaling pathway by targeting TLR4. Exp Ther Med. 2018; 16:3913–20. https://doi.org/10.3892/etm.2018.6692 [PubMed]
- 17. Wei R, Cao G, Deng Z, Su J, Cai L. miR-140-5p attenuates chemotherapeutic drug-induced cell death by regulating autophagy through inositol 1,4,5-trisphosphate kinase 2 (IP3k2) in human osteosarcoma cells. Biosci Rep. 2016; 36:e00392. https://doi.org/10.1042/BSR20160238 [PubMed]
- 18. Köhler S, Bauer S, Horn D, Robinson PN. Walking the interactome for prioritization of candidate disease genes. Am J Hum Genet. 2008; 82:949–58. https://doi.org/10.1016/j.ajhg.2008.02.013 [PubMed]
- 19. Zhang J, Suo Y, Liu M, Xu X. Identification of genes related to proliferative diabetic retinopathy through RWR algorithm based on protein-protein interaction network. Biochim Biophys Acta Mol Basis Dis. 2018; 1864:2369–2375. https://doi.org/10.1016/j.bbadis.2017.11.017 [PubMed]
- 20. Zhu L, Su F, Xu Y, Zou Q. Network-based method for mining novel HPV infection related genes using random walk with restart algorithm. Biochim Biophys Acta Mol Basis Dis. 2018; 1864:2376–83. https://doi.org/10.1016/j.bbadis.2017.11.021 [PubMed]
- 21. Jiang W, Zhan H, Jiao Y, Li S, Gao W. A novel lncRNA-miRNA-mRNA network analysis identified the hub lncRNA RP11-159F24.1 in the pathogenesis of papillary thyroid cancer. Cancer Med. 2018; 7:6290–98. https://doi.org/10.1002/cam4.1900 [PubMed]
- 22. McDermott JE, Diamond DL, Corley C, Rasmussen AL, Katze MG, Waters KM. Topological analysis of protein co-abundance networks identifies novel host targets important for HCV infection and pathogenesis. BMC Syst Biol. 2012; 6:28. https://doi.org/10.1186/1752-0509-6-28 [PubMed]