Abstract
BACKGROUND: N6-methyladenosine (m6A) is the most abundant mRNA modification. Whether m6A regulators can determine tumor aggressiveness and risk of immune evasion in pancreatic ductal adenocarcinoma (PDAC) remains unknown. METHODS: An integrated model named "m6Ascore" is constructed based on RNA-seq data of m6A regulators in PDAC. Association of m6Ascore and overall survival is validated across several different datasets. Overlaps of m6Ascore and established molecular classifications of PDAC is examined. Immune infiltration, enriched pathways, somatic copy number alterations (SCNAs), mutation profiles and response to immune checkpoint inhibitors are compared between m6Ascore-high and m6Ascore-low tumors. FINDINGS: m6Ascore is associated with dismal overall survival and increased tumor recurrence in PDAC as well as several other solid tumors including colorectal cancer and breast cancer. Basal-like (Squamous) PDAC has higher m6Ascore than that in the classical PDAC. Mechanism study showed m6Ascore-high tumors are characterized with reduced immune infiltration and T cells exhaustion. Meanwhile, m6Ascore is associated with genes regulating cachexia and chemoresistance in PDAC. Furthermore, distinct SCNAs patterns and mutation profiles of KRAS and TP53 are present in m6Ascore-high tumors, indicating immune evasion. m6Ascore-low tumors have higher response rates to immune checkpoint inhibitors (ICIs). INTERPRETATION: These findings indicate m6Ascore can predict aggressiveness and immune evasion in pancreatic cancer. This model has implications for pancreatic cancer prognosis and treatment response to ICIs. FUNDING: This work was supported in part by National Institutes of Health (NIH) grants to M. Li (R01 CA186338, R01 CA203108, R01 CA247234 and the William and Ella Owens Medical Research Foundation) and NIH/National Cancer Institute Q39 award P30CA225520 to Stephenson Cancer Center.