Next Article in Journal
The Timing and Effects of Low-Dose Ethanol Treatment on Acetaminophen-Induced Liver Injury
Next Article in Special Issue
AC016745.3 Regulates the Transcription of AR Target Genes by Antagonizing NONO
Previous Article in Journal
Efficient Isolation of Bacterial RNAs Using Silica-Based Materials Modified with Ionic Liquids
Previous Article in Special Issue
Titration of Androgen Signaling: How Basic Studies Have Informed Clinical Trials Using High-Dose Testosterone Therapy in Castrate-Resistant Prostate Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comprehensive Characterization of Androgen-Responsive circRNAs in Prostate Cancer

1
Shanghai Engineering Research Center of Industrial Microorganisms, School of Life Science, Fudan University, Shanghai 200433, China
2
Department of Urology, Shanghai Fourth People’s Hospital Affiliated to Tongji University School of Medicine, Shanghai 200434, China
*
Authors to whom correspondence should be addressed.
Life 2021, 11(10), 1096; https://doi.org/10.3390/life11101096
Submission received: 22 August 2021 / Revised: 20 September 2021 / Accepted: 21 September 2021 / Published: 15 October 2021
(This article belongs to the Special Issue Androgen Receptor and AR-Related Signaling in Health and Disease)

Abstract

:
The androgen receptor (AR) signaling pathway plays an important role in the initiation and progression of prostate cancer. Circular RNAs (circRNAs), the novel noncoding RNAs without 5′ to 3′ polarity or 3′ poly (A), play an important role in multiple diseases. However, the potential roles of androgen-responsive circRNAs in prostate cancer remain unclear. In this study, we identified 3237 androgen-responsive circRNAs and 1954 androgen-responsive mRNAs after dihydrotestosterone (DHT) stimulation using microarray. Among them, the expression of 1296 androgen-responsive circRNAs was consistent with that of their parent genes, and we thought AR might regulate the expression of these circRNAs at the transcriptional level. In addition, 1941 circRNAs expression was not consistent with their parent genes, and we speculated that AR may regulate the expression of those circRNAs at the posttranscriptional level through affecting alternative splicing. Analyzing the androgen-responsive circRNAs regulated at the posttranscriptional level, we identified two key RNA binding proteins (RBPs), WTAP and TNRC6, using the circInteractome database, which may play important role in the biogenesis of androgen-responsive circRNAs. Furthermore, we explored the potential biological functions and predicted the molecular mechanisms of two dysregulated circRNAs (circNFIA and circZNF561) in prostate cancer. In this study, we revealed that circNFIA was upregulated in prostate cancer tissues and plasma samples from patients with prostate cancer; circNFIA may play an oncogenic role in prostate cancer. In contrast, circZNF561 was downregulated and may act as a tumor suppressor in prostate cancer. Our results suggest that androgen-responsive circRNAs might regulate the progression of prostate cancer and could be novel diagnostic biomarkers.

1. Introduction

Prostate cancer is the second most frequently diagnosed cancer in men worldwide, and 1.4 million estimated cases were diagnosed in 2020 [1]. The incidence rates of prostate cancer increased in most countries, such as Sweden, Thailand, China, and Lithuania [2]. Prostate cancer is hormone-responsive, and castration leads to tumor regression in prostate cancer patients [3]. Androgen deprivation therapy (ADT) is the main treatment for advanced and metastatic prostate cancer [4]. However, almost all patients develop castration-resistant prostate cancer (CRPC) approximately 2–3 years after ADT, at which time, the serum PSA levels increase again and the disease rapidly worsens [5,6].
The androgen receptor (AR) signaling pathway plays a vital role in the diagnosis and treatment of prostate cancer. CRPC, previously defined as hormone-independent prostate cancer, is now known to still be androgen dependent. AR amplification [7], AR mutations [8], mutations in coactivators or corepressors [9,10,11], androgen-independent AR activation [12], and alternative androgen production [13,14] contribute to the progression of CRPC. Understanding the role of the AR signaling pathway in the progression of prostate cancer is important to develop future therapies.
With advances in high-throughput sequencing technologies and bioinformatics, a large number of circRNAs have been gradually discovered in mammalian cells [15]. CircRNAs, a novel type of noncoding RNAs with a covalent closed-loop structure, are derived from pre-mRNAs via back-splicing and have no 5′ cap or 3′ poly (A) [16,17]. CircRNAs are more stable and protected against RNA exonucleases in vivo. Accordingly, recent studies have suggested that circRNAs could be promising biomarkers in multiple diseases [18]. In addition, an increasing number of studies have found that circRNAs play important roles in diseases, such as diabetes [19], cardiovascular disease [20], chronic inflammatory diseases [21], and neurological diseases [22], especially in a variety of tumors. For example, circSMARCA5 acts as a tumor suppressor in liver cancer [23], malignant glioma [24], non-small cell lung cancer [25], and plays an oncogenic role in prostate cancer [26]. CircRNAs act as microRNA sponges to regulate the development and progression of disease. For example, CDR1as, also called ciRS-7, was reported to have approximately 70 conserved miR-7 binding sites and to act as a potential miR-7 sponge to contribute to several biological processes [27,28]. Recently, our group also found that circFOXO3 could promote prostate cancer progression by competitive binding of miR-29a-3p [29]. However, the biogenesis and biological functions of androgen-responsive circRNAs in prostate cancer still need to be further explored.
In this study, we identified circRNAs as diagnostic and therapeutic markers and explored the potential functions of circRNAs in prostate cancer. We identified a set of androgen-responsive circRNAs and two important regulatory factors, WTAP and TNRC6, which may involve in circRNA biogenesis. In addition, we explored the potential biological functions and predicted the molecular mechanism of two androgen-responsive circRNAs (circNFIA and circZNF561), which may act as promising biomarkers for the diagnosis and treatment of prostate cancer.

2. Materials and Methods

2.1. Human Prostate Cancer Samples

Prostate cancer patient tissue samples and matched adjacent normal samples were obtained from patients who underwent radical prostatectomy at Fudan University Shanghai Cancer Center. Serum samples from newly diagnosed prostate cancer patients (n = 26) and healthy donors (n = 19) who provided informed consent were obtained at Center named above. These patients did not receive any preoperative treatment. The study was approved by the Research Ethics Committee of Fudan University Shanghai Cancer Center. All the samples were collected and used for gene expression analysis by qPCR.

2.2. Microarray Analysis

LNCaP cells were starved for 3 days and treated with 100 nM DHT for 8 h. Then, the cells were collected, and microarray analysis was performed as previously described [30] (Shanghai Biotechnology Corporation, Shanghai, China). In brief, total RNA was amplified, labeled with a Low Input Quick Amp Labeling Kit and then purified using RNeasy mini kit. Each slide was hybridized with 1.65 μg Cy3-labeled RNA for 17 h and then scanned by an Agilent microarray scanner (Agilent technologies, Santa Clara, CA, USA). The data were extracted using Feature Extraction software 10.7 (Agilent Technologies). The raw data were normalized by the quantile algorithm and limma package in R (Version 3.6.2).

2.3. Cell Culture and Treatment

LNCaP, DU145, 22Rv1, PC-3, and WPMY-1 cells were purchased from the American Type Culture Collection (ATCC, Manassas, VA, USA). All the cells were cultured in RPMI-1640 or DMEM medium supplemented with 10% fetal bovine serum (ExCell Bio, Shanghai, China). The cells were maintained in an incubator at 37 °C containing 5% CO2. The androgen treatment assays were performed as described previously [31]. Briefly, LNCaP cells were starved for 3 days and then treated with DHT. The cells were collected and used in the following detection.

2.4. RNA extraction and qPCR

Total RNA was isolated from homogenized prostate cancer tissue or cultured cells using MagZol Reagent (Magen, Guangzhou, China) as previously described [29]. Total RNA was reverse transcribed using the PrimeScript RT reagent Kit (Takara, Japan) according to the manufacturer’s protocol. The relative expression of candidate genes was determined using AceQ qPCR SYBR Green Master Mix (Vazyme, Shanghai, China) on LightCycler® 480II (Roche, Basel, Switzerland). The primers used for the qPCR analysis are listed in Supplementary Table S1.

2.5. Transfection

Cells were seeded in 6-well plates until cell confluence reached 60%–80%. Liposomal cocktails with siRNA (GenePharma, Shanghai, China) (50 nM final) were generated with HilyMax (Dojindo Laboratories, Kumamoto, Japan) in Opti-MEM (Invitrogen, Carlsbad, CA, USA) according to the protocol. The transfected cells were incubated for 48 h and then used for gene expression analysis or other experiments. The sequences of the specific siRNAs are shown in Supplementary Table S1.

2.6. Cell Apoptosis

Cells were cultured in 6-well plates after transfection. Cell apoptosis was measured using the Annexin V-FITC apoptosis kit (Dojindo Laboratories, Kumamoto, Japan). After incubation at room temperature for 15 min, cell apoptosis was analyzed on a FACSCalibur flow cytometer (Becton, Dickinson and Company, Franklin Lakes, NJ, USA).

2.7. Cell Proliferation

Cell Counting Kit-8 (Dojindo Laboratories, Kumamoto, Japan) was used to assess cell proliferation according to the manufacturer’s protocol. Briefly, transfected cells were seeded in 96-well plates and cultured in 100 μL medium. 10 μL of CCK-8 solution was added to each well at 0, 24, 48, and 72 h and further incubated for 2 h to allow the colorimetric reaction to occur. The absorbance at 450 nm was measured by a Microplate Reader ELx808 (Bio-Tek, Winooski, VT, USA) to determine cell viability.

2.8. Cell Cycle

Transfected cells were cultured in 6-well plates for 48 h. Then, the cells were collected and treated with Triton X-100 (0.03%, Sigma, Darmstadt, Germany) and propidium iodide (PI, 50 ng/mL, Beyotime, Shanghai, China) for 15 min. Cell cycle analysis was performed by a FACSCalibur flow cytometer (Becton, Dickinson and Company, Franklin Lakes, NJ, USA).

2.9. Construction of ceRNA Network

The circBank (http://www.circbank.cn/ (accessed on 18 September 2021)) was used to predict miRNA binding sites. We found 59 miRNAs targeted to circNFIA and 26 miRNAs targeted to circZNF561. The networks were drawn using Cytoscape 3.0. Following this, we identified the target genes of miRNAs interacting with circNFIA or circZNF561 using the miRDB database. The target score ≥ 90 was chosen to construct circRNA-miRNA-mRNA networks using Cytoscape 3.0.

2.10. Statistical Analysis

The numerical data are presented as the mean ± standard deviation (SD) of at least three determinations. Statistical comparisons between groups of normalized data were performed using t-tests or Mann–Whitney U-tests according to the test conditions. A p < 0.05 was considered statistically significant with a 95% confidence level.

3. Results

3.1. Identification of Androgen-Responsive circRNAs in Prostate Cancer

To identify the differentially expressed mRNAs and circRNAs after DHT stimulation, we performed high-throughput microarray analysis and the expression of 88,751 circRNAs and 18,854 mRNAs was detected. Of these, 3237 circRNAs and 1954 mRNAs were differentially expressed. As shown in Figure 1A,B, 2295 circRNAs and 969 mRNAs were upregulated and 942 circRNAs and 985 mRNAs were downregulated after DHT stimulation. Interestingly, the number of upregulated circRNAs was approximately 2.5 times that of downregulated circRNAs, while the number of upregulated mRNAs was close to that of downregulated mRNAs. The result suggests that the AR signaling pathway may have differences regulating the expression of mRNAs and circRNAs and circRNAs may not be meaningless splicing byproducts.
We also analyzed the distribution of the circRNAs on the human chromosomes and we found that the number of upregulated circRNAs derived from each chromosome was greater than the number of downregulated circRNAs (Figure 1C). In addition, most circRNAs were derived from exons (CDS region, 3′UTR, and 5′UTR) and a few circRNAs were derived from introns, intergenic regions, and antisense strands (Figure 1D).

3.2. Verification of Androgen-Responsive circRNAs

To verify the expression of androgen-responsive circRNAs, we quantified part of the androgen-responsive circRNAs expression in LNCaP cells after treatment with different doses of DHT (0, 0.1, 1, 10, 100, and 1000 nM) by qPCR. Among them, 12 of 20 circRNAs (hsa-circ-0076151, hsa-circ-0076150, hsa-circ-0001055, hsa-circ-0001599, hsa-circ-0005035, hsa-circ-0005079, hsa-circ-0052505, hsa-circ-0076157, hsa-circ-0085656, hsa-circ-0072906, hsa-circ-0029943, and hsa-circ-0006404) were significantly upregulated and 8 of 20 circRNAs (hsa-circ-0086235, hsa-circ-0008326, hsa-circ-0056198, hsa-circ-0018744, hsa-circ-0043541, hsa-circ-0071616, hsa-circ-0073483, and hsa-circ-0012755) were significantly downregulated after DHT stimulation (Figure 2). The expression of those circRNAs were consistent with the microarray data.

3.3. Characteristics of Androgen-Responsive circRNA Expression

Both mRNAs and circRNAs are derived from pre-mRNA by alternative splicing. We analyzed the expression of circRNAs and their parent genes in the microarray data. The results showed that only 51.99% of the 2087 circRNAs derived from androgen-upregulated mRNAs were upregulated and 18.22% of 1158 circRNAs derived from androgen-downregulated mRNAs were downregulated (Figure 3A,B). The results showed that only part of the circRNAs derived from androgen-induced or -reduced genes were upregulated or downregulated after DHT treatment, which suggested that the circRNAs and their parent genes had different regulated patterns.
According to the expression of mRNAs and circRNAs after DHT treatment, we divided androgen-responsive circRNAs into six groups (Figure 3D). We focused on the circRNAs in groups a and d that were regulated by androgen at the transcriptional level and the circRNAs in groups b and e that were regulated by androgen at the posttranscriptional level. Among them, the expression of hsa-circ-0076151, hsa-circ-0001055, hsa-circ-0005035, hsa-circ-0005079, and hsa-circ-0073627 was increased after DHT treatment (group a) (Figure 2), while hsa-circ-0008326 expression was decreased after DHT treatment (group d) (Figure 2).

3.4. Androgen Regulated circRNAs at tRanscriptional Level

To further verify the expression levels of circRNAs and their parent genes, we detected the expression of the above six circRNAs (in groups a and d) and their parent genes at a series of time points after DHT treatment (Figure 4A). The results showed that the expression level of these six circRNAs was basically the same as that of their parent genes. The expression of hsa-circ-0076151, hsa-circ-0001055, hsa-circ-0005035, hsa-circ-0005079, hsa-circ-0073627, and their parent genes FKBP5, AFF3, IGF1R, MBOAT2, KCNN2 was increased after DHT treatment, while hsa-circ-0008326 and RFX3 expression was decreased.
To test the hypothesis that the expression of these six circRNAs were transcriptionally regulated by androgen, we predicted the AREs within the upstream 10 kb from the transcription start site using the hTFtarget database. As shown in Figure 4B, we found many AR binding sites in those genes’ promoters. Furthermore, a ChIP-qPCR assay was performed using anti-AR antibodies to pull down DNA fragments with AREs in LNCaP cells with or without DHT stimulation. Consistently, AR was directly associated with the FKBP5, AFF3, RFX3, IGF1R, MBOAT2, and KCNN2 promoters after DHT treatment (Figure 4C). The results showed that AR might directly regulate the expression of FKBP5, AFF3, RFX3, IGF1R, MBOAT2, KCNN2, and their corresponding circRNAs.

3.5. Androgen Posttranscriptional Regulated circRNA Expression

We identified a large number of differentially expressed circRNAs, but their host genes were not differentially expressed (group b and group e). We speculated that the androgen may affect the alternative splicing and regulate these circRNAs expression in the posttranscriptional level. To further explore the regulation of circRNA expression, we firstly identified the RBPs that interacted with the circRNAs in group b and group e using CircInteractome. Then, we calculated the coefficient μ, where
( μ = number   of   androgen - induced   circRNA   interacted   with   a   specific   RBP ( a RBP ) Total   number   of   androgen - induced   circRNA ( a ) number   of   androgen - reduced   circRNA   interacted   with   a   specific   RBP ( b RBP ) Total   number   of   androgen - induced   circRNA ( b ) )  
of each RBP. The higher the μ value was, the more likely it was that the RBPs interacted with androgen-induced circRNAs. Conversely, the smaller the μ value was, the more likely it was that the RBPs interacted with androgen-reduced circRNAs (Figure 5A). The results showed that WTAP, TAF15, FOX2, and RBPMS potentially regulated the androgen-downregulated circRNAs, while TNRC6, C17ORF85, AUF1, FXR1, MOV10, FXR2, and ZC3H7B potentially regulated the androgen-upregulated circRNAs. Among those RBPs, WTAP, and TNRC6 were the most significant (Figure 5A). Next, we detected the expression of circRNAs potentially associated with WTAP or TNRC6 after knockdown of WTAP or TNRC6, respectively. As shown in Figure 5B, after knockdown of WTAP, the expression of hsa_circ_0012079, hsa_circ_0087714, and hsa_circ_0018744 was significantly decreased. Similarly, after knockdown of TNRC6, the expression of hsa_circ_0090744, hsa_circ_0076815, hsa_circ_0049519, hsa_circ_0014100, and hsa_circ_0061590 was significantly decreased (Figure 5C). We initially revealed that WTAP and TNRC6 may be involved in those circRNA’s biogenesis.

3.6. The Expression of circNFIA and circZNF561 in Prostate Cancer

To explore the biological roles of androgen-responsive circRNAs in prostate cancer, we identified differentially expressed circRNAs in 30 pairs of prostate cancer tissue samples and corresponding adjacent normal tissues by qPCR. Two circRNAs, hsa-circ-0012755 (also called circNFIA) and hsa-circ-0049154 (also called circZNF561), were significantly dysregulated in prostate cancer. The results showed that circNFIA expression was increased and circZNF561 expression was decreased in prostate cancer tissue samples compared to adjacent normal tissues (Figure 6A,B). Furthermore, we also detected the expression of circNFIA and circZNF561 in prostate cancer cells and plasma samples. As shown in Figure 6C,D, the expression of circNFIA was significantly upregulated in prostate cancer cell lines (LNCaP, LNCaP-AI, DU145, 22Rv1, and PC-3 cells), while circZNF561 was downregulated in prostate cancer cells. The expression of circNFIA and circZNF561 in plasma samples was further detected by qPCR. Consistent with the above results, circNFIA was significantly overexpressed in the plasma samples of prostate cancer patients, while the expression of circZNF561 was decreased (Figure 6E,F). These results suggest that circNFIA and circZNF561 may play roles in the progression of prostate cancer.

3.7. The Biological Functions of circNFIA and circZNF561 in Prostate Cancer

Next, we knocked down circNFIA and circZNF561 expression in the DU145 and LNCaP cell lines and explored their possible biological role in prostate cancer. As shown in Figure 7A,B, sicircNFIA and sicircZNF561 significantly decreased the expression of circNFIA and circZNF561. The CCK-8 assay showed that knockdown of circNFIA significantly inhibited cell proliferation and sicircZNF561 promoted proliferation in prostate cancer (Figure 7C–F).
In addition, we also tested the effect of knocking down circNFIA and circZNF561 on the cell cycle. The result showed that knocking down circNFIA increased the percentage of cells in the G1 phase (49.5 to 66.62% in DU145, 60.94 to 67.86% in LNCaP) and decreased the percentage of cells in the S phase (39.5 to 27.72% in DU145, 25.49 to 21.53% in LNCaP) to arrest cell cycle progression (Figure 8A,B). In contrast, sicircZNF561 decreased the percentage of cells in the G1 phase and increased the percentage of cells in the S phase, these results suggested sicircZNF561 could promote the cell cycle (Figure 8C,D).
Similarly, apoptosis assays showed that sicircNFIA significantly promoted cell apoptosis in DU145 and LNCaP cells. Compared with the control, apoptotic cells increased from 6.79% (DU145) and 7.07% (LNCaP) to 8.46% (DU145) and 7.85% (LNCaP) (Figure 9A,B). After knocking down circZNF561, we detected the percentage of apoptosis cells was decreased (Figure 9C,D). Taken together, we found that circNFIA may play an oncogenic role and circZNF561 may act as a tumor suppressor in prostate cancer. These results indicated that circNFIA and circZNF561 can be used as candidate tumor biomarkers.

3.8. Construction of a circRNA-miRNA ceRNA Network

As mentioned in previous reports, endogenous circRNAs could act as efficient microRNA sponges to regulate the expression of protein-coding genes [29,32]. To explore the mechanism by which circNFIA and circZNF561 regulated prostate cancer progression, the putative candidate miRNAs binding to circNFIA or circZNF561 were predicted using circBank. As shown in Figure 10A,E, we found 59 miRNAs targeted to circNFIA and 26 miRNAs targeted to circZNF561. Next, we constructed a circRNA-mediated competing endogenous RNA (ceRNA) network. Using the miRDB database, we identified the target genes of miRNAs interacting with circNFIA or circZNF561. The networks were drawn using Cytoscape 3.0 (Figure 10B,F)
Moreover, we performed Gene Ontology (GO) and KEGG pathway analysis of potential target genes in the network using the DAVID Bioinformatics Database (Figure 10C,D,G,H). Our results showed that target genes within the circNFIA-related network were mainly enriched in the cGMP-PKG, Rap1, Ras, PI3K-Akt and MAPK signaling pathways and were associated with regulating cell migration, cell-cell adhesion and apoptotic processes. (Figure 10C,D). The circZNF561-related network was mainly enriched in the Rap1, Ras, MAPK, and PI3K-Akt signaling pathways and was associated with regulating protein ubiquitination, GTPase activity, apoptotic processes, and cell proliferation (Figure 10G,H).

4. Discussion

ADT is the main treatment for prostate cancer, especially advanced and metastatic prostate cancer. With AR mutation or amplification, a low concentration of androgen can still activate the AR signaling pathway in CRPC, causing rapid worsening of prostate cancer in patients. Therefore, our group and others focused on exploring the important regulator of the AR signaling pathway in prostate cancer. Androgen-responsive noncoding RNAs, including miRNAs [33] and lncRNAs [31], have been reported, while androgen-responsive circRNAs have rarely been studied. In this study, we identified androgen-responsive circRNAs through high-throughput microarray analysis. The expression of six circRNAs (hsa-circ-0076151, hsa-circ-0001055, hsa-circ-0005035, hsa-circ-0005079, hsa-circ-0008326, and hsa-circ-0073627) regulated at the transcriptional level was validated by qPCR. Importantly, ChIP-seq data showed that the AR peaks in these circRNA loci significantly increased after DHT treatment, suggesting that these circRNAs were direct targets of AR.
Recent studies have reported that the biogenesis of circRNAs is mainly regulated by cis-elements or trans-factors, that is, through complementary base pairing on the flanking sequence [15,16] or binding to RBPs special motifs [34,35]. In this study, we identified RBPs involved in the biogenesis of circRNA using CircInteractome. We found that WTAP mainly participated in the formation of androgen-reduced circRNAs, while TNRC6 mainly participated in androgen-induced circRNAs. To further explore the importance of WTAP/TNRC6 in the biogenesis of circRNAs, we knocked down WTAP or TNRC6 and detected candidate circRNAs expression. The expression of hsa_circ_0012079, hsa_circ_0087714, and hsa_circ_0018744 was significantly decreased after WTAP knockdown. Similarly, the expression of hsa_circ_0090744, hsa_circ_0076815, hsa_circ_0049519, hsa_circ_0014100, and hsa_circ_0061590 was significantly decreased after TNRC6 knockdown. These results indicated that WTAP and TNRC6 might both involve in circRNA processing.
CircRNAs were more stable and tissue-specific than linear RNAs. Recent studies have reported that circRNAs could act as promising noninvasive markers for cancer diagnosis and prognosis [36,37]. In addition, abnormal expression of circRNAs has been observed in multiple types of cancer. In this study, we found that two circRNAs, circNFIA and circZNF561, were significantly differentially expressed in prostate cancer cells, tissue samples and plasma samples from patients with prostate cancer. Our results showed that circNFIA was significantly upregulated in prostate cancer cells, prostate cancer tissues, and plasma samples, while circZNF561 was downregulated in prostate cancer cells, tissues and plasma samples. Therefore, circNFIA and circZNF561 are potential diagnostic markers in prostate cancer.
Many studies have reported that circRNA expression is related to the etiology and progression of cancer and is involved in regulating the key signaling pathway of cancers [38]. To characterize the biological roles of circNFIA and circZNF561, we knocked down circNFIA and circZNF561 and explored the effect on cell proliferation, cell cycle, and cell apoptosis. In this study, we revealed that circNFIA may play an oncogenic role, while circZNF561 may act as a tumor suppressor in prostate cancer.
As mentioned in previous reports, circRNAs can act as miRNA sponges to regulate target genes and affect tumor progression. Wei et al. reported that CDR1as inhibited the tumor suppressor miR-7 to promote cell proliferation in NSCLC [39]. Zhang et al. reported that circFOXO3 can act as a miR-138-5p and miR-432-5p sponge to promote the expression of NFAT5 and promote the progression of glioblastoma [40]. CircFOXO3 could also promote the occurrence and development of prostate cancer by competitively binding miR-29a-3p [29]. To explore the potential mechanism by which circNFIA and circZNF561 played a role in prostate cancer, we constructed circNFIA- or circZNF561-mediated ceRNA networks by using the circBank and miRDB databases. We demonstrate that this study provides useful information for exploring potential therapeutic and prognostic biomarkers for prostate cancer.

5. Conclusions

In conclusion, we performed a androgen-responsive circRNAs analysis in prostate cancer patient tissues and identified two important regulatory factors, WTAP and TNRC6, which may involve in circRNA biogenesis. In addition, we explored the potential biological functions and predicted the molecular mechanism of two androgen-responsive circRNAs (circNFIA and circZNF561) by computational analysis and experimental verification, which may act as promising biomarkers for the diagnosis and treatment of prostate cancer.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/life11101096/s1, Table S1: Sequence information of primers.

Author Contributions

Conceptualization, Z.K. and Y.L. (Yao Li); methodology, Z.K., Y.L. (Yali Lu), Y.L. (Yao Li), and Y.X.; software, X.W.; validation, Z.K., J.L., D.L., Y.L. (Yao Li), and Y.X.; formal analysis, D.L., Y.H., and C.W.; data curation, Z.K., Y.L. (Yao Li), and Y.X.; writing—original draft preparation, Z.K.; writing—review and editing, Y.L. (Yao Li), and Y.X.; visualization, Z.K. and Y.L. (Yali Lu); supervision, Y.L. (Yao Li) and Y.X.; project administration, Y.L. (Yao Li), and Y.X.; funding acquisition, Y.L. (Yao Li) and Y.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Shanghai Science and Technology Development Foundation (20ZR1404500), the China’s 13th Five Year Programs for the prevention and cure of great infectious diseases (2018ZX10731301-001-005), and Project SKLGE-1901, supported by the Open Research Funds of the State Key Laboratory of Genetic Engineering, Fudan University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef]
  2. Center, M.M.; Jemal, A.; Lortet-Tieulent, J.; Ward, E.; Ferlay, J.; Brawley, O.; Bray, F. International variation in prostate cancer incidence and mortality rates. Eur Urol. 2012, 61, 1079–1092. [Google Scholar] [CrossRef]
  3. Huggins, C.; Hodges, C.V. Studies on prostatic cancer. I. The effect of castration, of estrogen and androgen injection on serum phosphatases in metastatic carcinoma of the prostate. CA Cancer J. Clin. 1972, 22, 232–240. [Google Scholar] [CrossRef]
  4. Scher, H.I.; Liebertz, C.; Kelly, W.K.; Mazumdar, M.; Brett, C.; Schwartz, L.; Kolvenbag, G.; Shapiro, L.; Schwartz, M. Bicalutamide for advanced prostate cancer: The natural versus treated history of disease. J. Clin. Oncol. 1997, 15, 2928–2938. [Google Scholar] [CrossRef]
  5. Chandrasekar, T.; Yang, J.C.; Gao, A.C.; Evans, C.P. Mechanisms of resistance in castration-resistant prostate cancer (CRPC). Transl. Androl. Urol. 2015, 4, 365–380. [Google Scholar]
  6. Kelly, W.K.; Scher, H.I. Prostate specific antigen decline after antiandrogen withdrawal: The flutamide withdrawal syndrome. J. Urol. 1993, 149, 607–609. [Google Scholar] [CrossRef]
  7. Visakorpi, T.; Hyytinen, E.; Koivisto, P.; Tanner, M.; Keinanen, R.; Palmberg, C.; Palotie, A.; Tammela, T.; Isola, J.; Kallioniemi, O.P. In vivo amplification of the androgen receptor gene and progression of human prostate cancer. Nat. Genet. 1995, 9, 401–406. [Google Scholar] [CrossRef]
  8. Gregory, C.W.; Johnson, R.J.; Mohler, J.L.; French, F.S.; Wilson, E.M. Androgen receptor stabilization in recurrent prostate cancer is associated with hypersensitivity to low androgen. Cancer Res. 2001, 61, 2892–2898. [Google Scholar]
  9. Makridakis, N.M.; di Salle, E.; Reichardt, J.K. Biochemical and pharmacogenetic dissection of human steroid 5 alpha-reductase type II. Pharmacogenetics 2000, 10, 407–413. [Google Scholar] [CrossRef]
  10. Wang, Q.; Carroll, J.S.; Brown, M. Spatial and temporal recruitment of androgen receptor and its coactivators involves chromosomal looping and polymerase tracking. Mol. Cell 2005, 19, 631–642. [Google Scholar] [CrossRef]
  11. Chung, A.C.; Zhou, S.; Liao, L.; Tien, J.C.; Greenberg, N.M.; Xu, J. Genetic ablation of the amplified-in-breast cancer 1 inhibits spontaneous prostate cancer progression in mice. Cancer Res. 2007, 67, 5965–5975. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wang, Q.; Li, W.; Zhang, Y.; Yuan, X.; Xu, K.; Yu, J.; Chen, Z.; Beroukhim, R.; Wang, H.; Lupien, M.; et al. Androgen receptor regulates a distinct transcription program in androgen-independent prostate cancer. Cell 2009, 138, 245–256. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Geller, J.; Albert, J.; Loza, D. Steroid levels in cancer of the prostate—Markers of tumour differentiation and adequacy of anti-androgen therapy. J. Steroid. Biochem. 1979, 11, 631–636. [Google Scholar] [CrossRef]
  14. Mostaghel, E.A. Abiraterone in the treatment of metastatic castration-resistant prostate cancer. Cancer Manag. Res. 2014, 6, 39–51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Ivanov, A.; Memczak, S.; Wyler, E.; Torti, F.; Porath, H.T.; Orejuela, M.R.; Piechotta, M.; Levanon, E.Y.; Landthaler, M.; Dieterich, C.; et al. Analysis of intron sequences reveals hallmarks of circular RNA biogenesis in animals. Cell Rep. 2015, 10, 170–177. [Google Scholar] [CrossRef] [Green Version]
  16. Kelly, S.; Greenman, C.; Cook, P.R.; Papantonis, A. Exon Skipping Is Correlated with Exon Circularization. J. Mol. Biol. 2015, 427, 2414–2417. [Google Scholar] [CrossRef] [PubMed]
  17. Zhang, X.O.; Wang, H.B.; Zhang, Y.; Lu, X.; Chen, L.L.; Yang, L. Complementary sequence-mediated exon circularization. Cell 2014, 159, 134–147. [Google Scholar] [CrossRef] [Green Version]
  18. Chen, S.; Li, T.; Zhao, Q.; Xiao, B.; Guo, J. Using circular RNA hsa_circ_0000190 as a new biomarker in the diagnosis of gastric cancer. Clin. Chim. Acta. 2017, 466, 167–171. [Google Scholar] [CrossRef]
  19. Xu, H.; Guo, S.; Li, W.; Yu, P. The circular RNA Cdr1as, via miR-7 and its targets, regulates insulin transcription and secretion in islet cells. Sci. Rep. 2015, 5, 12453. [Google Scholar] [CrossRef]
  20. Wang, K.; Long, B.; Liu, F.; Wang, J.X.; Liu, C.Y.; Zhao, B.; Zhou, L.Y.; Sun, T.; Wang, M.; Yu, T.; et al. A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223. Eur. Heart J. 2016, 37, 2602–2611. [Google Scholar] [CrossRef] [PubMed]
  21. Miao, Q.; Zhong, Z.; Jiang, Z.; Lin, Y.; Ni, B.; Yang, W.; Tang, J. RNA-seq of circular RNAs identified circPTPN22 as a potential new activity indicator in systemic lupus erythematosus. Lupus 2019, 28, 520–528. [Google Scholar] [CrossRef] [PubMed]
  22. Lukiw, W.J. Circular RNA (circRNA) in Alzheimer’s disease (AD). Front. Genet. 2013, 4, 307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Yu, J.; Xu, Q.G.; Wang, Z.G.; Yang, Y.; Zhang, L.; Ma, J.Z.; Sun, S.H.; Yang, F.; Zhou, W.P. Circular RNA cSMARCA5 inhibits growth and metastasis in hepatocellular carcinoma. J. Hepatol. 2018, 68, 1214–1227. [Google Scholar] [CrossRef] [PubMed]
  24. Barbagallo, D.; Caponnetto, A.; Cirnigliaro, M.; Brex, D.; Barbagallo, C.; D’Angeli, F.; Morrone, A.; Caltabiano, R.; Barbagallo, G.M.; Ragusa, M.; et al. CircSMARCA5 Inhibits Migration of Glioblastoma Multiforme Cells by Regulating a Molecular Axis Involving Splicing Factors SRSF1/SRSF3/PTB. Int. J. Mol. Sci. 2018, 19, 480. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Wang, Y.; Li, H.; Lu, H.; Qin, Y. Circular RNA SMARCA5 inhibits the proliferation, migration, and invasion of non-small cell lung cancer by miR-19b-3p/HOXA9 axis. Onco. Targets Ther. 2019, 12, 7055–7065. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Kong, Z.; Wan, X.; Zhang, Y.; Zhang, P.; Zhang, Y.; Zhang, X.; Qi, X.; Wu, H.; Huang, J.; Li, Y. Androgen-responsive circular RNA circSMARCA5 is up-regulated and promotes cell proliferation in prostate cancer. Biochem. Biophys. Res. Commun. 2017, 493, 1217–1223. [Google Scholar] [CrossRef]
  27. Lou, J.; Hao, Y.; Lin, K.; Lyu, Y.; Chen, M.; Wang, H.; Zou, D.; Jiang, X.; Wang, R.; Jin, D.; et al. Circular RNA CDR1as disrupts the p53/MDM2 complex to inhibit Gliomagenesis. Mol. Cancer. 2020, 19, 138. [Google Scholar] [CrossRef]
  28. Hanniford, D.; Ulloa-Morales, A.; Karz, A.; Berzoti-Coelho, M.G.; Moubarak, R.S.; Sanchez-Sendra, B.; Kloetgen, A.; Davalos, V.; Imig, J.; Wu, P.; et al. Epigenetic Silencing of CDR1as Drives IGF2BP3-Mediated Melanoma Invasion and Metastasis. Cancer Cell 2020, 37, 55–70. [Google Scholar] [CrossRef]
  29. Kong, Z.; Wan, X.; Lu, Y.; Zhang, Y.; Huang, Y.; Xu, Y.; Liu, Y.; Zhao, P.; Xiang, X.; Li, L.; et al. Circular RNA circFOXO3 promotes prostate cancer progression through sponging miR-29a-3p. J. Cell Mol. Med. 2020, 24, 799–813. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, C.; Han, X.; Yang, L.; Fu, J.; Sun, C.; Huang, S.; Xiao, W.; Gao, Y.; Liang, Q.; Wang, X.; et al. Circular RNA circPPM1F modulates M1 macrophage activation and pancreatic islet inflammation in type 1 diabetes mellitus. Theranostics 2020, 10, 10908–10924. [Google Scholar] [CrossRef]
  31. Wan, X.; Huang, W.; Yang, S.; Zhang, Y.; Pu, H.; Fu, F.; Huang, Y.; Wu, H.; Li, T.; Li, Y. Identification of androgen-responsive lncRNAs as diagnostic and prognostic markers for prostate cancer. Oncotarget 2016, 7, 60503–60518. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Hansen, T.B.; Jensen, T.I.; Clausen, B.H.; Bramsen, J.B.; Finsen, B.; Damgaard, C.K.; Kjems, J. Natural RNA circles function as efficient microRNA sponges. Nature 2013, 495, 384–388. [Google Scholar] [CrossRef] [PubMed]
  33. Mo, W.; Zhang, J.; Li, X.; Meng, D.; Gao, Y.; Yang, S.; Wan, X.; Zhou, C.; Guo, F.; Huang, Y.; et al. Identification of novel AR-targeted microRNAs mediating androgen signalling through critical pathways to regulate cell viability in prostate cancer. PLoS ONE 2013, 8, e56592. [Google Scholar] [CrossRef] [Green Version]
  34. Conn, S.J.; Pillman, K.A.; Toubia, J.; Conn, V.M.; Salmanidis, M.; Phillips, C.A.; Roslan, S.; Schreiber, A.W.; Gregory, P.A.; Goodall, G.J. The RNA binding protein quaking regulates formation of circRNAs. Cell 2015, 160, 1125–1134. [Google Scholar] [CrossRef] [Green Version]
  35. Errichelli, L.; Dini, M.S.; Laneve, P.; Colantoni, A.; Legnini, I.; Capauto, D.; Rosa, A.; De Santis, R.; Scarfo, R.; Peruzzi, G.; et al. FUS affects circular RNA expression in murine embryonic stem cell-derived motor neurons. Nat. Commun. 2017, 8, 14741. [Google Scholar] [CrossRef]
  36. Li, Y.; Zheng, Q.; Bao, C.; Li, S.; Guo, W.; Zhao, J.; Chen, D.; Gu, J.; He, X.; Huang, S. Circular RNA is enriched and stable in exosomes: A promising biomarker for cancer diagnosis. Cell Res. 2015, 25, 981–984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Vo, J.N.; Cieslik, M.; Zhang, Y.; Shukla, S.; Xiao, L.; Zhang, Y.; Wu, Y.M.; Dhanasekaran, S.M.; Engelke, C.G.; Cao, X.; et al. The Landscape of Circular RNA in Cancer. Cell 2019, 176, 869–881. [Google Scholar] [CrossRef] [Green Version]
  38. Kristensen, L.S.; Hansen, T.B.; Veno, M.T.; Kjems, J. Circular RNAs in cancer: Opportunities and challenges in the field. Oncogene 2018, 37, 555–565. [Google Scholar] [CrossRef] [Green Version]
  39. Zhang, X.; Yang, D.; Wei, Y. Overexpressed CDR1as functions as an oncogene to promote the tumor progression via miR-7 in non-small-cell lung cancer. Onco. Targets Ther. 2018, 11, 3979–3987. [Google Scholar] [CrossRef] [Green Version]
  40. Zhang, S.; Liao, K.; Miao, Z.; Wang, Q.; Miao, Y.; Guo, Z.; Qiu, Y.; Chen, B.; Ren, L.; Wei, Z.; et al. CircFOXO3 promotes glioblastoma progression by acting as a competing endogenous RNA for NFAT5. Neuro. Oncol. 2019, 21, 1284–1296. [Google Scholar] [CrossRef]
Figure 1. Identification of androgen-responsive circRNAs in prostate cancer. (A) The heatmap for androgen-responsive mRNAs and circRNAs in LNCaP cells after 8 h of 100 nM DHT treatment. (B) The volcano plot for androgen-responsive mRNAs and circRNAs. (C) The statistical analysis of androgen-responsive circRNAs arranged by chromosome. (D) The classification of androgen-responsive circRNAs according to their genomic localization.
Figure 1. Identification of androgen-responsive circRNAs in prostate cancer. (A) The heatmap for androgen-responsive mRNAs and circRNAs in LNCaP cells after 8 h of 100 nM DHT treatment. (B) The volcano plot for androgen-responsive mRNAs and circRNAs. (C) The statistical analysis of androgen-responsive circRNAs arranged by chromosome. (D) The classification of androgen-responsive circRNAs according to their genomic localization.
Life 11 01096 g001
Figure 2. Verification of androgen-responsive circRNAs. The expression levels of 20 circRNAs in LNCaP after different doses of DHT treatment for 48 h were detected by qPCR. Data are mean ± SD. The p-value was calculated by the unpaired two-tailed Student’s t-test. * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 2. Verification of androgen-responsive circRNAs. The expression levels of 20 circRNAs in LNCaP after different doses of DHT treatment for 48 h were detected by qPCR. Data are mean ± SD. The p-value was calculated by the unpaired two-tailed Student’s t-test. * p < 0.05, ** p < 0.01, *** p < 0.001.
Life 11 01096 g002
Figure 3. Classification of androgen-responsive circRNA expression. (A,B) The Venn diagrams for androgen-induced genes encoding androgen-induced circRNAs and androgen-reduced genes encoding androgen-reduced circRNAs. (C) The schematic diagram for the biogenesis of mRNA and circRNA. (D) The classification of androgen-responsive circRNAs according to the expression of circRNAs and their parent genes.
Figure 3. Classification of androgen-responsive circRNA expression. (A,B) The Venn diagrams for androgen-induced genes encoding androgen-induced circRNAs and androgen-reduced genes encoding androgen-reduced circRNAs. (C) The schematic diagram for the biogenesis of mRNA and circRNA. (D) The classification of androgen-responsive circRNAs according to the expression of circRNAs and their parent genes.
Life 11 01096 g003
Figure 4. Androgen-regulated circRNAs at the transcriptional level. (A) qPCR analysis of the expression of hsa-circ-0076151, hsa-circ-0001055, hsa-circ-0005035, hsa-circ-0005079, hsa-circ-0008326, hsa-circ-0073627, and their parent genes at a series of time points after 10 nM DHT treatment. Data are mean ± SD, The p-value was calculated by the unpaired two-tailed Student’s t-test. * p < 0.05, ** p < 0.01, *** p < 0.001. # p < 0.05, ## p < 0.01, ### p < 0.001 (mRNAs). (B) Prediction of AREs within 10 kb upstream of the transcription start site of FKBP5, AFF3, RFX3, IGF1R, MBOAT2, and KCNN2 using hTFtarget. (C) ChIP-qPCR showed that AR peaks in these circRNA loci significantly increased after DHT treatment.
Figure 4. Androgen-regulated circRNAs at the transcriptional level. (A) qPCR analysis of the expression of hsa-circ-0076151, hsa-circ-0001055, hsa-circ-0005035, hsa-circ-0005079, hsa-circ-0008326, hsa-circ-0073627, and their parent genes at a series of time points after 10 nM DHT treatment. Data are mean ± SD, The p-value was calculated by the unpaired two-tailed Student’s t-test. * p < 0.05, ** p < 0.01, *** p < 0.001. # p < 0.05, ## p < 0.01, ### p < 0.001 (mRNAs). (B) Prediction of AREs within 10 kb upstream of the transcription start site of FKBP5, AFF3, RFX3, IGF1R, MBOAT2, and KCNN2 using hTFtarget. (C) ChIP-qPCR showed that AR peaks in these circRNA loci significantly increased after DHT treatment.
Life 11 01096 g004
Figure 5. Identification of RBPs involved in regulating circRNA biogenesis. (A) Calculation of the ratio of each RBP interacting with androgen-induced circRNAs and androgen-reduced circRNAs. (B) Analysis of the expression of circRNAs potentially interacting with WTAP after silencing WTAP. (C) Analysis of the expression of circRNAs potentially interacting with TNRC6 after siTNRC6. The p-value was calculated by the unpaired two-tailed Student’s t-test. * p < 0.05; ** p < 0.01; *** p < 0.001 (circRNAs).
Figure 5. Identification of RBPs involved in regulating circRNA biogenesis. (A) Calculation of the ratio of each RBP interacting with androgen-induced circRNAs and androgen-reduced circRNAs. (B) Analysis of the expression of circRNAs potentially interacting with WTAP after silencing WTAP. (C) Analysis of the expression of circRNAs potentially interacting with TNRC6 after siTNRC6. The p-value was calculated by the unpaired two-tailed Student’s t-test. * p < 0.05; ** p < 0.01; *** p < 0.001 (circRNAs).
Life 11 01096 g005
Figure 6. The expression of circNFIA and circZNF561 in prostate cancer. (A,B) qPCR analysis of the expression of circNFIA and circZNF561 in 30 pairs of prostate cancer tissue samples and corresponding adjacent tissue samples. (C,D) qPCR analysis of the expression of circNFIA and circZNF561 in prostate cancer cell lines. (E,F) qPCR analysis of the expression of circNFIA and circZNF561 in 26 prostate cancer plasma samples and 19 normal plasma samples. The p-value was calculated by the unpaired two-tailed Student’s t-test. * p < 0.05; ** p < 0.01; *** p < 0.001 (circRNAs).
Figure 6. The expression of circNFIA and circZNF561 in prostate cancer. (A,B) qPCR analysis of the expression of circNFIA and circZNF561 in 30 pairs of prostate cancer tissue samples and corresponding adjacent tissue samples. (C,D) qPCR analysis of the expression of circNFIA and circZNF561 in prostate cancer cell lines. (E,F) qPCR analysis of the expression of circNFIA and circZNF561 in 26 prostate cancer plasma samples and 19 normal plasma samples. The p-value was calculated by the unpaired two-tailed Student’s t-test. * p < 0.05; ** p < 0.01; *** p < 0.001 (circRNAs).
Life 11 01096 g006
Figure 7. The biological function of circNFIA and circZNF561. (A,B) siRNA was designed to specifically target circNFIA (A left) or circZNF561 (B left). The efficiency of sicircNFIA (A right) and circZNF561 (B right) was confirmed by qPCR. (CF) Cell proliferation analysis was performed with a CCK-8 assay in LNCaP and DU145 cells after knockdown of circNFIA and circZNF561. The p-value was calculated by the unpaired two-tailed Student’s t-test. * p < 0.05; ** p < 0.01; *** p < 0.001 (circRNAs).
Figure 7. The biological function of circNFIA and circZNF561. (A,B) siRNA was designed to specifically target circNFIA (A left) or circZNF561 (B left). The efficiency of sicircNFIA (A right) and circZNF561 (B right) was confirmed by qPCR. (CF) Cell proliferation analysis was performed with a CCK-8 assay in LNCaP and DU145 cells after knockdown of circNFIA and circZNF561. The p-value was calculated by the unpaired two-tailed Student’s t-test. * p < 0.05; ** p < 0.01; *** p < 0.001 (circRNAs).
Life 11 01096 g007
Figure 8. The biological function of circNFIA and circZNF561. (AD) Cell cycle assay was performed in DU145 and LNCaP cells after knockdown of circNFIA and circZNF561.
Figure 8. The biological function of circNFIA and circZNF561. (AD) Cell cycle assay was performed in DU145 and LNCaP cells after knockdown of circNFIA and circZNF561.
Life 11 01096 g008
Figure 9. The biological function of circNFIA and circZNF561. (AD) Cell apoptosis assays were performed in DU145 and LNCaP cells after knockdown of circNFIA and circZNF561. The data are presented as the mean ± SD (n = 3). The p-value was calculated by the unpaired two-tailed Student’s t-test. * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 9. The biological function of circNFIA and circZNF561. (AD) Cell apoptosis assays were performed in DU145 and LNCaP cells after knockdown of circNFIA and circZNF561. The data are presented as the mean ± SD (n = 3). The p-value was calculated by the unpaired two-tailed Student’s t-test. * p < 0.05; ** p < 0.01; *** p < 0.001.
Life 11 01096 g009
Figure 10. Construction of the circNFIA/circZNF561-mediated ceRNA network in prostate cancer. (A,E) Prediction of circNFIA- and circZNF561-targeted miRNAs by circBank. (B,F) Clusters of circNFIA- or circZNF561-mediated ceRNA networks in prostate cancer. The networks were drawn using Cytoscape 3.0. (C,D) GO and KEGG pathway analysis of target genes in circNFIA-mediated ceRNA networks. (G,H) GO and KEGG pathway analysis of target genes in circZNF561-mediated ceRNA networks.
Figure 10. Construction of the circNFIA/circZNF561-mediated ceRNA network in prostate cancer. (A,E) Prediction of circNFIA- and circZNF561-targeted miRNAs by circBank. (B,F) Clusters of circNFIA- or circZNF561-mediated ceRNA networks in prostate cancer. The networks were drawn using Cytoscape 3.0. (C,D) GO and KEGG pathway analysis of target genes in circNFIA-mediated ceRNA networks. (G,H) GO and KEGG pathway analysis of target genes in circZNF561-mediated ceRNA networks.
Life 11 01096 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kong, Z.; Lu, Y.; Wan, X.; Luo, J.; Li, D.; Huang, Y.; Wang, C.; Li, Y.; Xu, Y. Comprehensive Characterization of Androgen-Responsive circRNAs in Prostate Cancer. Life 2021, 11, 1096. https://doi.org/10.3390/life11101096

AMA Style

Kong Z, Lu Y, Wan X, Luo J, Li D, Huang Y, Wang C, Li Y, Xu Y. Comprehensive Characterization of Androgen-Responsive circRNAs in Prostate Cancer. Life. 2021; 11(10):1096. https://doi.org/10.3390/life11101096

Chicago/Turabian Style

Kong, Zhe, Yali Lu, Xuechao Wan, Jun Luo, Dujian Li, Yan Huang, Chenji Wang, Yao Li, and Yaoting Xu. 2021. "Comprehensive Characterization of Androgen-Responsive circRNAs in Prostate Cancer" Life 11, no. 10: 1096. https://doi.org/10.3390/life11101096

APA Style

Kong, Z., Lu, Y., Wan, X., Luo, J., Li, D., Huang, Y., Wang, C., Li, Y., & Xu, Y. (2021). Comprehensive Characterization of Androgen-Responsive circRNAs in Prostate Cancer. Life, 11(10), 1096. https://doi.org/10.3390/life11101096

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop