N6-Methylandenosine-Related lncRNAs Play an Important Role in the Prognosis and Immune Microenvironment of Pancreatic Ductal Adenocarcinoma

Background: Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive, fatal tumor. N6-methylandenosine (m6A) methylation is the major epigenetic modification of RNA including lncRNAs. The roles of m6A-related lncRNAs in PDAC have not been fully clarified. The aim of this study is to assess gene signatures and prognostic value of m6A-related lncRNAs in PDAC. Methods The Cancer Genome Atlas (TCGA) dataset and the International Cancer Genome Consortium (ICGC) dataset were explored to identify m6A-related lncRNAs. Univariate, least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression were performed to construct the m6A-related lncRNAs prognostic riskscore (m6A-LPR) model to predict the overall survival (OS) in the TCGA training cohort. Kaplan–Meier curve with log-rank test and receiver operating characteristic (ROC) curve were used to evaluate the prognostic value of the m6A-LPR. Furthermore, the robustness of the m6A-LPR was further validated in the ICGC cohort. Tumor immunity was evaluated using ESTIMATE and CIBERSORT algorithms. Results A total of 262 m6A-related lncRNAs were identified in two datasets. In the TCGA training cohort, 28 prognostic m6A-related lncRNAs were identified and the m6A-LPR including four m6A-related lncRNAs was constructed. The m6A-LPR was able to identify high-risk patients with significantly poorer OS and accurately predict OS in both the TCGA training cohort and the ICGC validation cohort. Analysis of tumor immunity revealed that high-risk group had remarkably lower stromal, immune, and ESTIMATE scores. Moreover, high-risk group was associated with significantly higher levels of plasma B cells and resting NK cells infiltration, and lower levels of infiltrating resting memory CD4 T cells, monocytes and resting mast cells. Conclusions Our study proposed a robust m6A-related prognostic signature of lncRNAs for predicting OS in PDAC, which provides some clues for further studies focusing on the mechanism process underlying m6A modification of lncRNAs.