Human natural killer cells in major histocompatibility complex class I deficiency

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

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 (log₂ 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.

B2M overexpression correlates with malignancy and immune signatures in human gliomas

Because of the limited treatment strategy of gliomas, the key of diagnosis and treatment is finding new molecular biomarkers. Here, we explored the potential of β2-microglobulin (B2M) to serve as a hopeful candidate for immunotherapy or diagnostic biomarker in gliomas. The genomic profiles, clinical characteristics, and immune signatures were analyzed based on TCGA and CGGA databases. We carried out the whole statistical analyses using R project. High B2M expression correlated with worse prognosis. Somatic mutations of gliomas with high B2M expression are associated with PTEN deletion and EGFR amplification. Isocitrate dehydrogenase (IDH) mutations accounted for 82% in gliomas with low B2M expression. In addition, B2M positively correlated with ESTIMATE scores, interacted with infiltrating immune and stromal cell types. B2M also suppressed anti-tumor immunity through immune related processes. Meanwhile, B2M was associated with immune checkpoint molecules and inflammatory activities. Finally, functional annotation of the identified B2M related genes verified that B2M was a potential candidate for immunotherapy. We confirmed that B2M played a critical role in tumor progression, patient prognosis and immunotherapy of gliomas.