Disease Markers
Volume 2019 (2019), Article ID 1372571, 13 pages
https://doi.org/10.1155/2019/1372571

Research Article

Bioinformatics Analysis of Potential Key Genes in Trastuzumab-Resistant Gastric Cancer

Guangda Yang,¹ Liumeng Jian,² Xiangan Lin,³ Aiyu Zhu,¹ and Guohua Wen¹

¹Department of Cancer Chemotherapy, Zengcheng District People’s Hospital of Guangzhou (BoJi-Affiliated Hospital of Sun Yat-sen University), Guangzhou, China
²Department of Neurology, Zengcheng District People’s Hospital of GuangZhou (BoJi-Affiliated Hospital of Sun Yat-sen University), Guangzhou, China
³Department of Cancer Chemotherapy, Sun Yat-sen Memorial Hospital of Sun Yat-sen University, Guangzhou, China

Correspondence should be addressed to Guangda Yang; yangguangda.123com@163.com

Received 6 July 2019; Revised 3 October 2019; Accepted 9 November 2019; Published 20 December 2019

Academic Editor: Michele Malaguarnera

Copyright © 2019 Guangda Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. This study was performed to identify genes related to acquired trastuzumab resistance in gastric cancer (GC) and to analyze their prognostic value. Methods. The gene expression profile GSE77346 was downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were obtained by using GEO2R. Functional and pathway enrichment was analyzed by using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Search Tool for the Retrieval of Interacting Genes (STRING), Cytoscape, and MCODE were then used to construct the protein-protein interaction (PPI) network and identify hub genes. Finally, the relationship between hub genes and overall survival (OS) was analyzed by using the online Kaplan-Meier plotter tool. Results. A total of 327 DEGs were screened and were mainly enriched in terms related to pathways in cancer, signaling pathways regulating stem cell pluripotency, HTLV-I infection, and ECM-receptor interactions. A PPI network was constructed, and 18 hub genes (including one upregulated gene and seventeen downregulated genes) were identified based on the degrees and MCODE scores of the PPI network. Finally, the expression of four hub genes (ERBB2, VIM, EGR1, and PSMB8) was found to be related to the prognosis of HER2-positive (HER2+) gastric cancer. However, the prognostic value of the other hub genes was controversial; interestingly, most of these genes were interferon- (IFN-) stimulated genes (ISGs). Conclusions. Overall, we propose that the four hub genes may be potential targets in trastuzumab-resistant gastric cancer and that ISGs may play a key role in promoting trastuzumab resistance in GC.

1. Introduction

Gastric cancer is the fifth most commonly diagnosed cancer and the third leading cause of cancer-related deaths [1]. The majority of gastric cancer cases are associated with lifestyle factors [2] and infectious agents, including the bacterium Helicobacter pylori [2, 3] and Epstein-Barr virus (EBV) [4, 5]. Although many biomarkers (including HER2, E-cadherin, fibroblast growth factor receptor, PD-L1, and TP53) have been studied as prognostic markers, the 5-year survival rate of gastric cancer remains low [6].

The human epidermal growth factor receptor-2 (HER-2) gene, a proto-oncogene mapped to chromosome 17 (17q12–q21), is frequently found to be amplified and/or overexpressed in gastric cancer [7]. Additionally, HER2 positivity is often associated with a worse prognosis [8, 9]. A phase III trial (the ToGA trial) confirmed that trastuzumab, a HER-2 monoclonal antibody, markedly improved the outcome of HER-2-positive (HER2+) gastric cancer patients [10]. However, a large proportion of patients developed resistance to trastuzumab after continuous treatment despite the effectiveness of this therapeutic [11]. Thus, there is an urgent need to explore the molecular mechanisms of trastuzumab resistance in gastric cancer and to identify effective biomarkers.

Bioinformatics analysis has been widely used to identify key genes in cancer. Interestingly, Piro et al. obtained the gene expression profiles of trastuzumab-sensitive and trastuzumab-resistant cell lines and found that fibroblast growth factor receptor 3 (FGFR3) was associated with trastuzumab resistance in gastric cancer [12]. In the present study, we aimed to further screen DEGs and predict their underlying function by utilizing the same data. More importantly, hub genes affecting trastuzumab resistance in GC patients were identified by a using protein-protein interaction (PPI) network, PPI network modules, and survival analyses.

2. Materials and Methods

2.1. Microarray Data

The microarray data for GSE77346 deposited by Piro et al. into the GEO database were obtained on the GPL10558 platform (Illumina HumanHT-12 v4.0 Expression BeadChip). The expression profiles are provided for five samples, including one sample of a trastuzumab-sensitive cell line (NCI-N87) and four samples of trastuzumab-resistant cell lines (N87-TR1, N87-TR2, N87-TR3, and N87-TR4).

2.2. Identification of DEGs

The web tool GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r/) was utilized to screen differentially expressed genes (DEGs) between trastuzumab-resistant and trastuzumab-sensitive gastric cancer cells. These DEGs were identified as important genes that may play an important role in the development of gastric cancer. The cutoff criterion were and .

2.3. Functional and Pathway Enrichment Analysis

We performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses by using the Database for Annotation, Visualization, and Integrated Discovery (DAVID), which is a comprehensive set of functional annotation tools. A value of <0.05 was set as the cutoff criterion.

2.4. PPI Network and Module Selection

The Search Tool for the Retrieval of Interacting Genes (STRING) database was used to analyze the PPI network of DEGs. Then, the results were visualized by Cytoscape software. The cutoff criterion for the combined score was >0.4. Subsequently, the degrees of genes in the PPI network and the MCODE plugin in Cytoscape were used to identify hub genes. Genes with or MCODE were identified as hub genes. Furthermore, MCODE was used to screen modules in the PPI network.

2.5. mRNA Expression of the Hub Genes

UALCAN (http://ualcan.path.uab.edu/index.html) [13], an interactive web portal, can be used to analyze the relative expression of a query gene across tumor and normal samples. The expression levels of the 18 hub genes were analyzed in gastric cancer and normal samples. The calculated value is shown.

2.6. Survival Analysis Based on the Hub Genes

The Kaplan-Meier plotter (KM plotter) tool (http://kmplot.com/analysis/index.php?p=background) was used to predict the prognostic value of the hub genes in gastric cancer patients [14]. The patients were divided into two groups according to the particular gene expression level (high vs. low expression). The OS of the two patient groups was then analyzed based on these categories. The hazard ratios (HRs) with 95% confidence intervals and log rank values are shown.

3. Results

3.1. Identification of DEGs

After data preprocessing, a total of 327 genes were identified, including 128 upregulated genes and 199 downregulated genes.

3.2. GO and KEGG Pathway Enrichment Analyses

We conducted GO and KEGG pathway enrichment analyses by using DAVID software. The top five GO terms of the DEGs are shown in Table 1. Regarding biological processes, DEGs were significantly involved in the type I interferon signaling pathway, negative regulation of viral genome replication, cell proliferation, positive regulation of the apoptotic process, and extracellular matrix organization. Regarding the cellular component, DEGs were enriched in the extracellular exosome, intermediate filament, cytoplasm, proteinaceous extracellular matrix, and basolateral plasma membrane. For molecular function, DEGs were enriched in structural molecule activity, calcium ion binding, protein homodimerization activity, 2-5-oligoadenylate synthetase activity, and fibronectin binding. The most enriched KEGG pathways included pathways in cancer, signaling pathways regulating stem cell pluripotency, HTLV-I infection, ECM-receptor interaction, and central carbon metabolism in cancer.

Table 1: Functional and pathway enrichment analyses of the DEGs in trastuzumab-resistant gastric cancer.

Category
Description
Count
value

GO function
BP
Type I interferon signaling pathway
13
1.04

BP
Negative regulation of viral genome replication
10
2.46

BP
Cell proliferation
23
5.99

BP
Positive regulation of apoptotic process
18
2.51

BP
Extracellular matrix organization
13
1.92

CC
Extracellular exosome
90
1.28

CC
Intermediate filament
11
2.26

CC
Cytoplasm
122
3.98

CC
Proteinaceous extracellular matrix
16
5.65

CC
Basolateral plasma membrane
13
5.98

MF
Structural molecule activity
17
8.93

MF
Calcium ion binding
27
4.52

MF
Protein homodimerization activity
26
1.31

MF
2-5-oligoadenylate synthetase activity
3
1.83

MF
Fibronectin binding
4
1.06

KEGG pathway

hsa05200
Pathways in cancer
18
5.47

hsa04550
Signaling pathways regulating pluripotency of stem cells
8
1.22

hsa05166
HTLV-I infection
11
1.56

hsa04512
ECM-receptor interaction
6
1.91

hsa05230
Central carbon metabolism in cancer
5
2.68

Note: top five terms were selected according to value. Abbreviation: DEGs: differentially expressed genes; GO: gene ontology; BP: biological process; CC: cellular component; MF: molecular function; HTLV-I: human T lymphocyte virus-I; ECM: extracellular matrix.

3.3. Construction of the PPI Network and Module Identification

The PPI network of the DEGs was constructed with 186 nodes and 480 edges by using the STRING database (Figure 1). was set as the cutoff criterion. The top 16 genes were ERBB2, OAS2, OASL, ISG15, OAS1, VIM, IFIT1, IFIT2, IFIT3, MX1, EGR1, IFI27, IFI44L, IFITM3, IFI44, and BST2 (Figure 2(a) and Table 2). Subsequently, a significant module with an MCODE was selected; this module had 16 nodes and 111 edges, and the included genes were PSMB8, OAS2, OASL, ISG15, OAS1, IFIT1, IFIT2, IFIT3, MX1, EGR1, IFI27, IFI44L, IFITM3, IFI44, BST2, and SAMD9 (Figure 2(b) and Table 2). Thus, 18 genes were identified as hub genes.

Figure 1: PPI network of the DEGs. Abbreviation: DEGs: differentially expressed genes.

(a)

(b)

Figure 2: The genes identified by degree (a) and MCODE score (b) in the PPI network.

Table 2: The hub genes in the PPI network.

Gene
Regulation
Degree
MCODE scores

ERBB2
Down
29
—

OAS2
Down
24
11.54

OASL
Down
23
11.54

ISG15
Down
23
11.54

OAS1
Down
23
11.54

VIM
Up
21
—

IFIT1
Down
19
11.54

IFIT2
Down
18
11.54

IFIT3
Down
18
11.54

MX1
Down
18
11.54

EGR1
Down
18
12.00

IFI27
Down
14
11.54

IFI44L
Down
14
12.00

IFITM3
Down
14
11.54

IFI44
Down
14
12.00

BST2
Down
14
11.54

PSMB8
Down
—
12.00

SAMD9
Down
—
10.00

3.4. mRNA Expression and Survival Analysis

UALCAN database was used to analyze the expression levels of the 18 hub genes. Compared with normal samples, primary gastric cancer samples had higher expression of ERBB2, PSMB8, IFI44, IFI44L, IFIT2, IFIT3, ISG15, OAS1, BST2, IFIT1, IFITM3, MX1, and OAS2. The expression of EGR1 was lower in primary gastric cancer samples than in normal samples. However, no significant difference in VIM, OASL, SAMD9, or IFI27 expression was observed between primary gastric cancer samples and normal samples (Figure S1).

Kaplan-Meier plotter was used to predict the prognostic value of the 18 hub genes. For all gastric cancer cases, our results showed that high ERBB2 expression was associated with the worse overall survival of GC patients, VIM, IFI44, IFIT2, and MX1 showed similar associations () (Figures 3(a) and 3(c) and Table 3). Additionally, low EGR1 expression was associated with the poorer overall survival of GC patients, and similar associations were found for PSMB8, ASMD9, BST2, IFI27, and IFIT1 () (Figures 4(a) and 4(c) and Table 3). For HER2-gastric cancer, the high expression of ERBB2, VIM, or IFI44 was associated with worse overall survival () (Table 3). In addition, the low expression of PSMB8, IFI44L, IFIT3, ISG15, OAS1, SAMD9, BST2, IFI27, IFIT1, or OAS2 was associated with the poorer overall survival of HER2-GC patients () (Table 3). However, the expression of EGR1, IFIT2, OASL, IFITM3, and MX1 was not associated with the overall survival of HER2-GC patients. For HER2+ gastric cancer, the high expression of VIM, IFI44, IFI44L, IFIT2, IFIT3, ISG15, OAS1, or OASL was associated with worse overall survival () (Figure 3(d) and Table 3). In addition, the low expression of ERBB2, EGR1, or PSMB8 was associated with the poorer overall survival of HER2+ GC patients () (Figures 3(b), 4(b), and 4(d) and Table 3). However, the expression of ASMD9, BST2, IFI27, IFIT1, IFITM3, MX1, and OAS2 was not associated with the overall survival of HER2+ GC patients (Table 3).

(a)

(b)

(c)

(d)

Figure 3: Kaplan-Meier curves depicting the overall survival of all patients with gastric cancer (GC) (a) or HER2-positive (HER2+) gastric cancer (b) with high or low expression of ERBB2 or VIM ((c) GC, (d) HER2+ GC).

Table 3: Overall survival analysis based on the expression of the hub genes.

Gene name
ID
OS (all)
OS (HER2-)
OS (HER2+)

HR (95% CI)
value
HR (95% CI)
value
HR (95% CI)
value

ERBB2
210930_s_at
1.37 (1.14-1.65)
0.00088
1.33 (1.06-1.66)
0.013
0.73 (0.55-0.97)
0.027

VIM
201426_s_at
1.29 (1.09-1.53)
0.0031
1.52 (1.21-1.9)
0.00027
1.4 (1.01-1.92)
0.04

EGR1
201694_s_at
0.63 (0.52-0.76)
8.5
0.78 (0.61-1)
0.052
0.6 (0.46-0.79)
0.00023

PSMB8
209040_s_at
0.45 (0.37-0.56)
6.6
0.39 (0.3-0.5)
1.7
0.62 (0.46-0.85)
0.0024

IFI44
214059_at
1.46 (1.22-1.74)
0.000028
1.39 (1.1-1.75)
0.0056
1.67 (1.28-2.19)
0.00015

IFI44L
204439_at
0.85 (0.71-1.03)
0.1
0.75 (0.58-0.98)
0.037
1.59 (1.16-2.19)
0.0037

IFIT2
217502_at
1.26 (1.04-1.53)
0.018
0.84 (0.67-1.05)
0.12
1.7 (1.24-2.33)
0.00075

IFIT3
204747_at
1.14 (0.95-1.36)
0.16
0.75 (0.59-0.95)
0.017
1.54 (1.15-2.07)
0.0037

ISG15
205483_s_at
0.86 (0.71-1.03)
0.094
0.73 (0.56-0.93)
0.012
1.5 (1.1-2.05)
0.01

OAS1
202869_at
1.15 (0.97-1.37)
0.1
0.79 (0.63-0.99)
0.039
1.39 (1.07-1.8)
0.013

OASL
210797_s_at
1.16 (0.98-1.37)
0.09
0.84 (0.67-1.05)
0.12
1.41 (1.04-1.92)
0.028

SAMD9
219691_at
0.73 (0.61-0.88)
0.00065
0.57 (0.44-0.74)
0.00002
1.23 (0.94-1.6)
0.14

BST2
201641_at
0.7 (0.58-0.85)
0.00027
0.55 (0.42-0.7)
1.9
1.28 (0.94-1.73)
0.12

IFI27
202411_at
0.7 (0.59-0.84)
0.000096
0.58 (0.46-0.73)
1.8
0.85 (0.63-1.13)
0.26

IFIT1
203153_at
0.73 (0.61-0.87)
0.00061
0.63 (0.49-0.81)
0.00025
1.23 (0.95-1.6)
0.12

IFITM3
212203x_at
1.09 (0.91-1.29)
0.36
0.89 (0.71-1.12)
0.32
1.27 (0.96-1.69)
0.096

MX1
202086_at
1.27 (1.05-1.54)
0.015
1.15 (0.9-1.47)
0.26
1.32 (0.99-1.76)
0.056

OAS2
204972_at
0.91 (0.77-1.08)
0.3
0.77 (0.61-0.96)
0.021
1.25 (0.93-1.67)
0.13

Abbreviation: OS: overall survival; HER2-: human epidermal growth factor receptor-2 negative; HER2+: human epidermal growth factor receptor-2 positive.

(a)

(b)

(c)

(d)

Figure 4: Kaplan-Meier curves depicting the overall survival of all patients with gastric cancer (GC) (a) or HER2-positive (HER2+) gastric cancer (b) with high or low expression of EGR1 or PSMB8 ((c) GC, (d) HER2+ GC).

We also analyzed the prognostic value of the expression of the 18 hub genes in gastric cancer with different clinicopathological characteristics, including gender, stage, differentiation, and treatment. As shown in Table S1, not all the hub genes had the prognostic value in gastric cancer with different parameters, but some had prognostic value in gastric cancer with specific clinicopathological characteristics. For example, the expression of VIM had prognostic value in almost all categories of gastric cancer. The expression of PSMB8, IFIT2, and IFIT1 had prognostic value in gastric cancer at different stages.

4. Discussion

In this study, a total of 327 DEGs were screened, including 128 upregulated genes and 199 downregulated genes. Eighteen genes were identified as hub genes, including one upregulated gene (VIM) and seventeen downregulated genes (ERBB2, PSMB8, OAS2, OASL, ISG15, OAS1, IFIT1, IFIT2, IFIT3, MX1, EGR1, IFI27, IFI44L, IFITM3, IFI44, BST2, and SAMD9). However, survival analysis based on the expression of these genes indicated that only one overexpressed gene (VIM) and three downregulated genes (ERBB2, EGR1, and PSMB8) were significantly associated with the poorer overall survival of HER2+ GC patients.

The data showed that the overexpression of ERBB2 (or HER2) was associated with the worse overall survival of GC patients. Some studies have confirmed that ERBB2 positivity is correlated with a worse prognosis [8, 9], but others have found no relationship between ERBB2 status and prognosis [15, 16]. Therefore, the relationship between ERBB2 status and the prognosis of GC patients remains controversial. However, ERBB2 was found to be downregulated in trastuzumab-resistant cells in the present study, and low ERBB2 expression was associated with a poorer prognosis of HER2+ GC patients. Interestingly, HER2 loss was observed in GC patients treated with trastuzumab in the clinic [17]. This phenomenon indicates that ERBB2 may play an important role in promoting resistance to trastuzumab. The present study found that the overexpression of VIM (vimentin) was correlated with a poorer prognosis of all GC patients, including HER2-GC patients and HER2+ GC patients. Importantly, high VIM expression was associated with a poorer prognosis of almost all categories of GC patients. VIM expression is required for epithelial-mesenchymal transition (EMT) [18]. Some studies have indicated that VIM overexpression is associated with a poorer prognosis among GC patients [19, 20]; this finding is consistent with our results and indicates that VIM may be involved in the development of trastuzumab resistance. The roles of proteasome subunit beta type-8 (PSMB8) and EGR1 in trastuzumab-resistant gastric cancer are controversial. We found that low expression of PSMB8 and EGR1 was associated with a poorer prognosis in all GC patients, including HER2-GC patients and HER2+ GC patients. However, the expression level of PSMB8 was higher in gastric cancer patients than in control patients. A previous study found that increased PSMB8 expression was associated with a lower survival rate of GC patients [21]. Similarly, increased expression of EGR1 was found to be significantly correlated with the depth of invasion and poorer survival of GC patients [22, 23]. Therefore, the roles of PSMB8 and EGR1 in trastuzumab-resistant gastric cancer need to be further investigated. More importantly, many drugs that target ERBB2, including lapatinib [24], trastuzumab [25], and pertuzumab [26], have been investigated. In addition, the ToGA trial confirmed that trastuzumab markedly improves the outcome of HER2+ gastric cancer patients [10]. Moreover, inhibitors of PSMB8 [27] or VIM [28] have been investigated. Carfilzomib, a novel, irreversible proteasome inhibitor, was shown to have potent activity against preclinical models of multiple myeloma [27]. Phenethyl isothiocyanate has been evaluated in trials studying the prevention and treatment of leukemia, lung cancer, tobacco use disorder, and lymphoproliferative disorders [28]. Together, ERBB2, VIM, and PSMB8 may be effective targets in gastric cancer, but more experimental investigations and clinical trials are needed.

IFI44, IFI44L, IFIT2, IFIT3, ISG15, OAS1, and OASL were identified in the subnetwork module. All of these genes are interferon- (IFN-) stimulated genes (ISGs) [29], which mediate the antiviral action of interferon. The interferon-induced protein 44 (IFI44) and interferon-induced protein 44-like (IFI44L) genes belong to the IFI44 family [30]. A previous study found that overexpression of IFI44L decreased doxorubicin chemoresistance and was associated with the better survival of hepatocellular carcinoma patients [31]. The IFN-induced protein with tetratricopeptide repeats (IFITs) family (including IFIT1, IFIT2, IFIT3/4, and IFIT5) is among hundreds of ISGs [32, 33]. Studies have shown that IFIT2 depletion induces cell migration and is associated with poor prognosis in patients with oral squamous cell carcinoma (OSCC) [34, 35]. Similarly, decreased IFIT2 expression predicted poor therapeutic outcomes of GC patients [36]. However, a previous study found that high IFIT3 expression could enhance the chemotherapeutic resistance of pancreatic ductal adenocarcinoma (PDAC) cells and was independently associated with the poor survival of PDAC patients [37].

Another previous study demonstrated that interferon-stimulated gene 15 (ISG15) downregulation could enhance cisplatin resistance via the DNA damage/repair pathway in A549/DDP cells [38]. However, ISG15 was found to be overexpressed in breast carcinoma, and ISG15 overexpression was associated with an unfavourable prognosis [39]. OAS1 and OASL belong to the 2,5-oligoadenylate synthetase (2-5OAS) gene family [40]. Many ISGs, including the OASL gene, have been tested for their antiviral specificity [41]. A previous study showed that OASL gene upregulation is involved in the inhibition of lung cancer cell proliferation and apoptosis [42]. Therefore, we propose that ISGs may play an important role in promoting GC resistance to trastuzumab, but the roles of these genes are controversial and need to be further investigated.

EBV plays an important role in gastric carcinogenesis [4, 5]. EBV latent membrane protein (LMP1) and latent membrane protein 2A (LMP2A) regulate the expression of the hub gene EGR1 [43, 44]. LMP2A has been found to suppress the expression of HER2 via the TWIST/YB-1 axis in EBV-associated gastric carcinoma [45]. Importantly, EBV and HER2 may exhibit crosstalk during human gastric carcinogenesis [46]. Interestingly, LMP1 can establish an antiviral state via the induction of ISGs, including OAS [47]. The EBV microRNA BART16 suppresses type I IFN signaling [48]. Thus, EBV may play a key role in gastric carcinogenesis, including in trastuzumab-acquired resistance, by regulating some of the identified hub genes or ISGs.

5. Conclusions

In summary, the current study identified genes related to acquired trastuzumab resistance in gastric cancer and analyzed their prognostic value. We found that four hub genes, including ERBB2, VIM, EGR1, and PSMB8, may participate in the development of chemoresistance to trastuzumab. However, the data in the present study were obtained by bioinformatics analysis, and the findings remain to be confirmed by further investigations. Therefore, more experiments are required to ascertain the clinical value of the identified genes as biomarkers and the underlying mechanism.

Data Availability

The data used to support the findings of this study are included within the article and the supplementary information files.

Conflicts of Interest

The authors declare no potential conflicts of interest.

Authors’ Contributions

Guangda Yang, Liumeng Jian, and Xiangan Lin contributed equally to this work.

Acknowledgments

We thank all the doctors in our department for their help with patient care so that we had sufficient time to conduct this study.

Supplementary Materials

Figure S1: the expression levels of ERBB2 (a), VIM (b), EGR1 (c), PSMB8 (d), IFI44 (e), IFI44L (f), IFIT2 (g), IFIT3 (h), ISG15 (i), OAS1 (j), OASL (k), SAMD9 (l), BST2 (m), IFI27 (n), IFIT1 (o), IFITM3 (p), MX1 (q), and OAS2 (r) in gastric cancer (UALCAN database). Red represents primary gastric cancer and blue represents normal samples. Table S1: survival analyses of the hub genes in gastric cancer with different parameters. (Supplementary Materials)

References

F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: A Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.

C. M. den Hoed and E. J. Kuipers, “Gastric cancer: how can we reduce the incidence of this disease?” Current Gastroenterology Reports, vol. 18, no. 7, p. 34, 2016.

L. Gao, A. Michel, M. N. Weck, V. Arndt, M. Pawlita, and H. Brenner, “Helicobacter pylori infection and gastric cancer risk: evaluation of 15 H. pylori proteins determined by novel multiplex serology,” Cancer Research, vol. 69, no. 15, pp. 6164–6170, 2009.

S. Akiba, C. Koriyama, R. Herrera-Goepfert, and Y. Eizuru, “Epstein-Barr virus associated gastric carcinoma: epidemiological and clinicopathological features,” Cancer Science, vol. 99, no. 2, pp. 195–201, 2008.

X. Chen, H. Chen, F. A. Castro, J. K. Hu, and H. Brenner, “Epstein–Barr virus infection and gastric cancer: a systematic review,” Medicine, vol. 94, no. 20, article e792, 2015.

N. Baniak, J.-L. Senger, S. Ahmed, S. C. Kanthan, and R. Kanthan, “Gastric biomarkers: a global review,” World Journal of Surgical Oncology, vol. 14, no. 1, p. 212, 2016.

N. Boku, “HER2-positive gastric cancer,” Gastric Cancer, vol. 17, no. 1, pp. 1–12, 2014.

Y. Geng, X. Chen, J. Qiu et al., “Human epidermal growth factor receptor-2 expression in primary and metastatic gastric cancer,” International Journal of Clinical Oncology, vol. 19, no. 2, pp. 303–311, 2014.

M. Bayrak, O. F. Olmez, E. Kurt et al., “Prognostic significance of c-erbB2 overexpression in patients with metastatic gastric cancer,” Clinical and Translational Oncology, vol. 15, no. 4, pp. 307–312, 2013.

Y. J. Bang, E. van Cutsem, A. Feyereislova et al., “Trastuzumab in combination with chemotherapy versus chemotherapy alone for treatment of HER2-positive advanced gastric or gastro-oesophageal junction cancer (ToGA): a phase 3, open-label, randomised controlled trial,” The Lancet, vol. 376, no. 9742, pp. 687–697, 2010.

R. Nahta and F. J. Esteva, “HER2 therapy: molecular mechanisms of trastuzumab resistance,” Breast Cancer Research, vol. 8, no. 6, p. 215, 2006.

G. Piro, C. Carbone, I. Cataldo et al., “An FGFR3 autocrine loop sustains acquired resistance to trastuzumab in gastric Cancer patients,” Clinical Cancer Research, vol. 22, no. 24, pp. 6164–6175, 2016.

D. S. Chandrashekar, B. Bashel, S. Balasubramanya et al., “UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses,” Neoplasia, vol. 19, no. 8, pp. 649–658, 2017.

Á. Nagy, A. Lánczky, O. Menyhárt, and B. Győrffy, “Validation of miRNA prognostic power in hepatocellular carcinoma using expression data of independent datasets,” Scientific Reports, vol. 8, no. 1, article 9227, 2018.

Y. Y. Janjigian, D. Werner, C. Pauligk et al., “Prognosis of metastatic gastric and gastroesophageal junction cancer by HER2 status: a European and USA international collaborative analysis,” Annals of Oncology, vol. 23, no. 10, pp. 2656–2662, 2012.

J. Matsubara, Y. Yamada, Y. Hirashima et al., “Impact of insulin-like growth factor type 1 receptor, epidermal growth factor receptor, and HER2 expressions on outcomes of patients with gastric cancer,” Clinical Cancer Research, vol. 14, no. 10, pp. 3022–3029, 2008.

F. Pietrantonio, M. Caporale, F. Morano et al., “HER2 loss in HER2-positive gastric or gastroesophageal cancer after trastuzumab therapy: implication for further clinical research,” International Journal of Cancer, vol. 139, no. 12, pp. 2859–2864, 2016.

A. M. Richardson, L. S. Havel, A. E. Koyen et al., “Vimentin is required for lung adenocarcinoma metastasis via heterotypic tumor cell-cancer-associated fibroblast interactions during collective invasion,” Clinical Cancer Research, vol. 24, no. 2, pp. 420–432, 2018.

M. Iwatsuki, K. Mimori, T. Fukagawa et al., “The clinical significance of vimentin-expressing gastric cancer cells in bone marrow,” Annals of Surgical Oncology, vol. 17, no. 9, pp. 2526–2533, 2010.

T. Utsunomiya, T. Yao, K. Masuda, and M. Tsuneyoshi, “Vimentin-positive adenocarcinomas of the stomach: co-expression of vimentin and cytokeratin,” Histopathology, vol. 29, no. 6, pp. 507–516, 1996.

C. H. Kwon, H. J. Park, Y. R. Choi et al., “PSMB8 and PBK as potential gastric cancer subtype-specific biomarkers associated with prognosis,” Oncotarget, vol. 7, no. 16, pp. 21454–21468, 2016.

T. Sun, H. Tian, Y.-G. Feng, Y.-Q. Zhu, and W.-Q. Zhang, “Egr-1 promotes cell proliferation and invasion by increasing β-catenin expression in gastric cancer,” Digestive Diseases and Sciences, vol. 58, 2012.

E. Myung, Y. L. Park, N. Kim et al., “Expression of early growth response-1 in human gastric cancer and its relationship with tumor cell behaviors and prognosis,” Pathology - Research and Practice, vol. 209, no. 11, pp. 692–699, 2013.

S. R. Johnston and A. Leary, “Lapatinib: a novel EGFR/HER2 tyrosine kinase inhibitor for cancer,” Drugs Today, vol. 42, no. 7, pp. 441–453, 2006.

I. Treish, R. Schwartz, and C. Lindley, “Pharmacology and therapeutic use of trastuzumab in breast cancer,” American Journal of Health-System Pharmacy, vol. 57, no. 22, pp. 2063–2076, 2000.

M. C. Franklin, K. D. Carey, F. F. Vajdos, D. J. Leahy, A. M. de Vos, and M. X. Sliwkowski, “Insights into ErbB signaling from the structure of the ErbB2-pertuzumab complex,” Cancer Cell, vol. 5, no. 4, pp. 317–328, 2004.

D. J. Kuhn, Q. Chen, P. M. Voorhees et al., “Potent activity of carfilzomib, a novel, irreversible inhibitor of the ubiquitin-proteasome pathway, against preclinical models of multiple myeloma,” Blood, vol. 110, no. 9, pp. 3281–3290, 2007.

L. Mi, B. L. Hood, N. A. Stewart et al., “Identification of potential protein targets of isothiocyanates by proteomics,” Chemical Research in Toxicology, vol. 24, no. 10, pp. 1735–1743, 2011.

P. Hertzog, S. Forster, and S. Samarajiwa, “Systems biology of interferon responses,” Journal of Interferon & Cytokine Research, vol. 31, no. 1, pp. 5–11, 2011.

I. C. McDowell, T. H. Modak, C. E. Lane, and M. Gomez-Chiarri, “Multi-species protein similarity clustering reveals novel expanded immune gene families in the eastern oyster Crassostrea virginica,” Fish & Shellfish Immunology, vol. 53, pp. 13–23, 2016.

W.-C. Huang, S.-L. Tung, Y.-L. Chen, P.-M. Chen, and P.-Y. Chu, “IFI44L is a novel tumor suppressor in human hepatocellular carcinoma affecting cancer stemness, metastasis, and drug resistance via regulating met/Src signaling pathway,” BMC Cancer, vol. 18, no. 1, p. 609, 2018.

X. Zhou, J. J. Michal, L. Zhang et al., “Interferon induced IFIT family genes in host antiviral defense,” International Journal of Biological Sciences, vol. 9, no. 2, pp. 200–208, 2013.

C. Wacher, M. Muller, M. J. Hofer et al., “Coordinated regulation and widespread cellular expression of interferon-stimulated genes (ISG) ISG-49, ISG-54, and ISG-56 in the central nervous system after infection with distinct viruses,” Journal of Virology, vol. 81, no. 2, pp. 860–871, 2006.

K. C. Lai, C. J. Liu, K. W. Chang, and T. C. Lee, “Depleting IFIT2 mediates atypical PKC signaling to enhance the migration and metastatic activity of oral squamous cell carcinoma cells,” Oncogene, vol. 32, no. 32, pp. 3686–3697, 2013.

K. C. Lai, K. W. Chang, C. J. Liu, S. Y. Kao, and T. C. Lee, “IFN-induced protein with tetratricopeptide repeats 2 inhibits migration activity and increases survival of oral squamous cell carcinoma,” Molecular Cancer Research, vol. 6, no. 9, pp. 1431–1439, 2008.

L. Chen, W. Zhai, X. Zheng et al., “Decreased IFIT2 expression promotes gastric cancer progression and predicts poor prognosis of the patients,” Cellular Physiology and Biochemistry, vol. 45, no. 1, pp. 15–25, 2018.

Y. Zhao, A. Altendorf-Hofmann, I. Pozios et al., “Elevated interferon-induced protein with tetratricopeptide repeats 3 (IFIT3) is a poor prognostic marker in pancreatic ductal adenocarcinoma,” Journal of Cancer Research and Clinical Oncology, vol. 143, no. 6, article 2351, pp. 1061–1068, 2017.

Y. Huo, Z. Zong, Q. Wang, Z. Zhang, and H. Deng, “ISG15 silencing increases cisplatin resistance via activating p53-mediated cell DNA repair,” Oncotarget, vol. 8, no. 64, pp. 107452–107461, 2017.

N. Bektas, E. Noetzel, J. Veeck et al., “The ubiquitin-like molecule interferon-stimulated gene 15 (ISG15) is a potential prognostic marker in human breast cancer,” Breast Cancer Research, vol. 10, no. 4, article R58, 2008.

S. Kakuta, S. Shibata, and Y. Iwakura, “Genomic structure of the mouse 2,5-oligoadenylate synthetase gene family,” Journal of Interferon & Cytokine Research, vol. 22, no. 9, pp. 981–993, 2002.

J. W. Schoggins, S. J. Wilson, M. Panis et al., “A diverse range of gene products are effectors of the type I interferon antiviral response,” Nature, vol. 472, no. 7344, pp. 481–485, 2011.

J. Lv, L. Wang, H. Shen, and X. Wang, “Regulatory roles of OASL in lung cancer cell sensitivity to Actinidia chinensis planch root extract (acRoots),” Cell Biology and Toxicology, vol. 34, no. 3, pp. 207–218, 2018.

M. Vockerodt, W. Wei, E. Nagy et al., “Suppression of the LMP2A target gene, EGR-1, protects Hodgkin’s lymphoma cells from entry to the EBV lytic cycle,” The Journal of Pathology, vol. 230, no. 4, pp. 399–409, 2013.

J. H. Kim, W. S. Kim, J. H. Kang, H. Y. Lim, Y. H. Ko, and C. Park, “Egr-1, a new downstream molecule of Epstein-Barr virus latent membrane protein 1,” FEBS Letters, vol. 581, no. 4, pp. 623–628, 2007.

Y. Zhang, X. Zhao, C. Tan et al., “Epstein-Barr virus latent membrane protein 2A suppresses the expression of HER2 via a pathway involving TWIST and YB-1 in Epstein-Barr virus-associated gastric carcinomas,” Oncotarget, vol. 6, no. 1, pp. 207–220, 2015.

F. S. Cyprian, N. al-Antary, and A.-E. al Moustafa, “HER-2/Epstein-Barr virus crosstalk in human gastric carcinogenesis: a novel concept of oncogene/oncovirus interaction,” Cell Adhesion & Migration, vol. 12, no. 1, pp. 1–4, 2017.

J. Zhang, S. C. Das, C. Kotalik, A. K. Pattnaik, and L. Zhang, “The latent membrane protein 1 of Epstein-Barr virus establishes an antiviral state via induction of interferon-stimulated genes,” The Journal of Biological Chemistry, vol. 279, no. 44, pp. 46335–46342, 2004.

M. J. G. Hooykaas, M. van Gent, J. A. Soppe et al., “EBV microRNA BART16 suppresses type I IFN signaling,” Journal of Immunology, vol. 198, no. 10, pp. 4062–4073, 2017.