Circular RNA circDLC1 Inhibits MMP1-Mediated Liver Cancer Progression via Interaction with HuR

doi:10.21203/rs.3.rs-74121/v1

Download PDF

Research

Circular RNA circDLC1 Inhibits MMP1-Mediated Liver Cancer Progression via Interaction with HuR

https://doi.org/10.21203/rs.3.rs-74121/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 01 Jan, 2021

Read the published version in Theranostics →

Version 1

posted 10 Sep, 2020

You are reading this latest preprint version

Background: N6-methyladenosine (m6A) modification has been demonstrated to be closely related to cancer progression. KIAA1429, a key component of the m6A methyltransferase complex, has recently been reported to promote hepatocellular carcinoma (HCC) progression by regulating the m6A methylation. However, the involvement of circular RNAs (circRNAs) in KIAA1429-mediated HCC progression is still unknown.

Methods: RNA sequencing (RNA-seq) and methylated RNA immunoprecipitation sequencing (m6A-seq) were utilized to identify KIAA1429-regulated circRNAs. The effects of circDLC1 on proliferation and metastasis of hepatoma cells were examined in vitro and in vivo. RT-qPCR was used to measure the expression of circDLC1 in HCC tissues and hepatoma cells. RNA FISH, RIP assays and biotin-labeled RNA pull-down were used to investigate the downstream effector of circDLC1. The downstream targets of circDLC1 were identified using RNA-seq.

Results: Our data demonstrated that circDLC1 was downregulated in HCC tissues and closely relevant to favorable prognosis. Overexpression of circDLC1 inhibited the proliferation and motility of hepatoma cells in vitro and in vivo, while silencing of circDLC1 played the opposite role. Mechanistic investigations revealed that circDLC1 could bind to RNA-binding protein HuR, which subsequently reduced the interaction between HuR and MMP1 mRNAs, thus inhibited the expression of MMP1, finally contributed to inhibition of HCC progression.

Conclusion: Our work suggests that circDLC1, a downstream target of KIAA1429, is a promising prognostic marker for HCC patients, and the circDLC1-HuR-MMP1 axis may serve as a potential therapeutic target for HCC treatment.

Cancer Biology

Hepatocellular carcinoma

KIAA1429

circular RNA

RNA-binding protein HuR

MMP1

Liver cancer is one of the most common malignant tumors and the fourth leading cause of cancer-related death worldwide (1). Hepatocellular carcinoma (HCC), which is the major form of liver cancer, consists of around 90% of primary liver cancer (2). The high rate of postsurgical recurrence and metastasis is the main cause of poor prognosis in HCC patients, with a five-year survival of less than 20% (3, 4). However, the molecular mechanisms underlying HCC growth and metastasis remain largely unknown.

N6-Methyladenosine (m6A) modification, methylation of adenosine at the N6 position, is a widespread and abundant modification in messenger RNA (mRNA) and non-coding RNAs (ncRNAs) in eukaryotes (5). Emerging studies revealed that m6A modification is tightly correlated with cancer-associated signaling pathways, including RNA stability, cell fate determination, transcription splicing, translation efficiency, and drug resistance (6-8). KIAA1429, a key component of the m6A methyltransferase complex, has been proved to be crucial for RNA m6A methylation (5). A series of studies have shown that KIAA1429 is aberrantly expressed in many types of cancer, such as breast cancer (9), gastric cancer (10) and HCC (11, 12) and regulates the proliferation and metastasis of cancer cells. Previous studies only focused on KIAA1429-mediated post-translational modification on mRNA. However, the modification of KIAA1429 on non-coding RNAs and its effect on downstream signaling pathway remain unknown.

Circular RNAs (circRNAs), known as covalently closed and single-stranded transcripts, have been proved to exist in many species (13). Recent studies have reported that circRNAs was aberrantly expressed in many types of cancer, such as HCC (14), gastric cancer (15), and pancreatic ductal adenocarcinoma (16). The main function of circRNAs in tumor is sponging microRNAs (miRNAs), encoding protein or binding protein (14-16). Interestingly, more recent studies revealed that methyltransferase-like-3 (METTL3), core complex of methyltransferase, is involved in circRNA biogenesis and function by regulating m6A modification (17, 18). Thus, the specific roles of KIAA1429 in mediating circRNAs pathway are urgently needed to explore.

In this study, we identified circDLC1 (derived from the exons 14, 15 and 16 of the DLC1 gene) as a KIAA1429-regulated circRNA, which was closely associated with favorable prognosis. In addition, circDLC1 was found to interact with RNA-binding protein HuR and block the interaction between HuR and MMP1 mRNAs. Our study identified a novel circRNA that was regulated by m6A modification and provided new insights into the pathogenesis of HCC.

Human samples

All HCC and paired adjacent normal tissue samples (n = 110) were obtained from deidentified patients who had undergone surgical resection in West China Hospital (Sichuan University, Chengdu, China). The protocols of this study were conducted in accordance with the Ethical Review Committees of Sichuan University. According to the policies of the committee, the written informed consent was provided from each patient.

Cell lines

Huh-7, Hep3B, HepG2 and SK-Hep1 cell lines were purchased from the Shanghai Cell Bank Type Culture Collection Committee (CBTCCC, Shanghai, China). SNU449 cell line was obtained from ATCC (Manassas, VA, USA). Short tandem repeat (STR) analysis were performed by third-party biology services with all cell lines to characterize the identities of the cells (Feiouer Biology Co., Ltd., Chengdu, China). All cells were cultured in DMEM (HyClone, Logan, UT, USA) containing 10% FBS (HyClone) at 37 °C and 5% CO₂. Actinomycin D was purchased from Sigma-Aldrich.

Cross-Linked RNA immunoprecipitation assays

0.3% formaldehyde was used to cross-link and glycine solution was used to quench the cells (Millipore). Cross-Linked RNA immunoprecipitation assays were performed according to the manufacturer’s protocol of Magna RIP^TM RNA-binding Protein Immunoprecipitation Kit (Millipore, Massachusetts, USA). In brief, magnetic beads were incubated with 5 μg of antibody against HuR (Abcam, Cambridge, USA) or DHX9 (Abcam, Cambridge, USA) and normal IgG (Millipore, Massachusetts, USA) for 30 min at 25 °C. The cell lysate and treated magnetic beads were incubated overnight at 4 °C. After 6 times washing, proteinase K digestion buffer was used to digest the RNA-protein complexes for 45 min at 55 °C. The enriched RNAs were extracted by phenol: chloroform: isoamyl alcohol (125 : 24 : 1). Then, the purified RNAs were subsequently analyzed by RT-qPCR to assess the enrichment of target RNAs to the target proteins.

Biotin-labeled RNA pull-down

0.3% formaldehyde was used to cross-link and glycine solution was used to quench the cells (Millipore). Sangon Biotech synthesized the biotin-labelled circDLC1 probes and control probes (Additional file1: Table S6). Biotin-labeled RNA pull-down was conducted according to the manufacturer’s instructions of the EZ-Magna ChIRP RNA Interactome Kit (Millipore). Briefly, 1/5 associated RNA-protein complexes were separated for RNA purification, and the remaining were subjected to protein purification. RNA was purified with RNeasy Mini Kit (Qiagen Inc.). RT-qPCR was used to identify and quantify the enrichment of circDLC1. On the other hand, 4/5 associated RNA-protein complexes were treated with RNase A, RNase H and DNase to elute proteins. The purified proteins were subsequently analyzed by Western blot.

Animal studies

All animal experiments were authorized by the Institutional Animal Care and Use Committee of Sichuan University. In this study, the BALB/c nude mice were all male and 6 weeks old, which were purchased from BEIJING HFK BIOSCIENCE (Beijing, China). All mice were fed under standard pathogen-free conditions. All surgical procedures were performed with sodium pentobarbital anesthesia. 5 × 10⁵ cells were injected subcutaneously into the right axilla of mice. Tumor volume was measured by a caliper weekly and calculated as length × width² × 0.52. For the liver orthotopic-implanted models, each liver of mice was injected with 1 × 10⁶ cells. After 6 weeks, mice were sacrificed, and the livers were visualized using the [email protected] Lumina II system (Caliper Life Sciences, Hopkinton, MA). All tumor samples were embedded in paraffin and confirmed by hematoxylin and eosin staining. Furthermore, 2 × 10⁶ cells were injected into mice to established the tail intravenous injection models through the tail vein. After 6 weeks, mice were sacrificed, and the lungs were visualized using the [email protected] Lumina II system (Caliper Life Sciences, Hopkinton, MA).

Statistical analysis

All statistical analyses were performed using GraphPad Prism 8 software (GraphPad Software). The Student t test, x2 test, the Wilcoxon signed-rank test, was conducted as appropriate. Kaplan-Meier method was applied to measure the survival curves , and the log-rank test was used to assess the differences. Correlations were assessed by Pearson correlation analysis. P values less than 0.05 indicated statistical significance. *P < 0.05, **P < 0.01, ***P < 0.001.

Supplementary tables are described in Additional file 1. Supplementary methods and figure legends are described in Additional file 2.

CircDLC1 is a KIAA1429-regulated circRNA and clinically relevant to HCC patient prognosis.

To identify circRNAs that are regulated by KIAA1429, we first characterized circular RNA transcripts using RNA-seq from KIAA1429 stable knockdown cells (shKIAA1429) and control cells (shCtrl). We detected 73 distinct circRNAs with 2-fold change (p<0.05). Among these differentially expressed circRNAs, 33 were upregulated and 40 were downregulated in shKIAA1429 cells compared with shCtrl cells (Figure 1A). Our previous work has revealed the target transcripts of KIAA1429-mediated m6A methylation by m6A-seq (12). To identify the circRNAs directly regulated by KIAA1429, we take intersection of RNA-seq and m6A-seq, 6 circRNAs (circDLC1, circCACUL1, circUBAP2, circSETD2, circIGF2BP2, circPICALM) were screened out (Figure 1B). Among them, the expression of circDLC1 showed the most remarkable change upon KIAA1429 silencing (Figure 1C), indicating that circDLC1 could be the main downstream target of KIAA1429. To verify the RNA-seq results, the expression of circDLC1 in shKIAA1429 and shCtrl hepatoma cells was examined by RT-qPCR. As shown in Figure 1D, the expression of circDLC1 was significantly increased in shKIAA1429 compared to shCtrl cells.

Next, we examined the expression of circDLC1 in paired HCC and adjacent normal tissue samples (n = 40) by RT-qPCR. The expression of circDLC1 was much lower in HCC tissues than paired adjacent normal tissues (Figure 1E and Additional file 6: Figure S1A). Furthermore, to investigate the correlation between the expression level of circDLC1 and HCC patient prognosis, the expression levels of circDLC1 in 110 HCC patients (including the previous 40 samples) were measured. According to the relative expression level of circDLC1, the patients were divided into circDLC1 high expression group (n = 63) and circDLC1 low expression group (n = 47). It was revealed that circDLC1 low expression group showed more advanced tumor stage (Additional file 1: Table S1), demonstrating that circDLC1 was associated with HCC progression. In univariate analysis, we found that the status of ascites, serum AFP level, TNM stage and BCLC stage, microvascular invasion, macrovascular invasion and circDLC1 expression were correlated with OS or RFS (Additional file 1: Table S2). In multivariate regression analysis, the status of ascites, microvascular invasion and low circDLC1 expression were shown to be independent risk factors for OS, while serum AFP level and low circDLC1 expression were independent risk factors for RFS (Additional file 1: Table S3). In addition, according to the Kaplan-Meier survival analyses, low expression of circDLC1 was associated with worse OS (P < 0.0001)and RFS (P = 0.0003) (Figure 1F, G).

On the other hand, the expression of KIAA1429 in paired HCC and adjacent normal tissue samples (n = 40) was examined by RT-qPCR and IHC. Notably, the expression of KIAA1429 was obviously upregulated in HCC tissues compared to paired adjacent normal tissues (Figure1H, I). Moreover, we revealed a negative correlation between the expression of KIAA1429 and circDLC1 in HCC tissues (Figure 1J). Collectively, these results indicated that circDLC1 was regulated by KIAA1429 and low circDLC1 expression predicted poor prognosis in HCC patients.

The characteristics of circDLC1

Next, we examined the physical circular structure of circDLC1. CircDLC1 is derived from exons 14, 15 and 16 of the DLC1 gene [CircBase ID: hsa_circ_0135718, termed circDLC1], and the result of Sanger sequencing proved the presence of back-splicing junction (Figure 2A). Also, circDLC1 was observed to resist the digestion by RNase R that specifically degraded linear RNAs but not circRNAs (Figure 2B). Owing to its circular structure, circDLC1 was more stable than mDLC1 with actinomycin D treatment (Figure 2C). Moreover, random hexamer and oligo (dT)18 primers were utilized to further confirm the circularity of circDLC1, circDLC1 was reduced in reverse-transcription efficiency by oligo-dT primers due to the lack of polyA tail (Figure 2D). These results indicated that the formation of circDLC1 was not due to genomic rearrangement. To investigate the localization of circDLC1, we performed subcellular RNA fractionation assays and fluorescence in situ hybridization (FISH), which demonstrated the cytoplasmic enrichment of circDLC1 (Figure 2E, F). Collectively, these findings demonstrated the circularity of circDLC1 and the localization of circDLC1 in hepatoma cells.

DHX9 can regulate the expression of circDLC1

Previous study revealed that the long flanking introns with inverted complementary sequences usually facilitated the biogenesis of circRNAs from their internal exons. (19). These complementary sequences can promote the biogenesis of circRNAs through enhancing back-splicing. By comparing the sequence of flanking introns of the circDLC1 gene, the highly reverse complementary sequences were observed between intron 13 and intron 16 (77% identity over 230 nucleotides, Additional file 6: Figure S1B). They were named as I13RC (reverse complementary sequences in intron 13) and I16RC (reverse complementary sequences in intron 16), respectively. Next, to confirm whether the biogenesis of circDLC1 was determined by I13RC and I16RC, we constructed a series of plasmids of circDLC1, including wild type (a 3754-nt region of the DLC1 gene, spanning from intron 13 to intron 16), mutant #1 (deletion of I13RC), mutant #2 (deletion of I16RC) and mutant #3 (deletion of I13RC and I16RC) (20) (Figure 3A). The result of RT-qPCR indicated that only the wild type plasmid could overexpress circDLC1 after transfection with five types of plasmids, which demonstrated that I13RC and I16RC were essential for the production of circDLC1 (Figure 3B). In addition, northern blot further confirmed this result (Figure 3C).

Previous study revealed that DExH-Box Helicase 9 (DHX9) is a nuclear RNA helicase that inhibits the pairing of complementary sequences of RNA (21). It was reported that DHX9 was involved in the biogenesis of circRNA by binding their flanking inverted complementary sequences (21, 22). To test whether DHX9 regulate the expression of circDLC1, we examined the expression of circDLC1 in hepatoma cell lines after silencing DHX9. In DHX9 knockdown cells, the expression of circDLC1 was obviously increased, while preDLC1 did not show significant changes (Figure 3D, E and Additional file 6: Figure S1C). Moreover, RNA immunoprecipitation assays (RIP) demonstrated that the I13RC and I16RC were significantly enriched in DHX9 immunoprecipitates (Figure 3F, G). Also, we measured the expression of DHX9 in HCC samples (n = 40). DHX9 was remarkably upregulated in HCC tissues (Figure 3H, I) and the expression of circDLC1 and DHX9 showed negative correlation in HCC tissues (Figure 3J). Collectively, our data indicated that DHX9 could bind to the flanking inverted complementary sequences of circDLC1 and inhibit the pairing of these sequences, which subsequently inhibited the production of circDLC1.

CircDLC1 inhibits proliferation and metastasis of hepatoma cells in vitro and in vivo

Six human hepatoma cell lines were chosen to measure the expression of circDLC1 by RT-qPCR. Among them, the lowest expression of circDLC1 was observed in HCC-LM9, moderate expression of circDLC1 was observed in SNU 449, HepG2 and Hep3B, and the highest expression of circDLC1 was observed in Huh7 and SK-Hep1(Additional file 7: Figure S2A). For functional studies, small interfering RNAs (siRNAs) that target the back-splice sequence of circDLC1 was used to knock down the expression of circDLC1 in Huh7 and SNU 449 cells, while no significant change was found in mDLC1 and DLC1 protein levels upon circDLC1 knockdown (Additional file 7: Figure S2B, C). Furthermore, we successfully constructed stable circDLC1-overexpressing cells with HCC-LM9 and SNU 449 cells and stable circDLC1 knockdown cells with SK-Hep1 cells (Additional file 7: Figure S2D-F). It was shown that knockdown of circDLC1 significantly promoted cell proliferation and cell cycle, whereas overexpression of circDLC1 inhibited cell proliferation and cell cycle (Figure 4A-D). Furthermore, silencing of circDLC1 enhanced the migratory and invasive capacities of hepatoma cells, while the motility of hepatoma cells was markedly impaired when circDLC1 was overexpressed (Figure 4E-H).

To further explore the effects of circDLC1 in vivo, circDLC1 knockdown SK-Hep1 cells and circDLC1-overexpressing HCC-LM9 cells were subcutaneously injected into nude mice. We observed that tumor volume was remarkably increased when circDLC1 was knocked down, whereas overexpression of circDLC1 showed an inverse result (Figure 5A-D). Next, the effect of circDLC1 on tumor metastasis in vivo was assessed by establishing liver orthotopic implanted models and tail intravenous injection models. Six weeks later, circDLC1 knockdown group showed higher fluorescence value of GFP than control group both in liver and lung. Haematoxylin eosin (H&E) staining validated the number of metastatic foci was also increased both in liver and lung. However, the liver and lung metastasis model exhibited opposite results when circDLC1 was overexpressed (Figure 5E-H). Collectively, these observations indicated that circDLC1 was a tumor suppressor, which inhibited proliferation and metastasis of hepatoma cells in vitro and in vivo.

CircDLC1 interacts with RNA-binding protein HuR

Many studies have demonstrated that circRNAs play regulatory roles via sponging miRNAs and binding proteins (15, 16, 22, 23). To confirm whether circDLC1 serves as “miRNA sponges”, RNA immunoprecipitation assays (RIP) were performed in SNU 449 and HCC-LM9 cells using a specific antibody against argonaute 2 (AGO2). The result showed that there was no significant difference between IgG group and AGO2 group, which indicated that circDLC1 may not act as “miRNA sponges” (Additional file 8: Figure S3A). Next, we conducted bioinformatics analysis to screen the possible binding protein for circDLC1 by using circinteractome (24), CSCD (25) and RBPDB (26). The intersection of these databases suggested RNA-binding protein HuR is the most possible binding protein for circDLC1 (Figure 6A).

To confirm the interaction between circDLC1 and HuR, we then examined the subcellular location of HuR. Subcellular protein fractionation assay revealed that HuR existed in both cytoplasm and nucleus (Additional file 8: Figure S3B). Given that circDLC1 was located in the cytoplasm, RIP assays were conducted using the anti-HuR antibody to verify whether circDLC1 could bind HuR. In HCC-LM9 and SNU 449 cells, the significant enrichment of circDLC1 was observed in HuR immunoprecipitates compared with IgG pellet (Figure 6B, C). Furthermore, we conducted biotin-labeled RNA pull-down using specific biotin-labeled circDLC1 probe (against the back-splice sequence) and control probe. The data showed that circDLC1 and HuR were prominently enriched in circDLC1 probe group in HCC-LM9 and SNU 449 cells (Figure 6D, E). Also, fluorescence in situ hybridization at a subcellular level confirmed interaction between circDLC1 and HuR (Figure 6F). Then we observed that overexpression of circDLC1 did not affect the HuR level significantly (Figure 6G, H). Taken together, our data demonstrated that circDLC1 may physically bind HuR without affecting HuR expression.

CircDLC1 inhibited the expression of MMP1 by reducing the interaction between HuR and MMP1 mRNAs

To better understand the tumor suppressor roles of circDLC1 on hepatoma cells, we performed RNA-sequencing to analyze the differentially expressed genes which affected by overexpressing circDLC1. Hierarchical clustering showed a total of 347 upregulated genes and 277 downregulated genes in circDLC1-overexpressing SNU 449 cells (Figure 7A). Pathway enrichment analysis and Gene Ontology (GO) analysis indicated that the extracellular matrix related pathways and genes showed the most significant changes (Figure 7B and Additional file 8: Figure S3C). Among the extracellular matrix-related genes, a series of matrix metallopeptidases (MMP1, MMP2, MMP3, MMP10) exhibited remarkable reduction. Then we verified the mRNA level of MMP1, MMP2, MMP3 and MMP10 in circDLC1-overexpressing cells (Figure 7C, D).

HuR is the omnipresent member of the Hu/ELAV (human/embryonic lethal abnormal vision) RBP family and is responsible for the stabilization and (or) translation of many target mRNAs, which typically bear U-rich 3 UTRs (27-29). Importantly, previous reports demonstrated that circRNA may impair the stabilization and translation of mRNAs by competitively binding with RNA-binding protein HuR (30, 31). Thus, we hypothesized that circDLC1 may act as “sponges” for HuR, leading to decreased expressions of downstream target mRNA and protein. To test this hypothesis, we tested the RNA levels of four candidate mRNA (MMP1, MMP2, MMP3, MMP10) after knocking down HuR with siRNA. The result showed that only MMP1 was downregulated (Figure 7E, F and Additional file 8: Figure S3D), indicating MMP1 is the downstream target of circDLC1-HuR axis. Then we verified that protein level of MMP1 was significantly decreased in circDLC1-overexpressing cells (Additional file 8: Figure S3E).

Furthermore, we predicted the binding possibility for HuR-circDLC1 and HuR-MMP1 mRNA by two softwares, catRAPID (32) and PRIdictor (33), which indicated that HuR held a high potential to bind with circDLC1 and 3’UTR of MMP1 mRNA (Additional file 9: Figure, S4A, B). More importantly, the prediction suggested that polypeptide at 126-177 of HuR possessed high binding potential both for circDLC1 and 3’ UTR of MMP1 mRNA (Additional file 9: Figure, S4C, D). Next, reduction of the interaction between HuR and MMP1 mRNA in circDLC1-overexpressing cells was confirmed by RIP assays (Figure 7G, H). In addition, we observed that the expression of MMP1 was significantly increased by circDLC1 knockdown and deletion of HuR markedly repressed the promotion (Figure 7I). Taken together, these observations demonstrated that circDLC1 inhibited expression of MMP1 by reducing the interaction between HuR and MMP1 mRNA.

CircDLC1 inhibits the metastasis in hepatoma cells through HuR-MMP1 axis

A previous study has demonstrated that MMP1 inhibited cell proliferation and invasion in HCC (34). Our data proved that MMP1 possessed the inhibition capacity of invasion and migration in hepatoma cells (Additional file 10: Figure S5A-C). Functionally, as revealed by scratch wound healing and transwell assays, knockdown of circDLC1 could promote the motility of hepatoma cells, while the promotion could be blocked by deletion of MMP1 (Figure 8A, B and Additional file 11: Figure S6A, B). In addition, we assessed the expression level of MMP1 in HCC samples (n = 40) by RT-qPCR and IHC. Notably, MMP1 was obviously upregulated in HCC tissues (Figure 8C and Additional file 11: Figure S6C). We also observed that the expression of MMP1 was negatively correlated with the expression of circDLC1 in HCC tissues (Figure 8D). Taken together, these data indicated that the metastasis of HCC can be inhibited by circDLC1 via HuR-MMP1 axis.

Mounting evidence reveal that m6A modification plays a critical role in different human cancers [6, 35]. As a primary component of the m6A methyltransferase complex, KIAA1429 has recently been demonstrated to promote HCC progression [11, 12]. More importantly, growing evidence suggested that the biogenesis and function of circRNAs could be regulated by m6A modification and m6A regulated-circRNAs may provide potential therapeutic targets [17, 18]. In our present work, we identified circDLC1 as a novel KIAA1429-regulated circRNA, which was clinically relevant to favorable prognosis of HCC patients. Moreover, we clarified the molecular mechanism underlying circDLC1-mediated inhibition of growth and metastasis in hepatoma cells. Our findings will facilitate the understanding of the function of m6A-regulated circRNAs and the assessment of m6A-regulated circRNAs as therapeutic targets for HCC.

Multiple studies have reported that the alternative splicing of RNA was regulated by m6A modification [6]. For example, METTL3 methylated the transited adenosine in codon 273 that resulted in overexpression of R273H mutated p53 protein by enhancing the preferential pre-mRNA splicing [36]. In addition, a previous study revealed that circRNAs were derived from their parental pre-mRNAs by different alternative splicing [37]. In our study, we found that the expression of preDLC1 was increased in KIAA1429 knockdown cells as well as the expression of circDLC1 (data not shown), suggesting that the regulatory role of KIAA1429 in the production of circDLC1 is implicated with the processing of preDLC1. The detailed mechanism of KIAA1429-mediated alternative splicing during the processing of preDLC1 into circDLC1 will be fully illustrated in our future work. On the other hand, increasing studies suggested that the production of circRNAs can be regulated by two or more RNA-binding proteins synergistically [37]. In our study, we only demonstrated that DHX9 was involved in the production of circDLC1 by binding to the flanking inverted complementary sequences, which is similar to the mechanism reported by previous studies [21, 22]. However, whether other RBPs are involved in the biogenesis of circDLC1 or not requires further exploration.

Previous studies demonstrated that HuR is responsible for the stabilization and (or) translation of many target mRNAs, which typically bear U-rich 3’ UTRs [27-29]. For instance, HuR has been reported to promote the progression of many cancers by enhancing the expression of multiple angiogenic factors [29]. Recent studies revealed that circRNAs may impair the stabilization and translation of mRNAs by competitively binding with HuR [30, 31]. Consistently, in the present study, we proved that circDLC1 could reduce the interaction between HuR and MMP1 mRNA, therefore inhibit the expression of MMP1 and HCC progression. Our future work will focus on illustrating the competitive binding sites between HuR and circDLC1, as well as the binding sites between HuR and 3’ UTR of MMP1 mRNA.

MMP1, a member of the zinc-dependent endopeptidase family, has been demonstrated to be closely relevant to migration and invasion in a number of cancers [38]. Interestingly, recent studies revealed that MMP1 could also promote cancer cell proliferation, including HCC [34, 39]. In this study, we confirmed that MMP1 is able to promote both the proliferation and metastasis of hepatoma cells. However, the detailed role of MMP1 in proliferation-related signaling pathway and the underlying mechanism remain to be investigated.

In summary, our study identified circDLC1 as a novel downstream effector of KIAA1429-mediated m6A modification and confirmed the critical role of m6A modification in regulating the biogenesis of circRNAs. Furthermore, we demonstrated that circDLC1 suppresses the growth and metastasis of HCC via the HuR-MMP1 axis. Our findings shed new light on the underlying mechanisms of HCC progression and suggest that circDLC1 may serve as a prognostic biomarker and potential therapeutic target of HCC patients.

m6A: N6-methyladenosine; CircRNAs: Circular RNAs; HCC: Hepatocellular carcinoma; RT-qPCR: Reverse transcription reaction and quantitative real-time PCR; RBPs: RNA binding proteins; HuR: ELAV-like protein 1; RIP: RNA-protein immunoprecipitation; IHC：Immunohistochemistry；MMP1：matrix metallopeptidase1.

Ethics approval and consent to participate

The protocols used in this study were approved by the Ethical Review Committees of Sichuan University, and written informed consent was provided from all the patients. The animal studies were authorized by the Animal Ethic Review Committees of the West China Hospital, China. All animal experiments were strictly implemented in compliance with the NIH Guide for the Care and Use of Laboratory Animals.

Consent for publication

All subjects have written informed consent.

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.

Competing interests

The authors have declared no conflict of interest.

Funding

This work was supported by grants from the National Key Technologies R&D Program (2018YFC1106800), the Natural Science Foundation of China (81972747, 81972204, 81872004, 81800564, 81770615, 81702327, 81700555 and 81672882), the Science and Technology Support Program of Sichuan Province (2019YFQ0001, 2018SZ0115, 2017SZ0003), the Science and Technology Program of Tibet Autonomous Region (XZ201801-GB-02), the 1.3.5 project for disciplines of excellence, West China Hospital, Sichuan University (ZYJC18008), the Natural Science Foundation of Guangdong Province (2019A1515011097), the Innovation Program of Shenzhen (Grant No. JCYJ20180508165208399), the Science and Technology Planning Project of Guangzhou (201904010089).

Authors’ contributions

Conceptualization, Hailing Liu, Yong Zeng, Kefei Yuan and Jiwei Huang; Data curation, Hailing Liu, Tian Lan, Hui Li and Lin Xu; Formal analysis, Hailing Liu, Tian Lan and Lin Xu; Funding acquisition, Yong Zeng, Kefei Yuan and Xiangzheng Chen; Investigation, Xing Chen, Haotian Liao and Jinpeng Du; Methodology, Hui Li, Jinju Wang and Xiangzheng Chen; Reagents Support, Xuefeng Li; Software, Haotian Liao and Yunshi Cai; Supervision, Yong Zeng, Kefei Yuan and Jiwei Huang; Original draft, Hailing Liu and Tian Lan; all the authors participated in approval of final manuscript.

Acknowledgements

We thank Yan Wang, Jinkui Pi and Xiangwei Ren from Core Facility of West China Hospital, and Yang Yang, Bo Su from the Animal Experimental Center of West China Hospital for their technique support. We are most grateful for West China Biobanks, Department of Clinical Research Management, West China Hospital, Sichuan University for their support of human tissue samples. We thank for prof. Yan Li (Zhongnan Hospital of Wuhan University, Wuhan, China) for kindly supplying HCCLM9 cells.

Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394-424.
Llovet JM, Zucman-Rossi J, Pikarsky E, Sangro B, Schwartz M, Sherman M, et al. Hepatocellular carcinoma. Nat Rev Dis Primers. 2016;2:16018.
Hlady RA, Robertson KD. A Three-Pronged Epigenetic Approach to the Treatment of Hepatocellular Carcinoma. Hepatology. 2018;68(4):1226-8.
Finn RS, Zhu AX, Farah W, Almasri J, Zaiem F, Prokop LJ, et al. Therapies for advanced stage hepatocellular carcinoma with macrovascular invasion or metastatic disease: A systematic review and meta-analysis. Hepatology. 2018;67(1):422-35.
Dominissini D, Moshitch-Moshkovitz S, Schwartz S, Salmon-Divon M, Ungar L, Osenberg S, et al. Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature. 2012;485(7397):201-6.
Huang H, Weng H, Chen J. m(6)A Modification in Coding and Non-coding RNAs: Roles and Therapeutic Implications in Cancer. Cancer Cell. 2020;37(3):270-88.
Lin S, Choe J, Du P, Triboulet R, Gregory RI. The m(6)A Methyltransferase METTL3 Promotes Translation in Human Cancer Cells. Mol Cell. 2016;62(3):335-45.
Li A, Chen YS, Ping XL, Yang X, Xiao W, Yang Y, et al. Cytoplasmic m(6)A reader YTHDF3 promotes mRNA translation. Cell Res. 2017;27(3):444-7.
Qian JY, Gao J, Sun X, Cao MD, Shi L, Xia TS, et al. KIAA1429 acts as an oncogenic factor in breast cancer by regulating CDK1 in an N6-methyladenosine-independent manner. Oncogene. 2019;38(33):6123-41.
Miao R, Dai CC, Mei L, Xu J, Sun SW, Xing YL, et al. KIAA1429 regulates cell proliferation by targeting c-Jun messenger RNA directly in gastric cancer. J Cell Physiol. 2020.
Cheng X, Li M, Rao X, Zhang W, Li X, Wang L, et al. KIAA1429 regulates the migration and invasion of hepatocellular carcinoma by altering m6A modification of ID2 mRNA. Onco Targets Ther. 2019;12:3421-8.
Lan T, Li H, Zhang D, Xu L, Liu H, Hao X, et al. KIAA1429 contributes to liver cancer progression through N6-methyladenosine-dependent post-transcriptional modification of GATA3. Mol Cancer. 2019;18(1):186.
Chen LL. The biogenesis and emerging roles of circular RNAs. Nat Rev Mol Cell Biol. 2016;17(4):205-11.
Liang WC, Wong CW, Liang PP, Shi M, Cao Y, Rao ST, et al. Translation of the circular RNA circbeta-catenin promotes liver cancer cell growth through activation of the Wnt pathway. Genome Biol. 2019;20(1):84.
Liu H, Liu Y, Bian Z, Zhang J, Zhang R, Chen X, et al. Circular RNA YAP1 inhibits the proliferation and invasion of gastric cancer cells by regulating the miR-367-5p/p27 (Kip1) axis. Mol Cancer. 2018;17(1):151.
Wong CH, Lou UK, Li Y, Chan SL, Tong JH, To KF, et al. CircFOXK2 Promotes Growth and Metastasis of Pancreatic Ductal Adenocarcinoma by Complexing with RNA-Binding Proteins and Sponging MiR-942. Cancer Res. 2020;80(11):2138-49.
Di Timoteo G, Dattilo D, Centron-Broco A, Colantoni A, Guarnacci M, Rossi F, et al. Modulation of circRNA Metabolism by m(6)A Modification. Cell Rep. 2020;31(6):107641.
Zhang L, Hou C, Chen C, Guo Y, Yuan W, Yin D, et al. The role of N(6)-methyladenosine (m(6)A) modification in the regulation of circRNAs. Mol Cancer. 2020;19(1):105.
Zhang XO, Wang HB, Zhang Y, Lu X, Chen LL, Yang L. Complementary sequence-mediated exon circularization. Cell. 2014;159(1):134-47.
Ashwal-Fluss R, Meyer M, Pamudurti NR, Ivanov A, Bartok O, Hanan M, et al. circRNA biogenesis competes with pre-mRNA splicing. Mol Cell. 2014;56(1):55-66.
Aktas T, Avsar Ilik I, Maticzka D, Bhardwaj V, Pessoa Rodrigues C, Mittler G, et al. DHX9 suppresses RNA processing defects originating from the Alu invasion of the human genome. Nature. 2017;544(7648):115-9.
Yu J, Xu QG, Wang ZG, Yang Y, Zhang L, Ma JZ, et al. Circular RNA cSMARCA5 inhibits growth and metastasis in hepatocellular carcinoma. J Hepatol. 2018;68(6):1214-27.
Zhu YJ, Zheng B, Luo GJ, Ma XK, Lu XY, Lin XM, et al. Circular RNAs negatively regulate cancer stem cells by physically binding FMRP against CCAR1 complex in hepatocellular carcinoma. Theranostics. 2019;9(12):3526-40.
Dudekula DB, Panda AC, Grammatikakis I, De S, Abdelmohsen K, Gorospe M. CircInteractome: A web tool for exploring circular RNAs and their interacting proteins and microRNAs. RNA Biol. 2016;13(1):34-42.
Xia S, Feng J, Chen K, Ma Y, Gong J, Cai F, et al. CSCD: a database for cancer-specific circular RNAs. Nucleic Acids Res. 2018;46(D1):D925-D9.
Cook KB, Kazan H, Zuberi K, Morris Q, Hughes TR. RBPDB: a database of RNA-binding specificities. Nucleic Acids Research. 2010;39(Database):D301-D8.
Lebedeva S, Jens M, Theil K, Schwanhausser B, Selbach M, Landthaler M, et al. Transcriptome-wide analysis of regulatory interactions of the RNA-binding protein HuR. Mol Cell. 2011;43(3):340-52.
Mukherjee N, Corcoran DL, Nusbaum JD, Reid DW, Georgiev S, Hafner M, et al. Integrative regulatory mapping indicates that the RNA-binding protein HuR couples pre-mRNA processing and mRNA stability. Mol Cell. 2011;43(3):327-39.
Nabors LB, Gillespie,G.Y., Harkins,L. and King,P.H. HuR, a RNA stability factor, is expressed in malignant brain tumors and binds to adenine- and uridine-rich elements within the 3' untranslated regions of cytokine and angiogenic factor mRNAs. Cancer Res. 2001(61,):2154–61.
Abdelmohsen K, Panda AC, Munk R, Grammatikakis I, Dudekula DB, De S, et al. Identification of HuR target circular RNAs uncovers suppression of PABPN1 translation by CircPABPN1. RNA Biol. 2017;14(3):361-9.
Liu B, Yang G, Wang X, Liu J, Lu Z, Wang Q, et al. CircBACH1 (hsa_circ_0061395) promotes hepatocellular carcinoma growth by regulating p27 repression via HuR. J Cell Physiol. 2020.
Agostini F, Zanzoni A, Klus P, Marchese D, Cirillo D, Tartaglia GG. catRAPID omics: a web server for large-scale prediction of protein-RNA interactions. Bioinformatics. 2013;29(22):2928-30.
Tuvshinjargal N, Lee W, Park B, Han K. PRIdictor: Protein-RNA Interaction predictor. Biosystems. 2016;139:17-22.
Kim E, Kim D, Lee J-S, Yoe J, Park J, Kim C-J, et al. Capicua suppresses hepatocellular carcinoma progression by controlling the ETV4-MMP1 axis. Hepatology. 2018;67(6):2287-301.

Download PDF

Journal Publication

published 01 Jan, 2021

Read the published version in Theranostics →

Version 1

posted 10 Sep, 2020

You are reading this latest preprint version

Circular RNA circDLC1 Inhibits MMP1-Mediated Liver Cancer Progression via Interaction with HuR

Status:

Journal Publication

Version 1

Abstract

Figures

Background

Materials and Methods

Results

Discussion

Conclusion

Abbreviations

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1