A potential EBV-related classifier is associated with the efficacy of immunotherapy in gastric cancer

Gastric cancer (GC) is a global health problem. As a complicated heterogeneous disease, The Cancer Genome Atlas (TCGA) Research Network recognized four subtypes of GC, including Epstein-Barr virus (EBV)-positive GC (EBVaGC), which accounts for approximately 9% of all GC cases. The response to immunotherapy in GC is limited, and whether EBV status is a predictor of immunotherapy remains controversial.

The differential gene expression analysis was utilized to compare the gene expression between the EBV-positive group and the EBV-negative group in the TCGA-Stomach Adenocarcinoma (TCGA-STAD) cohort. Weighted gene co-expression network analysis (WGCNA), protein-protein interaction (PPI) network analysis, and gene functional enrichment analysis were used to investigate the most pivotal hub genes and their roles. The "Estimation of Stromal and Immune cells in Malignant Tumours using Expression data" (ESTIMATE) and CIBERSORT algorithms were performed to infer the immune compositions of tissue samples. Furthermore, quantitative real-time polymerase chain reaction (RT-qPCR) and survival analysis were used to validate the expression and prognosis of these hub genes in Sun Yat-sen University Cancer Center (SYSUCC) cohort. A GC cohort that received anti-programmed cell death 1 (PD-1) therapy was used to analyze the predicted efficacy based on the expression of hub genes.

There is a total of 1,686 differentially expressed genes (DEGs) between the EBV-positive group and EBV-negative group, and WGCNA identified a yellow-colored module that was most related to EBV features. Functional enrichment analysis of 144 genes in this yellow module demonstrated that these genes primarily performed immune-related functions, and PPI network analysis through the CytoHubba plug-in identified 11 hub genes in the network. The RT-qPCR results in the SYSUCC cohort further validated that the hub genes were all increased in the EBV-positive group. Finally, we found that a potential classifier was associated with the efficacy of immunotherapy based on the expression of these 11 hub genes.

Our study identified several hub genes associated with EBV status that are closely related to the immune microenvironmental features of EBVaGC and may be used as molecular markers for predicting the immune response in GC.