M6A-mediated upregulation of circMDK promotes tumorigenesis and acts as a nanotherapeutic target in hepatocellular carcinoma
Molecular Cancer volume 21, Article number: 109 (2022)
Emerging evidence suggest the critical role of circular RNAs (circRNAs) in disease development especially in various cancers. However, the oncogenic role of circRNAs in hepatocellular carcinoma (HCC) is still largely unknown.
RNA sequencing was performed to identify significantly upregulated circRNAs in paired HCC tissues and non-tumor tissues. CCK-8 assay, colony formation, transwell, and xenograft mouse models were used to investigate the role of circRNAs in HCC proliferation and metastasis. Small interfering RNA (siRNA) was used to silence gene expression. RNA immunoprecipitation, biotin pull-down, RNA pull-down, luciferase reporter assay and western blot were used to explore the underlying molecular mechanisms.
Hsa_circ_0095868, derived from exon 5 of the MDK gene (named circMDK), was identified as a new oncogenic circRNA that was significantly upregulated in HCC. The upregulation of circMDK was associated with the modification of N6-methyladenosine (m6A) and poor survival in HCC patients. Mechanistically, circMDK sponged miR-346 and miR-874-3p to upregulate ATG16L1 (Autophagy Related 16 Like 1), resulting to the activation of PI3K/AKT/mTOR signaling pathway to promote cell proliferation, migration and invasion. Poly (β-amino esters) (PAEs) were synthesized to assist the delivery of circMDK siRNA (PAE-siRNA), which effectively inhibited tumor progression without obvious adverse effects in four liver tumor models including subcutaneous, metastatic, orthotopic and patient-derived xenograft (PDX) models.
CircMDK could serve as a potential tumor biomarker that promotes the progression of HCC via the miR-346/874-3p-ATG16L1 axis. The PAE-based delivery of siRNA improved the stability and efficiency of siRNA targeting circMDK. The PAE-siRNA nanoparticles effectively inhibited HCC proliferation and metastasis in vivo. Our current findings offer a promising nanotherapeutic strategy for the treatment of HCC.
Hepatocellular carcinoma (HCC) is the fifth most common cancer worldwide and the second leading cause of cancer death . Despite rapid advances in diagnosis, surgical techniques, targeted therapies and immunotherapies, the 5-year overall survival rate of HCC patients remains unsatisfactory largely due to recurrent distant metastases and resistance to antitumor drugs. The underlying molecular mechanisms of HCC tumorigenesis, metastasis, and resistance still remain unclear. Therefore, it is essential to identify prognostic biomarkers, explore the potential mechanisms of HCC tumorigenesis and progression, and develop targeted therapies.
Circular RNAs (circRNAs) are another class of endogenous non-coding RNAs that are generated from introns or exons through back-splicing. They possess a covalently closed circular loop without 5’ end caps and 3’ polyadenylated tails . Many studies have demonstrated that circRNAs can act as sponges for microRNAs (miRNAs) or bind to proteins, and altered circRNA levels can lead to aberrant expression of gene products . An increasing number of studies have shown that circRNAs are widely expressed in mammals and are expressed to varying degrees in various cancers, such as esophageal squamous cell carcinoma, urothelial carcinoma, and colorectal cancer [4,5,6]. In addition, some circRNAs have been reported to participate in the development of HCC, regulating cell proliferation, tumor metastasis, and glycolysis [7,8,9]. Recent researches elucidated the role of some circRNA in HCC. For example, circASAP1 promoted HCC growth and pulmonary metastasis by enhancing the expression of MAPK1 and CSF-1 via sponging miR-326 and miR-532-5p , while circDLC1 inhibited HCC development and progression via circDLC1-HuR-MMP1 axis, serving as a potential therapeutic target for HCC . CircRNA-SORE increases the sensitivity of HCC cells to sorafenib by inhibiting the Wnt/β-catenin pathway . However, it remains to be explored whether HCC-associated circRNAs play a tumorigenic role through other regulatory mechanisms.
N6-methyadenosine (m6A) modification is emerging as a new frontier of epigenetic regulation with multiple functions including splicing, export, protein translation, stability and tumorigenesis , which has drawn increasingly more attention. m6A RNA methylation is the most prevalent internal mRNA modification in mammals and is mediated by m6A methyltransferases (METTL3/14, WTAP, RBM15/15B, KIAA1429, VIRMA and ZC3H13, called “writers”), demethylases (FTO, ALKBH5 and ALKBH3, called “erasers”) or m6A binding proteins (YTHDC1/2, YTHDF1/3, IGF2BP1/3, HNRNP and eIF3, called “readers”) . Several members (e.g. METTL3, FTO and IGF2BPs) actively participate in diverse human cancers such as acute lymphoblastic leukemia , breast cancer  and endometrial cancer . A bridge between m6A and circRNAs has been established earlier. For instance, m6A-mediated circRNAs that primarily localized to the cytoplasm resulted in enhanced of mRNA stability or increased of protein translation [17,18,19]. For instance, METTL3-dependent m6A was involved in the DGCR8-mediated maturation of miR126 and pri-miR126 . In HCC cells, KIAA1429 could induce m6A methylation of GATA3, promoting malignant phenotypes . Nevertheless, the roles of m6A-modified circRNAs in HCC still need further investigations.
In recent years, the remarkable progress in genome research has allowed gene therapy to emerge as a new approach to treating cancer. RNA interference (RNAi), a method of regulating target genes, shows promise for the developing novel molecular therapeutic drugs that encode “non-drug” targets which are not suitable for conventional therapies. RNAi molecules can be divided into miRNA, small interfering RNA (siRNA) and short hairpin RNA (shRNA) . Among them, siRNAs show great potential for nucleic acid therapy because of their strong and specific RNAi triggering activity . SiRNAs are short double-stranded RNAs that enable sequence-specific gene silencing of complementary messenger RNAs (mRNAs), induce mRNA degradation and inhibit the production of target proteins. The siRNA-based therapies have emerged as a promising strategy for targeting a variety of diseases . However, the efficiency of gene silencing by naked siRNAs is very low because naked siRNA molecules are rapidly degraded by nucleases in the blood and undergo rapid renal clearance in vivo . In addition, the large size and negative charge of siRNA impede its penetration into cell membranes and prevent its accumulation within cells. Therefore, effective delivery is a critical issue in bringing siRNA to target cells and tissues. To deliver siRNA efficiently, various materials have been developed, including lipids, polymers, dendrimers, polymer micelles, and metal-core nanoparticles (NPs). Interestingly, poly (β-amino esters) (PAEs) offer prospects to circumvent the difficulty in effective delivery of siRNAs. PAEs are a significant class of biodegradable synthetic polymers that can be degraded by intracellular esterase, resulting in significantly improved biocompatibility . They have been investigated as gene and drug delivery vehicles and tissue engineering scaffolds considering their pH sensitivities, high biocompatibilities and structural diversity for flexible synthesis and gene loading.
In this study, we identified circMDK as an oncogenic circRNA in HCC by RNA sequencing (RNA-seq), and validated ATG16L1 as its target gene by sponging miR-346 and miR-874-3p, resulting in activation of PI3K/AKT/mTOR pathway. PAEs-mediated in vivo delivery of circMDK siRNA (PAE-siRNA) resulted in significant attenuation of tumor progression in four mouse models of HCC. Taken together, our study demonstrated that circMDK might serve as a potential biomarker for HCC, and provided a promising nanotherapeutic strategy for the treatment of HCC.
Material and methods
Patients and specimens
A total of 10 paired HCC and matched adjacent noncancerous liver (ANL) tissues were collected from patients who had undergone primary surgical resection in Mengchao Hepatobiliary Hospital of Fujian Medical University (Fujian, China). Another 35 paired HCC and ANL tissues were obtained from surgical resections of HCC patients without preoperative treatment at Guizhou Provincial People’s Hospital (Guizhou, China). All samples were identified by two pathologists independently and the human specimen collection was approved by the Ethics Committee of Mengchao Hepatobiliary Hospital and Guizhou Provincial People’s Hospital. Written informed consent was obtained from all patients according to the policies of the committee. The 10 paired samples were sent for RNA sequencing (RNA-seq) while the 35 paired samples were used for qRT-PCR and western blot analysis. All samples were stored at -80 °C until use. The detailed clinicopathological features are described in Supplementary Tables 1, 2 and 3.
Data were assessed by SPSS 18.0 and expressed as mean ± SD. The statistical significance of differences was evaluated by two-tailed Student’s ttest or two-way ANOVA. Statistical analyses were performed using Prism software (GraphPad Software). The correlations among circMDK, IGF2BP1 and ATG16L1 expression in HCC patients were calculated by Pearson correlation analysis and P value computed using the R language function. P values < 0.05 were considered statistically significant.
CircMDK is upregulated in HCC
To investigate the role of circRNAs in HCC tumorigenesis, we performed RNA-seq analyses of ribosomal RNA-depleted total RNA obtained from 10 paired HCC and ANL tissues. A total of 92,204 distinct circRNAs were detected, of which four circRNAs were upregulated and 38 were downregulated in HCC compared with ANL tissues (Supplementary Fig. 1A, Fig. 1A). Based on differential expression levels of dysregulated circRNAs, we found that hsa_circ_0095868 (termed circMDK in the remainder of the article) was significantly upregulated (fold change ≥ 10 and p < 0.05). To further verify our RNA-seq data, we detected the expression of circMDK in another 35 paired HCC samples by qRT-PCR (Fig. 1B). As shown in Fig. 1C, circMDK was significantly higher in 74.3% (26 of 35) of HCC tissues. The elevated expression of circMDK was significantly correlated with a higher TNM stage in HCC patients (Supplementary Table 3). We then analyzed the relationship between circMDK expression and prognosis of HCC patients in 35 paired HCC samples. Although the overall survival (OS) curve did not reach statistical significance (p > 0.05), the higher expression of circMDK conferred a poorer prognosis as compared to patients with lower circMDK expression. This lack of statistical significance could be due to the small sample size and short follow-up time; the 35 HCC patients were followed for only up to 30 months (Supplementary Fig. 1B). Moreover, the protein expression of MDK was higher in HCC tissues compared to normal liver tissues based on The Human Protein Atlas (THPA, https://www.proteinatlas.org) and its high expression was associated with a poorer 5-year survival probability (Supplementary Fig. 1C, D).
In addition, the expression levels of circMDK were confirmed in the normal liver cell line (LO2) and seven HCC cell lines including HCCLM3, SK-Hep-1, SMMC-7721, PLC/PRF/5, HepG2, Huh7 and Hep3B (Fig. 1D). Among the HCC cell lines, Hep3B and Huh7 cells showed the higher expression. Thus, we selected Hep3B and Huh7 cell lines to investigate the biological function and downstream regulatory mechanism of circMDK.
Characteristics of circMDK in HCC cells
CircMDK is generated from exon 5 of MDK gene with a length of 104 nt. The back-spliced junction site of circMDK was amplified using divergent primers and confirmed by Sanger sequencing (Fig. 1E). PCR analysis showed that circMDK could be amplified by divergent primers in cDNA reverse-transcribed from random hexamers, but not gDNA primers (Fig. 1F). Resistance to digestion with RNase R, a highly processive 3’ to 5’ exoribonuclease that digests linear RNAs, confirmed that circMDK harbored a closed loop structure (Fig. 1G, Supplementary Fig. 1E). Treatment with Actinomycin D showed that the half-life of circMDK transcript within 24 h was more stable in comparison to MDK (Fig. 1H). The results of nuclear and cytoplasmic fractionation revealed that circMDK was predominantly localized in the cytoplasm (Fig. 1I, Supplementary Fig. 1F), which was consistent with the results of RNA FISH assay (Supplementary Fig. 1G). Collectively, these results demonstrated that circMDK was abundant, circular and stable transcript and was significantly up-regulated in HCC.
Knockdown of circMDK impedes cell proliferation, migration and promotes cell apoptosis
To explore the role of circMDK in HCC progression in vitro, three siRNA specifically targeting the back-splicing region of circMDK (siRNA-#1, 2, 3) and overexpressing circMDK vector (OE-circMDK) were constructed to knockdown and overexpress circMDK. While siRNA-#1, 2, 3 all showed more than 50% knockdown efficiency in Huh7 and Hep3B cells, siRNA-#1 was chosen for the succeeding experiments due to its highest inhibitory effect (Supplementary Fig. 1H). Moreover, we found that siRNA-#1 could successfully knockdown circMDK expression but had no effect on mRNA expression of MDK in Hep3B cells (Supplementary Fig. 1I). CCK-8 assays and colony formation assays showed that silencing circMDK significantly reduced the proliferative capabilities of HCC cells while overexpressing circMDK promoted cell proliferation (Fig. 2A-D). Transwell assays showed that migration and invasion abilities of HCC cells were markedly reduced when circMDK was decreased (Fig. 2E). In contrast, increasing the level of circMDK increased the migration and invasion rate of cells (Fig. 2F). In addition, we also performed flow cytometry and found that knockdown of circMDK enhanced the apoptosis ratio of HCC cells while OE-circMDK led to opposite results (Fig. 2G). In summary, these results collectively indicated that silencing circMDK reduced the proliferation, migration and invasion, and promoted apoptosis in vitro.
The m6A modification of circMDK improves the transcriptome stability
Recent studies have suggested that m6A modification in mRNA and circRNAs is extremely widespread, and functionally modulates the eukaryotic transcriptome to influence RNA splicing, export, localization, translation and stability . To explore whether circMDK undergoes m6A modification, we first predicted the m6A sites in circMDK using an online bioinformatics tool m6A Avar (http://m6avar.renlab.org/) and found one RRACU m6A sequence motif at exon 5 site of circMDK (Fig. 3A). Next, we performed MazF to provide nucleotide-resolution quantification of m6A methylation sites. MazF toxin, an ACA-sequence-specific endoribonuclease, is sensitive to m6A sites, representing the first m6A-sensitive RNA cleavage enzyme. The methyl-sensitivity of MazF allowed simple analyses of both m6A demethylase and methyltransferase activity (Fig. 3B). As shown in Fig. 3C, the m6A level was higher in Hep3B and Huh7 cells compared with LO2 cells. We then detected the mRNA expression of Mett3, Mett14, Fto, IGF2BP1, IGF2BP2 and IGF2BP3 in Hep3B cells, and found that the most abundantly expressed was IGF2BP1, suggesting the IGF2BP1 protein was a m6A-binding protein of circMDK (Fig. 3D). Additionally, the expression levels of IGF2BP1 were determined by western blot and qRT-PCR, which indicated that IGF2BP1 was highly expressed in HCC tissues compared with normal tissues (Fig. 3E and F).
To further investigate whether the interaction between circMDK and IGF2BP1 exerts a biological function, RNA pull-down assays were performed (Fig. 3G). As shown in Fig. 3H, we conducted a biotin-labeled RNA pull-down using a specific biotin-labeled circMDK probe (against the back-splice sequence) and a control probe. The results of silver staining showed that circMDK and IGF2BP1 were prominently enriched in the circMDK probe group. Furthermore, the significant enrichment of circMDK was observed in IGF2BP1 immunoprecipitants by using the anti-IGF2BP1 antibody compared with the IgG pellet (Fig. 3I). In addition, we found that knockdown of IGF2BP1 resulted in the decreased expression of circMDK (Fig. 3J). After treatment with actinomycin D to block the newly synthesized circMDK RNA, silencing IGF2BP1 reduced the RNA stability of circMDK (Fig. 3K). The correlation analyses also suggested that the expression of circMDK was positively related to IGF2BP1 (Fig. 3L). In summary, these findings suggested that IGF2BP1 could bind to circMDK in vitro, and m6A modification enhanced the transcriptome stability of circMDK, which may partially account for the significant upregulation of circMDK in HCC.
CircMDK acts as a sponge for miR-346 and miR-874-3p
Considering the cytoplasmic distribution of circMDK, we hypothesized that circMDK might exert its effects by targeting miRNAs. According to the theory of competing endogenous RNA (ceRNA), circRNAs can share the same miRNAs with mRNA. Thus, we predicted the potential targets of circMDK using miRNA target prediction tools including miRDB, miRanda and TargetScan. We found that 17 candidate miRNAs were overlapped between three databases and selected the top 7 miRNAs for further analysis (Supplementary Fig. 2A). To verify our prediction, we designed a biotinylated-circMDK probe and confirmed the pull-down efficiency in HCC cells over-expressing circMDK (Fig. 4A). We found that miR-346 and miR-874-3p were abundantly pulled down by circMDK probe in Huh7 and Hep3B cells (Fig. 4B). To further consolidate the direct binding of miRNAs and circMDK, we conducted RIP assays using an anti-Ago2 antibody (Fig. 4C), which showed that circMDK was preferentially enriched in Ago2-RIPs compared with control IgG-RIPs (Fig. 4D). Furthermore, the Ago2-RIP assay also demonstrated that the enrichment of circMDK was much higher in miR-346 and miR-874-3p mimic groups than in the miR-Ctrl group (Fig. 4E and F). Subsequently, dual-luciferase reporter assays also showed that overexpression of miR-346 or miR-874-3p reduced the luciferase activity of the wildtype circMDK reporter gene (WT), but not the mutant circMDK vector (Mut) (Fig. 4G, Supplementary Fig. 2B). The results of FISH assay also confirmed that circMDK was colocalized with miR-346 and miR-874-3p in the cytoplasm (Supplementary Fig. 2C).
Moreover, biotin-labeled miRNA pulldown assays showed that the expression of circMDK was markedly elevated in HCC cells transfected with biotin-labeled miR-346 and miR-874-3p, compared with the control group (Fig. 4H). Meanwhile, circMDK silencing significantly increased the expression of miR-346 and miR-874-3p (Fig. 4I), whereas circMDK expression was obviously suppressed by overexpression of miR-346 or miR-874-3p (Fig. 4J). These results indicated that circMDK directly sponged miR-346 and miR-874-3p.
In addition, we investigated the involvement of miR-346 and miR-874-3p in the mechanism where circMDK promotes HCC development and progression. HCC cells were co-transfected with circMDK plasmids together with miR-346 or miR-874-3p mimics. We found that overexpression of miR-346 or miR-874-3p could inhibit cell growth and proliferation (Supplementary Fig. 3A, B), suppress migration and invasion (Supplementary Fig. 3C, D), and promote apoptosis (Supplementary Fig. 3E). These results suggested that the oncogenetic effects of circMDK were partially mediated by negative regulation of miR-346 and miR-874-3p.
CircMDK decoys miR-346 and miR-874-3p to upregulate their common target, ATG16L1
To further explore the regulatory role of miR-346 and miR-874-3p with circMDK, we first predicted the common target genes of miR-346 and miR-874-3p using MicroCosm, TargetScan, miRDB and miRanda (Fig. 5A). The intersection of these databases indicated ATG16L1 was the most possible target gene for circMDK (Fig. 5B). We found that the expression of ATG16L1 could be affected by modulating the expression of miR-346 and miR-874-3p (Fig. 5C and D). Luciferase reporter assays showed that overexpression of miR-346 and/or miR-874-3p repressed the luciferase activity in HCC cells transfected with the wildtype ATG16L1 3’-UTR reporter plasmid, whereas no obvious inhibition was observed in cells transfected with the corresponding mutant reporter plasmid (Fig. 5E). The predicted miRNA binding sites between ATG16L1 and miR-346/miR-874-3p were depicted in Supplementary Fig. 2D.
Moreover, biotin-labeled miRNA pulldown assays verified that ATG16L1 was the target gene of miR-346 and miR-874-3p (Fig. 5F). Further investigations showed that ATG16L1 was highly expressed in HCC tissues based on the results of western blot and qRT-PCR (Fig. 5G, H). Kaplan–Meier survival analysis was performed to report 120 months OS rates based on high or low ATG16L1 expression, and indicated that high expression of circMDK was associated with worse OS (Fig. 5I). We also found that there was a positive correlation between circMDK and ATG16L1 expression in HCC tissues (Fig. 5J). Subsequently, we used THPA (https://www.proteinatlas.org) proteomic data to verify the expression of ATG16L1 in human HCC tissues and normal liver tissues and the data showed that ATG16L1 was remarkably upregulated in HCC tissues (Fig. 5K). Further investigations showed that downregulation of circMDK decreased the expression of ATG16L, while its upregulation increased ATG16L expression both the mRNA and protein levels (Fig. 5L, M). Taken together, these findings indicated the existence of a circMDK-miR-346/miR-874-3p-ATG16L1 regulatory axis.
ATG16L1 is responsible for circMDK-mediated cell proliferation, migration, invasion and apoptosis
To explore the ability of circMDK to promote tumor progression in an ATG16L1-dependent manner, Huh7 and Hep3B cells were transfected with si-circMDK or ATG16L1-overexpression vector. In vitro functional experiments demonstrated that the restoration of ATG16L1 expression partially rescued the regulatory effects of circMDK knockdown on proliferation (Fig. 6A and B), migration (Fig. 6C), invasion (Fig. 6D) and apoptosis (Fig. 6E) in HCC cells. ATG16L1 has been reported to promote tumorigenesis in ovarian cancer by activating the PI3K/AKT/mTOR pathway . Therefore, we investigated whether modulating the circMDK-mir-346/miR-874-3p-ATG16L1 axis could activate PI3K/AKT/mTOR pathway. As shown in Fig. 6F, circMDK silencing decreased ATG16L1 protein levels and the phosphorylation levels of Akt (p-Akt) and PI3K (p-PI3K). In contrast, circMDK overexpression increased ATG16L1 protein levels and the levels of p-Akt and p-PI3K. In addition, the rescue experiments showed that ATG16L1 overexpression increased the levels of p-Akt and p-PI3K, whereas ATG16L1 knockdown had the opposite effects (Fig. 6G). Collectively, these findings demonstrated that circMDK was an oncogenic circRNA that promoted the cell growth, migration, invasion and inhibited apoptosis in HCC via the circMDK-mir-346/miR-874-3p-ATG16L1-PI3K/AKT/mTOR signaling pathway.
Synthesis and characterization of PAE-siRNA complex in vitro
Considering the critical role of circMDK upregulation in HCC, we sought to develop an efficient delivery vehicle of siRNA targeting circMDK as a potential therapy for HCC. As shown in the previous results, siRNA delivery via Lipo3000 can effectively knock down the expression of circMDK and inhibit its effects on cell proliferation, migration and invasion ability in vitro. We then sought to further verify the in vivo effect of circMDK knockdown on the proliferation and metastasis of HCC.
We synthesized a type of PAEs for targeted delivery of siRNA for circMDK. The branched PAEs were synthesized through Michael addition reaction of amines and acrylates using glycerol triacrylate as cross-linker (Fig. 7A). The obtained PAEs had a light-yellow appearance in solution (Fig. 7B) with a diameter of 92.57 nm and a positive zeta potential of 21.5 mV (Fig. 7C and D). Due to the protonation of amine groups derived positive charges in acidic environments, the circMDK siRNA adsorbed onto cationic PAEs via electrostatic interaction (Supplementary Fig. 4A). The morphology of the PAEs was further characterized by transmission electron microscope, which showed a well-defined core structure with uniformed size of ~ 100 nm (Fig. 7E). After binding of siRNA, the diameter of PAE-siRNA complex was increased to ~ 180 nm (Fig. 7C) and the zeta potential was reduced to -9.41 mV (Fig. 7D).
Furthermore, we explored the effect of PAE-siRNA complex on the stability of siRNA to serum nuclease by gel electrophoresis. Free siRNA and PAE-siRNA (100 nM siRNA equivalent) were incubated at 37℃ in DMEM medium with 10% FBS at multiple time points. As shown in Supplementary Fig. 4B, PAE-siRNA complex delayed the degradation of siRNA under FBS condition and the intensity of siRNA remained unchanged at least 6 h after incubation, suggesting that PAE-siRNA could protect naked siRNA from serum nuclease degradation. On the contrary, most of free siRNA rapidly degraded within 0.5 h in FBS incubation. This result indicated that PAE-siRNA complex effectively protected the siRNA against serum nuclease degradation.
To further study the cellular uptake of siRNA in vitro, we examined the internalization of free Cy5-siRNA and PAE-Cy5-siRNA (100 nM siRNA equivalent) in HepG2 cells. After incubating with different Cy5-siRNA formulations for 6 h, the PAE-Cy5-siRNA entered the HepG2 cells in much larger quantity than free Cy5-siRNA (Supplementary Fig. 4C), which suggested that PAE-siRNA complex could promote better internalization of siRNA and cellular uptake as compared to free siRNA.
The effects of knockdown of circMDK with PAE-siRNA in vitro
Firstly, the cytotoxicity assay was employed to test the safety of PAE-siRNA complex. With increasing siRNA concentrations (from 0 to 5 µM), the PAE-siRNA complex did not exhibit apparent toxicity on HepG2 cells compared with Lipofectamine 3000 transfection reagent-mediated siRNA delivery (Lipo3000-siRNA group) (Supplementary Fig. 4D). Subsequently, we performed qRT-PCR to assess the knockdown effect of Lipo3000-siRNA, PAE-siRNA and PAE-nsRNA in HepG2 cells. As shown in Supplementary Fig. 4E, the PAE-siRNA (100 nM siRNA equivalent) had a remarkably stronger inhibitory effect on circMDK mRNA expression, which was downregulated by about 60% compared to the untreated control group.
Furthermore, we found that the proliferative capacity of HepG2 cells treated with PAE-siRNA was significantly reduced based on CCK-8 assay (Supplementary Fig. 4F). Similarly, PAE-siRNA complex impeded the colony-formation ability in comparison with the nontreated and PAE-nsRNA group in HepG2 cells (Supplementary Fig. 4G). The results from the transwell chamber system showed that PAE-siRNA complex suppressed the migration and invasion of HepG2 cells (Supplementary Fig. 4H). To sum up, these results suggested that PAE-siRNA complex inhibited cell growth, migration, and invasion in vitro.
In vivo antitumor effects of PAE-siRNA complex
To explore the potential antitumor effects of PAE-siRNA complex in vivo, we constructed four HCC animal models, including subcutaneous model, metastatic model, PDX model, and orthotopic model.
Antitumor effects and targeting potential in subcutaneous tumor model
To investigate the antitumor effects and biodistribution of PAE-siRNA in vivo, we constructed subcutaneous models in Balb/c nude mice. The general experimental procedure is illustrated in Supplementary Fig. 5A. We firstly investigated the in vivo antitumoral efficacy of PAE-siRNA on Balb/c nude mice bearing subcutaneous HepG2 tumor-bearing models. Briefly, HepG2 tumor-bearing mice were randomly divided into 4 groups with peritumoral injection of saline, free siRNA, PAE-siRNA, and PAE-nsRNA (4 µg siRNA per mouse equivalent, 15µL siRNA of 20 µM stock), respectively. Body weights of mice in all groups were not significantly affected by any treatment measures, indicating no significant toxicity (Supplementary Fig. 5B). As shown in Supplementary Fig. 5C, the tumor growth of PAE-siRNA group was significantly inhibited compared with saline control, and free siRNA injection also exhibited slight effects on tumor suppression. The tumor sizes and tumor weight in PAE-siRNA group were the smallest among four groups and the inhibitory efficiency reached about 60%, calculated by tumor weight (Supplementary Fig. 5D, E).
Analysis of the haematoxylin and eosin (H&E) tumor sections suggested that PAE-siRNA caused more necrosis than in other groups (Supplementary Fig. 5F). The tumor proliferation and malignancy were evaluated according to the expression of Ki-67. Visually, tumor cells in the PAE-siRNA treated group exhibited decreased Ki-67 expression (Supplementary Fig. 5G). IF staining of the tumor sections for ATG16L1 (green fluorescence) indicated that the decreased expression of ATG16L1 induced by PAE-siRNA was more obvious than the other groups (Supplementary Fig. 5H). Moreover, the TUNEL staining (green fluorescence, Supplementary Fig. 6A), Bax (red fluorescence, Supplementary Fig. 6B), Bcl-2 (yellow fluorescence, Supplementary Fig. 6C) and Caspase-3 staining (pink fluorescence, Supplementary Fig. 6D) of tumor tissue sections all showed that PAE-siRNA induced more apoptosis in most cells. These data in subcutaneous models demonstrated our hypothesis that circMDK silencing inhibits the progression of subcutaneous liver tumors.
To furtherly explore the biodistribution of PAE-siRNA in vivo, we intravenously injected free Cy5-siRNA and PAE-Cy5-siRNA (40 µg siRNA equivalent, 150µL siRNA of 20 µM stock) into Balb/c nude mice bearing subcutaneous HepG2 tumors respectively. After 6–24 h of injection, the PAE-Cy5-siRNA effectively accumulated into the tumor tissues as indicated by a higher degree of fluorescence imaging than the free siRNA (Supplementary Fig. 6E). After 24 h injection, major organs as well as tumor tissues were eviscerated for ex vivo imaging. We observed that mice treated with PAE-Cy5-siRNA complex exhibited higher level of fluorescence intensity at tumor sites than those treated with free Cy5-siRNA (Supplementary Fig. 6F). In addition, the PAE-siRNA was mainly eliminated by the liver and the kidneys (Supplementary Fig. 6F). It is consistent with previous report that the liver can rapidly uptake, degrade, and eliminate NPs, while the kidneys could excrete nanoparticles through urine [29, 30]. The tumor size and weights were also significantly suppressed by PAE-Cy5-siRNA complex (Supplementary Fig. 6G, H). Taken together, these results suggested that PAE-siRNA could passively accumulate into tumor tissues through the enhanced permeation and retention effect, thereby minimizing non-specific tissue distribution.
Antitumor effects of PAE-siRNA complex in metastatic tumor model
We generated lung metastases models to explore the role of circMDK in regulating the progression of tumor metastasis in vivo (Supplementary Fig. 7A). No apparent change in body weight was observed during the period of treatment (Supplementary Fig. 7B).
After 2 weeks of tail vein injection of HepG2-luc cells, saline, free siRNA, PAE-siRNA and PAE-nsRNA (40 µg siRNA per mouse equivalent, 150µL siRNA of 20 µM stock) were injected into mice via tail vein, respectively. Subsequently, we obtained bioluminescence images (BLI) of mice to evaluate the tumor progression in vivo. As depicted in Supplementary Fig. 7C, the mice were imaged to visualize the lung metastases. We clearly observed that the reduced fluorescence intensity was weakest in PAE-siRNA group. Similarly, the results of BLI of lung and metastatic nodules in the lungs presented that PAE-siRNA complex decreased the number of metastatic lung nodules compared with the control group (Supplementary Fig. 7D, E). H&E staining of lung indicated that decreased tumor area was observed in the PAE-siRNA group than in other groups (Supplementary Fig. 7F). These results suggested that PAE-siRNA complex effectively suppressed tumor metastasis.
Antitumor effects of PAE-siRNA complex in PDX model
As depicted in Fig. 7F, we constructed HCC PDX tumor models in Balb/c nude mice from primary liver tumor of two patients (P1) to subcutaneous tumor of mice (P2 and P3). In order to verify the accuracy of PDX model establishment, we performed H&E staining and immunohistochemical staining of Hep, a liver cancer specific marker. These results showed that PDX tumor model maintained the structural characteristics and a degree of similarity with the primary tumor (Supplementary Fig. 8A, B). The experimental procedure for the PDX model of HCC is depicted in Fig. 7G. There were no remarkable changes of mouse body weight during the whole treatment period, indicating no apparent systemic toxicity for all treatment groups (Supplementary Fig. 8C).
We also tested the biodistribution of PAE-siRNA at various time points after treatment of free Cy5-siRNA and PAE-Cy5-siRNA (40 µg siRNA per mouse equivalent, 150µL siRNA of 20 µM stock). The fluorescence imaging indicated that PAE-Cy5-siRNA exhibited a better targeting effect on tumor sites than free Cy5-siRNA (Fig. 7H). Furthermore, the results of biodistribution were also verified in ex vivo tissues at 24 h (Fig. 7I).
At the end of the intervention, we observed that PAE-siRNA group resulted in a more significant decrease in tumor sizes and weight than those from the saline control, free siRNA and PAE-nsRNA group (Fig. 7J). As shown in Fig. 7K, H, E analysis of the tumor sections showed increased necrotic area in the PAE-siRNA group. The expression of ATG16L1 detected by IF (green fluorescence) was the weakest in PAE-siRNA group, which indicated that the knockdown of circMDK by PAE-siRNA suppressed the tumor growth through decreasing the expression of ATG16L1 in PDX model.
Antitumor effects of PAE-siRNA complex in orthotopic model
To further reveal the potential therapeutic effect of PAE-siRNA for in vivo HCC, a more relevant animal model, the orthotopic model was employed to investigate the efficacy of antitumor. We established orthotopic hepatic tumor models to mimic tumor microenvironment in vivo (Fig. 8A), and evaluated the antitumor effect of siRNA by systemic and various drug administration (40 µg siRNA per mouse equivalent, 150µL siRNA of 20 µM stock, Fig. 8B). No significant change of mouse body weight was observed during the whole experimental period (Fig. 8C), indicating that there was no apparent systemic toxicity caused by different siRNA treatments.
To further investigate the antitumor effects of PAE-siRNA complex in orthotopic model, we performed PET/CT scans using 18F-FDG as a tumor metabolic marker of mice. Cartoon illustration of PET/CT scanning process was shown in Fig. 8D. As depicted in Fig. 8E, we can clearly see that there was relatively lower accumulation of 18F-FDG in the liver region marked by the arrows of mice treated with PAE-siRNA than in other groups. Similarly, we observed that PAE-siRNA complex significantly suppressed the liver growth of mice compared to the saline group (Fig. 8F). H&E staining analysis of the liver sections presented that more necrotic areas in the PAE-siRNA group than in other groups were observed (Fig. 8G). In addition, the IF of ATG16L1 (green fluorescence) was obviously decreased in PAE-siRNA group (Fig. 8G), suggesting that PAE-siRNA complex reduced the expression of ATG16L1.
In vivo toxicity text of PAE-siRNA
To evaluate the toxicity and side effects of PAE-siRNA complex in vivo, the H&E-stained tissue slices of the major organs (heart, liver, spleen, lung and kidney) were analyzed after multiple dosing treatment. As shown in Supplementary Fig. 8D, no noticeable damage or pathological change were observed in the organ slices between the PAE-siRNA group and the control groups. In addition, there were no abnormal changes in ALT, AST, BUN, and CRE (Supplementary Table 5). These results revealed that PAE-siRNA possessed no evident side effects, and can be a safe agent for further in vivo applications.
The differential expression of various circRNAs has been reported in various diseases, especially in cancers. Their stability, abundance, conservation, and spatiotemporal specificity made them a hotspot of biomedical research in recent years . In this study, we found that circMDK expression was significantly increased in HCC tissues and cell lines. The increased expression of circMDK resulted in a poor prognosis, suggesting a critical tumor-promoting effect. In contrast to linear MDK, circMDK showed stable expression in HCC cells and was mainly localized in the cytoplasm, indicating a role in post-transcriptional gene regulation. In vitro functional assays showed that knockdown of circMDK promoted cell apoptosis and hindered proliferation, migration and invasion of HCC. The difference in expression and function of circMDK may be related to the tissue specificity of circRNA, which suggests the potential of circMDK as a promising prognostic biomarker to guide the development of personalized therapies for HCC patients.
The m6A modification is widespread throughout the transcriptome, accounting for about 0.2% to 0.6% of total adenosine of cellular RNA . M6A modification has been shown to participate in the development of various cancers. The fate of m6A-modified RNA depends on the function of the different proteins that identify them, which impact stability, translation efficiency, subcellular localization, and alternative polyadenylation [33,34,35]. For example, lncRNA GAS5-AS was reported to repress tumorigenesis and metastasis of cervical cancer by enhancing GAS5 stability and regulating m6A modifications of GAS5, which was dependent on ALKBH5 and YTHDF2. Furthermore, ALKBH5 and YTHDF2 are involved in the dynamic regulation of the balance between m6A methyltransferase and demethylase . Here, we identified a predicted m6A site in circMDK by a series of experiments and found that IGF2BP1 could bind to circMDK in vitro. Our findings showed that m6A modification of circMDK improved its RNA stability, which might partially account for the upregulation of circMDK in HCC. In addition to m6A modification, there might exist other mechanisms that are involved in the elevation of circMDK, such as DNA methylation, histone modification, and miRNA dysregulation, which deserve further exploration.
To date, more and more miRNAs have been found to take vital roles in the pathogenesis and progression of cancers . In this study, we identified circMDK as a sponge for miR-346 and miR-874-3p. MiR-346 has been reported to regulate a variety of cellular processes, including oncogenesis, inflammatory response and differentiation [38, 39]. In cervical cancer, miR-346 enhances Ago2 expression to regulate the activity of other miRNAs and migration and invasion . It is reported that miR-874-3p expression was abnormally reduced in a variety of cancers, including ovarian cancer , colorectal cancer  and HCC . Here, our findings revealed the significance of the interaction of circMDK and miR-346/miR-874-3p in HCC tumorigenesis.
In general, as a ceRNA, the function of circRNA depends on its target mRNAs. Using online databases, ATG16L1 was predicted as a common potential target of miR-346 and miR-874-3p, which was further confirmed by luciferase reporter and biotin-labeled miRNA pull-down assays. ATG16L1, a member of autophagy family, plays an important role in the occurrence and development of various diseases . ATG16L1 could be targeted by miR-142-3p, a novel miRNA that regulates autophagy in human colonic epithelial cells, suggesting the role of miRNA-mediated regulation of ATG16L1 in intestinal inflammation and Crohn disease . In our present study, we found that the expression of ATG16L1 was increased in HCC cells and tissues, implying that autophagy may promote HCC progression. Furthermore, multiple studies have reported that the classical PI3K/AKT/mTOR pathway plays a significant negative regulatory role in the formation of autophagosomes [46,47,48]. For instance, ATG16L1 was associated with tumor invasion and metastasis through activation of Akt signaling . Here, we demonstrated that circMDK up-regulates the expression of ATG16L1 by adsorption of miR-346 and miR-874-3p to activate PI3K/AKT/mTOR signaling pathway, and ultimately promoted proliferation, migration and invasion of HCC cells, as well as inhibition of apoptosis.
Additionally, we developed an siRNA-based NPs delivery in this study. The PAE-siRNA complex downregulated circMDK expression through siRNA silencing and facilitated entry into cells without additional chemical modifications. Applying biocompatible PAEs for siRNA delivery is beneficial: they can increase siRNA stability in serum, and enhance siRNA uptake by cells and tumor accumulation. This may be due to the smaller size of PAEs and positively charged surface. To our knowledge, this is the first study on the delivery of circRNA with PAEs in HCC. PAE-siRNA complex effectively down-regulated circMDK expression in HCC cells and exhibited adequate anti-tumor effects in four animal models (subcutaneous model, orthotopic model, lung metastasis model and PDX model). Notably, PAEs also showed no detectable toxicity in vitro or in vivo. This suggests the feasibility and safety of using PAEs as carriers to deliver circMDK siRNA for antitumor therapy. However, the detailed regulatory mechanism of PAE-siRNA complex in cell proliferation, migration, and invasion remains to be investigated in future works.
While we have revealed the oncogenic role of circMDK and its molecular mechanisms in HCC both in vitro and in vivo, there are several limitations that have to be discussed. First, while we have demonstrated the binding of circMDK with IGF2BP1 in vitro, our data do not further reveal the effect of inhibition of circMDK adenosine methylation on IGF2BP1 recruitment in this paper. For instance, how inhibiting the adenosine methylation of circMDK using a global inhibitor (e.g., Sinefungin) or by CRISPR technology could affect IGF2BP1 expression requires further investigations. In addition, the interaction between IGF2BP1 and miR-346/-874-3p was also not addressed in this study. Third, while we report here the overexpression of ATG16L1 in HCC cells and tissues, it is yet to be further verified by transmission electron microscopy or other assays whether the formation of autophagic vesicles is an important attribute in HCC tumorigenesis and whether such vesicles are formed following administration of our PAE-siRNA complex.
In conclusion, we demonstrated circMDK as an oncogenic circRNA and played an important role in tumorigenesis and progression of HCC by circMDK-miR-346/miR-874-3p-ATG16L1 axis. Our study is the first to demonstrate a direct and effective approach to suppressing HCC by PAE-siRNA complex targeting circRNAs, which suggests a promising nanotherapeutic option for HCC.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its Additional files.
Adjacent noncancerous liver
Haematoxylin and eosin
poly (β-amino esters)
short hairpin RNA
small interfering RNA
Wei L, et al. The emerging role of microRNAs and long noncoding RNAs in drug resistance of hepatocellular carcinoma. Mol Cancer. 2019;18:147. https://doi.org/10.1186/s12943-019-1086-z.
Huang A, Zheng H, Wu Z, Chen M, Huang Y. Circular RNA-protein interactions: functions, mechanisms, and identification. Theranostics. 2020;10:3503–17. https://doi.org/10.7150/thno.42174.
Zhou WY, et al. Circular RNA: metabolism, functions and interactions with proteins. Mol Cancer. 2020;19:172. https://doi.org/10.1186/s12943-020-01286-3.
Fan L, Cao Q, Liu J, Zhang J, Li B. Circular RNA profiling and its potential for esophageal squamous cell cancer diagnosis and prognosis. Mol Cancer. 2019;18:16. https://doi.org/10.1186/s12943-018-0936-4.
Wei WS, et al. LRPPRC regulates redox homeostasis via the circANKHD1/FOXM1 axis to enhance bladder urothelial carcinoma tumorigenesis. Redox Biol. 2021;48: 102201. https://doi.org/10.1016/j.redox.2021.102201.
Xu H, et al. CircRNA_0000392 promotes colorectal cancer progression through the miR-193a-5p/PIK3R3/AKT axis. J Exp Clin Cancer Res. 2020;39:283. https://doi.org/10.1186/s13046-020-01799-1.
Huang XY, et al. Exosomal circRNA-100338 promotes hepatocellular carcinoma metastasis via enhancing invasiveness and angiogenesis. J Exp Clin Cancer Res. 2020;39:20. https://doi.org/10.1186/s13046-020-1529-9.
Yu J, et al. Circular RNA cSMARCA5 inhibits growth and metastasis in hepatocellular carcinoma. J Hepatol. 2018;68:1214–27. https://doi.org/10.1016/j.jhep.2018.01.012.
Hu ZQ, et al. Circular RNA Sequencing Identifies CircASAP1 as a Key Regulator in Hepatocellular Carcinoma Metastasis. Hepatology. 2020;72:906–22. https://doi.org/10.1002/hep.31068.
Liu H, et al. Circular RNA circDLC1 inhibits MMP1-mediated liver cancer progression via interaction with HuR. Theranostics. 2021;11:1396–411. https://doi.org/10.7150/thno.53227.
Xu J, et al. N(6)-methyladenosine-modified CircRNA-SORE sustains sorafenib resistance in hepatocellular carcinoma by regulating β-catenin signaling. Mol Cancer. 2020;19:163. https://doi.org/10.1186/s12943-020-01281-8.
Lee Y, Choe J, Park OH, Kim YK. Molecular Mechanisms Driving mRNA Degradation by m(6)A Modification. Trends Genet. 2020;36:177–88. https://doi.org/10.1016/j.tig.2019.12.007.
Wang T, Kong S, Tao M, Ju S. The potential role of RNA N6-methyladenosine in Cancer progression. Mol Cancer. 2020;19:88. https://doi.org/10.1186/s12943-020-01204-7.
Vu LP, 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:1369–76. https://doi.org/10.1038/nm.4416.
Cai X, et al. HBXIP-elevated methyltransferase METTL3 promotes the progression of breast cancer via inhibiting tumor suppressor let-7g. Cancer Lett. 2018;415:11–9. https://doi.org/10.1016/j.canlet.2017.11.018.
Liu J, et al. m(6)A mRNA methylation regulates AKT activity to promote the proliferation and tumorigenicity of endometrial cancer. Nat Cell Biol. 2018;20:1074–83. https://doi.org/10.1038/s41556-018-0174-4.
Chen RX, et al. N(6)-methyladenosine modification of circNSUN2 facilitates cytoplasmic export and stabilizes HMGA2 to promote colorectal liver metastasis. Nat Commun. 2019;10:4695. https://doi.org/10.1038/s41467-019-12651-2.
Yang Y, et al. Extensive translation of circular RNAs driven by N(6)-methyladenosine. Cell Res. 2017;27:626–41. https://doi.org/10.1038/cr.2017.31.
Di Timoteo G, et al. Modulation of circRNA Metabolism by m(6)A Modification. Cell Rep. 2020;31: 107641. https://doi.org/10.1016/j.celrep.2020.107641.
Li X, et al. Inhibition of METTL3/m6A/miR126 promotes the migration and invasion of endometrial stromal cells in endometriosis†. Biol Reprod. 2021;105:1221–33. https://doi.org/10.1093/biolre/ioab152.
Lan T, et al. KIAA1429 contributes to liver cancer progression through N6-methyladenosine-dependent post-transcriptional modification of GATA3. Mol Cancer. 2019;18:186. https://doi.org/10.1186/s12943-019-1106-z.
Kotowska-Zimmer A, Pewinska M, Olejniczak M. Artificial miRNAs as therapeutic tools: Challenges and opportunities. Wiley Interdiscip Rev RNA. 2021;12:e1640. https://doi.org/10.1002/wrna.1640
Saw PE, Song EW. siRNA therapeutics: a clinical reality. Sci China Life Sci. 2020;63:485–500. https://doi.org/10.1007/s11427-018-9438-y.
Shah K, et al. Nanostructured Hyaluronic Acid-based Materials for the Delivery of siRNA. Curr Pharm Des. 2018;24:2678–91. https://doi.org/10.2174/1381612824666180807123705.
Mainini F, Eccles MR. Lipid and polymer-based nanoparticle siRNA delivery systems for cancer therapy. Molecules. 2020;25. https://doi.org/10.3390/molecules25112692.
Perni S, Preedy EC, Prokopovich P. Amplify antimicrobial photo dynamic therapy efficacy with poly-beta-amino esters (PBAEs). Sci Rep. 2021;11:7275. https://doi.org/10.1038/s41598-021-86773-3.
Wen J, et al. Zc3h13 Regulates Nuclear RNA m(6)A Methylation and Mouse Embryonic Stem Cell Self-Renewal. Mol Cell. 2018;69:1028-1038.e1026. https://doi.org/10.1016/j.molcel.2018.02.015.
Satyavarapu EM, Nath S, Mandal C. Desialylation of Atg5 by sialidase (Neu2) enhances autophagosome formation to induce anchorage-dependent cell death in ovarian cancer cells. Cell Death Discov. 2021;7:26. https://doi.org/10.1038/s41420-020-00391-y.
Rojas JM, et al. Time-course assessment of the aggregation and metabolization of magnetic nanoparticles. Acta Biomater. 2017;58:181–95. https://doi.org/10.1016/j.actbio.2017.05.047.
Wong S, et al. Fibrillar pharmacology of functionalized nanocellulose. Sci Rep. 2021;11:157. https://doi.org/10.1038/s41598-020-79592-5.
Zhang Z, Yang T, Xiao J. Circular RNAs: Promising Biomarkers for Human Diseases. EBioMedicine. 2018;34:267–74. https://doi.org/10.1016/j.ebiom.2018.07.036.
Zeng Y, et al. Refined RIP-seq protocol for epitranscriptome analysis with low input materials. PLoS Biol. 2018;16: e2006092. https://doi.org/10.1371/journal.pbio.2006092.
Wang X, et al. N(6)-methyladenosine Modulates Messenger RNA Translation Efficiency. Cell. 2015;161:1388–99. https://doi.org/10.1016/j.cell.2015.05.014.
Li J, et al. Dysregulated m6A-Related Regulators Are Associated With Tumor Metastasis and Poor Prognosis in Osteosarcoma. Front Oncol. 2020;10:769. https://doi.org/10.3389/fonc.2020.00769.
Kasowitz SD, et al. Nuclear m6A reader YTHDC1 regulates alternative polyadenylation and splicing during mouse oocyte development. PLoS Genet. 2018;14: e1007412. https://doi.org/10.1371/journal.pgen.1007412.
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:4909–21.
Zhang N, Hu X, Du Y, Du J. The role of miRNAs in colorectal cancer progression and chemoradiotherapy. Biomed Pharmacother. 2021;134: 111099. https://doi.org/10.1016/j.biopha.2020.111099.
Kempinska-Podhorodecka A, et al. Oncomir MicroRNA-346 Is Upregulated in Colons of Patients With Primary Sclerosing Cholangitis. Clin Transl Gastroenterol. 2020;11: e00112. https://doi.org/10.14309/ctg.0000000000000112.
Yang B, Dong R, Zhao H. Inhibition of microRNA-346 inhibits myocardial inflammation and apoptosis after myocardial infarction via targeting NFIB. Eur Rev Med Pharmacol Sci. 2020;24:11752–60. https://doi.org/10.26355/eurrev_202011_23827.
Guo J, Lv J, Liu M, Tang H. miR-346 Up-regulates Argonaute 2 (AGO2) Protein Expression to Augment the Activity of Other MicroRNAs (miRNAs) and Contributes to Cervical Cancer Cell Malignancy. J Biol Chem. 2015;290:30342–50. https://doi.org/10.1074/jbc.M115.691857.
Wang Y, et al. miR-874-3p mitigates cisplatin resistance through modulating NF-κB/inhibitor of apoptosis protein signaling pathway in epithelial ovarian cancer cells. Mol Cell Biochem. 2021. https://doi.org/10.1007/s11010-021-04271-6.
Wang Y, Chen H, Wei X. Circ_0007142 downregulates miR-874-3p-mediated GDPD5 on colorectal cancer cells. Eur J Clin Invest. 2021;51: e13541. https://doi.org/10.1111/eci.13541.
Leong KW, et al. miR-874-3p is down-regulated in hepatocellular carcinoma and negatively regulates PIN1 expression. Oncotarget. 2017;8:11343–55. https://doi.org/10.18632/oncotarget.14526.
Gammoh, N. The multifaceted functions of ATG16L1 in autophagy and related processes. J Cell Sci. 2020;133. https://doi.org/10.1242/jcs.249227.
Zhai Z, et al. Human autophagy gene ATG16L1 is post-transcriptionally regulated by MIR142-3p. Autophagy. 2014;10:468–79. https://doi.org/10.4161/auto.27553.
Chen X, et al. Ethanol Extract of Brucea javanica Seed Inhibit Triple-Negative Breast Cancer by Restraining Autophagy via PI3K/Akt/mTOR Pathway. Front Pharmacol. 2020;11:606. https://doi.org/10.3389/fphar.2020.00606.
Liu B, et al. Scoparone improves hepatic inflammation and autophagy in mice with nonalcoholic steatohepatitis by regulating the ROS/P38/Nrf2 axis and PI3K/AKT/mTOR pathway in macrophages. Biomed Pharmacother. 2020;125: 109895. https://doi.org/10.1016/j.biopha.2020.109895.
Qu L, et al. Glycyrrhizic acid ameliorates LPS-induced acute lung injury by regulating autophagy through the PI3K/AKT/mTOR pathway. Am J Transl Res. 2019;11:2042–55.
Devis-Jauregui L, Eritja N, Davis ML, Matias-Guiu X, Llobet-Navàs D. Autophagy in the physiological endometrium and cancer. Autophagy. 2021;17:1077–95. https://doi.org/10.1080/15548627.2020.1752548.
The authors appreciate the members of Zanxian Xia’s laboratory for their valuable comments and discussions. We would like to thank Kai Yuan, Fang Chen and Yongzhen Meng for their kind guidance on experiments. And we also appreciate Wangxi Hai, Ledu Zhou, Zhanwei Zhu and Cao Guo for their assistance in performing the PET/CT scans of mice and collection of HCC patients.
This work was supported by National Natural Science Foundation of China (U21A20384 and 82072293); Open Project Program of the State Key Laboratory of Proteomics (SKLP‐O201805); Key Research and Development Program of Hunan Province (2020SK2054); Central South University Graduate Research Innovation Project (CX20200380, 1053320192068, 2018zzts387, 1053320212227).
Ethics approval and consent to participate
The research protocol was reviewed and approved by the Human Research Ethics Committee of Mengchao Hepatobiliary Hospital and Guizhou Provincial People’s Hospital. And written informed consent was obtained from each patient included in the study.
Consent for publication
We have obtained patients’ consent for publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary materials and methods. References. Supplementary Table 1. Clinical information of 10 HCC patient samples used for RNA-sequencing. Supplementary Table 2. Clinical information of 35 HCC patients’ samples for validation of circMDK expression. Supplementary Table 3. Associations between the expression levels of circMDK and the clinicopathological characteristics of 35 HCC patients. Supplementary Table 4. Primers, siRNA, shRNA sequences and antibodies. Supplementary Table 5. Blood biochemistry of tumor-bearing mice treatedwith saline, PAE-nsRNA and PAE-siRNA. Abbreviations. Supplementary Figure 1. The correlation characterization of circMDK, MDK expression, and prognosis in HCC, (A, B, C, D, E, F, G). Supplementary Figure 2. The binding site of miR-346 and miR-874-3p with circMDK, (A, B, C, D). Supplementary Figure 3. MiR-346 and miR-874-3p are responsible for circMDK-mediated proliferation, migration, invasion and apoptosis, (A, B, C, D, E). Supplementary Figure 4. Characterization and effects of knockdown of circMDK with PAE-siRNA complex in vitro, (A, B, C, D, E, F, G, H). Supplementary Figure 5. The antitumor effects of PAE-siRNA complex in subcutaneous hepatic tumors, (A, B, C, D, E, F, G, H). Supplementary Figure 6. Biodistribution of PAE-siRNA complex in subcutaneous tumor model, (A, B, C, D, E, F, G, H). Supplementary Figure 7. Antitumor effects of PAE-siRNA complex in metastatic tumor model, (A, B, C, D, E). Supplementary Figure 8. Validation diagrams for PDX model construction and Histological observation of tissue sections, (A, B, C, D).
About this article
Cite this article
Du, A., Li, S., Zhou, Y. et al. M6A-mediated upregulation of circMDK promotes tumorigenesis and acts as a nanotherapeutic target in hepatocellular carcinoma. Mol Cancer 21, 109 (2022). https://doi.org/10.1186/s12943-022-01575-z
- Hepatocellular carcinoma (HCC)
- N6-methyladenosine (m6A)
- Poly (β-amino esters) (PAEs)
- Nanoparticles (NPs)