Abstract
Tumor molecular data sets are becoming increasingly complex, making it nearly impossible for humans alone to effectively analyze them. Here, we demonstrate the power of using machine learning (ML) to analyze a single-cell, spatial, and highly multiplexed proteomic data set from human pancreatic cancer and reveal underlying biological mechanisms that may contribute to clinical outcomes. We designed a multiplex immunohistochemistry antibody panel to compare T-cell functionality and spatial localization in resected tumors from treatment-naïve patients with localized pancreatic ductal adenocarcinoma (PDAC) with resected tumors from a second cohort of patients treated with neoadjuvant agonistic CD40 (anti-CD40) monoclonal antibody therapy. In total, nearly 2.5 million cells from 306 tissue regions collected from 29 patients across both cohorts were assayed, and over 1,000 tumor microenvironment (TME) features were quantified. We then trained ML models to accurately predict anti-CD40 treatment status and disease-free survival (DFS) following anti-CD40 therapy based on TME features. Through downstream interpretation of the ML models' predictions, we found anti-CD40 therapy reduced canonical aspects of T-cell exhaustion within the TME, as compared with treatment-naïve TMEs. Using automated clustering approaches, we found improved DFS following anti-CD40 therapy correlated with an increased presence of CD44+CD4+ Th1 cells located specifically within cellular neighborhoods characterized by increased T-cell proliferation, antigen experience, and cytotoxicity in immune aggregates. Overall, our results demonstrate the utility of ML in molecular cancer immunology applications, highlight the impact of anti-CD40 therapy on T cells within the TME, and identify potential candidate biomarkers of DFS for anti-CD40-treated patients with PDAC.