- Open Access
N6-methyladenosine-modified CircRNA-SORE sustains sorafenib resistance in hepatocellular carcinoma by regulating β-catenin signaling
Molecular Cancer volume 19, Article number: 163 (2020)
Background and aims
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.
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 . 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 . 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 epithelial-mesenchymal transition, the proliferation of cancer stem cells, and metabolic reprogramming, the exact cause is still elusive . 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, cis-regulation of transcription, and by acting as competing endogenous RNAs (ceRNAs) . Importantly, many circRNAs 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 , while circRHOT1 can inhibit HCC development and progression via recruiting TIP60 to initiate NR2F6 expression . Circ_0003418 increases the sensitivity of HCC cells to cisplatin by inhibiting the Wnt/β-catenin pathway . However, the roles of circRNAs in HCC sorafenib resistance remain unknown.
N6-methyladenosine (m6A) 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 . The m6A modification has been implicated in various cellular and physiological events, including carcinogenesis . 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 . 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 m6A modification of mRNA, which promotes cell proliferation . In HCC, m6A modification also regulates the progress of oncogenesis. SIRT1 promotes HCC by increasing the overall m6A modification level to attenuate the expression of the tumor suppressor, GNAO1 . However, the functions of m6A modification in HCC sorafenib resistance remain elusive. Moreover, the roles of m6A-modified 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 N6-methyladenosine (m6A) 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.
Materials and methods
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′-GAGTTGTTGCTGCTTGATGGAGTC-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 RiboBio (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.
The sh-circRNA-SORE sequence (same as the sicircRNA-SORE sequence) was cloned into GV248. The backbone of pLCDH-ciR (Geneseed, Guangzhou, China) was used to construct the circRNA-SORE over-expression plasmid (OE-circRNA-SORE). The linear sequence of circRNA-SORE together with the flanking circRNA-inducing sequence was amplified and subcloned into EcoRI and BamHI sites of pLCDH-ciR, and termed pLCDH-ciR-SORE. The primers used were as follows:
pLCDH-circRNA-SORE-M has a mutated circRNA-SORE m6A 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. Forty-eight hours after transfection, viruses were collected and concentrated by density gradient centrifugation, then frozen at − 80 °C for later use.
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 knock-down 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.
Our results show that circRNA-SORE is essential for sustaining sorafenib resistance and that silencing circRNA-SORE substantially increased the efficacy of sorafenib by inducing tumor cell apoptosis.
circRNA-SORE sustains sorafenib resistance by acting as an miRNA sponge for miR-103a-2-5p and miR-660-3p
circRNAs can function in numerous ways, including by sponging miRNAs, interacting with proteins, promoting gene transcription, and encoding peptides . Given that circRNA-SORE localizes to the cytoplasm, as described above, we hypothesized that miRNA sponge activity could be a possible mechanism for its functional effects. circRNA-miRNA interactions were predicted using Arraystar’s miRNA target prediction software based on TargetScan  & miRanda  (Fig. 3a, b). Among these predicted miRNAs, miR-96-3p, miR-103a-2-5p, miR-185-3p, miR-328-5p, miR-345-3p, miR-591, miR-609, and miR-660-3p were found to be associated with circRNA-SORE (Fig. 3c). By transfecting mimics of these miRNAs, we found that miR-103a-2-5p and miR-660-3p reduced sorafenib resistance in HepG2-SR cells (Fig. 3d, e). RNA-FISH for circRNA-SORE further revealed the colocalization of circRNA-SORE with miR-103a-2-5p and miR-660-3p (Fig. 3f). Consistently, biotinylated miRNA pull-down assays showed significant enrichment of circRNA-SORE with miR-103a-2-5p and miR-660-3p (Fig. 3g). Ago2 immunoprecipitation assays detected circRNA-SORE in the RNA-induced silencing complex (RISC) in sorafenib-resistant cells and not in wild-type cells (Fig. 3h). Moreover, according to the predicted miRNA binding sites on circRNA-SORE, luciferase reporter assays of wild-type and mutant linear forms of circRNA-SORE showed that miR-103a-2-5p and miR-660-3p could specifically target the wild-type linear form of circRNA-SORE resulting in decreased Renilla luciferase activity, but not the mutant 1 or 2 forms (Fig. 3i, j). These results suggested a miRNA sponging role of circRNA-SORE for miR-103a-2-5p and miR-660-3p in sorafenib-resistant HCC.
To elucidate the downstream mechanisms of circRNA-SORE, several recently reported sorafenib-resistance-related genes were screened in parental and sorafenib-resistant HCC cell lines. Among these, β-catenin was consistently up-regulated in sorafenib-resistant cell lines (Fig. 4a). Analysis of HCC patient RNA-seq expression data from The Cancer Genome Atlas (TCGA) project  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.
We further hypothesized that miR-103a-2-5p and miR-660-3p inhibit the expression of Wnt2b by targeting its 3′-untranslated region (3′-UTR). Reporter assays demonstrated that both miR-103a-2-5p and miR-660-3p decreased the luciferase activity of the wild-type reporter, while no significant changes were found using mutant reporters (Fig. 4j, k).
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.
circRNA-SORE is modulated by m6A RNA methylation
circRNA-SORE is likely generated through pre-mRNA back-splicing of the 7th and 8th exons of the TLE4 transcript. However, the detailed mechanisms controlling circRNA-SORE levels remain elusive. m6A is the most abundant base modification of RNA and it might play a role in circRNA biogenesis [19, 20]. In particular, m6A-containing circRNAs are reported to be subjected to endoribonucleolytic cleavage by YTHDF2-HRSP12-RNase P/MRP . The modification preferentially occurs in the consensus motif ‘RRm6ACH’ (R = G or A; H = A, C or U) . Using SRAMP and RMBase v2.0 to predict m6A sites, we found one m6A site in circRNA-SORE that is close to the junction region (Fig. 5a). It is possible that m6A modification plays important roles in the regulation of circRNA-SORE during the acquisition of sorafenib resistance.
In the complex mechanisms of m6A regulation, writers (m6A methyltransferases, such as METTL3 and METTL14) 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 . RNA pulldown assays indicated that circRNA-SORE interacted with m6A modifiers, such as METTL3, FTO and YTHDF1/2 (Fig. 5b). Importantly, m6A-specific immunoprecipitation assays showed increased m6A levels of circRNA-SORE in sorafenib-resistant HepG2 cells compared with wild-type controls, while m6A 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 m6A site in circRNA-SORE (TGAGTCCTACCTTGATGGAGTCTCT). The m6A 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 m6A demethylation via programmable RNA m6A editing by CRISPR-Cas9-ALKBH5 conjugates  further confirmed the existence of this m6A site on circRNA-SORE (Fig. 5h). Moreover, over-expression of circRNA-SORE with the m6A site mutated produced a lower m6A 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 post-transcriptional RNA metabolism, including mRNA localization, splicing, translation and degradation, which in turn regulates important biological processes . 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 m6A 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 m6A-mutated circRNA-SORE under actinomycin D treatment (Fig. 5n). These results indicate that the m6A 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 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.
Overall, our results show that circRNA-SORE can sequester 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, an increased level of N6-methyladenosine (m6A) at a specific circRNA-SORE site increased RNA stability and upregulated levels of circRNA-SORE in sorafenib-resistant cells (Fig. 6f).
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 . Moreover, the STORM trial  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 . 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) . 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 . First described in 1993, miRNAs are small noncoding RNA molecules that play crucial post-transcriptional regulatory roles. Importantly, miRNAs 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 . 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 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 . For example, the circRNA, CDR1as, harbors more than 60 conserved binding sites for miR-7 . Additionally, circRNAs may be associated with various miRNAs, as demonstrated for circHIPK3, which can bind to multiple miRNAs . 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 circMTO1 can act as a sponge of microRNA-9 to suppress HCC progression . Also, study of circ_0003418 showed that circRNAs can influence anti-tumor treatments such as cisplatin . Containing multiple binding sites for particular miRNAs, circRNAs have specific and efficient functionality. Furthermore, tissue- and stage-specific 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 . 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 . 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. m6A depends on m6A writers (METTL3, METTL14, KIAA1429, WTAP, RBM15 and ZC3H13), erasers (FTO and ALKBH5) and readers (YTHDC1, YTHDC2, YTHDF1, YTHDF2 and HNRNPC) . m6A modifications may play important roles in RNA production, stability and interactions in cancers . Metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), a conserved lncRNA, is highly methylated with m6A. Two of these m6A 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 . hnRNPA2B1 can recognize pri-miRNAs with m6A marks to promote interaction of DGR8 and pri-miRNAs and miRNA processing . However, in HCC, the roles of m6A modification are intricate and controversial. Knock-down of YTHDF2 suppresses proliferation of HCC cells ; however, YTHDF2 was also reported to suppress HCC tumors by targeting EGFR, IL11 and SRPINE2 . Indeed, the detailed roles of m6A modification in HCC need further investigation. Recently, m6A-modified circRNAs were found with cell-type-specific expression . m6A can recruit YTHDF3 and initiation factor, eIF4G2, to regulate protein synthesis from circRNAs . The change of m6A level by ALKBH5 and METTL3 affects circRNA biosynthesis in spermatogenesis by regulating enhanced splicing and promotes circRNA formation . YTHDF2–HRSP12–RNase P/MRP-mediated endoribonucleolytic cleavage is related to m6A-containing circRNA decay . However, the functions of m6A-modified circRNAs in HCC and sorafenib resistance were still elusive.
In this study, we identified a predicted m6A site in circRNA-SORE by a series of experiments. Our results show that the m6A level of circRNA-SORE is increased in sorafenib-resistant cells, and that the expression of circRNA-SORE is decreased when its m6A modification was inhibited. Recently, lncRNA GAS5-AS was found to enhance GAS5 stability by interacting with ALKBH5 and regulating m6A modifications of GAS5, which was dependent on ALKBH5 and YTHDF2 . In the present study, mechanistic dissection suggested that the m6A modification can stabilize circRNA-SORE. However, no significant differences in the expression level of m6A-related proteins were found between sorafenib-resistant and parental cells. The mechanism of increased circRNA-SORE m6A levels in sorafenib-resistant cells requires further investigation.
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 m6A 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.
Availability of data and materials
All data generated or analyzed during this study are included either in this article or in the supplementary information files.
Sorafenib Hepatocellular Carcinoma Assessment Randomized Protocol
competing endogenous RNAs
Acute lymphoblastic leukemia
Hepatitis B X-interacting protein
- a circular RNA up-regulated in sorafenib-resistant HCC cells:
Fluorescence in situ hybridization
RNA-induced silencing complex
The Cancer Genome Atlas
long non-coding RNAs
Metastasis-associated lung adenocarcinoma transcript 1
Forner A, Reig M, Bruix J. Hepatocellular carcinoma. Lancet. 2018;391(10127):1301–14.
Connell LC, Harding JJ, Abou-Alfa GK. Advanced hepatocellular Cancer: the current state of future research. Curr Treat Options in Oncol. 2016;17(8):43.
Xia S, Pan Y, Liang Y, Xu J, Cai X. The microenvironmental and metabolic aspects of sorafenib resistance in hepatocellular carcinoma. EBioMedicine. 2020;51:102610.
Lasda E, Parker R. Circular RNAs: diversity of form and function. RNA (New York). 2014;20(12):1829–42.
Zhang L, Zhou Q, Qiu Q, Hou L, Wu M, Li J, Li X, Lu B, Cheng X, Liu P, et al. CircPLEKHM3 acts as a tumor suppressor through regulation of the miR-9/BRCA1/DNAJB6/KLF4/AKT1 axis in ovarian cancer. Mol Cancer. 2019;18(1):144.
Li Z, Ruan Y, Zhang H, Shen Y, Li T, Xiao B. Tumor-suppressive circular RNAs: mechanisms underlying their suppression of tumor occurrence and use as therapeutic targets. Cancer Sci. 2019;110(12):3630-8.
Su Y, Lv X, Yin W, Zhou L, Hu Y, Zhou A, Qi F. CircRNA Cdr1as functions as a competitive endogenous RNA to promote hepatocellular carcinoma progression. Aging. 2019;11(19):8182–203.
Wang L, Long H, Zheng Q, Bo X, Xiao X, Li B. Circular RNA circRHOT1 promotes hepatocellular carcinoma progression by initiation of NR2F6 expression. Mol Cancer. 2019;18(1):119.
Chen H, Liu S, Li M, Huang P, Li X. circ_0003418 inhibits tumorigenesis and Cisplatin Chemoresistance through Wnt/β-catenin pathway in hepatocellular carcinoma. OncoTargets Ther. 2019;12:9539–49.
Lee Y, Choe J, Park OH, Kim YK. Molecular mechanisms driving mRNA degradation by m (6) a modification. Trends Genet. 2020;36(3):177–88.
Choe J, Lin S, Zhang W, Liu Q, Wang L, Ramirez-Moya J, Du P, Kim W, Tang S, Sliz P, et al. mRNA circularization by METTL3-eIF3h enhances translation and promotes oncogenesis. Nature. 2018;561(7724):556–60.
Vu LP, Pickering BF, Cheng Y, Zaccara S, Nguyen D, Minuesa G, Chou T, Chow A, Saletore Y, MacKay M, et al. The N (6)-methyladenosine (m (6) a)-forming enzyme METTL3 controls myeloid differentiation of normal hematopoietic and leukemia cells. Nat Med. 2017;23(11):1369–76.
Cai X, Wang X, Cao C, Gao Y, Zhang S, Yang Z, Liu Y, Zhang X, Zhang W, Ye L. HBXIP-elevated methyltransferase METTL3 promotes the progression of breast cancer via inhibiting tumor suppressor let-7g. Cancer Lett. 2018;415:11–9.
Liu X, Liu J, Xiao W, Zeng Q, Bo H, Zhu Y, Gong L, He D, Xing X, Li R, et al. SIRT1 regulates N (6) -methyladenosine RNA modification in hepatocarcinogenesis by inducing RANBP2-dependent FTO SUMOylation. Hepatology. 2020.
Lyu D, Huang S. The emerging role and clinical implication of human exonic circular RNA. RNA Biol. 2017;14(8):1000-6.
Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in drosophila. Genome Biol. 2003;5(1):R1.
Pasquinelli AE. MicroRNAs and their targets: recognition, regulation and an emerging reciprocal relationship. Nat Rev Genet. 2012;13(4):271–82.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98-102.
Yang Y, Fan X, Mao M, Song X, Wu P, Zhang Y, Jin Y, Yang Y, Chen LL, Wang Y, et al. Extensive translation of circular RNAs driven by N (6)-methyladenosine. Cell Res. 2017;27(5):626–41.
Chen RX, Chen X, Xia LP, Zhang JX, Pan ZZ, Ma XD, Han K, Chen JW, Judde JG, Deas O, et al. N (6)-methyladenosine modification of circNSUN2 facilitates cytoplasmic export and stabilizes HMGA2 to promote colorectal liver metastasis. Nat Commun. 2019;10(1):4695.
Csepany T, Lin A, Baldick CJ Jr, Beemon K. Sequence specificity of mRNA N6-adenosine methyltransferase. J Biol Chem. 1990;265(33):20117–22.
Fu Y, Dominissini D, Rechavi G, He C. Gene expression regulation mediated through reversible m (6) a RNA methylation. Nat Rev Genet. 2014;15(5):293–306.
Liu XM, Zhou J, Mao Y, Ji Q, Qian SB. Programmable RNA N (6)-methyladenosine editing by CRISPR-Cas9 conjugates. Nat Chem Biol. 2019;15(9):865–71.
Meyer KD, Jaffrey SR. The dynamic epitranscriptome: N6-methyladenosine and gene expression control. Nat Rev Mol Cell Biol. 2014;15(5):313–26.
Cheng AL, Kang YK, Chen Z, Tsao CJ, Qin S, Kim JS, Luo R, Feng J, Ye S, Yang TS, et al. Efficacy and safety of sorafenib in patients in the Asia-Pacific region with advanced hepatocellular carcinoma: a phase III randomised, double-blind, placebo-controlled trial. Lancet Oncol. 2009;10(1):25–34.
Llovet JM, Ricci S, Mazzaferro V, Hilgard P, Gane E, Blanc JF, de Oliveira AC, Santoro A, Raoul JL, Forner A, et al. Sorafenib in advanced hepatocellular carcinoma. N Engl J Med. 2008;359(4):378–90.
Hampton T. Cancer drug trials show modest benefit: drugs target liver, gastric, head and neck cancers. JAMA. 2007;298(3):273–5.
Bruix J, Takayama T, Mazzaferro V, Chau GY, Yang J, Kudo M, Cai J, Poon RT, Han KH, Tak WY, et al. Adjuvant sorafenib for hepatocellular carcinoma after resection or ablation (STORM): a phase 3, randomised, double-blind, placebo-controlled trial. Lancet Oncol. 2015;16(13):1344–54.
Chen J, Jin R, Zhao J, Liu J, Ying H, Yan H, Zhou S, Liang Y, Huang D, Liang X, et al. Potential molecular, cellular and microenvironmental mechanism of sorafenib resistance in hepatocellular carcinoma. Cancer Lett. 2015;367(1):1–11.
Xu J, Zheng L, Chen J, Sun Y, Lin H, Jin R-A, Tang M, Liang X, Cai X. Increasing AR by HIF-2α inhibitor (PT-2385) overcomes the side-effects of sorafenib by suppressing hepatocellular carcinoma invasion via alteration of pSTAT3, pAKT and pERK signals. Cell Death Dis. 2017;8:e3095.
Ji L, Lin Z, Wan Z, Xia S, Jiang S, Cen D, Cai L, Xu J, Cai X. miR-486-3p mediates hepatocellular carcinoma sorafenib resistance by targeting FGFR4 and EGFR. Cell Death Dis. 2020;11(4):250.
Lin Z, Xia S, Liang Y, Ji L, Pan Y, Jiang S, Wan Z, Tao L, Chen J, Lin C, et al. LXR activation potentiates sorafenib sensitivity in HCC by activating microRNA-378a transcription. Theranostics. 2020;10(19):8834–50.
Xu J, Lin H, Li G, Sun Y, Chen J, Shi L, Cai X, Chang C. The miR-367-3p increases Sorafenib chemotherapy efficacy to suppress hepatocellular carcinoma metastasis through altering the androgen receptor signals. EBioMedicine. 2016;12:55–67.
Xu J, Lin H, Li G, Sun Y, Shi L, Ma WL, Chen J, Cai X, Chang C. Sorafenib with ASC-J9((R)) synergistically suppresses the HCC progression via altering the pSTAT3-CCL2/Bcl2 signals. Int J Cancer. 2017;140(3):705–17.
Anastasiadou E, Jacob LS, Slack FJ. Non-coding RNA networks in cancer. Nat Rev Cancer. 2018;18(1):5–18.
Xia H, Ooi LL, Hui KM. MicroRNA-216a/217-induced epithelial-mesenchymal transition targets PTEN and SMAD7 to promote drug resistance and recurrence of liver cancer. Hepatology. 2013;58(2):629–41.
Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, Maier L, Mackowiak SD, Gregersen LH, Munschauer M, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495(7441):333–8.
Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature. 2013;495(7441):384–8.
Zheng Q, Bao C, Guo W, Li S, Chen J, Chen B, Luo Y, Lyu D, Li Y, Shi G, et al. Circular RNA profiling reveals an abundant circHIPK3 that regulates cell growth by sponging multiple miRNAs. Nat Commun. 2016;7:11215.
Bai N, Peng E, Qiu X, Lyu N, Zhang Z, Tao Y, Li X, Wang Z. circFBLIM1 act as a ceRNA to promote hepatocellular cancer progression by sponging miR-346. J Ex Clin Cancer Res. 2018;37(1):172.
Han D, Li J, Wang H, Su X, Hou J, Gu Y, Qian C, Lin Y, Liu X, Huang M, et al. Circular RNA circMTO1 acts as the sponge of microRNA-9 to suppress hepatocellular carcinoma progression. Hepatology. 2017;66(4):1151–64.
Qiu L, Xu H, Ji M, Shang D, Lu Z, Wu Y, Tu Z, Liu H. Circular RNAs in hepatocellular carcinoma: biomarkers, functions and mechanisms. Life Sci. 2019;231:116660.
Hentze MW, Preiss T. Circular RNAs: splicing's enigma variations. EMBO J. 2013;32(7):923–5.
Yang Y, Hsu PJ, Chen YS, Yang YG. Dynamic transcriptomic m (6) a decoration: writers, erasers, readers and functions in RNA metabolism. Cell Res. 2018;28(6):616–24.
Zhou C, Molinie B, Daneshvar K, Pondick JV, Wang J, Van Wittenberghe N, Xing Y, Giallourakis CC, Mullen AC. Genome-wide maps of m6A circRNAs identify widespread and cell-type-specific methylation patterns that are distinct from mRNAs. Cell Rep. 2017;20(9):2262–76.
Liu N, Dai Q, Zheng G, He C, Parisien M, Pan T. N (6)-methyladenosine-dependent RNA structural switches regulate RNA-protein interactions. Nature. 2015;518(7540):560–4.
Alarcón CR, Goodarzi H, Lee H, Liu X, Tavazoie S, Tavazoie SF. HNRNPA2B1 is a mediator of m (6) A-dependent nuclear RNA processing events. Cell. 2015;162(6):1299–308.
Yang Z, Li J, Feng G, Gao S, Wang Y, Zhang S, Liu Y, Ye L, Li Y, Zhang X. MicroRNA-145 modulates N (6)-Methyladenosine levels by targeting the 3′-Untranslated mRNA region of the N (6)-Methyladenosine binding YTH domain family 2 protein. J Biol Chem. 2017;292(9):3614–23.
Zhong L, Liao D, Zhang M, Zeng C, Li X, Zhang R, Ma H, Kang T. YTHDF2 suppresses cell proliferation and growth via destabilizing the EGFR mRNA in hepatocellular carcinoma. Cancer Lett. 2019;442:252–61.
Tang C, Xie Y, Yu T, Liu N, Wang Z, Woolsey RJ, Tang Y, Zhang X, Qin W, Zhang Y, et al. m (6) A-dependent biogenesis of circular RNAs in male germ cells. Cell Res. 2020;30(3):211–28.
Park OH, Ha H, Lee Y, Boo SH, Kwon DH, Song HK, Kim YK. Endoribonucleolytic Cleavage of m (6) A-Containing RNAs by RNase P/MRP Complex. Mol Cell. 2019;74(3):494–507.e498.
Wang X, Zhang J, Wang Y. Long noncoding RNA GAS5-AS1 suppresses growth and metastasis of cervical cancer by increasing GAS5 stability. Am J Transl Res. 2019;11(8):4909–21.
The authors would like to thank the Jeremy Allen, PhD, from Liwen Bianji, Edanz Group China (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.
This research was supported by the National Natural Science Foundation of China under Grant No. 81772546 (to Cai X), No. 81827804 (to Cai X) and No. 81902367 (to Xu J), the Zhejiang Provincial Natural Science Foundation of China under Grant No. LQ19H160026 (to Xu J), Zhejiang Clinical Research Center of Minimally Invasive Diagnosis and Treatment of Abdominal Diseases Grant No. 2018E50003 (to Cai X), Key Research and Development Project of Zhejiang Province No. 2018C03083 (to Cai X), and the Hepatobiliary and Pancreatic Cancer Research of Hubei Chen Xiaoping Science and Technology Development Foundation under Grant No. CXPJJH11900001–2019308 (to Xu J), No. CXPJJH11900009–03 (to Liang X) and No. CXPJJH11900001–2019209 (to Liang X).
Ethics approval and consent to participate
Ethical approval was obtained from the Sir RunRun Shaw Hospital Research Ethics Committee, and written informed consent was obtained from each patient.
Consent for publication
All authors have agreed to publish this manuscript.
The authors declared that they had no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Xu, J., Wan, Z., Tang, M. et al. N6-methyladenosine-modified CircRNA-SORE sustains sorafenib resistance in hepatocellular carcinoma by regulating β-catenin signaling. Mol Cancer 19, 163 (2020). https://doi.org/10.1186/s12943-020-01281-8
- Sorafenib resistance
- Hepatocellular carcinoma
- Circular RNA