Next Article in Journal
PARP Inhibitors and Myeloid Neoplasms: A Double-Edged Sword
Next Article in Special Issue
CircFISH: A Novel Method for the Simultaneous Imaging of Linear and Circular RNAs
Previous Article in Journal
Intraoperative HIFU Ablation of the Pancreas Using a Toroidal Transducer in a Porcine Model. The First Step towards a Clinical Treatment of Locally Advanced Pancreatic Cancer
Previous Article in Special Issue
The Role of Circular RNAs in DNA Damage Response and Repair
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Hypothesis

hsa_circ_0001275 Is One of a Number of circRNAs Dysregulated in Enzalutamide Resistant Prostate Cancer and Confers Enzalutamide Resistance In Vitro

1
Department of Histopathology and Morbid Anatomy, Trinity Translational Medicine Institute, Trinity College Dublin, D08 W9RT Dublin, Ireland
2
Department of Medical Oncology, St. James’s Hospital, D08 NHY1 Dublin, Ireland
3
Department of Medical Oncology, Tallaght University Hospital, D24 NR0A Dublin, Ireland
4
School of Medicine, Trinity Translational Medicine Institute, Trinity College Dublin, D08 W9RT Dublin, Ireland
5
Department of Medical Oncology, St. Vincent’s University Hospital, D04 YN26 Dublin, Ireland
6
CiberNed (Centro Investigacion Biomedica en Red Enfermedades Neurodegenerativas), IIBB, Rosello 161, 08036 Barcelona, Spain
7
Thoracic Oncology Research Group, Labmed Directorate, St. James’s Hospital, D08 RXOX Dublin, Ireland
8
School of Biological and Health Sciences, Technological University Dublin, D07 ADY7 Dublin, Ireland
9
Department of Histopathology, St. James’s Hospital, D08 X4RX Dublin, Ireland
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 25 November 2021 / Accepted: 7 December 2021 / Published: 20 December 2021
(This article belongs to the Special Issue Circular RNAs: New Insights into the Molecular Biology of Cancer)

Abstract

:

Simple Summary

Although newer generations of androgen deprivation therapy such as enzalutamide are providing hope, it is clinically challenging to deliver effective therapy to individuals with metastatic castrate-resistant prostate cancer. Between 20–40% of patients have intrinsic resistance to therapy and all patients will ultimately experience disease progression due to acquired resistance, which is a significant clinical dilemma. The aim of our study was to evaluate the role of circular RNAs (circRNAs) in enzalutamide-resistant prostate cancer as part of the effort to identify useful biomarkers for patient selection and potential new therapeutic targets. We confirmed that hsa_circ_0001275 was highly upregulated in an enzalutamide resistant cell line and demonstrated that its overexpression resulted in increased enzalutamide resistance. Our data showed that hsa_circ_0001275 was not expressed abundantly in patient plasma samples, however, a trend of expression was evident which paralleled disease activity indicating a possible association with enzalutamide resistance. Overall, we have provided evidence that hsa_circ_0001275 promotes enzalutamide resistance and thus may serve as a potential therapeutic target.

Abstract

Background: Enzalutamide is part of the treatment regimen for metastatic castration-resistant prostate cancer (MCRPC). However, both intrinsic and acquired resistance to the drug remain substantial clinical quandaries. circRNAs, a novel type of non-coding RNA, have been identified in a number of cancers including prostate cancer and have been associated with cancer development and progression. circRNAs have shown great potential as clinically useful blood-based ‘liquid biopsies’ and as therapeutic targets in prostate cancer. The aim of this study was to examine the role of circRNA transcripts in enzalutamide-resistant prostate cancer cells and assess their utility as biomarkers. Methods: An isogenic cell line model of enzalutamide resistance was subjected to circRNA microarray profiling. Several differentially expressed circRNAs, along with their putative parental genes were validated using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). circRNAs of interest were stably overexpressed in the control cell line and drug sensitivity was assessed using an ELISA-based proliferation assay. The candidate circRNA, hsa_circ_0001275, was measured in patient plasma samples using RT-droplet digital PCR (RT-ddPCR). Results: hsa_circ_0001275 and its parental gene, PLCL2, were significantly up-regulated in strongly resistant clones vs. control (p < 0.05). Overexpression of hsa_circ_0001275 in the control cell line resulted in increased resistance to enzalutamide (p < 0.05). While RT-ddPCR analysis of hsa_circ_0001275 expression in plasma samples of 44 clinical trial participants showed a trend that mirrored the stages of disease activity (as defined by PSA level), the association did not reach statistical significance. Conclusions: Our data suggest that increased levels of hsa_circ_0001275 contribute to enzalutamide resistance. hsa_circ_0001275 plasma expression showed a trend that mirrors the PSA level at specific disease time points, indicating that circRNAs mirror disease recurrence and burden and may be associated with enzalutamide resistance.

Graphical Abstract

1. Introduction

Prostate cancer is the most common invasive cancer diagnosed in men and the second leading cause of male cancer-associated mortality in the United States (US) [1]. The majority (79%) of patients with prostate cancer are diagnosed at a localised stage [2], which is potentially curable using radical prostatectomy, external beam radiotherapy (EBRT), or brachytherapy [3]. However, about 35% to 60% of men undergoing the initial curative effort will experience disease recurrence [4,5,6]. For these patients, androgen deprivation therapy (ADT) is the standard of treatment [7] but resistance inevitably develops resulting in castration-resistant prostate cancer (CRPC), which is currently incurable [8].
The treatment strategy for prostate cancer has evolved rapidly in the last decade due to greater availability and selection of therapeutic agents that improve overall survival (OS) of patients with CRPC such as second-generation antiandrogens (abiraterone and enzalutamide), chemotherapy (cabazitaxel), and radionuclide therapy (alpha-radium 223) [9,10,11,12]. The second-generation anti-androgens have stood at the forefront of this progress, however, despite the effectiveness of these new agents, 20% to 40% of patients will have intrinsic resistance to these therapies and all patients will eventually experience disease progression due to acquired resistance [9,13,14]. Additionally, recent phase III studies such as CHAARTED and LATITUDE have highlighted the importance of appropriate sequencing of treatment in optimizing clinical benefit and outcomes for patients with metastatic castrate sensitive prostate cancer (MCSPC) [15,16]. Even though an impressive therapeutic armamentarium for prostate cancer exists, there is a lack of effective clinical tools centred on patient selection making it difficult to provide the optimal therapy at each time point along the care pathway.
Our knowledge of the genomic and transcriptomic landscape of prostate and other cancers has broadened substantially due to improved sequencing technologies, in turn leading to the ability to study previously poorly defined biomarkers such as non-coding RNA (ncRNA) contributing to novel molecular biomarker discoveries and new therapeutic strategies [17,18]. ncRNA consists of several different classes including circular RNA (circRNA), microRNA (miRNA) and long ncRNA (lncRNA) [19]. circRNA is produced from the backsplicing of exons, introns, or both where the 3′ and 5′ ends are covalently closed into one continuous loop, which critically confers resistance to degradation by RNA exonucleases [20,21]. Thus, circRNA is present in the cytoplasm owing to its stability and can be partitioned into exosomes [22,23]. For these reasons, circRNAs are believed to be a novel stable plasma biomarker. Importantly, the triad of stability, abundance and evolutionary preservation between species indicate that circRNA may have a substantial regulatory function [24] and appear to have essential roles in cancer initiation, development and progression [25,26]. A number of circRNA biomarkers have been investigated in prostate cancer as summarised by Papatsirou et al. [27], however, there are few studies in the literature describing the role of circRNA in enzalutamide resistance in vitro [28,29]. For instance, Wu et al. demonstrated that the downregulation of circRNA17 resulted in the overexpression of AR-V7 leading to enzalutamide resistance in prostate cancer cell lines [28] and Greene et al. found that the downregulation of hsa_circ_0004870 may lead to enzalutamide resistance through the regulation of AR-V7 [29]. Enzalutamide competitively binds to the carboxy-terminal ligand domain (LBD) of AR [9], thus hindering AR translocation, AR cofactor recruitment, and AR binding to DNA [30]. Enzalutamide has been shown in previous phase 3 studies to prolong OS and progression-free survival (PFS) in patients who were chemotherapy naive [31] and in the cohort of patients who had chemotherapy [9].
Our group previously identified differential circRNAs expression in an enzalutamide resistance model, which consisted of age-matched LNCaP parental cells (Control) and two sub-lines displaying strong resistance (Clone 1) and weak resistance (Clone 9) [29]. In this current study, two highly expressed circRNAs (hsa_circ_0001275 and hsa_circ_0001721) in Clone 1 were investigated in vitro, and expression levels determined using RT-ddPCR in plasma samples of patients with MCRPC collected as part of the Phase II trial of radium-223 in combination with enzalutamide (CTRIAL-IE 13-21, NCT02225704) [32].

2. Materials and Methods

2.1. Cell Lines

Isogenic LNCaP clones were generated from single clones as described by Korpal et al. [33], which consisted of aged-matched LNCaP parental cells (Control) and two sub-lines displaying strong resistance (Clone 1) and weak resistance (Clone 9) to enzalutamide. These clones were gifted from Novartis (Dublin, Ireland). All clones were cultured in RPMI-1640 (Sigma-Aldrich, St. Louis, MO, USA) supplemented with 10% foetal bovine serum (FBS) (Sigma-Aldrich, St. Louis, MO, USA) and 1% penicillin, streptomycin, and amphotericin b (Sigma-Aldrich, St. Louis, MO, USA) at 37 °C in a humidified atmosphere containing 5% CO2. All cell lines were periodically examined for mycoplasma infection and experiments were undertaken at passages below 30. The STR profiles for the LNCaP Control, C1 and C9 cells were compared against the ATCC STR database (https://www.atcc.org/search-str-database, accessed on 19 November 2021). The matching criterion is based on an algorithm that compares the number of shared alleles between two cell line samples, expressed as a percentage. Cell lines with ≥80% match are considered to be related; derived from a common ancestor. Using the following STR criteria: (AMEL) + the following loci: D5S818, D13S317, D7S820, D16S539, vWA, TH01, TPOX, CSF1PO. The results for identity using the ATCC algorithm were as follows: LNCAP Control: 94% (Pass), LNCaP C1: 94% (Pass), LNCaP C9: 94% (Pass).

2.2. Total RNA Preparation and Reverse Transcription

The miRNeasy Mini Kit (Qiagen, Manchester, UK) was used to isolate RNA from cell lines according to the manufacturer’s instructions. RNA was isolated from 200 µL plasma using the miRNeasy Plasma/Serum Kit (Qiagen) according to the manufacturer’s instructions. The NanoDrop 1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) was used to quantify the RNA. cDNA was synthesized from either 1 µg (cell lines) or 0.1 µg (plasma) total RNA through random priming using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Waltham, MA, USA) according to the manufacturer’s instructions.

2.3. circRNA Enrichment Using RNase R

Degradation of linear RNA for circRNA enrichment by digestion using RNase R from Lucigen (LGC, Teddington, Middlesex, UK) was performed. Briefly, 2 µg RNA was digested using RNase R according to the protocol by Panda et al. [34] followed by immediate RNA isolation using the miRNeasy Kit (Qiagen) according to the manufacturer’s protocol. A control RNA reaction was also performed in parallel to the RNase R reaction but without Ribonuclease R. Both RNA with and without RNase R digestion were reverse transcribed to cDNA using High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems) and RT-qPCR was performed for circRNA analysis. Data were analysed according to the protocol by Panda et al. [34] (Supplementary Figure S1).

2.4. circRNA Microarray

The isogenic LNCaP clones (n = 3) were analysed using the Arraystar Human circRNA Array version 2.0 (Arraystar, Rockville, MD, USA). The RNA sample preparation and microarray hybridization were done according to the manufacturer’s protocol as described previously in Greene et al. [29].

2.5. RT-qPCR and RT-ddPCR

RT-qPCR was performed using SYBR Green Master Mix (Applied Biosystems) on the 7500 Fast Real-Time PCR System (Applied Biosystems) using the following cycling conditions; hold: 50 °C for 20 s (1 cycle), hold: 95 °C for 10 min (1 cycle), denaturation: 95 °C for 15 s and annealing: 60 °C for 30 s (40 cycles). PCR reactions were prepared in a 96-well plate with a final volume of 20 µL per reaction: 10 µL 2X SYBR Green Master Mix, 1 µL forward (FWD) primer (10 µM), 1 µL reverse (REV) primer (10 µM), 7 µL molecular grade water and 1 µL cDNA. Data were normalised to GAPDH and fold change (FC) was calculated using the ∆∆ Ct relative quantification method. RT-ddPCR was performed using EvaGreen Supermix (Bio-Rad, Hercules, CA, USA) on the QX200™ Droplet Digital PCR System (Bio-Rad). Each PCR reaction was first prepared in a 1.5 mL tube (Eppendorf, Hamburg, Germany) with a final volume of 20 µL containing 10 µL 2X EvaGreen Supermix, 1 µL custom FWD primer (10 µM), 1 µL custom REV primer (10 µM) and 6 µL molecular grade water and 2 µL plasma cDNA or 7 µL molecular grade water and 1 µL cell lines cDNA. The 20 µL PCR reaction mix and 70 μL Droplet Generation Oil for EvaGreen (Bio-Rad) was subjected to a droplet generation procedure according to the manufacturer’s instructions. Following completion of droplet generation, 40 μL droplet contents was transferred into ddPCR™ 96-well plates (Bio-Rad), sealed with Pierceable Foil Heat Seal (Bio-Rad) using a PX1™ PCR Plate Sealer (Bio-Rad) and thermal cycled immediately using a T100™ Thermal Cycler (Bio-Rad) under the following conditions; hold: 95 °C for 5 min (1 cycle), denaturation: 95 °C for 30 s and annealing: 60 °C for 60 s (40 cycles) then, a signal stabilization step at 4 °C for 5 min, 90 °C for 5 min and lastly holding at 4 °C. A ramp rate of 2 °C/s was used at every cycling step to avoid droplet instability during rapid temperature alterations and to ensure each individual droplet attained uniform temperature during cycling. The QX200 Droplet Reader (Bio-Rad) was used to count individual droplets based on fluorescence intensity following completion of thermal cycling and analysed by QuantaSoft™ Analysis Pro software (Bio-Rad). Reactions were performed in triplicate and RT-ddPCR results were averaged as absolute copy number (copies/20 μL). The divergent primers overlapping the splice junction of hsa_circ_0001275 and hsa_circ_0001721 were custom designed using Circular RNA Interactome (https://circinteractome.nia.nih.gov/, accessed on 17 April 2018) [35] and synthesised by Integrated DNA Technologies (IDT) (Coralville, Iowa, USA). The convergent primers of GAPDH, Phospholipase C Like 2 (PLCL2), Keratin 1 (KRT1), B-cell leukaemia/lymphoma 11B (BCL11B), Ribosomal L1 Domain Containing 1 (RSL1D1) and Vacuolar Protein Sorting 72 Homolog (VPS72) were custom designed using the in silico PCR software suite in the UCSC Genome Browser (https://www.genome.ucsc.edu/cgi-bin/hgPcr, accessed on 17 April 2018) and synthesised by IDT. Primer sequences are given in Supplementary Table S1.

2.6. Sanger Sequencing

The PCR product was sequenced by Source BioScience (Nottingham, UK). For PCR product preparation, 15 µL post-PCR reaction product was mixed with 6 µL ExoSAP-IT™ reagent (Applied Biosystems), incubated at 37 °C for 15 min to degrade remaining primers and nucleotides, followed by incubation at 80 °C for 15 min to inactivate ExoSAP-IT™ reagent. Following ExoSAP-IT™ treatment, the PCR products were prepared according to Source BioScience purified PCR product specifications.

2.7. Stable Overexpression of hsa_circ_0001275 and hsa_circ_0001721

pCD25-ciR vector containing GFP marker (Geneseed Biotech, Guangzhou, China) (Supplementary Figure S2) was used to overexpress circRNA of interest, hsa_circ_0001275 and hsa_circ_0001721. Primers to amplify both circRNAs were designed (Supplementary Table S2), and PCR products were sub-cloned into pCD25-ciR vector using the EcoRI and BamHI (Roche Diagnostics, Mannheim, Germany) restriction enzymes built into the primers for cloning purposes. LNCaP control cells were transfected with both an empty vector control (EVC) and the respective pCD25-ciR circRNA expressing plasmids using Lipofectamine® 3000 (Life Technologies Corporation, Carlsbad, CA, USA) according to the manufacturer’s protocol to generate the following three cell sub-lines: (i) EVC, (ii) hsa_circ_0001275 (1275 Stable), (iii) hsa_circ_0001721 (1721 Stable). To select for stable integration of the individual plasmids, transfected cells were selected by treatment with geneticin (G418—Invivogen, San Diego, CA, USA) at a concentration of 450 µg/mL (defined as the optimal dose from a G418 killing curve, Supplementary Figure S3) for more than 2 weeks. As the pCD25-ciR plasmid also expresses GFP (Supplementary Figure S2), cells were monitored regularly for enrichment of GFP-expressing cells. At this point expanded stable clones were then sorted using a BD FACSMelody™ (BD Biosciences, San Jose, CA, USA) flow cytometer for GFP expression to create the relevant EVC and circRNA overexpressing clones used in the analyses presented. To establish that the relevant circRNAs had indeed been stably selected for circRNA overexpression was confirmed using RT-qPCR.

2.8. BrdU Proliferation Assay

Cells were seeded at 1 × 10³/well in a final volume of 100 µL/well in a 96-well plate and treated with enzalutamide (Stratech Scientific Ltd., Suffolk, UK) at a range of concentrations (0–80 µM) for a period of seven days with dimethyl sulfoxide (DMSO) serving as the vehicle control. Cell proliferation was measured using a Cell Proliferation ELISA, BrdU (Roche Diagnostics Ltd., Sussex, UK) according to the manufacturer’s protocol. Absorbance was measured on a Vesamax tuneable microplate reader (Molecular Devices, San Jose, CA, USA) at 450 nm with the reference wavelength set to 690 nm. The absorbance value of drug-treated cells was normalized to cells treated with DMSO which was set to 100%.

2.9. Patient Samples

Plasma samples were obtained from 44 patients with MCRPC collected as part of the Phase II trial of enzalutamide in combination with radium-223 (CTRIAL-IE 13-21, NCT02225704) [32]. Docetaxel chemotherapy-naive patients with MCRPC who progressed on ADT were enrolled between June 2015 to July 2017 and the trial closed in August 2021. The treatment regimen consisted of enzalutamide (160 mg per day, oral) in combination with 6 cycles of radium-223 (55 kBq/kg monthly, intravenous) followed by enzalutamide alone until disease progression, unacceptable toxicity or withdrawal of consent. A full description of the patient characteristics in this cohort has previously been described [32]. Blood specimens were collected at baseline before the start of treatment, at month 4 of treatment then every 3 months (±1 month) until disease progression or study withdrawal. Blood samples were collected in K2-EDTA collection tubes and centrifuged at 2000× g, 4 °C for 15 min within 60 min of collection for plasma isolation and stored at −80 °C. The clinical component of this study was nested within a trial registered with Cancer Trials Ireland (CTRIAL-IE 13-21) and with ClinicalTrials.gov (ClinicalTrials.gov identifier: NCT02225704) before the first patient was enrolled. This study received ethical approval by the Clinical Research Ethics Committee of the Cork Teaching Hospitals on 21 July 2014, REF: ECM51 D1/07/14.

2.10. Statistical Analysis

The data were analysed using analysis of variance (ANOVA) and Student’s t-test using GraphPad Software 9.1.0 (San Diego, CA, USA) as appropriate. All results were summarised and presented as mean ± SEM followed by one-way/two-way ANOVA with Tukey’s or Friedman–Šídák multiple comparisons test as appropriate. A value of p < 0.05 was considered statistically significant.

3. Results

3.1. hsa_circ_0001275 and hsa_circ_0000129 Are Up-Regulated in Enzalutamide Resistant Cells

The top 5 up-regulated circRNAs—as ranked by highest FC from our previous enzalutamide resistant cell line panel circRNAs microarray data [29]—along with their associated parental genes were validated via RT-qPCR (Table 1).
Custom-designed divergent primers were used to amplify circRNAs across the backsplice junctions, instead of canonically spliced linear mRNA counterparts (Supplementary Table S1). In order to exclude the possibility of the linear form of RNA detection by the PCR primers, circRNA enrichment using RNase R for hsa_circ_0001275 was performed (Supplementary Figure S1). Additionally, the hsa_circ_0001275 and hsa_circ_0001721 RT-qPCR product was subjected to Sanger sequencing and demonstrated a match to the backsplice junction sequences from the circinteractome database (Supplementary Figure S4a,b).
Out of the five circRNAs validated, two were in line with the microarray data. hsa_circ_0001275 was significantly up-regulated in Clone 1 vs. Control (FC 5.01 ± SEM 0.74, p = 0.006) (Figure 1a). hsa_circ_0001275 is an antisense circRNA, located on chromosome 3 and its associated parental gene PLCL2 was also significantly up-regulated in Clone 1 vs. Control and Clone 9 vs. Control (FC 7.23 ± SEM 0.06, p = 0.002 and FC 8.33 ± SEM 1.19, p = 0.001, respectively) (Figure 1b). The second validated circRNA was an exonic circRNA, located on chromosome 1—hsa_circ_0000129 was up-regulated in Clone 1 vs. Control (FC 2.56 ± SEM 0.21, p = 0.004), however, its associated parental gene VPS72 was not dysregulated between the cell lines (FC 0.79 ± SEM 0.02, p = 0.31) (Figure 1c,d, respectively). The parental gene for hsa_circ_0001721, CDK14, was of interest and selected for validation as our group had previously shown that the circRNA was up-regulated in the highly resistant clone compared to control [29]. CDK14 was also up-regulated in Clone 1 vs. Control (FC 3.51 ± SEM 0.64, p = 0.009) (Figure 1e).
There were no differences in the expression of hsa_circ_0033144 and hsa_circ_0026462 in Clone 1 vs. Control (Supplementary Figure S5a,c), whereas hsa_circ_0000673 was significantly down-regulated in Clone 1 vs. Control (FC 0.31 ± SEM 0.04, p = 0.022) (Supplementary Figure S5e) in contrast to the microarray result. The associated parental genes for hsa_circ_0033144, hsa_circ_0026462 and hsa_circ_0000673 were not differentially expressed in Clone 1 vs. Control (Supplementary Figure S5b,d,f).

3.2. hsa_circ_0001275 Overexpression Contributed to Enzalutamide Resistance

We further studied the effects of hsa_circ_0001275 and hsa_circ_0001721 overexpression on enzalutamide resistance. Both circRNAs were selected on the basis that their parental gene was also up-regulated in the strongly resistant clone compared to control. The transfection efficiency was initially determined by subjective estimation through examination under EVOS™ Digital Colour Fluorescence Microscope (Life Technologies Corporation) (Figure 2a–d) and confirmed quantitatively using RT-qPCR (Figure 3a,b). hsa_circ_0001275 was significantly overexpressed in 1275 Stable vs. Control and vs. EVC (p < 0.05) (Figure 3a). The 1275 Stable cell line demonstrated a significantly increased resistance to enzalutamide at every concentration tested (p < 0.05) except for 80 µM compared to EVC (Figure 4a). The 1275 Stable cell line showed significantly increased resistance to enzalutamide at every concentration tested (p < 0.05) except for 2 µM and 80 µM compared to the Control cell line (Supplementary Figure S6a).
hsa_circ_0001721 was significantly overexpressed in 1721 Stable vs. Control and EVC (p < 0.05) (Figure 3b). The 1721 Stable cell line did not demonstrate enhanced resistance to enzalutamide compared to the EVC cell line (Figure 4b). The 1721 Stable cell line demonstrated significant resistance to enzalutamide at the 2, 5, 15 and 20 µM concentration range (p < 0.05) compared to the Control cell line (Supplementary Figure S6b).

3.3. circRNA Overexpression Did Not Significantly Change the Expression Level of Its Associated Parental Gene

The expression level of the parental gene that encodes its own specific circRNA was investigated in circRNA overexpression cell lines using RT-qPCR. There was no significant difference in the expression of PLCL2 in the cell lines that overexpressed hsa_circ_0001275 (p > 0.05, 1275 Stable vs. EVC and 1275 Stable vs. Control) (Figure 5a), nor CDK14 in 1721 Stable cell lines that overexpressed hsa_circ_0001721 (p > 0.05, 1721 Stable vs. EVC and 1721 Stable vs. Control) (Figure 5b).

3.4. Expression Pattern of hsa_circ_0001275 in Plasma of Patients with Prostate Cancer Treated with Enzalutamide

As our data demonstrated a potential role of hsa_circ_0001275 overexpression with enzalutamide resistance, we further evaluated this circRNA as a biomarker of enzalutamide resistance clinically via RT-ddPCR detection. hsa_circ_0001275 expression was quantified in the plasma of 44 patients with MCRPC on the radium-223 in combination with enzalutamide in metastatic castration-resistant prostate cancer phase II trial (CTRIAL-IE 13-21, NCT02225704) where patients received 6 cycles of radium-223 in combination with enzalutamide followed by enzalutamide alone until disease progression or treatment intolerance. The plasma from four different time points was selected according to the PSA level. This consisted of baseline (prior to treatment, TP1), PSA nadir (absolute lowest PSA level on treatment, TP2), first PSA progression (an increase of ≥25% and absolute increase of ≥2 ng/mL of the PSA level from the nadir, TP3) and second PSA progression (PSA level higher than the first PSA progression, TP4). All 44 patients had plasma available for both TP1 and TP2 time points, whereas only 21 patients had plasma for all four time points since not all patients had PSA progression at the time of this study. Assay optimisation was initially performed using optimized primer concentration and a custom-designed gBlock (Supplementary Table S3) as positive control under a range of temperatures (Supplementary Figure S7a). This step showed that absolute quantification can precisely detect the expression of both circRNAs as shown by the clear separation of positive droplets from the negative droplet background (Supplementary Figure S7b). Furthermore, the expression of hsa_circ_0001275 was also detected in the isogenic LNCaP clones where both circRNAs were significantly up-regulated in the strongly enzalutamide resistant clone compared to control (p < 0.05) further confirming our initial RT-qPCR results (Supplementary Figure S8). To estimate the accuracy between RT-ddPCR and RT-qPCR, we calculated the relative average copies/20 µL of transcript in Clone 1 to Control and Clone 9 to Control and compare the fold change in Clone 1 vs. Control and Clone 9 vs. Control for the two methodologies respectively. The results are presented in Supplementary Table S4. Overall, when compared across cell lines the methodologies have a similar pattern.
The expression of hsa_circ_0001275 in plasma samples showed a trend that mirrors the activity of the disease status according to the PSA level, however, this was not statistically significant (Figure 6a–d). All 44 patients had the first two time points available and their mean average copies/20 µL of hsa_circ_0001275 at TP2 (4 copies/20 µL) (plasma at time of PSA nadir) were lower than baseline TP1 (6 copies/20 µL) (prior to enzalutamide treatment) (Figure 6a). TP2 was the time point where the disease responded to enzalutamide treatment (Figure 6a) as shown by the mean PSA level which was significantly lower at TP2 (75 µg/L) compared to TP1 (10 µg/L) (p < 0.05) (Figure 6b). A subgroup analysis of 21 patients who had four available time points showed a similar trend where the mean average copies/20 µL of hsa_circ_0001275 was high at TP1 (9 copies/20 µL) then reduced at TP2 (6 copies/20 µL) but remained the same at TP3 (6 copies/20 µL) and increased again upon further PSA progression at TP4 (12 copies/20 µL) compared to TP3 (Figure 6c). The mean PSA at nadir, TP2 (20 µg/L) was significantly lower compared to TP1 (94 µg/L) and TP4 (186 µg/L) (p < 0.05), however, this was not significantly different compared to TP3 (90 µg/L) (Figure 6d). This indicated that the disease at TP3 is biochemically less resistant to enzalutamide compared to TP4.

4. Discussion

circRNAs were first discovered in the early 1990s [20], however, they were initially believed to represent by-products of splicing error until 2012 when a seminal paper described the identification of circRNAs in RNA-seq data generated from samples of paediatric acute lymphoblastic leukaemia suggesting that they may be more abundant and of more significant importance than previously perceived [41]. Recently, we and others showed that circRNAs may play important roles in numerous malignancies including prostate cancer [28,29]. However, the molecular function of circRNAs in drug-resistant prostate cancer and their potential as biomarkers to inform treatment strategies remain elusive. A previous study by Wu et al. demonstrated that hsa_circ_0001427 downregulation contributed to enzalutamide resistance, while Cao et al. showed that a circRNA derived from the AR was overexpressed when prostate cancer progressed to CRPC and could be detected in plasma [28,42]. Furthermore, our group identified a number of circRNAs that were differentially expressed in a model of enzalutamide resistance further supporting the notion that circRNAs play a role in enzalutamide resistance and potentially have utility as biomarkers to monitor for the emergence of enzalutamide resistance and prostate cancer progression [29].
From our initial array results, we examined the top 5 circRNAs ranked by highest FC in the strongly resistant cell line (p < 0.05, Clone 1 vs. Control) from our published circRNA microarray data [29] and validated them along with their associated parental gene via RT-qPCR. We found that hsa_circ_0001275 and hsa_circ_0000129 were up-regulated in the strongly resistant cell lines compared to Control (p < 0.05) validating our previous microarray results. Moreover, when examined, the associated parental gene PLCL2 from which hsa_circ_0001275 is derived, was also significantly up-regulated (p < 0.05, Clone 1 vs. Control) (Figure 1b). Interestingly, a previous study identified PLCL2 as part of a 23-gene expression panel that predicts metastatic lethal prostate cancer outcomes in patients with localised disease treated surgically [36] thus, further strengthening the possible association of this circRNA with prostate cancer. In contrast, VPS72 from which hsa_circ_0000129 is derived showed no significant difference in expression (Figure 1c).
Numerous studies have demonstrated that circRNAs are differentially expressed in various human diseases [29,41,43]. Despite this, there is a paucity of literature that has investigated circRNAs’ molecular function [44] and even fewer have investigated the role of circRNA in drug resistance [28]. In order to study the effects of circRNA on enzalutamide resistance, circRNAs of interest were overexpressed in the Control cell line that was sensitive to enzalutamide. hsa_circ_0001275 and hsa_circ_0001721 were selected for overexpression experiments instead of hsa_circ_0000129 since both associated parental genes were significantly up-regulated (p < 0.05, Clone 1 vs. Control) and known to be associated with invasion and metastasis (a known hallmark in cancer) [36,45], and hsa_circ_0001721 was previously validated by our group as being upregulated and associated with resistance to enzalutamide [29]. When stably overexpressed, hsa_circ_0001275 was shown to be associated with increased resistance to enzalutamide (p < 0.05) as determined using proliferation assays (Figure 4a). In contrast, the same proliferative effect was not demonstrated in the hsa_circ_0001721 overexpressing stable cell line suggesting that increased hsa_circ_0001721 does not contribute to enzalutamide resistance (Figure 4b). This indicates only certain up-regulated circRNAs may play a role in enzalutamide resistance by acting as a key regulator through competitive endogenous RNA (ceRNA) networks. Studies have shown a historical role of miRNA “sponging” in circRNA-miRNA–mRNA networks as playing important roles in both gene regulation and carcinogenesis [25]. For example, Xu et al. recently showed that dexamethasone-induced osteoblast growth inhibition can be reversed by downregulating hsa_circ_0001275 which mediated the miR-377/CDKN1B axis in an effort to discover new treatment for osteoporosis [46]. Furthermore, Formosa et al. previously demonstrated a role for miR-377 as a tumour suppressor miRNA targeting FZD4, a gene involved in epithelial-to-mesenchymal transition in prostate cancer [47]. The top five most likely miRNA binding sites for hsa_circ_0001275 were predicted using Arraystar’s miRNA prediction software, and four out of the five predicted miRNAs were hypothesized as tumour suppressor miRNAs in the literature (Supplementary Table S5). Thus, hsa_circ_0001275 may contribute to the development of enzalutamide resistance through these ceRNA networks, however, further investigations are required to fully delineate the circRNA-miRNA-mRNA network. The variability of circRNA expression is however, not dependent on its parental gene expression as shown by Vo et al. [48] and vice versa as shown in our data where expression of the circRNA parental gene was not significantly altered when its circRNA was overexpressed (Figure 5a,b). Thus far, studies in the literature have demonstrated that circRNAs are universally downregulated in proliferative cells across different cancer subtypes, which may be attributed to circRNA concentration dilution upon fast cell division as described by Bachmayr-Heyda et al. [26] and further supported by the findings from Gruner et al. [49] showing accumulation of circRNAs in non-proliferative cells such as ageing nervous tissue. Furthermore, Vo et al. [48] further demonstrated that by using a potent kinase inhibitor to slow down cell growth by decreasing cellular proliferation, the global abundance of circRNAs in LNCaP became elevated and this abundance is independent of the parental gene associated with the circRNA. Our data so far, showed that hsa_circ_0001275 was upregulated along with its parental gene PLCL2 in Clone 1. PLCL2 is associated with aggressive prostate cancer [36]. In the hsa_circ_0001275 (circ1275 Stable) overexpressed cell line, which was more resistant to enzalutamide, there was no difference in the expression of PLCL2 compared to control. PLCL2, therefore, does not contribute to the increase in resistance. As such, from the literature mentioned above, we believe that the enzalutamide resistance seen in our data was not due to the increase in malignancy/proliferation in general, however, the precise role of this circRNA on the enzalutamide mechanism of action is under investigation at this point in time.
Interest in circRNAs as candidate biomarkers in oncology has been established by their known dysregulation across different types of malignancy [50] and their ability to be detected in human bodily fluids such as saliva, blood and gastric fluid [51,52,53]. In order to further investigate the value of circRNA as a biomarker to inform drug resistance, we, therefore, investigated the level of hsa_circ_0001275 in plasma samples of clinical trial (CTRIAL-IE 13-21, NCT02225704) patients with prostate cancer on enzalutamide treatment via RT-ddPCR. These results suggested a trend that corresponded to the activity of the disease defined by the PSA level of patients on enzalutamide treatment, however, it did not reach statistical significance likely due to small cohorts and/or copy numbers were close to the limit of detection for RT-ddPCR. Previous studies have reported that circRNAs have tissue specificity [54] and thus may explain the low abundance of hsa_circ_0001275 in plasma especially if it is not actively secreted by tumour tissue. We did not examine whether or not the particular circRNAs studied were secreted from prostate cancer cells. Our data suggests that hsa_circ_0001275 was not expressed abundantly in plasma and currently limits its potential use as a liquid-based predictive biomarker of enzalutamide resistance. It may be that in the future increased sensitivity may be obtained using a panel of circRNAs.
There are some limitations in our study that should be acknowledged. One of the limitations of the current study is that the study used a single isogenic cell line model of enzalutamide resistance developed by Korpal et al. [33]. In the future, we aim to explore the overexpression of this circRNA in additional enzalutamide-sensitive prostate cancer cell lines to further delineate its association with drug resistance. An additional limitation is that the cohort of patients in the study was relatively small. After initial statistical analysis, we found that our circRNA data were positively skewed. Therefore, we performed the analysis using the Wilcoxon signed-rank test for patients with two time points (Figure 6a,b) and the Friedman test for those with four time points (Figure 6c,d). It has been argued that a sample size of more than five in each group is sufficient, however, the expected size of the difference across groups is also important. Our patient sample size was small which made it difficult to ascertain the statistical significance in particular to confidence interval and p-values. Additionally, despite utilising RT-ddPCR, low copy numbers of hsa_circ_0001275 were present in plasma samples, which may have contributed to large differences across groups. It could be argued that a pre-amplification step such as that used by Karousi et al. [55], may have helped in the detection of circRNAs, however, it must be noted that this approach can cause issues such as multiple PCR products for the same circRNA caused by rolling cDNA from the circRNA as exemplified in the study by Chen et al. [56]. Moving forward for such studies, it may be important to first evaluate their expression using an explant or organoid based protocol to determine whether or not individual circRNAs of interest can be detected within the conditioned media, which will pave the way for more robust studies to test for their presence in patient plasma. Second, we performed a time-point analysis in this study according to PSA levels only (Figure 6). PSA level is not the only surrogate maker to inform disease progression, however, increasing PSA is usually the initial indication of tumour regrowth which is then succeeded by progression of disease on imaging and ultimately symptomatic clinical presentations [57]. Thus, the time points selected may not be a solid reflection of disease progression and tumour burden. Third, this study involved combination treatment of six cycles of alpharadium and enzalutamide followed by enzalutamide maintenance until disease progression or treatment intolerance and the cohort was relatively small. Ideally, the results presented in this study should be examined in a larger study with patients on enzalutamide alone. Fourth, we did not study the level of circRNA expression in tissue samples as they were not available since tissue biopsy upon disease progression were not performed routinely as part of the trial.

5. Conclusions

We found hsa_circ_0001275 up-regulated in enzalutamide-resistant cell lines with data suggesting a protective role against enzalutamide in vitro. While our current research suggests that hsa_circ_0001275 plays a role in resistance to enzalutamide, the functional mechanisms underpinning this have yet to be elucidated for example effects on tumourigenicity or migration and these will form the basis of future studies. Our data also showed that hsa_circ_0001275 is not expressed abundantly in plasma. Although the trend of expression mirrors disease activity according to the PSA level suggest its association with enzalutamide resistance, further studies are required. The treatment landscape of prostate cancer has evolved rapidly in the last decade due to the greater availability and accessibility of therapeutic agents. Despite this arsenal of treatments, patients eventually progress due to drug resistance presenting a difficult conundrum. These novel circRNAs have potential in combating drug resistance and it seems certain that many that are involved remain undescribed. We believe this study provides beneficial groundwork for further exploration of circRNAs potential as biomarkers and therapeutic targets for drug-resistant prostate cancer.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/cancers13246383/s1, Figure S1: hsa_circ_0001275 enrichment using RNase R treatment, Table S1: Custom designed Primer sequence for GAPDH, circRNAs and its associated parental gene, Figure S2: circRNA expressing vector pCD25-ciR (Cat no. GS0106), Table S2: Custom designed Primer sequence for cloning of circRNAs, Figure S3: Representative images of LNCaP cells treated with a range of G418 concentration at the end of 7 days using bright field microscopy, Figure S4: Validation of (a) hsa_circ_0001275 and (b) hsa_circ_0001721 by Sanger sequencing, Figure S5: Validation of circRNAs and associated parental gene in the enzalutamide panel, Figure S6: The effect of enzalutamide on the proliferative rate of Control and cell line overexpressing (a) hsa_circ_0001275 (1275 Stable) and (b) hsa_circ_0001721 (1721 Stable), Table S3: Custom design gBlock as positive control for hsa_circ_0001275 and hsa_circ_000721, Figure S7: RT-ddPCR detection of circRNA, Figure S8: Absolute quantification of hsa_circ_0001275 in enzalutamide resistant panel, Table S4: (a) Relative average copies/20 µL of transcript in Clone 1 to Control and Clone 9 to Control using RT-ddPCR. (b) Fold change in Clone 1 vs. Control and Clone 9 vs. Control using RT-qPCR.; Table S5: Top 5 target miRNAs predicted for has_circ_0001275 using Arraystar’s miRNA prediction software based on Targetscan [1] and miRanda [2] with its hypothesized function and associated target gene.

Author Contributions

Conceptualization, M.C.J.L., S.P.F., A.-M.B., S.G.G. and J.G.; methodology, M.C.J.L., S.P.F., A.-M.B. and S.G.G.; validation, M.C.J.L. and S.G.G.; formal analysis, M.C.J.L., C.M. and P.P.; investigation, M.C.J.L.; resources, S.P.F. and R.S.M.; data curation, M.C.J.L.; writing—original draft preparation, M.C.J.L.; writing—review and editing, A.-M.B., S.P.F., J.G., R.S.M., S.G.G., C.M., K.R. and P.P.; visualization, M.C.J.L. and K.R.; supervision, S.P.F., R.S.M., A.-M.B., S.G.G. and O.S.; project administration, S.P.F.; funding acquisition, S.P.F. and R.S.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Prostate Cancer Foundation Young Investigator Award (S.P.F.) (Award number 11065), Movember (S.P.F.) (Award number 12161), Irish Cancer Society (S.P.F.) (Award number 13102) and Astellas.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Clinical Research Ethics Committee of the Cork Teaching Hospitals (protocol code: ECM5 (1) 01/07/14 and date of approval: 21 July 2014).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The microarray data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus [58] and are accessible through GEO Series accession number GSE118959 (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE118959, accessed on 31 October 2018).

Acknowledgments

The authors thank Eamon Breen from Flow Core, Trinity Translational Medicine Institute (Dublin, Ireland) for providing support in the cell sorting experiment. The authors thank Novartis (Dublin, Ireland) for providing the isogenic LNCaP clones. The authors thank Jake Gill from Accuscience at W91 XF79, Ireland for providing technical support on RT-ddPCR in this research project.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer statistics, 2020. CA Cancer J. Clin. 2020, 70, 7–30. [Google Scholar] [CrossRef]
  2. Noone, A.M.; Howlader, N.; Krapcho, M.; Miller, D.; Brest, A.; Yu, M.; Ruhl, J.; Tatalovich, Z.; Mariotto, A.; Lewis, D.R.; et al. SEER Cancer Statistics Review 1975–2015; National Cancer Institute: Bethesda, MD, USA, April 2018. [Google Scholar]
  3. Schulman, C.C.; Altwein, J.E.; Zlotta, A.R. Treatment options after failure of local curative treatments in prostate cancer: A controversial issue. BJU Int. 2000, 86, 1014–1022. [Google Scholar] [CrossRef] [PubMed]
  4. Pound, C.R.; Partin, A.W.; Epstein, J.I.; Walsh, P.C. Prostate-specific antigen after anatomic radical retropubic prostatectomy. Patterns of recurrence and cancer control. Urol. Clin. N. Am. 1997, 24, 395–406. [Google Scholar] [CrossRef]
  5. Moul, J.W. Prostate specific antigen only progression of prostate cancer. J. Urol. 2000, 163, 1632–1642. [Google Scholar] [CrossRef]
  6. Jemal, A.; Murray, T.; Ward, E.; Samuels, A.; Tiwari, R.C.; Ghafoor, A.; Feuer, E.J.; Thun, M.J. Cancer statistics, 2005. CA Cancer J. Clin. 2005, 55, 10–30. [Google Scholar] [CrossRef]
  7. Nelson, W.G.; De Marzo, A.M.; Isaacs, W.B. Prostate cancer. N. Engl. J. Med. 2003, 349, 366–381. [Google Scholar] [CrossRef]
  8. Heidenreich, A.; Bastian, P.J.; Bellmunt, J.; Bolla, M.; Joniau, S.; van der Kwast, T.; Mason, M.; Matveev, V.; Wiegel, T.; Zattoni, F.; et al. EAU guidelines on prostate cancer. Part II: Treatment of advanced, relapsing, and castration-resistant prostate cancer. Eur. Urol. 2014, 65, 467–479. [Google Scholar] [CrossRef]
  9. Scher, H.I.; Fizazi, K.; Saad, F.; Taplin, M.E.; Sternberg, C.N.; Miller, K.; de Wit, R.; Mulders, P.; Chi, K.N.; Shore, N.D.; et al. Increased survival with enzalutamide in prostate cancer after chemotherapy. N. Engl. J. Med. 2012, 367, 1187–1197. [Google Scholar] [CrossRef] [Green Version]
  10. Fizazi, K.; Scher, H.I.; Molina, A.; Logothetis, C.J.; Chi, K.N.; Jones, R.J.; Staffurth, J.N.; North, S.; Vogelzang, N.J.; Saad, F.; et al. Abiraterone acetate for treatment of metastatic castration-resistant prostate cancer: Final overall survival analysis of the COU-AA-301 randomised, double-blind, placebo-controlled phase 3 study. Lancet Oncol. 2012, 13, 983–992. [Google Scholar] [CrossRef]
  11. Parker, C.; Nilsson, S.; Heinrich, D.; Helle, S.I.; O’Sullivan, J.M.; Fossa, S.D.; Chodacki, A.; Wiechno, P.; Logue, J.; Seke, M.; et al. Alpha emitter radium-223 and survival in metastatic prostate cancer. N. Engl. J. Med. 2013, 369, 213–223. [Google Scholar] [CrossRef] [Green Version]
  12. de Bono, J.S.; Oudard, S.; Ozguroglu, M.; Hansen, S.; Machiels, J.P.; Kocak, I.; Gravis, G.; Bodrogi, I.; Mackenzie, M.J.; Shen, L.; et al. Prednisone plus cabazitaxel or mitoxantrone for metastatic castration-resistant prostate cancer progressing after docetaxel treatment: A randomised open-label trial. Lancet 2010, 376, 1147–1154. [Google Scholar] [CrossRef]
  13. Ryan, C.J.; Smith, M.R.; de Bono, J.S.; Molina, A.; Logothetis, C.J.; de Souza, P.; Fizazi, K.; Mainwaring, P.; Piulats, J.M.; Ng, S.; et al. Abiraterone in metastatic prostate cancer without previous chemotherapy. N. Engl. J. Med. 2013, 368, 138–148. [Google Scholar] [CrossRef] [Green Version]
  14. Antonarakis, E.S.; Lu, C.; Wang, H.; Luber, B.; Nakazawa, M.; Roeser, J.C.; Chen, Y.; Mohammad, T.A.; Chen, Y.; Fedor, H.L.; et al. AR-V7 and resistance to enzalutamide and abiraterone in prostate cancer. N. Engl. J. Med. 2014, 371, 1028–1038. [Google Scholar] [CrossRef] [Green Version]
  15. Sweeney, C.J.; Chen, Y.H.; Carducci, M.; Liu, G.; Jarrard, D.F.; Eisenberger, M.; Wong, Y.N.; Hahn, N.; Kohli, M.; Cooney, M.M.; et al. Chemohormonal Therapy in Metastatic Hormone-Sensitive Prostate Cancer. N. Engl. J. Med. 2015, 373, 737–746. [Google Scholar] [CrossRef] [PubMed]
  16. Fizazi, K.; Tran, N.; Fein, L.; Matsubara, N.; Rodriguez-Antolin, A.; Alekseev, B.Y.; Ozguroglu, M.; Ye, D.; Feyerabend, S.; Protheroe, A.; et al. Abiraterone plus Prednisone in Metastatic, Castration-Sensitive Prostate Cancer. N. Engl. J. Med. 2017, 377, 352–360. [Google Scholar] [CrossRef]
  17. Robinson, D.; Van Allen, E.M.; Wu, Y.M.; Schultz, N.; Lonigro, R.J.; Mosquera, J.M.; Montgomery, B.; Taplin, M.E.; Pritchard, C.C.; Attard, G.; et al. Integrative clinical genomics of advanced prostate cancer. Cell 2015, 161, 1215–1228. [Google Scholar] [CrossRef] [Green Version]
  18. Kamps, R.; Brandao, R.D.; Bosch, B.J.; Paulussen, A.D.; Xanthoulea, S.; Blok, M.J.; Romano, A. Next-Generation Sequencing in Oncology: Genetic Diagnosis, Risk Prediction and Cancer Classification. Int. J. Mol. Sci. 2017, 18, 308. [Google Scholar] [CrossRef]
  19. Bolton, E.M.; Tuzova, A.V.; Walsh, A.L.; Lynch, T.; Perry, A.S. Noncoding RNAs in prostate cancer: The long and the short of it. Clin. Cancer Res. 2014, 20, 35–43. [Google Scholar] [CrossRef] [Green Version]
  20. Jeck, W.R.; Sorrentino, J.A.; Wang, K.; Slevin, M.K.; Burd, C.E.; Liu, J.; Marzluff, W.F.; Sharpless, N.E. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA 2013, 19, 141–157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Jeck, W.R.; Sharpless, N.E. Detecting and characterizing circular RNAs. Nat. Biotechnol. 2014, 32, 453–461. [Google Scholar] [CrossRef] [PubMed]
  22. Lasda, E.; Parker, R. Circular RNAs Co-Precipitate with Extracellular Vesicles: A Possible Mechanism for circRNA Clearance. PLoS ONE 2016, 11, e0148407. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. 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]
  24. 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]
  25. Hansen, T.B.; Kjems, J.; Damgaard, C.K. Circular RNA and miR-7 in cancer. Cancer Res. 2013, 73, 5609–5612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Bachmayr-Heyda, A.; Reiner, A.T.; Auer, K.; Sukhbaatar, N.; Aust, S.; Bachleitner-Hofmann, T.; Mesteri, I.; Grunt, T.W.; Zeillinger, R.; Pils, D. Correlation of circular RNA abundance with proliferation—Exemplified with colorectal and ovarian cancer, idiopathic lung fibrosis, and normal human tissues. Sci. Rep. 2015, 5, 8057. [Google Scholar] [CrossRef]
  27. Papatsirou, M.; Artemaki, P.I.; Scorilas, A.; Kontos, C.K. The role of circular RNAs in therapy resistance of patients with solid tumors. Pers. Med. 2020, 17, 469–490. [Google Scholar] [CrossRef]
  28. Wu, G.; Sun, Y.; Xiang, Z.; Wang, K.; Liu, B.; Xiao, G.; Niu, Y.; Wu, D.; Chang, C. Preclinical study using circular RNA 17 and micro RNA 181c-5p to suppress the enzalutamide-resistant prostate cancer progression. Cell Death Dis. 2019, 10, 37. [Google Scholar] [CrossRef]
  29. Greene, J.; Baird, A.-M.; Casey, O.; Brady, L.; Blackshields, G.; Lim, M.; O’Brien, O.; Gray, S.G.; McDermott, R.; Finn, S.P. Circular RNAs are differentially expressed in prostate cancer and are potentially associated with resistance to enzalutamide. Sci. Rep. 2019, 9, 10739. [Google Scholar] [CrossRef] [Green Version]
  30. Monaghan, A.E.; McEwan, I.J. A sting in the tail: The N-terminal domain of the androgen receptor as a drug target. Asian J. Androl. 2016, 18, 687–694. [Google Scholar] [CrossRef]
  31. Beer, T.M.; Armstrong, A.J.; Rathkopf, D.E.; Loriot, Y.; Sternberg, C.N.; Higano, C.S.; Iversen, P.; Bhattacharya, S.; Carles, J.; Chowdhury, S.; et al. Enzalutamide in metastatic prostate cancer before chemotherapy. N. Engl. J. Med. 2014, 371, 424–433. [Google Scholar] [CrossRef] [Green Version]
  32. McDermott, R.S.; Greene, J.; McCaffrey, J.; Parker, I.; Helanova, S.; Baird, A.M.; Teiserskiene, A.; Lim, M.; Matthews, H.; Deignan, O.; et al. Radium-223 in combination with enzalutamide in metastatic castration-resistant prostate cancer: A multi-centre, phase II open-label study. Ther. Adv. Med. Oncol. 2021, 13, 17588359211042691. [Google Scholar] [CrossRef]
  33. Korpal, M.; Korn, J.M.; Gao, X.; Rakiec, D.P.; Ruddy, D.A.; Doshi, S.; Yuan, J.; Kovats, S.G.; Kim, S.; Cooke, V.G.; et al. An F876L mutation in androgen receptor confers genetic and phenotypic resistance to MDV3100 (enzalutamide). Cancer Discov. 2013, 3, 1030–1043. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Panda, A.C.; Gorospe, M. Detection and Analysis of Circular RNAs by RT-PCR. Bio Protoc. 2018, 8, e2775. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Dudekula, D.B.; Panda, A.C.; Grammatikakis, I.; De, S.; Abdelmohsen, K.; Gorospe, M. CircInteractome: A web tool for exploring circular RNAs and their interacting proteins and microRNAs. RNA Biol. 2016, 13, 34–42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Rubicz, R.; Zhao, S.; Wright, J.L.; Coleman, I.; Grasso, C.; Geybels, M.S.; Leonardson, A.; Kolb, S.; April, C.; Bibikova, M.; et al. Gene expression panel predicts metastatic-lethal prostate cancer outcomes in men diagnosed with clinically localized prostate cancer. Mol. Oncol. 2017, 11, 140–150. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Soudy, R.; Etayash, H.; Bahadorani, K.; Lavasanifar, A.; Kaur, K. Breast Cancer Targeting Peptide Binds Keratin 1: A New Molecular Marker for Targeted Drug Delivery to Breast Cancer. Mol. Pharm. 2017, 14, 593–604. [Google Scholar] [CrossRef]
  38. Mahapatra, S.; Klee, E.W.; Young, C.Y.; Sun, Z.; Jimenez, R.E.; Klee, G.G.; Tindall, D.J.; Donkena, K.V. Global methylation profiling for risk prediction of prostate cancer. Clin. Cancer Res. 2012, 18, 2882–2895. [Google Scholar] [CrossRef] [Green Version]
  39. Li, X.P.; Jiao, J.U.; Lu, L.I.; Zou, Q.; Zhu, S.; Zhang, Y. Overexpression of ribosomal L1 domain containing 1 is associated with an aggressive phenotype and a poor prognosis in patients with prostate cancer. Oncol. Lett. 2016, 11, 2839–2844. [Google Scholar] [CrossRef] [Green Version]
  40. Umemoto, T.; Yamato, M.; Ishihara, J.; Shiratsuchi, Y.; Utsumi, M.; Morita, Y.; Tsukui, H.; Terasawa, M.; Shibata, T.; Nishida, K.; et al. Integrin-alphavbeta3 regulates thrombopoietin-mediated maintenance of hematopoietic stem cells. Blood 2012, 119, 83–94. [Google Scholar] [CrossRef] [Green Version]
  41. Salzman, J.; Gawad, C.; Wang, P.L.; Lacayo, N.; Brown, P.O. Circular RNAs are the predominant transcript isoform from hundreds of human genes in diverse cell types. PLoS ONE 2012, 7, e30733. [Google Scholar] [CrossRef] [Green Version]
  42. Cao, S.; Ma, T.; Ungerleider, N.; Roberts, C.; Kobelski, M.; Jin, L.; Concha, M.; Wang, X.; Baddoo, M.; Nguyen, H.M.; et al. Circular RNAs add diversity to androgen receptor isoform repertoire in castration-resistant prostate cancer. Oncogene 2019, 38, 7060–7072. [Google Scholar] [CrossRef] [PubMed]
  43. Memczak, S.; Jens, M.; Elefsinioti, A.; Torti, F.; Krueger, J.; Rybak, A.; Maier, L.; Mackowiak, S.D.; Gregersen, L.H.; Munschauer, M.; et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 2013, 495, 333–338. [Google Scholar] [CrossRef] [PubMed]
  44. 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]
  45. Imawari, Y.; Mimoto, R.; Hirooka, S.; Morikawa, T.; Takeyama, H.; Yoshida, K. Downregulation of dual-specificity tyrosine-regulated kinase 2 promotes tumor cell proliferation and invasion by enhancing cyclin-dependent kinase 14 expression in breast cancer. Cancer Sci. 2018, 109, 363–372. [Google Scholar] [CrossRef] [Green Version]
  46. Xu, Y.; Sun, L.; Hu, J.; Xu, S. Knockdown of hsa_circ_0001275 reverses dexamethasone-induced osteoblast growth inhibition via mediation of miR-377/CDKN1B axis. PLoS ONE 2021, 16, e0252126. [Google Scholar] [CrossRef] [PubMed]
  47. Formosa, A.; Markert, E.K.; Lena, A.M.; Italiano, D.; Finazzi-Agro, E.; Levine, A.J.; Bernardini, S.; Garabadgiu, A.V.; Melino, G.; Candi, E. MicroRNAs, miR-154, miR-299-5p, miR-376a, miR-376c, miR-377, miR-381, miR-487b, miR-485-3p, miR-495 and miR-654-3p, mapped to the 14q32.31 locus, regulate proliferation, apoptosis, migration and invasion in metastatic prostate cancer cells. Oncogene 2014, 33, 5173–5182. [Google Scholar] [CrossRef] [PubMed]
  48. 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]
  49. Gruner, H.; Cortés-López, M.; Cooper, D.A.; Bauer, M.; Miura, P. CircRNA accumulation in the aging mouse brain. Sci. Rep. 2016, 6, 38907. [Google Scholar] [CrossRef]
  50. Qu, S.; Yang, X.; Li, X.; Wang, J.; Gao, Y.; Shang, R.; Sun, W.; Dou, K.; Li, H. Circular RNA: A new star of noncoding RNAs. Cancer Lett. 2015, 365, 141–148. [Google Scholar] [CrossRef]
  51. Bahn, J.H.; Zhang, Q.; Li, F.; Chan, T.M.; Lin, X.; Kim, Y.; Wong, D.T.; Xiao, X. The landscape of microRNA, Piwi-interacting RNA, and circular RNA in human saliva. Clin. Chem. 2015, 61, 221–230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Memczak, S.; Papavasileiou, P.; Peters, O.; Rajewsky, N. Identification and Characterization of Circular RNAs As a New Class of Putative Biomarkers in Human Blood. PLoS ONE 2015, 10, e0141214. [Google Scholar] [CrossRef]
  53. Shao, Y.; Li, J.; Lu, R.; Li, T.; Yang, Y.; Xiao, B.; Guo, J. Global circular RNA expression profile of human gastric cancer and its clinical significance. Cancer Med. 2017, 6, 1173–1180. [Google Scholar] [CrossRef] [PubMed]
  54. Xia, S.; Feng, J.; Lei, L.; Hu, J.; Xia, L.; Wang, J.; Xiang, Y.; Liu, L.; Zhong, S.; Han, L.; et al. Comprehensive characterization of tissue-specific circular RNAs in the human and mouse genomes. Brief. Bioinform. 2017, 18, 984–992. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Karousi, P.; Artemaki, P.I.; Sotiropoulou, C.D.; Christodoulou, S.; Scorilas, A.; Kontos, C.K. Identification of Two Novel Circular RNAs Deriving from BCL2L12 and Investigation of Their Potential Value as a Molecular Signature in Colorectal Cancer. Int. J. Mol. Sci. 2020, 21, 8867. [Google Scholar] [CrossRef] [PubMed]
  56. Chen, D.-F.; Zhang, L.-J.; Tan, K.; Jing, Q. Application of droplet digital PCR in quantitative detection of the cell-free circulating circRNAs. Biotechnol. Biotechnol. Equip. 2018, 32, 116–123. [Google Scholar] [CrossRef] [Green Version]
  57. Halabi, S.; Vogelzang, N.J.; Ou, S.S.; Owzar, K.; Archer, L.; Small, E.J. Progression-free survival as a predictor of overall survival in men with castrate-resistant prostate cancer. J. Clin. Oncol. 2009, 27, 2766–2771. [Google Scholar] [CrossRef] [Green Version]
  58. Edgar, R.; Domrachev, M.; Lash, A.E. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30, 207–210. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Validation of circRNAs and associated parental genes in the enzalutamide panel (a) hsa_circ_0001275, (b) PLCL2 parental gene, (c) hsa_circ_0000129, (d) VPS72 parental gene, (e) CDK14 parental gene for hsa_circ_0001721 (the circRNA expression has been previously presented [29]). Data graphed as the mean ± SEM (n = 3). Data were analysed using an ordinary one-way ANOVA followed by a Tukey’s post hoc test. (** p ≤ 0.01).
Figure 1. Validation of circRNAs and associated parental genes in the enzalutamide panel (a) hsa_circ_0001275, (b) PLCL2 parental gene, (c) hsa_circ_0000129, (d) VPS72 parental gene, (e) CDK14 parental gene for hsa_circ_0001721 (the circRNA expression has been previously presented [29]). Data graphed as the mean ± SEM (n = 3). Data were analysed using an ordinary one-way ANOVA followed by a Tukey’s post hoc test. (** p ≤ 0.01).
Cancers 13 06383 g001aCancers 13 06383 g001b
Figure 2. Representative images of Control cell line transfected with hsa_circ_0001275 using (a) bright field and (b) fluorescence microscopy and hsa_circ_0001721 using (c) bright field and (d) fluorescence microscopy. All scale bars are 400 µm. All images are 10× magnification.
Figure 2. Representative images of Control cell line transfected with hsa_circ_0001275 using (a) bright field and (b) fluorescence microscopy and hsa_circ_0001721 using (c) bright field and (d) fluorescence microscopy. All scale bars are 400 µm. All images are 10× magnification.
Cancers 13 06383 g002aCancers 13 06383 g002b
Figure 3. Relative expression of (a) hsa_circ_0001275 and (b) hsa_circ_0001721 in the stably transfected cell line compared to Control and EVC cell lines. Data graphed as the mean ± SEM (n = 3). Data were analysed using an ordinary one-way ANOVA followed by a Tukey’s post hoc test. (** p ≤ 0.01, *** p ≤ 0.001).
Figure 3. Relative expression of (a) hsa_circ_0001275 and (b) hsa_circ_0001721 in the stably transfected cell line compared to Control and EVC cell lines. Data graphed as the mean ± SEM (n = 3). Data were analysed using an ordinary one-way ANOVA followed by a Tukey’s post hoc test. (** p ≤ 0.01, *** p ≤ 0.001).
Cancers 13 06383 g003
Figure 4. The effect of enzalutamide on the proliferative rate of EVC and cell line overexpressing (a) hsa_circ_0001275 (1275 Stable) and (b) hsa_circ_0001721 (1721 Stable). Proliferation was measured using BrdU ELISA. Data graphed as the mean ± SEM (n = 3). Analysis was performed using a two-way ANOVA followed by a Šídák’s post hoc test. (* p ≤ 0.05, ** p ≤ 0.01, *** p ≤ 0.001). Dimethyl sulfoxide (DMSO).
Figure 4. The effect of enzalutamide on the proliferative rate of EVC and cell line overexpressing (a) hsa_circ_0001275 (1275 Stable) and (b) hsa_circ_0001721 (1721 Stable). Proliferation was measured using BrdU ELISA. Data graphed as the mean ± SEM (n = 3). Analysis was performed using a two-way ANOVA followed by a Šídák’s post hoc test. (* p ≤ 0.05, ** p ≤ 0.01, *** p ≤ 0.001). Dimethyl sulfoxide (DMSO).
Cancers 13 06383 g004
Figure 5. Relative expression of (a) PLCL2 in the stably transfected cell line which overexpressed hsa_circ_0001275 (1275 Stable) and (b) CDK14 in the stably transfected cell line which overexpressed hsa_circ_0001721 (1721 Stable). Data graphed as the mean ± SEM (n = 3). Data were analysed using an ordinary one-way ANOVA followed by a Tukey’s post hoc test.
Figure 5. Relative expression of (a) PLCL2 in the stably transfected cell line which overexpressed hsa_circ_0001275 (1275 Stable) and (b) CDK14 in the stably transfected cell line which overexpressed hsa_circ_0001721 (1721 Stable). Data graphed as the mean ± SEM (n = 3). Data were analysed using an ordinary one-way ANOVA followed by a Tukey’s post hoc test.
Cancers 13 06383 g005
Figure 6. Scatter plot of hsa_circ_0001275 level in plasma samples and PSA level. (a) hsa_circ_0001275 level in plasma samples and (b) PSA level of all 44 patients with first 2 available time points. Data graphed as mean ± SD (n = 44) and analysed using the Wilcoxon signed-rank test. (c) hsa_circ_0001275 level in plasma samples and (d) PSA level of 21 patients with all 4 time points available. Data graphed as the mean ± SD (n = 21) and analysed using one-way ANOVA with Friedman post hoc test. (* p ≤ 0.05, **** p ≤ 0.0001). Plasma at TP1 was selected as baseline prior to enzalutamide treatment, plasma at TP2 was selected during PSA nadir, plasma at TP3 was selected when PSA first progressed and plasma at TP4 was selected on second PSA progression.
Figure 6. Scatter plot of hsa_circ_0001275 level in plasma samples and PSA level. (a) hsa_circ_0001275 level in plasma samples and (b) PSA level of all 44 patients with first 2 available time points. Data graphed as mean ± SD (n = 44) and analysed using the Wilcoxon signed-rank test. (c) hsa_circ_0001275 level in plasma samples and (d) PSA level of 21 patients with all 4 time points available. Data graphed as the mean ± SD (n = 21) and analysed using one-way ANOVA with Friedman post hoc test. (* p ≤ 0.05, **** p ≤ 0.0001). Plasma at TP1 was selected as baseline prior to enzalutamide treatment, plasma at TP2 was selected during PSA nadir, plasma at TP3 was selected when PSA first progressed and plasma at TP4 was selected on second PSA progression.
Cancers 13 06383 g006aCancers 13 06383 g006b
Table 1. Top five up-regulated circRNAs in Clone 1 vs. Control based on FC with their associated parental gene (with hypothesised gene function) adapted from Greene et al. [29].
Table 1. Top five up-regulated circRNAs in Clone 1 vs. Control based on FC with their associated parental gene (with hypothesised gene function) adapted from Greene et al. [29].
circRNAFC/(p-Value)Parental GeneParental Gene Function
hsa_circ_00012755.8 (p = 0.047)PLCL2Part of a 23-gene signature, which predicted metastatic lethal prostate cancer outcomes [36].
hsa_circ_00264625.7 (p = 0.026)KRT1Target receptor overexpressed in breast cancer cells [37].
hsa_circ_00331445.2 (p = 0.012)BCL11BMethylation occurs in prostate cancer [38].
hsa_circ_00006734.2 (p = 0.038)RSL1D1Associated with an aggressive phenotype and a poor prognosis in patients with prostate cancer [39].
hsa_circ_00001293.9 (p = 0.039)VPS72May have a role in regulating haematopoietic stem cell activity [40].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lim, M.C.J.; Baird, A.-M.; Greene, J.; McNevin, C.; Ronan, K.; Podlesniy, P.; Sheils, O.; Gray, S.G.; McDermott, R.S.; Finn, S.P. hsa_circ_0001275 Is One of a Number of circRNAs Dysregulated in Enzalutamide Resistant Prostate Cancer and Confers Enzalutamide Resistance In Vitro. Cancers 2021, 13, 6383. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13246383

AMA Style

Lim MCJ, Baird A-M, Greene J, McNevin C, Ronan K, Podlesniy P, Sheils O, Gray SG, McDermott RS, Finn SP. hsa_circ_0001275 Is One of a Number of circRNAs Dysregulated in Enzalutamide Resistant Prostate Cancer and Confers Enzalutamide Resistance In Vitro. Cancers. 2021; 13(24):6383. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13246383

Chicago/Turabian Style

Lim, Marvin C. J., Anne-Marie Baird, John Greene, Ciara McNevin, Karine Ronan, Petar Podlesniy, Orla Sheils, Steven G. Gray, Ray S. McDermott, and Stephen P. Finn. 2021. "hsa_circ_0001275 Is One of a Number of circRNAs Dysregulated in Enzalutamide Resistant Prostate Cancer and Confers Enzalutamide Resistance In Vitro" Cancers 13, no. 24: 6383. https://0-doi-org.brum.beds.ac.uk/10.3390/cancers13246383

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