Skip to main content

circCYP24A1 promotes Docetaxel resistance in prostate Cancer by Upregulating ALDH1A3

Abstract

Background

Docetaxel (DTX) is the most widely prescribed first-line chemotherapy for advanced prostate cancer (PCa). Unfortunately, DTX resistance invariably emerges, leading to worse prognosis of PCa. Growing evidence has shown that circRNAs had complex spatiotemporal specificity during the tumor development and oncogenesis. This study was designed to investigate the biological functions and possible molecular mechanisms of circRNAs in DTX resistance of PCa.

Methods

circRNAs in established DTX-resistant DU145 cell line were identified by RNA sequencing. Biological function of circCYP24A1 was verified in vitro and in vivo. The potential role of circCYP24A1 in the development of DTX-resistant PCa was investigated via dual-luciferase reporter assays, RIP assays and RNA pull-down assays. Univariate and multivariate logistic regression analyses was used to predict DTX-chemotherapy response based on patients’ clinical and biological information.

Results

CircCYP24A1 was identified to be upregulated in DTX-resistant DU145 cells. Upregulated circCYP24A1 was found to suppress the DTX chemosensitivity in vitro and in vivo. Furthermore, we found that circCYP24A1 promoted DTX resistance in PCa via regulating ALDH1A3 expression by sponging miR-1301-3p and activating PI3K/AKT/mTOR signaling pathway. Statistical analyses elucidated that circCYP24A1 was an independent risk factor to predict DTX response (OR = 0.165; 95% CI: 0.038–0.723; P = 0.017).

Conclusions

This study demonstrated that circCYP24A played an essential role in DTX resistance in PCa, suggesting that circCYP24A1 could be a promising biomarker to predict DTX response and a potential therapeutic target in PCa patients resistant to DTX chemotherapy.

Background

Prostate cancer (PCa) has the greatest morbidity in aging men [1]. Notwithstanding benefitted from early diagnosis and locoregional therapy, 15% patients present with high-risk, local disease will progress to a metastatic stage and remain a pressing need for systemic therapy [2]. Attributed to the development of therapies for advanced PCa, docetaxel (DTX) is recommended as the first-line therapy applied in the administration of metastatic hormone-sensitive PCa (mHSPC) [3], as well as metastatic castration-resistant PCa (mCRPC) in combination with continuous androgen deprivation therapy [4,5,6].

Unfortunately, DTX resistance invariably emerges, leading to worse prognosis of PCa. In clinical practice, tumor volume (High vs. Low volume) on the conventional images was used to predict the benefits of DTX-based therapy for patients with mHSPC. However, this was just based on the subgroup analyses from two prospective trails [5, 7], and the definition of high or low volume of mHSPC remains controversial [8, 9]. Recently, genetic biomarkers have been shown to be promising in selecting mCRPC patients who might get benefits from PARP inhibitor and AKT inhibitor [10,11,12]. However, no genetic biomarkers have not been proven to be useful in predicting DTX response. Thus, development of biomarkers which is able to predict the benefits from DTX-based chemotherapy and novel strategies targeting signaling pathways involved in DTX resistance are urgently warranted in PCa patients.

Accumulating evidence descripts that the non-coding portion of the transcriptome is crucial for homoeostasis, metabolism, and cell fate [13, 14]. Among them, the expression of circular RNAs (circRNAs) shows complex spatiotemporal specificity during the development and oncogenesis of eukaryotes [15,16,17,18]. The closed-loop structure confers much more stability to circRNAs compared with linear RNAs [19]. Increasing number of researches have indicated the role of circRNAs in a variety of neoplastic processes, such as tumor development, tumor recurrence and chemotherapy resistance [20,21,22]. CircRNAs affect the oncogenesis by serving as miRNA sponges, which could mediate the progression of cancer via uregulating the downstream mRNA. For instance, the first characterized circRNA, ciRS-7, has been well verified to sponge miR-7 and inhibit the expression of miR-7 in diverse cancers [23].

In this work, we found that circCYP241, known as has_circ_0060927, was vital for generation and maintenance of DTX-resistant PCa. Our data illustrated that circCYP241 was highly expressed in DTX-resistant PCa cells. Down-regulation of circCYP241 notably increased DTX-induced cell death. Molecular mechanistic studies determined that circCYP241 bound with miR-1301-3p in the cytosol, and regulated ALDH1A3 expression through post-transcriptional manner. Clinical analysis elucidated that circCYP24A1 could be used to predict DTX response.

Methods

Cell culture and resistant cell line development

HEK293T, DU145 and 22Rv1 cell lines were purchased from the Cell Bank of Type Culture Collection, Chinese Academy of Science (Shanghai, China). All cell lines were determined for mycoplasma using MycAway™ Plus-Color One-Step Mycoplasma Detection Kit (YEASEN, Shanghai, China). The cells were incubated in RPMI-1640 medium (WISENT, Montreal, Canada) supplemented with 10% FBS (Gibco, Grand Island, NY) and 1% penicillin-streptomycin (WISENT), and maintained at 37 °C in a humidified 5% CO2 atmosphere.

The DTX-resistant DU145 (DU145-DR) sub-lines were developed by stepwise increased concentrations of DTX, essentially as described in the following. Parental DU145 cells were treated with DTX (MCE, Monmouth Junction, NJ) suspended in DMSO at an initial concentration of 2 nM for 2 days. The surviving cells were collected and seeded into new dish after one-cycle treatment. As DU145 cells shown resistance, the concentration of DTX was subsequently increased to 100 nM. Before increasing drug concentration or any experiment, cells were allowed to completely recover. Decreased death and increased proliferation indicated the adaptation of cells to DTX. The DTX-resistant sub-cell lines were continuously maintained in 10 nM DTX, and media containing DTX was freshened if necessary. The parental cells were cultured and passaged along with the construction of DU145-DR. Batches of cells were cryopreserved, and assays were conducted on same passages [24].

Tissue specimens and patients’ chemotherapy response assessment

PCa tissue specimens, including resistant and sensitive to DTX, were obtained from 70 patients receiving 6-cycles neoadjuvant treatment (goserelin plus DTX and prednisone) in an institutional board-approved research protocol (clinicaltrials.gov identification ID: NCT04869371) at Nanjing Drum Tower Hospital from Jan 2019 to Dec 2020. All patients received radical prostatectomy after neoadjuvant therapy. All patients provided written informed consent before treatment. The corresponding clinical, radiological and pathological characteristics were collected as well.

The whole-mount pathologic results were recognized as golden standard for assessing neoadjuvant treatment response. Specifically, the largest cross-section bidimensional in postoperative pathologic results shorter than 5 mm were defined as pathologic complete response (pCR) or minimal residual disease (MDR) [25]. In this study, DTX-sensitive was defined as patients with pCR or MDR and the other were regarded as DTX-resistant. This study was undertaken with the understanding and written consent of each subject, with approval of the Ethics Committee of Drum Tower Hospital, Medical School of Nanjing University (Nanjing, China), and in compliance with the principles of the 1964 Declaration of Helsinki.

RNA sequencing

Total RNA was purified from three generations of DU145-DR cells by TRIzol Reagent (Vazyme Biotech Co.,Ltd., Nanjing, China), and the RNA from aging-matched parental cells was set as control. Ribosomal RNA was diminished by RiboMinus Eukaryote kit (Qiagen, Valencia, CA) before establishing the library. Then, the Illumina HiSeq 2000 instrument (Illumina, San Diego, CA) was used to deep sequence the information from the RNA-seq library. The data of RNA-seq were first mapped to the hg38 human reference genome using TopHat2. After dropping out the sequences, which are along to the genome contiguously, the remaining reads were used to characterize circRNAs.

RNA isolation and quantitative realtime PCR (qRT-PCR) assays

TRIzol reagent (Vazyme) was utilized to purify RNA from cell lines following the manufacturers’ protocol. Cytoplasmic and nuclear RNA was extracted applied a Cytoplasmic & Nuclear RNA Purification Kit (Norgen Biotek, Thorold, ON, Canada) following the product description. PrimeScript® RT Reagent Kit with gDNA eraser (Takara, Kusatsu, Japan) were applied for producing cDNA from total RNA. qRT-PCR assays were conducted using SYBR® Premix ExTaq (TaKaRa), and data was collected by the StepOnePlus™ Real-Time PCR System (Applied Biosystems, CA). Otherwise, the expression of miRNAs was quantified by miRNA Universal SYBR qPCR Master Mix (Vazyme) and miRNA 1st Strand cDNA Synthesis Kit (Vazyme) was applied for the single-stranded cDNA synthesis. U6 snRNA was chosen as a loading control to normalize the expression of miRNAs and ACTB mRNA for mRNAs and circRNAs. All the primers for RNAs are listed in Additional file 8: Table S1.

RNA stability assays

Total RNA was treated with 3 U/μg RNase R (Epicentre Technologies, Madison, WI, USA) for 15 min at 37 °C, and RNA was purified for cDNA synthesis. Cells were seeded into dishes and added 5 μg/ml actinomycin D (ActD, MCE) to block transcription for the corresponding time period. The level of circCYP24A1 and CYP24A1 mRNA was analyzed through qRT-PCR.

RNA immunoprecipitation (RIP) assays

The RIP assays were performed with Magna RIP RNA-binding protein immunoprecipitation kit (Bersinbio) following the manufacturer’s protocol. The indicated cells (2X107) were lysed in complete RIP lysis buffer and incubated with human anti-Argonaute2 (AGO2, CST, USA) or non-specific anti-immunoglobulin G (IgG, Millipore) antibody-conjugated beads in rotation for 12 hours at 4 °C. The enriched RNA was obtained from the binding complexes by completely washed, eluted and purified. The enrichments of the target RNAs were detected by qRT-PCR.

RNA pull-down assays

The biotinylated RNA probes (circCYP24A1 and miR-1301-3p) were purchased from GenePharma Technology (Shanghai, China). RNA pull-down assay was conducted with an RNA Antisense Purification (RAP) Kit (BersinBio, Guangzhou, China) following the product manual. The RNA enrichment levels were analyzed by qRT-PCR and normalized to input group. The probes are listed in Additional file 9: Table S2.

Dual-luciferase reporter assay

The potential miR-1301-3p binding site of circCYP24A1 and its mutant were synthesized and subcloned into pGLO-miR plasmids respectively, named as pGLO-circCYP24A1-wild-type (WT) and pGLO-circCYP24A1 mutant (MUT). As the same, the 3′-UTR of ALDH1A3 containing the predicted miR-1301-3p-binding sites or mutant sites was amplified and cloned into pGLO-miR, named as pGLO-ALDH1A3-WT and pGLO-ALDH1A3-MUT. HEK293T cells were co-transfected with indicated luciferase plasmids and mimics. Firefly luciferase activity and Renilla luciferase activity were determined after transfected 48 h with a Dual-Luciferase reporter system (Vazyme), and Renilla luciferase activity was set as internal control.

Protein extraction and western blot

Cells were washed with pre-cooled wash buffer for 3 times and lysed in RIPA buffer (abcam, USA) containing protease inhibitor cocktail (Thermo Fisher Scientific, Waltham, MA) and phosphatase inhibitor cocktail (Thermo Fisher Scientific). Protein content was detected by bicinchoninic acid (BCA) protein assay kit (Vazyme) after centrifugation. The soluble fractions were mixed with 5 × loading buffer and boiled for 5 min. Protein samples were separated by SDS-PAGE and transferred to PVDF membranes (Bio-Rad Laboratories, Hercules, CA, USA). After blocked with 5% (w/v) non-fat milk, the membranes were incubated with appropriately diluted primary antibodies overnight at 4 °C. Horseradish peroxidase (HRP)-conjugated secondary antibodies were probed 1 h at 25 °C. The signals were visualized with ECL reagent (Vazyme), and the intensity values of protein bands were normalized to ACTB. The information about antibodies is showed in Additional file 10: Table S3.

Flow cytometry

Underwent different condition, cells were collected and stained with propidium iodide (PI) and APC-conjugated Annexin V as the product manual (MutiSciences, AP107-V01). Similarly, cell cycle could be detected based on the reagents from a PI/RNase Staining Kit (BD Biosciences). The proportion of apoptosis cells and cell cycle were then analyzed using flow cytometry (BD Biosciences).

CCK-8 assay

CCK-8 assay (Vazyme) was applied for assessment of cell survival. Cells (4 × 103) were seeded into 96-well plates and the DTX was added. CCK-8 working solution was added in after treating with DTX for 2 days. After incubating for 90 min at 37 °C, the absorbance at 450 nm was measured using a microplate reader (TECAN Group Ltd., Infinite M200pro, Mannedorf, CH/CHE, Switzerland).

Immunohistochemistry (IHC)

Clinical PCa specimens and xenografts were used for IHC following the product manual. Paraffin-embedded sections were deparaffinized through a standard procedure. HRP-conjugated species-specific antibodies were incubated in the sections after three times PBST washing. The final signals were detected using DAB Substrate Kit (Servicebio, Wuhan, China) according to the manufacturer’s protocol. For clinical PCa tissues, the staining intensity of ALDH1A3 were scored as 0 (low), 1 (moderate), and 2 (high), and the staining range of tumor was scored as 0 (0–25%), 1 (25–50%), 2 (50–75%), and 3 (75–100%) independently by two well-trained pathologists. The final score was the summation of staining intensity score and staining range score. Low expression was set as the summed scores score ≤ 3 while > 3 as high expression. Images were captured by applying a digital microscope camera (LEICA DM 1000, Nussloch, Germany). The details about antibodies are displayed in Additional file 10: Table S3.

TUNEL assay

To determine the apoptosis levels, xenografts slices of each mouse were prepared and stained using a TUNEL assay kit (Roche Applied Science, Mannheim, Germany) as the instructions. Images were captured by FV3000 confocal fluorescence microscope.

Fluorescence in situ hybridization (FISH)

Cy3-labeled circCYP24A1 (red) probe and FAM-labeled miR-1301-3p (green) probes were purchased from GenePharma. FISH assays were conducted with a FISH Kit (Bersinbio) following the operation manual. DAPI were applied to indicate nuclei. Images were captured by FV3000 confocal fluorescence microscope. The probes are listed in Additional file 11: Table S4.

Plasmid construction, small interfering RNA (siRNA), lentivirus and cell transfection

For constructing circCYP24A1 plasmids, circCYP24A1 sequence was synthesized and subcloned into pCDH-ciR by GenePharma. CircCYP24A1 siRNAs, ALDH1A3 siRNAs, miR-1301-3p mimics and inhibitor were obtained from GenePharma. SiRNAs and plasmids were transfected with Lipofectamine™ 3000 Reagent (Invitrogen, L3000015). Two circCYP24A1-knockdown lentiviruses, obtained from Genechem (Shanghai, China), were applied to transfect DU145-DR cells. Stable knockdown cells lines were obtained with 6 μg/mL puromycin (YEASEN) selection. All sequences were listed in Additional file 12: Table S5.

Animal study

Six-week-old male BALB/c athymic nude mice (n = 6 per group), purchased from Animal Core Facility of Nanjing Medical University (Nanjing, China), were keep in the specific pathogen-free environment at the Animal Laboratory Center of Nanjing University. Mice were subcutaneously injected with 2 × 106 DU145-DR cells stably transfected with corresponding lentivirus suspended in 100 μL of saline. Two weeks after establishment of mouse model, mice were intraperitoneally administrated with saline or 10 mg/kg DTX every 6 days. Four groups were assigned according to the treatment: group 1#, sh-circCYP24A1-transfected cells + DTX; group 2#, sh-circ CYP24A1-transfected cells + saline; group 3#, sh-NC-transfected cells + DTX; and group 4#, sh-NC-transfected cells + saline. Tumor diameters were measured every 3 days, and tumor volume was calculated as length × width2 × π / 6. Mice were sacrificed when tumor diameters longer than 1.5 cm. The tumor xenografts were harvested after sacrificed and fixed in paraformaldehyde for subsequent experiments. Animal experiments were approved and performed in accordance with the Institutional Animal Care and Use Committee of Drum Tower Hospital, Medical School of Nanjing University.

Statistical analysis

Mann-Whitney U test was used for continuous variables and the χ2 test for categorical variables. Student’s t test and one-way analysis of variance (ANOVA) were performed to calculate the significance of difference. To predict the chemotherapy response, univariate and multivariate logistic regression analyses were conducted for significant parameters. Receiver operating characteristic (ROC) analysis for the internal discrimination validation was done by SPSS software (IBM, USA). Corresponding area under the curve (AUC) with 95% confidence interval, sensitivity, specificity, and cut-off value for differentiation were derived. Linear correlation analyzes were used to assess correlations between circCYP24A1, miR-1301-3p and ALDH1A3 expression. The software SPSS 26.0 was employed for statistical analysis, of which all tests were two-sided with statistical significance set at P < 0.05 (*P < 0.05, **P < 0.01 and ***P < 0,001). All values were shown as the means ± standard deviation of three independent experiments.

Results

CircCYP24A1 elevates in DTX-resistant PCa

In order to explore the underlying mechanism of DTX resistance, DTX-resistant cell lines (DU145-DR) were established. Cells were considered to be DTX-resistant when they were able to remain 90–100% viable detected by CCK-8 assay in the presence of this 10-fold higher than IC50 of DTX (5.7 nM), and the IC50 of DU145-DR was 122.3 nM (Additional file 1: Fig. S1A). Moreover, DU145-DR cells showed resistance to DTX in vivo (Additional file 1: Fig. S1B, C). Additionally, flow cytometry analyses revealed that DU145-DR cells had a lower apoptosis rate (Additional file 1: Fig. S1D, F) and less G2/M phrase arrest (Additional file 1: Fig. S1E, G) than the parental cells after treated with 20 nM DTX. In general, DU145-DR cells revealed a significant resistance to DTX compared with the parental DU145 cells in vivo and in vitro.

RNA sequencing (RNA-seq) analysis was conducted to investigate the differential expression of circRNAs in DTX-resistant DU145 cells. With a qualifying standard of fold-change > 2 (or < − 2) and P value < 0.05, 107 upregulated and 67 downregulated circRNAs were identified (Fig. 1A, B; Additional file 2: Fig. S2A). Among these circRNAs, 94.1% were the products from back-splicing of protein-coding exons (Additional file 2: Fig. S2B). In order to amend the “Type I” RNA-seq error and accurately assess the variation, a heat map was generated to visually demonstrate the differential circRNAs expression with false discovery rate < 0.05 (Fig. 1C, D). Finally, 7 upregulated and 4 downregulated circRNAs were further confirmed in DTX-resistant PCa cells using qRT-PCR (Fig. 1E). Among them, has_circ_0060927 (Chr20:54157168–54,171,670) was found to be the most remarkably upregulated in DU145-DR cells. According to the CircBase database annotation, hsa_circ_0060927, with a length of 1106 bp, is back spliced from exons 3 to 11 of CYP24A1 gene and therefore is named as circCYP24A1. To further confirm the existence of circCYP24A1, the specific convergent and divergent primers were designed respectively. The result revealed that only the amplification from cDNA templates showed the existence of circCYP24A1, while both cDNA and gDNA templates could amplified by convergent primer of the CYP24A1 mRNA (Fig. 1F). Moreover, the back-splicing junction of circCYP24A1 was directly approved by first-generation sequencing analysis (Fig. 1G). Further, ActD and RNase R exonuclease assay were used to confirm the circular form of circCYP24A1, which showed more stable than linear mRNA (Fig. 1H, I). According to the RNA nuclear-cytoplasmic fractionation and FISH assay, circCYP24A1 located predominantly in the cytoplasm (Fig. 1J, K).

Fig. 1
figure 1

CircCYP24A1 was identified and characterized in DTX-resistant PCa. A Volcano plots visualizing 107 upregulated and 67 downregulated circRNAs in DU145-DR cells. B The chromosomal localization of those significantly expressed circRNAs. C Hot map visualizing 7 upregulated and 4 downregulated circRNAs in DU145-DR cells according to the log2 (fold-change) > 1 or < − 1 and q value < 0.05. D Basic information for the dysregulated circRNAs. E qRT-PCR analysis of the dysregulated circRNAs in DU145-DR and DU145 cells. F PCR and agarose gel electrophoresis analysis were used to verify the existence of linear and circular forms of CYP24A1 in cDNA and gDNA samples from DU145-DR cells. G Schematic illustration of circCYP24A1 back splicing junction validated by Sanger sequencing. H ActD treatment was used to evaluate the stability of circCYP24A1 and CYP24A1 mRNA in DU145-DR cells. I qRT-PCR analysis of the circCYP24A1 and CYP24A1 mRNA level after RNase R exonuclease treatment. J Nuclear-cytoplasmic fractionation assay was performed to observe the subcellular distribution of circCYP24A1 in DU145-DR cells. The ACTB mRNA and U6 snRNA were used as cytoplasmic and nuclear markers, respectively. K The circCYP24A1 probe labeled with Cy3 (red) was used to visualize the location of circCYP24A1 in DU145-DR cells via FISH assay, and the nuclei were stained with DAPI (blue). L FISH score of circCYP24A1 in DTX-sensitive (n = 30) and DTX-resistant (n = 40) PCa biopsy tissues. M Representative FISH images of circCYP24A1 in three DTX-resistant/−sensitive patients, respectively. The circCYP24A1 probe was labeled with Cy3 (red), and the nuclei were stained with DAPI (blue). All data are presented as the means ± SD of three independent experiments. n = 3, ***P < 0.001

Next, we verified the expression of circCYP24A1 in biopsy tissues from DTX sensitive (n = 30) and resistant (n = 40) patients. Interestingly, the circCYP24A1 expression was found to be significantly higher in DTX-resistant PCa specimens based on the FISH score (P = 0.002; Fig. 1L, M; Table 1). These findings revealed that circCYP24A1, elevated in DTX-resistant PCa, might regulate DTX sensitivity.

Table 1 Clinical and radiological characteristics of patients with sensitive or resistant to DTX

CircCYP24A1 promotes DTX resistance of PCa in vitro

To identify the potential roles of circCYP24A1 in DTX resistance, DU145-DR cells were transfected with shRNAs by lentivirus to knockdown circCYP24A1 (Fig. 2A, B). In addition, circCYP24A1 plasmid were transfected into DU145 and 22Rv1 cell lines for circCYP24A1 overexpression without significantly affecting the expression of CYP24A1 mRNA (Fig. 2C, D). Silence of circCYP24A1 evidently decreased IC50 value of DTX in DU145-DR cells, while the IC50 concentration of DTX was increased in DU145 and 22Rv1 cells with overexpression of circCYP24A1 (Fig. 2E-G). Furthermore, flow cytometry analysis was performed and the results revealed that knockdown circCYP24A1 caused the increased apoptosis rate (Fig. 2H) and more G2/M phrase arrest (Fig. 2K). Correspondingly, overexpression circCYP24A1 reduced cell apoptosis and G2/M phrase arrest under DTX treatment in both DU145 and 22Rv1 cells (Fig. 2I, J, L, M). These findings indicated that circCYP24A1 might partly promote the resistance of PCa cells to DTX.

Fig. 2
figure 2

CircCYP24A1 enhances DTX resistance in vitro. A Schematic illustration of two sh-circCYP24A1 RNAs targeted the back-splicing junction. B-D qRT-PCR analysis of circCYP24A1 and CYP24A1 mRNA expression in DU145-DR (B), DU145 (C) and 22Rv1 (D) cells transfected with indicated lentivirus. E-G CCK-8 assays were used to detected the DTX cytotoxicity after down−/up-regulation of circCYP24A1 in DU145-DR (E), DU145 (F) or 22Rv1 (G) cells. H-J The indicated cells were treated with DTX (20 nM for DU145-DR/DU145 and 3 nM for 22Rv1) or PBS (control) and subjected to Annexin V-APC and PI stanning to detect apoptotic rate by flow cytometry. K-M Cell cycle were detected by flow cytometry in indicated cells after treatment with DTX (20 nM for DU145-DR/DU145 and 3 nM for 22Rv1). All data are presented as the means ± SD of three independent experiments. n = 3, *P < 0.05, **P < 0.01, ***P < 0.001

CircCYP24A1 suppresses the DTX Chemosensitivity of PCa in vivo

To further verify the bio-function of circCYP24A1 in DTX chemosensitivity, the mouse xenograft model was established by subcutaneously injected DU145-DR cells transfected with shRNAs of circCYP24A1 into immunodeficient athymic mice. Two weeks after establishment, tumor-bearing nude mice were treated with DTX (10 mg/kg) every 6 days (Fig. 3A). The results revealed that silence of circCYP24A1 increased DTX efficacy, indicated by significantly smaller tumor appearance (Fig. 3B, C). In addition, lower volume of tumors was detected in the DTX/sh-circCYP24A1 group, compared with those in the other groups (Fig. 3D, E). Moreover, significantly increased cleaved caspase 3 expression and decreased Ki67 level were observed in DTX/sh-circCYP24A1 group compared with the DTX/sh-NC group (Fig. 3F). Correspondingly, the DTX/sh-circCYP24A1 group presented more TUNEL-positive nuclei, suggesting more apoptosis after DTX treatment (Fig. 3G). In general, our findings indicated that circCYP24A1 could suppress DTX chemosensitivity in vivo.

Fig. 3
figure 3

Knocking down circCYP24A1 suppressed the DTX chemoresistance of PCa in vivo. A Schematic diagram illustrates the xenograft mouse model treated by indicated concentration of DTX. B, C Images of the BALB/c nude mice (B) and the tumor xenografts (C) at the endpoint. D Tumor growth curves of tumor sizes are shown. Tumor volumes were measured every 3 days since the first DTX treatment. E Tumor volumes measured at the experimental endpoint. F Representative images of H&E and IHC (Cleaved caspase 3 and Ki67) staining of the xenograft tumors. G TUNEL staining of the xenograft tumors. All data are presented as the means ± SD of three independent experiments. n = 6, ***P < 0.001

CircCYP24A1 functions as a molecular sponge of miR-1301-3p in PCa cells

Growing evidence has shown that lots of circRNAs modulated the downstream mRNA expression by sponging miRNAs [26]. The intracellular distribution of circCYP24A1 suggested the ability of binding miRNAs. The cross-analysis of three miRNA databases (miRanda [27], RNAhybird [28] and Targetscan [29]) was used to predict the possible miRNAs which could bind to circCYP24A1. Eight candidate miRNAs were identified for the further validation (Fig. 4A). Next, the biotin-labeled circCYP24A1 probe was designed to indirectly capture the miRNAs in DU145-DR and DU145 cells. Our results revealed that circCYP24A1 was mostly enriched by the circCYP24A1 probe compared with the random sequence probe, confirming the specificity of the probe (Fig. 4B). Further, miR-2467-5p/− 1301-3p/−3135b were found to be the significantly enriched miRNAs in both DU145-DR and DU145 cells (Fig. 4C, D). Moreover, the levels of miR-1301-3p/−3135b were upregulated in DU145-DR cells after silencing circCYP24A1 (Fig. 4E). Interestingly, circCYP24A1 overexpression only reduced the level of miR-1301-3p in both DU145 and 22Rv1 cells (Fig. 4F, G). In addition, we confirmed that circCYP24A1 had potential effects only on miR-1301-3p expression among those 8 candidate miRNAs by up−/down-regulation of circCYP24A1. Furthermore, the RNA pull-down assays verified the interaction between biotin-labelled miR-1301-3p and circCYP24A1 (Fig. 4H, I). Dual-luciferase assay was used to confirm that miR-1301-3p could bind to circCYP24A1 directly. The full length of circCYP24A1 (wild-type, WT) was subcloned into pGLO-miR, and the mutation plasmid containing mutated potential binding region of miR-1301-3p was constructed. Luciferase activities were measured after transfected with miR-1301-3p and luciferase reporter plasmid, and a remarkable reduction of luciferase activities could be abolished by mutation (Fig. 4J). As previously reported, AGO2 was characterized as a potent component of RNA-induced silencing complex to participate the function of miRNA [30]. Thus, an anti-AGO2 RIP assay was performed to further testify the interaction between circCYP24A1 and miR-1301-3p (Fig. 4K). The ectopic expression of circCYP24A1 notably increased the enrichment of miR-1301-3p in AGO2 complex instead of IgG (Fig. 4L). Similarly, upregulation of miR-1301-3p dramatically increased enrichment of circCYP24A1, compared with that in the negative control (Fig. 4M). FISH assays visualized that circCYP24A1 and miR-1301-3p colocalized in the cytosol of DU145-DR cells (Fig. 4N). Moreover, miR-1301-3p levels were assessed by FISH score. Higher level of miR-1301-3p was found in specimen from those DTX-sensitive patients, compared with those from DTX-resistant patients (Fig. 4O, P; Table 1). In summary, these results suggested that miR-1301-3p might directly interact with circCYP24A1 in potential binding sites.

Fig. 4
figure 4

CircCYP24A1 serves as a sponge for miR-1301-1p. A Venn diagram showing the overlap of the eight candidate target miRNAs of circCYP24A1 predicted by miRanda, RNAhybird and TargetScan. B The efficiency of the circCYP24A1 probe in DU145-DR and DU145 cells. A random sequence probe served as a negative control. C-D The fold enrichment of the target miRNAs in pull-down complex of DU145-DR cells (C) and DU145 cells (D). E-G The expression of those candidate miRNAs after down−/up-regulation of circCYP24A1 in DU145-DR (E), 22Rv1 (F) and DU145 cells (G). H Schematic view of RNA pull down assays by biotin-labeled miRNA. I Biotin-labeled miRNA pull-down in DU145 and DU145-DR cells, and qRT-PCR results showing circCYP24A1 enrichment in miRNA probe. A random sequence probe served as a negative control. J The luciferase activities of the circCYP24A1 luciferase reporter in HEK293T cells transfected with miR-1301-3p mimics or mimic NC (upper) and the schematic of the wild-type (WT) and mutant (MUT) circCYP24A1 luciferase reporter vectors (bottom). K Schematic view of anti-AGO2 RIP assays. L-M. Anti-AGO2 RIP assays were performed in DU145 cells, and the enrichment of circCYP24A1 (L) and miR-1301-3p (M) was detected by qRT-PCR. N The colocalization of circCYP24A1 and miR-1301-3p in DU145-DR cells was detected by FISH assay. The circCYP24A1 probe was labeled with Cy3 (red), miR-1301-3p probes were labeled with FAM (green), and nuclei were stained with DAPI (blue). O FISH score of miR-1301-3p in DTX-sensitive (n = 30) and DTX-resistant (n = 40) PCa biopsy tissues. P Representative FISH images of miR-1301-3p in three DTX-resistant/−sensitive patients, respectively. The miR-1301-3p probe was labeled with FAM (green), and the nuclei were stained with DAPI (blue). All data are presented as the means ± SD of three independent experiments. n = 3, *P < 0.05, **P < 0.01, ***P < 0.001

CircCYP24A1 enhances the DTX-resistance of PCa cells through miR-1301-3p

To further verify the bio-functions of miR-1301-3p in DTX-resistant, corresponding mimics and inhibitor were used to regulate its expression in DU145-DR cells. The efficiency of miR-1301-3p mimics and inhibitor in cells was examined by qRT-PCR (Fig. 5A, B). The results of CCK-8 assays indicated that the miR-1301-3p inhibitor substantially attenuated DTX cytotoxicity and reversed the reduced IC50 value caused by circCYP24A1 downregulation (Fig. 5C). Similarly, flow cytometry analysis revealed that miR-1301-3p inhibitor could markedly decrease apoptosis rate (Fig. 5D, F) and less G2/M phrase arrest (Fig. 5E, G). Collectively, these data illustrated that miR-1301-3p promoted DTX sensitivity in DU145-DR cells and played a crucial function in the downstream of circCYP24A1.

Fig. 5
figure 5

CircCYP24A1 sponges miR-1301-3p to regulate DTX sensitivity in PCa cells. A The expression of miR-1301-3p in DU145-DR cells transfected with miR-1301-3p mimics. B The expression of miR-1301-3p in DU145-DR cells transfected with shcircCYP24A1 and miR-1301-3p inhibitor. C CCK-8 data for the IC50 of DU145-DR cells treated with 20 nM DTX after transfection with miR-1301-3p inhibitor combined with or without sh-circCYP24A1. D-E The indicated cells were treated with 20 nM DTX and subjected to detect apoptotic rate (D) or cell cycle (E) by flow cytometry. F Statistical analysis of the apoptotic rate. G Statistical analysis of the cell cycle. All data are presented as the means ± SD of three independent experiments. n = 3, *P < 0.05, **P < 0.01, ***P < 0.001

CircCYP24A1 upregulates ALDH1A3 expression by post-transcriptional mechanisms in DTX-resistant PCa cells

Cross-analysis between four bioinformatic database (PITA, miRanda, Targetscan and ENCORI) and mRNA-seq upregulated in DU145-DR cells were performed to explore the function of miR-1301-3p in DTX-resistant PCa cells. As the results, 16 potential target genes of miR-1301-3p were identified for the following investigations (Fig. 6A). ALDH1A3 was found to be the most possible downstream target of miR-1301-3p according to significant regulation of ALDH1A3 by miR-1301-3p (Fig. 6B, D). Thus, we supposed that ALDH1A3 and miR-1301-3p might form a ceRNA network in DTX-resistant PCa.

Fig. 6
figure 6

CircCYP24A1 upregulates ALDH1A3 expression by post-transcriptional mechanisms in DTX-resistant PCa. A Venn diagram visualizing the overlap of 13 genes which was predicted according to the cross-analysis of PITA, miRanda, ENCORI, TargetScan and the mRNA-seq results. B, D The expression of those candidate genes after transfected with miR-1301-3p mimic (B) or inhibitor (D). C The enrichment of ALDH1A3 mRNA was detected by biotinylated miRNA pull-down in DU145-DR and DU145 cells. E Anti-AGO2 RIP assays were performed in DU145 cells, and the enrichment of ALDH1A3 was assessed by qRT-PCR. F The luciferase activities were determined in HEK293T cells co-transfected with ALDH1A3 luciferase reporter vectors and miR-1301-3p mimics or mimic NC (bottom), and the schematic showed wild-type (WT) and mutant (MUT) ALDH1A3 luciferase reporter vectors (upper). G-L The RNA and protein levels of ALDH1A3 was detected by qRT-PCR or western blot in indicated cells, respectively. M IHC score of ALDH1A3 in DTX-sensitive (n = 30) and DTX-resistant (n = 40) PCa biopsy tissues. N Representative IHC images of ALDH1A3 in three DTX-resistant/−sensitive patients, respectively. All data are presented as the means ± SD of three independent experiments. n = 3, **P < 0.01. ***P < 0.001

To test the above hypothesis, biotin-labelled miR-1301-3p based RNA pull-down assays and anti-AGO2 RIP assays were performed to verify the direct interaction between miR-1301-3p and ALDH1A3 mRNA. As shown in Fig. 6C, biotin-labelled miR-1301-3p could directly pull down ALDH1A3 mRNA in both DU145-DR and DU145 cells. Likewise, the overexpression of miR-1301-3p augmented the enrichment of ALDH1A3 mRNA in AGO2 (Fig. 6E). Further, the wild-type and mutant ALDH1A3 3′-UTR were cloned into dual-luciferase reporter plasmids respectively. The luciferase activity was notably downregulated after ectopic expression of miR-1301-3p, which showed no influence on MUT-3′-UTR of ALDH1A3 plasmid (Fig. 6F). Moreover, down−/up- regulation of circCYP24A1 could modulate the mRNA and protein levels of ALDH1A3, which could be mostly rescued while silencing miR-1301-3p (Fig. 6G-L; Additional file 3: Fig. S3). The IHC staining results revealed that ALDH1A3 were upregulated induced DTX treatment and downregulated due to circCYP24A1 knockdown in tumor-bearing mice (Additional file 4: Fig. S4). Similar as circCYP24A1, the IHC score emphasized that the expression of ALDH1A3 were evidently elevated in DTX-resistant PCa patients compare with those sensitive to DTX (Fig. 6M, N). To summarize, the enhanced expression of ALDH1A3 was induced by circCYP24A1 via regulation of miR-1301-3p through post-transcriptional mechanisms and related to DTX resistance.

ALDH1A3 regulates the influence of circCYP24A1 on the DTX chemosensitivity of PCa cells via PI3K/AKT/mTOR axis

In the previous studies, ALDH1A3 has been verified to play a fundamental role in resistance of temozolomide [31]. Nonetheless, there has been none evidence supporting the association of ALDH1A3 and DTX resistance. To further confirm whether ALDH1A3 is involved in DTX resistance, cell cytotoxicity and IC50 assays were performed to verify this hypothesis. We found that knockdown of ALDH1A3, validated by western blot (Additional file 5: Fig. S5), induced decreased IC50 and mitigated DTX resistance in DU145-DR cells (Fig. 7A). Moreover, flow cytometry analysis uncovered that a notably increased apoptosis rate was induced by si-ALDH1A3 under the treatment of DTX in DU145-DR cells (Fig. 7B, C), and more G2/M phrase arrest was observed correspondingly (Fig. 7D, E).

Fig. 7
figure 7

ALDH1A3 promotes DTX resistance via activating PI3K/AKT/mTOR signaling pathway. A The IC50 values of DU145-DR cells were detected after treatment with 20 nM DTX and transfection with si-ALDH1A3. B-E Flow cytometry analysis were used to assess cell apoptosis (B, C) and cell cycle (D, E) in DU145-DR cells treated with 20 nM DTX after transfection of si-ALDH1A3. F-I The protein levels of PI3K and downstream molecules were determined in those cells co-transfected with siRNAs or overexpression plasmid and miR-1301-3p inhibitors by western blot. All data are presented as the means ± SD of three independent experiments. n = 3, **P < 0.01, ***P < 0.001

In accordance, PI3K/AKT signal pathway was reported to function momentously in various tumors to induce the resistance of CDK4/6 inhibitors or temozolomide and be involved in ALDH1A3-dependent tumor progression [32,33,34]. Obviously, our results illustrated that silence of ALDH1A3 induced a notably depressed activity of PI3K/AKT/mTOR pathway in DU145-DR cells (Fig. 7F). Likewise, knockdown of circCYP24A1 inactivated this pathway, which could be rescued by the knockdown of miR-1301-3p (Fig. 7G, H). In addition, we found that overexpressed circCYP24A1 enhanced the activity of PI3K/AKT/mTOR pathway (Fig. 7I). According to these findings, ectopic expression of ALDH1A3 might induce DTX resistance in PCa via PI3K/AKT/mTOR signaling pathway.

The clinical characteristics of the 70 included patients were listed in Additional file 13: Table S6, and the biological characteristics were based on the score of IHC standing (ALDH1A3) and FISH assay (circCYP24A1 and miR-1301-3p). Scatter plots illustrated a close relationship between the expression of circCYP24A1, miR-1301-3p and ALDH1A3, respectively (P < 0.05) (Additional file 6: Fig. S6 A-C). Meanwhile, the expression of circCYP24A1, miR-1301-3p and ALDH1A3 as well as the post-treatment prostate specific antigen (PSA) level were significantly associated with DTX chemotherapy response (P < 0.05 Table 1). In multivariable logistic regression analyses, both the expression of circCYP24A1 (OR = 0.165; 95% CI: 0.038–0.723; P = 0.017) and post-treatment PSA (OR = 0.033; 95% CI: 0.002–0.609; P = 0.022) could serve as independent risk factors of DTX response (Table 2, Additional file 7: Fig. S7,). In summary, these findings revealed that circCYP24A1, miR-1301-3p and ALDH1A3 were significantly associated with DTX resistance, and circCYP24A1 could be used to predict DTX response.

Table 2 Univariate and multivariate logistic regression analyses of clinical and biological parameters for the prediction of DTX response after chemotherapy

Discussion

DTX-based chemotherapy is recommended as the first-line therapy for mCRPC patients [3]. Since the publication of STAMPEDE and CHAARTED trials, DTX-mediated chemotherapy has become the first-line treatment for patients with mHSPC [5]. Even though DTX significantly prolongs the lifespan and improves the clinical outcomes, most patients develop DTX-resistance within 3 years [35]. Several molecular mechanisms regarding DTX resistance has been proposed. CircRNAs, a novel class of ncRNAs, had been suggested to mediate chemotherapy-resistance in multiple cancer. For instance, circRNA-SORE increases sorafenib resistance by stabilizing YBX1 in hepatocellular carcinoma [22]. CircAKT3 enhances resistance of cisplatin via miR-198/PI3KR1 axis in gastric cancer [36]. In present study, we identified that circCYP24A1 was upregulated in DU145-DR cells and DTX-resistant PCa tissues. Further experimental results indicated that circCYP24A1 promoted DTX resistance, suggesting that circCYP24A1 could serve as a potential DTX-resistant biomarker in PCa patients.

Recently, emerging evidence has indicated that the most common biological function of circRNAs was to play as molecular sponges of miRNAs and form ceRNA network to involve post-transcriptional regulation [37]. In this study, we found that circCYP24A1 could provide multiple potential binding sites for eight miRNAs according to the cross-analysis of miRNA databases. Subsequently, the data of RNA pull-down and RIP assays uncovered that miR-1301-3p was the most highly enriched candidate at circCYP24A1. Moreover, the rescue experiment revealed that suppression of DTX resistance induced by circCYP24A1-knockdown could be rescued mostly by miR-1301-3p inhibitor. Based on the above results, we confirmed that circCYP24A1 sponged miR-1301-3p and formed an activated ceRNA regulatory network to induce DTX-resistant in PCa.

It has been well described that, miRNAs exerted their biological functions mostly through post-transcriptional regulation by suppressing the translation or promoting the degradation of target mRNA. Khalid et al. reported miR-138 enhanced DTX resistance via K2/β-integrin axis in mCRPC [38] while Xu et al. showed miR-143 promoted DTX cytotoxicity through suppression of KRAS [39]. Several previous studies have suggested the anti-tumor activity of miR-1301-3p. Qiao et al. verified that miR-1301-3p suppressed tumor development of thyroid papillary cancer through decreasing expression of PCNA [40]. Peng et al. illustrated miR-1301-3p inhibited tumor proliferation through directly mediating ICT1 mRNA in breast cancer [41]. Nevertheless, there has been none evidence to verify the bio-functions of miR-1301-3p in DTX resistance. According to our results, ALDH1A3 and MXD4 were characterized as the potential target mRNA of miR-1301-3p. Meanwhile, we observed that ALDH1A3 mRNA levels had a higher fold-change than MDX4 after transfection of miR-1301-3p mimics or inhibitor, thus we set ALDH1A3 as the target mRNA for further investigation. Here, the results revealed that miR-1301-3p could directly interact with the ALDH1A3 mRNA and downregulated ALDH1A3 in DTX-resistant PCa.

ALDH1A3, a member of aldehyde dehydrogenase family 1, has been proven to convert acetaldehyde to acetate to produce pyruvate [42], nuclear acetyl-CoA [43, 44] and citrate [45]. Interestingly, few studies illustrated the connection between ALDH1A3 expression and DTX resistance. It has been reported that LoVo-1 and K1 cells with high expression of ALDH1A3 were resistant to the focal adhesion kinase autophosphorylation inhibitor Y15, while cancer cells with low ALDH1A3 level showed high sensitivity to Y15 [46]. A study raised by Kong. et al. displayed the pleural malignant mesotheliomas with dysregulation of ALDH1A3 were robustly resistant to cisplatin and pemetrexed [47]. In our study, a sequence of functional experiments data suggested that knockdown of ALDH1A3 decreased the resistance of DTX and increased the DTX-induced apoptosis. Moreover, our findings suggested that the abnormal expression of ALDH1A3 promoted DTX resistance via PI3K/AKT/mTOR axis.

From the multivariable analyses, the expression of circCYP24A1 and post-treatment PSA were found to be independent risk factors of DTX response. Though magnetic resonance imaging (MRI) is the standard imaging for detecting significant PCa, the parameters on MRI did not show good predicting performance compared with biomarkers such circRNA and PSA. As we know, PSA is the most common biomarker for the diagnosis of PCa. Also, post-treatment PSA has been reported to be effective to predict treatment response of neoadjuvant therapy [48, 49]. Zilli, T. et al. illustrated that nadir PSA before radiotherapy was a vitally independent risk factor for the prediction of the treatment response to neoadjuvant ADT [50]. However, nadir PSA could only be revealed after several cycles of treatment, which was not a valuable biomarker to determine the initial treatment strategies. Of interest, circCYP24A1 could be measured before the treatment, suggesting it to be a useful biomarker for decision making of the treatment options.

Conclusions

Our findings illustrated that circCYP24A participated in the generation and maintenance of DTX resistance potently, suggesting that circCYP24A1 could function as a crucial and promising biomarker to predict DTX treatment response and a therapeutic target in patients who are resistant to DTX chemotherapy (Fig. 8).

Fig. 8
figure 8

Schematic diagram for the mechanisms of circCYP24A1 functioning as a miRNA sponge to promote DTX resistance in PCa through circCYP24A1/miR-1301-3p/ALDH1A3/PI3K/AKT/mTOR signaling pathway

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

PCa:

Prostate cancer

DTX:

Docetaxel

mHPSC:

Metastatic hormone-sensitive prostate cancer

mCRPC:

Metastatic castration-resistant prostate cancer

ncRNAs:

Noncoding RNAs

circRNAs:

Circular RNAs

miRNAs:

MicroRNAs

siRNAs:

Small interfering RNAs

cDNA:

Complementary DNA

HRP:

Horseradish peroxidase

qRT-PCR:

Quantitative realtime PCR

CCK8:

Cell Counting Kit 8

TUNEL:

TdT-mediated dUTP nick end labeling

RIP:

RNA immunoprecipitation

IHC:

Immunohistochemistry

PI:

Propidium iodide

FISH:

Fluorescence in situ hybridization

ROC:

Receiver operating characteristic

AUC:

Area under the curve

ceRNA:

Competing endogenous RNA

MRI:

Magnetic resonance imaging

pCR:

Pathologic complete response

MDR:

Minimal residual disease

References

  1. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71:7–33.

    PubMed  Article  Google Scholar 

  2. Chang AJ, Autio KA, Roach M 3rd, Scher HI. High-risk prostate cancer-classification and therapy. Nat Rev Clin Oncol. 2014;11:308–23.

    PubMed  PubMed Central  Article  Google Scholar 

  3. Cornford P, van den Bergh RCN, Briers E, Van den Broeck T, Cumberbatch MG, De Santis M, et al. EAU-EANM-ESTRO-ESUR-SIOG guidelines on prostate Cancer. Part II-2020 update: treatment of relapsing and metastatic prostate Cancer. Eur Urol. 2021;79:263–82.

    CAS  PubMed  Article  Google Scholar 

  4. Tannock IF, de Wit R, Berry WR, Horti J, Pluzanska A, Chi KN, et al. Docetaxel plus prednisone or mitoxantrone plus prednisone for advanced prostate cancer. N Engl J Med. 2004;351:1502–12.

    CAS  PubMed  Article  Google Scholar 

  5. Kyriakopoulos CE, Chen YH, Carducci MA, Liu G, Jarrard DF, Hahn NM, et al. Chemohormonal therapy in metastatic hormone-sensitive prostate Cancer: long-term survival analysis of the randomized phase III E3805 CHAARTED trial. J Clin Oncol. 2018;36:1080–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  6. James ND, Sydes MR, Clarke NW, Mason MD, Dearnaley DP, Spears MR, et al. Addition of docetaxel, zoledronic acid, or both to first-line long-term hormone therapy in prostate cancer (STAMPEDE): survival results from an adaptive, multiarm, multistage, platform randomised controlled trial. Lancet. 2016;387:1163–77.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  7. Gravis G, Fizazi K, Joly F, Oudard S, Priou F, Esterni B, et al. Androgen-deprivation therapy alone or with docetaxel in non-castrate metastatic prostate cancer (GETUG-AFU 15): a randomised, open-label, phase 3 trial. Lancet Oncol. 2013;14:149–58.

    CAS  PubMed  Article  Google Scholar 

  8. Armstrong AJ, Szmulewitz RZ, Petrylak DP, Holzbeierlein J, Villers A, Azad A, et al. ARCHES: a randomized, phase III study of androgen deprivation therapy with Enzalutamide or placebo in men with metastatic hormone-sensitive prostate Cancer. J Clin Oncol. 2019;37:2974–86.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  9. Wenzel M, Wurnschimmel C, Nocera L, Colla Ruvolo C, Tian Z, Shariat SF, et al. Overall survival after systemic treatment in high-volume versus low-volume metastatic hormone-sensitive prostate Cancer: systematic review and network Meta-analysis. Eur Urol Focus. 2021;8:399-408.

  10. Hodgson DR, Dougherty BA, Lai Z, Fielding A, Grinsted L, Spencer S, et al. Candidate biomarkers of PARP inhibitor sensitivity in ovarian cancer beyond the BRCA genes. Br J Cancer. 2018;119:1401–9.

  11. Tutt ANJ, Garber JE, Kaufman B, Viale G, Fumagalli D, Rastogi P, et al. Adjuvant Olaparib for patients with BRCA1- or BRCA2-mutated breast Cancer. N Engl J Med. 2021;384:2394–405.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. Hyman DM, Smyth LM, Donoghue MTA, Westin SN, Bedard PL, Dean EJ, et al. AKT inhibition in solid tumors with AKT1 mutations. J Clin Oncol. 2017;35:2251–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. Amaral PP, Dinger ME, Mercer TR, Mattick JS. The eukaryotic genome as an RNA machine. Science. 2008;319:1787–9.

    CAS  PubMed  Article  Google Scholar 

  14. Mattick JS. RNA regulation: a new genetics? Nat Rev Genet. 2004;5:316–23.

    CAS  PubMed  Article  Google Scholar 

  15. Zheng Q, Bao C, Guo W, Li S, Chen J, Chen B, et al. Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nat Commun. 2016;7:11215.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  16. Yu CY, Li TC, Wu YY, Yeh CH, Chiang W, Chuang CY, et al. The circular RNA circBIRC6 participates in the molecular circuitry controlling human pluripotency. Nat Commun. 2017;8:1149.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  17. Chen N, Zhao G, Yan X, Lv Z, Yin H, Zhang S, et al. A novel FLI1 exonic circular RNA promotes metastasis in breast cancer by coordinately regulating TET1 and DNMT1. Genome Biol. 2018;19:218.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  18. Hsiao KY, Lin YC, Gupta SK, Chang N, Yen L, Sun HS, et al. Noncoding effects of circular RNA CCDC66 promote Colon Cancer growth and metastasis. Cancer Res. 2017;77:2339–50.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. Suzuki H, Zuo Y, Wang J, Zhang MQ, Malhotra A, Mayeda A. Characterization of RNase R-digested cellular RNA source that consists of lariat and circular RNAs from pre-mRNA splicing. Nucleic Acids Res. 2006;34:e63.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  20. Yang L, Chen Y, Liu N, Lu Y, Ma W, Yang Z, et al. CircMET promotes tumor proliferation by enhancing CDKN2A mRNA decay and upregulating SMAD3. Mol Cancer. 2022;21:23.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. Wang S, Su W, Zhong C, Yang T, Chen W, Chen G, et al. An eight-CircRNA assessment model for predicting biochemical recurrence in prostate Cancer. Front Cell Dev Biol. 2020;8:599494.

    PubMed  PubMed Central  Article  Google Scholar 

  22. Xu J, Ji L, Liang Y, Wan Z, Zheng W, Song X, et al. CircRNA-SORE mediates sorafenib resistance in hepatocellular carcinoma by stabilizing YBX1. Signal Transduct Target Ther. 2020;5:298.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  23. Weng W, Wei Q, Toden S, Yoshida K, Nagasaka T, Fujiwara T, et al. Circular RNA ciRS-7-a promising prognostic biomarker and a potential therapeutic target in colorectal Cancer. Clin Cancer Res. 2017;23:3918–28.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  24. O'Neill AJ, Prencipe M, Dowling C, Fan Y, Mulrane L, Gallagher WM, et al. Characterisation and manipulation of docetaxel resistant prostate cancer cell lines. Mol Cancer. 2011;10:126.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  25. McKay RR, Ye H, Xie W, Lis R, Calagua C, Zhang Z, et al. Evaluation of intense androgen deprivation before prostatectomy: a randomized phase II trial of Enzalutamide and Leuprolide with or without Abiraterone. J Clin Oncol. 2019;37:923–31.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  26. Kristensen LS, Andersen MS, Stagsted LVW, Ebbesen KK, Hansen TB, Kjems J. The biogenesis, biology and characterization of circular RNAs. Nat Rev Genet. 2019;20:675–91.

    CAS  PubMed  Article  Google Scholar 

  27. Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in Drosophila. Genome Biol. 2003;5:R1.

    PubMed  PubMed Central  Article  Google Scholar 

  28. Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R. Fast and effective prediction of microRNA/target duplexes. RNA. 2004;10:1507–17.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. Chiang HR, Schoenfeld LW, Ruby JG, Auyeung VC, Spies N, Baek D, et al. Mammalian microRNAs: experimental evaluation of novel and previously annotated genes. Genes Dev. 2010;24:992–1009.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. Ha M, Kim VN. Regulation of microRNA biogenesis. Nat Rev Mol Cell Biol. 2014;15:509–24.

    CAS  PubMed  Article  Google Scholar 

  31. Wu W, Wu Y, Mayer K, von Rosenstiel C, Schecker J, Baur S, et al. Lipid peroxidation plays an important role in chemotherapeutic effects of Temozolomide and the development of therapy resistance in human Glioblastoma. Transl Oncol. 2020;13:100748.

    PubMed  PubMed Central  Article  Google Scholar 

  32. Xia ZK, Wang W, Qiu JG, Shi XN, Li HJ, Chen R, et al. Discovery of a new CDK4/6 and PI3K/AKT multiple kinase inhibitor Aminoquinol for the treatment of hepatocellular carcinoma. Front Pharmacol. 2021;12:691769.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  33. Wei QT, Liu BY, Ji HY, Lan YF, Tang WH, Zhou J, et al. Exosome-mediated transfer of MIF confers temozolomide resistance by regulating TIMP3/PI3K/AKT axis in gliomas. Mol Ther Oncolytics. 2021;22:114–28.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. Nie S, Qian X, Shi M, Li H, Peng C, Ding X, et al. ALDH1A3 accelerates pancreatic Cancer metastasis by promoting glucose metabolism. Front Oncol. 2020;10:915.

    PubMed  PubMed Central  Article  Google Scholar 

  35. Petrylak DP, Tangen CM, Hussain MH, Lara PN Jr, Jones JA, Taplin ME, et al. Docetaxel and estramustine compared with mitoxantrone and prednisone for advanced refractory prostate cancer. N Engl J Med. 2004;351:1513–20.

    CAS  PubMed  Article  Google Scholar 

  36. Huang X, Li Z, Zhang Q, Wang W, Li B, Wang L, et al. Circular RNA AKT3 upregulates PIK3R1 to enhance cisplatin resistance in gastric cancer via miR-198 suppression. Mol Cancer. 2019;18:71.

    PubMed  PubMed Central  Article  Google Scholar 

  37. Han B, Chao J, Yao H. Circular RNA and its mechanisms in disease: from the bench to the clinic. Pharmacol Ther. 2018;187:31–44.

    CAS  PubMed  Article  Google Scholar 

  38. Sossey-Alaoui K, Plow EF. miR-138-mediated regulation of KINDLIN-2 expression modulates sensitivity to chemotherapeutics. Mol Cancer Res. 2016;14:228–38.

    CAS  PubMed  Article  Google Scholar 

  39. Xu B, Niu X, Zhang X, Tao J, Wu D, Wang Z, et al. miR-143 decreases prostate cancer cells proliferation and migration and enhances their sensitivity to docetaxel through suppression of KRAS. Mol Cell Biochem. 2011;350:207–13.

    CAS  PubMed  Article  Google Scholar 

  40. Qiao DH, He XM, Yang H, Zhou Y, Deng X, Cheng L, et al. miR-1301-3p suppresses tumor growth by downregulating PCNA in thyroid papillary cancer. Am J Otolaryngol. 2021;42:102920.

    CAS  PubMed  Article  Google Scholar 

  41. Peng X, Yan B, Shen Y. MiR-1301-3p inhibits human breast cancer cell proliferation by regulating cell cycle progression and apoptosis through directly targeting ICT1. Breast Cancer. 2018;25:742–52.

    PubMed  Article  Google Scholar 

  42. Sutendra G, Kinnaird A, Dromparis P, Paulin R, Stenson TH, Haromy A, et al. A nuclear pyruvate dehydrogenase complex is important for the generation of acetyl-CoA and histone acetylation. Cell. 2014;158:84–97.

    CAS  PubMed  Article  Google Scholar 

  43. Liu X, Cooper DE, Cluntun AA, Warmoes MO, Zhao S, Reid MA, et al. Acetate production from glucose and coupling to mitochondrial metabolism in mammals. Cell. 2018;175:502–13 e13.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. Mashimo T, Pichumani K, Vemireddy V, Hatanpaa KJ, Singh DK, Sirasanagandla S, et al. Acetate is a bioenergetic substrate for human glioblastoma and brain metastases. Cell. 2014;159:1603–14.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  45. Wellen KE, Hatzivassiliou G, Sachdeva UM, Bui TV, Cross JR, Thompson CB. ATP-citrate lyase links cellular metabolism to histone acetylation. Science. 2009;324:1076–80.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. Golubovskaya V, O'Brien S, Ho B, Heffler M, Conroy J, Hu Q, et al. Down-regulation of ALDH1A3, CD44 or MDR1 sensitizes resistant cancer cells to FAK autophosphorylation inhibitor Y15. J Cancer Res Clin Oncol. 2015;141:1613–31.

    CAS  PubMed  Article  Google Scholar 

  47. Kong B, Wu W, Cheng T, Schlitter AM, Qian C, Bruns P, et al. A subset of metastatic pancreatic ductal adenocarcinomas depends quantitatively on oncogenic Kras/Mek/Erk-induced hyperactive mTOR signalling. Gut. 2016;65:647–57.

    CAS  PubMed  Article  Google Scholar 

  48. Kang YJ, Jang WS, Kwon JK, Yoon CY, Lee JY, Ham WS, et al. Intermediate PSA half-life after neoadjuvant hormone therapy predicts reduced risk of castration-resistant prostate cancer development after radical prostatectomy. BMC Cancer. 2017;17:789.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  49. McGuire SE, Lee AK, Cerne JZ, Munsell MF, Levy LB, Kudchadker RJ, et al. PSA response to neoadjuvant androgen deprivation therapy is a strong independent predictor of survival in high-risk prostate cancer in the dose-escalated radiation therapy era. Int J Radiat Oncol Biol Phys. 2013;85:e39–46.

    CAS  PubMed  Article  Google Scholar 

  50. Zilli T, Dal Pra A, Kountouri M, Miralbell R. Prognostic value of biochemical response to neoadjuvant androgen deprivation before external beam radiotherapy for prostate cancer: a systematic review of the literature. Cancer Treat Rev. 2016;46:35–41.

    CAS  PubMed  Article  Google Scholar 

Download references

Acknowledgments

We are grateful to Dr. C. P. (Immunology and Reproduction Biology Laboratory & State Key Laboratory of Analytical Chemistry for Life Science, Medical School, Nanjing University) for providing technical assistance in animal study.

Funding

This study was supported by grants from the Project of Invigorating Health Care through Science, Technology and Education, Jiangsu Provincial Key Medical Discipline (Laboratory) (ZDXKB2016014). National Natural Science Foundation of China (82172639, 82102986).

Author information

Authors and Affiliations

Authors

Contributions

HY, HQ and LY performed experiments and data analysis; WZ, YY and QL provided technical assistance; HG, XQ and HQ provided scientific input and designed experiments; HY and LY wrote and final proofed the manuscript; MC designed and supervised the study; JH and JS verified the underlying data. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Xuefeng Qiu or Hongqian Guo.

Ethics declarations

Ethics approval and consent to participate

This study was undertaken with the understanding and written consent of each subject, with approval of the Ethics Committee of Drum Tower Hospital, Medical School of Nanjing University (Nanjing, China). Informed consent was obtained from all participants included in this study according to ethical committee regulations. Animal experiments were approved and performed in accordance with the Institutional Animal Care and Use Committee of Drum Tower Hospital, Medical School of Nanjing University.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

DTX-resistant cell lines (DU145-DR) were established. A. CCK-8 assays were performed to detected the DTX cytotoxicity in DU145 and DU145-DR cells. B. Tumor growth curves of tumor sizes are shown. Tumor volumes were measured every 3 days since the first DTX treatment. C. The weight growth curves of every group mouse. Mice weight were measured every 3 days during treatment. D, F. DU145 and DU145-DR cells were treated with DTX (20 nM) or PBS (control) and subjected to Annexin V-APC and PI stanning to detect apoptotic rate by flow cytometry. E, G. Cell cycle were detected by flow cytometry in indicated cells after treatment with DTX (20 nM). All data are presented as the means ± SD of three independent experiments. ***P < 0.001.

Additional file 2: Figure S2.

RNA sequencing analysis were conducted to investigate the differential expression of circRNAs in DTX-resistant DU145 cells. A. Hot map visualizing 107 upregulated and 67 downregulated circRNAs in DU145-DR cells according to the log2 (fold-change) > 1 or < − 1 and P value < 0.05. B. Pie chart showed the type of those abnormally expressed circRNAs.

Additional file 3: Figure S3.

The statistic diagram showed the quantification of the band intensity of ALDH1A3 protein level in in indicated cells. A. ALDH1A3 protein level after knocking down circCYP24A1 expression in DU145-DR cells. B. ALDH1A3 protein level after upregulated circCYP24A1 expression in DU145 and 22RV1 cells. C. ALDH1A3 protein level after transfected with shcircCYP24A1 and miR1301-3p inhibitor in DU145-DR cells.

Additional file 4: Figure S4.

Representative images of IHC (ALDH1A3) staining of the xenograft tumors.

Additional file 5: Figure S5.

Two ALDH1A3 siRNAs were used to knock-down ALDH1A3 expression. Western blot was performed to validate si-ALDH1A3 efficiency.

Additional file 6: Figure S6.

Scatter plot shows the correlation among the level of circCYP24A1, miR-1301-3p and ALDH1A3. A. CircCYP24A1 with miR-1301-3p. (R = − 0.45, P < 0.001) B. circCYP24A1 with ALDH1A3. (R = 0.39, P < 0.001) C. MiR-1301-3p with ALDH1A3. (R = − 0.33, P = 0.0047).

Additional file 7: Figure S7.

ROC analysis to compare the diagnostic accuracy of the identified parameters. PSA = prostate specific antigen; AUC = area under the curve; 95% confidence intervals is labeled in parentheses.

Additional file 8: Table S1.

Primers used for real-time PCR.

Additional file 9: Table S2.

Probes used for RNA pull-down.

Additional file 10: Table S3.

Primary antibodies used in this study.

Additional file 11: Table S4.

Probes used for FISH assay.

Additional file 12: Table S5.

SiRNA used for silencing target genes.

Additional file 13: Table S6.

Baseline characteristics of the included 70 high-risk PCa patients who received neoadjuvant therapy with DTX.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yin, H., Qin, H., Yang, L. et al. circCYP24A1 promotes Docetaxel resistance in prostate Cancer by Upregulating ALDH1A3. Biomark Res 10, 48 (2022). https://doi.org/10.1186/s40364-022-00393-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40364-022-00393-1

Keywords

  • Prostate cancer
  • circCYP24A1
  • miR-1301-3p
  • ALDH1A3
  • Docetaxel resistance