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.

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

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

Gene expression profile of immune cells from patients with nasopharyngeal carcinoma of stage T3-4N0]

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

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

Microarray Analysis of Novel Genes Involved in Nasopharyngeal Carcinoma

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

Cell Division Cycle-Associated Genes Are Potential Immune Regulators in Nasopharyngeal Carcinoma

[Analysis of Gene Expression Patterns of Periphery Lymphocytes in Patients with Nasopharyngeal Carcinoma].

Integrated Analysis of the Differential Cellular and EBV Mirna Expression Profiles in Microdissected Nasopharyngeal Carcinoma and Non-Cancerous Nasopharyngeal Tissues

Analysis of the Expression of Surface Receptors on NK Cells and NKG2D on Immunocytes in Peripheral Blood of Patients with Nasopharyngeal Carcinoma

Novel prognostic biomarkers in nasopharyngeal carcinoma unveiled by mega-data bioinformatics analysis

An individualized immune prognostic signature in nasopharyngeal carcinoma

Comprehensive analysis of gene expression and DNA methylation for human nasopharyngeal carcinoma

Identification of shared immune infiltration characteristic molecules in dermatomyositis and nasopharyngeal carcinoma using bioinformatics: Traits in dermatomyositis and nasopharyngeal cancer

Expression and Clinical Significance of Programmed Death Ligand 1 in Nasopharyngeal Carcinoma.

Comprehensive analysis of key genes associated with ceRNA networks in nasopharyngeal carcinoma based on bioinformatics analysis

Prognostic Value of Immune Score in Nasopharyngeal Carcinoma Using Digital Pathology

The identification of key genes in nasopharyngeal carcinoma by bioinformatics analysis of high-throughput data

Integrating iron metabolism-related gene signature to evaluate prognosis and immune infiltration in nasopharyngeal carcinoma