Abstract
Circular RNAs (circRNAs), a recently discovered non-coding RNA, have a number of functions including the regulation of miRNA expression. They have been detected in a number of malignancies including prostate cancer (PCa). The differential expression pattern of circRNAs associated with PCa and androgen receptor (AR) status was investigated in this study. circRNA profiling was performed using a high throughout microarray assay on a panel of prostate cell lines, which consisted of normal, benign, and malignant cells (n = 9). circRNAs were more commonly significantly up-regulated (p < 0.05) than downregulated in malignant cell lines (n = 3,409) vs. benign cell lines (n = 2,949). In a grouped analysis based on AR status, there were 2,127 down-regulated circRNAs in androgen independent cell lines compared to 2,236 in androgen dependent cell lines, thus identifying a potential circRNA signature reflective of androgen dependency. Through a bioinformatics approach, the parental genes associated with the top 10 differentially expressed circRNAs were identified such as hsa_circ_0064644, whose predicted parental gene target is RBMS3, and hsa_circ_0060539, whose predicted gene target is SDC4. Furthermore, we identified three circRNAs associated with the parental gene Caprin1 (hsa_circ_0021652, hsa_circ_0000288, and hsa_circ_0021647). Other studies have shown the importance of Caprin1 in PCa cell survival and drug resistance. Given the modified circRNA expression signatures identified here, these hypothesis generating results suggest that circRNAs may serve as potential putative diagnostic and predictive markers in PCa. However, further validation studies are required to assess the true potential of these markers in the clinical setting.
Introduction
Non-coding RNAs (ncRNAs), which include microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) (; ; ), play an important role in gene regulation (; ; ). circular RNAs (circRNAs), a recently discovered type of ncRNA, are generated from the backsplicing of exons, introns, or both to form exonic or intronic circRNAs (). circRNAs are covalently joined by their 3′- and 5′- ends, which are formed by back-splice events, thus presenting as closed continuous structures, making them highly stable and resistant to degradation (; ). circRNAs have many postulated functions such as the regulation of miRNA function by controlling the expression of miRNAs through a “sponging effect,” however this appears to apply to only a small number of circRNAs, such as ciRS-7 acting as a miR-7 sponge (). It has also been proposed that circRNAs may have a role in protein synthesis, with a number of circRNAs implicated in the translation of peptides, such as ZNF609 (). As research in this space increases, it is now clear that the circRNA-miRNA-mRNA network, plays an important role in both gene regulation and carcinogenesis (). circRNAs have also been detected in prostate cancer, which make them an attractive research target (). PCa growth and development is primarily dependent on the androgen receptor (AR), which is the target of therapeutic agents such as enzalutamide and abiraterone (). AR copy number gain and generation of variants such as AR-V7 are associated with the development of castration-resistant disease and drug resistance, however, despite significant advances; no suitable assay has become routinely available in the clinic to identify resistance mediators of AR targeting therapeutics (; ). Currently, PSA is used to screen men for PCa, yet when used in isolation, it is sub-optimal given issues with specificity and sensitivity (). Similarly, PSA cannot reliably detect early relapse or resistance to drug treatments, nonetheless earlier detection of castration-resistant disease would allow for more appropriate management of these patients and therefore significantly improve outcomes (). Thus, there is a need to identify new diagnostic and predictive biomarkers, as well as new therapeutic agents for use in the clinical setting. circRNAs may fulfill these roles, given their expression in PCa, in addition to studies showing circRNA such as circFOXO3 (hsa_circ_0006404) acting as an miRNA sponge in this disease ().
The aim of this study was to profile circRNAs in a panel of prostate cell lines and identify potential signatures associated with malignancy and androgen dependency. We propose that circRNAs have the potential to serve as useful biomarkers to improve diagnostic screening for PCa and/or identify men who are at risk of developing castration-resistant disease.
Materials and Methods
Cell Lines
A panel of malignant and benign prostate cell lines were used in this project, all of which were purchased from the ATCC (LGC standards, Middlesex, United Kingdom). The benign prostate cell lines (PWR-1E, RWPE-1) were cultured in Keratinocyte Serum Free Media (Thermo Fisher Scientific, CA, US), containing 0.05 mg/mL bovine pituitary extract (BPE) and 5 ng/mL Epidermal growth factor (EGF) (Thermo Fisher Scientific) and 1% Penicillin Streptomycin (P/S) (Merck KGaA, Darmstadt, Germany). The benign prostatic hyperplasia (BPH-1) and PCa cell lines (22Rv1, VCaP, LNCaP, DU145) were cultured in RPMI-1640 (Merck KGaA) supplemented with 10% Fetal Bovine Serum (FBS) (Merck KGaA) and 1% P/S. PC-3 was cultured in ATCC-formulated F-12K Medium (Thermo Fisher Scientific) supplemented with 10% FBS and 1% P/S.
RNA Preparation
Total RNA was prepared from cell lines (biological triplicates) according to manufacturer’s instructions. Briefly, cells were lysed directly with 2 mL TRIzol® Reagent (Thermo Fisher Scientific). For phase separation, 200 μL chloroform (Merck KGaA) was added to the cell lysate. The sample were vortexed vigorously for 15 s, incubated at RT for 5 min and centrifuged at 12,000 × g for 15 min at 4°C. The upper aqueous phase was transferred into a fresh tube and the RNA precipitated from the aqueous phase, by the addition of 400 μL isopropyl alcohol (Merck KGaA). The sample was mixed, incubated for 5 min at RT and centrifuged at 12,000 × g for 15 min at 4°C. One milliliter 75% EtOH (Merck KGaA) was added to pellet, and centrifuged at 12,000 × g for 10 min at 4°C. The EtOH was removed and a repeat pulse centrifuge was performed. Surplus EtOH was removed and the pellet air-dried. The sample was re-suspended in 30 μL molecular grade H2O before DNase treatment with Ambion® TURBOTM DNase (Thermo Fisher Scientific), and an additional RNA clean-up using standard EtOH precipitation.
CircRNA Microarray
The panel of cell lines (three biological replicates for each) were profiled using the Arraystar Human circRNA Array version 2.0 (Arraystar, MD, United States). The sample preparation and microarray hybridization were performed according to manufacturer’s instructions. Briefly, total RNA was digested with RNAse R (Epicenter, Illumina, CA, United States) to remove linear RNAs and enrich for circRNAs. The enriched circRNAs were amplified and transcribed into fluorescent cRNA using a random priming method with Arraystar Super RNA Labeling Kit (Arraystar). The labeled cRNAs were hybridized onto the Arraystar Human circRNA Array V2 (8 × 15 K). The array slides were washed and scanned on the Agilent Scanner G2505C. Agilent Feature Extraction software (version 11.0.1.1) was used to analyze acquired array images.
Microarray Data Analysis
Quantile normalization and subsequent data processing were executed using the R software package (). circRNAs with at least 4 out of 8 samples that were flagged as present or marginal (an attribute that denotes the quality of the entities) were considered as target circRNAs according to GeneSpring software’s definitions and instructions. circRNA and miRNA interactions were predicted with the Arraystar’s miRNA target prediction software based on TargetScan () and miRanda (). These target circRNAs were used for further differential analysis. A student’s paired t-test was used to identify significantly altered circRNAs. The false discovery rate (FDR) was applied to determine the threshold of p-value. An FDR of < 0.05 was used. Changes in expression were identified using unsupervised clustering analysis (euclidean distance measure and the “average” agglomeration method) and associated heat-maps generated.
Results
CircRNAs Are Differentially Expressed in Malignancy
The human circular RNA microarray (Arraystar) version 2.0 covers 13,617 previously discovered human circRNAs. Our cell line panel consisted of malignant (22Rv1, LNCaP, VCAP, DU145, and PC-3), benign (BPH-1), and normal (PWR1, RWPE-1) cell lines. AR expression in 22Rv1, LNCaP, and VCaP was confirmed using qPCR (Supplementary Methods and Supplementary Figure S1). In total 9,805 circRNAs were classified as present across the complete panel of cell lines. Initially, circRNA expression levels between cell lines based on malignancy was investigated. Differentially expressed circRNAs were examined by computing the fold change (FC) (i.e., the ratio of the group averages) for each circRNA between malignant, normal and benign cell lines. Associated parental genes for the top 10 altered circRNAs were predicted from circBase, which uses the co-ordinates from each detected circRNA to identify the associated gene from the UCSC genome database ().
CircRNAs were significantly differentially expressed between these groups (BPH vs. normal; malignant vs. BPH/normal) according to FC and p-value (Figure 1). We first examined the difference in expression of circRNAs between normal cell lines (PWR1, RWPE-1) and the benign cell line, BPH-1. CircRNAs were more often down-regulated in normal cell lines compared to the benign cell line (5,429 vs. 4,344). Next, we examined the expression of circRNAs between the normal/benign cells compared to malignant cells. In total, a higher number of circRNAs were down-regulated in normal/benign cells compared to malignant cells (5,300 vs. 4,457). The top 10 down and up-regulated circRNAs according to FC, along with their associated parental gene, are given in Tables 1, 2, respectively.
FIGURE 1
TABLE 1
| circRNA | FC | Chromosome | Parental gene |
| hsa_circ_0064644 | 5.1 | chr3 | RBMS3 () |
| hsa_circ_0082672 | 5.1 | chr7 | HIPK2 () |
| hsa_circ_0022382 | 3.3 | chr11 | FADS2 |
| hsa_circ_0049657 | 3.0 | chr19 | NFIX () |
| hsa_circ_0005402 | 3.0 | chr15 | ANXA2 () |
| hsa_circ_0082680 | 2.9 | chr7 | HIPK2 () |
| hsa_circ_0001187 | 2.5 | chr21 | DOPEY2 |
| hsa_circ_0082546 | 2.5 | chr7 | MTPN |
| hsa_circ_0022383 | 2.5 | chr11 | FADS2 |
| hsa_circ_0056731 | 2.4 | chr2 | CACNB4 |
Top 10 down-regulated circRNAs in PCa (malignant vs. normal/benign cell lines)*.
*References refer to the previous identification of these genes in cancer.
TABLE 2
| circRNA | FC | Chromosome | Parental gene |
| hsa_circ_0060539 | 7.8 | chr20 | SDC4 () |
| hsa_circ_0001589 | 6.3 | chr6 | HIST1H1D |
| hsa_circ_0026457 | 6.0 | chr12 | KRT5 () |
| hsa_circ_0043570 | 5.0 | chr17 | TNS4 () |
| hsa_circ_0060540 | 5.0 | chr20 | SDC4 () |
| hsa_circ_0008805 | 4.8 | chr17 | ARHGAP23 |
| hsa_circ_0084021 | 3.4 | chr8 | PLEKHA2 () |
| hsa_circ_0026358 | 3.2 | chr12 | KRT7 () |
| hsa_circ_0007058 | 3.2 | chr11 | TEAD1 () |
| hsa_circ_0004405 | 2.8 | chr5 | FAM169A |
Top 10 up-regulated circRNA in malignant vs. benign/normal cell lines*.
*References refer to the previous identification of these genes in cancer.
A number of circRNAs, that could play a role in PCa, were identified through their interaction with their associated parental genes. hsa_circ_0064644 was the highest down-regulated circRNA (FC 5, p = 0.001), located on chromosome 3 (Table 1). Its’ predicted parental gene target using circBase is RBMS3, which has previously been identified to suppress cell proliferation, migration, invasion, and angiogenesis in a number of different cancers (). RBMS3 has also been shown to act as a miRNA sponge in PCa via lncRNAs. Other genes of interest include HIPK2 (hsa_circ_0082672) and NFIX (hsa_circ_0049657), both known to be involved in the regulation of the AR (). hsa_circ_0005402, is linked to ANXA2, which has been suggested to play a role in epithelial-mesenchymal transition (EMT) in PCa () via transcriptional repression by ERG (). hsa_circ_0060539, located on chromosome 20, was the highest up-regulated circRNA. The predicted gene target for this circRNA is SDC4, a known active gene in PCa (). Similarly, KRT5 (hsa_circ_0026457) has been shown to be involved in the inhibition of PCa stem cell invasion, migration and proliferation, as has KRT7 (hsa_circ_0026358), which has been identified along with FAM129A, MME and SOD2 as part of a four-gene androgen regulated panel ().
CircRNAs Are Differentially Expressed According to Androgen Dependence
Lastly, the differing expression of circRNAs in relation to androgen sensitivity was examined, which we hypothesized would reflect the development of castration resistant PCa. Cell lines were stratified into two groups based on AR status: androgen dependent/castration-sensitive (22Rv1, LNCaP, VCaP) and androgen independent/castration resistant) cell lines (DU145, PC-3). Response to androgen ablation was demonstrated by confirming the presence or absence of AR-FL across cell lines (see Supplementary Figure S1). CircRNAs were more often down-regulated in cell lines that were androgen independent compared to androgen-dependent cell lines (5,067 down-regulated vs. 4,693 up-regulated, respectively) (Figure 2). The top 10 down and up-regulated circRNAs and their associated parental genes are given in Tables 3, 4, respectively. The up-regulated hsa_circ_0000361 (FC 19.2, p < 0.001) is an exonic circRNA located on chromosome 11 and is associated with the gene SIAE (Table 3). mRNA expression of SIAE has previously been shown to be reduced in patients with cancer (). We identified three circRNAs associated with the parental gene Caprin1 (hsa_circ_0021652, hsa_circ_0000288, and hsa_circ_0021647). Caprin1 plays an important role in PCa cell survival and its upregulation has been shown to be associated with drug resistance, and is increased in SPOP mutated cell lines ().
FIGURE 2
TABLE 3
| circRNA | FC | Chromosome | Parental gene |
| hsa_circ_0000367 | 19.2 | chr11 | SIAE () |
| hsa_circRNA_404686 | 18.5 | chr1 | GPR137B () |
| hsa_circ_0004183 | 17.7 | chr10 | FRMD4A () |
| hsa_circ_0001666 | 13.6 | chr6 | FAM120B |
| hsa_circRNA_400850 | 13.3 | chr11 | TRAPPC4 |
| hsa_circ_0021652 | 13.0 | chr11 | CAPRIN1 () |
| hsa_circ_0000288 | 12.7 | chr11 | CAPRIN1 () |
| hsa_circ_0021647 | 12.4 | chr11 | CAPRIN1 () |
| hsa_circ_0006220 | 11.4 | chr17 | TADA2A |
| hsa_circ_0043278 | 11.4 | chr17 | TADA2A |
Top 10 up-regulated circRNAs in androgen dependent vs. independent cell lines*.
*References refer to the previous identification of these genes in cancer.
TABLE 4
| circRNA | FC | Chromosome | Parental gene |
| hsa_circ_0000825 | 40.6 | chr18 | MTCL1 |
| hsa_circ_0084615 | 24.1 | chr8 | ASPH () |
| hsa_circ_0002083 | 18.2 | chr11 | ETS1 () |
| hsa_circ_0091934 | 11.7 | chrX | FLNA () |
| hsa_circ_0091894 | 10.3 | chrX | FLNA () |
| hsa_circ_0085923 | 9.1 | chr8 | PLEC () |
| hsa_circ_0007386 | 8.8 | chr2 | CRIM1 () |
| hsa_circ_0044468 | 8.6 | chr17 | ITGA3 () |
| hsa_circRNA_402986 | 7.9 | chr3 | PLOD2 () |
| hsa_circ_0025506 | 7.5 | chr12 | GPRC5A () |
Top 10 down-regulated circRNAs in androgen dependent vs. independent cell lines*.
*References refer to the previous identification of these genes in cancer.
The highest down-regulated circRNA by FC was hsa_circ_0000825 (FC 40, p < 0.0001) located on chromosome 18 and is associated with the gene MTCL1 (Table 4). This gene has previously been identified to play a role in the development of castration-resistant PCa (). Other parental genes of interest include FLNA (hsa_circ_0091934 and hsa_circ_0091894) located on the X chromosome, which plays a significant role in PCa development and progression, and has been highlighted as a possible biomarker for disease screening and detection ().
Discussion
CircRNAs are a novel type of ncRNA which have been identified across a range of cancers including PCa (; ). The most abundant class of ncRNA are lncRNAs which have been extensively investigated in PCa (). A number of lncRNAs have been identified as potential predictive and prognostic biomarkers including SChLAP1, UCA1 and PCAT14 (; ; ). Further work is exploring lncRNAs as therapeutic targets in PCa using small interfering RNAs (siRNAs), antisense oligonucleotides (ASOs) and small molecule inhibitors (). Therefore, there is significant interest in the role ncRNAs may play in the future management of PCa. The exact role of circRNAs in cancer has yet to be fully elucidated, however initial studies suggest circRNAs may play a role in cancer development, progression and the development of resistance to therapeutic agents (). CircRNAs have been identified as potential biomarkers in cancer due to their abundance, stability and disease specific activity (). The identification of circRNAs has been complicated by bioinformatics challenges, however the development of improved detection methodologies and statistical methods has aided circRNA research, particularly with improvements in reducing false-positive detection rates (). In PCa, a number of circRNAs have been identified which have been identified as potential biomarkers (). Our laboratory previously identified the circRNA, hsa_circ_0004870, to be associated with resistance to the androgen receptor inhibitor enzalutamide (). Other circRNAs of note that have been identified in PCa include circAMACR, circAURKA and circAR3 which is encoded by the AR gene ().
We identified a number of circRNAs that could represent a unique gene signature of aggressive PCa. Initial bioinformatics analysis identified a significant number of differentially expressed circRNAs between malignant vs. benign/normal cell lines. There was a trend for circRNAs to be down-regulated in normal compared to benign cells. This trend continued when we compared normal/benign cells to malignant cells. When we examined circRNA expression in cell lines according to androgen dependency, we found circRNAs again to be more often down-regulated in androgen independent cell lines, reflective of aggressive PCa.
Parental genes of interest were RBMS3 (hsa_circ_0064644) and SDC4 (hsa_circ_0060539). Interestingly, we identified three circRNAs associated with the parental gene Caprin1 (hsa_circ_0021652, hsa_circ_0000288, and hsa_circ_0021647) and two circRNAs associated with FLNA (hsa_circ_0091934 and hsa_circ_0091894). Many of these genes are currently under active investigation as biomarkers in PCa, and identifying their associated circRNAs may contribute to ongoing research in this field (; ; ; ).
There is evidence to suggest that circRNAs regulate miRNA function by competing for the pool of miRNA binding sites to influence the activities of miRNAs in regulating gene expression (). Interestingly, a number of miRNAs that have been detected in PCa are linked to the AR pathway (; ). miRNA binding sites can be predicted using a bioinformatics pipeline for each detected circRNA that identified an associated gene. For instance, for hsa_circ_0064644, the predicted miRNA is miR-15a, which has been show to play a role in PCa progression and metastases (). Therefore, predicted miRNAs could be used to identify circRNAs that play a role in gene regulation in PCa, and further our understanding of the importance of the circRNA-miRNA-mRNA network in this disease.
This study had several limitations. Firstly, we used a circRNA microarray platform to detect circRNAs compared to RNA-seq methods which may have higher sensitivity and accuracy to detect circRNAs. Secondly, this was a preliminary in vitro study which identified multiple potential genes and ongoing validation is required in cell lines and in vivo to confirm our findings here.
Conclusion
In conclusion, circRNAs were modified between benign and malignant cell lines, as well as androgen dependent and independent cell lines, which suggests a role for circRNAs in PCa initiation and progression. This study highlights the importance of circRNAs in PCa and their potential role in identifying the disease and monitoring response to treatment. The dysregulated signature identified here may prove useful for the development of a blood-based assay with both diagnostic and predictive value in this disease. However, this would require validation in larger patient cohorts.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.605686/full#supplementary-material
References
1
AgarwalV.BellG. W.NamJ.-W.BartelD. P. (2015). Predicting effective microRNA target sites in mammalian mRNAs.eLife4:e05005.
2
AntonarakisE. S.LuC.WangH.LuberB.NakazawaM.RoeserJ. C.et al (2014). AR-V7 and resistance to enzalutamide and abiraterone in prostate cancer.New Engl. J. Med.3711028–1038.
3
BahnJ. H.ZhangQ.LiF.ChanT. M.LinX.KimY.et al (2015). The landscape of microRNA, Piwi-interacting RNA, and circular RNA in human saliva.Clin. Chem.61221–230. 10.1373/clinchem.2014.230433
4
BarboroP.BenelliR.TosettiF.CostaD.CapaiaM.AstigianoS.et al (2020). Aspartate β-hydroxylase targeting in castration resistant prostate cancer modulates the NOTCH/HIF1α/GSK3β crosstalk.Carcinogenesis411246–1252. 10.1093/carcin/bgaa053
5
BarrettS. P.WangP. L.SalzmanJ. (2015). Circular RNA biogenesis can proceed through an exon-containing lariat precursor.eLife4:e07540.
6
BiiV. M.CollinsC. P.HocumJ. D.TrobridgeG. D. (2018). Replication-incompetent gammaretroviral and lentiviral vector-based insertional mutagenesis screens identify prostate cancer progression genes.Oncotarget915451–15463. 10.18632/oncotarget.24503
7
BuddW. T.Seashols-WilliamsS. J.ClarkG. C.WeaverD.CalvertV.PetricoinE.et al (2015). Dual Action of miR-125b as a tumor suppressor and OncomiR-22 promotes prostate cancer tumorigenesis.PLoS One10:e0142373. 10.1371/journal.pone.0142373
8
ChenS.HuangV.XuX.LivingstoneJ.SoaresF.JeonJ.et al (2019). Widespread and functional RNA circularization in localized prostate cancer.Cell176831–843.e22.
9
Di NunnoV.GattoL.SantoniM.CimadamoreA.Lopez-BeltranA.ChengL.et al (2018). Recent advances in liquid biopsy in patients with castration resistant prostate cancer.Front. Oncol.8:397. 10.3389/fonc.2018.00397
10
DongY.HeD.PengZ.PengW.ShiW.WangJ.et al (2017). Circular RNAs in cancer: an emerging key player.J. Hematol. Oncol.10:2.
11
DuH.WangX.DongR.HuD.XiongY. (2019). miR-601 inhibits proliferation, migration and invasion of prostate cancer stem cells by targeting KRT5 to inactivate the Wnt signaling pathway.Int. J. Clin. Exp. Pathol.124361–4379.
12
DuY.BiW.ZhangF.WuW.XiaS.LiuH. G. - (2015). protein-coupled receptor 137 accelerates proliferation of urinary bladder cancer cells in vitro.Biotechnol. Appl. Biochem.62855–860. 10.1002/bab.1329
13
EstellerM. (2011). Non-coding RNAs in human disease.Nat. Rev. Genet.12861–874.
14
FabrisL.CederY.ChinnaiyanA. M.JensterG. W.SorensenK. D.TomlinsS.et al (2016). The potential of MicroRNAs as prostate cancer biomarkers.Eur. Urol.70312–322.
15
Fotouhi GhiamA.TaebS.HuangX.HuangV.RayJ.ScarcelloS.et al (2017). Long non-coding RNA urothelial carcinoma associated 1 (UCA1) mediates radiation response in prostate cancer.Oncotarget84668–4689. 10.18632/oncotarget.13576
16
FuD.ShiY.LiuJ.-B.WuT.-M.JiaC.-Y.YangH.-Q.et al (2020). Targeting long Non-coding RNA to therapeutically regulate gene expression in cancer.Mol. Ther. Nucleic Acids21712–724. 10.1016/j.omtn.2020.07.005
17
GlazarP.PapavasileiouP.RajewskyN. (2014). CircBase: a database for circular RNAs.RNA201666–1670. 10.1261/rna.043687.113
18
GrabowskaM. M.ElliottA. D.DeGraffD. J.AndersonP. D.AnumanthanG.YamashitaH.et al (2014). NFI transcription factors interact with FOXA1 to regulate prostate-specific gene expression.Mol. Endocrinol.28949–964. 10.1210/me.2013-1213
19
GreeneJ.BairdA. M.BradyL.LimM.GrayS. G.McDermottR.et al (2017). Circular RNAs: biogenesis, function and role in human diseases.Front. Mol. Biosci.4:38. 10.3389/fmolb.2017.00038
20
GreeneJ.BairdA. M.CaseyO.BradyL.BlackshieldsG.LimM.et al (2019). Circular RNAs are differentially expressed in prostate cancer and are potentially associated with resistance to enzalutamide.Sci. Rep.9:10739.
21
GrinerN. B.YoungD.ChaudharyP.MohamedA. A.HuangW.ChenY.et al (2015). ERG oncoprotein inhibits ANXA2 expression and function in prostate cancer.Mol. Cancer Res.13368–379. 10.1158/1541-7786.mcr-14-0275-t
22
HansenT. B.KjemsJ.DamgaardC. K. (2013). Circular RNA and miR-7 in cancer.Cancer Res.735609–5612.
23
HeJ. Y.WeiX. H.LiS. J.LiuY.HuH. L.LiZ. Z.et al (2018). Adipocyte-derived IL-6 and leptin promote breast cancer metastasis via upregulation of Lysyl Hydroxylase-2 expression.Cell Commun. Signal.16:100.
24
HudsonB. D.HumN. R.ThomasC. B.KohlgruberA.SebastianA.ColletteN. M.et al (2015). SOST inhibits prostate cancer invasion.PLoS One10:e0142058. 10.1371/journal.pone.0142058
25
Imberg-KazdanK.HaS.GreenfieldA.PoultneyC. S.BonneauR.LoganS. K.et al (2013). A genome-wide RNA interference screen identifies new regulators of androgen receptor function in prostate cancer cells.Genome Res.23581–591. 10.1101/gr.144774.112
26
JeckW. R.SorrentinoJ. A.WangK.SlevinM. K.BurdC. E.LiuJ.et al (2013). Circular RNAs are abundant, conserved, and associated with ALU repeats.RNA19141–157. 10.1261/rna.035667.112
27
JiangZ.ZhangY.ChenX.WuP.ChenD. (2020). Long noncoding RNA RBMS3-AS3 acts as a microRNA-4534 sponge to inhibit the progression of prostate cancer by upregulating VASH1.Gene Ther.27143–156. 10.1038/s41434-019-0108-1
28
JinW.ChenF.WangK.SongY.FeiX.WuB. (2018). miR-15a/miR-16 cluster inhibits invasion of prostate cancer cells by suppressing TGF-β signaling pathway.Biomed. Pharmacother.104637–644. 10.1016/j.biopha.2018.05.041
29
JohnB.EnrightA. J.AravinA.TuschlT.SanderC.MarksD. S. (2004). Human MicroRNA targets.PLoS Biol.2:e363. 10.1371/journal.pbio.0020363
30
KnightJ. F.ShepherdC. J.RizzoS.BrewerD.JhavarS.DodsonA. R.et al (2008). TEAD1 and c-Cbl are novel prostate basal cell markers that correlate with poor clinical outcome in prostate cancer.Br. J. Cancer991849–1858. 10.1038/sj.bjc.6604774
31
KongZ.WanX.LuY.ZhangY.HuangY.XuY.et al (2020). Circular RNA circFOXO3 promotes prostate cancer progression through sponging miR-29a-3p.J. Cell. Mol. Med.24799–813. 10.1111/jcmm.14791
32
KurozumiA.GotoY.MatsushitaR.FukumotoI.KatoM.NishikawaR.et al (2016). Tumor-suppressive microRNA-223 inhibits cancer cell migration and invasion by targeting ITGA3/ITGB1 signaling in prostate cancer.Cancer Sci.10784–94. 10.1111/cas.12842
33
LeblancR.SahayD.HoussinA.Machuca-GayetI.PeyruchaudO. (2018). Autotaxin-β interaction with the cell surface via syndecan-4 impacts on cancer cell proliferation and metastasis.Oncotarget933170–33185. 10.18632/oncotarget.26039
34
LegniniI.Di TimoteoG.RossiF.MorlandoM.BrigantiF.SthandierO.et al (2017). Circ-ZNF609 Is a circular RNA that can be translated and functions in myogenesis.Mol. Cell6622–37.e9.
35
LimM. C. J.BairdA.-M.AirdJ.GreeneJ.KapoorD.GrayS. G.et al (2018). RNAs as candidate diagnostic and prognostic markers of prostate cancer—from cell line models to liquid biopsies.Diagnostics8:60. 10.3390/diagnostics8030060
36
LinH.-M.MahonK. L.SpielmanC.GurneyH.MallesaraG.StocklerM. R.et al (2017). Phase 2 study of circulating micicroRNA biomarkers in castration-resistant prostate cancer.Br. J. Cancer116:1002. 10.1038/bjc.2017.50
37
MartinyP. B.AlcobaD. D.NetoB. S.CarvalhoP. C.BrumI. S. A. (2018). Proteomic glimpse into the oncogenesis of prostate cancer.J. Appl. Biomed.16328–336. 10.1016/j.jab.2018.05.003
38
MehraR.UdagerA. M.AhearnT. U.CaoX.FengF. Y.LodaM.et al (2016). Overexpression of the long Non-coding RNA SChLAP1 independently predicts lethal prostate cancer.Eur. Urol.70549–552. 10.1016/j.eururo.2015.12.003
39
MisawaA.TakayamaK. I.InoueS. (2017). Long non-coding RNAs and prostate cancer.Cancer Sci.1082107–2114. 10.1111/cas.13352
40
MuharramG.SahgalP.KorpelaT.De FranceschiN.KaukonenR.ClarkK.et al (2014). Tensin-4-dependent MET stabilization is essential for survival and proliferation in carcinoma cells.Dev. Cell.29421–436. 10.1016/j.devcel.2014.03.024
41
OhS. Y.KangS. M.KangS. H.LeeH. J.KwonT. G.KimJ. W.et al (2020). Potential salivary mRNA biomarkers for early detection of oral cancer.J. Clin. Med.9:243. 10.3390/jcm9010243
42
O’SullivanJ. (2017). Controversies in PSA screening.Evid. Based Med.22:198. 10.1136/ebmed-2017-110858
43
R Core Team (2020). The R Project for Statistical Computing.Vienna: R Core Team.
44
RamnarineV. R.KobelevM.GibbE. A.NouriM.LinD.WangY.et al (2019). The evolution of long noncoding RNA acceptance in prostate cancer initiation, progression, and its clinical utility in disease management.Eur. Urol.76546–559. 10.1016/j.eururo.2019.07.040
45
RavipatyS.WuW.DalviA.TannaN.AndreaziJ.FrissT.et al (2017). Clinical validation of a serum protein panel (FLNA, FLNB and KRT19) for diagnosis of prostate cancer.J. Mol. Biomark. Diagn.8:323.
46
RodgersJ. J.McClureR.EpisM. R.CohenR. J.LeedmanP. J.HarveyJ. M.et al (2019). ETS1 induces transforming growth factor β signaling and promotes epithelial-to-mesenchymal transition in prostate cancer cells.J. Cell. Biochem.120848–860. 10.1002/jcb.27446
47
RomanelA.Gasi TandefeltD.ConteducaV.JayaramA.CasiraghiN.WetterskogD.et al (2015). Plasma AR and abiraterone-resistant prostate cancer.Sci. Transl. Med.7:312re10.
48
SalzmanJ.GawadC.WangP. L.LacayoN.BrownP. O. (2012). Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types.PLoS One7:e30733. 10.1371/journal.pone.0030733
49
SawadaY.KikugawaT.IioH.SakakibaraI.YoshidaS.IkedoA.et al (2020). GPRC5A facilitates cell proliferation through cell cycle regulation and correlates with bone metastasis in prostate cancer.Intern. J. Cancer J. Intern. Cancer1461369–1382. 10.1002/ijc.32554
50
ShiQ.ZhuY.MaJ.ChangK.DingD.BaiY.et al (2019). Prostate cancer-associated SPOP mutations enhance cancer cell survival and docetaxel resistance by upregulating Caprin1-dependent stress granule assembly.Mol. Cancer18:170.
51
ShihJ. W.WangL. Y.HungC. L.KungH. J.HsiehC. L. (2015). Non-Coding RNAs in castration-resistant prostate cancer: regulation of androgen receptor signaling and cancer metabolism.Intern. J. Mol. Sci.1628943–28978. 10.3390/ijms161226138
52
SuM.XiaoY.MaJ.TangY.TianB.ZhangY.et al (2019). Circular RNAs in cancer: emerging functions in hallmarks, stemness, resistance and roles as potential biomarkers.Mol. Cancer18:90.
53
SzaboL.SalzmanJ. (2016). Detecting circular RNAs: bioinformatic and experimental challenges.Nat. Rev. Genet.17679–692. 10.1038/nrg.2016.114
54
ThomasB. C.KayJ. D.MenonS.VowlerS. L.DawsonS. N.BucklowL. J.et al (2016). Whole blood mRNA in prostate cancer reveals a four-gene androgen regulated panel.Endocr. Relat. Cancer23797–812. 10.1530/erc-16-0287
55
TuckerD.ZhengW.ZhangD.-H.DongX. (2020). Circular RNA and its potential as prostate cancer biomarkers.World J. Clin. Oncol.11563–572. 10.5306/wjco.v11.i8.563
56
WarnerJ. R. (1999). The economics of ribosome biosynthesis in yeast.Trends Biochem. Sci.24437–440. 10.1016/s0968-0004(99)01460-7
57
WhiteN. M.ZhaoS. G.ZhangJ.RozyckiE. B.DangH. X.McFaddenS. D.et al (2017). Multi-institutional analysis shows that low PCAT-14 expression associates with poor outcomes in prostate cancer.Eur. Urol.71257–266. 10.1016/j.eururo.2016.07.012
58
YangN.WangL.LiuJ.LiuL.HuangJ.ChenX.et al (2018). MicroRNA-206 regulates the epithelial-mesenchymal transition and inhibits the invasion and metastasis of prostate cancer cells by targeting Annexin A2.Oncol. Lett.158295–8302.
59
ZhangC.XiongJ.YangQ.WangY.ShiH.TianQ.et al (2018). Profiling and bioinformatics analyses of differential circular RNA expression in prostate cancer cells.Future Sci. OA4:FSOA340-FSOA.
60
ZhengX.JiaB.LinX.HanJ.QiuX.ChuH.et al (2016). FRMD4A: a potential therapeutic target for the treatment of tongue squamous cell carcinoma.Int. J. Mol. Med.381443–1449. 10.3892/ijmm.2016.2745
61
ZhuL.XiP.-W.LiX.-X.SunX.ZhouW.-B.XiaT.-S.et al (2019). The RNA binding protein RBMS3 inhibits the metastasis of breast cancer by regulating Twist1 expression.J. Exper. Clin. Cancer Res.38:105.
Summary
Keywords
prostate cancer, circRNA, non-coding RNA, biomarkers, androgen signaling
Citation
Greene J, Baird A-M, Lim M, Flynn J, McNevin C, Brady L, Sheils O, Gray SG, McDermott R and Finn SP (2021) Differential CircRNA Expression Signatures May Serve as Potential Novel Biomarkers in Prostate Cancer. Front. Cell Dev. Biol. 9:605686. doi: 10.3389/fcell.2021.605686
Received
12 September 2020
Accepted
05 February 2021
Published
25 February 2021
Volume
9 - 2021
Edited by
Mingli Liu, Morehouse School of Medicine, United States
Reviewed by
Zoran Culig, Innsbruck Medical University, Austria; Yasser Ged, Johns Hopkins Medicine, United States
Updates
Copyright
© 2021 Greene, Baird, Lim, Flynn, McNevin, Brady, Sheils, Gray, McDermott and Finn.
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: John Greene, greenejo@tcd.ie
This article was submitted to Molecular Medicine, a section of the journal Frontiers in Cell and Developmental Biology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.