Yu-Peng Wu†Xiao-Dan Lin†Shao-Hao Chen†
Zhi-Bin Ke†Fei LinDong-Ning Chen
Xue-Yi Xue
Yong WeiQing-Shui ZhengYao-An Wen*
Ning Xu*- Department of Urology, The First Affiliated Hospital of Fujian Medical University, Fuzhou, China
Background: Prostate cancer (PCa) is one of the most common malignant tumors worldwide. Accumulating evidence has suggested that circular RNAs (circRNAs) are involved in the development and progression of various cancers, and they show great potential as novel biomarkers. However, the underlying mechanisms and specific functions of most circRNAs in PCa remain unknown. Here, we aimed to identify circRNAs with potential roles in PCa from the PCa expression profile.
Methods: We used data downloaded from the Gene Expression Omnibus to identify circRNAs that were differentially expressed between PCa samples and adjacent non-tumor samples. Relative expression levels of identified circRNAs were validated by quantitative real-time PCR. Micro (mi)RNA response elements were predicted by the CircInteractome database, and miRNA target genes were predicted by miRDB, miRTarBase, and TargetScan databases. Gene ontology (GO) enrichment analysis and pathway analysis revealed the potential biological and functional roles of these target genes. A circRNA–miRNA–mRNA interaction network was constructed by Cytoscape. The interaction between circRNAs and miRNAs in PCa was thoroughly reviewed in the PubMed. Finally, the mRNA expression of these genes was validated by the Cancer Genome Atlas (TCGA) and Gene Expression Profiling Interactive Analysis (GEPIA) databases. The expression of proteins encoded by these genes was further validated by the Human protein Atlas (HPA) database.
Results: A total of 60 circRNAs that were differentially expressed between PCa and healthy samples were screened, of which 15 were annotated. Three circRNAs (hsa_circ_0024353, hsa_circ_0085494, hsa_circ_0031408) certified the criteria were studied. The results of quantitative real-time PCR demonstrated that the expression of hsa_circ_0024353 was significantly downregulated in PC-3 cells when compared with RWPE-1 cells, while the expression of hsa_circ_0031408 and hsa_circ_0085494 was significantly upregulated in PC-3 cells when compared with RWPE-1 cells. GO and Kyoto Encyclopedia of Genes and Genomes analyses found that target genes were mainly enriched in metabolic processes and pathways involving phosphoinositide 3-kinase-Akt signaling, hypoxia-inducible factor-1 signaling, p53 signaling, and the cell cycle. A total of 11 miRNA target genes showing differential expression between PCa and healthy samples were selected, and their mRNA and protein expression were validated by GEPIA and HPA databases, respectively. Of these, PDE7B, DMRT2, and TGFBR3 were identified as potentially playing a role in PCa progression. Finally, three circRNA–miRNA–mRNA interaction axes were predicted by bioinformatics: hsa_circ_0024353–hsa-miR-940–PDE7B, hsa_circ_0024353–hsa-miR-1253–DMRT2, and hsa_circ_0085494–hsa-miR-330-3p–TGFBR3.
Conclusion: This study identified three circRNA–miRNA–mRNA interaction axes that might provide novel insights into the potential mechanisms underlying PCa development.
Introduction
Prostate cancer (PCa) is one of the most common cancers worldwide, and its incidence has gradually increased in many counties (Huang et al., 2018). For example, it has the seventh highest incidence and tenth highest mortality of all cancers in Chinese males (Zhang et al., 2017). Accumulating evidence suggests that at least 15% of PCa patients are diagnosed at an advanced stage (Cooperberg et al., 2010). Radical prostatectomy and androgen deprivation therapy (ADT) are the mainstay for patients with localized and advanced prostate cancer, respectively (Xue et al., 2018), while the standard therapy for metastatic PCa remains ADT. Nevertheless, metastatic PCa patients typically progress to the castration-resistant stage after 1 or 2 years of ADT (Josefsson et al., 2019). Therefore, as the means to cure and prevent PCa remain limited in the clinic, the identification of new biomarkers to monitor and intervene in PCa carcinogenesis is imperative.
Circular RNAs (circRNAs) are a class of endogenous non-coding RNAs that have covalent closed-loop structures with no 5′ cap or 3′ poly(A) tail. Their ring structure makes them more stable than mRNAs. They have many roles, including acting as micro (mi)RNA sponges (Thomson and Dinger, 2016), regulating alternative splicing (Feng J. et al., 2019), modulating parental genes (Li et al., 2015), sequestering proteins (Du et al., 2017), and functioning as scaffolds (Feng Z. et al., 2019). CircRNAs are proven essential elements in the initiation and progression of human disease, especially cancers (Wilusz and Sharp, 2013), and recent studies (Meng et al., 2017; Zhang et al., 2018; Verduci et al., 2019; Vo et al., 2019; Xie et al., 2019; Li S. et al., 2020) indicated that they function as potential biomarkers for cancer diagnosis and treatment.
In this study, we aimed to identify circRNAs that might be involved in PCa carcinogenesis by analyzing the expression profile of PCa. Data downloaded from the National Center of Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) were used to identify novel circRNAs that were differentially expressed in PCa. Gene ontology (GO) enrichment analysis and pathway analysis revealed the potential function of miRNA target genes. Finally, mRNA expression levels of these genes were validated by the Cancer Genome Atlas (TCGA) and Gene Expression Profiling Interactive Analysis (GEPIA) databases. The protein expression of these genes was further validated by the Human Protein Atlas (HPA) database.
Materials and Methods
CircRNA Screening
Differentially expressed circRNAs were screened and identified using the GEO2R online tool from the NCBI GEO database GSE140927. |LogFC| ≥ 1 and P < 0.05 were set as hub circRNA screened criteria.
Annotation of CircRNAs
We used circBase (Glažar et al., 2014) and circPrimer (Zhong et al., 2018) (version 1.2) software to annotate circRNAs of interest. The position, strand, genomic length, spliced length, and parental gene symbol were extracted from circBase. The composition of circRNA exons and introns was demonstrated by circPrimer. Coding sequences, untranslated regions, and start and stop codons were annotated by circPrimer.
Cell Culture
PC-3 and the human prostatic epithelial cell line RWPE-1 were purchased from the American Type Culture Collection (Manassas, VA). PC-3 cells were cultured in F-12K medium. The medium was supplemented with 10% fetal bovine serum (Invitrogen) and 5% penicillin–streptomycin (Invitrogen). RWPE-1 cells were cultured in keratinocyte serum-free medium supplemented with 50 mg/ml bovine pituitary extract, 5 ng/ml human epidermal growth factor, and 1% penicillin–streptomycin.
Quantitative Real-Time PCR (qRT-PCR)
Total RNA was isolated from cells and tissues using RNAiso Plus reagent (Takara, Kusatsu, Shiga, Japan). Reverse transcription (RT) was performed with random primers using the PrimeScriptTM RT reagent kit with gDNA Eraser (Takara) to detect the expression of hsa_circ_0024353, hsa_circ_0085494, and hsa_circ_0031408. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was amplified as an internal control. The DNA extraction kit (TIANGEN, Beijing, China) was used to extract genomic DNA from cells and tissues. PCR reactions were performed using Thermo Scientific DreamTaq PCR Master Mix (2×) (Thermo Fisher Scientific, Waltham, MA, United States). The amount of circRNAs and mRNA was quantified after the RT step using SYBR Green qPCR Master Mix (Takara) according to the manufacturer’s instructions. The primers used in this study were: 5′-TCCATCCTGCGAGCTCCTTG-3′ (F) and 5′-GCTGCATGGCACCTCTGTTC-3′ (R) for hsa_circ_0024353; 5′- GACTTGCAGGCTACGTTGAAGC-3′ (F) and 5′- CAAGTCCAGCTAGATCTGACACAAGAT-3′ (R) for hsa_circ_0085494; 5′- TCCTTTCTTGGCAACTGGAGGT-3′ (F), 5′-ACAAGTGAGGACAGCACGCA-3′ (R) for hsa_circ_0031408; and 5′-AGCCTCAAGATCATCAGCAATG-3′ (F) and 5′-ATGGACTGTGGTCATGAGTCCTT-3′ (R) for GAPDH.
Prediction of miRNA Response Elements (MREs) and miRNA Target Genes
We used the CircInteractome (Dudekula et al., 2016) database to predict MREs of circRNAs. Target genes of predicted miRNAs were predicted by miRDB (Chen and Wang, 2020)1, miRTarBase (Chou et al., 2018)2, and TargetScan (Agarwal et al., 2015)3 databases.
Construction of the CircRNA–miRNA–mRNA Interaction Network
We established a circRNA–miRNA–mRNA network to demonstrate the interactions occurring between mRNAs and circRNAs. The circRNA–miRNA potential correlation was predicted by the CircInteractome database, while miRNA-binding mRNAs were predicted by miRDB, miRTarBase, and TargetScan databases. To better understand key circRNAs, the circRNA–miRNA–mRNA interaction network was constructed using Cytoscape (Shannon et al., 2003) (version 3.6.1)4. The role of hub circRNAs was presented by Cytoscape.
GO and KEGG Analyses
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses including gene set analyses were performed using the ConsensusPathDB-human (Kamburov et al., 2013)5 database.
Validation of Target Genes by TCGA Database
The original mRNA expression profile used in the present study was downloaded from the TCGA database, and included 499 PCa tissues and 52 healthy prostatic tissues. Genes that were differentially expressed between tissues were screened for their interaction with predicted miRNA target genes.
Validation of Target Genes Expression by GEPIA and HPA Database
We used GEPIA (Tang et al., 2017)6 to validate the mRNA expression of target genes. Protein expression was validated by the HPA (Uhlen et al., 2017)7 database.
Results
Differential Expression of CircRNAs in Patients With PCa
A study flow diagram is presented in Figure 1. Expression profiles of the GSE140927 database (four PCa tissues and four paracancerous tissues of PCa) were screened, and GEO2R was used to identify differentially expressed circRNAs in patients with PCa. The top 250 RNAs, including 25 mRNAs, 165 lncRNAs, and 60 circRNAs, were obtained, and the 60 circRNAs were used for the following analysis.
Figure 1. Study flow diagram.
Selection of Screened CircRNAs
| LogFC| ≥ 1 and P < 0.05 were set as hub circRNA screened criteria. The top 15 circRNAs were identified as candidate circRNAs, including hsa_circ_0001407, hsa_circ_0050454, hsa_circ_0004529, hsa_circ_0024353, hsa_circ_0028880, hsa_circ_0044094, hsa_circ_0034546, hsa_circ_0045515, hsa_circ_0069230, hsa_circ_0025330, hsa_circ_0056631, hsa_circ_0082371, hsa_circ_0085494, hsa_circ_0069759, and hsa_circ_0031408. LogFC values of these circRNAs are shown in Table 1. We then used the circBase database to illustrate specific information about these circRNAs, including the position, strand, genomic length, spliced length, and gene symbol (Table 1). Only five circRNAs with genomic lengths less than 5 kb were included in this study: hsa_circ_0001407, hsa_circ_0024353, hsa_circ_0025330, hsa_circ_0085494, and hsa_circ_0031408. The corresponding gene symbol of each circRNA is shown in Table 1, and corresponding heatmaps are shown in Figure 2. No parental gene was found for hsa_circ_0001407, and hsa_circ_0025330 was excluded from further analysis because its spliced length exceeded 1 kb.
Table 1. Annotation of the top 15 differentially expressed circRNAs.
Figure 2. Heatmap of selected circRNAs. A total of five circRNAs were selected from four prostate cancer samples and four normal adjacent prostatic tissues samples. Red represents upregulation and green represents downregulation.
Annotation of Hub CircRNAs
The remaining three circRNAs (hsa_circ_0024353, hsa_circ_0085494, hsa_circ_0031408) were annotated by circPrimer. We next used circBase and circPrimer software to annotate these circRNAs. Corresponding gene symbols are APOC3 for hsa_circ_0024353 (Figure 3A), TMEM65 for hsa_circ_0085494 (Figure 3B), and STXBP6 for hsa_circ_0031408 (Figure 3C). hsa_circ_0024353 is formed by the circularization of exon 2 and exon 3 of the APOC3 gene, hsa_circ_0031408 is formed by the circularization of exon 3 of the STXBP6 gene, and hsa_circ_0085494 is formed by the circularization of exons 4–6 of the TMEM65 gene. mRNA and protein expression levels of these corresponding genes are shown in Figures 4A–C.
Figure 3. CircRNA annotation. The start codon, coding sequence, stop codon, parental gene symbol, position, and length are all presented in pattern graphs. (A) hsa_circ_0024353 is formed by the circularization of exon 2 and exon 3 of the APOC3 gene. (B) hsa_circ_0031408 is formed by the circularization of exon 3 of the STXBP6 gene. (C) hsa_circ_0085494 is formed by the circularization of exons 4–6 of the TMEM65 gene.
Figure 4. mRNA and protein expression of APOC3 (A), TMEM65 (B), and STXBP6 (C) genes of identified circRNAs.
Validation of the Expression of Identified CircRNAs in PCa
We used circPrimer (version 1.2) and Primer3 (version 4.1) software to design divergent primers for hub circRNAs. Primer diagrams are shown in Figures 5A–C. qRT-PCR was performed to validate the expression of identified circRNAs in RWPE-1 and PC-3 cells. The results of qRT-PCR demonstrated that the expression of hsa_circ_0024353 was significantly downregulated in PC-3 cells when compared with RWPE-1 cells, while the expression of hsa_circ_0031408 and hsa_circ_0085494 was significantly upregulated in PC-3 cells when compared with RWPE-1 cells (Figures 5A–C).
Figure 5. Pattern diagrams of divergent primers spanning back-splicing junctions and qRT-PCR validation of the expression of hsa_circ_0024353 (A), hsa_circ_0085494 (B), and hsa_circ_0031409 (C) in RWPE-1 and PC-3 cells. *P < 0.05, **P < 0.01, and ***P < 0.001.
MRE Prediction
The CircInteractome database predicted a total of 26, 10, and eight MREs for hsa_circ_0024353, hsa_circ_0031408, and hsa_circ_0085494, respectively (Table 2).
Table 2. Predicted miRNA response elements of circRNAs.
Prediction of miRNA Target Genes
Target genes of predicted miRNAs were predicted by miRDB, miRTarBase, and TargetScan databases, and interactions among these three databases were obtained. The selected genes were used in the construction of the circRNA–miRNA–mRNA interaction network.
Construction of the CircRNA–microRNA–mRNA Interaction Network
Hub circRNAs and mRNAs contain corresponding miRNA binding sites. Interactions between circRNAs and miRNAs were predicted by the Circular RNA Interactome online tool and combined with interactions identified between miRNAs and corresponding target genes. The circRNA–microRNA–mRNA interaction network was constructed using Cytoscape (version 3.6.1) software, and the roles of hub circRNAs are shown in Figures 6A,B. Previous publications describing interactions between circRNAs and miRNAs in PCa were thoroughly reviewed in the PubMed database. A total of 13 studies were identified describing 18 circRNA–miRNA interactions (Table 3).
Figure 6. The circRNA–miRNA-mRNA interaction network. (A) Overview of the entire network. (B) Specific relationship between circRNAs and interacting miRNAs in the network. Sizes of characters and circles represent the weight of circRNAs and miRNAs.
Table 3. PubMed literature review of circRNA–miRNA interactions.
GO and KEGG Analyses
GO and KEGG analyses of potential target genes of hsa_circ_0024353, hsa_circ_0031408, and hsa_circ_0085494 were performed by ConsensusPathDB.
hsa_circ_0024353
A total of 48 GO terms were shown to be enriched by ConsensusPathDB. These potential target genes were mainly enriched in the primary metabolic process, cellular metabolic process, organic substance metabolic process, and nitrogen compound metabolic process (Figure 7A). A total of 26 pathways were shown to be enriched by ConsensusPathDB. These potential target genes were mainly enriched in the phosphoinositide 3-kinase (PI3K)-Akt signaling pathway, cellular senescence, the cell cycle, the hypoxia-inducible factor (HIF)-1 signaling pathway, the T cell receptor signaling pathway, the p53 signaling pathway, and the phosphatidylinositol signaling system (Figure 7B).
Figure 7. GO and KEGG analyses of hsa_circ_0024353 (A,B), hsa_circ_0031408 (C,D), and hsa_circ_0085494 (E,F).
hsa_circ_0031408
A total of 15 GO terms were shown to be enriched by ConsensusPathDB. These potential target genes were mainly enriched in the primary metabolic process, cellular metabolic process, organic substance metabolic process, and the intracellular process (Figure 7C). A total of 16 pathways were shown to be enriched by ConsensusPathDB. These potential target genes were mainly enriched in viral carcinogenesis, purine metabolism, oocyte meiosis, and the cell cycle (Figure 7D).
hsa_circ_0085494
A total of 40 GO terms were shown to be enriched by ConsensusPathDB. These potential target genes were mainly enriched in cell communication, cellular metabolic process, and primary metabolic process (Figure 7E). A total of 17 pathways were shown to be enriched by ConsensusPathDB. These potential target genes were mainly enriched in the PI3K-Akt signaling pathway, miRNAs in cancer, and the AMP-activated protein kinase (AMPK) signaling pathway (Figure 7F).
Validation of miRNA Target Genes by TCGA Database
A total of 552 samples were downloaded from the TCGA database, including 52 normal adjacent tissues and 499 PCa samples. The clinicopathological characteristics of 367 patients with PCa from TCGA cohort were demonstrated in Table 4 (Patients with incomplete clinicopathological data were excluded). A comparison of these was used to identify a total of 1417 differentially expressed genes (| fold change| > 1.5 and P < 0.05). A total of 413 predicted differentially expressed mRNAs were also identified. Venn diagrams were used to select 11 overlapping genes, including SIAH3, EPHB1, CASKIN1, PDE7B, KIAA0408, HOXC4, SLC24A4, DMRT2, TGFBR3, PPP1R1C, and ONECUT2 (Figure 8A). The heatmap of these 11 genes was generated based on data from 552 samples downloaded from the TCGA database (Figure 8B). GO analyses of these 11 genes were also performed (Figure 8C).
Table 4. Clinicopathological characteristics of 367 patients with PCa from TCGA cohort.
Figure 8. Validation of miRNA-regulated genes by the TCGA database. (A) Interaction between predicted miRNA-regulated genes and differentially expressed genes between PCa and normal adjacent prostatic tissues selected from the TCGA database. (B) Heatmap generated from data of 552 samples downloaded from the TCGA database. (C) GO terms analyses of these 11 genes.
Validation of Expression in PCa Tissues by GEPIA and HPA Databases
The mRNA expression of the 11 differentially expressed genes was validated in the GEPIA database (Figure 9A). CASKIN1, HOXC4, and ONECUT2 were found to be upregulated in PCa compared with normal prostatic tissues. However, EPHB1, PDE7B, SLC24A4, DMRT2, TGFBR3, PPP1R1C, and KIAA0408 were downregulated in PCa tissues compared with normal prostatic tissues. SIAH3 showed no difference in expression between PCa tissues and normal prostatic tissues. To investigate the expression of proteins encoded by these genes, immunohistochemical assays performed in the HPA database were analyzed (Figure 9B). PDE7B and DMRT2 protein expression was found to be higher in PCa tissues than in normal prostatic tissues, which is inconsistent with mRNA findings from the GEPIA database. However, TGFBR3 protein expression was shown to be lower in PCa tissues than in normal prostatic tissues, which is consistent with GEPIA database results.
Figure 9. Validation of miRNA-regulated genes by GEPIA and HPA databases. (A) Differentially expressed genes identified in the GEPIA database. (B) Differentially expressed genes identified in the HPA database.
Discussion
CircRNAs are widely reported to be associated with a range of biological processes and to participate in many types of cancers (Meng et al., 2017; Bi et al., 2018). They have multiple functions, including acting as scaffolds in the assembly of protein complexes, modulating the expression of parental genes, and regulating RNA–protein interactions. One of their most important roles is functioning as miRNA sponges (Hansen et al., 2013). These are defined as competitive inhibitors of small RNAs in mammalian cells (Ebert et al., 2007), and circRNAs miRNA sponges regulate gene expression through their interactions (Han et al., 2017). Although previous studies (Chen et al., 2017; Meng et al., 2017; Kristensen et al., 2018) have reported circRNA miRNA sponge functions in many types of cancers, little is known about the role of circRNAs in PCa.
Kong et al. (2020) showed that hsa_circ_0006404 promoted PCa progression through sponging miR-29a-3p, while Cai et al. (2019) demonstrated that circHIPK3 overexpression accelerated the proliferation and invasion of PCa cells by regulating miRNA-339-3p. Additionally, Wang et al. (2019) revealed that circRNA ITCH suppressed PCa progression by upregulating HOXB13 expression through sponging miR-17-5p. In the present study, we used bioinformatics analyses to identify circRNAs that were differentially expressed between PCa tissues and corresponding normal adjacent prostatic tissues.
We screened and identified the top 250 RNAs from GSE140927 expression profiles, including 60 circRNAs. Of the top 15 circRNAs, only those with a genomic length less than 5 kb, a spliced length < 1 kb, and an identified parental gene symbol were selected for further analysis (hsa_circ_0024353, hsa_circ_0085494, and hsa_circ_0031408). The roles of these circRNAs in cancer, especially PCa, have not been explored.
CircRNAs have previously been reported to regulate gene expression by MREs (Hansen et al., 2013). We therefore used the CircInteractome database to predict 26, 10, and eight MREs for hsa_circ_0024353, hsa_circ_0031408, and hsa_circ_0085494, respectively, and used these to predict target genes by miRDB, miRTarBase, and TargetScan databases. GO and KEGG analyses showed that the predicted genes were mainly enriched in metabolic-related processes. Interactions among these databases enabled 11 genes to be selected which were mainly enriched in cell migration and enzyme-linked receptor protein signaling pathways.
We used Cytoscape to construct a circRNA–miRNA–mRNA interaction network which revealed that hsa_circ_0024353 interacts with the most dominant miRNAs such as hsa-miR-646, hsa-miR-940, hsa-miR-665, and hsa-miR-1253. Following GEPIA and HPA database validation of mRNA and protein expression of the 11 genes, PDE7B, DMRT2, and TGFBR3 were identified as novel genes regulated by the identified circRNAs and with potential roles in PCa. PDE7B is the target gene of hsa-miR-940 which can be regulated by hsa_circ_0024353; DMRT2 is the target gene of hsa-miR-1253 which can also be regulated by hsa_circ_0024353; and TGFBR3 is the target gene of hsa-miR-330-3p which can be regulated by hsa_circ_0085494.
The circRNAs identified in this study have not been reported previously. Therefore, we performed a literature review of circRNAs in PCa in the PubMed database, identifying 13 studies (Xia et al., 2018; Shen et al., 2019; Song et al., 2019; Wang et al., 2019; Wu et al., 2019; Hu and Guo, 2020; Kong et al., 2020; Li T. et al., 2020; Ou et al., 2020; Shan et al., 2020; Yan et al., 2020; Zheng et al., 2020a, b) of 18 circRNA–miRNA interactions. The miRNAs identified in this study have previously been shown to play important roles in PCa. Chen et al. (2018) found that long-non-coding (lnc)RNA GAS5 and ZFAS1 can be used as prognostic markers in PCa and that they are involved in translation targeted by miR-940 in PCa, while lncRNA FOXC2-AS1 was reported to promote proliferation and progression in PCa by targeting miR-1253/EZH2 (Chen et al., 2019). Liu et al. (2019) demonstrated that lncRNA LEF1-AS1 knockdown suppressed the initiation and development of PCa by targeting miR-330-5p. PDE7B, DMRT2, and TGFBR3 are involved in several types of tumor such as lung cancer (Liu et al., 2018), bladder cancer (Liu et al., 2013), and head and neck cancer (Fang et al., 2020). However, their roles in PCa have not been fully clarified. de Alexandre et al. (2015) suggested that PDE7B is involved in the regulation of cAMP and cyclic guanosine monophosphate levels associated with PCa. Nishan et al. (2019) demonstrated a role for DMRT2 in PCa patients after castration, while Sharifi et al. (2007) found that the downregulation of TGFBR3 is an important step in prostate tumorigenesis. Our findings can be used to predict three circRNA–miRNA–mRNA interaction axes that appear to be involved in PCa: hsa_circ_0024353–hsa-miR-940–PDE7B, hsa_circ_0024353–hsa-miR-1253–DMRT2, and hsa_circ_0085494–hsa-miR-330-3p–TGFBR3.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.
Author Contributions
Y-PW: conceptualization. X-DL, Z-BK, and S-HC: data curation. FL: formal analysis. X-YX and D-NC: investigation. Y-PW, NX, and YW: methodology. YW, Q-SZ, and X-YX: project administration. Q-SZ: visualization. Y-PW, Z-BK, X-DL, and S-HC: writing – original draft. Y-AW and NX: writing – review and editing. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Startup Fund for scientific research, Fujian Medical University (Grant No. 2017XQ2041), the Young and Middle-aged Talents Training Project of Fujian Provincial Health Commission (Grant No. 2019-ZQN-54), and the Joint Foundation of Fujian Province for Science and Technology Innovative Research Project (Grant No. 2017Y9093).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Abbreviations
ADT, androgen deprivation therapy; CDS, the coding sequence; circRNAs, circular RNAs; GEO, Gene Expression Omnibus; GO, Gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; MRE, microRNA response element; PCa, Prostate cancer; UTR, untranslated region.
Footnotes
- ^ http://mirdb.org/
- ^ http://mirtarbase.mbctu.edu.tw/php/index.php
- ^ http://www.targetscan.org/vert_72/
- ^ https://cytoscape.org/
- ^ http://cpdb.molgen.mpg.de/
- ^ http://gepia.cancer-pku.cn/
- ^ https://www.proteinatlas.org/
References
Agarwal, V., Bell, G. W., Nam, J.-W., and Bartel, D. P. (2015). Predicting effective microRNA target sites in mammalian mRNAs. eLife 4:e05005. doi: 10.7554/eLife.05005
Bi, W., Huang, J., Nie, C., Liu, B., He, G., Han, J., et al. (2018). CircRNA circRNA_102171 promotes papillary thyroid cancer progression through modulating CTNNBIP1-dependent activation of β-catenin pathway. J. Exp. Clin. Cancer Res. 37, 275–288. doi: 10.1186/s13046-018-0936-7
Cai, C., Zhi, Y., Wang, K., Zhang, P., Ji, Z., Xie, C., et al. (2019). CircHIPK3 overexpression accelerates the proliferation and invasion of prostate cancer cells through regulating miRNA-338-3p. Onco Targets Ther. 12, 3363–3372. doi: 10.2147/OTT.S196931
Chen, L., Zhang, S., Wu, J., Cui, J., Zhong, L., Zeng, L., et al. (2017). circRNA_100290 plays a role in oral cancer by functioning as a sponge of the miR-29 family. Oncogene 36, 4551–4561. doi: 10.1038/onc.2017.89
Chen, X., Yang, C., Xie, S., and Cheung, E. (2018). Long non-coding RNA GAS5 and ZFAS1 are prognostic markers involved in translation targeted by miR-940 in prostate cancer. Oncotarget 9, 1048–1062. doi: 10.18632/oncotarget.23254
Chen, Y., Gu, M., Liu, C., Wan, X., Shi, Q., Chen, Q., et al. (2019). Long noncoding RNA FOXC2-AS1 facilitates the proliferation and progression of prostate cancer via targeting miR-1253/EZH2. Gene 686, 37–42. doi: 10.1016/j.gene.2018.10.085
Chen, Y., and Wang, X. (2020). miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 48, D127–D131. doi: 10.1093/nar/gkz757
Chou, C.-H., Shrestha, S., Yang, C.-D., Chang, N.-W., Lin, Y.-L., Liao, K.-W., et al. (2018). miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 46, D296–D302. doi: 10.1093/nar/gkx1067
Cooperberg, M. R., Broering, J. M., and Carroll, P. R. (2010). Time trends and local variation in primary treatment of localized prostate cancer. J. Clin. Oncol. 28, 1117–1123. doi: 10.1200/JCO.2009.26.0133
de Alexandre, R. B., Horvath, A. D., Szarek, E., Manning, A. D., Leal, L. F., Kardauke, F., et al. (2015). Phosphodiesterase sequence variants may predispose to prostate cancer. Endocrine Relat. Cancer 22, 519–530. doi: 10.1530/erc-15-0134
Du, W. W., Zhang, C., Yang, W., Yong, T., Awan, F. M., and Yang, B. B. (2017). Identifying and Characterizing circRNA-Protein Interaction. Theranostics 7, 4183–4191. doi: 10.7150/thno.21299
Dudekula, D. B., Panda, A. C., Grammatikakis, I., De, S., Abdelmohsen, K., and Gorospe, M. (2016). CircInteractome: a web tool for exploring circular RNAs and their interacting proteins and microRNAs. RNA Biol. 13, 34–42. doi: 10.1080/15476286.2015.1128065
Ebert, M. S., Neilson, J. R., and Sharp, P. A. (2007). MicroRNA sponges: competitive inhibitors of small RNAs in mammalian cells. Nat. Methods 4, 721–726. doi: 10.1038/nmeth1079
Fang, W. Y., Kuo, Y. Z., Chang, J. Y., Hsiao, J. R., Kao, H. Y., Tsai, S. T., et al. (2020). The tumor suppressor TGFBR3 blocks lymph node metastasis in head and neck cancer. Cancers 12:1375. doi: 10.3390/cancers12061375
Feng, J., Chen, K., Dong, X., Xu, X., Jin, Y., Zhang, X., et al. (2019). Genome-wide identification of cancer-specific alternative splicing in circRNA. Mol. Cancer 18, 35–43. doi: 10.1186/s12943-019-0996-0
Feng, Z., Meng, S., Zhou, H., Xu, Z., Tang, Y., Li, P., et al. (2019). Functions and potential applications of circular RNAs in cancer stem cells. Front. Oncol. 9:500. doi: 10.3389/fonc.2019.00500
Glažar, P., Papavasileiou, P., and Rajewsky, N. (2014). circBase: a database for circular RNAs. RNA 20, 1666–1670. doi: 10.1261/rna.043687.113
Han, D., Li, J., Wang, H., Su, X., Hou, J., Gu, Y., et al. (2017). Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression. Hepatology 66, 1151–1164. doi: 10.1002/hep.29270
Hansen, T. B., Jensen, T. I., Clausen, B. H., Bramsen, J. B., Finsen, B., Damgaard, C. K., et al. (2013). Natural RNA circles function as efficient microRNA sponges. Nature 495, 384–388. doi: 10.1038/nature11993
Hu, Y., and Guo, B. (2020). Circ-MTO1 correlates with favorable prognosis and inhibits cell proliferation, invasion as well as miR-17-5p expression in prostate cancer. J. Clin. Lab. Anal. 34:e23086. doi: 10.1002/jcla.23086
Huang, E. Y., Chang, Y.-J., Huang, S.-P., Lin, V. C., Yu, C.-C., Huang, C.-Y., et al. (2018). A common regulatory variant in SLC35B4 influences the recurrence and survival of prostate cancer. J. Cell. Mol. Med. 22, 3661–3670. doi: 10.1111/jcmm.13649
Josefsson, A., Larsson, K., Freyhult, E., Damber, J.-E., and Welén, K. (2019). Gene expression alterations during development of castration-resistant prostate cancer are detected in circulating tumor cells. Cancers 12:39. doi: 10.3390/cancers12010039
Kamburov, A., Stelzl, U., Lehrach, H., and Herwig, R. (2013). The ConsensusPathDB interaction database: 2013 update. Nucleic Acids Res. 41, D793–D800. doi: 10.1093/nar/gks1055
Kong, Z., Wan, X., Lu, Y., Zhang, Y., Huang, Y., Xu, Y., et al. (2020). Circular RNA circFOXO3 promotes prostate cancer progression through sponging miR-29a-3p. J. Cell. Mol. Med. 24, 799–813. doi: 10.1111/jcmm.14791
Kristensen, L. S., Hansen, T. B., Venø, M. T., and Kjems, J. (2018). Circular RNAs in cancer: opportunities and challenges in the field. Oncogene 37, 555–565. doi: 10.1038/onc.2017.361
Li, S., Pei, Y., Wang, W., Liu, F., Zheng, K., and Zhang, X. (2020). Extracellular nanovesicles-transmitted circular RNA has_circ_0000190 suppresses osteosarcoma progression. J. Cell. Mol. Med. 24, 2202–2214. doi: 10.1111/jcmm.14877
Li, T., Sun, X., and Chen, L. (2020). Exosome circ_0044516 promotes prostate cancer cell proliferation and metastasis as a potential biomarker. J. Cell. Biochem. 121, 2118–2126. doi: 10.1002/jcb.28239
Li, Z., Huang, C., Bao, C., Chen, L., Lin, M., Wang, X., et al. (2015). Exon-intron circular RNAs regulate transcription in the nucleus. Nat. Struct. Mol. Biol. 22, 256–264. doi: 10.1038/nsmb.2959
Liu, C., Yang, Z., Deng, Z., Zhou, Y., Gong, Q., Zhao, R., et al. (2018). Upregulated lncRNA ADAMTS9-AS2 suppresses progression of lung cancer through inhibition of miR-223-3p and promotion of TGFBR3. IUBMB Life 70, 536–546. doi: 10.1002/iub.1752
Liu, D. C., Song, L. L., Liang, Q., Hao, L., Zhang, Z. G., and Han, C. H. (2019). Long noncoding RNA LEF1-AS1 silencing suppresses the initiation and development of prostate cancer by acting as a molecular sponge of miR-330-5p via LEF1 repression. J. Cell. Physiol. 234, 12727–12744. doi: 10.1002/jcp.27893
Liu, X. L., Xiao, K., Xue, B., Yang, D., Lei, Z., Shan, Y., et al. (2013). Dual role of TGFBR3 in bladder cancer. Oncol. Rep. 30, 1301–1308. doi: 10.3892/or.2013.2599
Meng, S., Zhou, H., Feng, Z., Xu, Z., Tang, Y., Li, P., et al. (2017). CircRNA: functions and properties of a novel potential biomarker for cancer. Mol. Cancer 16:94. doi: 10.1186/s12943-017-0663-2
Nishan, U., Rosa-Ribeiro, R., Cesar, C. L., and Carvalho, H. F. (2019). Transcription regulators are transiently expressed during the prostate gland adaptation to the hypoandrogenic environment. Histol. Histopathol. 34, 1025–1036. doi: 10.14670/hh-18-105
Ou, R., Lv, J., Zhang, Q., Lin, F., Zhu, L., Huang, F., et al. (2020). circAMOTL1 motivates AMOTL1 expression to facilitate cervical cancer growth. Mol. Ther. Nucleic Acids 19, 50–60. doi: 10.1016/j.omtn.2019.09.022
Shan, G., Shao, B., Liu, Q., Zeng, Y., Fu, C., Chen, A., et al. (2020). circFMN2 sponges miR-1238 to promote the expression of LIM-homeobox gene 2 in prostate cancer cells. Mol. Ther. Nucleic Acids 21, 133–146. doi: 10.1016/j.omtn.2020.05.008
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303
Sharifi, N., Hurt, E. M., Kawasaki, B. T., and Farrar, W. L. (2007). TGFBR3 loss and consequences in prostate cancer. Prostate 67, 301–311. doi: 10.1002/pros.20526
Shen, J., Chen, L., Cheng, J., Jin, X., Mu, Y., Li, Q., et al. (2019). Circular RNA sequencing reveals the molecular mechanism of the effects of acupuncture and moxibustion on endometrial receptivity in patients undergoing infertility treatment. Mol. Med. Rep. 20, 1959–1965. doi: 10.3892/mmr.2019.10386
Song, Z., Zhuo, Z., Ma, Z., Hou, C., Chen, G., and Xu, G. (2019). Hsa_Circ_0001206 is downregulated and inhibits cell proliferation, migration and invasion in prostate cancer. Artif. Cells Nanomed. Biotechnol. 47, 2449–2464. doi: 10.1080/21691401.2019.1626866
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 45, W98–W102. doi: 10.1093/nar/gkx247
Thomson, D. W., and Dinger, M. E. (2016). Endogenous microRNA sponges: evidence and controversy. Nat. Rev. Genet. 17, 272–283. doi: 10.1038/nrg.2016.20
Uhlen, M., Zhang, C., Lee, S., Sjöstedt, E., Fagerberg, L., Bidkhori, G., et al. (2017). A pathology atlas of the human cancer transcriptome. Science 357:eaan2507. doi: 10.1126/science.aan2507
Verduci, L., Strano, S., Yarden, Y., and Blandino, G. (2019). The circRNA-microRNA code: emerging implications for cancer diagnosis and treatment. Mol. Oncol. 13, 669–680. doi: 10.1002/1878-0261.12468
Vo, J. N., Cieslik, M., Zhang, Y., Shukla, S., Xiao, L., Zhang, Y., et al. (2019). The landscape of circular RNA in cancer. Cell 176, 869.e13–881.e123. doi: 10.1016/j.cell.2018.12.021
Wang, X., Wang, R., Wu, Z., and Bai, P. (2019). Circular RNA ITCH suppressed prostate cancer progression by increasing HOXB13 expression via spongy miR-17-5p. Cancer Cell Int. 19:328. doi: 10.1186/s12935-019-0994-8
Wilusz, J. E., and Sharp, P. A. (2013). Molecular biology. A circuitous route to noncoding RNA. Science 340, 440–441. doi: 10.1126/science.1238522
Wu, G., Sun, Y., Xiang, Z., Wang, K., Liu, B., Xiao, G., et al. (2019). Preclinical study using circular RNA 17 and micro RNA 181c-5p to suppress the enzalutamide-resistant prostate cancer progression. Cell Death Dis. 10:37. doi: 10.1038/s41419-018-1048-1
Xia, Q., Ding, T., Zhang, G., Li, Z., Zeng, L., Zhu, Y., et al. (2018). Circular RNA expression profiling identifies prostate cancer- specific circRNAs in prostate cancer. Cell. Physiol. Biochem. 50, 1903–1915. doi: 10.1159/000494870
Xie, J., Wang, S., Li, G., Zhao, X., Jiang, F., Liu, J., et al. (2019). circEPSTI1 regulates ovarian cancer progression via decoying miR-942. J. Cell. Mol. Med. 23, 3597–3602. doi: 10.1111/jcmm.14260
Xue, D., Lu, H., Xu, H.-Y., Zhou, C.-X., and He, X.-Z. (2018). Long noncoding RNA MALAT1 enhances the docetaxel resistance of prostate cancer cells via miR-145-5p-mediated regulation of AKAP12. J. Cell. Mol. Med. 22, 3223–3237. doi: 10.1111/jcmm.13604
Yan, Z., Xiao, Y., Chen, Y., and Luo, G. (2020). Screening and identification of epithelial-to-mesenchymal transition-related circRNA and miRNA in prostate cancer. Pathol. Res. Pract. 216:152784. doi: 10.1016/j.prp.2019.152784
Zhang, H.-D., Jiang, L.-H., Sun, D.-W., Hou, J.-C., and Ji, Z.-L. (2018). CircRNA: a novel type of biomarker for cancer. Breast Cancer 25, 1–7. doi: 10.1007/s12282-017-0793-9
Zhang, Y., Su, X., Kong, Z., Fu, F., Zhang, P., Wang, D., et al. (2017). An androgen reduced transcript of LncRNA GAS5 promoted prostate cancer proliferation. PLoS One 12:e0182305. doi: 10.1371/journal.pone.0182305
Zheng, Y., Chen, C. J., Lin, Z. Y., Li, J. X., Liu, J., Lin, F. J., et al. (2020a). Circ_KATNAL1 regulates prostate cancer cell growth and invasiveness through the miR-145-3p/WISP1 pathway. Biochem Cell Biol. 98, 396–404. doi: 10.1139/bcb-2019-0211
Zheng, Y., Li, J. X., Chen, C. J., Lin, Z. Y., Liu, J. X., and Lin, F. J. (2020b). Extracellular vesicle-derived circ_SLC19A1 promotes prostate cancer cell growth and invasion through the miR-497/septin 2 pathway. Cell Biol. Int. 44, 1037–1045. doi: 10.1002/cbin.11303
Keywords: prostate cancer, circRNA, bioinformatics analysis, circRNA–microRNA–mRNA interaction axis, signaling pathway
Citation: Wu Y-P, Lin X-D, Chen S-H, Ke Z-B, Lin F, Chen D-N, Xue X-Y, Wei Y, Zheng Q-S, Wen Y-A and Xu N (2020) Identification of Prostate Cancer-Related Circular RNA Through Bioinformatics Analysis. Front. Genet. 11:892. doi: 10.3389/fgene.2020.00892
Received: 01 May 2020; Accepted: 20 July 2020;
Published: 14 August 2020.
Edited by:
Mehdi Pirooznia, National Heart, Lung, and Blood Institute (NHLBI), United StatesReviewed by:
Mustafa Ozen, Baylor College of Medicine, United StatesNeeti Sharma, Symbiosis International University, India
Nayara Izabel Viana, University of São Paulo, Brazil
Copyright © 2020 Wu, Lin, Chen, Ke, Lin, Chen, Xue, Wei, Zheng, Wen and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ning Xu, drxun@fjmu.edu.cn; Yao-An Wen, wyaxyz@qq.com
†These authors have contributed equally to this work