Abstract:Background: The prognostic value of tumor-infiltrating immune cells has been widely studied in nasopharyngeal carcinoma (NPC). However, the role of tumor-infiltrating immune cells in the diagnosis of NPC has not been fully elucidated. Thus, tumor-infiltrating immune cell-related biomarkers in the diagnosis of NPC patients were explored in the current study. Methods: Gene expression profiles of NPC patients were downloaded from the Gene Expression Omnibus (GEO) database. Differentially infiltrating immune cells (DDICs) between NPC and control samples were analyzed by the CIBERSORT algorithm. Weighted gene coexpression network analysis (WGCNA) was performed to screen hub genes significantly correlated with DDIC. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of hub genes were performed with R package clusterProfiler. The diagnostic value of hub genes was evaluated by receiver operating characteristic (ROC) curves. RT-qPCR was conducted to validate the expression patterns of diagnostic markers in NPC and adjacent control tissues. The correlations between diagnostic markers and immunomodulators were analyzed using the TISIDB. The protein-protein interaction (PPI) network based on immunomodulators significantly associated with diagnostic biomarkers was constructed and visualized by STRING. The functional enrichment analysis of genes in the PPI network was analyzed by the WebGestalt online tool. Results: The abundances of memory B cells, plasma cells, follicular helper T cells, activated NK cells, M0 macrophages, M1 macrophages, M2 macrophages, resting mast cells, and activated mast cells were significantly different between NPC and control samples. Dark orange was identified as the hub module, with a total of 371 genes associated with memory B cells, plasma cells, and M0 and M1 macrophages defined as hub genes, which were enriched into immune-related biological processes and pathways. FCER2, KHDRBS2, and IGSF9 were considered diagnostic biomarkers with areas under ROC curves as 0.985, 0.978, and 0.975, respectively. Moreover, real-time reverse transcriptase-polymerase chain reaction (RT-qPCR) suggested that the expression patterns of FCER2, KHDRBS2, and IGSF9 were consistent with the results in GEO datasets. TISIDB analysis revealed that FCER2, KHDRBS2, and IGSF9 had a strong association with 8 immunoinhibitors (BTLA, CD160, CD96, LAG3, PDCD1, TIGIT, CD244, and TGFB1) and 11 immunostimulators (CD27, CD28, CD40LG, CD48, ICOS, KLRC1, KLRK1, TMIGD2, TNFRSF13C, CXCR4, and C10 or f54). The PPI network implied that these 19 immunomodulators had interactions with other 50 genes. WebGestalt analysis demonstrated that 69 genes in the PPI network were enriched into cytokine-cytokine receptor interaction, NF-kappa B signaling pathway, and pathways in cancer. Conclusion: Our study identified novel diagnostic biomarkers and revealed potential immune-related mechanisms in NPC. These findings enlighten the investigation of the molecular mechanisms of tumor-infiltrating immune cells regulating NPC.

Infiltrating B-cell Subtypes and Associated Hub Genes in Nasopharyngeal Carcinoma Identified from Integrated Single-Cell and Bulk-Rna Sequencing Data.

The role of tertiary lymphoid structure and B cells in nasopharyngeal carcinoma: Based on bioinformatics and experimental verification

Detailed analysis of inflammatory cell infiltration and the prognostic impact on nasopharyngeal carcinoma

[Immunocyte infiltration characteristics of gene expression profile in nasopharyngeal carcinoma and clinical significance]

Single-cell transcriptomics reveals regulators underlying immune cell diversity and immune subtypes associated with prognosis in nasopharyngeal carcinoma

Peripheral NK-cell immunophenotypes predict progression in nasopharyngeal carcinoma patients.

Comprehensive single-cell sequencing reveals the stromal dynamics and tumor-specific characteristics in the microenvironment of nasopharyngeal carcinoma

Integrated analysis of bulk and single-cell RNA sequencing reveals the interaction of PKP1 and tumor-infiltrating B cells and their therapeutic potential for nasopharyngeal carcinoma

scRNA+TCR+BCR-seq revealed the proportions and gene expression patterns of dual receptor T and B lymphocytes in NPC and NLH

Genomic Analysis of Nasopharyngeal Carcinoma Reveals TME-Based Subtypes

Prognosis-related Molecular Subtyping of Nasopharyngeal Carcinoma

Differentially Infiltrated Identification of Novel Diagnostic Biomarkers Associated with Immune Infiltration in Nasopharyngeal Carcinoma

Single cell RNA-seq reveals the landscape of tumor and infiltrating immune cells in nasopharyngeal carcinoma

PD-1 + CXCR5 − CD4 + Th-CXCL13 cell subset drives B cells into tertiary lymphoid structures of nasopharyngeal carcinoma

Spatial Heterogeneity of Immune Infiltration Predicts the Prognosis of Nasopharyngeal Carcinoma Patients

Multimodal profiling uncovers tertiary lymphoid structures as a critical determinant of immunotherapy response and prognosis in nasopharyngeal carcinoma

Tumour heterogeneity and intercellular networks of nasopharyngeal carcinoma at single cell resolution

Integrative analysis of BSG expression in NPC through immunohistochemistry and public high-throughput gene expression data

CD8+ Tumor-Infiltrating Lymphocyte Abundance Is a Positive Prognostic Indicator in Nasopharyngeal Cancer

Tumor cells PD-L1 expression as a favorable prognosis factor in nasopharyngeal carcinoma patients with pre-existing intratumor-infiltrating lymphocytes.

Investigation of differentially expressed genes in nasopharyngeal carcinoma by integrated bioinformatics analysis