Intrahepatic cholangiocarcinoma (ICC) is a malignant disease characterized by onset occult, rapid progression, high relapse rate, and high mortality. However, data on how the tumor microenvironment (TME) regulates ICC metastasis at the transcriptomic level remains unclear. This study aimed to explore the mechanisms and interactions between hepatocytes and ICC cells.
We analyzed the interplay between ICC and liver microenvironment through cytokine antibody array analysis. Then we investigated the role of N6-methyladenosine (m6A) modification and the downstream target in vitro, in vivo experiments, and in clinical specimens.
Our study demonstrated that cytokine CCL3, which is secreted by hepatocytes, promotes tumor metastasis by regulating m6A modification via vir-like m6A methyltransferase associated (VIRMA) in ICC cells. Moreover, immunohistochemical analyses showed that VIRMA correlated with poor outcomes in ICC patients. Finally, we confirmed both in vitro and in vivo that CCL3 could activate VIRMA and its critical downstream target SIRT1, which fuels tumor metastasis in ICC.
In conclusion, our results enhanced our understanding of the interaction between hepatocytes and ICC cells, and revealed the molecular mechanism of the CCL3/VIRMA/SIRT1 pathway via m6A-mediated regulation in ICC metastasis. These studies highlight potential targets for the diagnosis, treatment, and prognosis of ICC.
Intrahepatic cholangiocarcinoma (ICC) emanates from the epithelial lining of the intrahepatic biliary tree, and is a malignant disease characterized by onset occult, rapid progression, high relapse rate, and high mortality [1, 2]. As the second most common primary liver cancer, ICC accounts for about 3% of all gastrointestinal malignancies, and recent years have seen a rapid increase in ICC-related morbidity and mortality [3, 4]. Thus, there is an urgent need to explore the molecular mechanisms mediating the development of ICC, which will lay the basis for identifying novel diagnostic biomarkers and validating therapeutic targets.
Tumor metastasis is one of the most important indicators of tumor progression and a leading cause of mortality in cancer patients. Besides, the tumor microenvironment (TME) has been shown to play a key role in driving tumor development and metastasis . ICC is more prone to intrahepatic spreading and metastasis as compared with extrahepatic cholangiocarcinoma. The ICC cells break through the wall of the intrahepatic bile duct and invade the liver parenchyma. The ICC liver microenvironment is mainly composed of liver parenchymal cells, which account for 85%, and non-parenchymal cells . Previous data have demonstrated that liver cells can modulate the TME via the secretion of cytokines, chemokines, and other tumor-regulating factors which are vital for tumor metastasis [7,8,9]. Other studies have shown that various cytokines and chemokines such as CCL2 and CXCL12 in the liver microenvironment mediate many physiological and pathological processes in cholangiocarcinoma [10, 11].
Chemokines, a family of small proteins (8–10 kDa) that interact with specific G-protein–coupled chemokine receptors, are key components of the TME . The chemokines play an important role in processes controlled at epigenetic, transcriptional, and post-transcriptional regulation levels, which affect the occurrence and development of tumor cells . In addition, there has been an association between extracellular signals and N6-methyladenosine (m6A) modification [14, 15]. For example, TGF-β can directly modulate the m6A dynamics which trigger cancer proliferation and metastasis [16, 17].
The m6A methylation occurs on the sixth N of adenylate (A) in RNA, which is the most universal and abundant epigenetic modification in eukaryotic mRNA [18, 19]. The m6A RNA methylation is a dynamic and reversible process mainly mediated by the methyltransferase complex. The process is mediated by m6A demethylases and is related to diverse cellular functions by recognizing and binding the m6A reader proteins . The m6A methyltransferase complex consists of METTL3, METTL14, VIRMA (vir-like m6A methyltransferase associated), and WTAP, while the m6A demethylases comprise FTO and AlkBH5. m6A methyltransferase includes YTH domain-containing family proteins (YTHDFs, comprising YTHDF1/2/3) and insulin-like growth factor 2 binding proteins (IGF2BPs, comprising IGF2BP1/2/3), which can specifically recognize m6A modification sites and recruit related proteins to regulate the stability of mRNA, nuclear output or translation efficiency. Another study suggested that the m6A modification and m6A binding proteins participate in a range of physiological and pathological processes such as DNA damage , embryonic development , immune response , tumorigenesis, or tumor progression [24, 25]. Besides, m6A methylation plays an essential role in the promotion and metastasis of many types of malignancies such as gastric cancer , hepatocellular cancer , ovarian cancer , breast cancer , and colorectal cancer [30, 31], thus is a potential target for cancer therapy. For instance, VIRMA promotes tumor development in breast cancer , while METTL3-mediated m6A modification of HDGF mRNA modulates the progression of gastric cancer . To date, however, data on the biological roles and mechanisms of m6A modification in ICC remain unclear.
In this study, we analyzed the interplay between ICC and liver microenvironment and investigated the role of m6A modification in ICC. Our study demonstrated that CCL3, which is secreted by hepatocytes, promotes tumor metastasis by fueling m6A modification in ICC. Moreover, we showed that VIRMA was a major modulator in the m6A modification, which correlated with poor outcomes in ICC patients and promoted ICC metastasis. In addition, SIRT1 was shown to be a critical downstream target of VIRMA, which fuels tumor metastasis. Taken together, our data enhanced our understanding of the interaction between hepatocytes and ICC cells, and uncovered the molecular mechanism of VIRMA-mediated m6A modification in ICC metastasis.
Materials and methods
Tissue specimens and clinical data
Ethical approval for this study was provided by the Institutional Review Board of Sun Yat-sen Memorial Hospital. A total of 110 paraffin-embedded ICC specimens and corresponding adjacent normal tissues were obtained from Sun Yat-sen Memorial Hospital, Sun Yat-sen University between January 2011 and December 2020.
Immortalized normal liver epithelial cells (THLE3) and human normal bile duct epithelial cell lines (HIBEpiC) were obtained from the American Type Culture Collection (ATCC, Manassas, VA, USA). The cells were cultured as specified by respective manufacturers. Two intrahepatic cholangiocarcinoma cell lines (HuCCT1 and RBE) were acquired from the Cell Bank of the Chinese Academy of Sciences, Shanghai, China. The cell lines were grown at 37 °C under 5% CO2 in a 1640 medium supplemented with 10% fetal bovine serum (FBS, Hyclone, USA). The cells were seeded at a density of 5 × 105 cells in a T25 culture flask and passaged every 4–5 days.
Here, we employed both the transwell migration and invasion assays. The ICC cells were inoculated in the upper transwell chamber (Corning, USA) with an indicated density of 3–5 × 104 cells per well in the 24-well plates, with or without Matrigel (BD Biosciences, USA), and incubated for 24 h. Thereafter, unmigrated and uninvaded ICC cells were removed with cotton swabs. The ICC cells on the bottom surface were then fixed in methanol for 10 min, at room temperature, followed by staining with 0.5% crystal violet. The experiments were carried out three times and in triplicate.
Quantitative real-time PCR
Total RNA was extracted from tissues and cells using TRIzol Reagent (Invitrogen, USA), following the manufacturer’s instructions. Quantitative real-time PCR (qRT-PCR) was performed using PrimeScript RT reagent Kit and SYBR Premix Ex Taq (Takara, Japan), and employed primer sequences as shown in Additional file 1: Table S1. Results were normalized to GAPDH expression and calculated using the 2–ΔΔCt method. The experiments were repeated at least three times.
Full-length VIRMA (Ensembl: ENSG00000164944) cDNA was designed and synthesized by GenePharma (Shanghai, China), and then the cloned fragments were ligated into pcDNA3.1. Lipofectamine 3000 (Invitrogen, USA) was served as the role of transfecting the plasmids into the ICC cells, following the manufacturer’s protocol. Small hairpin RNA (shRNA) of VIRMA (shVIRMA) and negative control (shNC) were designed as shown in Additional file 2: Table S2 and synthesized by GenePharma (Shanghai, China). In order to acquire stable VIRMA-knockdown ICC cell lines, the cells were infected with lentiviral particles and then selected under 8 μg/ml puromycin pressure.
All animal experiments were approved by the Animal Research Committee of Sun Yat-sen University Cancer Center. Female Balb/c nude mice (Balb/c-nu, weighing ~ 15–20 g), were purchased from the Guangdong Medical Laboratory Animal Center (Guangzhou, China). To generate a subcutaneous xenograft tumor model, we injected 5 × 106 transfected HuCCT1 cells, which were suspended in 0.2 ml of PBS, into the flank of nude mice (five mice per group). On the other hand, an orthotopic xenograft tumor model was established by inoculating 3 × 106 transfected HuCCT1 cells in 50 μl of PBS into the liver of nude mice (five mice per group). For the treated groups, HuCCT1 cells were treated with CCL3 (100 ng/ml) for 6 h before inoculation. For subcutaneous injection, the intratumoral multiple-point injection of CCL3 (20 mg/kg) diluted in 50 μl PBS was performed every 5 days. The control groups were treated with PBS. Subcutaneous tumor size was measured twice a week. After in vivo fluorescence imaging, all the study mice were sacrificed with the in vivo imaging system (IVIS) spectrum after four weeks. The mice tumors and organs were dissected, photographed, weighed, and stained.
Additional methods and experimental details are described in Additional file 4.
Adjacent-cancer hepatocytes promote malignant progression of ICC
This study analyzed a total of 280 cholangiocarcinoma patient records with complete follow-up and clinical data from Sun Yat-sen Memorial Hospital of Sun Yat-sen University between 2009 and 2018 (Fig. 1). The data demonstrated that patients with ICC had a worse prognosis compared to those with ECC (Fig. 1a). In addition, compared to the ECC patients, patients with ICC had a higher chance of distant metastasis (ICC vs ECC: 39.2% vs 24.2%) and liver metastases (ICC vs ECC: 28.4% vs 12.9%) (Fig. 1b). Analysis of the SEER database showed the same trend in survival or prognosis of ICC and ECC patients. The SEER database analysis demonstrated that the median survival of the ICC patients was significantly shorter than that of the patients with ECC (Fig. 1a). Besides, ICC patients were significantly more susceptible to distant metastasis (ICC vs ECC: 35.5% vs 26.1%) (Fig. 1c). These findings showed that ICC is more prone to metastasis and has a worse prognosis than ECC. Thus, we hypothesized that these differences in disease development might be due to differential TME. Analysis of clinical data showed that hepatocytes, an important component of the liver microenvironment, played a crucial role in ICC metastasis. As shown in Fig. 1d, our data demonstrated that the ICC cells grew infiltratively along the intrahepatic bile duct, and then into adjacent liver parenchyma, which was in contact with hepatocytes.
To dissect the functional relationship between hepatocytes and ICC development and progression, we co-cultured the ICC cell lines RBE and HuCCT1 with hepatocytes (THLE3). The growth of the RBE and HuCCT1 cells after co-culture with THLE3 was measured using cell counting kit-8 (CCK-8) assays. The data showed that the co-culture with THLE3 significantly promoted the proliferation of RBE and HuCCT1 cells (Fig. 1e). Moreover, the transwell assay revealed that co-culture with THLE3 robustly enhanced the migratory and invasive ability of the ICC cells (Fig. 1f). In addition, we conducted cell cycle progression and cell apoptosis analysis by flow cytometry. As indicated in Fig. 1g, the co-culture with THLE3 resulted in a significant increase in S-phase cells and a reduction in cells at G0/1 phases. The apoptosis assays showed that the percentage of apoptotic cells was significantly higher in the co-culture group (Fig. 1h). Overall, these data indicated that hepatocytes promote the metastasis of ICC.
Hepatocyte-secreted CCL3 promotes migration and invasion of the ICC cells
To assess the key chemotactic molecule in hepatocytes that affect the metastasis of ICC cells, we employed the cytokine antibody assay to examine the level of cytokines secreted by hepatocytes after co-culture with the ICC cells. As shown in Fig. 2a and Additional file 3: Fig. S1, our results showed a high secretion of cytokine CCL3 and CCL23 in hepatocytes after co-culture with RBE, an ICC cell. We then used qRT-PCR and ELISA to verify the prominent hypersecretion of cytokine CCL3 in supernatants collected from the THLE3 cell line after co-culture with the ICC cells (Fig. 2b and d), while a slightly elevated level of CCL23 secretion was observed in Fig. 2c. Next, to evaluate the role of different cytokines in the regulation of migration and invasion ability of the ICC cells, we conducted transwell migration assays and matrigel invasion assays of the RBE cells and HuCCT1 cells stimulated with CCL3 and CCL23. Our findings showed that CCL3 significantly enhanced the migration and invasion ability of HuCCT1 and RBE cells compared with the control group (Fig. 2e–f). On the contrary, there was no significant difference between the control and the treatment group with CCL23, which indicated that CCL3 was a critical factor in ICC infiltration and metastasis.
The expression of CCL3 in the hepatocytes was significantly upregulated after co-culture with the ICC cells, while there was no significant correlation in the CCL3 expression in ICC cells after co-culture with THLE3 (Fig. 2g, h). Moreover, immunohistochemical staining was performed to profile the localization of CCL3 protein in ICC tumor and adjacent liver tissues. The data showed that the CCL3 was mainly expressed in the cytoplasm of hepatocytes and intercellular matrix adjacent to the tumor tissue, with cancer cells and distant hepatocytes staining negative (Fig. 2i). Functional enrichment analysis showed that genes that interacted with CCL3 were enriched for functions related to RNA metabolism (Fig. 2j). These results suggested that CCL3 is a key chemokine in the interplay between hepatocytes and ICC cells, and could influence RNA metabolic processes such as RNA modification, RNA synthesis, RNA cleavage, and RNA degradation in ICC.
CCL3 promotes ICC tumor migration and invasion by regulating m6A modification via VIRMA
Our dot blot assays demonstrated that m6A was among the genes that were shown to interact with CCL3. Global m6A abundance was detected in the mRNA of RBE cells co-cultured with THLE3 and stimulated with CCL3. As shown in Fig. 3a, b, the results indicated that the m6A in total RNA was upregulated in ICC cells after co-culture with THLE3 and stimulation with CCL3. In addition, the mRNA and protein expression profile of VIRMA, one of the major components of the m6A methyltransferase complex, was upregulated in all m6A-related writers and readers (Fig. 3c).
Further analysis revealed that VIRMA expression was significantly upregulated in ICC cells treated with 100 ng/ml CCL3 for 6 h (Fig. 3d, e). It was shown that CCL3 could not significantly increase VIRMA expression in the ICC cells treated with CCR5 receptor antagonist (Fig. 3f). These results suggested that the binding of CCL3 to the CCR5 receptor on the membrane of cancer cells promotes ICC migration and invasion by regulating VIRMA, an m6A methyltransferase. We then conducted transwell migration and invasion assays of the RBE and HuCCT1 cells with different treatments to assess the role of CCL3/VIRMA pathway on the migratory and invasive ability of the cells. As demonstrated in Fig. 3g, human recombinant cytokine CCL3 effectively promoted migration and invasion of ICC cells, a phenomenon that was inhibited by knockdown of VIRMA, thus confirming the pro-metastatic effects of CCL3/VIRMA pathway in ICC.
VIRMA promotes ICC development and progression in vitro
To examine the roles of VIRMA in the CCA cells, we profiled the VIRMA expression in CCA cell lines using qRT-PCR and Western blot assays. The results illustrated significant upregulation of the VIRMA mRNA and protein levels in ICC cell lines (RBE and HuCCT1) compared with human intrahepatic biliary epithelial cells (HIBEpiC) (Fig. 4a). Thereafter, we performed VIRMA overexpression and knockdown studies in RBE and HuCCT1 cells. The ICC cells were transfected with control siRNA (siNC) or siRNA targeting VIRMA (siVIRMA-1 and siVIRMA-2). On the other hand, the overexpression cell lines were constructed as oeVIRMA while the matched control cell lines were named as vector. The VIRMA expression in the cells was confirmed with qRT-PCR and Western blot (Additional file 3: Fig. S2).
The growth of the RBE and HuCCT1 cells after knockdown or overexpression of the VIRMA was measured using cell counting kit-8 (CCK-8) and colony formation assays. As shown in Fig. 4b–d, there was a declining trend in cell proliferation in RBE and HuCCT1 cells with VIRMA knockdown while overexpression resulted in significantly increased cell proliferation. Moreover, the transwell migration assay demonstrated that downregulation of VIRMA impeded the migratory ability of the ICC cells, while upregulation of VIRMA significantly elevated the migration of the ICC cells. Similarly, the transwell invasion assay showed that the invasive ability of HuCCT1 and RBE cells was significantly suppressed in response to the downregulation of VIRMA, while it was enhanced by VIRMA upregulation (Fig. 4e, f). We then conducted wound healing assays to display the slow migration of cells with VIRMA knockdown, which was in contrast to the increase observed in VIRMA overexpression cells (Fig. 4g, h). In the meantime, we detected the expression of the Epithelial-Mesenchymal Transition (EMT) marker, angiogenic marker vascular endothelial growth factor (VEGF), and proliferation marker Cyclin D1 in ICC cells with VIRMA knockdown or overexpression. As displayed in Fig. 4i, j, the mRNA and protein levels of N-cadherin (mesenchymal marker), Vimentin (mesenchymal marker), VEGF, and Cyclin D1 were significantly positively correlated with VIRMA expression, while there was a prominent inverse correlation between the E-cadherin (epithelial marker) expression and VIRMA expression. These results implied that VIRMA promotes ICC proliferation, invasion, and metastasis in vitro.
Furthermore, we performed cell cycle progression and cell apoptosis analyses. As demonstrated in Fig. 5a, b, VIRMA-deficient cells had a partial cell cycle arrest at the G0/1 transition but a significant increase in the G0/1 phase and a reduction in cells at S and G2-M phases. The apoptosis examination detected by flow cytometry (Fig. 5c, d) and TUNEL assay (Fig. 5e, f) indicated that the percentage of apoptotic cells was significantly elevated in cells transfected with siVIRMA. There appeared an obvious decrease in Bcl-2 expression, the antiapoptotic protein, in ICC cells with down-expression of VIRMA (Fig. 5g, h). From this, we can deduce that downregulation of VIRMA induces ICC cell cycle arrest and promotes cell apoptosis.
VIRMA promotes ICC proliferation and metastasis in vivo and is correlated with poor prognosis of ICC patients
We demonstrated that the downregulation of VIRMA inhibits the proliferation, invasion, and metastasis of ICC cells in vitro. Here, we analyzed the effect of VIRMA in the development of ICC in mice. Stable cells with modified VIRMA expression were subcutaneously injected into the BABL/c nude mice. First, we investigated the effect of VIRMA on the proliferation of ICC. As shown in Fig. 6a–d, VIRMA knockdown (HuCCT1-shVIRMA) led to significantly low tumor growth rate and tumor weight compared with the other groups (HuCCT1: 0.75 ± 0.05 g, HuCCT1-shNC: 0.77 ± 0.02 g, HuCCT1-shVIRMA: 0.28 ± 0.01 g). The H&E staining demonstrated that VIRMA knockdown significantly reduced the tumor density and differentiation degree of the ICC. On the other hand, immunohistochemical staining showed that the expression of tumor progression markers such as N-cadherin, Vimentin, VEGF, Ki67, and Cyclin D1 was decreased along with the downregulation of VIRMA, except E-cadherin exhibited a contrary tendency (Fig. 6j).
In parallel, we analyzed the expression of VIRMA in ICC patients using the TCGA dataset and showed that the VIRMA expression was remarkably increased in CCA tissues (Fig. 6g). Moreover, immunohistochemical staining in 110 human ICC patients and normal adjacent tissues revealed that VIRMA expression was significantly higher in ICC cancer compared to the adjacent counterparts (Fig. 6e, f). In addition, we conducted a Kaplan–Meier survival analysis to investigate the correlation between VIRMA expression and ICC patient prognosis (Fig. 6h, i). The results demonstrated that ICC patients with high VIRMA expression had poorer overall survival (OS) and disease-free survival (DFS) (*P < 0.05). Taken together, these findings suggested that the higher expression of VIRMA was associated with adverse prognosis in ICC patients.
SIRT1 was selected as a downstream target of VIRMA-mediated m6A modification
To identify the target genes of VIRMA, MeRIP-seq and RNA-seq were performed in stable VIRMA downregulation ICC cells. According to the results of the MeRIP-seq, it was found that more than 80% of m6A binding sites were located in the protein-coding region (CDS region) and were enriched in the 5′ untranslated region (5′ UTR) as well as 3′ untranslated region (3′ UTR). Further analysis showed that the proportion of specific m6A binding sites in the 3′UTR was reduced in the siVIRMA group as compared with that in the siNC group and this was consistent with findings that VIRMA was mainly involved in regulating m6A modification in the 3′UTR region of mRNA (Fig. 7a).
The results of the Motif for the m6A binding site indicated that the most significant m6A Motif was GATG (Additional file 3: Fig. S3c), which provided a basis for subsequent screening of the downstream targets. Network enrichment analysis was then conducted according to MeRIP-seq results and it was identified that the cell cycle pathway, as well as the Foxo signaling pathway, were significantly associated with the oncogenic activities of VIRMA (Fig. 7b). Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) enrichment analysis were carried out using gene set enrichment analysis (GSEA) (Additional file 3: Fig. S3). Results of the RNA-seq showed that 1188 transcripts were markedly downregulated on VIRMA knockdown. Further, MeRIP-seq presented that m6A peaks of 1011 transcripts exhibited decreased abundance. In addition, it was found that 2072 genes were significantly correlated with VIRMA in TCGA dataset (Fig. 7c).
Intriguingly, results of the transcriptome-sequencing, MeRIP-seq in the genes selected in the TCGA dataset revealed that 16 transcripts were overlapped (Fig. 7d). Further, it was noted that SIRT1 showed the most pronounced difference and was reported to be close to tumor progression among the candidate genes. The m6A usually happens in RRACH (R = G or A, H = A, C or U) consensus sequence, results of MeRIP-seq in the present study showed that an m6A peak was detected around the stop codon of SIRT1 mRNA in non-target control siRNA cells and were all decreased upon VIRMA knockdown (Fig. 7g). The SIRT1 mRNA was also significantly downregulated in VIRMA-knockdown ICC cells both at mRNA and protein levels (Fig. 7e, f). Furthermore, the RNA pull-down and RIP assays were also conducted to verify the positive association between SIRT1 and VIRMA (Fig. 7h, i). The correlation analysis between SIRT1 and VIRMA based on the TCGA ICC dataset in Fig. 7j reached a similar conclusion. Therefore, a preliminary conclusion was drawn that SIRT1 was selected for further studies as a candidate target of VIRMA-mediated m6A modification.
Downregulation of SIRT1 repressed tumor progression of ICC cells in vitro
As a kind of histone deacetylase, SIRT1 may regulate the chromatin acetylation process. Therefore, we observed that the H3K4 acetylation (H3K4ac) protein level was decreased after SIRT1 knockdown in both RBE and HuCCT1 cells (Fig. 8a). Cell viability was determined by CCK8 assays, from which we observed that downregulation of SIRT1 notably inhibited tumor cell growth in ICC (Fig. 8b). Moreover, we conducted cell cycle progression analysis and cell apoptosis assay of RBE and HuCCT1 cells after SIRT1 knockdown. As illustrated in Fig. 8c, d, depletion of SIRT1 induced tumor cell arrest in the G0/1 phase. And it also markedly increased cell apoptosis in ICC (Fig. 8e, f). Transwell assays (Fig. 8g, h) and wound healing assays (Fig. 8i, j) were performed through SIRT1 knockdown to evaluate the role of SIRT1 in the regulation of migration and invasion ability in the ICC cells. Results of the current study showed that the downregulation of SIRT1 significantly decreased the migration and invasion ability of HuCCT1 and RBE cells compared with the control group.
CCL3/VIRMA/SIRT1 pathway accelerates the ICC malignant process
It was evident that the CCL3 can abolish the effect of VIRMA knockdown in SIRT1 expression (Fig. 9a, b and Additional file 3: Fig. S4) during the rescue experiments in vitro. Therefore, SIRT1 is a downstream target of CCL3/VIRMA in ICC which served as an oncogene to promote migration and invasion of ICC cells.
To explore the effect of the CCL3/VIRMA/SIRT1 pathway on the ICC malignant process, we constructed subcutaneous and orthotopic tumor models in nude mice.
We first injected different cells (HuCCT1 + PBS, HuCCT1 + CCL3, HuCCT1-shNC + CCL3, and HuCCT1-shVIRMA + CCL3) subcutaneously into nude mice. As shown in Fig. 9c–f, tumor growth of HuCCT1 ICC cells was significantly faster than that in the PBS group after CCL3 treatment, while HuCCT1 cells with VIRMA knockdown grew slower even after receiving CCL3 stimulation. Our data showed that the tumor weight in the CCL3 treatment group was significantly higher than that in the PBS group, while the VIRMA knockdown group was not affected by CCL3 treatment, and the tumor weight was significantly lower than that in the control group (HuCCT1 + PBS: 0.81 ± 0.02 g, HuCCT1 + CCL3: 1.45 ± 0.08 g, HuCCT1-shNC + CCL3: 1.53 ± 0.04 g, and HuCCT1-shVIRMA + CCL3: 0.27 ± 0.02 g). The immunohistochemistry analyses confirmed that the expression of VIRMA, SIRT1, and Ki67 was upregulated in the tumor tissues with CCL3 stimulation, and it was also reversed in the tumor tissues with VIRMA down-expression (Fig. 9k). Together, these data demonstrated that CCL3 promotes the growth of ICC cells by regulating VIRMA/SIRT1 in ICC subcutaneous tumor mice model.
We further showed that CCL3 promoted the metastasis of intrahepatic bile duct cancer cells by regulating VIRMA/SIRT1 pathway in the orthotopic ICC tumor model. The liver metastasis was significantly enhanced with CCL3 stimulation, whereas the progression was significantly inhibited by VIRMA knockdown (Fig. 9g–j). Furthermore, IHC staining of the nude mice tissues (Fig. 9l) also verified that SIRT1 was positively correlated with VIRMA expression. With CCL3 stimulation, the Ki67 was up-regulated in the tumor tissues and it was also reversed in the tumor tissues with VIRMA down-expression. Consistent with the above observations, it was evident that the CCL3 had positive effects on tumor growth and liver metastasis of ICC and the progression could be abolished through VIRMA knockdown.
In conclusion, the results of the current study confirmed that CCL3/VIRMA/SIRT1 pathway promotes ICC proliferation and tumor metastasis in vivo.
Metastasis is the main clinical challenge in ICC. Tumor metastasis is affected by the TME, a unique environment that arises because of tumor progression. Most studies on the mechanism of ICC metastasis have focused on the analysis of the ICC cells, but only a few have assessed the role of other cells in TME. The liver cancer microenvironment is composed of cells such as immune cells, fibroblasts, and liver parenchymal cells as well as extracellular matrix components . Previous studies have demonstrated that ICC development is affected by the liver microenvironment, which promotes tumor proliferation, infiltration, and metastasis by releasing cytokines and regulating the transcription and translation of tumor-related genes [34, 35].
In this study, we showed that the co-culture of the ICC cells with hepatocytes significantly enhanced the migration and invasion ability of the cells. Moreover, to clarify the molecular interaction between hepatocytes and ICC cells, we performed a cytokine microarray analysis and showed that CCL3 is the main chemokine that is involved in the interaction. CCL3, also referred to as macrophage inflammatory protein-1 (MIP-1α), belongs to the C–C chemokine family. Studies have shown that CCL3 contributes to the development and progression of various malignancies, such as esophageal cancer, breast cancer, osteosarcoma, leukemia, and multiple myeloma . As expected, our data demonstrated that recombinant human CCL3 obviously promoted migration and invasion of the ICC cells.
Thereafter, we analyzed the mechanisms of CCL3 activity in the communication between hepatocytes and ICC. Functional enrichment analysis showed that the genes that interacted with CCL3 were enriched for functions related to RNA metabolism. We further demonstrated that chemical modification occurs at the RNA level, which affects gene-regulatory networks at the epigenetic level. Previous studies demonstrated an association between extracellular signals such as cytokines or chemokines and m6A modification [14, 15]. The m6A RNA modification has been identified as a source of epigenetic regulation in DNA modification. The m6A modification accounts for more than 80% of RNA modification and is a reversible and dynamic process regulated by m6A methyltransferase, demethylases, and readers . m6A methylation of mRNA has been reported to play significant roles in many RNA functions, such as RNA splicing, RNA processing, translation regulation, and RNA decay [18,19,20, 38, 39]. However, data on the role of m6A modification in ICC remains limited. In this study, there was an increased m6A modification in ICC cells co-cultured with hepatocytes, which was potentially affecting cancer metastasis. We showed that CCL3 was the key factor mediating the interaction between the liver microenvironment and ICC cells. Subsequently, we characterized methyltransferase VIRMA and showed that it is the aberrant m6A modification in ICC cells co-cultured with hepatocytes.
VIRMA is an important member of the m6A methyltransferase complex and is mainly located in the nucleosome and cytoplasm. As a key regulator of m6A methylation, VIRMA was shown to form a complex with WTAP/METTL3/METTL14 heterodimer to catalyze m6A formation [25, 40]. Recent studies have shown that VIRMA is upregulated in several malignancies such as breast cancer , lung cancer , germ cell tumor , and liver cancer , and is involved in tumor proliferation and metastasis. So far, however, there are limited studies that have evaluated the biological function and mechanism of VIRMA-mediated m6A modification in ICC. In this study, we showed that VIRMA was significantly upregulated in ICC tissues, which correlated with poor prognosis and thus regarded as a prognostic factor for tumor recurrence in ICC. Moreover, our in vitro and in vivo assays demonstrated that VIRMA promotes malignant development and progression of ICC, and its expression showcased a significant correlation with various metastatic, angiogenic and proliferation gene markers. These findings suggested that VIRMA may function as a tumor-promoting gene and might be a potential clinical prognostic biomarker and therapeutic target for ICC.
Analysis of the RNA-Seq and MeRIP-Seq data identified SIRT1 (Sirtuin type 1), as a crucial downstream target of VIRMA. Our RNA pull-down and RIP assays demonstrated that SIRT1 was positively regulated by VIRMA and modified by VIRMA-mediated m6A methylation. Thereafter, we conducted GO and KEGG enrichment analyses to explore the major biological functions and downstream signaling pathways of the genes from the RNA-Seq and MeRIP-Seq analyses, respectively. SIRTI, a nicotinamide adenine dinucleotide (NAD+)-dependent histone deacetylase, plays important roles in oncogenesis via transcription, translation, and post-translational modifications . SIRT1 was shown to be a potential tumor promoter based on its role in negatively regulating many tumors suppressors [45, 46]. In sync with previous studies, our findings showed that SIRT1 facilitates the proliferation and migration of tumor cells in ICC. As previously demonstrated, SIRT1 modulates various pathways such as the p53 signaling pathway [47, 48] and FoxO signaling pathway  to control cell cycle, metabolism, proliferation, differentiation, and epigenetics. These data agree with that obtained by GO and KEGG enrichment analysis of MeRIP-seq. Therefore, CCL3/VIRMA/SIRT1 axis may be one of the important mechanisms underlying tumor evolution and metastasis and is a potential new therapeutic target for ICC.
In our study, we only employed the m6A dot blot assay to quantify the m6A levels of ICC, and showed that SIRT1 is a downstream target of VIRMA. Thus, there is a need for further studies which would take the following variables into account: (1) examine the m6A levels of ICC using a colorimetric strategy or liquid chromatography-mass spectrometry (LC–MS), (2) detect whether VIRMA functions independently of its m6A catalytic activity in cancer progression, and (3) develop a peptide inhibitor to target VIRMA domain and explore whether it may be beneficial in the treatment of ICC.
In conclusion, our study highlights the critical role of VIRMA-mediated m6A modification in ICC progression and metastasis. Mechanistically, we demonstrated that CCL3 is secreted by hepatocytes and may promote metastasis of ICC cells by regulating m6A methylation. The regulation of m6A methylation is mediated by VIRMA, which epigenetically promotes SIRT1 expression through an m6A methylation-dependent mechanism. Our results suggested that the interaction between hepatocytes and ICC cells might offer a possible interventional target for ICC. Besides, the m6A modification on tumor metastasis will contribute to further studies that would explore molecular mechanisms and identify efficient treatment strategies against ICC.
Availability of data and materials
The dataset supporting the conclusions of this article is included within the article and its additional files.
Razumilava N, Gores GJ. Cholangiocarcinoma. Lancet. 2014;383:2168–79.
Chen H, et al. RNA N6-methyladenosine methyltransferase METTL3 facilitates colorectal cancer by activating the m6A-GLUT1-mTORC1 axis and is a therapeutic target. Gastroenterology. 2021;160:1284-1300.e16.
Gentilini A, Pastore M, Marra F, Raggi C. The role of stroma in cholangiocarcinoma: the intriguing interplay between fibroblastic component, immune cell subsets and tumor epithelium. Int J Mol Sci. 2018;19:2885.
Korbecki J, Grochans S, Gutowska I, Barczak K, Baranowska-Bosiacka I. CC chemokines in a tumor: a review of pro-cancer and anti-cancer properties of receptors CCR5, CCR6, CCR7, CCR8, CCR9, and CCR10 Ligands. Int J Mol Sci. 2020;21:7619.
Miranda-Gonçalves V, et al. The component of the m6A writer complex VIRMA is implicated in aggressive tumor phenotype, DNA damage response and cisplatin resistance in germ cell tumors. J Exp Clin Cancer Res. 2021;40:268.
The authors would like to thank the patients who participated in our study.
This work was supported by the National Natural Science Foundation of China (Grant Nos. 82002505, 82103142, and 81874057), and the Medical Science and Technology Foundation of Guangdong Province (Grant No. A2022371).
Present address: Medical Research Center, Sun Yat-Sen Memorial Hospital, Sun Yat-Sen University, Guangzhou, 510120, China
Present address: Department of Pancreato-Biliary Surgery, Sun Yat-Sen Memorial Hospital, Sun Yat-Sen University, Guangzhou, 510120, China
Shurui Zhou, Kege Yang, Shaojie Chen and Guoda Lian contributed equally to this work
Authors and Affiliations
Guangdong Provincial Key Laboratory of Malignant Tumor Epigenetics and Gene Regulation, Sun Yat-Sen Memorial Hospital, Sun Yat-Sen University, Guangzhou, 510120, China
YL made substantial contributions to the conception and design of the work. SZ and KY were the major contributors to experimental validation and manuscript writing. SC and GL revised the manuscript. YH, HY, YZ, and KH analyzed and interpreted the data. DY and HL were responsible for project administration. All authors read and approved the final manuscript.
The study was conducted according to the guidelines of the Declaration of Helsinki. The collection of the clinical sample was approved by the Institutional Review Board of Sun Yat-sen Memorial Hospital, Sun Yat-sen University under the Grant Number: SYSKY-2022-492-01. The in vivo study was approved by the Animal Experimentation Ethics Committee of Sun Yat-sen University Cancer Center, Sun Yat-sen University under Grant Number: L102012020080F. Ethical review and approval were not applicable in bioinformatics and in vitro analysis performed in this study.
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.
Functional enrichment analysis of cytokine antibody array. Figure S2. Validation of VIRMA knockdown and overexpression in ICC cells. (a) qRT-PCR (up) and western blot (down) analysis of VIRMA expression in VIRMA-knockdown ICC cells. siVIRMA-1 and siVIRMA-2 indicated VIRMA knockdown in both RBE and HuCCT1 cells. (b) qRT-PCR (up) and western blot (down) analysis of VIRMA expression in VIRMA-overexpression ICC cells. oeVIRMA indicated RBE cells and HuCCT1 cells transfected with an overexpressing plasmid. The relative quantification was calculated using the 2−ΔΔCt method and normalized based on GAPDH. Figure S3. MeRIP-seq and RNA-seq analysis. (a) GO analysis of down-regulated genes identified for MeRIP-seq. (b) Motif analysis. (c) The peak annotation for m6A signal enrichment in MeRIP-seq. (d) GO analysis of down-regulated genes identified for RNA-seq. (e) KEGG analysis of down-regulated genes identified for RNA-seq. (f) GO analysis of up-regulated genes identified for RNA-seq. (g) KEGG analysis of up-regulated genes identified for RNA-seq. Figure S4. Quantification of rescue experiments.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Zhou, S., Yang, K., Chen, S. et al. CCL3 secreted by hepatocytes promotes the metastasis of intrahepatic cholangiocarcinoma by VIRMA-mediated N6-methyladenosine (m6A) modification.
J Transl Med21, 43 (2023). https://doi.org/10.1186/s12967-023-03897-y