Circ_0035483 Functions as a Tumor Promoter in Renal Cell Carcinoma via the miR-31-5p-Mediated HMGA1 Upregulation

Abstract

Background

Renal cell carcinoma (RCC) that originates from the proximal renal tubules is the most common cancer of the human kidney. Increasing circRNA/miRNA/mRNA networks have been found in RCC regulation. This study will explore the regulatory relation of circular RNA (circRNA) circ_0035483, microRNA-31-5p (miR-31-5p) and high mobility group A1 (HMGA1).

Methods

The levels of circ_0035483, miR-31-5p and HMGA1 were measured by real-time polymerase chain reaction (qRT-PCR) or Western blot. Cell proliferation was determined using Cell Counting Kit-8 (CCK-8) and colony formation assays. Cell migration and invasion were assessed by transwell assay. HMGA1 and epithelial–mesenchymal transition (EMT)-related protein levels were quantified using Western blot. Glycolytic metabolism was evaluated by glucose consumption and lactate production. The interaction between targets was confirmed via dual-luciferase reporter, RNA immunoprecipitation (RIP) and RNA pull-down assays. In vivo experiment was performed through the establishment of xenograft models in mice.

Results

Circ_0035483 expression was upregulated in RCC tissues and cells. The inhibitory effects on RCC cell proliferation, migration, invasion, EMT and glycolysis were induced after circ_0035483 was downregulated. MiR-31-5p was identified as a target of circ_0035483 and miR-31-5p upregulation was related to the function of circ_0035483 knockdown in RCC cells. Additionally, miR-31-5p targeted HMGA1 and inhibited the malignant behaviors of RCC cells by negatively regulating HMGA1. Moreover, HMGA1 expression was regulated by circ_0035483 via targeting miR-31-5p. Circ_0035483 also affected tumor growth in vivo by relying on the miR-31-5p/HMGA1 axis.

Conclusion

These findings clarified that the tumor-promoting function of circ_0035483 in RCC was partly achieved by regulating the miR-31-5p/HMGA1 axis.

Keywords:

• circ_0035483
• renal cell carcinoma; RCC
• miR-31-5p
• HMGA1

Introduction

Renal cell carcinoma (RCC) accounts for more than 90% of cancers occurred in the genitourinary system, and it is a heterogeneous disease containing different kinds of histological subtypes.¹ Due to the asymptomatic property in the early phase, the targeted diagnosis has been emerging as an effective strategy for RCC screening.² Tumor metastasis of RCC is commonly found in human anatomical sites, such as lung, bone, liver, and brain.³ The current management for RCC therapy has been rapidly developed, especially the targeted therapy.⁴ The molecular biomarkers are pivotal for improving the clinical diagnosis and optimizing the therapeutic outcomes.⁵

Noncoding circular RNAs (circRNAs) can serve as the molecular sponges of microRNAs (miRNAs) in cancer cells.⁶ MiRNAs also can control the expression of downstream targets by binding to the 3ʹ-untranslated regions (3ʹ-UTRs) of messenger RNAs (mRNAs).⁷ Moreover, circRNAs regulate the gene expression through sponging miRNAs to play their regulatory roles in various kinds of cancers.⁸ For example, Mao et al have reported that circRNA_000554 blocked the epithelial-mesenchymal transition (EMT) of breast cancer cells via targeting miR-182 to regulate the ZFP36 level.⁹ Yin et al reported that circ_0046263 facilitated the malignant behaviors of nasopharyngeal carcinoma by upregulating the expression of IGFBP3 via binding to miR-133a-5p.¹⁰ Knockdown of circ_0109320 was shown to repress cell proliferation and invasion by the miR-595/E2F7 axis in non-small cell lung cancer.¹¹

A recent RCC research indicated that circ_0035483 was upregulated in RCC and its promoting effect on tumor growth was associated with the miR-335/CCNB1 axis.¹² The roles of circ_0035483 in cell migration, invasion and metabolic process of RCC are unreported, and its other molecular mechanisms need to be investigated. MicroRNA-31-5p (miR-31-5p) was downregulated in RCC and played as a suppressive role by targeting FLOT1 or CDK1.^13,14 Dong et al found that high mobility group A1 (HMGA1) was overexpressed in RCC and generated the oncogenic function in cell progression of RCC.¹⁵

This study hypothesizes that circ_0035483 targets miR-31-5p and miR-31-5p targets HMGA1, and further explores the regulatory potential of circ_0035483 on HMGA1 by targeting miR-31-5p. The key point of this study is to discover the circ_0035483/miR-31-5p/HMGA1 axis in the development of RCC.

Materials and Methods

Tissue Collection

60 patients with RCC have participated in the current study and no therapy was performed for them before the surgical resection at Central Hospital of Baoji. RCC tissues (n=60) and coterminous normal tissues (n=60, more than 3 cm from RCC tissues) were carefully collected during the surgery. Instantly, these fresh tissues were immerged in RNAlater^® RNA Stabilization Agent (Qiagen, Hilden, Germany) for stabilizing RNA and then saved in liquid nitrogen. 60 patients have provided the written informed consent forms for this research. Moreover, all operating procedures followed the Declaration of Helsinki involving human subjects and got ratification from the Ethics Committee of Central Hospital of Baoji.

Cell Culture

Human RCC cell lines 786-O and CaKi-1 were chosen for RCC research in our study, using human kidney tube cell line HK-2 as the normal control. Three cell lines were all purchased from American Type Culture Collection (ATCC, Manassas, VA, USA). Cell medium was prepared by mixing Roswell Park Memorial Institute-1640 (RPMI-1640; Gibco, Carlsbad, CA, USA), fetal bovine serum (FBS; 10%, Gibco) and penicillin/streptomycin solution (1%, Gibco). Cells were cultured using the nutrient medium in a humidified incubator with 5% CO₂ at 37°C, using.

Cell Transfection

Small interfering RNAs (si-circ_0035483#1, si-circ_0035483#2, and negative control si-NC), miRNA mimics (miR-31-5p and miR-NC), miRNA inhibitors (anti-miR-31-5p and anti-miR-NC), and lentiviral short hairpin RNA (shRNA) vectors (sh-circ_0035483 and sh-NC) were bought from GenePharma (Shanghai, China). In addition, an empty pcDNA expression vector was acquired from Invitrogen (Carlsbad, CA, USA) and the recombinant pcDNA-HMGA1 vector was constructed for HMGA1 overexpression. After 786-O and CaKi-1 cell density reached 60%, the mixture of oligonucleotides or vectors, Lipofectamine™ 3000 Reagent (Invitrogen) and Opti-MEM™ I Reduced Serum Medium (Gibco) was added following the operating instruction for users.

RNA Acquisition and Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR)

Total RNA was extracted from the tissues samples and cells using TRIzol™ Reagent (Invitrogen), and the reverse transcription was performed via QuantiTect Reverse Transcription Kit (Qiagen) according to the manufacturers’ guidelines. Then QuantiTect SYBR^® Green PCR Kit (Qiagen) was used for RNA quantification on ABI 7500 Real-Time PCR system (Applied Biosystems, Foster City, CA, USA). The sequences for primers were shown as below: forward, 5ʹ-TCATCCTTCTAAGGTAGCAGCAG-3ʹ and reverse, 5ʹ-CTTCGACGCTGTAGCCTTGT-3ʹ (circ_0035483); forward, 5ʹ-ACACTCCAGCTGGGAGGCAAGATGCTGGC-3ʹ and reverse, 5ʹ-TGGTGTCGTGGAGTCG-3ʹ (miR-31-5p); forward, 5ʹ-CATCCGCATTTGCTACCAGC-3ʹ and reverse, 5ʹ-GGGCTCCTTCTGACTCCCTA-3ʹ (HMGA1); forward: 5ʹ-CGGGAAATCGTGCGTGACAT-3ʹ and reverse: 5ʹ-AACCGCTCATTGCCAATGGT-3ʹ (β-actin); forward: 5ʹ-CTCGCTTCGGCAGCACA-3ʹ and reverse: 5ʹ-ACGCTTCACGAATTTGCGT-3ʹ (U6). β-actin and U6 acted as the housekeeping genes for circ_0035483 or HMGA1 and miR-31-5p. The comparative cycle threshold (2^−∆∆Ct) method was applied for the analysis of the relative expression.

Cell Counting Kit-8 (CCK-8) Assay

The absorbance (450 nm) was measured at 0 h, 24 h, 48 h and 72 h to assess the proliferative ability of transfected 786-O and CaKi-1 cells. After cell transfection, 10 μL CCK-8 solution (Dojindo, Kumamoto, Japan) was respectively added into each well in the 96-well plates, followed by 4 h-incubation at 37°C and the detection under the microplate reader (Bio-Rad, Hercules, CA, USA).

Colony Formation Assay

Colony formation assay was also used for cell proliferation analysis. 786-O and CaKi-1 cells (200 cells/well) were transplanted into a 6-well plate after transfection for 48 h. Two weeks later, these cell colonies were fastened in 4% paraformaldehyde (Invitrogen) and dyed in 0.1% crystal violet (Solarbio, Beijing, China). Then the staining fluid was discarded and the colonies were washed in Phosphate Buffered Saline (PBS; Gibco), followed by counting the colonies manually.

Transwell Assay

The transwell chamber (Corning Inc., Corning, NY, USA) was applied to perform the migration assay and transwell chamber pre-coated with matrigel (Corning Inc.) was used in the invasion assay. Briefly, 786-O and CaKi-1 cells in serum-free RPMI-1640 medium were respectively seeded into the top chamber and 10% FBS+RPMI-1640 medium was pipetted into the bottom chamber. After cell culture for 24 h in the normal condition, cells in the top chamber passed into the bottom chamber were considered as the migratory and invasive cells. 4% paraformaldehyde (Invitrogen) and 0.1% crystal violet (Solarbio) were used for cell fixation and staining. Eventually, cell images were acquired with 100 × magnification under the inverted microscope (Olympus, Tokyo, Japan), and cell number of each group was counted.

Western Blot

Protein extraction was conducted by Radioimmunoprecipitation assay (RIPA) Lysis Buffer for tissues and cells (Solarbio). Protein products were electrophoresed on sodium dodecyl sulfate polyacrylamide gel (Invitrogen) and transferred onto the Polyvinylidene Fluoride (PVDF) membranes (Invitrogen), then the membranes were soaked in 5% nonfat dry milk (Cell Signaling Technology, CST; Boston, MA, USA) for protein blocking. After the incubation of primary antibodies at 4°C overnight and secondary antibody (CST, #7074, 1:3000) at room temperature for 1 h, the protein blots were discerned using ECL Chemiluminescent Substrate Reagent Kit (Invitrogen) and the expression analysis was performed by Image Lab software (Bio-Rad). Primary antibodies used in this study were listed as follows: E-cadherin (CST, #3195, 1:1000), Vimentin (CST, #5741, 1:1000), Snail (CST, #3879, 1:1000), HMGA1 (CST, #7777, 1:1000), β-actin (CST, #4970, 1:1000).

Detection of Glucose Consumption and Lactate Production

Glucose Assay Kit-WST^® and Lactate Assay Kit-WST^® (Dojindo) were respectively applied for detecting the consumption of glucose and the production of lactate. At 48 h post-transfection, cell supernatant was pipetted into a new tube and diluted with ultrapure water. After cell seeding into the 96-well plates, the levels of glucose and lactate were determined as per the users’ manuals.

Dual-Luciferase Reporter Assay

The sequence of circ_0035483 (wild-type, WT) containing the binding sites for miR-31-5p was inserted into the pmirGLO vector (Promega, Madison, WI, USA) to construct the recombinant luciferase plasmidscirc_0035483-WT. These binding sites of miR-31-5p were mutated in the circ_0035483 sequence, and the mutant-type (MUT) sequence was used for the construction of circ_0035483-MUT. HMGA1 3ʹUTR WT and MUT sequences were also cloned into the pmirGLO vector to obtain the luciferase plasmids HMGA1 3ʹUTR-WT and HMGA1 3ʹUTR-MUT. 786-O and CaKi-1 cells were plated into the well-plates, then co-transfected with the above each plasmid and miR-31-5p or miR-NC. The firefly luciferase activity (normalized by the renilla luciferase activity) was detected using the Dual-Luciferase Reporter Assay System (Promega), referring to the producer’s specification.

RNA Immunoprecipitation (RIP) Assay

786-O and CaKi-1 cells were harvested and incubated with antibody-conjugated magnetic beads according to the Dynabeads™ Protein A Immunoprecipitation Kit (Invitrogen). Total RNA was isolated from the magnetic beads of Anti-Argonaute-2 (Anti-Ago2) or Anti-immunoglobulin G (Anti-IgG) group, then the qRT-PCR was carried out for the examination of circ_0035483, miR-31-5p and HMGA1 levels.

Biotin-Coupled miR-31-5p Capture

786-O and CaKi-1 cells were transfected with the biotin-miR-31-5p mimic (Bio-miR-31-5p) or Bio-miR-NC, and cell lysates were incubated with the magnetic beads in Pierce™ Magnetic RNA-Protein Pull-Down Kit (Thermo Fisher Scientific, Waltham, MA, USA). The enrichment of circ_0035483 was quantified by qRT-PCR after the extraction of total RNA.

In vivo Assay

Twelve five-week-old male BALB/c nude mice were provided by Vital River Laboratory Animal Technology Co. Ltd (Beijing, China) and randomly divided into two groups. Circ_0035483 was stably knocked down by the construction of 786-O cell strain with stable expression of sh-circ_0035483, with sh-NC as the negative control group. 5 × 10⁶ cells were subcutaneously injected into the right flank of mice back. Tumor volume (length × width² × 0.5) was measured weekly and all those mice were euthanatized by CO₂ asphyxia method (displacing 30% air in the cage each minute) after 35 days. Tumors were excised from mice and weighed on the electronic scale. Then the expression levels of circ_0035483, miR-31-5p and HMGA1 in tumor tissues were assayed by qRT-PCR or Western blot. This animal assay was performed at Central Hospital of Baoji in accordance with the Guidelines of Laboratory Animals of National Institutes of Health (NIH) and approved by Animal Ethical Committee of Central Hospital of Baoji.

Statistical Analysis

Our data were revealed as the mean ± standard error of mean (SEM). The linear relationship was analyzed using Pearson’s correlation coefficient. SPSS 24.0 was used for the statistical analysis in this study. Student’s t-test or one-way analysis of variance (ANOVA) followed by Tukey’s test were performed for analyzing the difference of two groups or multiple groups. P<0.05 was in the presence of significant difference.

Results

Circ_0035483 Level Was Overexpressed in RCC Tissues and Cells

The expression detection of circ_0035483 was performed by qRT-PCR in tissues and cells. Contrasted with normal tissues, circ_0035483 expression was significantly upregulated in RCC tissues (Figure 1A). Also, the level of circ_0035483 was higher in 786-O and CaKi-1 cells than that in normal HK-2 cells (Figure 1B). The dysregulation of circ_0035483 implied its potential association with the cellular regulation of RCC.

Figure 1 Circ_0035483 level was overexpressed in RCC tissues and cells. (A and B) The relative circ_0035483 expression was examined by qRT-PCR in RCC tissues and cells (786-O and Caki-1), as well as the normal tissues and HK-2 cells. *P < 0.05.

Silencing Circ_0035483 Retarded Proliferation, Migration, Invasion and Glycolysis of RCC Cells

After transfection of si-circ_0035483#1 or si-circ_0035483#2 in 786-O and CaKi-1 cells, we observed the prominent downregulation of circ_0035483 expression relative to transfection of si-NC (Figure 2A). As the consequences of circ_0035483 knockdown, the inhibitory effects were found on cell proliferation by performing CCK-8 (Figure 2B–C) and colony formation (Figure 2D) assays, as well as cell migration and invasion by transwell assay (Figure 2E–F). Inhibition of E-cadherin and upregulation of Vimentin or Snail are the common indications for EMT.¹⁶ E-cadherin protein expression was increased but Vimentin and Snail protein levels were reduced in the si-circ_0035483#1 and si-circ_0035483#2 groups compared to si-NC group (Figure 2G), showing that EMT process was blocked by knockdown of circ_0035483. In addition, the low circ_0035483 expression induced the repression of glucose consumption (Figure 2H) and lactate production (Figure 2I) in 786-O and CaKi-1 cells. These experimental results indicated that the silence of circ_0035483 suppressed the malignant behaviors in RCC cells.

Figure 2 Silencing circ_0035483 retarded proliferation, migration, invasion and glycolysis of RCC cells. Si-NC, si-circ_0035483#1 or si-circ_0035483#2 was transfected into 786-O and CaKi-1 cells. (A) The qRT-PCR was used to measure the expression level of circ_0035483. (B–D) CCK-8 assay (B and C) and colony formation assay (D) were performed for detecting the proliferative ability. (E and F) Transwell assay was used for migration/invasion detection. (G) Western blot was used for the detection of EMT markers. (H and I) Glucose consumption and lactate production were determined for the glycolytic evaluation. *P < 0.05.

MiR-31-5p Was a miRNA Target of Circ_0035483

The miRNA targets of circ_0035483 were predicted by bioinformatic softwares circinteractome, circBank and Starbase3.0. Venn diagram was used for the further screening of miRNAs. As depicted in Figure 3A, three miRNAs (miR-31-5p, miR-432-5p, miR-582-5p) were mutually predicted by three softwares and considered as the candidate targets of circ_0035483. The qRT-PCR exhibited miR-31-5p and miR-432-5p were downregulated in RCC tissues compared to normal tissues while miR-582-3p expression was unchanged in RCC tissues contrasted with normal tissues (Figure 3B). The expression change of miR-31-5p was more significant than miR-432-5p, and we then further explored the target relation between circ_0035483 and miR-31-5p. Circinteractome has shown the binding sites between the sequences of circ_0035483 and miR-31-5p as the presentation in Figure 3C. Overexpression of miR-31-5p by miR-31-5p mimic (Figure 3D) resulted in the inhibition of the luciferase activity of circ_0035483-WT group, while luciferase activity of circ_0035483-MUT group was unaffected in 786-O and CaKi-1 cells (Figure 3E–F). In RIP assay, the higher enrichment of circ_0035483 and miR-31-5p was also observed in Anti-Ago group than those in Anti-IgG group (Figure 3G–H). Moreover, RNA pull-down assay indicated that circ_0035483 was largely pulled down by miR-31-5p (Figure 3I). Also, miR-31-5p level was obviously promoted after the introduction of si-circ_0035483#2 in 786-O and CaKi-1 cells in contrast with si-NC group (Figure 3J). There was a downregulating level of miR-31-5p in RCC tissues and cells compared to the normal tissues and cells (Figure 3K–L). The above data identified miR-31-5p as a miRNA target of circ_0035483.

Figure 3 MiR-31-5p was a miRNA target of circ_0035483. (A) The mutual miRNAs were screened by Veen diagram after the target prediction of circ_0035483 using circinteractome, circBank and Starbase3.0. (B) The detection of miR-31-5p, miR-432-5p, miR-582-5p was carried out by qRT-PCR in RCC tissues. (C) The binding sites of miR-31-5p were presented in the sequence of circ_0035483 by circinteractome. (D) The transfection efficiency of miR-31-5p mimic was analyzed using qRT-PCR. (E–I) The binding between miR-31-5p and circ_0035483 was performed using the dual-luciferase reporter assay (E and F), RIP (G and H) and RNA pull-down assay (I). (J) The level of miR-31-5p was assayed via qRT-PCR in si-NC or si-circ_0035483#2 transfected 786-O and CaKi-1 cells. (K and L) The qRT-PCR was applied for detecting the miR-31-5p expression in RCC tissues (K) and cells (L). *P < 0.05.

Circ_0035483 Affected the Malignant Behaviors of RCC Cells via Negatively Regulating miR-31-5p

The expression inhibition for miR-31-5p was achieved by anti-miR-31-5p transfection in 786-O and CaKi-1 cells (Figure 4A). The relation of miR-31-5p and circ_0035483 was explored in regulating RCC cellular processes. As expected, si-circ_0035483#2-induced suppressive influences on cell proliferation (Figure 4B–D), migration and invasion (Figure 4E–F) were countervailed by the transfection of anti-miR-31-5p. The alleviative phenomenons mediated by anti-miR-31-5p for si-circ_0035483#2 were also found in EMT process (Figure 4G–H) and glycolytic metabolism (Figure 4I–J) of 786-O and CaKi-1 cells. Therefore, the effects of circ_0035483 on RCC cells were associated with the negative regulation of miR-31-5p.

Figure 4 Circ_0035483 affected the malignant behaviors of RCC cells via negatively regulating miR-31-5p. (A) The miR-31-5p expression was quantified by qRT-PCR after transfection of anti-miR-NC or anti-miR-31-5p. (B–H) The cellular analysis was performed by CCK-8 assay (B and C) and colony formation assay (D) for cell proliferation, transwell assay for cell migration (E) and invasion (F), Western blot for EMT protein detection (G and H) in si-NC, si-circ_0035483#2, si-circ_0035483#2+anti-miR-NC or si-circ_0035483#2+anti-miR-31-5p. (I and J) Glycolysis was analyzed by examining the glucose consumption (I) and lactate production (J) in the above four groups. *P < 0.05.

HMGA1 Was a Downstream Target of miR-31-5p

Starbase3.0 has manifested that HMGA1 3ʹUTR sequence contained the miR-31-5p-binding sites (Figure 5A). Subsequently, the dual-luciferase reporter assay (Figure 5B–C) and RIP assay (Figure 5D–E) validated the combination between HMGA1 and miR-31-5p in 786-O and CaKi-1 cells. Our qRT-PCR and Western blot results further suggested the negative modulation of HMGA1 by miR-31-5p because of the HMGA1 upregulation by miR-31-5p inhibition and HMGA1 downregulation by miR-31-5p overexpression (Figure 5F–G). HMGA1 mRNA and protein levels were found to be evidently enhanced in RCC tissues (Figure 5H–I) and cells (Figure 5J–K). Taken together, HMGA1 was a target in the downstream of miR-31-5p.

Figure 5 HMGA1 was a downstream target of miR-31-5p. (A) Starbase3.0 showed the binding sites of miR-31-5p in HMGA1 3ʹUTR sequence. (B–E) Dual-luciferase reporter assay (B and C) and RIP assay (D and E) were conducted for the analysis of interaction between miR-31-5p and HMGA1. (F and G) The qRT-PCR and Western blot were adopted for assessing the effects of miR-31-5p inhibition or overexpression on HMGA1 mRNA and protein levels in 786-O and CaKi-1 cells. (H–K) HMGA1 expression was assayed by mRNA level via qRT-PCR and protein level via Western blot in RCC/Normal tissues (H and I) and RCC/HK-2 cells (J and K). *P < 0.05.

MiR-31-5p Downregulated HMGA1 to Impede the Progression of RCC

The qRT-PCR and Western blot showed that the mRNA and protein levels of HMGA1 were increased in 786-O and CaKi-1 cells transfected with pcDNA-HMGA1 compared to cells transfected with pcDNA (Figure 6A–B). Functional analysis revealed that miR-31-5p overexpression led to the distinct inhibition in cell proliferation (Figure 6C-E), migration (Figure 6F) and invasion (Figure 6G), which was ameliorated following the transfection of pcDNA-HMGA1. In the EMT process, the overexpression of HMGA1 also rescued the miR-31-5p-induced upregulation of E-cadherin and reduction of Vimentin or Snail (Figure 6H–I). Similarly, the inhibitory influences of miR-31-5p on glucose consumption and lactate production were reversed after HMGA1 expression was promoted (Figure 6J–K). MiR-31-5p acted as a repressor in RCC via targeting HMGA1.

Figure 6 MiR-31-5p downregulated HMGA1 to impede the progression of RCC. (A and B) The mRNA and protein levels of HMGA1 were measured by qRT-PCR and Western blot in 786-O and CaKi-1 cells transfected with pcDNA or pcDNA-HMGA1. (C–I) Functional experiments were conducted by cell proliferation using CCK-8 assay (C and D) and colony formation assay (E), cell migration (F) and invasion (G) using transwell assay, EMT protein detection (H and I) using Western blot after transfection of miR-NC, miR-31-5p, miR-31-5p+pcDNA or miR-31-5p+pcDNA-HMGA1. (J and K) The assessment of glycolytic metabolism was carried out by glucose consumption (J) and lactate production (K) in four transfection groups. *P < 0.05.

Circ_0035483 Could Regulate HMGA1 by Serving as a Sponge of miR-31-5p

To analyze the potential association of circ_0035483 and HMGA1 via miR-31-5p, HMGA1 expression was detected after 786-O and CaKi-1 cells were transfected with si-circ_0035483#2, si-circ_0035483#2+anti-miR-31-5p or the matched controls. As illustrated in Figure 7A–B, anti-miR-31-5p transfection evidently abolished the HMGA1 expression inhibition caused by knockdown of circ_0035483. Pearson’s correlation coefficient has shown the negative relation between miR-31-5p and circ_0035483 (r=−0.465, P=0.0002) or HMGA1 (r=−0.331, P=0.0099), as well as the positive correlation between circ_0035483 and HMGA1 (r=0.313, P=0.0148) (Figure 7C–E). Collectively, HMGA1 could be regulated by circ_0035483 via targeting miR-31-5p.

Figure 7 Circ_0035483 could regulate HMGA1 by serving as a sponge of miR-31-5p (A and B) After 786-O and CaKi-1 cells were transfected with si-circ_0035483#2, si-circ_0035483#2+anti-miR-31-5p or the corresponding controls, the detection of HMGA1 mRNA and protein expression was performed using qRT-PCR and Western blot. (C–E) The linear correlations among circ_0035483, miR-31-5p and HMGA1 expression in RCC tissues were analyzed by Pearson’s correlation coefficient. *P < 0.05.

RCC Tumor Growth in vivo Was Reduced by Downregulation of Circ_0035483 via the miR-31-5p/HMGA1 Axis

Xenograft tumor assay demonstrated that tumor volume (Figure 8A) and weight (Figure 8B) were decreased in mice of sh-circ_0035483 group by comparison to that of sh-NC group. After the RNA isolation and qRT-PCR analysis, we noticed that sh-circ_0035483 triggered the repression of circ_0035483 expression (Figure 8C) and an opposite effect on miR-31-5p level (Figure 8D) in contrast with sh-NC group. In addition, HMGA1 mRNA and protein levels were reduced by the downregulation of circ_0035483 in mice (Figure 8E-F). In vivo, downregulating circ_0035483 could inhibit tumor growth of RCC via regulating the miR-31-5p/HMGA1 axis.

Figure 8 RCC tumor growth in vivo was reduced by downregulation of circ_0035483 via the miR-31-5p/HMGA1 axis. (A and B) Tumor growth was assessed by measuring tumor volume (A) and weight (B) in xenograft model of sh-NC or sh-circ_0035483 group. (C and D) The qRT-PCR was applied to analyze the relative levels of circ_0035483 (C) and miR-31-5p (D) in tumor tissues. (E and F) The quantification of HMGA1 was performed by qRT-PCR and Western blot in tumor tissues. *P < 0.05.

Discussion

RCC is one of the most lethal urological malignancies with high incidence and mortality rates.¹⁷ It is important to know the pathological mechanisms and discover biomolecules for RCC. Circ_0035483 was found to inhibit the malignant development of RCC by targeting the miR-31-5p/HMGA1 axis, which disclosed a specific functional mechanism of circ_0035483 and the potential of circ_0035483 as a biological target in RCC.

CircRNAs have been reported as crucial regulators in diverse human cancers, such as breast cancer, gastric cancer, lung cancer, and colorectal cancer.^18–21 Cellular processes of RCC were reported to be affected by different circRNAs. For instance, overexpression of circ_0001368 inhibited RCC cell growth and invasion by regulating the levels of miR-492 and LATS2;²² hsa-circ-0072309 acted as an anti-tumor molecule via targeting miR-100 to inactivate the downstream signal pathways;²³ circ-AKT3 inhibited the migratory and invasive abilities of RCC cells via acting on the miR-296-3p/E-cadherin network.²⁴ In addition to the tumor-inhibitory roles, circRNAs also can exert the tumor-promoting effect on RCC. Li et al discovered that hsa_circ_0054537 contributed to the proliferation and migration in RCC cells by sponging miR-130a-3p to downregulate the cMet expression.²⁵ Downregulation of circ_001895 or circ_000926 impeded cell growth and invasion through targeting the miR-296-5p/SOX12 or miR-411/CDH2 axis in RCC cells.^26,27 In this report, the experimental results denoted that RCC cell proliferation, migration, invasion and EMT process were all suppressed with the inhibition of circ_0035483 expression.

Glycolysis is a frequent metabolic pathway in the developing process of tumor cells, accompanied by the heavy consumption of glucose and production of lactate under aerobic condition.²⁸ Glycolytic activity is increased in proliferating cells and glycolysis is also used as a target site in the therapy of kidney cancer.²⁹ Carcinogenic circRNAs have been found to promote glycolysis in cancer cells, such as circMAGI3 in non-small cell lung cancer,³⁰ circMAT2B in gastric cancer³¹ and circDENND4C in colorectal cancer.³² However, the association of circRNAs and glycolysis in RCC is not fully explicit. The downregulation of circ_0035483 suppressed the glucose consumption and lactate production in RCC cells here, showing that circ_0035483 overexpression facilitated the glycolytic metabolism in RCC.

Regarding the researches of circRNAs as the sponges of miRNAs, Jin et al have declared that circ_0000523 retarded the evolution of colorectal cancer via miR-31³³ and circ_0001564 served for a miR-29c-3p sponge to mediate the tumorigenicity in osteosarcoma.³⁴ After the online prediction for the miRNA targets of circ_0035483 and Veen diagram analysis, three miRNAs were selected as the candidates and miR-31-5p was used for the final research subject after the expression analysis in RCC. Our data validated the combination between circ_0035483 and miR-31-5p, and the negative regulation of circ_0035483 on miR-31-5p in RCC cells. Moreover, the revert experiments manifested that the regulation of circ_0035483 in RCC cells was achieved by targeting miR-31-5p at least in part. Meanwhile, HMGA1 was demonstrated to be a target of miR-31-5p in RCC. The anti-cancer effects of miR-31-5p in RCC cell proliferation, migration, invasion and glycolysis were dependent on inhibiting the expression of HMGA1. More interestingly, circ_0035483 could regulate the level of HMGA1 by targeting miR-31-5p. The previous studies have identified the circRNA/miRNA/mRNA networks in many cancers.^35–37 Not only in vitro, the circ_0035483/miR-31-5p/HMGA1 axis was also found in regulating RCC tumor growth in vivo.

Conclusion

Altogether, circ_0035483 contributed to the cellular malignant behaviors in RCC via the action on the miR-31-5p/HMGA1 axis. These findings indicated that circ_0035483 might have great clinical values as a useful biomarker for the early diagnosis and treatment of RCC patients.

Highlights

1. Silence of circ_0035483 suppresses RCC cell proliferation, migration, invasion and glycolytic metabolism.

2. Circ_0035483/miR-31-5p and miR-31-5p/HMGA1 axes are identified in the regulation of RCC cell malignant behaviors.

3. RCC tumor growth in vivo was regulated by circ_0035483/miR-31-5p/HMGA1 axis.

Ethics Approval and Consent Participate

Written informed consent was obtained from patients with approval by the Institutional Review Board in Central Hospital of Baoji.

Disclosure

The authors declare that they have no financial conflicts of interest.