N6-methyladenosine-modified CircRNA-SORE sustains sorafenib resistance in hepatocellular carcinoma by regulating β-catenin signaling

Accumulating evidence suggests that the primary and acquired resistance of hepatocellular carcinoma (HCC) to sorafenib is mediated by multiple molecular, cellular, and microenvironmental mechanisms. Understanding these mechanisms will enhance the likelihood of effective sorafenib therapy. In vitro and in vivo experiments were performed and clinical samples and online databases were acquired for clinical investigation. In this study, we found that a circular RNA, circRNA-SORE, which is up-regulated in sorafenib-resistant HCC cells, was necessary for the maintenance of sorafenib resistance, and that silencing circRNA-SORE substantially increased the efficacy of sorafenib-induced apoptosis. Mechanistic studies determined that circRNA-SORE sequestered miR-103a-2-5p and miR-660-3p by acting as a microRNA sponge, thereby competitively activating the Wnt/β-catenin pathway and inducing sorafenib resistance. The increased level of circRNA-SORE in sorafenib-resistant cells resulted from increased RNA stability. This was caused by an increased level of N6-methyladenosine (m6A) at a specific adenosine in circRNA-SORE. In vivo delivery of circRNA-SORE interfering RNA by local short hairpin RNA lentivirus injection substantially enhanced sorafenib efficacy in animal models. This work indicates a novel mechanism for maintaining sorafenib resistance and is a proof-of-concept study for targeting circRNA-SORE in sorafenib-treated HCC patients as a novel pharmaceutical intervention for advanced HCC.


Introduction
Hepatocellular carcinoma (HCC) is the most common cancer and the leading cause of cancer-associated mortality worldwide. Surgical resection is the first choice treatment for HCC, but not for advanced HCC patients [1]. For unresectable HCC, treatment options include transplantation, ablation, transarterial chemoembolization, targeted therapies and immunotherapies. Sorafenib, a multikinase inhibitor, is a first-line targeted drug approved by the US Food and Drug Administration for advanced HCC [2]. According to the results of the Sorafenib Hepatocellular Carcinoma Assessment Randomized Protocol (SHARP) trial, sorafenib can prolong survival of HCC patients. However, drug resistance limits its efficacy. Although some mechanisms have been reported for sorafenib resistance, such as epithelialmesenchymal transition, the proliferation of cancer stem cells, and metabolic reprogramming, the exact cause is still elusive [3]. Therefore, understanding the underlying molecular basis of HCC sorafenib resistance and developing mechanism-based therapies are urgently needed.
Circular RNAs (circRNAs) have a circular configuration through a typical 5′ to 3′-phosphodiester bond and are recognized as a class of functional non-coding RNAs (ncRNAs). CircRNAs regulate biological processes by mediating alternative splicing (AS) of RNAs, cisregulation of transcription, and by acting as competing endogenous RNAs (ceRNAs) [4]. Importantly, many cir-cRNAs are involved in cell proliferation, differentiation, apoptosis and invasion during tumor progression [5,6]. Recent studies show that the circRNA, CDR1as, can accelerate the proliferation and migration of HCC cells by promoting the expression of AFP via sponging miR-1270 [7], while circRHOT1 can inhibit HCC development and progression via recruiting TIP60 to initiate NR2F6 expression [8]. Circ_0003418 increases the sensitivity of HCC cells to cisplatin by inhibiting the Wnt/β-catenin pathway [9]. However, the roles of circRNAs in HCC sorafenib resistance remain unknown.
N6-methyladenosine (m 6 A) is the most prevalent internal modification associated with eukaryotic mRNAs and ncRNAs and it influences many steps of mRNA metabolism, including splicing, export, translation, and stability [10]. The m 6 A modification has been implicated in various cellular and physiological events, including carcinogenesis [11]. In acute lymphoblastic leukemia (ALL), levels of METTL3 (a component of the methyltransferase complex that catalyzes adenosine methylation) increase in AML patients and it plays an oncogenic role by inhibiting cell differentiation and apoptosis, and promoting cell proliferation through increased c-MYC, BCL2, and PTEN translation [12]. In breast cancer, mammalian hepatitis B X-interacting protein (HBXIP) enhances the expression of METTL3 by suppressing the tumor suppressor, let-7 g, and forms a positive feedback loop to enhance the level of HBXIP by facilitating m 6 A modification of mRNA, which promotes cell proliferation [13]. In HCC, m 6 A modification also regulates the progress of oncogenesis. SIRT1 promotes HCC by increasing the overall m 6 A modification level to attenuate the expression of the tumor suppressor, GNAO1 [14]. However, the functions of m 6 A modification in HCC sorafenib resistance remain elusive. Moreover, the roles of m 6 Amodified circRNAs in HCC sorafenib resistance also need further investigation.
In the present study, we found that circRNA-SORE (also named circRNA_104,797 and circ_0087293) was up-regulated in sorafenib-resistant HCC cells, and was necessary for the maintenance of sorafenib resistance. By acting as a ceRNA and sequestering miR-103a-2-5p and miR-660-3p, circRNA-SORE competitively activates the Wnt/β-catenin pathway and promotes sorafenib resistance. The increased levels of circRNA-SORE in sorafenib-resistant HCC are due to its increased stability resulting from increased N 6 -methyladenosine (m 6 A) levels of a specific adenosine in circRNA-SORE. These results provide a novel mechanism for maintaining sorafenib resistance and demonstrate a proof-of-concept for targeting circRNA-SORE in sorafenib-treated HCC patients as a novel pharmaceutical intervention for advanced HCC.

RNA pull-down
For biotinylated microRNA (miRNA) pull-down assays, biotin-miRNA-transfected HepG2-SR cell lysates were incubated with MyOne™ Dynabead® Streptavidin C1 (Invitrogen, Grand Island, NY, USA) at 4°C for 1 h. For biotinylated probe pull-down assays, washed MyOne™ Dynabead® Streptavidin C1 (Invitrogen, USA) was incubated with cell lysate at 4°C for 1 h for preclearance. The 3′ biotin-labeled circRNA-SORE probe (5′-GAGT TGTTGCTGCTTGATGGAGTC-3′-biotin) was synthesized by Tsingke Biotech and incubated with the beads at room temperature for 10 min for immobilization. Then, the biotinylated beads were incubated with HepG2-SR cell lysate at 4°C overnight. The biotinylated beads were magnetically separated and washed five times. For qRT-PCR assays, Trizol reagent was used to extract total RNA from the beads.
Silencing and overexpression of circRNA-SORE circRNA-SORE siRNA targeting the junction region of the circRNA-SORE sequence was synthesized by Ribo-Bio (Guangzhou, China). Transfection of plasmids, siRNA, or miRNA mimics or inhibitors (RiboBio, Guangzhou, China) was performed using Lipofectamine® 3000 (Invitrogen, USA) according to the manufacturer's protocol. The sequences of siRNAs against specific targets are listed in Table S1.
pLCDH-circRNA-SORE-M has a mutated circRNA-SORE m 6 A site. The sequences of mutant constructs are listed in Table S2. GV248-shcircRNA-104797, pLCDH-ciR, pLCDH-ciR-SORE, the psAX2 packaging plasmid, and the pMD2G envelope plasmid, were transfected into 293 T cells using the standard calcium phosphate transfection method to obtain lentiviral supernatant. Fortyeight hours after transfection, viruses were collected and concentrated by density gradient centrifugation, then frozen at − 80°C for later use.

Tissue samples
A total of 92 randomly selected HCC clinical samples were collected from the Department of General Surgery, Sir Run-Run Shaw Hospital, Zhejiang University, starting in February 2006, with all patients providing informed consent for the use of their tissues for scientific research. Another 40 tissue samples were obtained from randomly selected sorafenib-treated HCC patients to detect circRNA-SORE expression and for β-catenin immunohistochemistry staining. We reviewed pathology records to identify samples with confirmed HCC. Low and high circRNA-SORE expression levels were cut off by median expression values. The current study conformed to the principles of the Declaration of Helsinki and was approved by the Institutional Review Board of the Sir Run-Run Shaw Hospital.
Additional methods information can be found in Supplementary Materials.

Elevated levels of circRNA-SORE in sorafenib-resistant HCC
To characterize HCC sorafenib resistance in vitro, we established four HCC sorafenib-resistant (SR) cell lines (HepG2-SR, SKhep1-SR, Huh7-SR and LM3-SR) and evaluated their resistance using the Real-time Cell Analysis xCELLigence System. As shown in Fig. 1a, sorafenib-resistant cell lines were less susceptible to sorafenib-killing than wild-type cell lines and, therefore, exhibited a faster growth rate with sorafenib treatment. We previously found that circRNA-SORE (a circular RNA up-regulated in sorafenib-resistant HCC cells) was upregulated in sorafenib-resistant HCC (unpublished observation) and circRNA-SORE was consistently upregulated in the sorafenib-resistant cell lines (Fig. 1b). Moreover, we confirmed its upregulation in sorafenib-resistant cell-derived xenograft (CDX) and PDX models (Fig. 1c), indicating the importance of circRNA-SORE in the induction and maintenance of HCC sorafenib resistance.
According to the UCSC Genome Browser, circRNA-SORE is generated from back-splicing of the 7th and 8th exons of the TLE4 gene (Fig. 1d). We treated cells with actinomycin D (an RNA synthesis inhibitor) to confirm that circRNA-SORE was more stable than its host TLE4 mRNA (Fig. 1e). circRNA-SORE was resistant to RNase R (a linear RNA degrader), but linear TLE4 mRNA was not (Fig. 1f). We further validated the nature of circRNA-SORE by RT-PCR using convergent and divergent primers, which yielded a specific divergent band from cDNA, but not gDNA (Fig. 1g). Sanger sequencing of the PCR product showed the circRNA-SORE sequence and the sequence in circBase to be the same (Fig. 1g). The functions of a circRNA can be influenced by its cellular localization. Fluorescence in situ hybridization (FISH) for circRNA-SORE demonstrated its cytoplasmic localization (Fig. 1h).
Overall, we identified a circular RNA related to HCC sorafenib resistance, circRNA-SORE, that was localized in the cytoplasm and up-regulated in sorafenib-resistant cells.

circRNA-SORE is vital for sustaining sorafenib resistance
To clarify the function of circRNA-SORE in sorafenib resistance, we knocked down circRNA-SORE using siRNA targeting the junction region of the circRNA sequence. Importantly, sorafenib resistance was significantly compromised by circRNA-SORE knockdown in the HCC sorafenib-resistant cell lines (Fig. 2a,  b). Colony formation assays showed similar results (Fig. 2c), indicating a critical role for circRNA-SORE in HCC sorafenib resistance. Real-time cell analyses showed significantly reduced viability of circRNA-SORE knock-down sorafenib-resistant cells, but not of wild-type cells (Fig. 2d), indicating that circRNA-SORE likely functions in a specific sorafenib-resistant context. JC-1 assays further revealed increased apoptosis in sorafenib-resistant cells upon circRNA-SORE knock-down (Fig. 2e). Importantly, HCC patients with lower circRNA-SORE expression levels had better recurrence free survival and overall survival than patients with higher circRNA-SORE expression levels ( Fig. 2f, g). Thus circRNA-SORE is potentially important for HCC prognosis.
To elucidate the downstream mechanisms of circRNA-SORE, several recently reported sorafenibresistance-related genes were screened in parental and sorafenib-resistant HCC cell lines. Among these, βcatenin was consistently up-regulated in sorafenibresistant cell lines (Fig. 4a). Analysis of HCC patient RNA-seq expression data from The Cancer Genome Atlas (TCGA) project [18] indicated increased CTNNB1 expression in tumor compared to normal liver tissue (Fig. 4b). Importantly, survival analysis showed that patients with higher CTNNB1 expression had worse overall survival compared with those with lower CTNNB1 expression ( Fig. 4c with group cutoffs at 75%/25%). Furthermore, silencing circRNA-SORE in three HCC sorafenib-resistant cell lines showed that Wnt2b/β-catenin signals may be closely linked with circRNA-SORE function because Wnt2b and βcatenin were significantly down-regulated upon circRNA-SORE silencing (Fig. 4d). Moreover, we measured circRNA-SORE expression in our collection of HCC tumor samples and found that the immunohistochemistry score of β-catenin was significantly correlated with the expression of circRNA-SORE (Fig. 4e), indicating that circRNA-SORE could regulate Wnt2b/ β-catenin signals. Consistently, miR-103a-2-5p and miR-660-3p alone and in combination significantly decreased the expression of Wnt2b and β-catenin (Fig. 4f), and miR-103a-2-5p and miR-660-3p in combination suppressed the transcription of β-catenin downstream signals (Fig. 4g). To further assess the role of β-catenin in mediating circRNA-SORE functions, we performed rescue assays by over-expressing and silencing β-catenin. These assays showed that βcatenin could partially rescue the effect of circRNA-SORE (Fig. 4h, i), indicating that circRNA-SORE sustains sorafenib resistance by regulating Wnt2b/β-catenin signals.
Our results show that circRNA-SORE can specifically sponge miR-103a-2-5p and miR-660-3p to act as a ceRNA, thereby competitively activating the Wnt/βcatenin pathway to promote sorafenib resistance.
(See figure on previous page.) Fig. 3 circRNA-SORE sustains sorafenib resistance by acting as an miRNA sponge for miR-103a-2-5p and miR-660-3p. a Predicted miRNA interactions with circRNA-SORE using Arraystar's miRNA target prediction software based on TargetScan & miRanda. b Map of circRNA-SORE and its predicted miRNAs interactions. c qPCR analysis for predicted miRNAs pulled-down by circRNA-SORE. d CCK-8 assays for HepG2-SR cells transfected with miRNA mimics and controls. e CCK-8 assays for HepG2-SR and SKhep1-SR cells transfected with miR-103a-2-5p, miR-660-3p mimics and controls with or without sorafenib treatment. f The co-localization of circRNA-SORE (red) and miR-103a-2-5p/miR-660-3p (green) was observed by fluorescence in situ hybridization (FISH) in HepG2-SR cells. Cell nuclei were counterstained with DAPI (blue). Scale bar: 5 μm. g Left: qPCR analysis of circRNA-SORE in the pull-down lysates of biotinylated miR-103a-2-5p or miR-660-3p. Right: workflow for the biotinylated miRNA pull-down assays. h Ago2 radioimmunoprecipitation assays for circRNA-SORE in wildtype HCC cells and sorafenib-resistant cells. i-j Dual Renilla and firefly luciferase reporter assays using the linear form of wild-type and mutant circRNA-SORE in HepG2-SR cells transfected with NC or miR-103a-2-5p/miR-660-3p miRNA mimics. Three different independent experiments with three technical repetitions were performed. Data are expressed as the mean ± SEM. Statistical analyses used Student's t-test, and p < 0.05 was considered statistically significant. * p < 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001. NS; not statistically significant  [19,20]. In particular, m 6 Acontaining circRNAs are reported to be subjected to endoribonucleolytic cleavage by YTHDF2-HRSP12-RNase P/MRP [10]. The modification preferentially occurs in the consensus motif 'RRm 6 ACH' (R = G or A; H = A, C or U) [21]. Using SRAMP and RMBase v2.0 to predict m 6 A sites, we found one m 6 A site in circRNA-SORE that is close to the junction region (Fig. 5a). It is possible that m 6 A modification plays important roles in the regulation of circRNA-SORE during the acquisition of sorafenib resistance.
In the complex mechanisms of m 6 A regulation, writers (m6A methyltransferases, such as METTL3 and METT L14) and erasers (m6A demethylases, such as FTO and ALKBH5) determine the distribution of m6A on RNA, whereas readers (such as YTHDF1 and YTHDF2) mediate m6A-dependent functions [22]. RNA pulldown assays indicated that circRNA-SORE interacted with m 6 A modifiers, such as METTL3, FTO and YTHDF1/2 (Fig. 5b). Importantly, m 6 A-specific immunoprecipitation assays showed increased m 6 A levels of circRNA-SORE in sorafenib-resistant HepG2 cells compared with wild-type controls, while m 6 A levels were reduced when METTL3 or METTL14 were silenced and increased when FTO was silenced in HepG2-SR cells (Fig. 5c-f). To further verify its existence, we designed a set of sequence-specific morpholino antisense oligos (MAOs) targeting this m 6 A site in circRNA-SORE (TGAGTCCT ACCTTGATGGAGTCTCT). The m 6 A levels were significantly decreased in HepG2-SR cells treated with an MAO targeting circRNA-SORE (MAO-SORE) compared with control MAO (MAO-NC)-treated cells (Fig. 5g). Site-directed m 6 A demethylation via programmable RNA m 6 A editing by CRISPR-Cas9-ALKBH5 conjugates [23] further confirmed the existence of this m 6 A site on circRNA-SORE (Fig. 5h). Moreover, over-expression of circRNA-SORE with the m 6 A site mutated produced a lower m 6 A level than over-expression of wild-type circRNA-SORE in HepG2 cells (Fig. 5i). These results indicated that circRNA-SORE is regulated by m6A modification during the acquisition of sorafenib resistance.
m6A modification is implicated in all aspects of posttranscriptional RNA metabolism, including mRNA localization, splicing, translation and degradation, which in turn regulates important biological processes [24]. We found that the levels of circRNA-SORE were reduced by silencing adenosine methyltransferase METTL3/14 (Fig. 5j). Notably, MAO-SORE could reduce circRNA-SORE levels in sorafenib-resistant HCC cells, but not in wild-type HCC cells (Fig. 5k, l), further indicating that circRNA-SORE levels were promoted by the m 6 A modification. Intriguingly, RNA stability assays demonstrated that MAO-SORE reduced the stability of circRNA-SORE compared with MAO-NC in HepG2-SR cells (Fig. 5m), while wild-type circRNA-SORE was more stable than m 6 A-mutated circRNA-SORE under actinomycin D treatment (Fig. 5n). These results indicate that the m 6 A modification regulates circRNA-SORE expression by increasing its stability.

Targeting circRNA-SORE in vivo retards sorafenib-resistant HCC
To further support the in vitro findings and to discover potential clinical applications, we employed in vivo sorafenib-resistant models. Orthotopic implantations of SKhep1-SR cells expressing sh-circRNA-SORE were significantly more sensitive to sorafenib treatment than controls, while the two groups showed similar responses to solvent treatment (Fig. 6a). Moreover, local injection with sh-circRNA-SORE lentivirus around the subcutaneous implantation site in sorafenib-resistant LM3-CDX (See figure on previous page.) Fig. 4 miR-103a-2-5p and miR-660-3p regulate Wnt2b/β-catenin signaling. a Western blotting for β-catenin and βtubulin in three sorafenib-resistant cell lines compared with their parental control cell lines. b Box-plot (left) and gene expression level (right) of CTNNB1 expression in TCGA HCC tumor and matched TCGA normal liver tissues along with GTEx data. c Overall survival analysis of HCC patients (N = 182) from TCGA project with high or low CTNNB1 expression levels (defined by RNA sequencing with group cutoff in 75%/25% quartile). d Western blot for Wnt2b, β-catenin and β-tubulin in three sorafenib-resistant cell lines with or without si-circRNA-SORE. e Correlation of relative RNA expression of circRNA-SORE and protein expression of CTNNB1 by immunohistochemistry. The slides were reviewed and scored by an experienced pathologist without the knowledge of patient outcome. f Western blot for Wnt2b, β-catenin, β-tubulin with or without miR-103a-2-5p/miR-660-3p mimics in HepG2-SR. g qPCR analysis for downstream genes of β-catenin. h Western blot for β-catenin and β-tubulin (upper) and CCK8 assays (lower) for HepG2-SR cells transfected with si-circRNA-SORE and over-expressing β-catenin. i Western blot for β-catenin and β-tubulin (upper) and CCK8 assays (lower) for HepG2 cells transfected with si-β-catenin and over-expressing circRNA-SORE. j-k Renilla and firefly dual luciferase reporter assay using Wnt2b 3′-UTR constructs containing wild-type or mutant seed sequences in HepG2-SR cells transfected with negative control (NC) or miR-103a-2-5p/miR-660-3p miRNA mimics. Three different independent experiments with three technical repetitions were performed. Data are expressed as the mean ± SEM. Statistical analyses used Student's t-test, and p < 0.05 was considered statistically significant. * p < 0.05 and ** p < 0.01. NS; not statistically significant models significantly increased the sensitivity of HCC cells to sorafenib treatment ( Fig. 6b-d). Similar to in vitro results, immunohistochemistry also showed impaired expression of Wnt2b and β-catenin in the sh-circRNA-SORE group (Fig. 6e). These findings support targeting of circRNA-SORE in sorafenib-treated HCC patients as a novel pharmaceutical intervention for advanced HCC.

Discussion
Many HCC patients are diagnosed too late for surgery because early stages of HCC present no clear symptoms. Other interventions, including molecular targeted therapy, are limited by a lack of efficacy. Unfortunately, patients who qualify for liver resection or transplantation have a high incidence of recurrence and metastases. Therefore, there is an urgent unmet medical need for the development of life-prolonging therapies. As the first US Food and Drug Administration-approved molecular targeted drug, sorafenib provided a 3 month prolongation of the median overall survival time [25,26] and transient and limited efficacy of sorafenib was frequently reported in HCC patients. Adverse effects, such as rash, diarrhea, high blood pressure, and hand-foot syndrome, also limit high-dosage use of sorafenib [27]. Moreover, the STORM trial [28] suggested that HCC patients who accepted radical treatment (resection or ablation) could not benefit from sorafenib treatment in the adjuvant setting, posing more challenges to its clinical application. Sorafenib resistance prompts the need for new therapies to overcome resistance [29]. Multiple mechanisms underlying impaired sensitivity to sorafenib in HCC have been investigated, including Wnt/β-catenin, TGFβ, Ras/ MEK/ERK, PI3K/Akt, TNFα/NF-κB, and JAK/STAT pathways, autophagy, epithelial-mesenchymal transition, cancer stem cells, tumor microenvironment, and epigenetic regulation (involving miR-222, miR-494, miR-21 and miR122) [3]. To overcome sorafenib resistance and lower its onset concentration, efforts were made to develop combined therapies [30][31][32][33][34]. Yet, the overall outcomes of liver cancer are still far from satisfactory. In the present study, sorafenib-resistant cell lines and animal models were developed to simulate the sorafenib resistance in HCC patients.
Non-coding RNAs are functional RNAs transcribed from the genome from which proteins cannot usually be translated. They mainly include miRNAs, long non-coding RNAs (lncRNAs) and circRNAs. Recent studies indicate that ncRNAs play vital roles in diverse biological and pathological processes, including cancer [35]. First described in 1993, miRNAs are small noncoding RNA molecules that play crucial post-transcriptional regulatory roles. Importantly, miR-NAs are involved in the development of sorafenib resistance through complementary base pairing with mRNAs, predominantly in the 3′-UTR. For example, miRNA-216a/217-induced epithelial-mesenchymal transition promotes sorafenib resistance and liver cancer recurrence by targeting PTEN and SMAD7 [36]. However, its small linear structure has inherent limitations, including extensive but limited effect, poor specificity, instability, and off-target effects. In contrast to miRNAs, circRNAs have a covalently-closed loop structure with neither 5′ to 3′ polarity nor a polyadenylated tail, which makes them more stable than their linear counterparts and more resistant to RNase R degradation.
circRNAs were occasionally identified more than 20 years ago and were thought to be of low abundance and to result from alternative splicing errors during transcription. Using high-throughput sequencing and novel (See figure on previous page.) Fig. 5 circRNA-SORE is modulated by m 6 A RNA methylation. a Predicted m 6 A site in circRNA-SORE from overlapping results of a sequence-based N6-methyladenosine (m 6 A) modification site predictor (SRAMP) and RMBase v2.0. b RNA pulldown assays and western blotting for YTHDF1/2, METTL3, and FTO using a circRNA-SORE probe. c Flow diagram of m6A-specific immunoprecipitation (MeRIP) assays. d MeRIP assays for m 6 A-modified circRNA-SORE in HepG2-SR and HepG2-WT cell lines. e MeRIP assays for m 6    computational approaches, circRNAs derived from exons or introns were determined to be widespread and diverse endogenous eukaryotic ncRNAs participating in various normal and disease-related processes [37,38].
Recent findings indicate that circRNAs function through RNA interactions, protein interactions, or by serving as transcription or splicing regulators [15]. For example, the circRNA, CDR1as, harbors more than 60 conserved binding sites for miR-7 [38]. Additionally, circRNAs may be associated with various miRNAs, as demonstrated for circHIPK3, which can bind to multiple miRNAs [39].
Recently, particular functions of circRNAs have been revealed in HCC; circFBLIM1 and circ-FOXP1 can act as ceRNAs to promote HCC progression [40,41], while cir-cMTO1 can act as a sponge of microRNA-9 to suppress HCC progression [41]. Also, study of circ_0003418 showed that circRNAs can influence anti-tumor treatments such as cisplatin [9]. Containing multiple binding sites for particular miRNAs, circRNAs have specific and efficient functionality. Furthermore, tissue-and stagespecific expression makes circRNAs potential targets for clinical intervention. Also, because of their RNA stability, circRNAs are valued as candidates for non-invasive biomarkers. circ_0005075, circ_0016788, ciRS7, circ_ 0128298, circ_0091579, and circ-CDYL have been identified as potential diagnostic biomarkers for HCC [42]. However, the relationship between circRNAs and resistance of HCC to sorafenib has not been reported. Thus, we aimed to fully understand the circRNA landscape when HCC patients acquired sorafenib resistance during treatment. We identified consistently up-regulated circRNA-SORE in sorafenib resistant cell lines and CDX and PDX models that was critical for the maintenance of sorafenib resistance. Importantly, orthotopic in vitro developed sorafenib-resistant CDX model and subcutaneous in vivo developed sorafenib-resistant CDX model were applied in the present study, showing that circRNA-SORE silencing could effectively reverse the acquired sorafenib resistance and retard tumor progression. In particular, we found specific mechanisms by which circRNA-SORE could specifically bind miR-103a-2-5p and miR-660-3p to act as a miRNA sponge to competitively activate Wnt/β-catenin and induce sorafenib resistance, thereby identifying a potential biomarker for prediction of sorafenib resistance and a promising therapeutic target for HCC. Derived from the 7th and 8th exons of TLE4, circRNA-SORE is generated into a circular form by back-splicing [43]. It was recently reported that exon-derived circRNAs are predominantly located in the cytoplasm. In our study, both nucleocytoplasmic fractionation and circRNA-FISH assays revealed that circRNA-SORE located predominantly in the cytoplasm of HCC cells. m6A modification was first discovered in the 1970s but the development of RNA-sequencing techniques and clarification of the proteins involved make a deeper understanding of the process possible. m 6 A depends on m 6 A writers (METTL3, METTL14, KIAA1429, WTAP, RBM15 and ZC3H13), erasers (FTO and ALKBH5) and readers (YTHDC1, YTHDC2, YTHDF1, YTHDF2 and HNRNPC) [44]. m 6 A modifications may play important roles in RNA production, stability and interactions in cancers [45]. Metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), a conserved lncRNA, is highly methylated with m 6 A. Two of these m 6 A residues can prevent RNA local secondary structure formation and enhance the recognition and binding of hnRNPC to a U5-tract in the MALAT1 hairpin through an "m6A switch" mechanism [46]. hnRNPA2B1 can recognize pri-miRNAs with m 6 A marks to promote interaction of DGR8 and pri-miRNAs and miRNA processing [47]. However, in HCC, the roles of m 6 A modification are intricate and controversial. Knock-down of YTHDF2 suppresses proliferation of HCC cells [48]; however, YTHD F2 was also reported to suppress HCC tumors by targeting EGFR, IL11 and SRPINE2 [49]. Indeed, the detailed roles of m 6 A modification in HCC need further investigation. Recently, m 6 A-modified circRNAs were found with cell-type-specific expression [45]. m 6 A can recruit YTHDF3 and initiation factor, eIF4G2, to regulate protein synthesis from circRNAs [19]. The change of m 6 A level by ALKBH5 and METTL3 affects circRNA (See figure on previous page.) Fig. 6 Targeting circRNA-SORE in vivo retards sorafenib-resistant HCC. a In vivo luminescent imaging of orthotopically implanted SKhep1-SR-luc or SKhep1-SR-luc-sh-cricRNA-SORE cells in 4-6 week-old male BALB/c nude mice treated with solvent or sorafenib (30 mg/kg/day) for 4 weeks from the fourth week after implantation. Luminescence intensity ranges from low (blue) to high (red). Tumor burdens were quantified by total photon flux (p/s). b-d Flow Chart of LM3-CDX-SR model construction. The first CDX generation was constructed in 4-6 week-old male BALB/c nude mice and treated with sorafenib (30 mg/kg daily, oral gavage). Twelve weeks later, the most resistant xenograft was disaggregated and implanted subcutaneously into 4-6 week-old BALB/c nude mice as the second SR-CDX. Four weeks after implantation, the second SR-CDX mice were treated with sorafenib (30 mg/kg daily, oral gavage) and locally injected with sh-circRNA-SORE lentivirus or its negative control (twice a week for 2 weeks). Mice were euthanized on the sixth week and tumors were isolated and the tumor weights were measured. The subcutaneous tumor size was measured and recorded every 2 days using the Vernier caliper as follows: tumor volume (mm 3 ) = (L × W 2 )/2, where L is the long axis and W the short axis. e Immunohistochemistry for Wnt2b and β-catenin comparing sh-circRNA-SORE and control CDX. Scale bar, 50 μm. f Hypothetical model for circRNA-SORE function in HCC. circRNA-SORE acts as a sponge for miR-103a-2-5p and miR-660-3p to promote HCC sorafenib resistance by regulating the Wnt2b/β-catenin signaling pathway biosynthesis in spermatogenesis by regulating enhanced splicing and promotes circRNA formation [50]. YTHD F2-HRSP12-RNase P/MRP-mediated endoribonucleolytic cleavage is related to m 6 A-containing circRNA decay [51]. However, the functions of m 6 A-modified circRNAs in HCC and sorafenib resistance were still elusive.
In this study, we identified a predicted m 6 A site in circRNA-SORE by a series of experiments. Our results show that the m 6 A level of circRNA-SORE is increased in sorafenib-resistant cells, and that the expression of circRNA-SORE is decreased when its m 6 A modification was inhibited. Recently, lncRNA GAS5-AS was found to enhance GAS5 stability by interacting with ALKBH5 and regulating m 6 A modifications of GAS5, which was dependent on ALKBH5 and YTHDF2 [52]. In the present study, mechanistic dissection suggested that the m 6 A modification can stabilize circRNA-SORE. However, no significant differences in the expression level of m 6 A-related proteins were found between sorafenibresistant and parental cells. The mechanism of increased circRNA-SORE m 6 A levels in sorafenib-resistant cells requires further investigation.

Conclusion
To conclude, upregulated circRNA-SORE was necessary for the maintenance of sorafenib resistance, and silencing circRNA-SORE substantially increased the efficacy of sorafenib-induced apoptosis. Mechanistic studies revealed that circRNA-SORE sequestered miR-103a-2-5p and miR-660-3p by acting as a miRNA sponge, thereby competitively activating the Wnt/β-catenin pathway and inducing sorafenib resistance. Importantly, the increased levels of circRNA-SORE in sorafenib-resistant cells resulted from increased RNA stability, which was triggered by increased levels of m 6 A at specific sites in circRNA-SORE. In vivo delivery of circRNA-SORE RNAi by local short hairpin RNA lentivirus injection substantially enhanced sorafenib efficacy in animal models. This work describes a novel mechanism of sorafenib resistance maintenance and is a proof-of-concept demonstration for targeting circRNA-SORE in sorafenib-treated HCC patients as a novel pharmaceutical intervention for advanced HCC.