High Throughput Transcriptome Data Analysis and Computational Verification Reveal Immunotherapy Biomarkers of Compound Kushen Injection for Treating Triple-Negative Breast Cancer

BACKGROUND: Although notable therapeutic and prognostic benefits of compound kushen injection (CKI) have been found when it was used alone or in combination with chemotherapy or radiotherapy for triple-negative breast cancer (TNBC) treatment, the effects of CKI on TNBC microenvironment remain largely unclear. This study aims to construct and validate a predictive immunotherapy signature of CKI on TNBC. METHODS: The UPLC-Q-TOF-MS technology was firstly used to investigate major constituents of CKI. RNA sequencing data of CKI-perturbed TNBC cells were analyzed to detect differential expression genes (DEGs), and the GSVA algorithm was applied to explore significantly changed pathways regulated by CKI. Additionally, the ssGSEA algorithm was used to quantify immune cell abundance in TNBC patients, and these patients were classified into distinct immune infiltration subgroups by unsupervised clustering. Then, prognosis-related genes were screened from DEGs among these subgroups and were further overlapped with the DEGs regulated by CKI. Finally, a predictive immunotherapy signature of CKI on TNBC was constructed based on the LASSO regression algorithm to predict mortality risks of TNBC patients, and the signature was also validated in another TNBC cohort. RESULTS: Twenty-three chemical components in CKI were identified by UPLC-Q-TOF-MS analysis. A total of 3692 DEGs were detected in CKI-treated versus control groups, and CKI significantly activated biological processes associated with activation of T, natural killer and natural killer T cells. Three immune cell infiltration subgroups with 1593 DEGs were identified in TNBC patients. Then, two genes that can be down-regulated by CKI with hazard ratio (HR) > 1 and 26 genes that can be up-regulated by CKI with HR < 1 were selected as key immune- and prognosis-related genes regulated by CKI. Lastly, a five-gene prognostic signature comprising two risky genes (MARVELD2 and DYNC2I2) that can be down-regulated by CKI and three protective genes (RASSF2, FERMT3 and RASSF5) that can be up-regulated by CKI was developed, and it showed a good performance in both training and test sets. CONCLUSIONS: This study proposes a predictive immunotherapy signature of CKI on TNBC, which would provide more evidence for survival prediction and treatment guidance in TNBC as well as a paradigm for exploring immunotherapy biomarkers in compound medicines.