Construction of a cancer-associated fibroblasts-related long non-coding RNA signature to predict prognosis and immune landscape in pancreatic adenocarcinoma.

Abstract

Background: Cancer-associated fibroblasts (CAFs) are an essential cell population in the pancreatic cancer tumor microenvironment and are extensively involved in drug resistance and immune evasion mechanisms. Long non-coding RNAs (lncRNAs) are involved in pancreatic cancer evolution and regulate the biological behavior mediated by CAFs. However, there is a lack of understanding of the prognostic signatures of CAFs-associated lncRNAs in pancreatic cancer patients. Methods: Transcriptomic and clinical data for pancreatic adenocarcinoma (PAAD) and the corresponding mutation data were obtained from The Cancer Genome Atlas database. lncRNAs associated with CAFs were obtained using co-expression analysis. lncRNAs were screened by Cox regression analysis using least absolute shrinkage and selection operator (LASSO) algorithm for constructing predictive signature. According to the prognostic model, PAAD patients were divided into high-risk and low-risk groups. Kaplan-Meier analysis was used for survival validation of the model in the training and validation groups. Clinicopathological parameter correlation analysis, univariate and multivariate Cox regression, time-dependent receiver operating characteristic (ROC) curves, and nomogram were performed to evaluate the model. The gene set variation analysis (GSVA) and gene ontology (GO) analyses were used to explore differences in the biological behavior of the risk groups. Furthermore, single-sample gene set enrichment analysis (ssGSEA), tumor mutation burden (TMB), ESTIMATE algorithm, and a series of immune correlation analyses were performed to investigate the relationship between predictive signature and the tumor immune microenvironment and screen for potential responders to immune checkpoint inhibitors. Finally, drug sensitivity analyses were used to explore potentially effective drugs in high- and low-risk groups. Results: The signature was constructed with seven CAFs-related lncRNAs (AP005233.2, AC090114.2, DCST1-AS1, AC092171.5, AC002401.4, AC025048.4, and CASC8) that independently predicted the prognosis of PAAD patients. Additionally, the high-risk group of the model had higher TMB levels than the low-risk group. Immune correlation analysis showed that most immune cells, including CD8+ T cells, were negatively correlated with the model risk scores. ssGSEA and ESTIMATE analyses further indicated that the low-risk group had a higher status of immune cell infiltration. Meanwhile, the mRNA of most immune checkpoint genes, including PD1 and CTLA4, were highly expressed in the low-risk group, suggesting that this population may be “hot immune tumors” and have a higher sensitivity to immune checkpoint inhibitors (ICIs). Finally, the predicted half-maximal inhibitory concentrations of some chemical and targeted drugs differ between high- and low-risk groups, providing a basis for treatment selection. Conclusion: Our findings provide promising insights into lncRNAs associated with CAFs in PAAD and provide a personalized tool for predicting patient prognosis and immune microenvironmental landscape.