CircGPR137B/miR-4739/FTO feedback loop suppresses tumorigenesis and metastasis of hepatocellular carcinoma

Background Emerging evidence indicates that circular RNAs (circRNAs) and m6A RNA methylation participate in the pathogenesis and metastasis of multiple malignancies including hepatocellular carcinoma (HCC). However, it remains undocumented how circRNAs form a feedback loop with the m6A modification contributing to HCC. Methods A novel hsa_circ_0017114 (circGPR137B) was identified from three pairs of primary HCC and adjacent normal tissues by circRNA expression profiling. The association of circGPR137B and miR-4739 with clinicopathological parameters and prognosis in patients with HCC was analyzed by RT-qPCR, fluorescence in situ hybridization and TCGA cohorts. The role of circGPR137B in HCC was estimated in vitro and in vivo. RT-qPCR, western blot, m6A dot blot, RIP, MeRIP and dual-luciferase reporter assays were used to validate the reciprocal regulation of the feedback loop among circGPR137B, miR-4739 and m6A demethylase FTO. Meanwhile, the expression, function and prognosis of FTO in HCC were investigated by RT-qPCR, western blot, TCGA and rescue experiments. Results We identified a new dramatically downregulated circGPR137B in HCC tissues, and found that downregulation of circGPR137B or upregulation of miR-4739 was associated with poor prognosis in patients with HCC. Ectopic expression of circGPR137B strikingly repressed the proliferation, colony formation and invasion, whereas knockdown of circGPR137B harbored the opposite effects. Moreover, restored expression of circGPR137B inhibited tumor growth and lung metastasis in vivo. Further investigations showed that circGPR137B, co-localized with miR-4739 in the cytoplasm, acted as a sponge for miR-4739 to upregulate its target FTO, which mediated m6A demethylation of circGPR137B and promoted its expression. Thus, a feedback loop comprising circGPR137B/miR-4739/FTO axis was formed. FTO suppressed cell growth and indicated favorable survival in patients with HCC. Conclusion Our results demonstrate that circGPR137B inhibits HCC tumorigenesis and metastasis through the circGPR137B/miR-4739/FTO feedback loop. This positive feedback mechanism executed by functional coupling between a circRNA sponge and an m6A modification event suggests a model for epigenetics. Supplementary Information The online version contains supplementary material available at 10.1186/s12943-022-01619-4.


Introduction
Hepatocellular carcinoma (HCC) is one of the most frequent pernicious tumors with higher incidence and the second cancer-related mortality in China [1] and the fifth mortality in the United States [2]. Despite the application of various novel treatments, HCC cases still have poor prognoses due to the malignant infiltration and metastasis [3]. The dysregulation of noncoding RNAs (ncRNAs) has been implicated in the pathogenesis of HCC [4,5]. Therefore, the identification of candidate ncRNAs may offer new diagnostic and therapeutic targets for HCC.
Furthermore, m 6 A modification of circNSUN2 promotes colorectal liver metastasis [35] and m 6 A-modified circ-SORE sustains sorafenib resistance in HCC [36]. METTL14-mediated m 6 A modification of circORC5 suppresses gastric cancer progression [37]. However, it remains unclear how a circRNA forms a feedback loop with m 6 A modification contributing to HCC. We herein identified a novel circ_0017114 (circGPR137B), and found that low expression of circG-PR137B or high expression of miR-4739 was associated with poor survival and tumor recurrence in patients with HCC. CircGPR137B curbed HCC tumorigenesis and metastasis through the circGPR137B/miR-4739/ FTO feedback loop. Our findings imply that circG-PR137B is a potential predictive biomarker and therapeutic target for HCC.

Clinical data
The HCC clinical data, as well as the expression of miR-NAs, FTO and YTHDC2 were collected from The Cancer Genome Atlas (TCGA) dataset (https:// genome-cancer. ucsc. edu). A tissue microarray (TMA) including 87 paired HCC tissue samples (Lot No. XT16-029) was purchased from Shanghai Outdo Biotech Company. Ten paired frozen HCC tissue samples were stored in liquid nitrogen in our laboratory. Protocols were approved by the Ethics Committee of Zhong Shan Hospital. The specimens were classified based on the TNM staging system, and diagnosed by two independent pathologists.

CircRNA microarray analysis
Total RNA from 3 paired HCC tissues was quantified using the NanoDrop ND-1000. The sample preparation and microarray hybridization were conducted based on Arraystar's standard protocols. CircRNA microarray analysis was performed as previously reported [10].

Fluorescence in situ hybridization (FISH)
The expression and cellular localization of circGPR137B or miR-4739 in HCC tissue samples were detected by FISH analysis. Digoxin-labeled probe sequence for circG-PR137B (5′-GGC TTT GAA AAT CAC TCT GTG AAC ATA GCA -3′) and Biotin-labeled probe sequence for miR-4739 (5′-AGG GCC CCT CCG CTC CTC CTC CCT T-3′) were synthesized for FISH analysis. The specific description of FISH analysis was performed as previously reported [10]. The analysis software Image-pro plus 6.0 (Media Cybernetics, Rockville, MD, USA) was used to analyze immunofluorescence accumulation optical density (IOD).

Quantitative real-time PCR (RT-qPCR)
Total RNA was extracted using TRIzol, reverse transcription was performed using M-MLV and cDNA amplification using the SYBR Green Master Mix kit (Takara, Otsu, Japan). Total RNA was isolated using a High Pure miRNA isolation kit (Roche) and RT-PCR using a TaqMan Micro-RNA Reverse Transcription kit (Life Technologies). The nuclear and cytoplasmic fractions were isolated using NE-PER Nuclear and Cytoplasmic Extraction Reagents (Thermo Scientific). The primers used were listed in Supplementary Table S1.

Actinomycin D and RNase R treatment
Transcription was prevented by the addition of 2 mg/ ml Actinomycin D or DMSO (Sigma-Aldrich, St. Louis, MO, USA) as the negative control. Total RNA (2 μg) was incubated for 30 min at 37 °C with 3 U/μg of RNase R (Epicentre Technologies, Madison, WI, USA).

RNA immunoprecipitation (RIP) and m 6 A RIP (MeRIP)
RIP assay was performed in HepG2 and Hep3B cell lines using a Magna RIP RNA-binding protein Immunoprecipitation Kit (Millipore) according to the manufacturer's instructions. Antibodies for RIP assays against Ago2 and IgG were purchased from Abcam (ab5072, Rabbit polyclonal antibody, Cambridge, MA, USA). Anti-m 6 A antibody (A-1801-020, Epigentek) and Magna MeRIP m6A Kit (17-10,499, Millipore).were used for the MeRIP assay.

In vivo tumorigenesis assay
BALB/c (nu/nu) nude mice (male, 6-8 week) were purchased from Shanghai SIPPR-BK Laboratory Animal Co. Ltd. (Shanghai, China). All the animals were handled according to institutional guidelines, and approved by the Animal Ethics Committee of Zhong Shan Hospital. HepG2 cells stably transfected with pLV-circGPR137B and empty vectors (pLV-circ-Luci (2A) puro，GenePharma，Shanghai, China) were labeled by luciferase. The mice were subcutaneously inoculated with 6 × 10 7 HepG2 cells stably transfected with pLV-circGPR137B and empty vectors or Sk-hep-1 cells stably transfected with pLV-sh-circGPR137B/sh-NC. The tumor weight and size were measured every other day, and the tumor volume was calculated based on the formula: length × width 2 /2.

Orthotopic implantation liver tumor model
To establish an orthotopic implantation liver tumor model, 6 × 10 7 HepG2 cells stably transfected with pLV-circGPR137B and empty vectors which were labelled by luciferase were injected orthotopically into the left liver lobe of nude mice. Each mouse was injected with 100ul (6× 10 6 cells). Each group had 5 mice. After 4 weeks, the luciferase signaling was collected for pLV-circG-PR137B and control groups. When mice were sacrificed, the liver weight of the mice was encoded.

Caudal vein pulmonary metastasis model
HepG2 cells stably transfected with circGPR137B or pcDNA3.1 were cultured in complete medium. When the cells were 70% confluent, the medium was replaced with fresh medium to remove dead and detached cells. Subsequently, 6 × 10 7 cells were used to inject into mice via the tail vein. Each mouse was injected with 100ul (6× 10 6 cells). Each group had 5 mice. The progression of pulmonary metastasis was investigated for 4 weeks.

Hematoxylin and eosin (HE) staining
Mice tumor tissues were harvested, fixed in 4% paraformaldehyde, and preserved at normal atmospheric temperature. The liver tissues were sliced into 5 μm sections and stained with HE for the histological studies.

Immunohistochemistry (IHC) analysis
IHC analysis of Ki-67 levels was conducted on paraffin slides of mouse liver and metastatic lung tumor tissues using anti-Ki-67 antibodies (Abcam, Cambridge, UK) and anti-FTO (#ab126605, Abcam). The detailed description of IHC was conducted as previously reported [5].

Statistical analysis
Statistical analysis was executed with GraphPad Prism 7 (La Jolla, CA, USA). In brief, the values are expressed as the mean ± standard deviation (SD). Student's test and analysis of variance were used for comparisons between groups. Kaplan-Meier analysis was used to evaluate the association of circGPR137B, miR-4739 or FTO with HCC prognosis. Pearson Correlation Analysis was used to analyze the correlation of circGPR137B with miR-4739. The categorical data were analyzed by chi-square of Fisher's exact tests. A P-value of < 0.05 was considered statistically significant.

Downregulation of circGPR137B is associated with poor prognosis in patients with HCC
To screen the differentially expressed circRNAs in HCC, we conducted a circRNA expression profiling in 3 paired HCC tissue samples. The heatmap indicated 6524 upregulated circRNAs and 5986 downregulated circRNAs in HCC tissues as compared with the adjacent normal tissues (Fig. 1A). According to the criteria of P < 0.05 and FC > 1.5, 96 upregulated circRNAs and 51 downregulated circRNAs were further identified (Fig. 1B). The volcano plot showed that hsa_circ_0017114 derived from a linear RNA GPR137B, displayed a dramatical downregulation in HCC (P = 0.0098; FC = 15.5613) and was nominated as circGPR137B (Fig. 1C). The circRNA profiling (Fig. 1D) and RT-qPCR (Fig. 1E) indicated a decreased circG-PR137B expression in HCC as compared with the adjacent normal tissues. FISH analysis further validated this result in 87 paired HCC tissues (P = 0.0113; Fig. 1F, G).
According to circGPR137B expression and the prognostic data, we obtained a cutoff value of circGPR137B, and divided the cases into circGPR137B-high and circGPR137B-low groups. We found that circGPR137B harbored no relationship with the clinicopathological characteristics in HCC (Supplementary Table S2). The cases as well as the early stage ones but not the stage III cases ( Supplementary Fig. S1) with circGPR137B-low expression possessed a poorer survival as compared with those with circGPR137B-high expression (Fig. 1H, I). A Cox proportional hazard model was established, and univariate and multivariate analysis unveiled that TNM stage rather than circGPR137B expression was an independent prognostic factor in HCC (Supplementary Table S3).

The characteristics of circGPR137B in HCC cells
hsa_circ__0017114 (chr1:236332005-236,372,209), derived from exon 1, 7 regions within G protein-coupled receptor 137B (GPR137B) locus, is located on chromosome 14q1(q42.3), and nominated as circGPR137B, whose genomic sequence is 40,204 nt and spliced length is 1537 nt. (Fig. 2A). The mature sequence for circGPR137B is shown in supplementary data. When exposed to RNase R treatment for 2 h, GPR137B expression level was extremely lowered, but circGPR137B exhibited no significant alternation duo to its resistance to RNase R in HepG2 and Hep3B cells (Fig. 2B). We then investigated the stability of circG-PR137B after exposure to the actinomycin D, an inhibitor of the transcription at indicated time points. We found that the transcription half-life of circGPR137B was much longer than that of linear GPR137B in HepG2 and Hep3B cells (Fig. 2C). The cytoplasmic or nuclear localization of circGPR137B and GPR137B was determined by qPCR and FISH analyses, which indicated that, circGPR137B was predominantly localized in the cytoplasm of HCC tissues in comparison with the linear GPR137B (Fig. 2D, E) .

circGPR137B inhibits in vitro cell growth and invasion
We estimated the expression of circGPR137B in multiple HCC cell lines by qRT-PCR analysis, which indicated that circGPR137B harbored a higher expression in SK-hep-1 cell line, but a lower expression in HepG2 and Hep3B cell lines in comparison with the normal liver tissues (Fig. 3A). Then, the overexpression efficiency of circGPR137B in HepG2 and Hep3B cell lines as well as the silencing efficiency of sh-circGPR137B in SK-hep-1 cell line was defined by qRT-PCR analysis (Fig. 3B). Further investigations indicated that, ectopic circGPR137B expression repressed the cell viability (Fig. 3C), colony formation (Fig. 3D) and cell invasive capabilities (Fig. 3E) in HepG2 and Hep3B cells, whereas knockdown of circG-PR137B harbored the opposite effects ( Fig. 3C-E).

MiR-4739 has a negative correlation with circGPR137B expression and poor prognosis in HCC
Depending on the circRNA expression profiling and miRbase, circGPR137B was identified to possess the potential to bind with 5 miRNAs (miR-1249-5p, miR-3916, miR-6760-5p, miR-214-3p and miR-4739) (Fig. 4A) and their binding sites can be shown in Fig. 4B. We analyzed the luciferase activity of circG-PR137B in hepG2 cells after treatment with 5 miRNA mimics and found that the luciferase activities of circGPR137B were markedly reduced by miR-4739 rather than the other miRNAs (Fig. 4C). FISH analysis indicated that miR-4739 expression levels were significantly increased (Fig. 4D, E) and had a negative correlation with circGPR137B ( Fig. 4F), but harbored no correlation with GPR137B in HCC (Supplementary Fig. S2). The consistent result was validated in 220 HCC tissue samples in the TCGA cohort (Fig. 4G). Then, a cutoff value of miR-4739 (0.586) was obtained and categorized the cases into miR-4739-high and miR-4739-low groups. Increased miR-4739 expression was linked to the pathological stage (P = 0.029) and tumor size (P = 0.018) in HCC (Supplementary Table  S4). Kaplan-Meier analysis revealed that, patients with high miR-4739 expression harbored poorer survival as compared with those with low miR-4739 expression (Fig. 4H). The advanced cases as well as early stage ones with high miR-4739 expression harbored no differences in overall survival and tumor recurrence as compared with those with low miR-4739 expression ( Supplementary Fig. S3). Univariate and multivariate analyses demonstrated that high expression of miR-4739 was an independent prognostic factor of tumor recurrence rather than poor survival in patients with HCC (Supplementary Tables S5, 6).

CircGPR137B acts as a sponge for miR-4739 in HCC
The binding sites between WT circGPR137B 3'UTR and miR-4739 can be demonstrated in Fig. 5A. The expression levels of miR-4739 were reduced by ectopic expression of circGPR137B in HepG2 and Hep3B cells (Fig. 5B), but miR-4739 mimics had no effects on circGPR137B expression ( Supplementary Fig.S4). PRL-TK-pMIR-Luc vector including WT or Mut circGPR137B 3'UTR was co-transfected with miR-4739 mimics into HepG2 and Hep3B cells, and the result indicated that, miR-4739 mimics lowered the luciferase activity of WT circGPR137B 3'UTR, but exerted no effect on the luciferase activity of Mut circGPR137B 3'UTR and Mut miR-4739 also had no impact on that of WT circGPR137B 3'UTR in HepG2 and Hep3B cells (Fig. 5C). Moreover, the Ago2 occupancy  Table  S7). Further RIP assay implementation indicated that the enrichment of endogenous circGPR137B and miR-4739 pulled-down from Ago2-expressed HepG2 and Hep3B was markedly elevated in Ago2 pellet as compared with the input control (Fig. 5D). In addition, a cytoplasmic co-localization between circGPR137B and miR-4739 was shown in HepG2 cells by FISH analysis (Fig. 5E). The cotransfection of circGPR137B lentiviruses with miR-4739 mimics in HepG2 and Hep3B cells indicated that miR-4739 promoted the cell proliferation and invasion, and counteracted tumor suppressive effects of circGPR137B (Fig. 5F, G).

FTO is regulated by circGPR137B/miR-4739 axis and indicates a favorable prognosis in HCC
According to the analyses of TargetScan7.1, the targets of miR-4739 were screened and m 6 A demethylase FTO might be a target of miR-4739 (Fig. 6A). It was confirmed that miR-4739 mimics could reduce the luciferase activity of WT FTO 3'UTR, but exerted no effect on Mut FTO 3'UTR and Mut miR-4739 had no impact on the luciferase activity of WT FTO 3'UTR as compared with miR-NC group in HepG2 and Hep3B cells (Fig. 6B). Furthermore, qRT-PCR (Fig. 6C) and western blot analysis (Fig. 6D) indicated that miR-4739 mimics markedly decreased expression of FTO, p21, p27, cleaved caspase-3/9 and E-cadherin and increased expression of N-cadherin and vimentin. miR-4739 also reversed circGPR137B-induced upregulation of FTO, p21, p27, cleaved caspase-3/9 and E-cadherin and downregulation of N-cadherin and vimentin in HepG2 and Hep3B cells. Then, after the co-transfection of FTO plasmids and miR-4739 mimics into HepG2 and Hep3B cells for 48 h, we found that ectopic expression of FTO prevented the cell proliferation, colony formation, and cell invasion, and reversed tumor-promoting effects caused by miR-4739 ( Fig. 6E-G).
According to the RT-qPCR and western blot analyses, the expression of FTO was dramatically decreased in 10 pared HCC tissue samples as compared with the adjacent normal tissues (Fig. 7A, B). This result was validated in 50 paired and 371 unpaired HCC tissues (Fig. 7C). Then, we found that decreased expression of FTO was associated with age and distant metastasis in HCC (Supplementary  Table S8). Kaplan-Meier analysis showed that the cases with high FTO expression indicated favorable survival as compared with those with low FTO expression (Fig. 7D). Univariate and multivariate analyses uncovered that FTO low expression was an independent prognostic factor of poor survival in HCC (Supplementary Table S9).

FTO mediates m 6 A modification of circGPR137B to form a positive feedback loop with circGPR137B
It has been shown that m 6 A modification of circRNA is implicated in cancer [35][36][37]. We supposed that FTOmediated m 6 A modification of circGPR137B could repress HCC progression. MeRIP and m 6 A dot blot indicated that FTO overexpression reduced the total m 6 A levels (Fig. 8A, B) and circGPR137B m 6 A levels (Fig. 8C), but increased circGPR137B RNA levels in HepG2 and Hep3B cells (Fig. 8D). Knockdown of FTO could decrease circGPR137B expression (Fig. 8E). Then, knockdown of circGPR127B increased miR-4739 expression and downregulated FTO (Fig. 8F), and the miR-4739 inhibitor increased FTO expression in the SK-hep-1 cell line (Fig. 8G). RIP assays validated that FTO could bind to circGPR137B in HCC cells (Fig. 8H). However, FTO had no effects on GPR137B m 6 A levels and could not bind to GPR137B (Supplementary Fig. S5). Further functional assays demonstrated that knockdown of circG-PR137B promoted cell proliferation and invasion and counteracted FTO-induced antitumor effects in HepG2 and Hep3B cells (Fig. 8I, J).

CircGPR137B represses in vivo HCC tumorigenesis
Luciferase-labeled HepG2 cells stably transfected with circGPR137B or pcDNA3.1 were injected into the nude mice. Luciferase signals were examined using ex vivo imaging, and the result indicated that the luciferase signals of the tumor tissues in circGPR137B group were slighter than those in the control group (Fig. 9A). The tumors were removed from the mice when the mice was sacrificed, and those in the circGPR137B group seemed smaller than those in the control group. Statistical analysis indicated a decreased tumor size and weight in the circGPR137B group as compared with the control group (Fig. 9B, C). HE staining indicated a reduced tumor cell number and IHC analysis showed lowered Ki-67 and FTO levels in the circGPR137B group as compared with the control group (Fig. 9D). qRT-PCR analysis further showed that overexpression of circGPR137B decreased miR-4739 expression and increased FTO expression in the circGPR137B group as compared with the control group (Fig. 9E). Knockdown of circGPR137B in vivo exhibited a promoting effect on tumor growth (Fig. 9F). IHC analysis indicated elevated Ki-67 levels and lowered FTO levels in the sh-circGPR137B group as compared with the control group (Fig. 9G). qRT-PCR analysis further showed that knockdown of circGPR137B increased miR-4739 expression and decreased FTO expression in the sh-circGPR137B group as compared with the control group (Fig. 9H).

CircGPR137B suppresses in vivo lung metastasis of HCC
Luciferase-labeled HepG2 cells stably transfected with circGPR137B or CON were injected into the mice to construct orthotopic liver tumor models. Luciferase signals were examined using ex vivo imaging, and the result demonstrated liver tumor growth with peritoneal metastasis. The luciferase signals were obviously fainter in the circGPR137B group than those in the control group (Fig. 10A). The size of the liver tumors in the circG-PR137B group appeared smaller than those in the control group, and a decreased liver weight was shown in circGPR137B group as compared with the control group (Fig. 10B), but the body weight in mice had no difference between circGPR137B and control groups (Supplementary Fig. S6). Besides, HepG2 cells stably transfected with circGPR137B or CON were injected into the caudal vein to establish pulmonary metastasis models. we found that the formation rate of the metastatic lung tumors was 60.0% (3/5) in the circGPR137B group and 100.0% (5/5) in the control group, and a decreased lung weight was shown in the circGPR137B group as compared with the control group (Fig. 10C). HE staining indicated a diminished cell number in orthotopic liver tumors and metastatic lung tumor cells and IHC analysis exhibited a dwindling Ki-67 index and elevated FTO expression in liver tumors and metastatic lung tumor tissues in the circGPR137B group as compared with the control group (Fig. 10D). qRT-PCR indicated that overexpression of circGPR137B decreased miR-4739 expression and increased FTO expression in orthotopic liver tumors as compared with the control group (Fig. 10E).

Discussion
It has been shown that upregulated circPPP1R12A [18] and downregulated circYAP1 [38] are related to the tumor metastasis and TNM stage, and can represent independent prognostic factors in patients with HCC. Herein, we identified a novel differentially-expressed circGPR137B by a circRNA profiling and validated its downregulation in HCC tissues. The downregulation of circGPR137B was associated with poor HCC survival Fig. 10 CircGPR137B inhibits in vivo lung metastasis of HCC. A Ex vivo imaging of the luciferase signals in orthotopic liver tumors with peritoneal infiltration by stable transfection with circGPR137B or control group into luciferase-labeled HepG2 cells. B Schematic representation of the comparison of orthotopic liver tumors, and statistical analysis of the total liver weight between circGPR137B and control groups. C Schematic representation of the comparison of metastatic pulmonary tumors, and statistical analysis of the total lung weight between circGPR137B and control groups. D HE staining of the tumor tissue cells from orthotopic liver tumors and metastatic pulmonary tumors in circGPR137B and control groups and IHC analysis of the Ki-67 proliferation index and FTO expression levels in orthotopic liver tumors and metastatic pulmonary tumors between circGPR137B and control groups. E. RT-qPCR analysis of the expression levels of circGPR137B, miR-4739 and FTO in orthotopic liver tumors between circGPR137B and control groups as well as early-stage cases, and might offer a potential marker for early detection of HCC.
CircRNAs play a dual role in cancer. circPPP1R12A, circRNA-100,338 or circSLC3A2, for example, acts as an oncogenic factor [18,23,24], but circTRIM33-12, circ_0051443 or circYAP1 acts as a tumor suppressor in HCC [22,25,38]. Herein, we investigated the role of circGPR137B in HCC cells and found that circG-PR137B repressed cell proliferation, colony formation, and cell invasion as well as liver tumor infiltration and pulmonary metastasis, while silencing circGPR137B harbored the opposite effects. These findings suggested that circGPR137B might be a tumor suppressor in HCC.
It is known that circRNAs can act as miRNA sponges to regulate tumor progression [26][27][28]. Here, we identified circGPR137B-specific binding with miR-4739 in HCC cells. MiR-4739 regulates osteogenic and adipocytic differentiation of bone marrow stromal cells [39] and pleural fibrosis [40], and act as a potential biomarker for diabetes [41]. We here found that increased miR-4739 expression was associated with pathological stage and tumor size, and acted as an independent prognostic factor of tumor recurrence in HCC. MiR-4739 can also be sponged by lncRNA VPS9D1-AS1 to prompt prostate cancer tumorigenesis [42]. We validated that circG-PR137B could act as a miR-4739 sponge, and miR-4739 reversed the tumor-suppressive effects of circGPR137B, Fig. 11 Schematic illustration of circGPR137B suppressing proliferation and metastasis of liver cancer through a positive feedback loop. CircGPR137B induces FTO protein expression by sponging miR-4739 which binds to FTO mRNA in order to inhibit its expression. FTO protein enters into cellular nuclear to mediate m6A demethylation of circGPR137B and promotes its expression. Thus, a positive feedback forms on circGPR137B production, resulting in strongly repressing liver cancer cell proliferation, colony formation, invasion, and lung migration indicating that, circGPR137B might act as a sponge for miR-4739 to inhibit HCC progression.
Increasing data show that FTO as an m 6 A demethylase plays an oncogenic role in acute myeloid leukemia [43], breast cancer [44] and melanoma [45], but restrains ovarian cancer stem cell self-renewal [46]. We here found that FTO was identified as a direct target of miR-4739 and indicated a favorable prognosis in HCC. hsa_circ_0072309 promotes tumorigenesis and invasion in non-small cell lung carcinoma by regulating miR-607/FTO axis [47]. However, our findings indicated that circGPR137B acted as a tumor suppressor in HCC by sponging miR-4739/ FTO axis. Modification of circRNAs m 6 A can regulate cancer progression [35][36][37]. We further confirmed that FTO modulated m 6 A-dependent modification of circG-PR137B and demonstrated that circGPR137B inhibited the tumorigenesis and metastasis of HCC through the circGPR137B/miR-4739/FTO feedback loop (Fig. 11).

Conclusion
Altogether, our findings demonstrate that, downregulation of circGPR137B or upregulation of miR-4739 is associated with poor prognosis in HCC. Restored circG-PR137B suppresses the malignancy of HCC through the circGPR137B/miR-4739/FTO feedback loop. Our study offers direct evidence for a feedback loop formed by cir-cRNA and m6A modification, producing new insights into epigenetics.