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.

Identification of Four Key Biomarkers and Small Molecule Drugs Innasopharyngeal Carcinoma by Weighted Gene Co-expression Network Analysis

Identification of four key biomarkers and small molecule drugs in nasopharyngeal carcinoma by weighted gene co-expression network analysis

Identification of Prognostic Biomarkers for Response to Radiotherapy by DNA Microarray in Nasopharyngeal Carcinoma Patients

Identifying key genes and small molecule compounds for nasopharyngeal carcinoma by various bioinformatic analysis

Bioinformatic identification of candidate biomarkers and related transcription factors in nasopharyngeal carcinoma

Identification of Novel Kinase-Transcription Factor-mRNA-miRNA Regulatory Network in Nasopharyngeal Carcinoma by Bioinformatics Analysis

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

Identification of key pathways and genes in nasopharyngeal carcinoma using bioinformatics analysis

Identification of Key Genes in Nasopharyngeal Carcinoma Based on Bioinformatics Analysis

Screening Key Genes and Biological Pathways in Nasopharyngeal Carcinoma by Integrated Bioinformatics Analysis

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

Identification of novel biomarkers and small-molecule compounds for nasopharyngeal carcinoma with metastasis

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

Bioinformatics analysis identifies hub genes and pathways in nasopharyngeal carcinoma

Identification of markers for predicting prognosis and endocrine metabolism in nasopharyngeal carcinoma by miRNA-mRNA network mining and machine learning

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

Functional analysis of the nasopharyngeal carcinoma primary tumor-associated gene interaction network

[Key Prediction Genes of Nasopharyngeal Carcinoma:Screening Based on Systematic Bioinformatics and Validation by Cell Experiments]

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

Microarray Analysis of Novel Genes Involved in Nasopharyngeal Carcinoma

A practical method to screen and identify functioning biomarkers in nasopharyngeal carcinoma