Skip to main content

A metabolism-related 4-lncRNA prognostic signature and corresponding mechanisms in intrahepatic cholangiocarcinoma



Long non-coding RNA (lncRNA) plays a critical role in the malignant progression of intrahepatic cholangiocarcinoma (iCCA). This study aimed to establish a 4-lncRNA prognostic signature and explore corresponding potential mechanisms in patients with iCCA.


The original lncRNA-seq and clinical data were collected from the TCGA and GEO databases. Overlapping and differentially expressed lncRNAs (DE-lncRNAs) were further identified from transcriptome data. Univariate regression analysis was performed to screen survival-related DE-lncRNAs, which were further selected to develop an optimal signature to predict prognosis using multivariate regression analysis. The Kaplan-Meier survival curve visualized the discrimination of the signature on overall survival (OS). The area under the curve (AUC) and C-index were used to verify the predictive accuracy of the signature. Combined with clinical data, multivariate survival analysis was used to reveal the independent predictive capability of the signature. In addition, a prognostic nomogram was constructed. Finally, the common target genes of 4 lncRNAs were predicted by the co-expression method, and the corresponding functions were annotated by GO and KEGG enrichment analysis. Gene set enrichment analysis (GSEA) was also performed to explore the potential mechanism of the signature. Quantitative real-time PCR was used to evaluated the expression of 4 lncRNAs in an independent cohort.


We identified and constructed a 4-lncRNA (AC138430.1, AGAP2-AS1, AP001783.1, and AP005233.2) prognostic signature using regression analysis, and it had the capability to independently predict prognosis. The AUCs were 0.952, 0.909, and 0.882 at 1, 2, and 3 years, respectively, and the C-index was 0.808, which showed good predictive capability. Subsequently, combined with clinical data, we constructed a nomogram with good clinical application. Finally, 252 target genes of all four lncRNAs were identified by the co-expression method, and functional enrichment analysis showed that the signature was strongly correlated with metabolism-related mechanisms in tumourigenesis. The same results were also validated via GSEA.


We demonstrated that a metabolism-related 4-lncRNA prognostic signature could be a novel biomarker and deeply explored the target genes and potential mechanism. This study will provide a promising therapeutic strategy for patients with intrahepatic cholangiocarcinoma.

Peer Review reports


Intrahepatic cholangiocarcinoma (iCCA) is a malignant hepatobiliary tumour originating from the intrahepatic bile ducts and accounts for 10–20% of all bile duct malignancies [1]. Over the past few decades, the incidence of iCCA has been reported to be steadily increasing in most of the world, and it frequently occurs in patients with underlying liver disease [1, 2]. Because of its high malignancy, hidden nature, and aggressive nature, iCCA lacks obvious symptoms in the early stages, which makes diagnosis and treatment difficult for clinicians [2, 3]. Although the prognosis of iCCA patients has improved in recent years with the enhancement of surgical techniques, the application of adjuvant therapies, and the gradual promotion of emerging immunotherapies and targeted therapies for specific mutation sites [4, 5], the long-term prognosis of iCCA patients still needs to be strengthened. Therefore, there is an urgent need to explore new prognostic indicators and reveal potential mechanisms to better understand disease progression.

lncRNAs are greater than 200 bp in length and do not have a protein-coding function [6]. With the increasing exploration and development of cancerous molecular pathology, numerous studies have proven that lncRNAs are functional RNAs that can play a corresponding role similar to that of oncogenes or tumour suppressor genes [7] and are not, as previously thought, “transcriptional noise” (non-functional RNA). Recently, researchers have become very interested in their role in tumour invasion, metastasis, and signalling pathways, and previous studies have demonstrated the involvement of lncRNAs as key molecules in the malignant biology of various cancers [8,9,10]. At present, some studies have described the relevant roles of fractional lncRNAs in cholangiocarcinoma [11,12,13], but the role of more lncRNAs in the malignant progression of intrahepatic cholangiocarcinoma still needs to be further identified and analysed.

In this study, we identified and constructed a metabolism-related lncRNA prognostic signature and incorporated it into the nomogram for OS prediction. The relevant target genes of the signature and their potential mechanisms were also explored. The findings may provide new prognostic biomarkers that reveal a novel perspective for the individualized treatment of patients with iCCA.

Materials and methods

Data obtaining and pre-processing

The CHOL RNA-seq dataset and clinical data were downloaded from the TCGA database ( Another dataset, GSE107943, was collected from the GEO database ( Then, the lncRNA matrix was annotated and extracted from the RNA-seq dataset using the “Perl” language and “GenomicTools” package.

Differential expression analysis

The differentially expressed genes and DE-lncRNAs were analysed by using the “edgeR” package with the threshold FDR < 0.05 coupled with |log2foldchange(FC)| > 2 in two datasets [14]. Overlapping DE-lncRNAs were extracted and exhibited using the “Venn” package [15], and the heatmap and volcano map were plotted to illustrate the up- (red) and downregulated (green) DE-lncRNAs using the “ggplot2” package.

Construction and validation of the prognostic signature

The survival-related DE-lncRNAs were screened using univariate regression analysis and were further incorporated into the multivariable regression analysis to obtain independent prognostic DE-lncRNAs for constructing a prognostic signature. p < 0.05 was deemed the threshold. The coefficients and expression levels of DE-lncRNAs were obtained to calculate the risk score. Based on the median risk score, we divided all patients into high- and low-risk groups. The Kaplan-Meier (KM) survival curve described the predictive power of indicators using the “survival” and “survminer” packages. The AUC of the receiver operating characteristic (ROC) curve and C-index were applied to verify the predictive accuracy using the “timeROC” and “rms” packages. The clinical value of the signature was analysed. To reveal the potential mechanism of the prognostic signature, we performed gene set enrichment analysis (GSEA).

Validation of lncRNAs using quantitative real-time PCR

Tumour tissues with pathological diagnosis of intrahepatic cholangiocarcinoma and normal tissues were prospective collected from the Chinese PLA general Hospital. qRT-PCR was performed to validate the gene expression changes of all four lncRNAs. TRIzol reagent (Ambion) was used to extracted total RNA; NanoPhotometer® C40 Touch (IMPLEN) was used to assess the RNA purity based on the ratio of OD260/280 and 260/230; Eppendorf Mastercycler® was used to perform reverse transcription of qualified RNA to single-stranded complementary DNA according to the manufacturer’s instructions; QuanStudio™ 5 Real-Time PCR instrument was used to implement real-time quantification; β-actin was used as internal reference; finally, we recorded cycle threshold (Ct) and calculated the relative expression of four lncRNAs using the 2−ΔΔCt method. Primers sequences of four lncRNAs and β-actin were shown in Supplementary Table S1.

Co-expression prediction of target genes and functional annotation

We calculated the Pearson’s correlation coefficients between four DE-lncRNA expression profiles and protein-coding genes (PCGs) to determine the co-expression relationship, where PCGs with |Pearson’s correlation coefficient| > 0.5 were considered lncRNA-associated PCGs. Then, we took the intersection of PCGs and differentially expressed genes, and the overlapping PCGs were further functionally annotated using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis in the “clusterProfiler” package [16].

Development and assessment of the nomogram

To reveal the independent prognostic value of the signature, univariate and multivariate regression analyses were implemented on the signature and clinical data involving age, sex, race, T stage, N stage, M stage, and AJCC stage. Based on all independently predictive variables, we constructed a prediction nomogram using the “rms” and “foreign” packages. The predictivity of the nomogram was validated using AUC and calibration curves.

Statistical analysis

R version 4.0.2 software and its resource packages were used for statistical analysis and related visualization graphics. Pearson’s correlation coefficients were calculated between the expression profiles of four DE-lncRNAs and PCGs to determine the co-expression relationship among them, and |coef| > 0.5 was considered to have relevance. Principal component analysis (PCA) was implemented to efficiently downscale high-dimensional lncRNA sequencing data to improve the ability for data identification and was implemented to analyse the capability to distinguish patients with high or low risk based on a risk model, differentially expressed lncRNAs, and whole lncRNA sequences. The function of the prognostic signature was annotated by GSEA. Statistical significance of all statistical tests implemented in this study was determined as two-sided p < 0.05.


Differential expression analysis

Meeting the screening criteria, a total of 1851 DE-lncRNAs were identified from the TCGA dataset, of which 1334 were upregulated and 518 were downregulated (Fig. 1a-b), while a total of 362 DE-lncRNAs were identified from the GEO dataset, of which 197 were upregulated and 165 were downregulated (Fig. 1c-d). Finally, 214 overlapping DE-lncRNAs were extracted from the intersection of DE-lncRNAs (Fig. 1e).

Fig. 1

Differentially expressed analysis. a-b Heatmap and volcano plot in the TCGA dataset. c-d Heatmap and volcano plot in the GEO dataset. e Venn diagram of the differentially expressed lncRNA intersection. f Survival-related lncRNAs were obtained using univariate regression analysis

Construction and validation of the prognostic signature

Univariate Cox regression analysis showed that six lncRNAs were significantly related to OS in all DE-lncRNAs (p < 0.05, Fig. 1f). Subsequently, an optimal prognostic signature involving 4 lncRNAs (AC138430.1, AGAP2-AS1, AP001783.1, and AP005233.2) was constructed using multivariate Cox regression analysis, of which the lncRNAs AGAP2-AS1 and AP005233.2 were upregulated, whereas the lncRNAs AC138430.1 and AP001783.1 were downregulated in tumour tissues. The risk coefficients suggested that all 4 lncRNAs were risk factors for iCCA with |coef| > 0 (Table 1).

Table 1 The Four lncRNAs were identified using the multivariate analysis

We extracted the coefficients and expression of 4 lncRNAs and calculated risk scores with the following equation: risk score = (0.459) *AC138430.1 + (0.004) *AGAP2-AS1 + (0.096) *AP001783.1 + (0.010) *AP005233.2. Depending on the median risk score, we divided all patients into high- and low-risk groups (n = 16 and 17, respectively). Corresponding survival condition plots showed that the high-risk group had higher mortality than the low-risk group (Fig. 2a-b), and the related gene expression patterns of four lncRNAs in the two risk groups are shown in Fig. 2c. The Kaplan-Meier survival curve revealed that patients with high risk scores had a significantly shorter OS (Fig. 2d). The AUCs were 0.952, 0.909, and 0.882 at 1, 2, and 3 years, respectively (Fig. 2e), and the C-index was 0.808, which demonstrated that this signature performed well as a predictor of prognosis. The PCA results showed a clear distribution of high and low risk on both sides based on the risk model, which could well classify patients into high- and low-risk groups compared to the other two approaches (Fig. 3a-c). qRT-PCR analysis showed the similar trends of four lncRNAs expression compared to differentially expression analysis (Supplementary Figure S1), which revealed all four lncRNAs were differentially expressed between the tumour and normal tissue in an independent cohort and may involve in the tumorigenesis of iCCA.

Fig. 2

a-c Survival condition plots and heatmap of 4 IRGs. d Kaplan-Meier survival curve. e Time-dependent ROC curves for predicting OS at 1, 2, and 3 years

Fig. 3

Principal components analysis. a-c Distribution of high and low risk based on all lncRNAs, differentially expressed lncRNAs, and the risk model

Co-expression prediction and functional enrichment analysis

We used the co-expression method to predict the highly correlated PCGs of four DE-lncRNAs. A total of 252 differentially expressed PCGs were identified (Fig. 4a), which had a high correlation with all four lncRNAs (coef> 0.5, p < 0.05). Alluvial plot of PCGs illustrating associations with lncRNAs and survival state (Fig. 4b). The list of differentially expressed genes and PCGs is shown in Supplementary Tables S23. The GO and KEGG analyses showed that the function of PCGs was highly enriched in metabolism-related processes (Fig. 4b-c). The results showed that the four DE-lncRNAs were extensively involved in the metabolism-related processes of tumourigenesis.

Fig. 4

a Venn diagram of the intersection of co-expressed genes. b Association of mRNAs with lncRNAs and regulated types. c GO enrichment analysis. d KEGG pathway enrichment analysis

Gene set enrichment analysis

Based on these results, the prognostic signature has been shown to have good predictive power for OS. To further explore the potential mechanisms of the prognostic signature, we performed GSEA in the high-risk group. As expected, the enriched results showed that the signature had a strong correlation with metabolism in tumourigenesis (p < 0.05), which was similar to the pathways of co-expression target genes. The key enrichment pathways of the prognostic signature were “nicotinate and nicotinamide metabolism”, “butanoate metabolism”, “regulation of autophagy”, and “pantothenate and COA biosynthesis” (Fig. 5a). The detailed GSEA results are shown in Supplementary Figure S2.

Fig. 5

a Four crucially enriched pathways were displayed by using gene set enrichment analysis. b-c Forest plot of univariate and multivariate regression analyses

Independent prognosis analysis

The clinical data containing 33 patients with follow-up were collected from the TCGA database, as shown in Table 2. This prognostic signature was verified to be an independent prognostic indicator using Cox regression analyses, as shown in Fig. 5b-c, and the risk score and N stage were significantly correlated with the OS of patients in both univariate and multivariate regression analyses (p < 0.05). The results demonstrated that they had independent predictability.

Table 2 Clinicopathologic characteristics of patients with intrahepatic cholangiocarcinoma

Construction and verification of the nomogram

Next, the above independent predictors were integrated into a nomogram for predicting 1-, 2-, and 3-year survival rates (Fig. 6a). Based on the nomogram score, the 1-, 2-, and 3-year survival rates of patients could be well predicted according to their nomogram scores. The AUCs at 1, 2, and 3 years were 0.912, 0.920, and 0.923, respectively, in the nomogram prediction (Fig. 6b). The multivariable ROC curve demonstrated that the risk score had the best predictive capability (Fig. 6c). In addition, the calibration curves showed good agreement between the predicted survival rate and the actual survival rate at 1, 2, and 3 years (Fig. 6d-f). All the results suggested that the nomogram had good predictive accuracy and clinical application prospects.

Fig. 6

a Nomogram for predicting OS at 1, 2, and 3 years. b Time-dependent ROC curves of nomogram prediction at 1, 3, and 5 years. c The AUC for the risk score and clinical features. d-f Calibration curves showing the probability of 1-, 2-, and 3-year OS between the nomogram prediction and practical observation in the TCGA dataset


Like most malignant tumours, the diagnosis and prognosis of patients with intrahepatic cholangiocarcinoma remain relatively poor, even though various treatments have been gradually improved and optimally applied [1]. In recent years, exploring tumour biomarkers has become an increasingly popular field in tumour therapy, which has greatly aroused the interest of researchers [5]. Meanwhile, with the extensive development of sequencing technology and the deepening of transcriptomics, lncRNAs have gradually entered the perspective of researchers [6, 17, 18]. Recently, many studies have demonstrated that lncRNAs are not “transcriptional noise”, as originally thought, but that they play roles as oncogenes or tumour suppressor genes [7] and act as competing endogenous RNAs (ceRNAs) to regulate key target genes in tumour development [19, 20]. Previous studies have identified a variety of lncRNAs involved in tumour metabolism, autophagy, and immune-related processes [21,22,23,24], while many researchers have developed various predictive models based on prognosis-related lncRNAs for application [25,26,27]. We are very interested in pathways related to tumour metabolism and would like to develop a predictive model for intrahepatic cholangiocarcinoma. Zhang et al.’s study identified a ceRNA regulation network containing 25 lncRNAs that were highly involved in bone metabolism-related biological processes and corresponding pathways [28]. Huang et al. also identified 9 lncRNAs that may be mainly involved in metabolism-related pathways of lung squamous cell carcinoma [29]. However, fewer lncRNAs have been developed for metabolism-related processes in intrahepatic cholangiocarcinoma, and the identification of more prognostic lncRNAs and related models is urgently needed to provide novel therapeutic targets in intrahepatic cholangiocarcinoma.

In this study, four survival-related and metabolism-related lncRNAs (AC138430.1, AGAP2-AS1, AP001783.1, and AP005233.2) were identified for prognosis prediction in patients with intrahepatic cholangiocarcinoma based on transcriptome data from public databases, and were verified in an independent cohort via quantitative real-time PCR. It shows these four lncRNAs may highly involve in the tumorigenesis of iCCA.

Recently, lncRNA AGAP2-AS1 (AGAP2 antisense RNA 1) has been reported to play crucial roles in some tumours, such as colorectal cancer [30], pancreatic cancer [31], and breast cancer [32]. It is extensively involved in important processes of tumourigenesis and development and participates in the expression process of target genes as a ceRNA. However, to our knowledge, high-quality studies to explore the potential mechanism of AGAP2-AS1 in intrahepatic cholangiocarcinoma are lacking. Hence, studying the correlation of AGAP2-AS1 with tumourigenesis in intrahepatic cholangiocarcinoma is significantly meaningful. The lncRNA AP005233.2 has only been identified as a key lncRNA in pathway crosstalk of lung adenocarcinoma [33]. Consequently, the mechanism of AP005233.2 in intrahepatic cholangiocarcinoma or other tumours also needs further investigation. Notably, no study has explored the mechanisms of lncRNA AC138430.1 and AP001783.1 in tumours. Therefore, their pivotal roles in intrahepatic cholangiocarcinoma are also worth revealing in future studies for the first time. In addition, according to whole genome sequencing, we identified the potential PCGs of the four lncRNAs using the co-expression method, and these expression states of PCGs may be upregulated or downregulated via the lncRNA-mRNA regulation network. Next, we performed GO and KEGG analyses to reveal a potential mechanism. The results showed that these four lncRNAs were extensively involved in metabolic processes. Previous studies have demonstrated that metabolism is highly related to tumourigenesis and metastasis of intrahepatic cholangiocarcinoma [34, 35], such as the fatty acid synthesis pathway [36]. Zhang et al. revealed that fatty acid-related metabolism plays a crucial role in the tumourigenesis of intrahepatic cholangiocarcinoma [37]. Based on four metabolism-related lncRNAs, we developed a prognostic signature with independent predictive capability using Cox regression analysis. The AUCs were 0.952, 0.909, and 0.882 at 1, 2, and 3 years, respectively, and the C-index was 0.808, which showed that our prognostic signature had good predictive accuracy. The multi-ROC curve showed that the signature had better predictive capability than other clinicopathologic characteristics, such as AJCC stage, so we could predict the patient’s prognosis using this signature. The PCA results showed that the signature including four lncRNAs could well distinguish high- and low-risk patients with intrahepatic cholangiocarcinoma to guide clinical grouping. Meanwhile, as we expected, GSEA also revealed that the signature containing four lncRNAs was highly related to metabolic processes, such as the nicotinate and nicotinamide metabolism pathway; therefore, these lncRNAs would have key effectiveness in tumour metabolism. Finally, we combined the two independent prognostic factors (N stage and risk score) to construct a nomogram for clinical application. The corresponding results also indicated that this nomogram had good discrimination capability. The calibration curves showed that there were good discriminative and calibration capabilities for this nomogram.

In addition, some limitations must also be noted in this study. First, all the data in our study were obtained from public databases to perform the retrospective analysis, and selection bias is inevitable. Secondly, our study included the number of prospective samples may be insufficient, thus, further large prospective cohort studies must be implemented to confirm the effectiveness of our prognostic signature. Then, due to the different sequencing methods in the TCGA and GEO databases, the standardization process was difficult to unify in our study, so we only completed the differential expression analysis for screening the overlapping DE-lncRNAs. But notably, our study validated the differential expression of these four DE-lncRNAs in an independent cohort using qRT-PCR, so them could be considered as a valuable predictive factor.


We identified and validated a 4-lncRNA (AC138430.1, AGAP2-AS1, AP001783.1, and AP005233.2) prognostic signature that had a good predictive capability for prognosis in patients with intrahepatic cholangiocarcinoma. Meanwhile, the signature was highly related to metabolic pathways, which helps researchers deeply understand the correlation of metabolism with tumourigenesis. Hopefully, it will provide a new perspective for exploring biomarkers in the tumour metabolic microenvironment of intrahepatic cholangiocarcinoma.

Availability of data and materials

The datasets downloaded for supporting the results of this article are publicly available at the TCGA ( and GEO (



American Joint Committee on Cancer


The area under curve


Overall survival


Gene Expression Omnibus


Gene Ontology


Intrahepatic cholangiocarcinoma


Kyoto Encyclopedia of Genes and Genomes


Receiver operating characteristic


The Cancer Genome Atlas


Protein-coding genes


Principal component analysis


Gene set enrichment analysis


  1. 1.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Zhang H, Yang T, Wu M, Shen F. Intrahepatic cholangiocarcinoma: epidemiology, risk factors, diagnosis and surgical management. Cancer Lett. 2016;379(2):198–205.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Razumilava N, Gores GJ. Cholangiocarcinoma. Lancet. 2014;383(9935):2168–79.

    Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Kelley RK, Bridgewater J, Gores GJ, Zhu AX. Systemic therapies for intrahepatic cholangiocarcinoma. J Hepatol. 2020;72(2):353–63.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Lamarca A, Barriuso J, McNamara MG, Valle JW. Molecular targeted therapies: ready for “prime time” in biliary tract cancer. J Hepatol. 2020;73(1):170–85.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Chandra Gupta S, Nandan Tripathi Y. Potential of long non-coding RNAs in cancer patients: From biomarkers to therapeutic targets. Int J Cancer. 2017;140(9):1955–67.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Wu P, Mo Y, Peng M, Tang T, Zhong Y, Deng X, et al. Emerging role of tumor-related functional peptides encoded by lncRNA and circRNA. Mol Cancer. 2020;19(1):22.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Goodall GJ, Wickramasinghe VO. RNA in cancer. Nat Rev Cancer. 2021;21(1):22-36.

  9. 9.

    Chi Y, et al. Long Non-Coding RNA in the Pathogenesis of Cancers. Cells. 2019;8(9):1015.

  10. 10.

    Fang Y, Fullwood MJ. Roles, functions, and mechanisms of long non-coding RNAs in cancer. Genomics Proteomics Bioinformatics. 2016;14(1):42–54.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Shen S, Wang J, Zheng B, Tao Y, Li M, Wang Y, et al. LINC01714 enhances gemcitabine sensitivity by modulating FOXO3 phosphorylation in cholangiocarcinoma. Mol Ther Nucleic Acids. 2020;19:446–57.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Wan T, Wang H, Gou M, Si H, Wang Z, Yan H, et al. LncRNA HEIH promotes cell proliferation, migration and invasion in cholangiocarcinoma by modulating miR-98-5p/HECTD4. Biomed Pharmacother. 2020;125:109916.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Li Z, Jiang X, Huang L, Li J, Ji D, Xu Y, et al. Up-regulation of ZFAS1 indicates dismal prognosis for cholangiocarcinoma and promotes proliferation and metastasis by modulating USF1 via miR-296-5p. J Cell Mol Med. 2019;23(12):8258–68.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Chen H, Boutros PC. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinformatics. 2011;12(1):35.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16(5):284–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Wardell CP, Fujita M, Yamada T, Simbolo M, Fassan M, Karlic R, et al. Genomic characterization of biliary tract cancers identifies driver genes and predisposing mutations. J Hepatol. 2018;68(5):959–69.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Xue R, et al. Genomic and transcriptomic profiling of combined hepatocellular and intrahepatic cholangiocarcinoma reveals distinct molecular subtypes. Cancer Cell. 2019;35(6):932–947.e8.

    CAS  Article  Google Scholar 

  19. 19.

    Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014;505(7483):344–52.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Xu W, Yu S, Xiong J, Long J, Zheng Y, Sang X. CeRNA regulatory network-based analysis to study the roles of noncoding RNAs in the pathogenesis of intrahepatic cholangiocellular carcinoma. Aging (Albany NY). 2020;12(2):1047–86.

    CAS  Article  Google Scholar 

  21. 21.

    Huang F, Chen W, Peng J, Li Y, Zhuang Y, Zhu Z, et al. LncRNA PVT1 triggers Cyto-protective autophagy and promotes pancreatic ductal adenocarcinoma development via the miR-20a-5p/ULK1 Axis. Mol Cancer. 2018;17(1):98.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Li Y, Jiang T, Zhou W, Li J, Li X, Wang Q, et al. Pan-cancer characterization of immune-related lncRNAs identifies potential oncogenic biomarkers. Nat Commun. 2020;11(1):1000.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Sun J, et al. Identification of tumor immune infiltration-associated lncRNAs for improving prognosis and immunotherapy response of patients with non-small cell lung cancer. J Immunother Cancer. 2020;8(1):e000110.

  24. 24.

    Meng X, Feng C, Fang E, Feng J, Zhao X. Combined analysis of RNA-sequence and microarray data reveals effective metabolism-based prognostic signature for neuroblastoma. J Cell Mol Med. 2020;24(18):10367–81.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Kong W, Wang X, Zuo X, Mao Z, Cheng Y, Chen W. Development and validation of an immune-related lncRNA signature for predicting the prognosis of hepatocellular carcinoma. Front Genet. 2020;11:1037.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Li W, Chen QF, Huang T, Wu P, Shen L, Huang ZL. Identification and validation of a prognostic lncRNA signature for hepatocellular carcinoma. Front Oncol. 2020;10:780.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Liu Y, Wang L, Liu H, Li C, He J. The prognostic significance of metabolic syndrome and a related six-lncRNA signature in esophageal squamous cell carcinoma. Front Oncol. 2020;10:61.

    Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Zhang M, Cheng L, Zhang Y. Characterization of dysregulated lncRNA-associated ceRNA network reveals novel lncRNAs with ceRNA activity as epigenetic diagnostic biomarkers for osteoporosis risk. Front Cell Dev Biol. 2020;8:184.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Huang G, Huang Q, Xie Z, Zhou H, Cao J, Shi L, et al. A nine-long non-coding RNA signature for prognosis prediction of patients with lung squamous cell carcinoma. Cancer Biomark. 2019;26(3):239–47.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Hong S, Yan Z, Song YM, Bi MM, Li S. LncRNA AGAP2-AS1 augments cell viability and mobility, and confers gemcitabine resistance by inhibiting miR-497 in colorectal cancer. Aging (Albany NY). 2020;12(6):5183–94.

    CAS  Article  Google Scholar 

  31. 31.

    Hui B, Ji H, Xu Y, Wang J, Ma Z, Zhang C, et al. RREB1-induced upregulation of the lncRNA AGAP2-AS1 regulates the proliferation and migration of pancreatic cancer partly through suppressing ANKRD1 and ANGPTL4. Cell Death Dis. 2019;10(3):207.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Han J, et al. MSC-induced lncRNA AGAP2-AS1 promotes stemness and trastuzumab resistance through regulating CPT1 expression and fatty acid oxidation in breast cancer. Oncogene. 2021;40(4):833-47.

  33. 33.

    Qi G, Kong W, Mou X, Wang S. A new method for excavating feature lncRNA in lung adenocarcinoma based on pathway crosstalk analysis. J Cell Biochem. 2019;120(6):9034–46.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Jia X, Lu S, Zeng Z, Liu Q, Dong Z, Chen Y, et al. Characterization of gut microbiota, bile acid metabolism, and cytokines in intrahepatic cholangiocarcinoma. Hepatology. 2020;71(3):893–906.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Chaisaingmongkol J, et al. Common molecular subtypes among asian hepatocellular carcinoma and cholangiocarcinoma. Cancer Cell. 2017;32(1):57–70.e3.

    CAS  Article  Google Scholar 

  36. 36.

    Nie J, Zhang J, Wang L, Lu L, Yuan Q, An F, et al. Adipocytes promote cholangiocarcinoma metastasis through fatty acid binding protein 4. J Exp Clin Cancer Res. 2017;36(1):183.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Zhang B, Zhou BH, Xiao M, Li H, Guo L, Wang MX, et al. KDM5C represses FASN-mediated lipid metabolism to exert tumor suppressor activity in intrahepatic cholangiocarcinoma. Front Oncol. 2020;10:1025.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We acknowledge TCGA and GEO databases for providing their platforms and contributors for uploading their meaningful datasets.


This study was supported by funds from National Key Sci-Tech Special Project of China (2018ZX10302204003) and Medical Data Research and Development Project of PLA General Hospital (2018MBD-019).

Author information




WBZ, ZZW and FW contributed equally to this work, and were considered as the co-first authors. WBZ, ZZW and FW participate in design and conception of this study. FW and LCL downloaded the data form corresponding databases. WBZ performed the bioinformatics analysis. All the authors participate in writing the manuscript, while RL and MGH revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Rong Liu or Minggen Hu.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the ethics committee of the Chinese PLA general Hospital. All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. Written informed consent was obtained from all patients.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Fig. S1.

A-D The results of quantitative real-time PCR showed that relative expression level of 4 key lncRNAs between tumor and normal tissue. * P < 0.05; ** P < 0.01; *** P < 0.001. Supplementary Fig. S2 A-D Four detailed gene set enrichment pathways in high-risk group. Table S1. Primers used for quantitative real-time PCR. Table S2. Correlation analasisi of protein coding gene with each lncRNA. Table S3. Differentially expressed protein coding gene.

Rights and permissions

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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zou, W., Wang, Z., Wang, F. et al. A metabolism-related 4-lncRNA prognostic signature and corresponding mechanisms in intrahepatic cholangiocarcinoma. BMC Cancer 21, 608 (2021).

Download citation


  • Long non-coding RNA
  • Intrahepatic cholangiocarcinoma
  • Overall survival
  • Signature
  • Nomogram