Epigenetic dysregulation-mediated COL12A1 upregulation predicts worse outcome in intrahepatic cholangiocarcinoma patients
Clinical Epigenetics volume 15, Article number: 13 (2023)
Collagen type XII alpha 1 chain (COL12A1) is associated with human cancer progression. Nevertheless, the expression pattern and the function of COL12A1 in intrahepatic cholangiocarcinoma (iCCA) remain unknown. The present study was performed to assess the role of COL12A1 in iCCA.
A total of 1669 genes, differentially expressed between iCCA and nontumor liver tissue samples, were identified as potential tumor-specific biomarkers for iCCA patients. Of these, COL12A1 was significantly upregulated in clinical iCCA tissue samples and correlated with epithelial–mesenchymal transition gene set enrichment score and advanced tumor stage in clinical iCCA. COL12A1-high expression was associated with the poor prognoses of iCCA patients (n = 421) from four independent cohorts. Promoter hypermethylation-induced downregulation of miR-424-5p resulted in COL12A1 upregulation in clinical iCCA. Experimental knockout of COL12A1 inhibited the proliferation, invasiveness and growth of iCCA cells. MiR-424-5p had a therapeutic potential in iCCA via directly targeting COL12A1.
Promoter hypermethylation-induced miR-424-5p downregulation contributes to COL12A1 upregulation in iCCA. COL12A1 is a promising druggable target for epigenetic therapy of iCCA.
Intrahepatic cholangiocarcinoma (iCCA) is a highly aggressive malignancy with dismal prognosis [1,2,3]. Its incidence and mortality are growing worldwide in recent years . Due to lack of disease-specific biomarker for early detection of iCCA, most iCCA patients are typically diagnosed at an advanced stage when the currently available systemic therapies have limited effectiveness . Therefore, identification of tumor-specific biomarkers with druggable targets is urgently needed to discover new therapeutic agents for personalized therapy of iCCA.
Genome and transcriptome sequencing had provided valuable insights into the molecular etiologies of iCCA in the past decade [6,7,8,9]. Accordingly, several molecular inhibitors binding to fibroblast growth factor receptor (FGFR) 2 or isocitrate dehydrogenase 1 or 2 (IDH1 or 2) were developed for management of advanced iCCA carrying FGFR2 fusion or IDH1/2 mutation . Recent research provided new insights into molecular pathogenesis of iCCA at proteomic scale . Nevertheless, there is still lack of research works to reveal common hallmarks and tumor-specific molecular events in iCCA from different omics datasets.
To reveal common molecular signatures involved in iCCA development, we performed an integrative analysis of multiple omics datasets including liver cancer samples from different centers. Furthermore, among the identified potential tumor-specific biomarkers, we investigated the expression pattern, clinical relevance, biological function and upstream regulatory mechanism of collagen type XII alpha 1 chain (COL12A1) in iCCA.
Identification of potential tumor-specific mRNA biomarkers for iCCA patients
Differential analysis showed that 2696 human genes were aberrantly expressed across three independent iCCA datasets (Fig. 1A). Of these, 1027 genes were also aberrantly expressed in hepatocellular carcinoma (HCC) samples from TCGA dataset, and the remaining 1669 genes were identified as potential iCCA-specific mRNA biomarkers (Fig. 1A). Weighted correlation network analysis (WGCNA) showed that the expression profiles of these potential tumor-specific biomarkers clustered into 6 gene expression modules (Fig. 1B). Among these, we are interested in the brown module (MEbrown) as it was positively correlated with necrosis detected in iCCA tumor and inversely correlated overall survival time of iCCA patients from E-MTAB-6389 cohort (Fig. 1C). Cluster analysis of a total of 76 genes from brown module showed that, except for 7 mRNA biomarkers (LRRC25, ZDHHC14, OVGP1, MRPS28, MTMR4, GRTP1 and L2HGDH) downregulated in iCCA, the remaining detected mRNA biomarkers were markedly upregulated in iCCA tissue samples relative to nontumor liver tissue samples from E-MTAB-6389 cohort (Fig. 1D).
Enrichment analysis showed that the identified mRNA biomarkers are significantly enriched in 7 main biological processes including natural killer cell mediated cytotoxicity, aminoglycan catabolic process, fatty acid derivative biosynthetic process, specification of animal organ identity, chondroitin sulfate metabolic process, amino sugar metabolic process and negative regulation of reproductive process (Fig. 1E). Volcano plot showing that COL12A1, IL8, CD248, ANTXR1 and THBS2 were the top five upregulated genes from brown module (Fig. 1F). Given iCCA is a collagen-rich malignancy, and the expression fold change of COL12A1 ranked first among the identified differentially expressed genes from brown module (Fig. 1F). We subsequently investigated and validated the expression profile, clinical significance, upstream regulatory mechanism and the biological role of COL12A1 in iCCA.
The upregulation of COL12A1 in clinical iCCA tissues
Reanalysis of multiple transcriptomic datasets revealed that COL12A1 mRNA was significantly upregulated in clinical iCCA samples relative to the adjacent nontumor liver tissue or HCC tissue samples (Fig. 2A, Additional file 1: Fig. S1A). Consistent with the findings in omics datasets, immunoblotting and RT-PCR confirmed that COL12A1 mRNA and protein were significantly overexpressed in clinical iCCA tissue specimens (Fig. 2B). In addition, COL12A1 mRNA and protein were upregulated in the indicated human iCCA cells in comparison with normal intrahepatic bile duct cells (Fig. 2C). Furthermore, COL12A1 mRNA and protein expression was similar between clinical HCC tissue and nontumor liver tissue samples (Fig. 2B, Additional file 1: Fig. S1B), whereas COL12A1 protein was downregulated in human HCC cell lines relative to normal liver cell line HL7702 (Fig. S1C). IHC and immunofluorescence staining showed that COL12A1 protein was expressed in both cytoplasm and nucleus of iCCA cells, as well as the extracellular matrix of clinical iCCA tissue (Fig. 2D, E).
COL12A1 upregulation associated with poor clinical outcomes in iCCA patients
COL12A1 expression was markedly increased in iCCA patients at the advanced stage from TCGA, OEP001105 or LVC1202 cohorts (Additional file 7: Table S1, Fig. 3A, B). Survival analyses showed that iCCA patients with high COL12A1 mRNA expression had a significant worse prognosis than iCCA patients with low COL12A1 mRNA expression from TCGA, E-MTAB-6389 or OEP001105 cohort (Fig. 3C), and the corresponding optimal cutoff value for dividing iCCA patients into COL12A1 mRNA-high/low expression (log2 expression) group was 11.375 in TCGA cohort, 6.02 in E-MTAB-6389 cohort and 5.255 in OEP001105 cohort, respectively. The unfavorable impact of COL12A1 upregulation on the prognosis of iCCA patients was also observed at the translational level in two independent iCCA cohorts (Fig. 3D, Additional file 2: Fig. S2), and the optimal cutoff point for stratifying iCCA patients into two groups was -0.6726799 (log2 ratio of unique peptides) in OEP001105 cohort and 3 (IHC score of COL12A1 staining) in LVC1202 cohort, respectively. The median overall survival time was significantly shorter in iCCA patients with COL12A1-high expression than iCCA patient with COL12A1-low expression from TCGA, OEP001105, E-MTAB-6389 and LVC1202 cohorts (18.55 (10.16–29.54) months vs. 53.45 (29.62–75.53) months, p = 0.014, Fig. 3E). Furthermore, the recurrence risk of iCCA was higher in COL12A1-high expression (COL12A1 staining score ≥ 4) group (Fig. 3F). Additionally, multivariate analyses determined that COL12A1-high expression was an independent risk factor associated with the poor prognosis of iCCA patients (Fig. 3G).
COL12A1 promoter methylation pattern in iCCA and the corresponding nontumor samples
To investigate the upstream mechanism of COL12A1 upregulation in iCCA, we first deciphered the methylation pattern of COL12A1 promoter/enhancer regions by analysis of two independent DNA methylation datasets of clinical iCCA (Additional file 3: Fig. S3). The methylation degree of COL12A1 enhancer region was similar between iCCA tissues and the adjacent nontumor liver tissues in both TCGA and GSE156299 datasets (Fig. 4A, B); the methylation degree of COL12A1 promoter region was significantly increased in iCCA tissues relative to the paired nontumor liver tissues in TCGA dataset (Fig. 4A), while no significant difference was observed in the samples included in GSE156299 dataset (Fig. 4B). In addition, the methylation degree of several regions in COL12A1 locus including cg03564793, cg14375912, cg04074140 and cg19132213 was significantly increased in iCCA tissue samples compared with the adjacent nontumor samples (Additional file 4: Fig. S4), indicating that the methylation pattern of COL12A1 promoter is not the cause for its upregulation in human iCCA samples.
MiR-424-5p directly binding to COL12A1 mRNA in iCCA cells
TargetScanHuman database (version 8.0) prediction revealed that 7 candidate miRNAs (miR-424-5p, -15a-5p, -497-5p, -15b-5p, -6838-5p, -16-5p and miR-195-5p) can bind to COL12A1 3'UTR (Additional file 5: Fig. S5). Both sample-paired and unpaired analyses determined that miR-424-5p expression was significantly decreased in clinical iCCA relative to the adjacent nontumor samples (Fig. 4C and Additional file 6: Fig. S6A). Except for miR-6838-5p that was not detected in TCGA dataset, the remaining candidate miRNAs expression levels were similar between clinical iCCA and the adjacent nontumor samples (Fig. 4C and Additional file 6: Fig. S6A). Correlation analyses showed that miR-424-5p was significantly negatively correlated with COL12A1 mRNA expression in clinical iCCA samples from TCGA (Additional file 6: Fig. S6B, C). Additionally, miR-424-5p was markedly downregulated in human iCCA cells compared with normal intrahepatic bile duct cell HIBEpiC (Additional file 6: Fig. S6D).
Immunoblotting showed that overexpression of miR-424-5p suppressed COL12A1 expression in HuCCT1 and RBE cells (Fig. 4D). Luciferase reporter assays showed that luciferase activity was markedly reduced or increased in HEK-293T cells co-transfected with miR-424-5p plus wild-type COL12A1 3'UTR or mutated 3'UTR of COL12A1, respectively (Fig. 4E, F). RT-PCR showed that overexpression of miR-424-5p downregulated COL12A1 mRNA expression in HuCCT1 and RBE cells (Fig. 4G). In addition, the downregulation of miR-424-5p in clinical iCCA and the inverse association of COL12A1 mRNA with miR-424-5p expression were validated in our internal iCCA dataset (Fig. 4H). These findings suggest that miR-424-5p can result in translational inhibition or degradation of COL12A1 by directly targeting COL12A1 mRNA.
Promoter hypermethylation-mediated miR-424-5p downregulation in clinical iCCA
To decipher the upstream regulatory mechanism for miR-424-5p downregulation in clinical iCCA, we tested the promoter methylation pattern of miR-424 in 10 pairs of clinical iCCA samples by targeted bisulfite sequencing. The results showed that CpG sites methylation degree of miR-424 promoter area was markedly increased in iCCA compared with the corresponding paired nontumor samples (Fig. 5A). The promoter methylation degree of miR-424 was significantly inversely associated with miR-424-5p expression in iCCA (Fig. 5B). In addition, DNA methyltransferase writers (DNMTs) were inversely correlated with miR-424 family members at the transcriptional level (Fig. 5C, D). These results indicate that DNMTs-mediated promoter CpG sites hypermethylation results in miR-424-5p downregulation in clinical iCCA.
We treated iCCA cells with decitabine to assess whether DNA demethylating agents can inhibit COL12A1 expression via inducing miR-424 promoter demethylation. Immunoblotting showed that COL12A1 expression was decreased in HuCCT1 and CCLP1 cells after being treated with the indicated dose of decitabine (Fig. 5E). Targeted bisulfite sequencing showed that CpG sites methylation degree of miR-424 promoter area was decreased after decitabine treatment in HuCCT1 and CCLP1 cells (Fig. 5F). RT-PCR determined that miR-424-5p expression was increased whereas COL12A1 mRNA expression was decreased in HuCCT1 and CCLP1 cells treated with decitabine (Fig. 5G, H). Moreover, COL12A1 expression levels were lower in subcutaneous xenograft tumors from the decitabine treatment group compared with the controls (Fig. 5I). Collectively, these findings showing that DNMTs-mediated miR-424-5p downregulation resulted in COL12A1 upregulation in iCCA cells.
COL12A1 promoting proliferation and growth of iCCA cells in vitro and in vivo
To assess the impact of COL12A1 on iCCA cell progression, we silenced COL12A1 expression in HuCCT1 and RBE cells using COL12A1 CRISPR/Cas9 KO plasmid, respectively (Fig. 6A). In cell viability assays, COL12A1 knockout weakened HuCCT1 and RBE cells proliferation (Fig. 6B). In colony-forming assays, COL12A1 knockout inhibited tumor cells growth by attenuating colony-forming ability of iCCA cells (Fig. 6C). In subcutaneous xenograft models, COL12A1 knockout inhibited tumor growth and proliferation (Fig. 6D, E).
COL12A1 positively correlated with the epithelial–mesenchymal transition (EMT) program in clinical iCCA
Like the EMT-associated marker vimentin (VIM), COL12A1 was significantly positively correlated with EMT-gene set enrichment score in 3 independent iCCA datasets (Fig. 7A, B). Immunoblotting showed that COL12A1 knockout increased the epithelial marker E-cadherin expression whereas reduced mesenchymal markers (N-cadherin and vimentin) expression in HuCCT1 and RBE cells (Fig. 7C). These findings indicate that COL12A1 upregulation plays a key role in maintaining the EMT program of iCCA cells.
The therapeutic value of miR-424-5p in iCCA by sponging COL12A1
MiR-424-5p overexpression can inhibit the growth, invasion and migration of iCCA cell in vitro . However, it is still unclear whether miR-424-5p has antitumor potential in vivo. To evaluate the treatment effect of miR-424-5p in iCCA, we treated the tumor-bearing mice by tail vein administration of miR-424-5p agonist (AgomiR, 20 nmol per mouse). In the xenograft model, tumor growth was markedly inhibited in AgomiR treatment group relative to the controls (Fig. 8A). COL12A1 expression was decreased in the AgomiR group compared with the controls (Fig. 8B).
In the present study, we identified 1669 potential iCCA-specific biomarkers by reanalysis of multiple transcriptomic datasets of iCCA and HCC. These results will contribute to more precisely revealing tumor-specific molecular signatures involved in iCCA initiation and/or progression and discovering potential therapeutic targets for iCCA patients. The omics data provided valuable resources for determining the common and distinct events at genomic, epigenetic, transcriptomic and proteomic scale in iCCA [6,7,8,9, 11, 13, 14]. By integrative analysis of transcriptomics and clinical profiles of iCCA patients, we are interested in aberrantly expressed genes clustered in gene expression module correlated with necrosis that occurred in tumor and overall survival of iCCA patients. Based on the results from enrichment analyses and the marked collagen-rich tumor microenvironment in iCCA, we subsequently investigated the role of COL12A1 in iCCA.
COL12A1 is aberrantly upregulated in human malignancies including gastric cancer, colon cancer, colorectal cancer and pancreatic adenocarcinoma, and its upregulation correlated with tumor progression and poor prognosis [15,16,17,18]. Consistent with previous findings, the present research showed that COL12A1 is significantly upregulated in clinical iCCA relative to the nontumor liver tissues or HCC samples. Higher expression of COL12A1 was associated with tumor progression and worse survival in iCCA patients. Moreover, our loss-of-function experiments demonstrated that COL12A1 upregulation promotes iCCA cells proliferation and growth in vitro and in vivo. Thus, COL12A1 is a promising biomarker for the prognosis prediction of iCCA patients.
EMT program plays a pivotal role in the growth and metastasis of cancers including iCCA [5, 19]. In the present study, we first revealed that, like vimentin (a classic mesenchymal marker), COL12A1 is positively correlated with EMT-gene set enrichment scores in clinical iCCA samples, indicating that COL12A1 mRNA may be served as a new EMT maker for clinical iCCA. Additionally, we found that COL12A1 can regulate EMT biomarkers expression at the protein level in human iCCA cells, which suggests that COL12A1 regulates the growth and metastasis process of iCCA cells possibly via modulating the EMT program.
Accumulating evidence suggests that epigenetic dysregulation modulates liver cancer development by regulating tumor-associated gene expression [20,21,22,23]. Unlike promoter/enhancer demethylation-mediated upregulation of ITPR3, MCM2 and NUP37 in hepatocellular carcinoma [20, 21]. In this study, we first revealed that DNA methylation degree of COL12A1 promoter/enhancer was not demethylated in iCCA relative to the adjacent nontumor samples, indicating that COL12A1 upregulation is independent of the methylation level of its promoter/enhancer region in clinical iCCA.
MicroRNA dysregulation plays a pivotal role in human diseases, particularly in cancers [12, 24]. In the present study, we first revealed that miR-424-5p is an upstream regulator of COL12A1 in iCCA. It has been demonstrated that miR-424-5p was downregulated in several kinds of cancer and had antitumor functions in breast cancer and iCCA cells [12, 25]. Consistent with the previous finding, our results also showed that miR-424-5p was markedly downregulated in clinical iCCA tissues relative to the paired nontumor liver tissues.
MiR-424-5p can regulate Hippo/YAP1 and Wnt/ꞵ-catenin signaling via targeting circACTN4 and YAP1 in iCCA cells , showing that miR-424-5p plays a crucial role in iCCA initiation and/or progression and may have therapeutic potential via inhibiting Hippo/YAP1 and Wnt/ꞵ-catenin signaling in iCCA. In our study, we observed that miR-424-5p had therapeutic potential in the xenograft iCCA model via targeting COL12A1. These results elucidated the pivotal role of miR-424-5p in iCCA progression and provided preclinical evidence for further translational studies to investigate the therapeutic effect of miR-424-5p in iCCA patients.
There is still a lack of research to decipher the epigenetic regulatory mechanism of miR-424-5p downregulation in human malignancies. In this study, we first demonstrated that promoter CpG sites of miR-424 were aberrantly hypermethylated in iCCA, which induced miR-424-5p downregulation in clinical iCCA. Furthermore, DNMTs-dictated aberrant promoter hypermethylation silenced miR-424-5p in iCCA cells. Although we first revealed the epigenetic mechanism of COL12A1 upregulation in clinical iCCA and determined the protumor function of COL12A1 in human iCCA cells, there are also some limitations in our present study. First, we only established HuCCT1 cell-derived xenograft mouse models to assess the role of COL12A1 in vivo, there may be heterogeneities among different iCCA models such as patient-derived xenograft and hydrodynamic model of iCCA. Second, we did not evaluate the therapeutic effect of miR-424-5p mimic in iCCA patients via targeting COL12A1 since there is still a lack of public evidence to make a clinical trial. Future studies are needed to explore the detailed mechanisms downstream of COL12A1 in iCCA progression and develop potential inhibitor-targeting COL12A1 mRNA or protein for in vivo use.
COL12A1 is upregulated in iCCA by promoter CpG sites hypermethylation-induced downregulation of miR-424-5p, which promotes iCCA progression. MiR-424-5p has therapeutic potential in iCCA by targeting COL12A1. Our research determined the pivotal role of COL12A1 in iCCA progression and indicated that COL12A1 is a promising druggable target for epigenetic therapy of iCCA. Therefore, it is reasonable for use of miR-424-5p mimic or potential inhibitor-targeting COL12A1 to manage iCCA progression or postoperative recurrence. This hypothesis should be tested in future translational studies enrolling iCCA patients.
Public omics datasets
Public bulk-RNA seq datasets (TCGA-CHOL , TCGA-LIHC , GSE76297 , GSE26566 , E-MTAB-6389 , GSE32879 , E-MTAB-4171 , OEP001105  and OEP000321 ), including gene expression profile of 565 intrahepatic cholangiocarcinoma (iCCA), 607 hepatocellular carcinoma (HCC), 7 combined hepatocellular carcinoma (CHC), 5 focal nodular hyperplasia (FNH) and 327 adjacent nontumor liver tissue specimens included in The Cancer Genome Atlas (TCGA), ArrayExpress, Gene Expression Omnibus (GEO), and National Omics Data Encyclopedia (NODE) were obtained and processed using R package TCGAbiolinks (version 2.20.1), oligo (version 1.56.0) and GEOquery (version 2.60.0), respectively [31,32,33]. Differentially expressed genes between nontumor liver tissue and iCCA or HCC samples were inferred using R package limma (version 3.50.3) based on |logFC| > 0.1 and p value < 0.05 .
Re-analyses of DNA methylation datasets GSE156299  and TCGA-CHOL  were performed to decipher the methylation patterns of COL12A1 gene in iCCA and nontumor liver tissues. The raw data of DNA methylation were processed using minifi (version 1.38.0) as previously described [21, 36]. Weighted gene co-expression network analysis (WGCNA) was performed for weighted correlation network analysis of gene expression patterns in iCCA .
Single-sample gene set enrichment was performed using R package GSVA (version1.42.0) and GSEABase (version 1.56.0), and epithelia-to-mesenchymal transition (EMT) gene sets (GOBP_EPITHELIAL_TO_MESENCHYMAL_TRANSITION) in Molecular Signatures Database c5.all.v7.5.1.symbols.gmt was applied to calculate the EMT-gene set enrichment score in each sample .
Patients and clinical tissue specimen
A total of 10 iCCA, 10 HCC, and the paired nontumor liver tissue samples were prospectively collected from treatment-naive patients with iCCA or HCC post radical surgical treatment from May to September 2020 at our hospital. The clinicopathological characteristics of 10 iCCA/HCC patients are summarized in Additional file 8: Table S2. After snap-frozen in liquid nitrogen, the collected tissue samples were stored at − 80 refrigerator for subsequent molecular biology experiments as we previously described . An external cohort (LVC1202, SUPERBIOTEK, Shanghai, CN), including 60 primary iCCA patients without receiving neoadjuvant treatment prior to primary surgery, was obtained for survival analysis. The TNM staging of each iCCA patient after surgical resection was performed in accordance with the 8th edition of the American Joint Committee on Cancer (AJCC) TNM classification system . The pathological diagnosis of iCCA or HCC was performed according to World Health Organization criteria . Written informed consent was collected from each patient enrolled in the present study. The experimental protocol (5826Z/20, 0531B/22) for using clinical tissue specimens was reviewed and approved by Institutional Ethics Committee in our hospital.
Cell lines, cell culture, stable cell line establishment and transient infection
RBE cell line was obtained from Yuchi Biological Technology Co., Ltd. (Shanghai, CN). HuCCT1 (JCRB0425) and HCCC9810 cell lines were kindly provided by the Key Laboratory of Carcinogenesis and Cancer Invasion (Fudan University), Ministry of Education, Shanghai, CN. CCLP1 and normal human intrahepatic biliary cell (HIBEpiC) lines were kindly provided by Key Laboratory of Combined Multiorgan Transplantation, National Health Commission, Hangzhou, Zhejiang, CN. HEK-293T cell line was purchased from the Cell Bank of Chinese Academy of Sciences (Shanghai, CN). RBE cell line came from a 37-year-old yellow race woman with primary iCCA. HCCC9810 cell line came from a 60-year-old yellow race woman with primary iCCA. CCLP1 cell line came from a 48-year-old Caucasian female with primary iCCA. HuCCT1 cell line came from a 56-year-old yellow race man with primary iCCA. HIBEpiC cells were isolated from normal liver tissue. All cell lines were tested by short tandem repeat typing. Furthermore, all cells were cultured in RPMI-1640 medium (Thermo Fisher Scientific, USA) containing 10% fetal bovine serum (FBS) (Thermo Fisher Scientific, USA), 100 U/mL penicillin (Cienry, Huzhou, CN) and 100 μg/mL streptomycin (Cienry, Huzhou, CN) at 37 °C in 5% CO2 atmosphere. Human hepatoma cell lines (including Hep3B, Bel-7402, HCC-LM3 and Huh7) and normal liver cell line HL7702 were cultured as we previously described . Additionally, all of the cell lines had been tested regularly for mycoplasma contamination as previously described .
COL12A1 CRISPR/Cas9 KO plasmid (sc-402361) was purchased from Santa Cruz Biotechnology (USA). 2 × 105 HuCCT1 or RBE cells (per well) in the logarithmic phase were seeded in a 6-well tissue culture plate and incubated with antibiotic-free standard growth medium (3 mL per well). Upon cells growing to a 50–70% confluency after initial seeding, 2 μg (0.1 μg/μL) of plasmid DNA diluted in 130 μL of Plasmid Transfection Medium (sc-108062, Santa Cruz Biotechnology, USA) and 10 μL of UltraCruz® Transfection Reagent (sc-395739, Santa Cruz Biotechnology, USA) diluted in 140 μL of Plasmid Transfection Medium, were transfected into HuCCT1 and RBE cells following the manufacturer’s protocol. After incubation for 48 h, cells with successful transfection of COL12A1 KO plasmid were determined by detecting the green fluorescent protein (GFP) using fluorescent microscopy (Leica DMi8, Germany). Then, after cell sorting for GFP by FACS (MoFloAstrios EQ, Beckman), the harvested cells were cultured in a standard growth medium for 3–4 weeks and the selected single-cell colony with COL12A1 knockout was then used for the downstream experiments.
The mimic of miRNAs including miR-424-5p (miR10001341-1-5), miR-424-3p (miR10004749-1-5), miR-497-5p (miR10002820-1-5), miR-195-5p (miR10000461-1-5) and the negative control (miR-NC, miR1N0000001-1-10) were purchased from RIBOBIO Co. Ltd. (Guangzhou, CN). Transient transfection with the mimic of miRNAs was performed using riboFECT CP Transfection Kit (C10511-05, RiboBio Co. Ltd., Guangzhou, CN) in accordance with the manufacturer's protocol.
DNA methylation inhibitor decitabine (S1200) was purchased from Selleck (Shanghai, CN). 3 × 105 HuCCT1 or CCLP1 cells (per well) were seeded in a 6-well cell culture plate and incubated for 48 h in RPMI-1640 containing 10% FBS and 10/20 μmol/mL decitabine.
Protein extraction and immunoblotting
Protein extraction from fresh clinical tissue samples or human cell lines and immunoblotting procedures were performed as previously described . Samples containing 20–40 μg proteins were subjected to SDS-PAGE (8–10% polyacrylamide) and blotted on polyvinylidene difluoride (PVDF) membranes. After blocking with 5% skim milk, the membranes were first incubated with COL12A1 antibody (1:2000, Sigma-Aldrich, #HPA009143), Vimentin antibody (1:1000, Abcam, #ab92547), E-cadherin antibody (1:1000, Cell Signaling Technology, #3195), or N-cadherin antibody (1:1000, Cell Signaling Technology, #13116), and then incubated with Goat Anti-Rabbit IgG H&L (1:5000, Abcam, #ab205718). GAPDH (1:5000, Abcam, #ab181602) was used as a loading control.
RNA isolation and real-time PCR
Total RNA isolation from clinical tissue samples or cell lines was performed using TRIzol™ Plus RNA Purification Kit (12183555, Invitrogen™, Thermo Fisher Scientific, USA) following the instruction for the user. Real-time PCR was performed as previously described . The primers of the human COL12A1 and β-actin gene (GCD0259363) were purchased from GeneChem (Shanghai, CN). The primer sequences of COL12A1 and β-actin are as follows: COL12A1, forward 5′-AGGTCGGATGACGGGAAGA-3′ and reverse 5′-GCGGACATTCAAGGTGCTG-3′; β-actin, forward 5′-CGACAGATGCAGAACGAGA-3′ and reverse 5′-GACCCTGGATGTGACAGCTC-3′.
MicroRNAs (miRNAs) in clinical tissue samples and human cell lines were extracted using an ultrapure miRNA isolation kit (5080576001, Roche) following the instructions for the user. The cDNA template was produced using Bulge-Loop miRNA qRT-PCR Starter Kit (C10211-2, RIBOBIO, Guangzhou, CN) following the instructions for the user. U6 snRNA was used as the endogenous control, and miRNAs primers (MQPS0001272-1-200, MQPS0001271-1-200, MQPS0000002-1-200) were purchased from RiboBio Co. Ltd. (Guangzhou, CN). A final volume of 20 μL mixture containing Bulge-Loop miRNA qRT-PCR Starter Kit (C10211-1, RIBOBIO, Guangzhou, CN) reagents, RT-PCR Grade Water (AM9935, Applied Biosystems, Thermo Fisher Scientific, USA) and cDNA template was amplified using ABI QuantStudio-5 Real-Time PCR System (Applied Biosystems, Thermo Fisher Scientific, USA) in accordance with the manual for users. Reaction conditions were set as 95 °C for 10 min, 40 cycles of denaturation at 95 °C for 2 s, annealing at 60 °C for 20 s and extension at 70 °C for 10 s. The ΔΔCT method was used to determine the relative miRNA expression level in each sample.
Hematoxylin and eosin staining and immunohistochemistry
Immunohistochemical (IHC) or Hematoxylin and eosin (H&E) staining was performed as previously described [2, 21]. Paraffin-embedded liver tissue sections (3 μm thickness) were incubated with anti-COL12A1 (1:200, Sigma #HPA009143) or anti-Ki67 (1:300, Abcam, #ab16667), respectively. The COL12A1 staining score was used to divide iCCA patients into COL12A1-low and COL12A1-high expression groups based on an optimal threshold.
A total of 2 × 103 iCCA cells, seeded on coverslips coated in a 24-well cell culture plate, were incubated in RPMI-1640 medium containing 10% FBS and allowed to adhere for 6–8 h. After 3 times washing using phosphate-buffered saline (PBS), cells were fixed with 4% paraformaldehyde for 10 min at room temperature, then washed 3 times with ice-cold PBS, permeabilized with 0.1% Triton-X100 in PBS for 10 min, blocked with 1% bovine serum albumin (BSA) and 22.52 mg/mL glycine in PBS supplemented with 0.1% Tween 20 (PBST) at room temperature for 30 min, and finally incubated with anti-COL12A1 antibody (1:50) overnight at 4 °C in the dark. After 3 times washing using PBS, cells were then incubated with fluorescence-labeled secondary antibody (1:200, Abcam, #ab150079) for 1 h in the dark, washed 3 times with PBS, followed by incubation with 0.1 μg/mL 4′,6-diamidino-2-phenylindole (DAPI) at room temperature for 1 min in the dark, and finally rinsed 3 times with PBS in the dark. Coverslips were mounted with UltraCruz® Aqueous Mounting Medium (sc-24941, Santa Cruz) and sealed with nail oil. Finally, cells on the coverslip were microphotographed using a Leica fluorescence microscope (TCS SP8, Germany).
Cell viability assay
We assessed iCCA cell viability using CCK8 assays (Beyotime, Shanghai, CN) as previously described . Briefly, 1 × 103 iCCA cells (per well) were seeded into 96-well plates and incubated for 24 h and 48 h, respectively. Then, CCK8 (1:10) combined with fresh growth medium was added to each well. After incubation for 1 h, the absorbance at 450 nm was tested as previously described .
Cell colony-forming assay
1 × 103 iCCA cells in the logarithmic phase were seeded into a 6-cm cell culture dish and incubated for 2 weeks. After 3 times washing with PBS, the cell colonies were stained with 0.05% crystal violet (Sigma-Aldrich, USA) for 30 min, then photographed and counted using ImageJ (version.1.8.0), which were performed at room temperature.
Dual-luciferase reporter assays
The sense and antisense fragments of the human COL12A1 3'UTR including either the wild-type (GOSE0310921) or the mutated sites (GOSE0310922) binding to miR-424-5p were synthesized separately by GeneChem (Shanghai, CN). Then the respective sense and corresponding antisense strands were hybridized and then cloned into the GV272-Firefly-Luciferase vector at the XhoI and KpnI sites. 2 × 104 HEK-293 T cells (per well) were seeded into a 96-well cell culture plate, and co-transfected with either the GV272 vector (Control), GV272-COL12A1-3'UTR (wild-type), or the GV272-COL12A1-Mut-3'UTR, and the miR-424-5p mimic and CV045-TK promoter-Renilla-Luciferase plasmid (GeneChem, Shanghai, CN) using LipofectamineTM3000 transfection reagent (Invitrogen, Thermo Fisher Scientific, USA) following the instructions for the user. After incubation for 48 h at 37 °C in 5% CO2, the luciferase activity of cell lysate was tested using a Dual-Glo luciferase assay kit (E1910, Promega). Luminescence was measured by a SpectraMax i3x Multi-model Microplate reader and the SoftMax Pro7 (Molecular Devices, USA).
Genomic DNA isolation and targeted bisulfite sequencing
The extraction of genomic DNA in fresh-frozen iCCA and the paired nontumor liver tissue samples (n = 10), as well as human iCCA cells, was performed as we previously described . The concentration of each extracted DNA sample was fluorometrically tested using Qubit 4.0 (Invitrogen, Thermo Fisher Scientific). Bisulfite conversion and next-generation sequencing (NGS) library preparation were performed as previously described . According to the miRNA promoter recognition methods described by Annalisa Marsico et al. , we identified miR-424 promoter locus (chrX: 133683002–133684302) and designed primers for miR-424 promoter region using MethPrimer (version 2.0.0) [43, 44]. The putative primers for miR-424 promoter region (Additional file 9: Table S3) were validated and synthesized by BGI genomics (Shenzhen, CN). A total of 1 µg genomic DNA was subjected to bisulfite DNA conversion with ZYMO EZ DNA Methylation-Gold Kit (#D5006, Zymo Research, Irvine, CA, USA) accordingly to the manufacturer’s instructions. Then, a total of 500 ng bisulfite-converted DNA sample was used for PCR amplification using KAPA HiFi HotStart Uracil+ ReadyMix PCR Kit (#KK2801, Kapa Biosystems, Wilmington, MA, USA) combined with bisulfite-specific primers according to the standard conditions.
After the quality control using Labchip GX Touch (PerkinElmer, Germany), one-twentieth of the elution PCR reaction products were subjected to second round PCR for ligation of barcoded sequencing adapter with T4 DNA ligase. Thereafter, all samples were quantified on an Agilent 2100 (Agilent, USA) and adjusted for the same copy number for subsequent sequencing, and the sequencing was performed on a NovaSeq6000 system (Illumina, USA) using the NovaSeq 6000 S4 Reagent Kit v1.5 (300 cycles) in collaboration with Yunbios (Shanghai, CN).
The basic statistics on the quality of the raw reads was performed using The FastQC tool (http://www.bioinformatics.babraham.ac.uk/projects/fastq). The adapters and low-quality reads were removed using Trimmomatic v0.36 with default parameter . The alignment of the bisulfite sequence to the human reference genome (hg19) was performed using the BSMAP v2.7.3 with parameters ‘-n 0 -g 0 -v 0.08 -m 50 -× 1000’ . The unique mapped reads were kept for the subsequent analyses and the sequence depths of methylated cytosines were greater than 20X.
Experimental animal model
The animal experiments were performed as previously described . The experimental protocol (1615/22) for all animal experiments was reviewed and approved by the Institutional Animal Care & Use Committee (IACUC). Five-week-old male NSG mice used in the study were purchased from Model Animal Research Centre of Nanjing University, Nanjing, CN.
To determine the direct impact of COL12A1 on iCCA progression, 4 × 106 HuCCT1 cells with or without COL12A1 knockout, resuspended in 200 μL PBS, were subcutaneously implanted into the flanks of 5-week-old male NSG mice. After observations for 8 weeks for tumor formation, mice were killed according to IACUC protocol. The collected tumor lesions were measured and weighted by a digital caliper and an electronic balance, respectively, and then fixed tumor tissue with 4% paraformaldehyde for use in downstream experiments. The tumor volume was calculated based on the formula as previously described .
Decitabine drug treatment model: HuCCT1 cells (3 × 106) were implanted into the flanks of 5-week-old male NSG mice. Upon tumor reaching 25 to 50 mm3, mice were randomly divided into two groups with each group including 5 mice. Then, mice were treated with decitabine (2 mg/kg, 5 times weekly) or PBS by tail vein injection.
MiR-424-5p agonist (modified miR-424-5p mimic) administration model: After implantation of 3 × 106 HuCCT1 cells into the flanks of 5-week-old male NSG mice, tumor development was monitored 3 times weekly at 3wks after implantation. Upon tumor reaching 130–200 mm3, mice were randomly divided into AgoMIR (20 nmol, i.v. 2 times weekly) group or control (100ul PBS, i.v. 2 times weekly) group, with each group consisting of 3 mice. AgoMIR (miR40001341-4-5) was purchased from RIBOBIO Co. Ltd. (Guangzhou, CN).
All statistical analyses were performed with R v.4.1.0 (R Foundation for Statistical Computing, Vienna, Austria) or with SPSS v.22.0 (IBM Corp., Armonk, NY, USA). Unless otherwise indicated, data are expressed as means ± standard error of the mean (SEM). A comparison of two groups was made using a two-tailed Student’s t test or Wilcoxon test. A comparison of multiple groups was made using a one-way analysis of variance (ANOVA). Correlation analysis between two variables was assessed by Pearson correlation analysis. The association between COL12A1 expression and the clinicopathological features of iCCA patients was assessed by chi-square test or Fisher's exact test. Survival analysis was performed with the Kaplan–Meier method, and the statistical difference was determined by log-rank test. The “surv_cutpoint ()” function in R package survminer (https://cran.r-project.org/web/packages/survminer/index.html) was applied to determine the optimal cutoff point of COL12A1 mRNA/protein expression level corresponding to the survival of iCCA patients. Cox proportional hazards regression model was applied for multivariate survival analysis.
Availability of data and materials
All data are included in the main text or online supplementary materials.
Collagen type XII alpha 1 chain
The Cancer Genome Atlas
Gene Expression Omnibus
Focal nodular hyperplasia
Combined hepatocellular carcinoma
Intraductal tubulopapillary neoplasm of biliary duct
Intraductal papillary neoplasm of biliary duct
Bridgewater J, Galle PR, Khan SA, Llovet JM, Park JW, Patel T, et al. Guidelines for the diagnosis and management of intrahepatic cholangiocarcinoma. J Hepatol. 2014;60(6):1268–89.
Tang Z, Yang Y, Zhang J, Fu W, Lin Y, Su G, et al. Quantitative proteomic analysis and evaluation of the potential prognostic biomarkers in cholangiocarcinoma. J Cancer. 2019;10(17):3985–99.
Tang Z, Yang Y, Wang X, Meng W, Li X. Meta-analysis of the diagnostic value of Wisteria floribunda agglutinin-sialylated mucin1 and the prognostic role of mucin1 in human cholangiocarcinoma. BMJ Open. 2019;9(1):e021693.
Sirica AE, Gores GJ, Groopman JD, Selaru FM, Strazzabosco M, Wei Wang X, et al. Intrahepatic cholangiocarcinoma: continuing challenges and translational advances. Hepatology. 2019;69(4):1803–15.
Banales JM, Marin JJG, Lamarca A, Rodrigues PM, Khan SA, Roberts LR, et al. Cholangiocarcinoma 2020: the next horizon in mechanisms and management. Nat Rev Gastroenterol Hepatol. 2020;17(9):557–88.
Farshidfar F, Zheng S, Gingras MC, Newton Y, Shih J, Robertson AG, et al. Integrative genomic analysis of cholangiocarcinoma identifies distinct IDH-mutant molecular profiles. Cell Rep. 2017;19(13):2878–80.
Jiao Y, Pawlik TM, Anders RA, Selaru FM, Streppel MM, Lucas DJ, et al. Exome sequencing identifies frequent inactivating mutations in BAP1, ARID1A and PBRM1 in intrahepatic cholangiocarcinomas. Nat Genet. 2013;45(12):1470–3.
Chaisaingmongkol J, Budhu A, Dang H, Rabibhadana S, Pupacdi B, Kwon SM, et al. Common molecular subtypes among Asian hepatocellular carcinoma and cholangiocarcinoma. Cancer Cell. 2017;32(1):57-70.e3.
Job S, Rapoud D, Dos Santos A, Gonzalez P, Desterke C, Pascal G, et al. Identification of four immune subtypes characterized by distinct composition and functions of tumor microenvironment in intrahepatic cholangiocarcinoma. Hepatology. 2020;72(3):965–81.
Kelley RK, Bridgewater J, Gores GJ, Zhu AX. Systemic therapies for intrahepatic cholangiocarcinoma. J Hepatol. 2020;72(2):353–63.
Dong L, Lu D, Chen R, Lin Y, Zhu H, Zhang Z, et al. Proteogenomic characterization identifies clinically relevant subgroups of intrahepatic cholangiocarcinoma. Cancer Cell. 2022;40(1):70-87.e15.
Chen Q, Wang H, Li Z, Li F, Liang L, Zou Y, et al. Circular RNA ACTN4 promotes intrahepatic cholangiocarcinoma progression by recruiting YBX1 to initiate FZD7 transcription. J Hepatol. 2022;76(1):135–47.
Zhang M, Yang H, Wan L, Wang Z, Wang H, Ge C, et al. Single-cell transcriptomic architecture and intercellular crosstalk of human intrahepatic cholangiocarcinoma. J Hepatol. 2020;73(5):1118–30.
Sia D, Hoshida Y, Villanueva A, Roayaie S, Ferrer J, Tabak B, et al. Integrative molecular analysis of intrahepatic cholangiocarcinoma reveals 2 classes that have different outcomes. Gastroenterology. 2013;144(4):829–40.
Chivu-Economescu M, Necula LG, Matei L, Dragu D, Bleotu C, Sorop A, et al. Collagen family and other matrix remodeling proteins identified by bioinformatics analysis as hub genes involved in gastric cancer progression and prognosis. Int J Mol Sci. 2022;23(6):3214.
Mikula M, Rubel T, Karczmarski J, Goryca K, Dadlez M, Ostrowski J. Integrating proteomic and transcriptomic high-throughput surveys for search of new biomarkers of colon tumors. Funct Integr Genomics. 2011;11(2):215–24.
Wu Y, Xu Y. Integrated bioinformatics analysis of expression and gene regulation network of COL12A1 in colorectal cancer. Cancer Med. 2020;9(13):4743–55.
Li Y, Su Z, Wei B, Qin M, Liang Z. Bioinformatics analysis identified MMP14 and COL12A1 as immune-related biomarkers associated with pancreatic adenocarcinoma prognosis. Math Biosci Eng. 2021;18(5):5921–42.
Dongre A, Weinberg RA. New insights into the mechanisms of epithelial–mesenchymal transition and implications for cancer. Nat Rev Mol Cell Biol. 2019;20(2):69–84.
Guerra MT, Florentino RM, Franca A, Lima Filho AC, Dos Santos ML, Fonseca RC, et al. Expression of the type 3 InsP3 receptor is a final common event in the development of hepatocellular carcinoma. Gut. 2019;68(9):1676–87.
Tang Z, Yang Y, Chen W, Li E, Liang T. Demethylation at enhancer upregulates MCM2 and NUP37 expression predicting poor survival in hepatocellular carcinoma patients. J Transl Med. 2022;20(1):49.
Hernandez-Meza G, von Felden J, Gonzalez-Kozlova EE, Garcia-Lezana T, Peix J, Portela A, et al. DNA methylation profiling of human hepatocarcinogenesis. Hepatology. 2021;74(1):183–99.
Wang Q, Liang N, Yang T, Li Y, Li J, Huang Q, et al. DNMT1-mediated methylation of BEX1 regulates stemness and tumorigenicity in liver cancer. J Hepatol. 2021;75(5):1142–53.
Rupaimoole R, Slack FJ. MicroRNA therapeutics: towards a new era for the management of cancer and other diseases. Nat Rev Drug Discov. 2017;16(3):203–22.
Dastmalchi N, Hosseinpourfeizi MA, Khojasteh SMB, Baradaran B, Safaralizadeh R. Tumor suppressive activity of miR-424-5p in breast cancer cells through targeting PD-L1 and modulating PTEN/PI3K/AKT/mTOR signaling pathway. Life Sci. 2020;259:118239.
Cancer Genome Atlas Research Network. Electronic address wbe, Cancer Genome Atlas Research N. Comprehensive and integrative genomic characterization of hepatocellular carcinoma. Cell. 2017;169(7):1327-1341 e23.
Andersen JB, Spee B, Blechacz BR, Avital I, Komuta M, Barbour A, et al. Genomic and genetic characterization of cholangiocarcinoma identifies therapeutic targets for tyrosine kinase inhibitors. Gastroenterology. 2012;142(4):1021-1031 e15.
Oishi N, Kumar MR, Roessler S, Ji J, Forgues M, Budhu A, et al. Transcriptomic profiling reveals hepatic stem-like gene signatures and interplay of miR-200c and epithelial–mesenchymal transition in intrahepatic cholangiocarcinoma. Hepatology. 2012;56(5):1792–803.
Nojima M, Matsui T, Tamori A, Kubo S, Shirabe K, Kimura K, et al. Global, cancer-specific microRNA cluster hypomethylation was functionally associated with the development of non-B non-C hepatocellular carcinoma. Mol Cancer. 2016;15(1):31.
Gao Q, Zhu H, Dong L, Shi W, Chen R, Song Z, et al. Integrated proteogenomic characterization of HBV-related hepatocellular carcinoma. Cell. 2019;179(2):561-577 e22.
Colaprico A, Silva TC, Olsen C, Garofano L, Cava C, Garolini D, et al. TCGAbiolinks: an R/Bioconductor package for integrative analysis of TCGA data. Nucleic Acids Res. 2016;44(8):e71.
Carvalho BS, Irizarry RA. A framework for oligonucleotide microarray preprocessing. Bioinformatics. 2010;26(19):2363–7.
Sean D, Meltzer PS. GEOquery: a bridge between the gene expression omnibus (GEO) and BioConductor. Bioinformatics. 2007;23(14):1846–7.
Ritchie ME, Phipson B, Wu D, Hu YF, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
Goeppert B, Stichel D, Toth R, Fritzsche S, Loeffler MA, Schlitter AM, et al. Integrative analysis reveals early and distinct genetic and epigenetic changes in intraductal papillary and tubulopapillary cholangiocarcinogenesis. Gut. 2022;71(2):391–401.
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, et al. Minfi: a flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics. 2014;30(10):1363–9.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9:1–13.
Hanzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinform. 2013;14:7.
Chun YS, Pawlik TM, Vauthey JN. 8th edition of the AJCC cancer staging manual: pancreas and hepatobiliary cancers. Ann Surg Oncol. 2018;25(4):845–7.
Wittekind C. Pitfalls in the classification of liver tumors. Pathologe. 2006;27(4):289–93.
Pu W, Wang C, Chen S, Zhao D, Zhou Y, Ma Y, et al. Targeted bisulfite sequencing identified a panel of DNA methylation-based biomarkers for esophageal squamous cell carcinoma (ESCC). Clin Epigenetics. 2017;9:129.
Marsico A, Huska MR, Lasserre J, Hu HY, Vucicevic D, Musahl A, et al. PROmiRNA: a new miRNA promoter recognition method uncovers the complex regulation of intronic miRNAs. Genome Biol. 2013;14(8):1–23.
Li LC, Dahiya R. MethPrimer: designing primers for methylation PCRs. Bioinformatics. 2002;18(11):1427–31.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.
Li YX, Li W. BSMAP: whole genome bisulfite sequence MAPping program. BMC Bioinform. 2009;10:1–9.
We thank Dr. Zhangfu Yang (Fudan University, CN) for kindly providing HuCCT1 and HCCC9810 cell lines, Dr. Xiaode Feng (Zhejiang University, CN) for kindly providing CCLP1 and HIBEpiC cell lines, and the epigenetic research team working at Yunbios (Shanghai, CN) for their help on targeted bisulfite sequencing data analyses, respectively.
This study was supported by the National Key Research and Development Program of China (2019YFC1316000) and National Natural Science Foundation of China (82188102).
Ethics approval and consent to participate
The experimental protocol (5826Z/20, 0531B/22) for using clinical liver tissue specimens was reviewed and approved by the Institutional Ethics Committee of the First Affiliated Hospital of Zhejiang University (Hangzhou, CN). The experimental protocol (1615/22) for all animal experiments was reviewed and approved by the Institutional Animal Care & Use Committee of the First Affiliated Hospital, Zhejiang University School of Medicine (Hangzhou, CN). Written informed consents were collected from each patient enrolled in this study.
Consent for publication
All authors agreed to publish this manuscript.
We have no conflict of interest to declare in this study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original online version of this article was revised: The corrected supplementary material for Figure S4 has been uploaded.
COL12A1 mRNA expression is upregulated in clinical iCCA. (A) Boxplot to show COL12A1 mRNA expression levels in clinical iCCA, CHC, FNH and adjacent nontumor tissues samples from GSE32879 dataset. ns, p ≥ 0.05, **** p < 0.00001, using two-tailed unpaired Student's t test. Boxplot to show COL12A1 mRNA expression levels in clinical iCCA and HCC tissue samples form TCGA, GSE76297 or our in-house dataset, respectively. * p < 0.05; **** p < 0.00001, using two-tailed paired/unpaired Student's t test accordingly. Gallbladder or extrahepatic carcinoma tissue samples (n = 4) in TCGA cholangiocarcinoma dataset were excluded from our analysis. (B) Boxplot to show COL12A1 mRNA expression levels in clinical HCC and nontumor liver tissue samples from TCGA, E-MTAB-4171, GSE76297 or OEP000321 dataset, respectively. ns, p ≥ 0.05, using two-tailed paired/unpaired Student's t test accordingly. (C) COL12A1 protein levels in whole cell lysate of human HCC cell lines (including Hep3B, Bel-7402, HCC-LM3, and Huh7) and normal liver cell line HL7702 were evaluated by immunoblotting. Bar plots showing the relative expression level of COL12A1 protein in the indicated cell lines. Experiments were in triplicates. ****p < 0.0001, based on one-way analysis of variance (ANOVA). iCCA, intrahepatic cholangiocarcinoma; HCC, hepatocellular carcinoma; CHC, combined hepatocellular cholangiocarcinoma; FNH, focal nodular hyperplasia; TCGA, The Cancer Genome Atlas.
Forest plot to show univariate survival analysis results of iCCA patients (n = 60). iCCA, intrahepatic cholangiocarcinoma.
The overall methylation profile of COL12A1 gene in iCCA and adjacent nontumor tissue samples from TCGA (A) and GSE156299 (B) dataset, respectively. TCGA, The Cancer Genome Atlas; iCCA, intrahepatic cholangiocarcinoma.
Deciphering DNA methylation pattern of COL12A1 promoter region in iCCA and adjacent nontumor tissue samples. (A) Boxplots to show the methylation degree of CpG islands at COL12A1 enhancer region in iCCA and paired adjacent nontumor tissues samples from TCGA dataset. ns, p ≥ 0.05, using two-tailed paired Student's t test. (B) Boxplots to show the methylation degree of CpG islands at COL12A1 enhancer region in iCCA, intraductal tubulopapillary neoplasm of biliary duct (ITPN), intraductal papillary neoplasm of biliary duct (IPNB), and adjacent nontumor tissues samples from GSE156299 dataset. ns, p ≥ 0.05; *p < 0.05; **p < 0.01, using two-tailed unpaired Wilcox test. TCGA, The Cancer Genome Atlas; iCCA, intrahepatic cholangiocarcinoma.
Candidate miRNAs (miR-15a-5p, -497-5p, -15b-5p, -6838-5p, -16-5p, -195-5p and miR-424-5p) binding to wild-type COL12A1 3'UTR regions were predicted by TargetScanHuman database (version 8.0).
The expression patterns of the candidate miRNAs in clinical iCCA and adjacent nontumor liver tissue samples from TCGA dataset. (A) Boxplots showing that miRNAs (miR-424-5p, -424-3p, -497-5p, -15b-5p, -15a-5p, -16-5p, -195-5p and miR-16-5p) in iCCA and nontumor liver tissue samples from TCGA dataset. ns, p ≥ 0.05. **p < 0.01; ***p < 0.001, using two-tailed unpaired Student’s t test. (B-C) Heatmap shows the association of COL12A1 with miRNAs (miR-424-5p, -424-3p, -497-5p, -15b-5p, -15a-5p, -16-5p, -195-5p and miR-16-5p) in iCCA and nontumor samples from TCGA dataset. Numbers in heatmap denote Pearson coefficient (B) and statistical significance (C). Gallbladder or extrahepatic carcinoma tissue samples (n = 4) in TCGA cholangiocarcinoma dataset were excluded from our analysis. (D) MiR-424-5p in human iCCA (HuCCT1 and CCLP1 cell line) and normal human intrahepatic biliary cell (HIBEpiC) lines was determined by RT-PCR, respectively. RT-PCR, real-time PCR; TCGA, The Cancer Genome Atlas; iCCA, intrahepatic cholangiocarcinoma.
Correlation of COL12A1 expression with clinicopathological features of 60 iCCA patients.
The clinicopathological characteristics of 10 iCCA/HCC patients for COL12A1 expression profile analysis.
The designed primers of miR-424 promoter for targeted bisulfite sequencing.
About this article
Cite this article
Tang, Z., Yang, Y., Zhang, Q. et al. Epigenetic dysregulation-mediated COL12A1 upregulation predicts worse outcome in intrahepatic cholangiocarcinoma patients. Clin Epigenet 15, 13 (2023). https://doi.org/10.1186/s13148-022-01413-5