- Open Access
Hsa_circ_0003998 promotes epithelial to mesenchymal transition of hepatocellular carcinoma by sponging miR-143-3p and PCBP1
Journal of Experimental & Clinical Cancer Research volume 39, Article number: 114 (2020)
Circular RNAs (circRNAs) play a critical regulatory role in cancer progression. However, the underlying mechanisms of circRNAs in hepatocellular carcinoma (HCC) metastasis remain mostly unknown.
Has_circ_0003998 (circ0003998) was identified by RNAs sequencing in HCC patients with /without portal vein tumor thrombus (PVTT) metastasis. The expression level of circ0003998 was further detected by in situ hybridization on tissues microarray (ISH-TMA) and qRT-PCR in 25 HCC patients with PVTT metastasis. Moreover, the 25 HCC patients with PVTT metastasis and 50 HCC patients without PVTT metastasis were recruited together to analyze the correlation between circ0003998 expression and HCC clinical characteristics. Transwell, migration and CCK8 assays, as well as nude mice model of lung or liver metastasis were used to evaluate the role of circ0003998 in epithelial to mesenchymal transition (EMT) in HCC. The regulatory mechanisms of circ0003998 in miR-143-3p and PCBP1 were determined by dual-luciferase reporter assay, nuclear-cytoplasmic fractionation, fluorescent in situ hybridization, RNA pull- down, microRNA sequence, western blot and RNA immunoprecipitation.
Compared with adjacent normal liver tissues (ANL), circ0003998 expression was significantly upregulated in PVTT tissues and HCC tissues, and its expression correlates with the aggressive characteristics of HCC patients. Further assays suggested that circ0003998 promoted EMT of HCC both in vitro and in vivo. Mechanistically, our data indicated that circ0003998 may act as a ceRNA (competing endogenous RNA) of microRNA-143-3p to relieve the repressive effect on EMT-related stimulator, FOSL2; meanwhile, circ0003998 could bind with PCBP1-poly(rC) binding protein 1 (PCBP1) to increase the expression level of EMT-related genes, CD44v6.
Circ0003998 promotes EMT of HCC by circ0003998/miR-143-3p/FOSL2 axis and circ0003998 /PCBP1/CD44v6 axis.
Hepatocellular carcinoma (HCC) is one of the most common and lethal tumors worldwide due to metastases and recurrence [1, 2]. Epithelial to mesenchymal transition (EMT) is a complex biological process that plays key role in tumor metastases, and it disrupts the intercellular junctions, polarity, order, and consistency of the cells, leading to tumor recurrence . More and more studies have elucidated the indispensable role of EMT in metastatic dissemination of HCC; and it highlighted the need for further study of EMT-related molecular mechanisms in HCC [3, 4].
Covalently closed circular RNAs (circRNAs) are connected by the back-splicing of exons (3′ and 5′ ends) or introns of precursor mRNAs . Recent studies have shown the critical functions of circRNAs, such as microRNAs (miRNAs) or RNA-binding proteins (RBPs) sponging, proteins translation, transcription, and splicing modulation [6, 7]. As reported, circRNAs perform better as diagnostic and therapeutic targets than linear transcripts by virtue of its structural stability, species conservation, and cell/tissue-specificity [8, 9]. More importantly, abnormal expression of non-coding RNAs is closely related to EMT in the metastatic process of cancers [10, 11].
In the present study, we investigated the expression profile of circRNAs in HCC patients with or without portal vein tumor thrombus (PVTT) metastasis using RNA-sequencing (RNA-seq). And we characterized a novel circRNAs, circ0003998 (circBase ID: hsa_circ_0003998) that derived from ADP ribosylation factor guanine nucleotide exchange factor 2 (ARFGEF2) and located at chr20:47570092–47,580,435. We further showed that circ0003998 plays a key role in HCC metastasis through circ0003998/miR-143-3p/FOSL2 axis and circ0003998 /PCBP1/CD44v6 axis.
Materials and methods
Patient tissues and cell lines
The study recruited 25 HCC patients with PVTT metastasis (cohort 1) and 50 HCC patients without PVTT metastasis (cohort 2) from the Eastern Hepatobiliary Surgery Hospital (Shanghai, China). The patients with a history of preoperative chemoradiotherapy were excluded. HCC tissues, PVTT tissues, and corresponding adjacent normal liver tissues (ANL) tissues were collected from cohort 1; and another HCC tissues and corresponding ANL tissues were collected from cohort 2. The tissues were verified by two pathologists independently. The study complied with the Declaration of Helsinki and was approved by the Human Ethics Committee of the Eastern Hepatobiliary Surgery Hospital (Shanghai, China). Moreover, all patients in this study provided written informed consent. HCC cell lines (HepG2, HuH7, Hep3B, and PLC/PRF/5) were purchased from the American Type Culture Collection (ATCC, USA). The MHCC97H cells, hepatocyte cell line L02, and HEK293T cells were obtained from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). No cell lines were contaminated by other cells, such as Hela, as shown by short tandem repeat (STR) data (Biowing Applied Biotechnology Co., Ltd., Shanghai, China) (Additional file 1). All cells were maintained in humidified incubators (37 °C) with 5% CO2 and were cultured in Dulbecco’s modified Eagle’s medium (SIGMA, USA) with 10% fetal bovine serum (GIBCO, BRL).
Total RNA isolation, RNA-seq, and RNase R treatment
Total RNA were isolated from the cultured cells and fresh tissues by Trizol reagent (Invitrogen, USA). Total RNA was extracted by RNAprep pure Tissue Kit (TIANGEN), and rRNA was depleted by GeneRead rRNA Depletion Kit (QIAGEN). The mRNA was purified from total RNA by oligo (dT) magnetic beads, and fragmented into 200–500 bp; the cleaved RNA fragments were reverse-transcribed into cDNA, and enriched by PCR to create the final complementary DNA (cDNA) libraries. The harvested target bands were quantified by Agilent 2100 and then subjected to deep sequencing with the Illumina HiSeq 2000. HCC cells were mixed with RNase R (3 U/μg, Epicentre, Madison, WI) at 37 °C for 15 min, and then qRT-PCR was used to assess the expression stability of circ0003998 as compared to ARFGEF2 mRNA.
Reverse transcription reaction and quantitative real-time PCR
The first-strand cDNA was synthesized by the PrimeScript RT reagent kit (Takara Bio Inc., China), and the reverse transcription of miRNAs was generated by miRNAs First Strand cDNA Synthesis Poly A Tailing Kit (Sangon Bio Inc., China). SYBR Premix Ex TaqII (Takara Bio Inc., China) was used in Quantitative real-time polymerase chain reaction (qRT-PCR). β-actin and U6 were used as endogenous references for mRNA and miRNAs respectively. The comparative 2-ΔΔCt method was used to calculate the relative fold-change of target expression. The sequences of the primers in this study were showed in Additional file 2. All primers were designed and purchased from Sangon (Shanghai, China).
Tissue microarray (TMA) and in situ hybridization (ISH)
After dewaxed in xylene and rehydrated with gradient alcohol, the TMAs were digested by protease K (20 μg/mL) at 37 °C. 3% methanol-H2O2 was used to block endogenous peroxidase. Next, TMAs were hybridized with specific digoxin-labeled circ0003998 probe (8 ng/μL, Digoxin- 5′-Digoxin-GGCCTCCTGCAACTTTAATGGCAGATG-Digoxin-3′, Service Biotechnology, Wu han, China) at 37 °C overnight. TMAs were incubated in BSA serum at room temperature for 30 min and then incubated with mouse anti-digoxin labeled peroxidase (anti-DIG-HRP) at 37 °C for 40 min. The TMAs were stained with freshly prepared diaminoaniline (DAB) solution, and the nucleus was stained with Harris hematoxylin stain.
Fluorescence in situ hybridization (FISH)
FISH assay was used to find the intracellular location of circ0003998. Probes of FISH assay for circ0003998, human U6, and human18S were synthesized by RiboBio Co., Ltd. (Guangzhou, China). Briefly, cells were rinsed with phosphate buffer saline (PBS) and fixed in 4% formaldehyde solution for 10 min at room temperature, and then incubated with 0.5% Triton X-100 solution for 5 min at 4 °C. After pre-hybridization for 30 min, the cells were hybridized with fluorescence-hybridization probes overnight in the dark at 37 °C. Then, cells were washed three times in 4X/2X/1X SSC solution in the dark at 42 °C for 5 min respectively, and laser scanning confocal microscopy was used to visualize the images at 400 × magnification.
Cytoplasmic and nuclear RNA isolation was detected by the PARIS™ Kit (Invitrogen, USA). Briefly, after washed in PBS, the cultured cells were resuspended in cold cell fraction buffer and then incubated on ice for 10 min. The cells samples were centrifuged at 4 °C to separate the supernatant containing cytoplasmic fractions from the pellet containing nuclear fraction. The nuclear pellet was rewashed by the cell fraction buffer and was lysed by the disruption buffer. Both the cytoplasmic fraction and the nuclear fraction were respectively mixed with 2X lysis/ binding solution, followed by 100% ethanol. After centrifugation and washing, the RNA of nuclear and cytoplasmic were respectively obtained with the eluting solution.
Establishment of stable cell lines over expressing hsa_circ_0003998
The circ0003998 expressing lentivirus vector (HBLV-circ0003998-GFP-PURO, termed as OE-circ-vector) and its control vector (HBLV-GFP-PURO, termed as NC-circ-vector) were supplied by Hanbio (Shanghai, China). HCC cells were infected with the lentivirus vectorfoloowed by selection with 2 μg/mL puromycin for 1 week.
Construction of siRNA, plasmids, miRNA mimic, and transient transfection
The small interfering RNAs (siRNAs) of circ0003998 (si-circ0003998) and the corresponding negative control (siRNA-NC) were synthesized by Hanbio (Shanghai, China). The PCBP1 expressing plasmid (OE-PCBP1), FOSL2 expressing plasmid (OE-FOSL2) and the corresponding negative plasmid vector (vector) were provided by Hanbio (Shanghai, China). The siRNAs of PCBP1 (si-PCBP1), siRNAs of FOSL2 (si-FOSL2) were synthesized by GenePharma (Shanghai, China). The microRNA-143-3p mimics (miR-143-3p) and corresponding negative control (miR-143-3p-NC) were synthesized by GenePharma (Suzhou, China). After digestion with Xho I/Xba I, circ0003998 cDNA (full-length: 304 bp) from HepG2 was amplified by qRT-PCR which cloned into pmirGLO Dual-Luciferase miRNA Target Expression Vector (Promega Corporation, USA) to synthesize wild type-plasmid for the luciferase reporter assay (Luci-circ0003998-WT). The corresponding mutant type-plasmid (Luci-circ000 3998-MT) contained with mutated miRNA binding sites and was synthesized by Genscript (Nanjing, China). Lipofectamine 3000 kit (Invitrogen) was used in transient transfection according to the manufacturer’s instructions.
The cells were seeded in six-well plates and the center of each well was straightly scratched using a 10-μL plastic pipette tip. After floating cells were gently removed with PBS, cells in the six well plates continued to be cultured for 48 h. The wound-healing process were monitored by inverted light microscope.
Cell counting Kit-8 assay
Cell proliferation was assessed by Cell Counting Kit-8 kit (CCK-8, Dojindo Chemical Laboratory, Kumamoto, Japan). HCC cells (3 × 103/well) were seeded in 96-well plates with 6 replicates. 10 μL/well of CCK8 solution and 90 μL/well of fresh medium were mixed and added to each well at 0, 24, 48, 72, and 96 h, respectively, After incubation for 1.5 h in cell incubator, the absorbance of the medium at 450 nm was measured.
Migration assay and invasion assay
24-well transwell migration champers (8 μm size, Corning-3422, USA) and pre-biocoated matrigel transwell invasion chambers (8 μm size, Corning-354,480, USA) were used for migration assay and invasion assay respectively. Briefly, a total of 5 × 104 HuH7, 8 × 104 MHCC97H, and 8 × 104 HepG2 cells were resuspended in 200 μL serum-free DMEM medium and were seeded into the inner chambers, respectively. 600 μL DMEM medium containing 10% FBS as the attractant was loaded to the bottom chambers. For migration assay, the cells were incubated for 24 h, and for invasion assay, the cells were incubated for 48 h. After incubation, the cells that migrated or invaded through the pores were fixed with 4% paraformaldehyde, and stained with 0.1% crystal violet, and then counted at least in five random fields under 200 × microscope.
Luciferase reporter assay
HCC cells or HEK293T cells (3 × 105) were seeded in a 24-well plate overnight. 250 ng of Luci- 0003998-WT/Luci-circ0003998-MT were co-transfected into the cells with 100 nM miR-143-3p mimic/miR-143-3p mimic-NC using Lipofectamine 3000. In 48 h after transfection, firefly and Renilla luciferase fluorescence were detected by the Dual-Luciferase Reporter Assay kit (Promega, USA). The relative luciferase activities were quantified by Firefly/Renilla fluorescence.
RNA immunoprecipitation (RIP) assay
Magna RIP RNA Binding Protein Immunoprecipitation Kit (Millipore, Billerica, MA, USA) was used in RIP assay according to the manufacturer’s protocol. HepG2-OE cells (2 × 107) were lysed in RIP lysis buffer. 5 μg of human anti-targeted protein antibody or negative anti-IgG antibodies was incubated with magnetic beads for 2 h. Then, 100 μL RIP lysate was incubated with bead-antibody complex in 900 μL RIP Immunoprecipitation buffer (shaking at 4 °C overnight). After the beads were incubated with proteinase K buffer for 30 min at 55 °C, the immunoprecipitated RNAs were finally extracted to further detect the expression levels of circ0003998 by qRT-PCR.
HepG2 cells (2 × 107) were washed in cold PBS and then lysed in a lysis buffer. Biotin-labeled circ0003998 probe (5′-Biotin-aaaTTCCAGTTCTCTGGCCTCCTGCAACTTTAATGGCAGATG TGACTACAT-3′) and control probe (5′-Biotin-aaaATGTAGTCACATCTGCCATTAAAGTTGC AGGAGGCCAGAGAACTGGAA-3′) were synthesized by CloudSeq Biotech Inc. (Shanghai, China). Biotin-coupled probes were bounded on magnetic beads and then incubated with lysates of HepG2-OE cells overnight at 4 °C. After purification, the pull-down product was analyzed by miRNA sequencing and protein spectrum.
In vivo tumorigenesis and metastasis assays
The animal experiments met the demands of laboratory animal welfare and ethics, and approved by the Institutional Animal Care and Use Committee of Naval Medical University (Shanghai, China). The nude male BALB/c mice (four weeks old) were purchased from Animal Center of Naval Medical University and housed under pathogen-free conditions with standard pellet diet and water. Twenty mice were randomly divided into four groups for the construction of two model. MHCC97H-OE cells (3 × 106) or MHCC97H-NC cells (3 × 106) were injected into the tail vein to establish the lung metastatic model. MHCC97H-OE cells (5 × 106) or MHCC97H-NC cells (5 × 106) were injected into the spleen to establish liver metastatic model. IVIS@ LuminaII system was used to measure the fluorescence value of GFP at excitation wave-length of 488 nm, which monitor the metastases of two types of mice models. Mice were sacrificed with CO2; their lungs and livers were stained with H&E and re-examined microscopically for the development of metastatic foci.
HCC cells were lysed in a RIPA buffer (Beyotime Biotechnology, Nantong, China). BCA Protein Assay kit (Beyotime Biotechnology, Nantong, China) was used to detect protein concentration. The cell lysates were separated by polyacrylamide gel electrophoresis (4–20% SurePAGE Gel, Genscript, Nanjing, China), and then transferred onto to polyvinylidene difluoride (PVDF) membrane (Bio-Rad, CA, USA). After blocking with 5% nonfat milk in Tris-buffered saline-Tween (TBST) buffer at room temperature for 90 min, the membranes were incubated with the primary antibody at 4 °C overnight. After wash with TBST buffer (10 min × 3), the membranes were incubated with the secondary antibody at room temperature for 1 h. The membranes were examined using the enhanced chemiluminescence kit (Beyotime Biotechnology, Nantong, China). Anti-E-cadherin, anti-Vimentin, anti-N-cadherin, anti-Slug and anti-Snail antibodies were purchased from Cell Signaling Technologies (CST#9782; Boston, MA, USA). Anti-β-actin (ab227387), anti-CD44v6 (ab78960) and anti-FOSL2 (ab124830) antibodies were purchased from Abcam (Cambridge, UK).
SPSS 19.0 software and Prism version 7.0 were used for statistical analyses. Data from three independent experiments were shown as mean ± standard deviation (SD), followed by Student’s t-test, one-way analysis of variance (ANOVA) and chi-square tests. The correlations were analyzed by Pearson correlation. The result with *p-value< 0.05, **p-value < 0.01, or ***p-value < 0.001 was considered as statistically significant.
Confirmation of circ0003998 characteristics
We run the RNA-sequencing of HCC samples from patients with PVTT metastasis (PVTT group, n = 3) and without PVTT metastasis (HCC group, n = 3). Based on fold change value > 2.0, we identified 13 upregulated, and 8 down regulated circRNAs in the PVTT group as depicted by the heat map (Fig. 1a). Circ0003998 derived from exons (6 to 7) of ARFGEF2 gene were significantly upregulated in the PVTT group (p = 0.000000273) than in the HCC group; and therefore, circ0003998 was selected as the candidate circRNA.
To further identify the existence of circ00039998, we amplified circ0003998 by the designed divergent primers and inserted the corresponding PCR products into the T vector for Sanger sequencing. As observed, the corresponding result was consistent with the back-spliced region of circ0003998 provided in circBASE (Fig. 1b). In addition, we designed the convergent primers for the amplification of linear mRNA from exon 6 of ARFGEF2, in order to rule out other potential factors that could induce such head-to-tail splicing [10, 12]. And the results showed, in HCC cells, as opposed to that the linear transcripts could be produced from both cDNA and genomic DNA (gDNA) by convergent primers, circ0003998 was only detected from cDNA instead of gDNA by divergent primers (Fig. 1c). Subsequently, we confirmed the circular characteristics of circ0003998, as the linear transcripts of ARFGEF2 in HepG2 cells could be degraded by the RNase R treatment while circ0003998 could resist to the RNase R treatment (Fig. 1d).
We next detected the expression level of circ0003998 in 25 HCC patients with PVTT metastasis by both qRT-PCR and TMAs-ISH. The results indicated that circ0003998 was significantly upregulated in HCC and PVTT tissues as compared to ANL tissues of HCC patients (cohort 1) (Fig. 1e-g). We further identified that such high expression levels of circ0003998 was correlated with advanced TNM stage (P = 0.029) and high serum levels of AFP (P = 0.016) in 75 HCC patients (25 samples from cohort 1 and 50 samples from cohort 2) (Table 1). Moreover, the expression level of circ0003998 was also higher in HCC cells than in L02 cells (Fig. 1h).
Circ0003998 promotes the EMT of HCC in vitro and in vivo
To explore the biological function of circ0003998 in HCC, we detected the expression profiles of mRNA in 3 paired HCC-OE cells and HCC-NC cells by transcription sequencing (mRNA-seq). Compared with HCC-NC cells, 1027 mRNAs were differentially expressed in HCC-OE cells, among which, expression of 431 mRNAs was upregulated and expression of 596 mRNAs was down regulated, as shown by the volcano plot (Fig. 1i). Gene Ontology (GO) analysis indicated that these mRNAs are enriched in EMT pathways, such as migration, invasion, and focal adhesion (Fig. 1j). Hence, we focused our attention on the role of circ0003998 in EMT of HCC.
We infected MHCC97H and HepG2 cells (MHCC97H-OE/NC and HepG2-OE/NC) with OE/NC -circ-vector and transfected si/NC-circ0003998 into HuH7 and HepG2 cells (HuH7-si/NC and HepG2-si/NC) (Fig. S1a). Firstly, we found that over expression of circ0003998 inhibits the protein level of E-cadherin and enhances the protein level of Vimentin, N-cadherin, Snail and Slug; while down expression of circ0003998 showed the opposite results (Fig. 1k, Fig. S1b). Secondly, we explored the role of circ0003998 in EMT of HCC by in vitro cell phenotypic assays. Transwell assay indicated that over expression of circ0003998 and knockdown of circ0003998 expression respectively increases and decreases cell metastasis (Fig. 2a-d). Wound healing assay showed that the mobility of HCC cells is significantly enhanced when circ003998 was over expressed, but markedly impaired when circ0003998 expression was knocked down (Fig. 2e-h). Furthermore, the results of CCK8 assay demonstrated that ectopic expression of circ0003998 promotes the proliferation of HCC cells, whereas the down regulation of circ0003998 restrained the proliferation of HCC cells (Fig. 2i-j). Thirdly, we applied nude mice model with lung metastatic or with liver metastatic to study the role of circ0003998 in EMT of HCC in vivo. By monitoring the fluorescence value of GFP, we found that circ0003998 significantly accelerated the metastasis of HCC into both lung and liver. (Fig. 2k-m). Furthermore, haematoxylin and eosin (H&E) staining suggested that over expression of circ0003998 remarkably promoted metastasis of HCC into both lung and liver (Fig. 2n). And the number of metastatic foci was also increased with the over expression of circ0003998 in liver metastatic mice model (Fig. 2o-p). These results indicated that circ0003998 could play an oncogenic role in the EMT in HCC both in vitro and in vivo.
Circ0003998 function as a sponge for miR-143-3p
Given that circRNAs could function as sponges for miRNAs and the latter was mostly located in the cytoplasm, we firstly tried to explore the subcellular localization of circ0003998 in the HCC cells. By FISH assay, we found circ0003998 in both cytoplasm and nucleus (Fig. 3a), which was further verified by quantitatively detecting the expression of circ0003998 in the nuclear and cytoplasmic fractions in HCC cell lines (Fig. 3b). Considering that miRNAs silence the expression of targeted genes by binding with argonaute 2 (AGO2) , secondly, we run RIP assay in HepG2 cells using AGO2 antibodies. We found that circ0003998 was significantly enriched in the presence of AGO2 antibodies, suggesting that circ0003998 could act as the platform for the binding between miRNAs and AGO2 (Fig. 3c). Taken together, our data demonstrated the potentiality of circ0003998 to sponge miRNAs.
Following, we run the pull-down assay using specific probes against circ0003998 to screen the specific miRNA that possibly sponged by circ0003998. The precipitation complexes obtained from HCC cells were analyzed by miRNAs sequencing, and the result indicated that miR-143-3p is significantly enriched and it might be the potential targets of circ0003998 (Fig. 3d). We further confirmed this result by FISH assay that most circ0003998 (red) and miR-143-3p (green) are co-located in the cytoplasm in HCC cells (Fig. 3e, Fig. S1c). Additionally, dual-luciferase reporter assay in both HEK293T cells and HepG2 cells also indicated the consistent result that miR-143-3p mimic could inhibit the luciferase activity of wide-type circ0003998 vector (circ0003998-WT vector) as opposed to mutant-type circ0003998 vector (circ0003998-Mut vector) (Fig. 3f-h).
Next, we detected the correlation between the expression of circ0003998 and miR-143 -3p. The expression levels of miR-143-3p were down regulated in HCC and PVTT tissues as compared to ANL tissues (cohort 1), respectively (Fig. 3i). Moreover, the expression level of miR-143-3p showed no significant changes after over expression or silencing of circ0003998 in HepG2 cells (Fig. 3j), while the expression level of circ0003998 showed no significant changes after transfection with miR-143-3p mimic or inhibitor (Fig. 3k). Additionally, the functional rescue experiments indicated that miR-143-3p mimics could abolish the invasion-promoting effects of circ0003998 over expression, whereas miR-143-3p inhibitors could reverse the invasion- suppressing effects of silencing circ0003998 (Fig. 3l-o). These findings suggested that circ0003998 directly sponges with miR-143-3p as the ceRNA in HCC, and circ0003998 is neither induced into the degradation of miR-143-3p nor regulated by miR-143-3p.
Circ0003998 promotes the EMT in HCC cells through miR-143-3p/FOSL2 pathway
Given above results, we further screened the targeted genes of miR-143-3p which play the oncogenetic role in the EMT of HCC and might be released due to the sponge between circ0003998 and miR-143-3p. Firstly, we predicted the possible target genes of miR-143-3p in TargetScan (http://www.targetscan.org/vert_72/) and miRanda (http://miranda.org.uk/) (Additional file 3). Secondly, we analyzed the mRNA expression profile in both HCC cells (circ0003998-OE vs. circ0003998-NC) and HCC tissues (with/without PVTT metastasis) by mRNA-seq. According to the prediction of target genes of miR-143-3p and the results of mRNA-seq, we screened 9 differentially expressed genes (NOB1, MSI-2, FOSL2, ITGA6, MAPK7, TRL2, S100PBP, SFXN1, and EFS) (Fig. 4a).
Among the 9 candidate genes, the mRNA expression of FOSL2, ITGA6, MAPK7 were significantly down regulated by miR-143-3p mimics and upregulated by miR-143-3p inhibitors in HCC cells (Fig. 4b). Then the mRNA expression of these three targets (FOSL2, ITGA6, MAPK7) were detected after over expressing or silencing circ0003998 and we found that the mRNA expression of FOSL2 is significantly enhanced in MHCC97H-OE and HepG2-OE cells whereas decreased in HuH7-si and HepG2-si cells (Fig. 4c). Moreover, the protein expression of FOSL2 altered with the same trend as the mRNA levels under the regulation of circ0003998 (Fig. 4d). Meaningfully, FOSL2 has been identified as the promoter of EMT on TGF-β/Smad3 signaling pathway in advanced cancers [14,15,16]. Thus we presumed that circ0003998 promoted the EMT of HCC by protecting FOSL2 from down regulating by miR-143-3p, which was verified by means of the following assays.
Firstly, we predicted the binding sites of miR-143-3p with FOSL2 by TargetScan and miRanda (Fig. 4e). Secondly, we identified that circ0003998 increased the expression of FOSL2 in HepG2 cells, which was reversed by exogenous miR-143-3p mimics; while the silence of circ0003998 suppressed the expression of FOSL2 in HepG2 cells, which was retarded by miR-143-3p inhibitors (Fig. 4f-g). Thirdly, the functional rescue experiments consistently showed that FOSL2 is the downstream target gene of circ0003998 /miR-143- 3p axis to promote migration of HCC (Fig. S1d). Finally, we detected the the correlation of expression between circ0003998 and FOSL2 in HCC tissues (cohort 1). The expression levels of FOSL2 were highly upregulated in HCC tissues and PVTT tissues as compared to ANL tissues by both ISH-TMAs and qRT-PCR analysis (Fig. 4h-j). And by Pearson correlation analysis, the expression levels between FOSL2 and circ0003998 showed a positive association in both HCC tissues and PVTT tissues, respectively (Fig. 4k). Thus, our study suggested that circ0003998 functioned as a ceRNA by adsorption of miR-143-3p to release FOSL2, which promoted EMT of HCC.
Circ0003998 binds to poly(rC) binding protein 1 (PCBP1)
Results from above study revealed that circ0003998 played stimulating role in EMT of HCC by sponging with miR-143-3p. Meanwhile, considering that the interaction with RNA binding proteins (RBPs) is an important part of circRNA function , it remained to reveal whether circ0003998 could bind with RBPs to participate in the EMT of HCC. Firstly we performed pull-down assay and protein spectrum to identified the specific RBP bound by circ0003998 (Fig. 5a). And poly(rC) binding protein 1 (PCBP1) was selected as the binding protein of circ0003998 as it was not only precipitated by circ0003998-specific probes but also received the highest prediction score in catRAPID (Additional file 4). Secondly, RIP experiments in HepG2 cells verified the results of pull down assay as it showed that circ0003998 is considerably enriched by specific antibodies against PCBP1 compared with IgG (Fig. 5b). Thirdly, FISH assay in HepG2 cells also confirmed that circ0003998 (red) and PCBP1 (green) were co-located in the nucleus and cytoplasm, especially in the nuclear membrane (Fig. 5c, Fig. S1c).
Further, we predicted the potential binding sites between circ0003998 and PCBP1. Firstly, we analyzed the nucleotides sequence of circ0003998 by online catRAPID and PRIdictor (http://www.rna-society.org/rnainter/php/PRIdictor.php), which indicated that nucleotides at 220–300 with a high potential to bind with protein (Fig. 5d). Secondly, we predicted that nucleotides at 220–280 of circ0003998 and polypeptide at 31–82 of PCBP1 could be the potential binding sites of each other by online catRAPID (Fig. 5e). Thirdly, we analyzed that the polypeptide at 31–82 of PCBP1 were belonged to its K homology 1 (KH1) and we modeled the 3D structure of KH1 (Fig. 5f). Moreover, the results in functional rescue assays also consistently confirmed circ0003998 could bind with PCBP1, for PCBP1 could reverse the role of circ0003998 on the migration ability in HCC cells. (Fig. 5g-h).
Circ0003998 promotes the EMT in HCC cell through PCBP1/CD44v6 pathway
Since PCBP1 is recognized as a tumor suppressor in previous studies [18,19,20], we hypothesized that function of PCBP1 could be inhibited when bind to circ0003998, which thereby promote EMT of HCC. Firstly, Gene Transcription Regulation Database (GTRD, http://gtrd.biouml.org/) were used to predict the possible target genes of PCBP1 (Additional file 5), which were significantly enriched in the EMT-related pathway by GO analysis (Fig. 5i). Secondly, we measured the expression levels of five target genes of PCBP1 (CD44v6, STAT3, MAPK1, ITGB1, DAB2) based on the results of according to the analysis of mRNA sequencing in HCC and tissues and results of previous studies (Additional file 6, Fig. 4a). Among the five potential target genes, the mRNA expression levels of CD44v6 were significantly upregulated in HepG2-OE and MHCC97H-OE cells, while down regulated in HuH7-si and HepG2-si cells (Fig. 5j-k). In addition, previous study reported that PCBP1 regulated the expression of DAB2 at translation level . However, the protein level of DAB2 showed no significant changes with circ003998 upregulation or silencing (Fig. S1e).
Given above results and the promoting-role of CD44v6 in EMT in cancers [22, 23], we further studied that whether circ0003998 promoted the EMT of HCC by protecting the CD44v6 from down regulation by PCBP1. Firstly, we found that the protein levels of CD44v6 was significantly upregulated in HepG2-OE and MHCC97H-OE cells, while it was down regulated in HuH7-si and HepG2-si cells (Fig. 5l). Next, rescue assay showed that PCBP1 could abolish circ000398 regulation on both the mRNA and proteins expression level of CD44v6 in HCC cells (Fig. 5m-n). Then the mRNA levels of CD44v6 were highly upregulated in 25 pairs of HCC and PVTT tissues as compared to ANL tissues (Fig. 5o). Pearson correlation analysis indicated that the expression level of CD44v6 were positively associated with that of circ0003998 both in HCC and PVTT tissues (Fig. 5p). These findings further demonstrated that circ0003998 inhibited HCC metastasis through circ0003998/PCBP1/CD44v6 axis.
More and more circRNAs are discovered using high-throughput sequencing techniques and bioinformatics analysis. Previous study has showed that circ0003998 could promote cell proliferation and invasion by targeting miR-326 in non-small cell lung cancer . In this study, our data for the first time identified that circ0003998 could promote EMT in HCC by in vitro functional assays and in vivo mice models; and meanwhile circ0003998 was upregulated in HCC tissues and PVTT tissues, which were correlated with advanced TNM stage and high serum levels of AFP. Furthermore, our previous study has confirmed that expression of circ0003998 is also upregulated in plasma in HCC, which is significantly correlated with poor overall survival (OS) in HCC patients; and therefore, circ0003998 could be used as independently prognostic factors for poor OS in HCC .
Mechanically, we identified that circ0003998 could sponge miR-143-3p, the suppressor of proliferation, invasion, and metastasis in multiple cancers [26,27,28] as well as the inhibitor of the EMT-related genes [29,30,31]. Furthermore, there was no significant correlation between the expression of miR-143-3p and circ0003998. Previous reports showed that ceRNAs could regulate the activity of miRNAs; however, whether ceRNAs regulate the expression levels of miRNAs might rely on the cellular context [32, 33]. Our data further showed that circ0003998 could sponge with miR-143-3p to increase the expression of FOSL2 and thereby promoted the EMT in HCC and it has been consistently showed that FOSL2 is the target gene of miR-143-3p in osteosarcoma in previous study .
The tertiary structures of circRNAs result in higher protein adsorbing capacity than those of linear RNA sequences. As such, the circRNA-interacting RBPs serves as an essential molecular action mode in genesis, translation, transcriptional regulation of target genes, and extracellular transport . PCBP1 as one of the the RBPs inhibited the tumor formation and metastasis by translation silencing, mRNA alternative splicing, or transcription of carcinogenic genes . Meaningfully, PCBP1 is reported to participate in EMT pathways in cancer, especially in the TGF-β pathway [37,38,39]. And our data showed that circ0003998 bound to PCBP1 to increase the expression of CD44v6 and thereby promoted the EMT in HCC, which was consistent with former studies that PCBP1 regulated alternative splicing of the CD44v6 and thereby inhibited the invasion in HCC .
Previous study showed that PCBP1 regulated DAB2 in the translation process . However, there were no significant changes in the protein level of DAB2 with circ0003998 regulating in our data. The sponging affinity of circRNAs with RBPs was not the same as that between circRNAs and miRNAs, and circRNAs absorbing RBPs might be more complicated than expected . Thus further research is necessary to understand whether circ0003998 inhibits the translation regulating-role of PCBP1 or whether circ0003998 induces the degradation of PCBP1.
We identified that circ0003998 expression is remarkably upregulated in HCC, PVTT tissues, and which is correlated with advanced TNM stage and high serum levels of AFP. Circ0003998 locates in both cytoplasm and nucleus and it regulates EMT of HCC by both circ0003998/miR-143 -3p/FOSL2 axis and circ0003998/PCBP1/CD44v6 axis (Fig. 6).
Availability of data and materials
The data used and analyzed during the current study are available from the corresponding author on reasonable request.
Portal vein tumor thrombus
PCBP1-poly(rC) binding protein 1
Epithelial to mesenchymal transition
ADP ribosylation factor guanine nucleotide exchange factor 2
Adjacent non-cancerous liver tissues
short tandem repeat
Quantitative polymerase chain reaction
polymerase chain reaction
Fluorescence in situ hybridization
Haematoxylin and eosin
Phosphate buffer saline
Small interfering RNAs
Cell counting kit-8 kit
One-way analysis of variance
FOS like 2, AP-1 transcription factor subunit
Variant isoform 6 of CD44
Over expression of circ0003998
Silencing of circ0003998
Villanueva A. hepatocellular Carcinoma. N Engl J Med. 2019;380:1450–62..
Budhu A, Forgues M, Ye QH, Jia HL, He P, Zanetti KA, Kammula US, Chen Y, Qin LX, Tang ZY, et al. Prediction of venous metastases, recurrence, and prognosis in hepatocellular carcinoma based on a unique immune response signature of the liver microenvironment. Cancer Cell. 2006;10:99–111..
Dongre A, Weinberg RA. New insights into the mechanisms of epithelial-mesenchymal transition and implications for cancer. Nat Rev Mol Cell Biol. 2019;20:69–84.
Piera-Velazquez S, Jimenez SA. Endothelial to Mesenchymal transition: role in physiology and in the pathogenesis of human diseases. Physiol Rev. 2019;99:1281–324.
Li X, Yang L, Chen LL. The biogenesis, functions, and challenges of circular RNAs. Mol Cell. 2018;71:428–42.
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.
Vo JN, Cieslik M, Zhang Y, Shukla S, Xiao L, Zhang Y, Wu YM, Dhanasekaran SM, Engelke CG, Cao X, et al. The landscape of circular RNA in Cancer. Cell. 2019;176:869–81 e13.
Wong CM, Tsang FH, Ng IO. Non-coding RNAs in hepatocellular carcinoma: molecular functions and pathological implications. Nat Rev Gastroenterol Hepatol. 2018;15:137–51.
Fang J, Hong H, Xue X, Zhu X, Jiang L, Qin M, Liang H, Gao L. A novel circular RNA, circFAT1(e2), inhibits gastric cancer progression by targeting miR-548g in the cytoplasm and interacting with YBX1 in the nucleus. Cancer Lett. 2019;442:222–32.
Zheng F, Wang M, Li Y, Huang C, Tao D, Xie F, Zhang H, Sun J, Zhang C, Gu C, et al. CircNR3C1 inhibits proliferation of bladder cancer cells by sponging miR-27a-3p and downregulating cyclin D1 expression. Cancer Lett. 2019;460:139–51.
Wei S, Zheng Y, Jiang Y, Li X, Geng J, Shen Y, Li Q, Wang X, Zhao C, Chen Y, et al. The circRNA circPTPRA suppresses epithelial-mesenchymal transitioning and metastasis of NSCLC cells by sponging miR-96-5p. EBioMedicine. 2019;44:182–93.
Jeck WR, Sharpless NE. Detecting and characterizing circular RNAs. Nat Biotechnol. 2014;32:453–61.
Yu J, Xu QG, Wang ZG, Yang Y, Zhang L, Ma JZ, Sun SH, Yang F, Zhou WP. Circular RNA cSMARCA5 inhibits growth and metastasis in hepatocellular carcinoma. J Hepatol. 2018;68:1214–27.
Li Z, Niu H, Qin Q, Yang S, Wang Q, Yu C, Wei Z, Jin Z, Wang X, Yang A, et al. lncRNA UCA1 mediates resistance to Cisplatin by regulating the miR-143/FOSL2-signaling pathway in ovarian Cancer. Mol Ther Nucleic Acids. 2019;17:92–101.
Sun L, Guo Z, Sun J, Li J, Dong Z, Zhang Y, Chen J, Kan Q, Yu Z. MiR-133a acts as an anti-oncogene in hepatocellular carcinoma by inhibiting FOSL2 through TGF-beta/Smad3 signaling pathway. Biomed Pharmacother. 2018;107:168–76.
Wang J, Sun D, Wang Y, Ren F, Pang S, Wang D, Xu S. FOSL2 positively regulates TGF-beta1 signalling in non-small cell lung cancer. PLoS One. 2014;9:e112150.
Zang J, Lu D, Xu A. The interaction of circRNAs and RNA binding proteins: An important part of circRNA maintenance and function. J Neurosci Res. 2020;98(1).
Zhang Y, Meng L, Xiao L, Liu R, Li Z, Wang YL. The RNA-binding protein PCBP1 functions as a tumor suppressor in prostate Cancer by inhibiting mitogen activated protein kinase 1. Cell Physiol Biochem. 2018;48:1747–54.
Jiang P, Li Z, Tian F, Li X, Yang J. Fyn/heterogeneous nuclear ribonucleoprotein E1 signaling regulates pancreatic cancer metastasis by affecting the alternative splicing of integrin beta1. Int J Oncol. 2017;51:169–83.
Wang X, Guo J, Che X, Jia R. PCBP1 inhibits the expression of oncogenic STAT3 isoform by targeting alternative splicing of STAT3 exon 23. Int J Biol Sci. 2019;15:1177–86.
Chaudhury A, Hussey GS, Ray PS, Jin G, Fox PL, Howe PH. TGF-beta-mediated phosphorylation of hnRNP E1 induces EMT via transcript-selective translational induction of Dab2 and ILEI. Nat Cell Biol. 2010;12:286–93.
Branco da Cunha C, Klumpers DD, Koshy ST, Weaver JC, Chaudhuri O, Seruca R, Carneiro F, Granja PL, Mooney DJ. CD44 alternative splicing in gastric cancer cells is regulated by culture dimensionality and matrix stiffness. Biomaterials. 2016;98:152–62.
Hernandez JR, Kim JJ, Verdone JE, Liu X, Torga G, Pienta KJ, Mooney SM. Alternative CD44 splicing identifies epithelial prostate cancer cells from the mesenchymal counterparts. Med Oncol. 2015;32:159.
Yu W, Jiang H, Zhang H, Li J. Hsa_circ_0003998 promotes cell proliferation and invasion by targeting miR-326 in non-small cell lung cancer. Onco Targets Ther. 2018;11:5569–77.
Qiao GL, Chen L, Jiang WH, Yang C, Yang CM, Song LN, Chen Y, Yan HL, Ma LJ. Hsa_circ_0003998 may be used as a new biomarker for the diagnosis and prognosis of hepatocellular carcinoma. Onco Targets Ther. 2019;12:5849–60.
Jin YP, Hu YP, Wu XS, Wu YS, Ye YY, Li HF, Liu YC, Jiang L, Liu FT, Zhang YJ, et al. miR-143-3p targeting of ITGA6 suppresses tumour growth and angiogenesis by downregulating PLGF expression via the PI3K/AKT pathway in gallbladder carcinoma. Cell Death Dis. 2018;9:182.
Xia C, Yang Y, Kong F, Kong Q, Shan C. MiR-143-3p inhibits the proliferation, cell migration and invasion of human breast cancer cells by modulating the expression of MAPK7. Biochimie. 2018;147:98–104.
Karimi L, Mansoori B, Shanebandi D, Mohammadi A, Aghapour M, Baradaran B. Function of microRNA-143 in different signal pathways in cancer: new insights into cancer therapy. Biomed Pharmacother. 2017;91:121–31.
Panza E, Ercolano G, De Cicco P, Armogida C, Scognamiglio G, Botti G, Cirino G, Ianaro A. MicroRNA-143-3p inhibits growth and invasiveness of melanoma cells by targeting cyclooxygenase-2 and inversely correlates with malignant melanoma progression. Biochem Pharmacol. 2018;156:52–9.
Guo L, Fu J, Sun S, Zhu M, Zhang L, Niu H, Chen Z, Zhang Y, Guo L, Wang S. MicroRNA-143-3p inhibits colorectal cancer metastases by targeting ITGA6 and ASAP3. Cancer Sci. 2019;110:805–16.
He Z, Yi J, Liu X, Chen J, Han S, Jin L, Chen L, Song H. MiR-143-3p functions as a tumor suppressor by regulating cell proliferation, invasion and epithelial-mesenchymal transition by targeting QKI-5 in esophageal squamous cell carcinoma. Mol Cancer. 2016;15:51.
Yang R, Xing L, Zheng X, Sun Y, Wang X, Chen J. The circRNA circAGFG1 acts as a sponge of miR-195-5p to promote triple-negative breast cancer progression through regulating CCNE1 expression. Mol Cancer. 2019;18:4.
Zhang PF, Wei CY, Huang XY, Peng R, Yang X, Lu JC, Zhang C, Gao C, Cai JB, Gao PT, et al. Circular RNA circTRIM33-12 acts as the sponge of MicroRNA-191 to suppress hepatocellular carcinoma progression. Mol Cancer. 2019;18:105.
Sun X, Dai G, Yu L, Hu Q, Chen J, Guo W. MiR-143-3p inhibits the proliferation, migration and invasion in osteosarcoma by targeting FOSL2. Sci Rep. 2018;8:606.
Patop IL, Wust S, Kadener S. Past, present, and future of circRNAs. EMBO J. 2019;38:e100836.
Meng Q, Rayala SK, Gururaj AE, Talukder AH, O'Malley BW, Kumar R. Signaling dependent and coordinated regulation of transcription, splicing, and translation resides in a single coregulator, PCBP1. Proc Natl Acad Sci. 2007;104:5866–71.
Brown AS, Mohanty BK, Howe PH. Identification and characterization of an hnRNP E1 translational silencing motif. Nucleic Acids Res. 2016;44:5892–907.
Zhang ZZ, Shen ZY, Shen YY, Zhao EH, Wang M, Wang CJ, Cao H, Xu J. HOTAIR long noncoding RNA promotes gastric Cancer metastasis through suppression of poly r(C)-binding protein (PCBP) 1. Mol Cancer Ther. 2015;14:1162–70.
Tripathi V, Sixt KM, Gao S, Xu X, Huang J, Weigert R, Zhou M, Zhang YE. Direct regulation of alternative splicing by SMAD3 through PCBP1 is essential to the tumor-promoting role of TGF-beta. Mol Cell. 2016;64:549–64.
Zhang T, Huang XH, Dong L, Hu D, Ge C, Zhan YQ, Xu WX, Yu M, Li W, Wang X, et al. PCBP-1 regulates alternative splicing of the CD44 gene and inhibits invasion in human hepatoma cell line HepG2 cells. Mol Cancer. 2010;9:72.
Li Q, Wang Y, Wu S, Zhou Z, Ding X, Shi R, Thorne RF, Zhang XD, Hu W, Wu M. CircACC1 Regulates Assembly and Activation of AMPK Complex under Metabolic Stress. Cell Metab. 2019;2:157-73.
We thank Department of The Third Department of Hepatic Surgery, Eastern Hepatobiliary Surgery Hospital, Naval military Medical University, Shanghai, China for providing the HCC tissue samples and related anonymous clinical data.
This work was supported by the National Natural Science Foundation of China (NO.81672335, 81672350 and 81872225), and the “medical professionals cross fund” of Shanghai Jiaotong University (NO. YG2016ZD10).
Ethics approval and consent to participate
All experimental protocols were approved by the Ethics Committee of Hepatic Surgery, Eastern Hepatobiliary Surgery Hospital and performed in accordance with the relevant guidelines and regulations. Written informed consents were obtained from all patients. All animal experiments were approved by the Institutional Animal Care and Use Committee of Naval Medical University (Shanghai, China).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The identification of HCC cells in this study by short tandem repeat.
The primer details and siRNA sequence used in this study.
Predicted target genes for hsa-miR-143-3p in TargetScan and miRanda.
Prediction score of RBPs binding to circ0003998 using the online catRAPID algorithm.
Predicted target genes of PCBP1 in Gene Transcription Regulation Database.
Tumor promoter in HCC which are targets of PCBP1.
(a) QRT-PCR detected the relative expression level of circ0003998 after over expressing or silencing circ0003998 in HCC cells. (b) The relative protein expression level of Snail and Slug in the HCC cells with circ0003998 over expression and silencing. (c) FISH was performed to observe the cellular location of circ0003998 (red), miR-143-3p (green) and PCBP1 (green) in HuH7 cells (magnification, × 400, scale bar, 25 μm). (d) Rescue transwell assay was performed after trans-fection with indicated vectors, OE-FOSL2 or si-FOSL2 (magnification, × 100, scale bar, 100 μm). (e) The relative protein expression level of DAB2 in HCC cells with circ0003998 over expression and silencing. (f) The relative mRNA expression of circ003998, FOSL2 and CD44v6 in HCC and PVTT tissues. *p-value< 0.05, **p-value< 0.01, ***p-value< 0.001.
About this article
Cite this article
Song, Ln., Qiao, Gl., Yu, J. et al. Hsa_circ_0003998 promotes epithelial to mesenchymal transition of hepatocellular carcinoma by sponging miR-143-3p and PCBP1. J Exp Clin Cancer Res 39, 114 (2020). https://doi.org/10.1186/s13046-020-01576-0
- Hepatocellular carcinoma