OncoTargets and Therapy (Jun 2021)

Identification of Enzalutamide Resistance-Related circRNA-miRNA-mRNA Regulatory Networks in Patients with Prostate Cancer

  • Yu J,
  • Sun S,
  • Mao W,
  • Xu B,
  • Chen M

Journal volume & issue
Vol. Volume 14
pp. 3833 – 3848

Abstract

Read online

JunJie Yu,1,2 Si Sun,1,2 WeiPu Mao,1,2 Bin Xu,3 Ming Chen3– 5 1Surgical Research Center, Institute of Urology, School of Medicine, Southeast University, Nanjing, People’s Republic of China; 2Department of Medical College, Southeast University, Nanjing, Jiangsu, People’s Republic of China; 3Department of Urology, Affiliated Zhongda Hospital of Southeast University, Nanjing, People’s Republic of China; 4Institute of Urology, Southeastern University, Nanjing, People’s Republic of China; 5Department of Urology, Affiliated Lishui People’s Hospital of Southeast University, Nanjing, People’s Republic of ChinaCorrespondence: Ming ChenDepartment of Urology, Southeast University Zhongda Hospital, No. 87 Dingjiaqiao, Hunan Road, Gulou District, Nanjing, 210009, People’s Republic of ChinaTel/Fax +86-13913009977Email [email protected]: This study aimed to identify enzalutamide resistant-related (EnzR-related) circRNAs and to characterize and validate circRNA-miRNA-mRNA ceRNA regulatory network and corresponding prognostic signature of prostate cancer patients.Methods: We obtained circRNA expression microarray from the Gene Expression Omnibus (GEO) database and performed differential expression analysis to identify EnzR-related circRNAs using the limma package. The miRNA and mRNA expression profiling were downloaded and performed differential expression analysis, then overlapped with predicted candidates. Next, we established circRNA-miRNA-mRNA ceRNA network and PPI network utilized Cytoscape software and STRING database, respectively. In addition, univariate and Lasso Cox regression analyses were applied to generate a prognostic signature. Receiver operating characteristic (ROC) curves and Kaplan–Meier analysis were used to evaluate the reliability and sensitivity of the signature. Ultimately, we chose hsa_circ_0047641 to validate the feasibility of the EnzR-related ceRNA regulatory pathway using qRT-PCR, CCK8 and Transwell assays.Results: We identified 13 EnzR-related circRNAs and constructed a ceRNA regulatory network that contained two downregulated circRNAs (has-circ-00000919 and has-circ-0000036) and two upregulated circRNAs (has-circ-0047641 and has-circ-0068697), and their sponged 6 miRNAs and 167 targeted mRNAs. Subsequently, these targeted mRNAs were performed to implement PPI analysis and to identify 10 Hub genes. Functional enrichment analysis provided new ways to seek potential biological functions. Besides, we established a prognostic signature of PCa patients based on 8 prognostic-associated mRNAs. We confirmed that the survival rates of PCa patients with high-risk subgroup were slightly lower than those with low-risk subgroup in the TCGA dataset (p< 0.001), and ROC curves revealed that the AUC value for prognostic signature was 0.816. Finally, the functional analysis suggested that knockdown of hsa_circ_0047641 could inhibit the progression of PCa and could reverse Enz-resistance in vitro.Conclusion: We identified 13 EnzR-related circRNAs, and constructed and confirmed that EnzR-related circRNA-miRNA-mRNA ceRNA network and corresponding prognostic signature could be a useful prognostic biomarker and therapeutic target.Keywords: ceRNA, circRNA, prognostic signature, prostate cancer, The Cancer Genome Atlas, Gene Expression Omnibus

Keywords