Abstract
Pancreatic cancer is a highly aggressive disease with poor prognosis. N6-methyladenosine (m6A) is critical for post-transcriptional modification of messenger RNA (mRNA) and long non-coding RNA (lncRNA). However, the m6A-associated lncRNAs (m6A-lncRNA) and their values in predicting clinical outcomes and immune microenvironmental status in pancreatic cancer patients remain largely unexplored. This study aimed to evaluate the importance of m6A-lncRNA and established a m6A-lncRNA signature for predicting immunotherapeutic response and prognosis of pancreatic cancer. The m6A-lncRNA co-expression networks were constructed using data from the TCGA and GTEx database. Based on the least absolute shrinkage and selection operator (LASSO) analysis, we constructed an 8 m6A-lncRNA signature risk model, and selection operator (LASSO) analysis, and stratified patients into the high- and low-risk groups with significant difference in overall survival (OS) (HR = 2.68, 95% CI = 1.74-4.14, P < 0.0001). Patients in the high-risk group showed significantly reduced OS compared to patients in the low-risk group (P < 0.001). The clinical characteristics and m6A-lncRNA risk scores were used to construct a nomogram which accurately predicted the OS in pancreatic cancer. TIMER 2.0 were used to investigate tumor immune infiltrating cells and its relationship with pancreatic cancer. CIBERSORT analysis revealed increased higher infiltration proportions of M0 and M2 macrophages, and lower infiltration of naive B cell, CD8(+) T cell and Treg cells in the high-risk group. Compared to the low-risk group, functional annotation using ssGSEA showed that T cell infiltration and the differential immune-related check-point genes are expressed at low level in the high-risk group (P < 0.05). In summary, our study constructed a novel m6A-associated lncRNAs signature to predict immunotherapeutic responses and provided a novel nomogram for the prognosis prediction of pancreatic cancer.