Single-cell and bulk RNA sequencing identifies T cell marker genes score to predict the prognosis of pancreatic ductal adenocarcinoma

Abstract Pancreatic ductal adenocarcinoma (PDAC) is one of the lethal malignancies, with limited biomarkers identified to predict its prognosis and treatment response of immune checkpoint blockade (ICB). This study aimed to explore the predictive ability of T cell marker genes score (TMGS) to predict their overall survival (OS) and treatment response to ICB by integrating single-cell RNA sequencing (scRNA- seq ) and bulk RNA- seq data. Multi-omics data of PDAC were applied in this study. The uniform manifold approximation and projection (UMAP) was utilized for dimensionality reduction and cluster identification. The non-negative matrix factorization (NMF) algorithm was applied to molecular subtypes clustering. The Least Absolute Shrinkage and Selection Operator (LASSO)-Cox regression was adopted for TMGS construction. The prognosis, biological characteristics, mutation profile, and immune function status between different groups were compared. Two molecular subtypes were identified via NMF: proliferative PDAC (C1) and immune PDAC (C2). Distinct prognoses and biological characteristics were observed between them. TMGS was developed based on 10 T cell marker genes (TMGs) through LASSO-Cox regression. TMGS is an independent prognostic factor of OS in PDAC. Enrichment analysis indicated that cell cycle and cell proliferation-related pathways are significantly enriched in the high-TMGS group. Besides, high-TMGS is related to more frequent KRAS , TP53 , and CDKN2A germline mutations than the low-TMGS group. Furthermore, high-TMGS is significantly associated with attenuated antitumor immunity and reduced immune cell infiltration compared to the low-TMGS group. However, high TMGS is correlated to higher tumor mutation burden (TMB), a low expression level of inhibitory immune checkpoint molecules, and a low immune dysfunction score, thus having a higher ICB response rate. On the contrary, low TMGS is related to a favorable response rate to chemotherapeutic agents and targeted therapy. By combining scRNA- seq and bulk RNA- seq data, we identified a novel biomarker, TMGS, which has remarkable performance in predicting the prognosis and guiding the treatment pattern for patients with PDAC.