Emerging Roles Of hsa-circ-0046600 Targeting The miR-640/HIF-1α Signalling Pathway In The Progression Of HCC

Introduction

Hepatocellular carcinoma (HCC) is a malignant tumour with a high incidence and risk of death worldwide and is related to obesity, alcohol abuse, dyslipidaemia, metabolic syndrome and hepatitis.^1,2 Currently, although several studies have shown that many genes are involved in HCC and many treatment strategies have been developed for this disease, the 5-year survival rate of patients is still very low.^3,4 Therefore, it is still necessary to explore the mechanism and therapeutic targets in further research on HCC.

Circular RNAs (circRNAs) are non-coding RNAs with a length of approximately 100 bp to 4 kb and are found in almost all major eukaryotes (such as animals, plants and fungi), including in the blood, different subcellular compartments (including nuclei and cytoplasms), and extracellular exosomes.⁵ At present, most studies focus on the mode of action of circRNAs (which is similar to that of long non-coding RNAs (lncRNAs)), as microRNA (miRNA) sponges affect the negative regulatory effect of miRNAs on mRNAs, thereby regulating gene expression.⁶

In this study, we discovered a functionally unknown circRNA, hsa-circ-0046600, which is associated with HCC, by using the bioinformatics software circBase (http://www.circbase.org/).⁷ hsa-circ-00446600 is a circRNA consisting of the anti-shear of exons 4, 5, and 6 of its maternal gene B3GNTL1 and has a length of 227 bp. We found that hsa-circ-0046600 is highly expressed in HCC and is involved in regulating the migration of HCC cells. As hsa-circ-0046600 is mainly expressed in the cytoplasm, we mainly studied the competitive endogenous RNA (ceRNA) mechanism of hsa-circ-0046600. Studies have shown that hsa-circ-004660 regulates the functions of its target genes, such as HIF, by acting as a sponge for miR-640. These results provide an important basis for the use of hsa-circ-0046600 as a potential target for the diagnosis and treatment of HCC.

Materials And Methods

Human Tissue Samples And Cell Culture

A total of 20 paired HCC tissue and adjacent normal tissue samples were obtained from patients who underwent surgical procedures at HENAN PROVINCIAL PEOPLE’S HOSPITAL of Zhengzhou University. All of the patients provided written consent, and approval was obtained from the Ethics Committee of HENAN PROVINCIAL PEOPLE’S HOSPITAL of Zhengzhou University. The patient consent was written informed consent, and that this was conducted in accordance with the Declaration of Helsinki. The HL-7702 (L-02), SK-HEP-1, SMMC-7721, HepG2, HuH7, and 293T cell lines were obtained from the Cell Bank of Chinese Academy of Sciences (Shanghai, China). HepG2, SMMC-7721, SK-HEP-1, and HuH7 cells were cultured in Dulbecco’s modified Eagle medium (Gibco, USA), and L-02 and SMMC-7721 cells were cultured in RPMI 1640 medium (Gibco, USA) supplemented with 10% foetal bovine serum (Gibco, USA), 50 U/mL penicillin, and 50 μg/mL streptomycin (Gibco, USA). The cells were maintained at 37°C in a humidified incubator at 5% CO₂.

Small Interfering RNA (siRNA) And miRNA Transfection

hsa-circ-0046600 siRNA and miR-640 mimics were all purchased from RiboBio (Guangzhou, China). SK-HEP-1 and HepG2 cells were seeded 24 h before transfection. The siRNA and miRNA mimics were transfected using Lipofectamine 3000 (Invitrogen, USA) with serum-free medium. After 48 h, the cells were collected for testing.

Quantitative Reverse Transcription PCR (qRT-PCR)

Total RNA was extracted from cell lines or the 20 paired HCC tissue and adjacent normal tissue samples using TRIzol Reagent (Invitrogen, USA). A total of 1 µg of total RNA was reverse transcribed into cDNA using a Thermo Fisher Scientific RevertAid First Star and cDNA Synthesis Kit (Thermo Fisher Scientific, USA) (the circRNA primer was a random hexamer primer; the miRNA primer was synthesized by RIBO). Then, qRT-PCR amplification of the cDNA was performed using PowerUp^TM SYBR^TM Green Master Mix (Applied Biosystems, USA). The sequences of the miRNA primers synthesized by RiboBio and of the circRNA and mRNA primers were as follows:

hsa-circ-0046600, 5`-GCTCCTAACCCAGTCGGATAC-3` (forward), 5`-CCGAATCCAAAAAGCAAAGGTA-3` (reverse);

B3GNTL1 mRNA, 5`-TCGGATACGCTAAAAATCAAG-3` (forward), 5`-CTGGGTTAGGAGCTGCTCCG-3` (reverse); and

β-actin, 5`-CCTGTACGCCAACACAGTGC-3` (forward), 5`-ATACTCCTGCTTGCTGATCC-3` (reverse).

Fluorescence In Situ Hybridization (FISH)

Hybridization was performed overnight with hsa-circ-0046600 probes, sections were washed, and DAPI stain was added and incubated for 8 min in the dark. The sections were observed under a Nikon upright fluorescence microscope, and images were acquired. The nuclei stained with DAPI appeared blue under UV excitation and contained hsa-circ-0046600, which was positively labelled with Cy3 (red emission). The sequence of the hsa-circ-0046600 probe used for FISH was 5ʹ-TTTAG CGTAT CCGAC TGGGT TAGGA GCTGC-3ʹ.

Wound Scratch Assay

SK-HEP-1 and HepG2 cells were seeded into 6-well plates and transfected on the second day. Cells were then incubated at 5% CO₂/37°C until they fused into a monolayer, and the cell monolayer was scraped along a straight line with a 200 µl pipette tip to create a gap. The migration distance was calculated using an inverted microscope at 0 h, 24 h, and 48 h.

Western Blot Analysis

When cells reached 80% confluence, the medium was discarded, and the cells were washed with pre-chilled 1× phosphate-buffered saline (PBS). Then, 320 μL of cell lysis buffer (RIPA buffer supplemented with 3.2 µL PMSF) was added to the cells to extract cellular protein. After a 30-min incubation on ice, the cells were scraped into a 1.5-mL centrifuge tube and subjected to centrifugation at 4°C for 15 min at 12,000 rpm. Protein quantification was performed using a NanoDrop ND-1000. Proteins were separated by sodium dodecyl sulfate polyacrylamide gel electrophoresis and electrotransferred to a polyvinylidene fluoride membrane. The membrane was incubated with a primary antibody at 4°C overnight, washed extensively with 1% TBST, and incubated with a secondary antibody conjugated to horseradish peroxidase (1:1000) at room temperature for 3 h. Immunolabelling was visualized using an electrochemiluminescence system. The primary antibodies included anti-N-cadherin (dilution 1:1000, AF4039) and anti-HIF-1α (dilution 1:200, AF1009) antibodies, which were purchased from Affinity (USA), and an anti-GAPDH antibody (dilution 1:2000, AF7021).

Acquisition Of hsa-circ-0046600, miR-640 And Target mRNAs

The GSE97332 dataset in the Gene Expression Omnibus (GEO) database was used to screen for differentially expressed circRNAs in HCC.⁸ R software was used to screen circRNAs with |logFC |> 2 and adjusted P-value <0.05. CircInteractome (https://circinteractome.nia.nih.gov/), TargetScan 7.2 (http://www.targetscan.org/vert_72/), miRDB (http://www.mirdb.org/) and mirDIP (http://ophid.utoronto.ca/mirDIP/index_confirm.jsp) were used to predict the assembly sites of circRNA-miRNA or miRNA-mRNA.^9–11 For the meta-analysis, GEO datasets were screened to verify the expression of miR-640 in HCC (dataset screening keywords: HCC and (microRNAs OR miRNA) and Homo sapiens). R software was used to screen mRNAs with |logFC |> 1 and adjusted P-value <0.05 in The Cancer Genome Atlas (TCGA)-HCC data.¹²

Functional Enrichment Analysis And Protein-Protein Interaction (PPI) Network Construction

To understand the underlying biological processes and pathways of the aberrantly expressed mRNAs, Gene Ontology (GO) biological process term and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were conducted using the Cytoscape plug-in ClueGO and Metascape (Database for Annotation, Visualization, and Integrated Discovery; http://metascape.org).¹³ The GO enrichment analysis was based on a P-value < 0.05 and a kappa score of 0.3, which were applied as the thresholds to divide the tree into term clusters. Significant enrichment results were visualized using Cytoscape software 3.3.0. The STRING database (http://string-db.org/)¹⁴ was used as background data to construct a PPI network using proteins with an interaction score greater than 0.4.

Statistical Analysis

Statistical analysis was performed using SPSS (version 25.0), and data are recorded as the mean ± standard deviation (SD). Differences between the two groups were analysed using Student’s t-test. Differences between the three groups were analysed using one-way analysis of variance (ANOVA). Pearson’s correlation analysis was used to study the correlation between circRNAs and miRNAs. P < 0.05 was considered statistically significant.

Results

hsa-circ-0046600 Expression Is Upregulated In HCC And Associated With Clinicopathological Variables

The GSE97332 dataset showed that hsa-circ-0046600 was highly expressed in HCC tissues (Figure 1A). The qRT-PCR results showed that the expression level of hsa-circ-0046600 was upregulated by approximately 1.35–4.04-fold in the tumour cell line compared to that in the normal liver cell line (L-02) (Figure 1B). hsa-circ-0046600 was highly expressed in 90% of the HCC tissue samples relative to the corresponding adjacent tissue samples (P = 0.003; Figure 1C and D). In addition, the expression of hsa-circ-0046600 was significantly upregulated in patients with tumours ≥5 cm (P = 0.003) and stage III/IV (P = 0.028). High hsa-circ-0046600 expression was significantly associated with pathological vascular invasion (P = 0.019). However, tumour differentiation, patient age and sex were irrelevant (all P > 0.05, Table 1). We used FISH to locate the expression of hsa-circ-0046600 in HCC and adjacent tissue samples and randomly selected 5 fields of view to count. The results showed that the expression in the tumour group was approximately 1.83 times higher than that in the adjacent normal tissue group.

Table 1 Correlation Between hsa-circ-0046600 Expression And The Clinicopathological Features Of HCC Patients

Figure 1 hsa-circ-0046600 expression in HCC. (A) The top 20 differentially expressed circRNAs in the GSE97332 dataset. In circBase, hsa-circ-102272 is named hsa-circ-0046600. (B) The expression levels of hsa-circ-0046600 in HepG2, SK-HEP-1, HuH7, and SMMC-7721 cells were compared with the expression level of hsa-circ-0046600 in L-02 cells: the fold change was 4.04, 1.88, 1.35, and 1.38, respectively. (C) The expression level of hsa-circ-0046600 was upregulated in 20 HCC tissue samples compared with matched noncancerous tissue samples (fold change≈1.7). (D) hsa-circ-0046600 expression was upregulated in 90% (18/20) of liver tissue samples. (E) The localization of hsa-circ-0046600 was assessed in HCC and adjacent tissues. The results showed that hsa-circ-0046600 was highly expressed in tumour tissue (fold change ≈1.9). Data are reported as the mean ± SD. *P-value < 0.05, **P-value < 0.01, ***P-value < 0.001.

hsa-circ-0046600 Promotes The Migration Of HCC Cells

To reveal the effect of hsa-circ-0046600 in HCC, SK-HEP-1 and HepG2 cells were transfected with an hsa-circ-0046600 siRNA inhibitor or the corresponding negative control (NC). The qRT-PCR results indicated that the siRNA affected only hsa-circ-0046600 and not B3GNTL1 mRNA (Figure 2A). A wound scratch test showed that hsa-circ-0046600 knockdown inhibited the migration of SK-HEP-1 cells (siCirc0046600-A group, the migration distance was 40.33% relative to the NC group, P=0.002; siCirc0046600-B group, the migration distance was 42.33% relative to the NC group, P=0.002) and HepG2 cells (siCirc0046600-A group, the migration distance was 67.33% relative to the NC group, P=0.001; siCirc0046600-B group, the migration distance was 88.25% relative to the NC group, P=0.045) (Figure 2B–E). In SK-HEP-1 and HepG2 cells, Western blot analysis showed that the relative expression level of N-cadherin was significantly lower in the siCirc0046600-A and siCirc0046600-B groups than in the NC group (P<0.05), indicating that hsa-circ-0046600 expression is positively correlated with N-cadherin expression (Figure 2F–H). GAPDH was used as an internal control. Taken together, these results indicate that hsa-circ-0046600 is overexpressed in HCC cells and promotes the migration of these cells.

Figure 2 The role of hsa-circ-0046600 in HCC migration. (A) SiCirc0046600 had a specific effect on hsa-circ-0046600, and the corresponding mRNA was not affected. (B) The migration rates of SK-HEP-1 cells were determined by a wound scratch assay. The migratory ability of the experimental group was lower than that of the control group. (C) Quantification of relative migration in SK-HEP-1 cells transfected with NC, siCirc0046600-A, siCirc0046600-B. (D) The migration rates of HepG2 cells were determined by a wound scratch assay. The migratory ability of the experimental group was lower than that of the control group. (E) Quantification of relative migration in HepG2 cells transfected with NC, siCirc0046600-A, siCirc0046600-B. (F) Western blot analysis showed that N-cadherin and HIF-1α expression was significantly reduced in the experimental group. (G, H) Quantification of relative protein in HepG2 and SK-HEP-1 cells transfected with NC, siCirc0046600-A, siCirc0046600-B. Data are reported as the mean ± standard deviation. *P-value < 0.05, **P-value < 0.01, ***P-value < 0.001.

Features Of hsa-circ-0046600

To confirm the circular characteristics of hsa-circ-0046600, we conducted experiments with 3 different approaches. Sanger sequencing confirmed that hsa-circ-0046600 is a head-to-tail splice product of exon 4 and exon 6 of the PCR product (Figure 3A). Next, we used actinomycin D to inhibit the transcription of hsa-circ-0046600 and B3GNTL1 mRNA in HepG2 cells and then measured the half-lives of these molecules. The results showed that hsa-circ-0046600 was more stable than B3GNTL1 mRNA (Figure 3B). Unlike mRNAs, circRNAs do not have a poly-A tail; therefore, we used random hexamer or oligo(dT) 18 primers for reverse transcription RNA experiments in HepG2 cells. When the oligo(dT) 18 primers were used, the relative expression of hsa-circ-0046600 was significantly lower than that measured with the random hexamer primers, while B3GNTL1 mRNA expression was not, confirming that hsa-0046600 has no poly-A tail (Figure 3C). Furthermore, FISH for hsa-circ-0046600 in HepG2 cells showed that hsa-circ-0046600 has a primarily cytoplasmic distribution (Figure 3D).

Figure 3 Features of hsa-circ-0046600. (A) The genomic locus of the B3GNTL1 gene and hsa-circ-0046600 is shown. The reverse splicing product of exon 4, exon 5, and exon 6 of B3GNTL1 was confirmed by Sanger sequencing. (B) Total RNA was extracted from HepG2 cells treated with actinomycin D for 0, 4, 8, 12, and 24 h, and the levels of hsa-circ-0046600 and B3GNTL1 mRNA were evaluated by qRT-PCR. (C) qRT-PCR analysis of hsa-circ-0046600 and B3GNTL1 mRNA in HepG2 cells was performed with different reverse transcription primers. (D) A FISH assay was performed to determine the subcellular localization of hsa-circ-0046600. Data are reported as the mean ± SD. **P-value < 0.01, ***P-value < 0.001.

hsa-circ-0046600 Acts As A Molecular Sponge Of miR-640

To explore the ceRNA mechanism of hsa-circ-0046600, we used CircInteractome and miRDB to predict miRNAs targeted by circRNAs. miR-640 scored high in both prediction methods (the context+ score percentile was 97, the site type was 8 mer-1a). Subsequently, reduced luciferase activity was observed in cells co-transfected with the hsa-circ-0046600 WT/miR-640 construct; however, the luciferase activity in cells transfected with the hsa-circ-0046600-Mut construct was not affected (Figure 4A). miR-640 expression was upregulated in HepG2 and SK-HEP-1 cells after transfection with si-circRNA-0046600 (Figure 4B). Next, we conducted a meta-analysis of seven sets of microarrays containing miR-640 in the GEO database. A total of 403 HCC samples and 134 control samples were included. The heterogeneity between the datasets was small (P = 0.198, I 2 = 30.2%), and the combined standard mean difference (SMD) calculated by the fixed effects model was −0.30 (95% confidence interval (CI): −0.61- −0.08, P = 0.008). There were no significant publication biases between the datasets (Figure 4C). These results indicate that miR-640 is significantly downregulated in HCC. In addition, qRT-PCR and a correlation analysis showed a negative correlation between the expression of hsa-circ-0046600 and miR-640 in 20 HCC tissue samples (r = −0.5689, P = 0.0089) (Figure 4D). The above results indicate that hsa-circ-0046600 acts by targeting miR-640 in HCC.

Figure 4 Correlation of miR-640 with hsa-circ-0046600 and the expression of miR-640 in HCC. (A) Relative luciferase activity was reduced after the co-transfection of hsa-circ-0046600-WT and miR-640 mimics. The predicted binding sites, sequences and models of hsa-circ-0046600-Mut are shown. (B) siCircRNA0046600 transfection caused a corresponding change in miR-640 in HepG2 cells. (C) The meta-analysis of miR-640 and HCC gene microarrays in the GEO database is shown. (D) miR-640 is negatively correlated with hsa-circ-0046600 in HCC. Data are reported as the mean ± SD. **P-value < 0.01, ***P-value < 0.001.

miR-640 Affects The Migration Of HCC Cells

HepG2 and SK-HEP-1 cells were seeded into 6-well plates, transfected with miR-640 mimics and harvested or scraped after transfection for the indicated times. Nontransfected cells served as a control. Western blotting showed that when the expression level of miR-640 was increased, the expression of N-cadherin was downregulated compared with that in the control group (P < 0.05, Figure 5A). Wound scratch tests showed that after 48 h, the increased expression of miR-640 reduced the migration of SK-HEP-1 cells (the migration distance was 75.10% relative to the NC group, P=0.0003) and HepG2 cells (the migration distance was 58.00% relative to the NC group, P = 0.034, Figure 5B–E).

Figure 5 miR-640 inhibits the migration of HepG2 and SK-HEP-1 cells. (A) After transfection of miR-640 mimics, the expression level of N-cadherin in HepG2 and SK-HEP-1 cells was lower than that in NC group. (B, C) Quantification of relative protein in HepG2 and SK-HEP-1 cells transfected with NC, miR-640 mimics. (D) The migration rates of HepG2 cells were determined by a wound scratch assay. The migratory ability of miR-640 mimics group was lower than that of the control group. (E) Quantification of relative migration in HepG2 cells transfected with NC, miR-640 mimics. (F) The migration rates of SK-HEP-1 cells were determined by a wound scratch assay. The migratory ability of miR-640 mimics group was lower than that of the control group. (G) Quantification of relative migration in SK-HEP-1 cells transfected with NC, miR-640 mimics. *P-value < 0.05, **P-value < 0.01, ***P-value < 0.001.

Prediction Of miR-640 Target Genes

An analysis of the TCGA database using R software revealed 879 differentially expressed mRNAs in HCC. Target gene prediction of miR-640 was performed using three prediction platforms, and a total of 64 mRNAs were predicted to be repeated (Figure 6A). By intersecting these 64 mRNAs with the 879 differentially expressed mRNAs, we obtained a total of 37 miR-640-targeted mRNAs differentially expressed between HCC and control samples, of which 35 mRNAs were significantly upregulated in HCC (Figure 6B).

Figure 6 Functional enrichment analysis and PPI network of miR-640 target genes. (A) The intersection of miR-640 target genes was predicted with bioinformatics software. (B) The intersection of the differentially expressed genes and the miR-640 target genes in the TCGA liver hepatocellular carcinoma (LIHC) dataset is shown, and 35 of the 37 genes were highly expressed in HCC. (C) GO and KEGG analyses of the 37 genes are shown. (D) The PPI network of the 37 genes is shown.

Functional Enrichment Analysis And PPI Network Of The Target Genes

Table 2 shows the results of GO enrichment and KEGG pathway analyses by Metascape for the 37 target genes. Representative pathways revealed by the GO analysis included four items: neural crest cell development, protein kinase activity regulation, response to alcohol, and mRNA metabolic process regulation. The KEGG analysis showed that the Rap1 signalling pathway was the main representative pathway (Table 2, Figure 6C). Among the pathways included, the Rap1 signalling pathway, protein kinase regulation and mRNA metabolic process regulation are the most relevant to tumours. The PPI network indicated that of the 37 mRNAs, the downstream proteins of 9 mRNAs formed 8 PPI pairs (Figure 6D). Upon combining this result with those from existing studies on miR-640, we verified HIF-1α as the downstream protein of hsa-circ-0046600. Consistent with the predicted results, the experimental results showed that in HepG2 and SK-HEP-1 cells, HIF-1α expression was significantly downregulated in the si-circ-004600 group compared to the NC group (P<0.05, Figure 2D).

Table 2 GO Enrichment Terms And KEGG Pathway Enrichment Analysis Of The Target Genes

Discussion

As a large portion of the non-coding RNA family, circRNAs have generated strong interest in the past few years.¹⁵ Second-generation sequencing techniques have demonstrated that circRNAs are abundant non-coding RNAs that are ubiquitous throughout eukaryotes.¹⁶ Similar to lncRNAs and miRNAs, circRNAs function as transcriptional competitors or enhancers of parental genes and translatable transcripts, interacting with RNA-binding proteins (RBPs) and even translating directly into proteins.¹⁷ Because of their ubiquitous existence and diversity, they are particularly important in cancer research. At present, an increasing number of circRNAs related to HCC, such as circMTO1, cSMARCA5 and Cul2, have been described.^18–20 However, there are still many unknown HCC-related circRNAs that have not been identified. We performed a bioinformatic analysis and explored the mechanism of the newly identified circRNA, hsa-circ-0046600, in HCC. First, we compared the microarray data of Cao et al, who reported that hsa-circ-004660 (hsa-circ-102272) was among the top 10 differentially expressed genes. According to our data, hsa-circRNA-0046600 is upregulated in most hepatoma cell lines, particularly in the HepG2 and SK-HEP-1 cell lines. Therefore, to obtain clearer results, we chose to verify the biological function and ceRNA mechanism of hsa-circ-0046600 in HepG2 and SK-HEP-1 cells. In addition, hsa-circRNA-0046600 expression was significantly upregulated in 90% of our tested HCC tissue samples compared with that in the paired adjacent normal tissue samples, as measured by qRT-PCR, and was positively correlated with tumour size, stage and pathological vascular invasion. Next, to verify whether there was a functional difference associated with hsa-circ-0046600, an hsa-circ-0046600 siRNA was constructed, and its specificity was verified. In vitro experiments showed that the migration of HCC cells after siRNA transfection and the expression of the migration-related protein N-cadherin were reduced. The effect of hsa-circ-0046600 on migration was not as obvious in HepG2 cells as it was in SK-HEP-1 cells. After repeated verification, we believe that the difference may be due to the poor migration of HepG2 cells themselves, resulting in a diminishing migration distance after 48 h. In addition, we demonstrated the proliferation and apoptosis of these two cell types after siRNA regulation. There was no significant difference between the experimental and control groups. However, hsa-circ-0046600 is closely related to tumour size, which contradicts the experimental results. We believe that this discrepancy may exist because tumour proliferation is a complex process in which hsa-circ-0046600 has little effect despite its involvement.

Unlike linear RNAs, circRNAs are produced by a covalent linkage between the 5ʹ and 3ʹ ends of an RNA molecule.²¹ circRNAs can consist of a single exon, multiple exons or multiple exons and introns, and they lack 5ʹ and 3ʹ ends and the circular structure of the poly-A tail.¹⁹ Therefore, circRNAs can avoid degradation by exonucleases and RNase R, making these molecules conservative and stable. According to the sequence of hsa-circ-0046600, we found that hsa-circ-0046600 consists of three exons from its parent gene, B3GNTL1. Due to the nature of circRNAs, hsa-circ-0046600 degraded much more slowly than did B3GNTL1 mRNA following treatment with actinomycin D, and hsa-circ-0046600 could not undergo reverse transcription with oligo(dT) 18 primers. We conclusively demonstrated that the molecules detected in the experiment were circRNAs instead of mRNAs. In addition, we also determined by FISH that hsa-circ-0046600 is present mainly in the cytoplasm. Therefore, we speculate that hsa-circ-0046600 can act as a sponge.

The ceRNA hypothesis posits that mRNAs, lncRNAs and circRNAs can release miRNA-targeted mRNA transcripts by competing with base-complementary miRNA response elements (using a sponge mechanism).²² Using bioinformatic software, we predicted that miR-640 is the miRNA most likely to bind to hsa-circ-0046600. Furthermore, the combination of hsa-circ-0046600-miR-640 was confirmed with the transfection of siCirc0046600 and a dual-luciferase reporter assay. However, the difference in the expression of miR-640 was not confirmed in HCC; therefore, we performed a meta-analysis of relevant microarrays to confirm that miRNA-640 is expressed at a higher level in normal tissue than in HCC tissue. In summary, we determined that the hsa-circ-0046600-miR-640 pathway plays an important role in the migration of HCC cells.

Next, we predicted the possible downstream mRNA binding partners of miR-640 and compared our results with those from the TCGA HCC dataset. This approach confirmed that the expression of 35 of the 37 closely related HCC mRNAs was upregulated in HCC tissue compared with normal tissue, and this result was in accordance with the expected results. Subsequently, we performed a functional analysis of the 37 mRNAs and constructed a PPI network and confirmed that there are 9 core mRNAs in the hsa-circ-0046600-miR-640 pathway. Yi-Chun Zhu et al verified that miR-640 bound to the 3ʹ-untranslated region (UTR) of HIF-1α mRNA and then inhibited the expression of HIF-1α with dual-luciferase reporter and Western blot experiments. Based on the functional enrichment analysis and the PPI network, as well as the results described by Yi-Chun Zhu et al, we selected HIF-1α as the downstream protein of hsa-circ-0046600 for verification by Western blot.²³ We used Western blotting to confirm that HIF-1α expression was downregulated after transfection with siCirc-0046600. The above results demonstrate that the hsa-circ-0046600-miR-640 combination can affect the migration of HCC cells through a variety of mechanisms.

Conclusion

Our study indicates that hsa-circ-0046600 is highly expressed in HCC tissues and cells and is closely related to tumour size, TNM stage, and pathological vascular invasion in HCC patients. hsa-circ-0046600 promotes the expression of HIF-1α by acting as a sponge of miR-640, further affecting the migration of HCC cells. hsa-circ-0046600 may be a potential new target for the diagnosis and treatment of HCC.