Abstract
Background: Anti-programmed cell death 1/programmed cell death ligand 1 (PD1/PDL1) therapy is an important part of comprehensive cancer therapy. However, many patients suffer from non-response to therapy. Tumor neoantigen burden (TNB) and cancer stemness play essential roles in the responsiveness to therapy. Therefore, the identification of drug candidates for anti-PD1/PDL1 therapy remains an unmet need. Methods: Three anti-PD1/PDL1 therapy cohorts were obtained from GEO database and published literatures. Cancer immune characteristics were analyzed using CIBERSORTX, GSVA, and ESTIMATE. WGCNA was employed to identify the gene modules correlated with cancer TNB and stemness. A machine-learning method was used to construct the immunotherapy resistance score (TSIRS). Pharmacogenomic analysis was conducted to explore the potential alternative drugs for anti-PD1/PDL1 therapy resistant patients. CCK-8 assay, EdU assay and wound healing assay were used to validate the effect of the predicted drug on cancer cells. Results: The therapy response and non-response cancer groups have different microenvironment features. TSIRS was developed based on tumor neoantigen and stemness. TSIRS can effectively predict the outcomes of patients with anti-PD1/PDL1 therapy in training, validation and meta cohorts. Meanwhile, TSIRS can reflect the characteristics of tumor microenvironment during anti-PD1/PDL1 therapy. PF-4708671 is identified as a potential alternative drug for patients with resistance to anti-PD1/PDL1 therapy. It possesses significant inhibitive effect on the proliferation and migration of BGC-823 cells. Conclusion: TSIRS is an effective tool in the identification of candidate patients who will be benefit from anti-PD1/PDL1 therapy. Small molecule drug PF-4708671 has the potential to be used in anti-PD1/PDL1 therapy resistant patients.