CircRNA has_circ_0001806 promotes hepatocellular carcinoma progression via the miR-193a-5p/MMP16 pathway

Abstract

Reportedly, circular RNAs (circRNAs) are crucial regulators in cancer progression. Nonetheless, the molecular mechanism of circRNAs in hepatocellular carcinoma (HCC) has not been fully clarified. Gene expression omnibus (GEO) database was employed to screen out the differentially expressed circRNAs in HCC. qRT-PCR and western blot were executed to detect circ_0001806 expression, miR-193a-5p expression, and MMP16 mRNA and protein expressions in HCC. The effect of circ_0001806 on HCC was analyzed by the CCK-8 method and Transwell experiment. RIP assay, pull-down experiment, and dual-luciferase reporter gene experiment were applied to validate the targeting relationships among circ_0001806, miR-193a-5p, and MMP16. Circ_0001806 was up-modulated in HCC tissues and cell lines. Knockdown of circ_0001806 impeded the multiplication, migration, and invasion of HCC cells. Circ_0001806 could up-regulate MMP16 expression through repressing miR-193a-5p, thereby facilitating the malignant biological behaviors of HCC. Circ_0001806 promoted HCC progression by regulating miR-193a-5p/MMP16 axis.

Circ_0001806; Hepatocellular carcinoma; MMP16; miR-193a-5p

Introduction

Liver cancer has become the sixth most frequently diagnosed cancer and the fourth leading cause of cancer-related death in the world (¹). Hepatocellular carcinoma (HCC) accounts for 75-85% of cases of primary liver cancer. Since early HCC is usually asymptomatic, HCC patients are usually at an advanced stage when diagnosed and have missed the opportunity for curative treatment (liver resection or liver transplantation) (²,³). The high recurrence rate after radical resection and chemoresistance of HCC are major challenges in the treatment of HCC (³,⁴). To improve the diagnosis, prevention, and treatment of HCC, in-depth research into the molecular mechanism of this disease is therefore of great importance.

Circular RNAs (circRNAs) are RNA transcripts with a covalently closed circular structure that lack a 5′ cap and 3′ tail. The sequence of circRNAs is conserved, the structure of circRNAs is stable, and the expression of circRNA is tissue-specific (⁵). These characteristics endow circRNAs with many potential functions, such as acting as a microRNA (miRNA) sponge, or combining with proteins to form RNA-protein complexes to modulate the function of proteins (⁶). Accumulating research has demonstrated that circRNAs are related to the pathogenesis of tumors (⁷- ¹⁰). Importantly, circRNA can function as competitive endogenous RNA (ceRNA) to regulate the expression level of miRNAs and target genes (⁶). For instance, circSLC8A1 is reported to directly interact with miR-130b or miR-494, acting as an miRNA sponge to modulate PTEN expression and downstream signaling pathways, thereby restraining bladder cancer progression (¹¹). Nevertheless, the role and mechanism of circRNA in HCC have not been fully clarified.

In this work, we identified a circRNA, has_circ_0001806, that was significantly up-modulated in HCC tissues based on the circRNA microarray data. Further research revealed that circ_0001806 was highly expressed in HCC tissues and cell lines. Circ_0001806 knockdown impeded the growth, migration, and invasion of HCC cells. Additionally, circ_0001806 worked as a ceRNA for miR-193a-5p to up-modulate matrix metallopeptidase 16 (MMP16).

Material and Methods

Tissues specimens

This study was approved by the Ethics Committee of the Affiliated Hospital of Jianghan University. Tissue specimens were collected from patients who were diagnosed with HCC between May 2014 and May 2019 at the Affiliated Hospital of Jianghan University, and stored in liquid nitrogen at -196°C. Written informed consent was signed by each patient. Two HCC microarray datasets (GSE94508 and GSE97332) were downloaded from the Gene Expression Omnibus (GEO) database.

Cell lines, cell culture, and transfection

Normal hepatic cell line (L02) and HCC cell lines (Huh-7, HepG2, SMMC-7721, Bel-7402, and HL-7702) were procured from the Cell Bank of the Chinese Academy of Sciences (China). Cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM) (Thermo Fisher, USA) supplemented with 10% fetal bovine serum (FBS) (Gibco, USA), 100 μg/mL streptomycin and 100 U/mL penicillin (Invitrogen, USA) at 37°C in 5% CO_2. In a 6-well plate, the HCC cells were planted at 3×10⁵ cells/well. The cells were transfected with the negative control of siRNA (si-NC, 5′-CCAGUUUACCUAACGCAAUTT-3′); siRNA targeting circ_0001806 (si-circ_0001806#1: 5′-GTCTCCCAGTGCTCCAGACAA-3′; si-circ_0001806#2: 5′-TCCCAGTGCTCCAGACAATGA-3′; si-circ_0001806#3: 5′-TCTCCCAGTGCTCCAGACAAT-3′); negative control of miR-193a-5p mimic (mimic-NC, 5′-UCACAACCUCCUAGAAAGAGUAGA-3′); miR-193a-5p mimic (miR-193a-5p, 5′- UGGGUCUUUGCGGGCGAGAUGA-3′); miR-193a-5p inhibitor (anti-miR-193a-5p, 5′-UCAUCUCGCCCGCAAAGACCCA-3′); and negative control of miR-193a-5p inhibitor (anti-miR-control, 5′-CAGUACUUUUGUGUAGUACAA-3′) using Lipofectamine^® 2000 (Invitrogen), when the cell confluence reached 90%. The aforementioned oligonucleotides were available from GenePharma Co., Ltd. (China).

Quantitative real-time polymerase chain reaction (qRT-PCR)

Total RNA in tissues and cells was isolated by TRIzol reagent (Invitrogen), and the purity and concentration of RNA were measured by NanoDrop 2000 (Thermo Fisher Scientific Inc.). With the ReverTra Ace qPCR RT Package (Toyobo, Japan), RNA was transcribed into cDNA. With cDNA as template, qRT-PCR was conducted with THUNDERBIRD SYBR qPCR Mix (Toyobo). The gene expression was normalized by GAPDH/U6 expression using the 2^−ΔΔCt method. The primers designed for this research are displayed in Table 1. To validate the circular structure of circ_0001806, total RNA (5 μg) was incubated at 37°C for 15 min with 3 U/μg of RNase R (Epicentre Biotechnologies, China), and subsequently qRT-PCR was performed, with GAPDH as the control of linear RNA.

Table 1
Primer sequences used in this study.

Gene	Sequence
circ_0001806	Forward: 5-GTGATCTGAAAGGGCCAGAG-3
	Reverse: 5-TCCACATCACCCTTCACCTT-3
MMP16	Forward: 5-TAGAGCGTGCAGATAATGACAAGGA-3
	Reverse: 5-TGAACTGCTAGCCTCTGGATTTGA-3
miR-193a-5p	Forward: 5-TATATGGGTCTTTGCGGGCG-3
	Reverse: 5-GTGCAGGGTCCGAGGT-3
GAPDH	Forward: 5-GGAGCGAGATCCCTCCAAAAT-3
	Reverse: 5-GGCTGTTGTCATACTTCTCATGG-3
U6	Forward: 5-ATTGGAACGATACAGAGAAGATT-3
	Reverse: 5-GGAACGCTTCACGAATTTG-3

Cell counting kit-8 (CCK-8) experiment

Huh-7 and HepG2 cells were seeded onto 96-well plates (1.0×10³ cells/well) and cultured for 0, 24, 48, and 72 h. Then, 10 μL of CCK-8 solution (Dojindo, Japan) was added to each well. Cells were incubated for 1 h, and the absorbance (at 450 nm wavelength) of each well was examined by a microplate reader (Bio-Rad, USA).

Migration and invasion experiments

Approximately 1×10⁴ transfected cells were suspended in 200 μL of serum-free medium and positioned in the upper compartment of each Transwell chamber (8-μm pore size, USA). A medium containing 10% FBS was added to the lower compartment as the chemoattractant. Cells were incubated at 37°C with 5% CO₂ for 48 h for the invasion experiment and 24 h for the migration experiment. Then, cells in the upper compartment were removed with cotton swabs and cells on the lower surface of the filter were fixed with methanol and stained with 0.1% crystal violet. Within 3 random areas, the number of the cells was counted under a microscope. Matrigel (BD Biosciences, USA) was used to cover the filter in the invasion assay, but was not used in the migration assay.

Dual-luciferase reporter experiment

HEK293T cells were planted onto 96-well plates (1×10⁴ cells/well), cultured at 37°C for 12 h, and then co-transfected with wild type (WT) or mutant type (MUT) luciferase reporter vectors and miR-193a-5p mimics or mimic-NC using Lipofectamine^® 2000 (Invitrogen). After 48 h, the Firefly and Renilla luciferase activities were tested by a Dual-Luciferase Reporter Assay kit (Promega, USA).

RNA immunoprecipitation (RIP) experiment

The RIP experiment was executed with an EZ-Magna RIP kit (EMD Millipore, USA). RIP lysis buffer plus cocktail (Roche Diagnostics, China) was applied to lyse Huh-7 and HepG2 cells. Supernatants were incubated with antibody anti-AGO2 or IgG (EMD Millipore) along with protein A/G magnetic beads. Next, proteinase K and DNase (Beyotime, China) were used to remove proteins and DNA from the mixture. Subsequently, an RNeasy MinElute Cleanup Kit (Qiagen, China) was used to purify the RNA, and the purified RNA was analyzed through qRT-PCR.

RNA pull-down experiment

Huh-7 and HepG2 cells were transfected with biotinylated miR-193a-5p. Then, the cells were rinsed with ice-cold phosphate buffer saline (PBS), harvested, and incubated on ice for 10 min in a lysis buffer (Sigma, USA). Afterward, the lysates were centrifuged at 12,000 g for 30 min at 4°C. The supernatant was incubated with M-280 streptavidin magnetic beads (Sigma). Then, the complex was eluted, the RNA in the complex was extracted, and qRT-PCR was executed to examine RNA levels.

Western blot

Cells were lysed with RIPA buffer (Beyotime, China) and protein was harvested and quantified by a BCA kit (Beyotime). SDS-PAGE was employed to dissolve the protein extractions, then the protein was transferred to PVDF membranes (Sigma-Aldrich). Next, the membranes were blocked in 5% skim milk at room temperature for 1 h. The membranes were incubated with anti-MMP16 antibody (ab73877, 1:500, Abcam, UK), anti-MMP2 antibody (ab86607, 1:500, Abcam), anti-MMP9 antibody (ab38898, 1:500, Abcam), and anti-GAPDH antibody (ab8245, 1:3000, Abcam) for 12 h at 4°C. Next, the membranes were incubated with the secondary antibody (1:5000, Abcam) for 1 h at room temperature. An ECL Plus reagent (Life Technologies, USA) was added onto the membranes and the signals were detected using a chemiluminescence detection device (Bio-Rad).

Statistical analysis

All the experiments were executed at least in triplicate. Statistical analysis was done with SPSS v22.0 (SPSS, Inc., USA). All data are reported as means±SD. To make the comparisons between two groups and among multiple groups, t-test and one-way ANOVA were used, respectively. The correlation was evaluated by Pearson's correlation coefficient. P<0.05 signified statistical significance.

Results

Circ_0001806 was up-modulated in HCC tissues and cell lines

First of all, we analyzed the microarray data of GSE94508 and GSE97332 downloaded from the GEO database, with |log2(FC)|>1.5 and P<0.05 as criteria, and then obtained the differentially expressed circRNAs (DEcircRNAs) in HCC and in normal tissue samples with GEO2R tool (Figure 1A). There were 242 DEcircRNAs in GSE94508 and 439 DEcircRNAs in GSE97332. Fifteen circRNAs were differentially expressed in HCC samples in two datasets (Table 2). For instance, as shown in the heat map, 5 up-regulated circRNAs were displayed, including circ_0001806 (Figure 1B and C). Importantly, RNase R could degrade linear RNA GAPDH rather than circ_0001806, which indicated that circ_0001806 possessed a loop structure (Figure 1D). Furthermore, circ_0001806 was predominantly localized in the cytoplasm of L02 cells, suggesting it could probably function as a ceRNA (Figure 1E). Subsequently, qRT-PCR was adopted to examine circ_0001806 expression in 49 pairs of HCC tissues and non-tumor tissues. The data demonstrated that circ_0001806 expression in HCC tissues was higher than that in normal tissues adjacent to cancer (Figure 1F). We also measured circ_0001806 expression in HCC cell lines (Huh-7, HepG2, SMMC-7721, Bel-7402, and HL-7702) and it was significantly up-modulated compared with L02 cells (Figure 1G). These data revealed that circ_0001806 was overexpressed in HCC. We hypothesized that circ_0001806 might be a vital regulator in HCC progression.

Figure 1
circ_0001806 is up-modulated in hepatocellular carcinoma (HCC). A, 15 circRNAs were differentially expressed in two microarray datasets (GSE94508, GSE97332) containing the circRNA expression profiles of HCC tissues and normal liver tissues. Heat maps showing the up-modulated circRNAs in HCC tissues compared to the matched non-tumor tissues GSE94508 (B) and GSE97332 (C). D, The levels of circ_0001806 were examined by qRT-PCR after treatment with RNase R. E, The subcellular distribution of circ_0001806 was analyzed by qRT-PCR after cellular RNA fractionation. GAPDH was used as the control for cytoplasm and U6 was used as the control for nucleus. F, The expression levels of circ_0001806 in 49 paired HCC tissues and matched adjacent normal tissues were examined by qRT-PCR. G, The expression levels of circ_0001806 in human HCC cell lines (Huh-7, HepG2, SMMC-7721, Bel-7402, and HL-7702) and the normal hepatic cell line (L02) were examined by qRT-PCR. Data are reported as means±SD. ***P<0.001 (ANOVA).

Table 2
The intersected circRNAs form differentially expressed circRNAs in 2 gene sets.

Gene symbol	GSE94508		GSE97332
	LogFC	P value	LogFC	P value
hsa_circ_0067934	1.831013	0.0223	2.6748404	0.00315
hsa_circ_0043438	1.554496	7.38E-09	2.0000682	0.0397
hsa_circ_0014879	1.548521	7.18E-09	1.5874528	0.0302
hsa_circ_0000981	1.541864	1.13E-08	1.6756029	0.0152
hsa_circ_0001806	1.539452	2.25E-09	2.3077778	0.000246
hsa_circ_0079480	1.524299	2.47E-06	-1.714286	4.57E-06
hsa_circ_0009581	-1.547747	0.000224	2.5098621	0.000521
hsa_circ_0084429	-1.566135	0.00311	1.7969988	0.000559
hsa_circ_0027774	-1.595124	0.000121	1.6627344	0.000397
hsa_circ_0056548	-1.759187	0.0000021	-1.66516	0.000954
hsa_circ_0007762	-1.951053	0.0000012	-1.553882	0.00238
hsa_circ_0078279	-2.060631	0.0000212	-1.501581	0.00155
hsa_circ_0059859	-2.064442	0.0000086	1.5601697	0.000159
hsa_circ_0055033	-2.843088	2.48E-06	2.6107158	0.0000107
hsa_circ_0004913	-6.093206	8.96E-13	-3.585576	0.0000169

• Names in red indicate upregulated circRNAs in both gene sets; values in blue indicate the logFC value of circRNAs in both gene sets.

Knockdown circ_0001806 repressed the growth and migration of HCC

Based on the above data, among the HCC cell lines, circ_0001806 expression was relatively higher in Huh-7 and HepG2 cells. Therefore, we transfected Huh-7 and HepG2 cells with 3 siRNAs targeting circ_0001806 (si-circ_0001806#1, si-circ_0001806#2, and si-circ_0001806#3), and si-circ_0001806#1 showed the best knockdown efficiency. Therefore, si-circ_0001806#1 was chosen for the follow-up experiments (Figure 2A). To probe the role of circ_0001806 in HCC, CCK-8 and Transwell methods were utilized to examine the multiplication, migration, and invasion of HCC cells. Relative to the control group, the proliferation of cells in circ_0001806 knockdown group was remarkably restrained (Figure 2B). Additionally, the migration and invasion of HCC cells in circ_0001806 knockdown group were also significantly repressed compared with those of the control group (Figure 2C and D). These data implied that knocking down circ_0001806 suppressed the proliferation, migration, and invasion of HCC cells, suggesting that circ_0001806 promoted HCC progression.

Figure 2
Knockdown of circ_0001806 inhibits multiplication, migration, and invasion of HCC cells. A, Three siRNAs targeting circ_0001806 (si-circ_0001806#1, si-circ_0001806#2, and si-circ_0001806#3) were transfected into Huh-7 and HepG2 cells and the expression level of circ_0001806 was detected by qRT-PCR. B, CCK-8 method was used to detect the multiplication of Huh-7 and HepG2 cells after transfection. C and D, Transwell experiment was employed to detect the migration and invasion of Huh-7 and HepG2 cells after transfection (scale bar 20 μm). Data are reported as means±SD. *P<0.05, **P<0.01, ***P<0.001; ns, not statistically significant (ANOVA and t-test).

Circ_0001806 targeted miR-193a-5p in HCC cells

To determine the downstream mechanism of circ_0001806 in HCC progression, the starBase database (http://starbase.sysu.edu.cn/) was searched, and it was predicted that miR-193a-5p was a potential downstream target of circ_0001806 (Figure 3A, Supplementary Figure S1). To verify the targeting relationship between circ_0001806 and miR-193a-5p, we co-transfected circ_0001806 wild-type plasmid (circ_0001806-WT) or mutant plasmid (circ_0001806-MUT) with miR-193a-5p mimic or control (miR-control) into HEK293T cells, respectively. The luciferase reporter gene assay revealed that, compared with the control group, miR-193a-5p mimics significantly repressed the luciferase activity of cells transfected with circ_0001806-WT plasmid, but exerted no significant effect on the luciferase activity of cells transfected with circ_0001806-MUT plasmid (Figure 3B). RIP and pull-down data revealed that circ_0001806 could bind with miR-193a-5p in both Huh-7 and HepG2 cells (Figure 3C and D). Moreover, after knocking down circ_0001806, miR-193a-5p expression in Huh-7 and HepG2 cells was significantly augmented (Figure 3E). In addition, data revealed that miR-193a-5p was significantly under-expressed in HCC tissues compared with normal tissues adjacent to cancer (Figure 3F). Pearson's correlation analysis showed that miR-193a-5p expression was negatively correlated with circ_0001806 expression (R²=0.4279, P<0.001) (Figure 3G). These data suggested that circ_0001806 adsorbed miR-193a-5p and negatively modulated its expression in HCC cells.

Figure 3
miR-193a-5p is the downstream target of circ_0001806. A, Bioinformatics analysis predicted that the sequence of miR-193a-5p matched the sequences of circ_0001806, and circ_0001806-WT (wild-type) and circ_0001806-MUT (mutant) luciferase reporter gene vectors were constructed. B, miR-193a-5p mimics (or miR control) were co-transfected with circ_0001806-WT (or circ_0001806-MUT) into HEK293T cells. After 48 h, the luciferase activity of the cells in each group was measured. C, The complex containing circ_0001806 and miR-193a-5p in Huh-7 and HepG2 cells was immunoprecipitated by anti-AGO2 using RIP assay. D, RNA pull-down assay was carried out to verify the interactions between circ_0001806 and miR-193a-5p. E, qRT-PCR was used to detect the expression of miR-193a-5p in Huh-7 and HepG2 cells after circ_0001806 knockdown. F, qRT-PCR was used to examine the expression levels of miR-193a-5p in 49 paired hepatocellular carcinoma (HCC) tissues and matched adjacent normal tissues. G, Pearson correlation analysis was utilized to evaluate the correlation between the expressions of circ_0001806 and miR-193a-5p in HCC tissues. Data are reported as means±SD. ***P<0.001 (t-test).

MMP16 was one of the direct targets of miR-193a-5p in HCC

Next, bioinformatics analysis was conducted with starBase, TargetScan (www.targetscan.org), and miRDB (www.mirdb.org) to identify candidate mRNAs targeted by miR-193a-5p, and analysis revealed that MMP16 3′UTR and miR-193a-5p had complementary binding sites (Figure 4A and B). GEPIA database (http://gepia.cancer-pku.cn/) indicated that the overall survival of patients with high MMP16 expression was significantly shorter than that of patients with low MMP16 expression (Supplementary Figure S2). The data of the luciferase reporter gene experiment showed that the luciferase activity of the cells co-transfected with the luciferase reporter plasmid MMP16-WT and miR-193a-5p mimic was significantly diminished, while the luciferase activity of the cells transfected with MMP16-MUT was unchanged (Figure 4C). Furthermore, anti-miR-193a-5p was transfected into Huh-7 and HepG2 cell lines to construct an miR-193a-5p under-expression model (Figure 4D). qRT-PCR showed that after knocking down circ_0001806, MMP16 mRNA expression in Huh-7 and HepG2 cells was significantly decreased, while after miR-193a-5p expression was inhibited, and MMP16 mRNA expression in Huh-7 and HepG2 cells was remarkably augmented (Figure 4E and F). Additionally, western blot data were consistent with the results of qRT-PCR (Figure 4G and H). The findings indicated that MMP16 was a direct target of miR-193a-5p, and circ_0001806 could positively regulate MMP16 expression in HCC cells.

Figure 4
MMP16 is the direct target of miR-193a-5p. A, Bioinformatics analysis was performed with starBase, TargetScan, and miRDB to predict the candidate mRNAs, which could be targeted by miR-193a-5p. B, Bioinformatics analysis predicted that the sequence of MMP16 3UTR matched the sequences of miR-193a-5p, and MMP16-WT (wild-type) and MMP16-MUT (mutant) luciferase reporter gene vectors were constructed. C, The miR-193a-5p mimics (or miR control) were co-transfected with MMP16-WT (or MMP16-MUT) into HEK293T cells, and the luciferase activity of the cells in each group was measured. D, Huh-7 and HepG2 cells were transfected with anti-miR-193a-5p and negative control (NC), and the expression level of miR-193a-5p was detected by qRT-PCR. E, qRT-PCR was used to detect the expression of MMP16 mRNA in Huh-7 and HepG2 cells after circ_0001806 knockdown. F, qRT-PCR was used to detect the expression of MMP16 mRNA in Huh-7 and HepG2 cells transfected with anti-miR-193a-5p. G, Western blot was used to detect the expression of MMP16 protein in Huh-7 and HepG2 cells after circ_0001806 knockdown (G) and in Huh-7 and HepG2 cells transfected with anti-miR-193a-5p (H). Data are reported as means±SD. **P<0.01, ***P<0.001 (t-test).

Circ_0001806 enhanced the multiplication, migration, and invasion of HCC cells through the miR-193a-5p/MMP16 axis

To assess whether circ_0001806 affected HCC progression by modulating the miR-193a-5p/MMP16 axis, compensation experiments were performed. Anti-miR-193a-5p or negative control (anti-miR-control) was co-transfected into Huh-7 and HepG2 cell lines transfected with si-circ_0001806#1 (Figure 5A). Inhibition of miR-193a-5p partly counteracted the inhibitory effect of circ_0001806 knockdown on MMP16 expression in Huh-7 and HepG2 cells (Figure 5B and C). Western blot results also showed that circ_0001806 knockdown could inhibit the expression of MMP2 and MMP9, but inhibition of miR-193a-5p partly counteracted this effect (Figure 5C). CCK-8 experiment data revealed that knocking down circ_0001806 repressed the multiplication of Huh-7 and HepG2 cells, while down-modulating miR-193a-5p partly counteracted this effect (Figure 5D). Additionally, the data of the Transwell experiment suggested that knocking down circ_0001806 impeded the migration and invasion of Huh-7 and HepG2 cells, while this effect could be reversed by the inhibition of miR-193a-5p (Figure 5E and F). Collectively, it was concluded that circ_0001806 enhanced the multiplication, migration, and invasion of HCC cells via the miR-193a-5p/MMP16 axis.

Figure 5
Circ_0001806 enhances the multiplication and invasion of hepatocellular carcinoma (HCC) cells through miR-193a-5p/MMP16 axis. Huh-7 and HepG2 cells were transfected with si-NC + anti-miR-Con, si-circ_0001806#1 + anti-miR-Con, or si-circ_0001806#1 + anti-miR-193a-5p. A, qRT-PCR was used to detect the expression of miR-193a-5p in Huh-7 and HepG2 cells. B, qRT-PCR was used to detect the expression of MMP16 mRNA in Huh-7 and HepG2 cells. C, Western blot was used to detect the expression of MMP16, MMP2, and MMP9 protein in Huh-7 and HepG2 cells. D, CCK-8 method was used to detect the multiplication of Huh-7 and HepG2 cells after transfection. The transwell experiment was employed to detect the migration and invasion of Huh-7 cells (E) and HepG2 cells (F) after transfection (scale bar 20 μm). Data are reported as means±SD. *P<0.05, **P<0.01, and ***P<0.001 (ANOVA).

Discussion

HCC is one of the most deadly malignancies (¹²,¹³). Clarifying the mechanism of HCC progression is crucial to develop novel biomarkers and therapy targets for this disease (¹⁴,¹⁵). Several circRNAs are known to play roles in the tumorigenesis and disease progression of HCC (¹⁶- ¹⁹). For instance, circ-ADD3 is under-expressed in HCC and it represses the metastasis of HCC by modulating EZH2 (¹⁷). Circ_0000517 is up-modulated in HCC and it contributes to HCC progression by modulating MAPK and Ras pathways (¹⁸). The expression of circ-PRKCI is associated with poor prognosis in HCC patients, and circ-PRKCI could inhibit cell apoptosis and promote cell invasion (¹⁹). Furthermore, knocking down circ_0000092 restrains HCC cell multiplication, migration, invasion, and angiogenesis by down-modulating HN1 (²⁰). In this work, we focused on the expression characteristics, biological function, and underlying mechanism of circ_0001806 in HCC. Circ_0001806 is reported to be remarkably up-modulated in colorectal cancer tissues and it is linked to TNM staging, depth of invasion, lymphatic metastasis, and distant metastasis (²¹). In this work, we demonstrated that circ_0001806 was also overexpressed in HCC tissues and cell lines. Additionally, circ_0001806 knockdown repressed the multiplication, migration and invasion of HCC cells. These findings indicated that circ_0001806 facilitated HCC progression.

Some circRNAs work as ceRNAs to regulate miRNAs, and this mechanism is common in cancer biology (²²,²³). For instance, circ_0001649 can function as ceRNA to absorb miR-127-5p, miR-612, and miR-4688 and activate SHPRH to repress HCC progression (²⁴). Down-regulating circ_0005986 enhances cell multiplication by increasing miR-129-5p and decreasing Notch1 expression to facilitate cell cycle progression of HCC cells (²⁵). Reportedly, miR-193a-5p is under-expressed in HCC, and miR-193a-5p suppresses the malignancy of HCC cells (²⁶,²⁷). This suggests that miR-193a-5p is a tumor suppressor in HCC. In this work, circ_0001806 was found to target miR-193a-5p and repress its expression level in HCC cells, which suggested that the abnormal high expression of circ_0001806 contributed to the under-expression of miR-193a-5p in HCC. Also, we demonstrated that inhibiting miR-193a-5p partially counteracted the effects of circ_0001806 knockdown on HCC cells, and this finding indicated that circ_0001806 partly exerts its biological function via regulating miR-193a-5p in HCC cells.

MMP16 is a membrane-type metalloprotease that facilitates cancer progression via multiple mechanisms (²⁸,²⁹). MMP16 can degrade the components of extracellular matrix and facilitate the migration and invasion of cancer cells (²⁹). Reportedly, MMP16 is abnormally expressed in diverse cancers such as melanoma, lung cancer, endometrial cancer, prostate cancer, and others (²⁹- ³³). For instance, MMP16 is overexpressed in gastric cancer and is linked to unfavorable prognosis (³³). Moreover, MMP16 is validated as the target gene of diverse miRNAs, such as miR-132, miR-328-3p, miR-33b, and miR-26a (³⁴- ³⁷). These studies suggest that miRNAs are crucial regulators of MMP16 expression. In HCC, it is reported that MMP16 overexpression enhances HCC metastasis by modulating epithelial-mesenchymal transition process (³⁸). In this work, we demonstrated that MMP16 was identified as a target gene of miR-193a-5p, and circ_0001806 positively modulated MMP16 expression by repressing miR-193a-5p. These data support that the ceRNA network consisting of circ_0001806, miR-193a-5p, and MMP16 is a novel mechanism of HCC progression. MMP9 and MMP2 are also crucial regulators of migration and invasion in HCC (³⁹,⁴⁰). Interestingly, we also demonstrated that circ_0001806 / miR-193a-5p axis not only regulated the expression of MMP16 in HCC cells, but also regulated the expression of MMP2 and MMP9. Our data suggested that circ_0001806 and miR-193a-5p are important regulators of matrix metalloprotease members, but the detailed mechanisms remain to be deciphered in following studies.

In summary, our research showed that circ_0001806 was up-modulated in HCC tissues and cells. Circ_0001806 knockdown repressed the multiplication and metastatic potential of HCC cells. Mechanistically, circ_0001806 participated in HCC progression through the miR-193a-5p/MMP16 axis. Our data suggested that circ_0001806 is a promising potential biomarker and therapy target for HCC.

References

• ¹
  Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 2018; 68: 394-424, doi: 10.3322/caac.21492.
  » https://doi.org/10.3322/caac.21492
• ²
  Lee M, Ko H, Yun M. Cancer metabolism as a mechanism of treatment resistance and potential therapeutic target in hepatocellular carcinoma. Yonsei Med J 2018; 59: 1143-1149, doi: 10.3349/ymj.2018.59.10.1143.
  » https://doi.org/10.3349/ymj.2018.59.10.1143
• ³
  Li Y, Li G, Chen X, Huang H, Liao L, Yuan T, et al. A novel lncRNA NONHSAT053785 acts as an independent risk factor for intrahepatic metastasis of hepatocellular carcinoma. Onco Targets Ther 2020; 13: 5455-5466, doi: 10.2147/OTT.S254455.
  » https://doi.org/10.2147/OTT.S254455
• ⁴
  Fornari F, Gramantieri L, Callegari E, Shankaraiah RC, Piscaglia F, Negrini M, et al. MicroRNAs in animal models of HCC. Cancers (Basel) 2019; 11: 1906, doi: 10.3390/cancers11121906.
  » https://doi.org/10.3390/cancers11121906
• ⁵
  Zhang Z, Yang T, Xiao J. Circular RNAs: promising biomarkers for human diseases. EBioMedicine 2018; 34: 267-274, doi: 10.1016/j.ebiom.2018.07.036.
  » https://doi.org/10.1016/j.ebiom.2018.07.036
• ⁶
  Li J, Yang J, Zhou P, Le Y, Zhou C, Wang S, et al. Circular RNAs in cancer: novel insights into origins, properties, functions and implications. Am J Cancer Res 2015; 5: 472-480.
• ⁷
  Liu J, Li D, Luo H, Zhu X. Circular RNAs: the star molecules in cancer. Mol Aspects Med 2019; 70: 141-152, doi: 10.1016/j.mam.2019.10.006.
  » https://doi.org/10.1016/j.mam.2019.10.006
• ⁸
  Zong L, Sun Q, Zhang H, Chen Z, Deng Y, Li D, et al. Increased expression of circRNA_102231 in lung cancer and its clinical significance. Biomed Pharmacother 2018; 102: 639-644, doi: 10.1016/j.biopha.2018.03.084.
  » https://doi.org/10.1016/j.biopha.2018.03.084
• ⁹
  Bi W, Huang J, Nie C, Liu B, He G, Han J, et al. CircRNA circRNA_102171 promotes papillary thyroid cancer progression through modulating CTNNBIP1-dependent activation of β-catenin pathway. J Exp Clin Cancer Res 2018; 37: 275, doi: 10.1186/s13046-018-0936-7.
  » https://doi.org/10.1186/s13046-018-0936-7
• ¹⁰
  Jiang Z, Shen L, Wang S, Wu S, Hu Y, Guo J, et al. Hsa_circ_0028502 and hsa_circ_0076251 are potential novel biomarkers for hepatocellular carcinoma. Cancer Med 2019; 8: 7278-7287, doi: 10.1002/cam4.2584.
  » https://doi.org/10.1002/cam4.2584
• ¹¹
  Lu Q, Liu T, Feng H, Yang R, Zhao X, Chen W, et al. Circular RNA circSLC8A1 acts as a sponge of miR-130b/miR-494 in suppressing bladder cancer progression via regulating PTEN. Mol Cancer 2019; 18: 111, doi: 10.1186/s12943-019-1040-0.
  » https://doi.org/10.1186/s12943-019-1040-0
• ¹²
  Tang A, Hallouch O, Chernyak V, Kamaya A, Sirlin CB. Epidemiology of hepatocellular carcinoma: target population for surveillance and diagnosis. Abdom Radiol (NY) 2018; 43: 13-25.
• ¹³
  Yang JD, Hainaut P, Gores GJ, Amadou A, Plymoth A, Roberts LR. A global view of hepatocellular carcinoma: trends, risk, prevention and management. Nat Rev Gastroenterol Hepatol 2019; 16: 589-604, doi: 10.1038/s41575-019-0186-y.
  » https://doi.org/10.1038/s41575-019-0186-y
• ¹⁴
  Schulze K, Imbeaud S, Letouzé E, Alexandrov LB, Calderaro J, Rebouissou S, et al. Exome sequencing of hepatocellular carcinomas identifies new mutational signatures and potential therapeutic targets. Nat Genet 2015; 47: 505-511, doi: 10.1038/ng.3252.
  » https://doi.org/10.1038/ng.3252
• ¹⁵
  Harding JJ, Nandakumar S, Armenia J, Khalil DN, Albano M, Ly M, et al. Prospective genotyping of hepatocellular carcinoma: clinical implications of next-generation sequencing for matching patients to targeted and immune therapies. Clin Cancer Res 2019; 25: 2116-2126, doi: 10.1158/1078-0432.CCR-18-2293.
  » https://doi.org/10.1158/1078-0432.CCR-18-2293
• ¹⁶
  Fu L, Jiang Z, Li T, Hu Y, Guo J. Circular RNAs in hepatocellular carcinoma: Functions and implications. Cancer Med 2018; 7: 3101-3109, doi: 10.1002/cam4.1574.
  » https://doi.org/10.1002/cam4.1574
• ¹⁷
  Sun S, Wang W, Luo X, Li Y, Liu B, Li X, et al. Circular RNA circ-ADD3 inhibits hepatocellular carcinoma metastasis through facilitating EZH2 degradation via CDK1-mediated ubiquitination. Am J Cancer Res 2019; 9: 1695-1707.
• ¹⁸
  Wang X, Wang X, Li W, Zhang Q, Chen J, Chen T. Up-regulation of hsa_circ_0000517 predicts adverse prognosis of hepatocellular carcinoma. Front Oncol 2019; 9: 1105, doi: 10.3389/fonc.2019.01105.
  » https://doi.org/10.3389/fonc.2019.01105
• ¹⁹
  Qi SX, Sun H, Liu H, Yu J, Jiang ZY, Yan P. Role and mechanism of circ-PRKCI in hepatocellular carcinoma. World J Gastroenterol 2019; 25: 1964-1974, doi: 10.3748/wjg.v25.i16.1964.
  » https://doi.org/10.3748/wjg.v25.i16.1964
• ²⁰
  Pu J, Wang J, Li W, Lu Y, Wu X, Long X, et al. hsa_circ_0000092 promotes hepatocellular carcinoma progression through up-regulating HN1 expression by binding to microRNA-338-3p. J Cell Mol Med 2020, doi: 10.1111/jcmm.15010.
  » https://doi.org/10.1111/jcmm.15010
• ²¹
  Sun J, Liu J, Zhu Q, Xu F, Kang L, Shi X. Hsa_circ_0001806 acts as a ceRNA to facilitate the stemness of colorectal cancer cells by increasing COL1A1. Onco Targets Ther 2020; 13: 6315-6327, doi: 10.2147/OTT.S255485.
  » https://doi.org/10.2147/OTT.S255485
• ²²
  Jeck WR, Sharpless NE. Detecting and characterizing circular RNAs. Nat Biotechnol 2014; 32: 453-461, doi: 10.1038/nbt.2890.
  » https://doi.org/10.1038/nbt.2890
• ²³
  Arnaiz E, Sole C, Manterola L, Iparraguirre L, Otaegui D, Lawrie CH. CircRNAs and cancer: biomarkers and master regulators. Semin Cancer Biol 2019; 58: 90-99, doi: 10.1016/j.semcancer.2018.12.002.
  » https://doi.org/10.1016/j.semcancer.2018.12.002
• ²⁴
  Su Y, Xu C, Liu Y, Hu Y, Wu H. Circular RNA hsa_circ_0001649 inhibits hepatocellular carcinoma progression via multiple miRNAs sponge. Aging (Albany NY) 2019; 11: 3362-3375, doi: 10.18632/aging.101988.
  » https://doi.org/10.18632/aging.101988
• ²⁵
  Fu L, Chen Q, Yao T, Li T, Ying S, Hu Y, et al. Hsa_circ_0005986 inhibits carcinogenesis by acting as a miR-129-5p sponge and is used as a novel biomarker for hepatocellular carcinoma. Oncotarget 2017; 8: 43878-43888, doi: 10.18632/oncotarget.16709.
  » https://doi.org/10.18632/oncotarget.16709
• ²⁶
  Li P, Xiao Z, Luo J, Zhang Y, Lin L. MiR-139-5p, miR-940 and miR-193a-5p inhibit the growth of hepatocellular carcinoma by targeting SPOCK1. J Cell Mol Med 2019; 23: 2475-2488, doi: 10.1111/jcmm.14121.
  » https://doi.org/10.1111/jcmm.14121
• ²⁷
  Wang JT, Wang ZH. Role of miR-193a-5p in the multiplication and apoptosis of hepatocellular carcinoma. Eur Rev Med Pharmacol Sci 2018; 22: 7233-7239.
• ²⁸
  Sounni NE, Noel A. Membrane type-matrix metalloproteinases and tumor progression. Biochimie 2005; 87: 329-342, doi: 10.1016/j.biochi.2004.07.012.
  » https://doi.org/10.1016/j.biochi.2004.07.012
• ²⁹
  Tatti O, Gucciardo E, Pekkonen P, Holopainen T, Louhimo R, Repo P, et al. MMP16 mediates a proteolytic switch to promote cell-cell adhesion, collagen alignment, and lymphatic invasion in melanoma. Cancer Res 2015; 75: 2083-2094, doi: 10.1158/0008-5472.CAN-14-1923.
  » https://doi.org/10.1158/0008-5472.CAN-14-1923
• ³⁰
  Hadchouel A, Decobert F, Franco-Montoya ML, Halphen I, Jarreau PH, Boucherat O, et al. Matrix metalloproteinase gene polymorphisms and bronchopulmonary dysplasia: identification of MMP16 as a new player in lung development. PLoS One 2008; 3: e3188, doi: 10.1371/journal.pone.0003188.
  » https://doi.org/10.1371/journal.pone.0003188
• ³¹
  Rak B, Mehlich D, Garbicz F, Domosud Z, Paskal W, Marczewska JM, et al. Post-transcriptional regulation of MMP16 and TIMP2 expression via miR-382, miR-410 and miR-200b in endometrial cancer. Cancer Genomics Proteomics 2017; 14: 389-401.
• ³²
  Lin HY, Cheng CH, Chen DT, Chen YA, Park JY. Coexpression and expression quantitative trait loci analyses of the angiogenesis gene-gene interaction network in prostate cancer. Transl Cancer Res 2016; 5: S951-S963, doi: 10.21037/tcr.2016.10.55.
  » https://doi.org/10.21037/tcr.2016.10.55
• ³³
  Cao L, Chen C, Zhu H, Gu X, Deng D, Tian X, et al. MMP16 is a marker of poor prognosis in gastric cancer promoting multiplication and invasion. Oncotarget 2016; 7: 51865-51874, doi: 10.18632/oncotarget.10177.
  » https://doi.org/10.18632/oncotarget.10177
• ³⁴
  Wang H, Li XT, Wu C, Wu ZW, Li YY, Yang TQ, et al. miR-132 can inhibit glioma cells invasion and migration by target MMP16 in vitro Onco Targets Ther 2015; 8: 3211-3218, doi: 10.2147/OTT.S79282.
  » https://doi.org/10.2147/OTT.S79282
• ³⁵
  Zhang M, Zhang J, Zhou Q. Elevated expression of microRNA-328-3p suppresses aggressive malignant behaviors via targeting matrix metalloprotease 16 in osteosarcoma. Onco Targets Ther 2019; 12: 2063-2070, doi: 10.2147/OTT.S195022.
  » https://doi.org/10.2147/OTT.S195022
• ³⁶
  Luo Y, Wang Q, Teng L, Zhang J, Song J, Bo W, et al. LncRNA DANCR promotes multiplication and metastasis in pancreatic cancer by regulating miRNA-33b. FEBS Open Bio 2020; 10: 18-27, doi: 10.1002/2211-5463.12732.
  » https://doi.org/10.1002/2211-5463.12732
• ³⁷
  Zheng W, Li ZY, Zhao DL, Li XL, Liu R. microRNA-26a directly targeting MMP14 and MMP16 inhibits the cancer cell multiplication, migration and invasion in cutaneous squamous cell carcinoma. Cancer Manag Res 2020; 12: 7087-7095, doi: 10.2147/CMAR.S265775.
  » https://doi.org/10.2147/CMAR.S265775
• ³⁸
  Shen Z, Wang X, Yu X, Zhang Y, Qin L. MMP16 promotes tumor metastasis and indicates poor prognosis in hepatocellular carcinoma. Oncotarget 2017; 8: 72197-72204, doi: 10.18632/oncotarget.20060.
  » https://doi.org/10.18632/oncotarget.20060
• ³⁹
  Alaseem A, Alhazzani K, Dondapati P, Alobid S, Bishayee A, Rathinavelu A. Matrix metalloproteinases: a challenging paradigm of cancer management. Semin Cancer Biol 2019; 56: 100-115, doi: 10.1016/j.semcancer.2017.11.008.
  » https://doi.org/10.1016/j.semcancer.2017.11.008
• ⁴⁰
  Lin XL, Li K, Yang Z, Chen B, Zhang T. Dulcitol suppresses proliferation and migration of hepatocellular carcinoma via regulating SIRT1/p53 pathway. Phytomedicine 2020; 66: 153112, doi: 10.1016/j.phymed.2019.153112.
  » https://doi.org/10.1016/j.phymed.2019.153112

Supplementary Material

Click here to view [pdf].

Publication Dates

• Publication in this collection
  29 Oct 2021
• Date of issue
  2021

History

• Received
  17 May 2021
• Accepted
  2 Sept 2021