Comprehensive single-cell transcriptomic and proteomic analysis reveals NK cell exhaustion and unique tumor cell evolutionary trajectory in non-keratinizing nasopharyngeal carcinoma

Abstract

Abstract Background Nonkeratinizing nasopharyngeal carcinoma (NK-NPC) has a strong association with Epstein-Barr virus (EBV) infection. The role of NK cells and the tumor cell evolutionary trajectory in NK-NPC remain unclear. In this study, we aim to investigate the function of NK cell and the evolutionary trajectory of tumor cells in NK-NPC by single-cell transcriptomic analysis, proteomics and immunohistochemistry. Methods NK-NPC (n = 3) and normal nasopharyngeal mucosa cases (n = 3) were collected for proteomic analysis. Single-cell transcriptomic data of NK-NPC (n = 10) and nasopharyngeal lymphatic hyperplasia (NLH, n = 3) were obtained from Gene Expression Omnibus (GSE162025 and GSE150825). Quality control, dimension reduction and clustering were based on Seurat software (v4.0.2) process and batch effects were removed by harmony (v0.1.1) software. Normal cells of nasopharyngeal mucosa and tumor cells of NK-NPC were identified using copykat software (v1.0.8). Cell-cell interactions were explored using CellChat software (v1.4.0). Tumor cell evolutionary trajectory analysis was performed using SCORPIUS software (v1.0.8). Protein and gene function enrichment analyses were performed using clusterProfiler software (v4.2.2). Results A total of 161 differentially expressed proteins were obtained between NK-NPC (n = 3) and normal nasopharyngeal mucosa (n = 3) by proteomics (log2 fold change > 0.5 and P value < 0.05). Most of proteins associated with the nature killer cell mediated cytotoxicity pathway were downregulated in the NK-NPC group. In single cell transcriptomics, we identified three NK cell subsets (NK1-3), among which NK cell exhaustion was identified in the NK3 subset with high ZNF683 expression (a signature of tissue-resident NK cell) in NK-NPC. We demonstrated the presence of this ZNF683 + NK cell subset in NK-NPC but not in NLH. We also performed immunohistochemical experiments with TIGIT and LAG3 to confirm NK cell exhaustion in NK-NPC. Moreover, the trajectory analysis revealed that the evolutionary trajectory of NK-NPC tumor cells was associated with the status of EBV infection (active or latent). The analysis of cell-cell interactions uncovered a complex network of cellular interactions in NK-NPC. Conclusions This study revealed that the NK cell exhaustion might be induced by upregulation of inhibitory receptors on the surface of NK cells in NK-NPC. Treatments for the reversal of NK cell exhaustion may be a promising strategy for NK-NPC. Meanwhile, we identified a unique evolutionary trajectory of tumor cells with active status of EBV-infection in NK-NPC for the first time. Our study may provide new immunotherapeutic targets and new sight of evolutionary trajectory involving tumor genesis, development and metastasis in NK-NPC.