N6-methyladenosine methylation regulates the tumor microenvironment of Epstein-Barr virus-associated gastric cancer

Abstract

BACKGROUND N6-methyladenosine (m6A) methylation modification exists in Epstein-Barr virus (EBV) primary infection, latency, and lytic reactivation. It also modifies EBV latent genes and lytic genes. EBV-associated gastric cancer (EBVaGC) is a distinctive molecular subtype of GC. We hypothesized EBV and m6A methylation regulators interact with each other in EBVaGC to differentiate it from other types of GC. AIM To investigate the mechanisms of m6A methylation regulators in EBVaGC to determine the differentiating factors from other types of GC. METHODS First, The Cancer Gene Atlas and Gene Expression Omnibus databases were used to analyze the expression pattern of m6A methylation regulators between EBVaGC and EBV-negative GC (EBVnGC). Second, we identified Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment of m6A-related differentially expressed genes. We quantified the relative abundance of immune cells and inflammatory factors in the tumor microenvironment (TME). Finally, cell counting kit-8 cell proliferation test, transwell test, and flow cytometry were used to verify the effect of insulin-like growth factor binding protein 1 (IGFBP1) in EBVaGC cell lines. RESULTS m6A methylation regulators were involved in the occurrence and development of EBVaGC. Compared with EBVnGC, the expression levels of m6A methylation regulators Wilms tumor 1-associated protein, RNA binding motif protein 15B, CBL proto-oncogene like 1, leucine rich pentatricopeptide repeat containing, heterogeneous nuclear ribonucleoprotein A2B1, IGFBP1, and insulin-like growth factor 2 binding protein 1 were significantly downregulated in EBVaGC (P < 0.05). The overall survival rate of EBVaGC patients with a lower expression level of IGFBP1 was significantly higher (P = 0.046). GO and KEGG functional enrichment analyses showed that the immunity pathways were significantly activated and rich in immune cell infiltration in EBVaGC. Compared with EBVnGC, the infiltration of activated CD4+ T cells, activated CD8+ T cells, monocytes, activated dendritic cells, and plasmacytoid dendritic cells were significantly upregulated in EBVaGC (P < 0.001). In EBVaGC, the expression level of proinflammatory factors interleukin (IL)-17, IL-21, and interferon-γ and immunosuppressive factor IL-10 were significantly increased (P < 0.05). In vitro experiments demonstrated that the expression level of IGFBP1 was significantly lower in an EBVaGC cell line (SNU719) than in an EBVnGC cell line (AGS) (P < 0.05). IGFBP1 overexpression significantly attenuated proliferation and migration and promoted the apoptosis levels in SNU719. Interfering IGFBP1 significantly promoted proliferation and migration and attenuated the apoptosis levels in AGS. CONCLUSION m6A regulators could remodel the TME of EBVaGC, which is classified as an immune-inflamed phenotype and referred to as a “hot” tumor. Among these regulators, we demonstrated that IGFBP1 affected proliferation, migration, and apoptosis.