Identification of a Prognosis-related CeRNA Network in Cholangiocarcinoma and Potentially Therapeutic Molecules Using a Bioinformatic Approach and Molecular Docking

Overview

Journal Sci Rep

Specialty Science

Date 2022 Sep 28

PMID 36171401

Authors

Xiaoling Gao

Wenhao Zhang

Yanjuan Jia

Hui Xu

Yuchen Zhu

Xiong Pei

Affiliations

Soon will be listed here.

Abstract

Cholangiocarcinoma (CCA) is a highly malignant disease with a poor prognosis, and mechanisms of initiation and development are not well characterized. It is long noncoding RNAs (lncRNAs) acting as miRNA decoys to regulate cancer-related RNAs in competing endogenous RNA (ceRNA) networks that suggest a possible molecular mechanism in CCA. The current study aims to find potential prognosis biomarkers and small molecule therapeutic targets based on the construction of a CCA prognosis-related ceRNA network. A transcriptome dataset for CCA was downloaded from the TCGA database. Differentially expressed lncRNAs (DElncRNAs), DEmiRNAs and DEmRNAs were identified based on the differential expression and a DEceRNA network was constructed using predicted miRNA-lncRNA and miRNA-mRNA interactions. Heat maps, PCA analysis, and Pathway enrichment analysis and GO enrichment analysis were conducted. The prognostic risk model and molecular docking were constructed based on identified key ceRNA networks. A DElncRNA-miRNA-mRNAs network consisting of 434 lncRNA-miRNA pairs and 284 miRNA-mRNA pairs with 200 lncRNAs, 21 miRNAs, and 245 mRNAs was constructed. There were three lncRNAs (AC090772.1, LINC00519, and THAP7-AS1) and their downstream mRNAs (MECOM, MBNL3, RCN2) screened out as prognostic factors in CAA. Three key networks (LINC00519/ hsa-mir-22/ MECOM, THAP7-AS1/hsa-mir-155/MBNL3, and THAP7-AS1/hsa-mir-155/RCN2) were identified based on binding sites prediction and survival analysis. A prognostic risk model was established with a good predictive ability (AUC = 0.66-0.83). Four anticancer small molecules, MECOM and 17-alpha-estradiol (-7.1 kcal/mol), RCN2 and emodin (-8.3 kcal/mol), RCN2 and alpha-tocopherol (-5.6 kcal/mol), and MBNL3 and 17-beta-estradiol (-7.1 kcal/mol) were identified. Based on the DEceRNA network and Kaplan-Meier survival analysis, we identified three important ceRNA networks associated with the poor prognosis of CCA. Four anti-cancer small molecules were screened out by computer-assisted drug screening as potential small molecules for the treatment of CCA. This study provides theoretical support for the development of ceRNA network-based drugs to improve the prognosis of CCA.

Citing Articles

Insights into the computer-aided drug design and discovery based on anthraquinone scaffold for cancer treatment: A systematic review.

Chua H, Moshawih S, Kifli N, Goh H, Ming L PLoS One. 2024; 19(5):e0301396.

PMID: 38776291 PMC: 11111074. DOI: 10.1371/journal.pone.0301396.

Identification and validation of a novel anoikis-related long non-coding RNA signature for pancreatic adenocarcinoma to predict the prognosis and immune response.

Jiang Y, Ye Y, Huang Y, Wu Y, Wang G, Gui Z J Cancer Res Clin Oncol. 2023; 149(16):15069-15083.

PMID: 37620430 DOI: 10.1007/s00432-023-05285-x.

References

Shaw A, Bauer T, de Marinis F, Felip E, Goto Y, Liu G . First-Line Lorlatinib or Crizotinib in Advanced -Positive Lung Cancer. N Engl J Med. 2020; 383(21):2018-2029. DOI: 10.1056/NEJMoa2027187. View

Fabris L, Sato K, Alpini G, Strazzabosco M . The Tumor Microenvironment in Cholangiocarcinoma Progression. Hepatology. 2020; 73 Suppl 1:75-85. PMC: 7714713. DOI: 10.1002/hep.31410. View

Conte F, Fiscon G, Chiara M, Colombo T, Farina L, Paci P . Role of the long non-coding RNA PVT1 in the dysregulation of the ceRNA-ceRNA network in human breast cancer. PLoS One. 2017; 12(2):e0171661. PMC: 5302781. DOI: 10.1371/journal.pone.0171661. View

Liu H, Zou Y, Zhu W, Sen-Liu , Zhang G, Ma R . lncRNA THAP7-AS1, transcriptionally activated by SP1 and post-transcriptionally stabilized by METTL3-mediated m6A modification, exerts oncogenic properties by improving CUL4B entry into the nucleus. Cell Death Differ. 2021; 29(3):627-641. PMC: 8901790. DOI: 10.1038/s41418-021-00879-9. View

Claperon A, Mergey M, Nguyen Ho-Bouldoires T, Vignjevic D, Wendum D, Chretien Y . EGF/EGFR axis contributes to the progression of cholangiocarcinoma through the induction of an epithelial-mesenchymal transition. J Hepatol. 2014; 61(2):325-32. DOI: 10.1016/j.jhep.2014.03.033. View

Li J, Miao B, Wang S, Dong W, Xu H, Si C . Hiplot: a comprehensive and easy-to-use web service for boosting publication-ready biomedical data visualization. Brief Bioinform. 2022; 23(4). DOI: 10.1093/bib/bbac261. View

Wang X, Hu K, Zhang Y, Yang C, Yao H . Comprehensive analysis of aberrantly expressed profiles of lncRNAs, miRNAs and mRNAs with associated ceRNA network in cholangiocarcinoma. Cancer Biomark. 2018; 23(4):549-559. DOI: 10.3233/CBM-181684. View

Yin Y, Yao S, Li C, Ma Y, Kang Z, Zhang J . Systematic analysis using a bioinformatics strategy identifies SFTA1P and LINC00519 as potential prognostic biomarkers for lung squamous cell carcinoma. Am J Transl Res. 2021; 13(1):168-182. PMC: 7847518. View

Pan S, Hu Y, Hu M, Xu Y, Chen M, Du C . [Corrigendum] S100A8 facilitates cholangiocarcinoma metastasis via upregulation of VEGF through TLR4/NF‑κB pathway activation. Int J Oncol. 2020; 56(4):1046. PMC: 7050978. DOI: 10.3892/ijo.2020.4977. View

10.

Statello L, Guo C, Chen L, Huarte M . Gene regulation by long non-coding RNAs and its biological functions. Nat Rev Mol Cell Biol. 2020; 22(2):96-118. PMC: 7754182. DOI: 10.1038/s41580-020-00315-9. View

11.

He Y, Xu R, Peng L, Hu X . Tumor infiltrating lymphocytes associated competitive endogenous RNA networks as predictors of outcome in hepatic carcinoma based on WGCNA analysis. PLoS One. 2021; 16(7):e0254829. PMC: 8321115. DOI: 10.1371/journal.pone.0254829. View

12.

Banales J, Cardinale V, Carpino G, Marzioni M, Andersen J, Invernizzi P . Expert consensus document: Cholangiocarcinoma: current knowledge and future perspectives consensus statement from the European Network for the Study of Cholangiocarcinoma (ENS-CCA). Nat Rev Gastroenterol Hepatol. 2016; 13(5):261-80. DOI: 10.1038/nrgastro.2016.51. View

13.

Kang Z, Guo L, Zhu Z, Qu R . Identification of prognostic factors for intrahepatic cholangiocarcinoma using long non-coding RNAs-associated ceRNA network. Cancer Cell Int. 2020; 20:315. PMC: 7364620. DOI: 10.1186/s12935-020-01388-4. View

14.

Tanaka M, Shibahara J, Ishikawa S, Ushiku T, Morikawa T, Shinozaki-Ushiku A . EVI1 expression is associated with aggressive behavior in intrahepatic cholangiocarcinoma. Virchows Arch. 2018; 474(1):39-46. DOI: 10.1007/s00428-018-2476-0. View

15.

Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z . clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation (Camb). 2021; 2(3):100141. PMC: 8454663. DOI: 10.1016/j.xinn.2021.100141. View

16.

John B, Enright A, Aravin A, Tuschl T, Sander C, Marks D . Human MicroRNA targets. PLoS Biol. 2004; 2(11):e363. PMC: 521178. DOI: 10.1371/journal.pbio.0020363. View

17.

Feng L, Wang J, Zhang J, Diao J, He L, Fu C . Comprehensive Analysis of E3 Ubiquitin Ligases Reveals Ring Finger Protein 223 as a Novel Oncogene Activated by KLF4 in Pancreatic Cancer. Front Cell Dev Biol. 2021; 9:738709. PMC: 8551701. DOI: 10.3389/fcell.2021.738709. View

18.

Fang S, Dong L, Liu L, Guo J, Zhao L, Zhang J . HERB: a high-throughput experiment- and reference-guided database of traditional Chinese medicine. Nucleic Acids Res. 2020; 49(D1):D1197-D1206. PMC: 7779036. DOI: 10.1093/nar/gkaa1063. View

19.

Sanaei M, Kavoosi F, Dehghani F . Comparative Analysis of the Effects of 17-Beta Estradiol on Proliferation, and Apoptosis in Hepatocellular Carcinoma Hep G2 and LCL-PI 11 Cell Lines. Asian Pac J Cancer Prev. 2018; 19(9):2637-2641. PMC: 6249463. DOI: 10.22034/APJCP.2018.19.9.2637. View

20.

Poliseno L, Salmena L, Zhang J, Carver B, Haveman W, Pandolfi P . A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature. 2010; 465(7301):1033-8. PMC: 3206313. DOI: 10.1038/nature09144. View