A minimal gene set characterizes TIL specific for diverse tumor antigens across different cancer types
Zhen Zeng, Tianbei Zhang, Jiajia Zhang, Shuai Li, Sydney Connor, Boyang Zhang, Yimin Zhao, J. L. Wilson, Dipika Singh, Rima M. Kulikauskas, Candice D. Church, Thomas H. Pulliam, Saumya Jani, Paul Nghiem, Suzanne L. Topalian, Patrick M. Forde, Drew M. Pardoll, Hongkai Ji, Kellie N. Smith
Abstract
Identifying tumor-specific T cell clones that mediate immunotherapy responses remains challenging. Mutation-associated neoantigen (MANA) -specific CD8+ tumor-infiltrating lymphocytes (TIL) have been shown to express high levels of CXCL13 and CD39 (ENTPD1), and low IL-7 receptor (IL7R) levels in many cancer types, but their collective relevance to T cell functionality has not been established. Here we present an integrative tool to identify MANA-specific TIL using weighted expression levels of these three genes in lung cancer and melanoma single-cell RNAseq datasets. Our three-gene "MANAscore" algorithm outperforms other RNAseq-based algorithms in identifying validated neoantigen-specific CD8+ clones, and accurately identifies TILs that recognize other classes of tumor antigens, including cancer testis antigens, endogenous retroviruses and viral oncogenes. Most of these TIL are characterized by a tissue resident memory gene expression program. Putative tumor-reactive cells (pTRC) identified via MANAscore in anti-PD-1-treated lung tumors had higher expression of checkpoint and cytotoxicity-related genes relative to putative non-tumor-reactive cells. pTRC in pathologically responding tumors showed distinguished gene expression patterns and trajectories. Collectively, we show that MANAscore is a robust tool that can greatly enrich candidate tumor-specific T cells and be used to understand the functional programming of tumor-reactive TIL.