- Open Access
Hsa_circ_0004296 inhibits metastasis of prostate cancer by interacting with EIF4A3 to prevent nuclear export of ETS1 mRNA
Journal of Experimental & Clinical Cancer Research volume 40, Article number: 336 (2021)
Circular RNAs (circRNAs) have been shown to play vital biological functions in various tumors, including prostate cancer (PCa). However, the roles of circRNAs in the metastasis of PCa remain unclear. In the present study, differentially expressed circRNAs associated with PCa metastasis were screened using high-throughput RNA sequencing, from which hsa_circ_0004296 was identified.
Quantitative real-time PCR (qRT-PCR) was used to detect the expression of circ_0004296 in PCa tissues and adjacent normal tissues as well as in blood and urine. Gain and loss of function experiments were performed to investigate the function of circ_0004296 in PCa. Bioinformatics analyses, RNA pull-down assay, and mass spectrometry were conducted to identify RNA-binding proteins. RNA immunoprecipitation and RNA and protein nuclear-cytoplasmic fractionation were performed to investigate the underlying mechanism. A xenograft mouse model was used to analyze the effect of circ_0004296 on PCa growth and metastasis in vivo.
The expression of circ_0004296 was decreased in PCa tissues, blood, and urine, which was negatively associated with metastasis. Furthermore, gain and loss of function experiments in vitro and in vivo showed that circ_0004296 inhibited the proliferation, migration, invasion, and epithelial-mesenchymal transition of PCa cells. Mechanistically, circ_0004296 regulated host gene ETS1 expression at the post-transcriptional level. EIF4A3 was identified and confirmed as the downstream binding protein of circ_0004296. EIF4A3 expression was significantly upregulated in PCa tissues and associated with PCa metastasis. Silencing EIF4A3 suppressed PCa cell proliferation, migration, invasion, and EMT.
Circ_0004296 overexpression efficiently inhibited ETS1 mRNA nuclear export by promoting EIF4A3 retention in the nucleus, leading to the downregulation of ETS1 expression and suppression of PCa metastasis; thus, circ_0004296 might be a potential biomarker and therapeutic target for patients with PCa.
Prostate cancer (PCa) is the second most common malignancy in men and causes around 300,000 deaths worldwide annually . Metastasis is the leading cause of death from PCa . Median survival of patients with newly diagnosed metastasis is approximately 42 months , while that of patients with progressive metastatic PCa is approximately 27 months . Therefore, it is important to further clarify the molecular mechanism of PCa metastasis to identify novel biomarkers and develop effective therapeutic targets.
Circular RNAs (circRNAs) are a class of covalently closed circular noncoding RNAs without 5′ and 3′ ends . Although circRNAs were initially considered to be splicing-associated noise, with the advance in high-throughput sequencing methods, a large number of circRNA candidates were found to have potential biological functions. Accumulated evidence has confirmed abnormal expression of circRNAs in multiple tumors, including PCa, wherein they function as oncogenes or tumor suppressor genes . CircRNAs are involved in the malignant behavior of tumor cells, such as proliferation, apoptosis, migration, and invasion [7,8,9]. At the molecular level, circRNAs play potentially important roles in tumorigenesis and tumor progression through competitive endogenous RNAs (ceRNAs), encoded proteins, interaction with proteins, and modulation of transcription [10,11,12,13]. Recent emerging evidence also suggests that circRNAs plays important roles in PCa, with the findings indicating their involvement in tumor cell proliferation, invasion, migration, and epithelial-mesenchymal transition (EMT) [14, 15]. However, the functions and mechanisms of most circRNAs remain largely unexplored in PCa metastasis.
The EMT process has been found to play a crucial role in PCa metastasis, wherein epithelial cells acquire mesenchymal phenotype along with loss of intercellular adhesion and conversion to migratory and invasive cells . It is well-known that the induction and maintenance of EMT are closely controlled by signaling pathways and transcription factors, which involve a complex regulatory network of genes and noncoding RNAs . Recent studies have also implicated circRNAs in cancer development, wherein they affect tumor cell migration, invasion, and cancer metastasis by regulating EMT [14, 18]. Previous studies have shown that the transcription factor ETS1 is involved in tumor progression and metastasis through EMT induction in various cancers, including hepatocellular carcinoma and colon cancer [19, 20].
In the present study, we identified circ_0004296 by RNA sequencing and qPCR and evaluated its role in PCa. The expression of circ_0004296 was decreased in PCa tissues, blood, and urine, and it was negatively associated with metastasis. Furthermore, gain and loss of function experiments in vitro and in vivo showed that circ_0004296 inhibited the proliferation, migration, invasion, and EMT of PCa cells. Mechanistically, circ_0004296 was localized in the nucleus, interacted with the RNA-binding protein (RBP) EIF4A3 to prevent nuclear export of host gene ETS1 mRNA, and subsequently inhibited its expression. Our study implied that circ_0004296 may serve as a specific biomarker and therapeutic target for patients with metastatic PCa (mPCa).
Materials and methods
All patients provided written informed consent. The present study was approved by the Ethics Committee of Shanghai Tenth People’s Hospital (approval number: SHSY-IEC-2014RES-36). Clinical pathological data were also collected through the hospital medical record system. Five pairs of PCa tissues and matched local metastatic lymph node tissues were obtained from radical resection of PCa and dissection of enlarged lymph nodes. Forty pairs of PCa tissues and matched adjacent normal tissues were obtained from radical resection of PCa. Forty-six urine specimens and 39 blood samples were obtained from patients with benign prostatic hyperplasia, localized PCa, and metastatic PCa. Tissue microarrays included 359 specimens from patients with PCa. These cases had corresponding follow-up information for an average of 32 months. All cases were confirmed by clinical and pathological diagnosis.
Whole-transcriptome sequencing (RNA-seq)
Five pairs of PCa tissues and matched local metastatic lymph node tissues were used for RNA-seq. Total RNA was isolated using the QubitRNA Assay Kit (cat. Q32852; Life Technologies, USA) following the manufacturer’s protocol. RNA-seq and data analysis were performed by Oebiotech (Shanghai OEbiotech Co., Ltd., Shanghai, China).
PCa cell lines (PC3, DU145, LNCaP, and 22RV-1) and the normal human prostate epithelial cell line RWPE-1 were purchased from the cell library of Shanghai Chinese Academy of Sciences (Shanghai, China). PCa cell lines were cultured in RPMI-1640 medium, supplemented with 10% fetal bovine serum (FBS), penicillin (100 U/ml), and streptomycin (100 μg/ml). RWPE-1 cells were cultured in Defined Keratinocyte SFM (1X) medium (Gibco, USA).
Quantitative real-time PCR (qRT-PCR)
Total RNA for PCa tissues and cells were isolated using TRIzol as described previously . Total RNA for urine and blood samples were isolated by magnetic bead nucleic acid extraction kit following the manufacturer’s instructions (Newskybio, China). According to the manufacturer’s instructions, RNA was reverse transcribed into cDNA using the HiScript® III 1st Strand cDNA Synthesis Kit (Vazyme, China), the qPCR was completed using ChamQ SYBR qPCR Master Mix Kit (Vazyme, China). The primers (Supplementary Table S1) were synthesized by Sangon (Shanghai, China).
PC3 and DU145 cells were transfected using Lipofectamine® 3000 (Invitrogen; Thermo Fisher Scientific, Inc., USA) following the manufacturer’s instructions. Plasmid and lentivirus expression vector constructs for overexpressing circ_0004296 were designed and synthesized by Zuorun (Shanghai, China) (Supplementary Table S1). Plasmid constructs for overexpressing ETS1 were designed and synthesized by Qihe (Shanghai, China). Small interfering RNAs (siRNAs) for the indicated gene sequences (Supplementary Table S1) were designed and synthesized by Ibsbio (Shanghai, China). DU145 cells were subsequently infected and selected with puromycin. The luciferase was then incorporated into DU145- circ_0004296 stable cells.
Cell proliferation assays
Cell proliferation was assessed by the Cell Counting Kit-8 (CCK8) assay and expressed as colony formation as described previously . Three independent replicates were set in each experiment; and all experiments were conducted in triplicate.
Wound healing assay
Cell migration was measured by the wound healing assay as described previously . Three independent replicates were set in each experiment; and all experiments were conducted in triplicate.
Cell migration and invasion were assessed by the Transwell migration and invasion assay as described previously . Three independent replicates were set in each experiment; and all experiments were conducted in triplicate.
Western blotting assay
Total protein was lysed from cells using radioimmunoprecipitation assay buffer (RIPA buffer; Beyotime, Shanghai, China). Protein quantification and western blotting assays were performed as described previously . The following primary antibodies were used: anti-E-cadherin (WL01482, Wanleibio, China), anti-ETS1 (#14069, Cell Signaling Technology (CST), USA), anti-N-cadherin (WL01047, Wanleibio, China), anti-Snail (#3879, CST), anti-EIF4A3 (ab18057, Abcam, USA), and anti-Vimentin (ab8069, Abcam).
Protein expression was evaluated by immunofluorescence staining (IF) of PCa tissues and cells. Briefly, 5-μm paraffin-embedded cross-sections of tissues or cells uniformly grown on slides were fixed with 4% paraformaldehyde, permeabilized with 0.1% Triton-100, and blocked with 5% bovine serum albumin (BSA). Subsequently, tissues or cells were incubated with anti-E-cadherin (WL01482, Wanleibio), anti-EIF4A3 (ab18057, Abcam, USA), and anti-Vimentin (ab8069, Abcam) overnight at 4 °C. Then, the cells were treated with the corresponding secondary antibody after washing for three times. DAPI (Beyotime) staining was used for nuclear localization. Images were captured with a confocal microscope (Leica Microsystems, Mannheim, Germany).
All animal experiments were approved by the Animal Care and Use Committee of Shanghai Tenth People’s Hospital of Tongji University (approval number: SHDSYY-2014-3028). BABL/c male nude mice were housed under specific pathogen-free conditions. To develop a subcutaneous xenograft mouse model, twelve 4-week-old male nude mice were randomly divided into two groups, which were subcutaneously injected with 5 × 106 DU145-vector or DU145-circ_0004296 cells. Tumor size was measured using a vernier caliper twice a week. Tumor volume was calculated as follows: 0.5 × length × width2. Five weeks later, mice were sacrificed, and all dissected tumor xenografts were weighed and subjected to immunohistochemistry (IHC). IHC was performed as described previously . To develop an orthotopic xenograft mouse model, 8-week-old male nude mice received intraprostatic injections of 2 × 105 DU145-vector or DU145-circ_0004296 cells in 10 μL of PBS, with six mice in each group. To develop a pulmonary metastatic model, 1 × 106 DU145-vector or DU145-circ_0004296 cells were injected into the tail vein of 8-week-old male nude mice, with six mice in each group. The intraprostatic and metastatic tumors were observed weekly using the in vivo imaging system (IVIS). Mice were sacrificed when the weight loss was ≥10%, and the lung metastatic foci were examined microscopically by H&E staining.
Fractionation of nuclear and cytoplasmic RNA
Nuclear and cytoplasmic RNA were extracted using the PARIS™ kit (Invitrogen, USA) following the manufacturer’s protocol. The qRT-PCR was used to detect RNA abundance in different cell fractions. The efficiency of nuclear and cytoplasmic RNA isolation was controlled by qRT-PCR using U6 and actin, respectively.
Fractionation of nuclear and cytoplasmic protein
Nuclear and cytoplasmic protein fractionation assays were performed with nuclear and cytoplasmic protein extraction kit (Beyotime) following the manufacturer’s protocol. The purity of the subcellular fractions were detected using anti-GAPDH (ab18602, Abcam) and anti-Histone 3 (13108-1-AP, Proteintech, China) antibodies by western blotting assay.
RNA pull-down assay
RNA pull-down assay for circ_0004296 and circ_0004296-binding proteins was performed using RNA pull-down kit (Bersinbio, Guangzhou, China). The biotinylated circ_0004296 probe used to pull down for circ_0004296 and RBPs was synthesized by GenePharma (Shanghai, China). DU145 cells were harvested and lysed with lysis buffer. The probes were incubated with streptavidin-coated magnetic beads to generate probe-coupled magnetic beads. The cell extract (800 μL) was incubated with magnetic beads coupled with the circ_0004296 probe or a negative control (NC) probe. Proteins pulled down by the probe were collected using a protein elution buffer. Mass spectrometry was used to identify differentially interacted proteins. Mass spectrometry and data analysis were performed by Wayenbio (Shanghai, China). The differential abundance of circ_0004296-binding proteins was confirmed by western blotting assay.
RNA fluorescence in situ hybridization analysis (RNA-FISH)
RNA-FISH experiments were completed using a Ribo™ FISH kit (Ribobio, China) following the manufacturer’s instructions. The FISH probe of circ_0004296 was synthesized by GenePharma (Supplementary Table S1). An 18S RNA probe was used as a positive control (Ribobio).
RNA immunoprecipitation assay (RIP)
RIP was completed with an EZ-Magna RIP kit (Millipore, Billerica, MA, USA). DU145 cells were harvested and lysed with lysis buffer. The cell extracts were incubated with magnetic beads coated with anti-EIF4A3 antibodies (ab180573, Abcam) or control (IgG, Beyotime). Subsequently, total RNA was extracted from the magnetic beads. qRT-PCR was used to determine the presence of ETS mRNA and circ_0004296.
Data were analyzed and plotted using SPSS 25.0 (IBM, SPSS, Chicago, IL, USA) and GraphPad Prism 7 (GraphPad Prism Software Inc., San Diego, CA), which were shown as mean ± standard deviation (SD). Student’s t-test was used to evaluate the difference between two groups, while ANOVA test was used for more than two groups. All experiments were conducted in triplicate. Linear correlations were analyzed by Pearson’s correlation coefficient. The disease-free survival and biochemical recurrence-free survival was analyzed by Kaplan-Meier method and the log-rank test. P < 0.05 was considered to be statistically significant.
Circ_0004296 is negatively associated with PCa metastasis
To investigate whether circRNAs are involved in PCa metastasis, five pairs of PCa tissues and matched local metastatic lymph node tissues were used for RNA-seq. The results showed a total of 22,676 circRNAs. Among these circRNAs, 208 upregulated circRNAs and 230 downregulated circRNAs were identified to show significant differential expression (fold change > 2.0 and P < 0.05) (Fig. 1A). Next, we designed divergent primers to detect six screened circRNA candidates in PCa cell lines, prostate normal epithelial cells, PCa tissues, and adjacent normal tissues (Supplementary Table S1). We chose circRNA_04265 (circBase ID: hsa_circ_0004296) as the candidate circRNA for further research, which is derived from back-splicing of exons 4, 5, 6, and 7 of host gene ETS1, with a length of 648 bp. The reasons were as followings: First, the circRNAs that are smaller than 2000 bp and annotated by the circBase database is selected. Then, circ_0004296 had relatively high expression abundance (Figure S1A); compared to that in prostate normal cells, circ_0004296 expression was significantly downregulated in all PCa cell lines (Fig. 1B and S1B-D); its host gene ETS1 was abnormally high expression and is an oncogenic gene in PCa, which suggests that the low expression of circ_0004296 may be functional. Consistently, compared to that in the adjacent normal tissue (n = 40), circ_0004296 was reduced in PCa tissues (Fig. 1C). Circ_0004296 was detected stably in urine (Fig. 1D) and blood plasma (Fig. 1E) from patients with benign prostatic hyperplasia, localized PCa, and metastatic PCa. Furthermore, the expression of circ_0004296 was negatively associated with distant metastasis (Fig. 1D and E). According to the consistent expression observed in PCa cell lines and PCa tissues by qRT-PCR, downregulated circ_0004296 might be involved in the occurrence and metastasis of PCa. We then designed divergent and convergent primers to amplify circ_0004296 and linear ETS1 mRNA by using cDNA and genomic DNA (gDNA) from PCa cells. Northern blotting showed that circ_0004296 can only be detected in the cDNA group, but not in the gDNA group (Fig. 1F). Sanger sequencing confirmed the back-spliced junction of circ_0004296 (Fig. 1G). To further confirm the stably expressed circ_0004296, PCa cells were treated with Rnase R and actinomycin D. The results of qRT-PCR showed that circ_0004296 had better stability than linear mRNA (Fig. 1H and I). These results suggested that circ_0004296 may play a suppressive role in PCa metastasis and serves as a noninvasive liquid biopsy biomarker.
Circ_0004296 inhibited PCa cell proliferation, migration, invasion, and EMT in vitro
To investigate the function of circ_0004296 in PCa cells, the recombinant pLenO-GTP-circ_0004296-overexpression plasmid was constructed. The integrity of the linear circ_0004296 RNA sequence carried on the plasmid was confirmed by northern blotting and Sanger sequencing (Fig. 2A and S1E, Supplementary file 1). Overexpression of circ_0004296 significantly increased the level of circ_0004296 but had no effect on the linear mRNA level in PC3 and DU145 cells (Fig. 2B). The proliferation assay, scratch wound healing, and Transwell assay revealed that overexpression of circ_0004296 markedly inhibited proliferation, migration, and invasion of PC3 and DU145 cells compared to the negative control (Fig. 2C-H).
It is known that EMT plays a critical role in mediating the metastasis of human tumors, including PCa . Next, IF and western blotting assays were performed to examine the impact of circ_0004296 on the expression of epithelial and mesenchymal markers. The results showed that circ_0004296 overexpression remarkably increased epithelial marker (E-cadherin) expression but reduced mesenchymal marker (N-cadherin, Vimentin and Snail) expression in PCa cells (Fig. 2I and J).
We also designed three siRNAs targeting the back-spliced junction (si-circ_0004296). The qRT-PCR results showed that siRNA#1 had the highest knockdown efficiency for circ_0004296 and had no effect on the linear mRNA level (Figure S2A and B). Furthermore, we found that circ_0004296 knockdown significantly promoted the proliferation, migration, and invasion of PCa cells (Figure S2C-H). Moreover, circ_0004296 knockdown remarkably decreased epithelial marker (E-cadherin) expression but increased mesenchymal marker (N-cadherin, Vimentin and Snail) expression in PCa cells (Figure S2I). Taken together, these findings demonstrated that circ_0004296 inhibited PCa cell proliferation, migration, invasion, and EMT.
Circ_0004296 regulated host gene ETS1 expression through interaction with the RBP EIF4A3
First, the intracellular localization of circ_0004296 was detected in PCa cells by the RNA-FISH assay and the nuclear and cytoplasmic fractionation assay. The results showed the predominant nuclear distribution of circ_0004296 (Fig. 3A and B). Overexpression of circ_0004296 in PCa cell lines remarkably decreased the ETS1 protein level, but had no significant impact on the expression of ETS1 mRNA (Fig. 3C). Additionally, silencing circ_0004296 remarkably increased the ETS1 protein level, but had no obvious effect on the level of ETS1 mRNA (Fig. 3D). These results suggested that circ_0004296 regulated host gene ETS1 expression at the post-transcriptional level. Recent studies have shown the significance of circRNA-protein interaction . Some RBPs have been identified to have an interplay with circRNAs through the regulation of gene expression in cancers [10, 15]. Given that circ_0004296 is mainly distributed in the nucleus, we investigated whether circ_0004296 regulated host gene ETS1 expression through interaction with an RBP by performing the RNA pull-down assay and mass spectrometry analysis (Fig. 3E, F and G). The specific proteins that interacted with circ_0004296 are presented in Supplementary Table S2. From the prediction results of CircInteractome, Starbase, and CircFunBase, we selected EIF4A3 as the downstream binding protein of circ_0004296 (Fig. 3F). RNA pull-down and western blotting assays confirmed the interaction of circ_0004296 with EIF4A3 (Fig. 3G). Silencing EIF4A3 remarkably decreased ETS1 protein and mRNA levels (Fig. 3H-J). These results suggested that circ_0004296 regulated host gene ETS1 expression by interacting with EIF4A3. To determine the potential functions of EIF4A3 in PCa, GSEA was performed to link the published molecular signature databases for PCa metastasis to high versus low EIF4A3 expression in TCGA databases (Gene Set: TOMLINS_METASTASIS_UP and SUNG_METASTASIS_STROMA_UP). GSEA showed that cell metastasis was significantly enriched in the PCa group, which strongly suggested that EIF4A3 was closely related to PCa metastasis (Fig. 4A). We also found that EIF4A3 was upregulated in PCa tissues and cells (Fig. 4B). Silencing EIF4A3 suppressed PCa cell proliferation, migration, invasion, and EMT (Fig. 4C-I). Kaplan-Meier analysis showed that EIF4A3 upregulation was associated with shorter biochemical recurrence-free survival of 359 patients with PCa (log-rank test, p = 0.017, Fig. 4J). Moreover, EIF4A3 expression was significantly upregulated in PCa tissues compared to that in normal samples from the TCGA database (Fig. 4K). Subsequently, Kaplan-Meier analysis suggested that EIF4A3 upregulation was associated with shorter disease-free survival of 492 patients with PCa from the TCGA database (log-rank test, p = 0.011, Fig. 4L). These results supported that EIF4A3 protein mediated circ_0004296 inhibition of PCa cell growth and metastasis. To investigate the function of ETS1 in PCa cells, three siRNAs targeting ETS1 were designed and the recombinant ETS1-overexpression plasmid was constructed. The results of western bolt showed that siRNA#1 had the highest knockdown efficiency for ETS1 (Figure S3B) and ETS1 expression was upregulated (Figure S3C). The proliferation assay and Transwell assay revealed that silencing ETS1 suppressed PCa cell proliferation and migration (Figure S3D, E). The proliferation assay and Transwell assay revealed that overexpression of ETS1 promoted proliferation and migration of PC3 and DU145 cells compared to the negative control (Figure S3F, G). Kaplan-Meier analysis showed that ETS1 upregulation was associated with shorter biochemical recurrence-free survival of 343 patients with PCa (log-rank test, p < 0.001, Figure S3H). The rescue experiment was performed to assess the effect of the circ_0004296/EIF4A3/ETS1 axis on PCa cells. The results showed that overexpression of ETS1 in pLenO-GTP-circ_0004296-overexpression cells could restore the function of circ_0004296-mediating inhibition of proliferation and migration (Fig. 5A and B).
Taken together, these findings demonstrated that circ_0004296 inhibited PCa cell proliferation, migration, invasion, and EMT through the EIF4A3/ETS1 axis.
Circ_0004296 inhibited nuclear export of ETS1 mRNA through binding to EIF4A3
Furthermore, GSEA supported that the mRNA splicing and mRNA export functional signatures were significantly enriched for the high versus low EIF4A3 expression in TCGA databases (Fig. 6A); this finding strongly suggested that EIF4A3 was closely related to ETS1 expression at the post-transcriptional level. To investigate the potential regulatory relationship between circ_0004296 and EIF4A3, the expression of EIF4A3 was measured after silencing or overexpression of circ_0004296. The results showed that silencing or overexpression of circ_0004296 had no effect on the expression of EIF4A3 in PC3 and DU145 cells (Fig. 6B-D). Additionally, the knockdown of EIF4A3 had no significant impact on circ_0004296 expression in PCa cells (Figure S3A). Previous studies have shown that the RBP EIF4A3, which is mainly located in the nucleus and is the core component of the exon junction complex (EJC), can facilitate nucleocytoplasmic transport of mRNA [27, 28]. Thus, we speculated that circ_0004296 might be involved in the regulation of nuclear ETS1 mRNA export in PCa cells. We therefore subjected circ_0004296-overexpressing cells to nuclear and cytoplasmic fractionation. The nuclear and cytoplasmic levels of EIF4A3 protein were detected by western blotting assay. The nuclear and cytoplasmic levels of ETS1 mRNA were determined by qRT-PCR. In circ_0004296-upregulated cells, the proportions of nuclear EIF4A3 protein and nuclear ETS1 mRNA were increased (Fig. 6E-H). RIP assay further demonstrated that EIF4A3 interacted with circ_0004296 and ETS1 mRNA in DU145 cells (Fig. 6I and J); this finding is consistent with the prediction of EIF4A3 binding to ETS1 mRNA by Starbase from HITS-CLIP seq . These data suggested that circ_0004296 efficiently inhibited ETS1 mRNA nuclear export by promoting EIF4A3 retention in the nucleus.
Circ_0004296 inhibited PCa growth and metastasis in vivo
To further examine the effect of circ_0004296 in vivo, we constructed luciferase-labeled stably DU145-overexpressing circ_0004296 and luciferase-labeled DU145-vector cells (Figure S3I and J). The results of the subcutaneous xenograft mouse model showed that overexpression of circ_0004296 inhibited PCa growth in vivo (Fig. 6A and B). The xenograft tumor was subjected to IHC assay. The results showed that the circ_0004296 overexpression downregulated the Ki67 protein level (Fig. 7C). The results of the orthotopic xenograft mouse model showed that overexpression of circ_0004296 suppressed PCa growth and pelvic peritoneal invasion (Fig. 7D). The results of IHC assay showed that circ_0004296 overexpression remarkably downregulated the mesenchymal marker (Vimentin) protein level but upregulated the epithelial marker (E-cadherin) protein level (Fig. 7E). Additionally, IHC showed that overexpression of circ_0004296 remarkably reduced the ETS1 protein level (Fig. 7F). IF was performed to confirm the change in the subcellular localization of EIF4A3 protein in xenograft tumor, and the results showed that the expression of EIF4A3 protein from the circ_0004296 overexpression group was increased in the nucleus as compared to that in the control group (Fig. 7G). Furthermore, The results of the in vivo tail vein injection model showed that overexpression of circ_0004296 suppressed metastases of PCa (Fig. 7H). The longevity of the mice post-injection further indicated circ_0004296 inhibited PCa metastasis in the intravenous xenograft assay (Fig. 7I). Taken together, these findings in vivo also demonstrated that circ_0004296 inhibited PCa cell growth and metastasis through the EIF4A3/ETS1 axis.
The present study first demonstrated that circ_0004296 was downregulated in PCa tissues and cell lines. We also confirmed that circ_0004296 was stably downregulated in the blood plasma and urine of patients with metastatic PCa, which predicted distant metastasis. Further in vitro and in vivo experiments confirmed that circ_0004296 inhibited the proliferation, migration, invasion, and EMT of PCa cells. Mechanistically, we showed that circ_0004296 was mainly located in the nucleus, and it interacted with RBP EIF4A3 and regulated the host gene ETS1 expression to inhibit the function of PCa cells. Furthermore, circ_0004296 efficiently inhibited nuclear export of ETS1 mRNA by promoting EIF4A3 retention in the nucleus (Fig. 8). Taken together, our study suggested that circ_0004296 might serve as a noninvasive liquid biopsy biomarker and a novel therapeutic target in human PCa.
Increasing evidence suggests that circRNAs are aberrantly expressed in multiple human tumors, and they play vital roles in tumor growth, tumor progression, metastasis, and drug resistance . Previous studies have shown that circ-AMOTL1L is expressed at a low level in PCa tissues and inhibit invasion and migration of PCa cells . In contrast, circ0005276 and circ-AGO2 are highly expressed in PCa and play a role in promoting cancer cell proliferation and migration [15, 29]. However, the biological functions of most circRNAs in PCa are currently unknown. In the present study, we screened 443 differentially expressed circRNAs in metastatic lymph node tissues and PCa tissues from the same patient. From the screening process, we selected a new circRNA, namely circ_0004296, whose host gene plays a key role in various tumors . exoRBase, a web-accessible database (http://www.exoRBase.org), is a repository of circular RNA (circRNA), long non-coding RNA (lncRNA) and messenger RNA (mRNA) derived from RNA-seq data analyses of human blood exosomes, indicating that circulating circRNAs may be diagnostic and prognostic biomarkers . CircSHKBP1 was abundant in serum exosomes and was similarly expressed in the serum and tumors of gastric cancer patients, suggested that circSHKBP1 is a stable biomarker for gastric cancer diagnosis and prognosis . Circ-CCAC1 plays a vital role in cholangiocarcinoma tumorigenesis and metastasis, and exosomal circ-CCAC1 may be an important biomarker/therapeutic target for cholangiocarcinoma . In this study, the downregulation of circ_0004296 in PCa cell lines and clinical samples suggested that circ_0004296 may play a suppressive role in PCa metastasis and serves as a noninvasive liquid biopsy biomarker.
It is known that EMT is important for the progression and metastasis of cancer cells. The primary PCa epithelial cells gain migratory and invasive properties by undergoing EMT . Importantly, our further functional studies revealed that overexpression of circ_0004296 remarkably inhibited proliferation, migration, invasion, and EMT of PC3 and DU145 cells. In line with this result, circ_0004296 knockdown significantly enhanced proliferation, migration, invasion, and EMT of PCa cells. Consistent with these findings, the in vivo experiment demonstrated that overexpression of circ_0004296 markedly inhibited the growth, metastasis, and EMT of PCa cells. Thus, these results strongly supported the suppressive role of circ_0004296 in PCa metastasis.
Although the role of circRNAs in cancers has been widely reported, the downstream regulatory mechanisms remain largely obscure. To date, most circRNAs have been reported to play an important role by acting as ceRNAs to inhibit miRNAs [8, 9]. However, in the present study, circ_0004296 was found to be mainly located in the nucleus of PCa cells, indicating circ_0004296 might not be an miRNA sponge or translated into peptides. At present, only few studies have reported on the mechanism of circRNAs located in the nucleus. Interestingly, circRNAs produced by precursor mRNA back-splicing of their host genes have potential regulatory relationship with their host genes. Circ0005276 plays a crucial role in regulating host gene XIAP expression at the transcriptional level by interacting with the RBP FUS in PCa . Recently, circPOK was found to bind the RNA/DNA-binding proteins ILF2 and ILF3 mainly in the nucleus and potentiate their regulatory activities in mediating both mRNA transcription and stability . It was demonstrated that circSMARCA5 maily located in the nucleus can suppress the expression of host gene SMARCA5 by binding to its parent gene locus to form an R-loop, which results in transcriptional pausing . In our present study, we found that circ_0004296 was mainly located in the nucleus and regulated the host gene ETS1 expression at the post-transcriptional level. Previous studies have shown that the transcription factor ETS1 is involved in tumor progression and metastasis through the transcriptional regulation of EMT-related genes in several cancers [19, 20], including PCa . Recent studies have reported that circRNAs may play crucial roles in regulating gene expression at the transcriptional level by interacting with RBPs [15, 34, 37]. Next, we screened the proteins that specifically bind to circ_0004296 by performing RNA pull-down and mass spectrometry analyses, and bioinformatics analysis was performed to predict RBPs that may bind to circ_0004296. RNA pull-down and western blotting assays confirmed that circ_0004296 interacted with the RBP EIF4A3, a core component of the EJC complex . EIF4A3 has been reported to promote the production of a majority of circRNAs [38,39,40]. On the contrary, a recent study found that EIF4A3 could inhibit the expression of circ_100290 in gastric cancer . In the present study, knockdown of EIF4A3 had no significant impact on circ_0004296 expression in PCa cells. Silencing EIF4A3 remarkably decreased ETS1 mRNA and protein levels. Moreover, RIP assay showed that EIF4A3 interacted with circ_0004296 and ETS1 mRNA. The EJC complex is a critical chromatin remodeler and regulates gene expression [27, 28, 42]. It is known that EIF4A3 plays essential roles in mRNA splicing, nuclear mRNA export, subcellular mRNA localization, translation efficiency, and nonsense-mediated mRNA decay . Previous research showed that circSEMA5A increased the expression of host gene SEMA5A by recruiting EIF4A3 to enhance mRNA stability of SEMA5A, and thereby promoted BC angiogenesis . Additionally, hsa_circ_0030042 promoted the cytoplasmic localization of EIF4A3 and inhibits FOXO1 and beclin1 expression by obstructing the interaction between EIF4A3 and their mRNAs . Emerging studies have demonstrated that through direct interaction, circRNAs can modulate the nuclear localization of proteins as well as their corresponding biological functions. Circ-Amotl1 can promote the nuclear localization of c-Myc, STAT3, PDK1, and AKT1 by direct interaction with them . Through its different regions, circ-Dnmt1 can bind to p53 and Auf1, and facilitate their nuclear localization . Circ-Foxo3 is predominantly localized in the cytoplasm and interacts with ID1, E2F1, HIF1α, and FAK, leading to their cytoplasm retention . Through the interaction with HuR, protein, circAGO2 facilitates its translocation from the nucleus to the cytoplasm . Furthermore, circ_0004296-overexpressing cells were subjected to nuclear and cytoplasmic fractionation. circ_0004296-overexpressing cells showed retention of EIF4A3 protein and ETS1 mRNA in the nucleus as compare to the control group. Thus, our study revealed a novel role of the circRNA circ_0004296 in PCa, wherein it directly interacted with EIF4A3 and suppressed nuclear export of ETS1 mRNA. Although EIF4A3 has been reported to be closely related to tumorigenesis, such as tumorigenesis of colorectal cancer  and pancreatic adenocarcinoma , its role in PCa has not yet been elucidated. Our present study found that EIF4A3 was highly expressed in PCa, and knockdown of EIF4A3 significantly suppressed proliferation, migration, invasion, and EMT of PCa cells. Thus, our results illustrated for the first time that EIF4A3 contributes to PCa development.
The upstream regulatory mechanisms of circRNA downregulation should thus be addressed in future research. The specific motif sites of EIF4A3 binding to circ_0004296 and ETS1 mRNA have not been elucidated and require further studies.
Overall, the present study demonstrated that circ_0004296 was downregulated in PCa and inhibited cancer metastasis by suppressing EMT. We also confirmed circ_0004296 was stably downregulated in the blood plasma and urine of patients with PCa, thus indicating its potential for precision targeted treatment. Mechanistically, EIF4A3 protein mediated circ_0004296-induced inhibition of host gene ETS1 expression at the post-transcriptional level. Furthermore, our research revealed a novel role of the circRNA circ_0004296, wherein it directly interacted with EIF4A3, suppressed nuclear export of ETS1 mRNA, and subsequently suppressed EMT in PCa. Thus, the present study provides novel insights into the potential roles of the circ_0004296/EIF4A3/ETS1 axis in the therapeutic management of PCa.
Availability of data and materials
The mass spectrometry analysis data are included in the manuscript. All other data are available from the corresponding authors upon reasonable request.
Real-time quantitative PCR
Cell Counting Kit-8
In vivo imaging system
RNA fluorescence in situ hybridization
Small interfering RNA
Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136(5):E359–86.
Center MM, Jemal A, Lortet-Tieulent J, Ward E, Ferlay J, Brawley O, et al. International variation in prostate cancer incidence and mortality rates. Eur Urol. 2012;61(6):1079–92.
James ND, Spears MR, Clarke NW, Dearnaley DP, De Bono JS, Gale J, et al. Survival with newly diagnosed metastatic prostate cancer in the “docetaxel era”: data from 917 patients in the control arm of the STAMPEDE trial (MRC PR08, CRUK/06/019). Eur Urol. 2015;67(6):1028–38.
Chowdhury S, Bjartell A, Lumen N, Maroto P, Paiss T, Gomez-Veiga F, et al. Real-world outcomes in first-line treatment of metastatic castration-resistant prostate cancer: the prostate cancer registry. Target Oncol. 2020;15(3):301–15.
Chen LL, Yang L. Regulation of circRNA biogenesis. RNA Biol. 2015;12(4):381–8.
Shang Q, Yang Z, Jia R, Ge S. The novel roles of circRNAs in human cancer. Mol Cancer. 2019;18(1):6.
Zheng X, Chen L, Zhou Y, Wang Q, Zheng Z, Xu B, et al. A novel protein encoded by a circular RNA circPPP1R12A promotes tumor pathogenesis and metastasis of colon cancer via hippo-YAP signaling. Mol Cancer. 2019;18(1):47.
Ding C, Wu Z, You H, Ge H, Zheng S, Lin Y, et al. CircNFIX promotes progression of glioma through regulating miR-378e/RPN2 axis. J Exp Clin Cancer Res. 2019;38(1):506.
Yu J, Xu QG, Wang ZG, Yang Y, Zhang L, Ma JZ, et al. Circular RNA cSMARCA5 inhibits growth and metastasis in hepatocellular carcinoma. J Hepatol. 2018;68(6):1214–27.
Abdelmohsen K, Panda AC, Munk R, Grammatikakis I, Dudekula DB, De S, et al. Identification of HuR target circular RNAs uncovers suppression of PABPN1 translation by CircPABPN1. RNA Biol. 2017;14(3):361–9.
Du WW, Yang W, Chen Y, Wu ZK, Foster FS, Yang Z, et al. Foxo3 circular RNA promotes cardiac senescence by modulating multiple factors associated with stress and senescence responses. Eur Heart J. 2017;38(18):1402–12.
Du WW, Yang W, Liu E, Yang Z, Dhaliwal P, Yang BB. Foxo3 circular RNA retards cell cycle progression via forming ternary complexes with p21 and CDK2. Nucleic Acids Res. 2016;44(6):2846–58.
Du WW, Fang L, Yang W, Wu N, Awan FM, Yang Z, et al. Induction of tumor apoptosis through a circular RNA enhancing Foxo3 activity. Cell Death Differ. 2017;24(2):357–70.
Yang Z, Qu CB, Zhang Y, Zhang WF, Wang DD, Gao CC, et al. Dysregulation of p53-RBM25-mediated circAMOTL1L biogenesis contributes to prostate cancer progression through the circAMOTL1L-miR-193a-5p-Pcdha pathway. Oncogene. 2019;38(14):2516–32.
Feng Y, Yang Y, Zhao X, Fan Y, Zhou L, Rong J, et al. Circular RNA circ0005276 promotes the proliferation and migration of prostate cancer cells by interacting with FUS to transcriptionally activate XIAP. Cell Death Dis. 2019;10(11):792.
Mak P, Leav I, Pursell B, Bae D, Yang X, Taglienti CA, et al. ERbeta impedes prostate cancer EMT by destabilizing HIF-1alpha and inhibiting VEGF-mediated snail nuclear localization: implications for Gleason grading. Cancer Cell. 2010;17(4):319–32.
De Craene B, Berx G. Regulatory networks defining EMT during cancer initiation and progression. Nat Rev Cancer. 2013;13(2):97–110.
Chen X, Chen RX, Wei WS, Li YH, Feng ZH, Tan L, et al. PRMT5 circular RNA promotes metastasis of urothelial carcinoma of the bladder through sponging miR-30c to induce epithelial-mesenchymal transition. Clin Cancer Res. 2018;24(24):6319–30.
Yalim-Camci I, Balcik-Ercin P, Cetin M, Odabas G, Tokay N, Sayan AE, et al. ETS1 is coexpressed with ZEB2 and mediates ZEB2-induced epithelial-mesenchymal transition in human tumors. Mol Carcinog. 2019;58(6):1068–81.
Peng C, Li Z, Niu Z, Niu W, Xu Z, Gao H, et al. Norcantharidin suppresses colon cancer cell epithelial-mesenchymal transition by inhibiting the alphavbeta6-ERK-Ets1 signaling pathway. Sci Rep. 2016;6:20500.
Zhang W, Mao S, Shi D, Zhang J, Zhang Z, Guo Y, et al. MicroRNA-153 decreases tryptophan catabolism and inhibits angiogenesis in bladder cancer by targeting Indoleamine 2,3-dioxygenase 1. Front Oncol. 2019;9:619.
Zhang W, Zhang J, Zhang Z, Guo Y, Wu Y, Wang R, et al. Overexpression of Indoleamine 2,3-dioxygenase 1 promotes epithelial-mesenchymal transition by activation of the IL-6/STAT3/PD-L1 pathway in bladder cancer. Transl Oncol. 2019;12(3):485–92.
Dudekula DB, Panda AC, 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(1):34–42.
Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42(Database issue):D92–7.
Meng X, Hu D, Zhang P, Chen Q, Chen M. CircFunBase: a database for functional circular RNAs. Database (Oxford). 2019;2019:baz003.
Du WW, Zhang C, Yang W, Yong T, Awan FM, Yang BB. Identifying and characterizing circRNA-protein interaction. Theranostics. 2017;7(17):4183–91.
Meyers C. Hmong children and their families: consideration of cultural influences in assessment. Am J Occup Ther. 1992;46(8):737–44.
Chan CC, Dostie J, Diem MD, Feng W, Mann M, Rappsilber J, et al. eIF4A3 is a novel component of the exon junction complex. RNA. 2004;10(2):200–9.
Chen Y, Yang F, Fang E, Xiao W, Mei H, Li H, et al. Circular RNA circAGO2 drives cancer progression through facilitating HuR-repressed functions of AGO2-miRNA complexes. Cell Death Differ. 2019;26(7):1346–64.
Dittmer J. The role of the transcription factor Ets1 in carcinoma. Semin Cancer Biol. 2015;35:20–38.
Li S, Li Y, Chen B, Zhao J, Yu S, Tang Y, et al. exoRBase: a database of circRNA, lncRNA and mRNA in human blood exosomes. Nucleic Acids Res. 2018;46(D1):D106–D12.
Xie M, Yu T, Jing X, Ma L, Fan Y, Yang F, et al. Exosomal circSHKBP1 promotes gastric cancer progression via regulating the miR-582-3p/HUR/VEGF axis and suppressing HSP90 degradation. Mol Cancer. 2020;19(1):112.
Xu Y, Leng K, Yao Y, Kang P, Liao G, Han Y, et al. A circular RNA, cholangiocarcinoma-associated circular RNA 1, contributes to cholangiocarcinoma progression, induces angiogenesis, and disrupts vascular endothelial barriers. Hepatology. 2021;73(4):1419–35.
Guarnerio J, Zhang Y, Cheloni G, Panella R, Mae Katon J, Simpson M, et al. Intragenic antagonistic roles of protein and circRNA in tumorigenesis. Cell Res. 2019;29(8):628–40.
Xu X, Zhang J, Tian Y, Gao Y, Dong X, Chen W, et al. CircRNA inhibits DNA damage repair by interacting with host gene. Mol Cancer. 2020;19(1):128.
Rodgers JJ, McClure R, Epis MR, Cohen RJ, Leedman PJ, Harvey JM, et al. ETS1 induces transforming growth factor beta signaling and promotes epithelial-to-mesenchymal transition in prostate cancer cells. J Cell Biochem. 2019;120(1):848–60.
Li Z, Huang C, Bao C, Chen L, Lin M, Wang X, et al. Exon-intron circular RNAs regulate transcription in the nucleus. Nat Struct Mol Biol. 2015;22(3):256–64.
Jiang Z, Tai Q, Xie X, Hou Z, Liu W, Yu Z, et al. EIF4A3-induced circ_0084615 contributes to the progression of colorectal cancer via miR-599/ONECUT2 pathway. J Exp Clin Cancer Res. 2021;40(1):227.
Meng L, Zheng Y, Liu S, Ju Y, Ren S, Sang Y, et al. ZEB1 represses biogenesis of circ-DOCK5 to facilitate metastasis in esophageal squamous cell carcinoma via a positive feedback loop with TGF-beta. Cancer Lett. 2021;519:117–29.
Zhang C, Han X, Yang L, Fu J, Sun C, Huang S, et al. Circular RNA circPPM1F modulates M1 macrophage activation and pancreatic islet inflammation in type 1 diabetes mellitus. Theranostics. 2020;10(24):10908–24.
Wang G, Sun D, Li W, Xin Y. CircRNA_100290 promotes GC cell proliferation and invasion via the miR-29b-3p/ITGA11 axis and is regulated by EIF4A3. Cancer Cell Int. 2021;21(1):324.
Giorgi C, Yeo GW, Stone ME, Katz DB, Burge C, Turrigiano G, et al. The EJC factor eIF4AIII modulates synaptic strength and neuronal protein expression. Cell. 2007;130(1):179–91.
Wang L, Li H, Qiao Q, Ge Y, Ma L, Wang Q. Circular RNA circSEMA5A promotes bladder cancer progression by upregulating ENO1 and SEMA5A expression. Aging (Albany NY). 2020;12(21):21674–86.
Yu F, Zhang Y, Wang Z, Gong W, Zhang C. Hsa_circ_0030042 regulates abnormal autophagy and protects atherosclerotic plaque stability by targeting eIF4A3. Theranostics. 2021;11(11):5404–17.
Yang Q, Du WW, Wu N, Yang W, Awan FM, Fang L, et al. A circular RNA promotes tumorigenesis by inducing c-myc nuclear translocation. Cell Death Differ. 2017;24(9):1609–20.
Du WW, Yang W, Li X, Awan FM, Yang Z, Fang L, et al. A circular RNA circ-DNMT1 enhances breast cancer progression by activating autophagy. Oncogene. 2018;37(44):5829–42.
Wu H, Qin W, Lu S, Wang X, Zhang J, Sun T, et al. Long noncoding RNA ZFAS1 promoting small nucleolar RNA-mediated 2′-O-methylation via NOP58 recruitment in colorectal cancer. Mol Cancer. 2020;19(1):95.
Li Q, Lei C, Lu C, Wang J, Gao M, Gao W. LINC01232 exerts oncogenic activities in pancreatic adenocarcinoma via regulation of TM9SF2. Cell Death Dis. 2019;10(10):698.
We thank Yan-yan He and Huijun Fu (pathologists at the Department of Pathology, Shanghai Tenth People’s Hospital, Tongji University, China) for their pathology expertise and technical assistance. We thank Meiyu Bao, Shaobo Xue, and Sujing Qiang (Central Laboratory of Shanghai Tenth People’s Hospital) for assistance in experiments. We thank International Science Editing (http://www.internationalscienceediting.com) for editing this manuscript.
This study was funded by the Shanghai Science Committee Foundation (grant number 19411967700) and Natural Science Foundation of China (grant numbers 81472389, 81802554, and 81672549).
Ethics approval and consent to participate
The present study was approved by the ethical committee of Shanghai Tenth People’s Hospital of Tongji University.
Consent for publication
Informed consent was obtained from all participants included in the study.
The authors declare that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The Sanger sequencing of the linear circ_0004296 RNA sequence carried on the plasmid.
Primers, siRNAs, probes, and plasmid used in the present study.
The specific RBPs interacted with circ_0004296.
A. The expression abundance of six screened circRNA candidates. B. qRT-PCR showed the expression of circ_0001708 (circRNA_19382) in PCa cell lines. C. qRT-PCR showed the expression of circ_0002842 (circRNA_22335) in PCa cell lines. D. qRT-PCR showed the expression of circ_0004390 (circRNA_00938) in PCa cell lines. E. Northern blotting assay confirmed the integrity of the linear circ_0004296 RNA sequence carried on the plasmid.
The circ_0004296 knockdown promoted the proliferation, migration, invasion, and EMT induction of PCa cells. A. Schematic illustration of designed siRNAs (siRNA#1, siRNA#2, and siRNA#3) targeting the back-splice junction of circ_0004296. B. The expression of circ_0004296 and ETS1 was detected by qRT-PCR in PC3 and DU145 cells transfected with siRNA-circ_0004296 or si-NC. C, D. CCK8 and colony formation assays were performed to measure the viability of PCa cells after circ_0004296 knockdown. E-H. The migration and invasion abilities of PCa cells were analyzed by wound healing and Transwell assays for migration and invasion after circ_0004296 knockdown. Scale bars =100um. I. Western blotting assay was performed to detected the expression of EMT-related proteins, namely E-cadherin, N-cadherin, Vimentin, and Snail, after circ_0004296 knockdown. J. Protein profiling was used to screen specific proteins bound by circ_0004296 in PCa cells. *P < 0.05, **P < 0.01, ***P < 0.001.
A. Knockdown of EIF4A3 had no effect on circ_0004296 expression. B. The efficiency of ETS1 knockdown at protein level. C. The efficiency of ETS1 overexpression at protein level. D. Colony formation assays were used to measure viability of PCa cells after ETS1 knockdown. E. The migration abilities of PCa cells were analyzed by Transwell assays for migration after ETS1 knockdown. Scale bars =100um. F. Colony formation assays were used to measure viability of PCa cells after ETS1 overexpression. G. The migration abilities of PCa cells were analyzed by Transwell assays for migration after ETS1 overexpression. Scale bars =100um. H. Kaplan-Meier analysis for different ETS1 protein expression in PCa patients. I, J. Luciferase-labeled stably DU145-overexpressing circ_0004296 and luciferase-labeled DU145-vector cells were constructed. Scale bars =200um. *P < 0.05, **P < 0.01, ***P < 0.001.
About this article
Cite this article
Mao, S., Zhang, W., Yang, F. et al. Hsa_circ_0004296 inhibits metastasis of prostate cancer by interacting with EIF4A3 to prevent nuclear export of ETS1 mRNA. J Exp Clin Cancer Res 40, 336 (2021). https://doi.org/10.1186/s13046-021-02138-8
- Prostate cancer