A metabolism-related 4-lncRNA prognostic signature and corresponding mechanisms in intrahepatic cholangiocarcinoma
BMC Cancer volume 21, Article number: 608 (2021)
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.
Intrahepatic cholangiocarcinoma (iCCA) is a malignant hepatobiliary tumour originating from the intrahepatic bile ducts and accounts for 10–20% of all bile duct malignancies . 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 . 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  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 (https://portal.gdc.cancer.gov/). Another dataset, GSE107943, was collected from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). 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 . Overlapping DE-lncRNAs were extracted and exhibited using the “Venn” package , 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 .
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.
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).
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).
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.
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 S2–3. 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.
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.
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.
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.
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 . In recent years, exploring tumour biomarkers has become an increasingly popular field in tumour therapy, which has greatly aroused the interest of researchers . 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  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 . Huang et al. also identified 9 lncRNAs that may be mainly involved in metabolism-related pathways of lung squamous cell carcinoma . 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 , pancreatic cancer , and breast cancer . 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 . 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 . Zhang et al. revealed that fatty acid-related metabolism plays a crucial role in the tumourigenesis of intrahepatic cholangiocarcinoma . 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.
American Joint Committee on Cancer
The area under curve
Gene Expression Omnibus
Kyoto Encyclopedia of Genes and Genomes
Receiver operating characteristic
The Cancer Genome Atlas
Principal component analysis
Gene set enrichment analysis
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. https://doi.org/10.1038/s41575-020-0310-z.
Zhang H, Yang T, Wu M, Shen F. Intrahepatic cholangiocarcinoma: epidemiology, risk factors, diagnosis and surgical management. Cancer Lett. 2016;379(2):198–205. https://doi.org/10.1016/j.canlet.2015.09.008.
Razumilava N, Gores GJ. Cholangiocarcinoma. Lancet. 2014;383(9935):2168–79. https://doi.org/10.1016/S0140-6736(13)61903-0.
Kelley RK, Bridgewater J, Gores GJ, Zhu AX. Systemic therapies for intrahepatic cholangiocarcinoma. J Hepatol. 2020;72(2):353–63. https://doi.org/10.1016/j.jhep.2019.10.009.
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. https://doi.org/10.1016/j.jhep.2020.03.007.
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. https://doi.org/10.1002/ijc.30546.
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. https://doi.org/10.1186/s12943-020-1147-3.
Goodall GJ, Wickramasinghe VO. RNA in cancer. Nat Rev Cancer. 2021;21(1):22-36. https://doi.org/10.1038/s41568-020-00306-0.
Chi Y, et al. Long Non-Coding RNA in the Pathogenesis of Cancers. Cells. 2019;8(9):1015. https://doi.org/10.3390/cells8091015.
Fang Y, Fullwood MJ. Roles, functions, and mechanisms of long non-coding RNAs in cancer. Genomics Proteomics Bioinformatics. 2016;14(1):42–54. https://doi.org/10.1016/j.gpb.2015.09.006.
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. https://doi.org/10.1016/j.omtn.2019.11.028.
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. https://doi.org/10.1016/j.biopha.2020.109916.
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. https://doi.org/10.1111/jcmm.14698.
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. https://doi.org/10.1093/bioinformatics/btp616.
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. https://doi.org/10.1186/1471-2105-12-35.
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. https://doi.org/10.1089/omi.2011.0118.
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. https://doi.org/10.1016/j.jhep.2018.01.009.
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.
Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014;505(7483):344–52. https://doi.org/10.1038/nature12986.
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. https://doi.org/10.18632/aging.102634.
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. https://doi.org/10.1186/s12943-018-0845-6.
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. https://doi.org/10.1038/s41467-020-14802-2.
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. https://doi.org/10.1136/jitc-2019-000110.
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. https://doi.org/10.1111/jcmm.15650.
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. https://doi.org/10.3389/fgene.2020.01037.
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. https://doi.org/10.3389/fonc.2020.00780.
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. https://doi.org/10.3389/fonc.2020.00061.
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. https://doi.org/10.3389/fcell.2020.00184.
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. https://doi.org/10.3233/CBM-182275.
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. https://doi.org/10.18632/aging.102940.
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. https://doi.org/10.1038/s41419-019-1384-9.
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. https://doi.org/10.1038/s41388-020-01574-8.
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. https://doi.org/10.1002/jcb.28177.
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. https://doi.org/10.1002/hep.30852.
Chaisaingmongkol J, et al. Common molecular subtypes among asian hepatocellular carcinoma and cholangiocarcinoma. Cancer Cell. 2017;32(1):57–70.e3.
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. https://doi.org/10.1186/s13046-017-0641-y.
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. https://doi.org/10.3389/fonc.2020.01025.
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).
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
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
About this article
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). https://doi.org/10.1186/s12885-021-08322-5
- Long non-coding RNA
- Intrahepatic cholangiocarcinoma
- Overall survival