Skip to main content

Advertisement

Prediction of key genes and pathways involved in trastuzumab-resistant gastric cancer

Abstract

Background

Trastuzumab has been prevailingly accepted as a beneficial treatment for gastric cancer (GC) by targeting human epidermal growth factor receptor 2 (HER2)-positive. However, the therapeutic resistance of trastuzumab remains a major obstacle, restricting the therapeutic efficacy. Therefore, identifying potential key genes and pathways is crucial to maximize the overall clinical benefits.

Methods

The gene expression profile GSE77346 was retrieved to identify the differentially expressed genes (DEGs) associated with the trastuzumab resistance in GC. Next, the DEGs were annotated by the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The DEGs-coded protein-protein interaction (PPI) networks and the prognostic values of the 20 hub genes were determined. Correlation of the hub genes were analyzed in The Cancer Genome Atlas. The prognostic values of hub genes were further validated by Kaplan-Meier (KM) plotter.

Results

A total of 849 DEGs were identified, with 374 in upregulation and 475 in downregulation. Epithelium development was the most significantly enriched term in biological processes while membrane-bounded vesicle was in cellular compartments and cell adhesion molecular binding was in molecular functions. Pathways in cancer and ECM-receptor interaction were the most significantly enriched for all DEGs. Among the PPI networks, 20 hub genes were defined, including CD44 molecule (CD44), HER-2, and cadherin 1 (CDH1). Six hub genes were associated with favorable OS while eight were associated with poor OS. Mechanistically, 2′-5′-oligoadenylate synthetase 1, 3 (OAS1, OAS3) and CDH1 featured high degrees and strong correlations with other hub genes.

Conclusions

This bioinformatics analysis identified key genes and pathways for potential targets and survival predictors for trastuzumab treatment in GC.

Background

Gastric cancer (GC) remains one of the leading common causes for cancer-related mortality and major global heath challenges [1,2,3,4]. Despite the incidence declining in industrialized nations, most new cases are occurred in South America, East Asia, and Eastern Europe [2, 5]. Surgery is the primary treatment for resectable GC [6]. However, the dissection extent of lymph node (D1, D2) remains controversial [3]. Kang et al. reported 46.5% patients who underwent curative surgery experienced recurrence, and half of the recurrence occurred in less than 3 years [7]. In the Dutch Gastric Cancer Group (DGCG) trial, 65% curative resected patients experienced recurrence with 30% overall survival (OS) for D1 and 35% for D2 [8]. Consistently, the Medical Research Council (MRC) trial reported a 34% 5-year OS [9]. Noteworthy, the inclusion of targeted drugs, such as angiogenesis inhibitors (ramucirumab) and epidermal growth factor receptor (EGFR) antibodies (nimotuzumab), have shown encouraging therapeutic benefits in GC patients [10, 11].

Trastuzumab, a monoclonal antibody targeting epidermal growth factor receptor 2 (HER2) in breast cancer [12], was also among the promising therapeutic management to the GC patients with HER2-positive [13, 14]. It eliminated the activity of HER2 receptor and weakened subsequent multiple signaling pathways [15]. The first randomized prospect trial had shown that a triplet regimen of trastuzumab, cisplatin, and a fluoropyrimidine significantly improved the median OS of GC with HER2 overexpression or amplification [13]. In fact, secondary resistance was acquired within a median of two therapeutic cycles [16]. Until now, the resistance to trastuzumab in GC remains a major obstacle with limited clinical benefits. Efficient biomarkers and underlying mechanism are yet to be fully elucidated.

Hereby, potential biomarkers and pathways associated with trastuzumab resistance were investigated in GC cell lines by the gene expression profile, GSE77346 [17], from the Genetic Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/). The prognostic values of the biomarkers and potential mechanisms were assessed.

Methods

Gene expression profile from GEO database

The gene expression profile, GSE77346, was retrieved from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/) [18]. The profile was generated by GPL10558, Illumina Human 48 K gene chips (Illumina HumanHT-12 V4.0 Expression BeadChip). The GSE77346 dataset consisted of one trastuzumab-sensitive NCI-N87 cell line and four trastuzumab-resistant cell lines (N87-TR1, N87-TR2, N87-TR3, N87-TR4). Briefly, all the cell lines were maintained in Roswell Park Memorial Institute (RPMI) 1640 medium with 10% heat-inactivated FBS. The green fluorescent protein (GFP) +/luciferase+ NCI-N87 cell lines were harvested and injected into the gastric walls of a nude mice. The tumor-bearing mice were received 20 mg/kg trastuzumab i.p. twice per week when the resulting tumors were detectable (Living Image Software program, Xenogen). The trastuzumab treatments were stopped when the tumors were relapsed. By repeated GFP flow cytometric sorting (FACSAria II sorter, Becton Dickinson), four trastuzumab-resistant cell lines were established [17]. Next, total RNA was retrieved by TRIzol reagent (Ambion, Warrington, UK). The synthesis of biotinylated cRNA (Illumina TotalPrep RNA Amplification Kit, Ambion) and the hybridization (Human HT-12 V4 BeadChip) were performed according to the manufacturer protocols. Probe intensity was obtained and normalized by the Illumina GenomeStudio software (Genome Studio V2011.1) [17]. The gene expression profiles GSE13861, including 84 samples (65 tumors and 19 normal tissues), were used for investigation of mRNAs expression of the hub genes between tumor and normal tissues (Illumina Human V3) [19]. For external validation on gene expression profiles with other target drugs, we further included GSE19043 and GSE95414. GSE19043 contained 21 samples from DiFi and GTL-16 cell lines, of which biological triplicates of DiFi cells with gefitinib (EGFR inhibition) and DMSO (control) were used in this study for validation. The platform was GPL5104, Sentrix HumanRef-8 v2 Expression BeadChip [20]. GSE95414 contained one parental NCI-N87 cell line and one trastuzumab-DM1 (T-DM1, trastuzumab emtansine)-resistant cell line. T-DM1 is designed to achieve a combinational therapy of trastuzumab and DM1 (a potent microtubule-disrupting drug, a maytansine derivative) [21]. The RNA was processed by Human Transcriptome Array 2.0 arrays (Affymetrix, GPL17586). Given the absence of biological replicates, the fold change between the T-DM1-resistant cell line and parental cell line was used for investigation (original study of GSE95414 is not yet published).

Data processing on DEGs

The differentially expressed genes (DEGs) between the trastuzumab-resistant cell lines and sensitive control were identified by the GEO2R analytical tool [22]. Benjamini and Hochberg method was used for false discovery rate (FDR). The cut-off values of DEGs were defined as adj.p value < 0.05 and log2 fold change (log FC) > 2 or < − 2. The DEG expression data were processed for a bidirectional hierarchical clustering plot (FunRich, http://www.funrich.org) [23].

Gene ontology and pathway analysis of DEGs

The Database for Annotation, Visualization, and Integrated Discovery (DAVID, http://david.abcc.ncifcrf.gov/) was employed for the gene ontology (GO) consortium reference, including biological processes (BP), cellular components (CC), and molecular functions (MF) [24, 25]. In addition, DAVID was also employed for pathway enrichment annotations with the data resources from Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/) pathway enrichment analysis [24, 26].

Protein-protein interaction (PPI) networks and module analysis

The interaction networks of the DEG-coded proteins were determined by the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING, http://www.string-db.org/) [27]. Node degree ≥ 5 was defined as the cut-off values for further PPI networks visualization by Cytoscape software (version 3.6.0; http://www.cytoscape.org/) [28]. The Molecular Complex Detection (MCODE) program embedded in Cytoscape was used to subcluster the PPI networks with predefined cutoff criterions (max. depth = 100, node score = 0.2 and k-score = 2) [29]. Hub genes were defined by the degree value (paired connections between each node). In addition, the betweenness centrality (defining the fraction of shortest paths involved in a given node) of the hub genes were also added.

Survival analysis of the hub genes

Kaplan-Meier (KM) plotter enables comprehensive analysis of the prognostic values among lists of genes in various cancers based on multiple genomic profiles, including GSE14210, GSE15459, GSE22377, GSE29272, GSE51105, and GSE62254 [30]. The prognostic values of overall survivals (OS) for hub genes were displayed with the hazard ratios (HR) and log-rank p values.

Hub genes correlation in TCGA

The gene expression profiling interactive analysis (GEPIA, http://gepia.cancer-pku.cn) was established for customized genomic analysis based on The Cancer Genome Atlas (TCGA) database [31]. The top 20 hub genes were extracted for interactive networks based on paired gene correlations of the stomach adenocarcinoma (STAD) cohort in TCGA (Pearson correlation coefficients). In addition, the mRNA expressions of the hub genes were also investigated between tumor and normal tissues.

Moreover, the stage-specific expression of each hub gene was also generated by GEPIA. The mRNA expressions of the hub genes of TCGA (STAD) were also retrieved from the Xena system, University of California, Santa Cruz (UCSC) for prognostic analysis [32].

Statistical analysis

Generally, p value < 0.05 was defined as cut-off criterion and considered statistically significant in all cases. SPSS 17.0 (Chicago, IL, USA) and Prism 5.0 (GraphPad Software, San Diego, CA) were used for statistical analysis and illustration.

Results

Identification of DEGs and heat map clustering

A total of 849 DEGs were identified to be associated with trastuzumab resistance, with 374 genes upregulated and 475 downregulated (Fig. 1). A bidirectional hierarchical clustering heat map of the DEGs was illustrated (Fig. 2).

Fig. 1
figure1

Volcano plot of the differentially expressed genes (DEGs) involved in trastuzumab-resistant gastric cancer (GC) with respect to control. The negative log10-adjusted p values (y-axis) were plotted against log2 fold change (log2FC) (x-axis). DEGs were identified by GEO2R. The threshold for significance was|log2FC| > 2 and adjusted p value < 0.05. Red, upregulated DEGs; green, downregulated DEGs

Fig. 2
figure2

Heat map for the DEGs in trastuzumab-resistant GC cell lines. The bidirectional hierarchical clustering heat map was generated by FunRich software. The expression values were all processed by log2 fold change in prior to the heat map construction. Blue represents downregulation; red represents upregulation

GO enrichment analysis

The GO enrichment analysis was conducted by the DAVID tool. A total of 193 BP terms significantly enriched, including epithelium development/cell surface receptor signaling pathway/locomotion (Table 1). A total of 23 CC terms were significantly enriched, including membrane-bounded vesicle/extracellular region part/extracellular vesicle (Table 1). A total of nine MF terms were significantly enriched, including top-ranked cell adhesion molecular binding/glycoprotein binding/growth factor binding (Table 1). Specifically, in each term, top ranked 10 most significantly enriched gene-ontologies of upregulated and downregulated DEGs were compared (Fig. 3). In BP term, nervous system development and response to type I interferon were significantly enriched in up/downregulated DEGs, respectively (Fig. 3a). In CC term, proteinaceous extracellular matrix and extracellular region part were significantly enriched in up/down regulated DEGs, respectively (Fig. 3b). In MF term, protein dimerization activity and cell adhesion molecule binding were significantly enriched in up/down regulated DEGs, respectively (Fig. 3c).

Table 1 Gene ontology analysis of the DEGs
Fig. 3
figure3

Gene ontology (GO) enrichment of the DEGs involved in trastuzumab resistance. a Biological function (BF) enrichment in up/downregulated DEGs. b Cellular component (CC) enrichment in up/downregulated DEGs. c Molecular function enrichment in up/downregulated DEGs

KEGG pathways analysis

Noteworthy, only two significant signaling pathways were identified in KEGG pathway analysis with cut-off values (p < 0.05, FDR < 0.05): pathways in cancer (hsa05200) and ECM-receptor interaction (hsa04512) (Table 2). The top ten enriched signaling pathways in upregulated and downregulated DEGs were illustrated, respectively (Fig. 4). Of note, no significant pathway was identified in upregulated set, and only one, the pathways in cancer (hsa5200), was identified as significant in downregulated set.

Table 2 KEGG pathway enrichment analysis
Fig. 4
figure4

Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the DEGs involved in trastuzumab resistance. a KEGG pathways in upregulated DEGs. b KEGG pathways in downregulated DEGs

PPI network and modules

Next, the PPI networks were initially obtained by the STRING database and visualized by Cytoscape with degrees of each nodes ≥ 5. A total of 291 nodes 1883 edges were included in the PPI networks (Fig. 5). The top 20 hub genes with highest degrees were determined, including CD44 molecule (CD44), erb-b2 receptor tyrosine kinase 2 (HER2), cadherin 1 (CDH1), 2′-5′-oligoadenylate synthetase 1–3 (OAS1–3), 2′-5′-oligoadenylate synthetase-like (OASL), ISG15 ubiquitin-like modifier (ISG15), bone morphogenetic protein 4 (BMP4), signal transducer and activator of transcription 1 (STAT1), early growth response 1 (EGR1), cyclin D1 (CCND1), vimentin (VIM), Wnt family member 5A (WNT5A), KIT proto-oncogene receptor tyrosine kinase (KIT), bone morphogenetic protein 2 (BMP2), interferon regulatory factor 9 (IRF9), MX dynamin-like GTPase 1 (MX1), FYN proto-oncogene, Src family tyrosine kinase (FYN), and HECT and RLD domain containing E3 ubiquitin protein ligase family member 6 (HERC6) (Fig. 5, Table 3). In addition, the top scored three modules were determined by MCODE in Cytoscape, with KEGG enrichment results (Fig. 6). Furthermore, the siRNAs of the hub genes were summarized (Additional file 1: Table S1) [33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51].

Fig. 5
figure5

Protein-protein interaction (PPI) networks of the DEGs. Red nodes represented upregulated genes and blue nodes represented downregulated genes (nodes ≥ 5). The interaction between genes was illustrated by lines

Table 3 Hub genes in the PPI networks
Fig. 6
figure6

The most scored three modules with KEGG enrichment results. a Module-1. b KEGG analysis of module 1. c Module 2. d KEGG analysis of module 2. e Module 3. f KEGG analysis of module 3. Red nodes represented upregulated genes while blue nodes represented downregulated genes

Prognostic analysis and mRNA expression of hub genes

The prognostic values of the hub genes were assessed by the KM plotter in GC. High HER2, CDH1, OAS1, OAS3, ISG15, BMP4, CCND1, and WNT5A expression levels were associated with poor OS, whereas high CD44, STAT1, EGR1, VIM, KIT, and FYN expression levels were associated with favorable OS. OAS2, OASL, BMP2, IRF9, MX1, and HERC6 were not significantly associated with OS (Fig. 7). The mRNAs expression of CD44, HER2, CDH1, OAS1, OAS2, OAS3, OASL, ISG15, STAT1, CCND1, and WNT5A were significantly upregulated in tumor while only KIT was significantly downregulated in tumor (TCGA) compared to normal (TCGA normal + GTEx normal) (Fig. 8a). Next, we further compared the mRNA expression of the hub genes between tumor (TCGA) and normal (TCGA) by the data retrieved from the Xena system. In fact, the results from the Xena (TCGA tumor vs TCGA normal) were different from GEPIA (TCGA tumor vs TCGA normal + GTEx normal). Only five hub genes (STAT1, OAS3, OAS2, CDH1, ISG15) significantly exhibited upregulation and two hub genes (KIT and EGR1) exhibited downregulation according to the thresholds (adj.p value < 0.05 and |logFC| > 1) (Additional file 2: Table S2). Interestingly, the gene with the most significant logFC value is KIT (logFC = − 2.11514), whereas the gene with the most significant adj.p value is STAT1 (adj.p value = 8.99E−12).

Fig. 7
figure7

Survival plots of the prognostic values (overall survival) of hub genes involved in trastuzumab-resistant GC. The survival values of the hub genes were generated by the Kaplan-Meier (KM) plotter. The expressions of hub genes were dichotomized by optimal cutoff values. Patients number = 593. p values were calculated by log rank method

Fig. 8
figure8

The mRNA expression and gene-gene correlation of the hub genes associated with trastuzumab-resistant GC. a The mRNAs expression of hub genes in tumor and normal tissues in TCGA, red: tumor, blue: normal. b the STAD of TCGA was calculated with Pearson’s correlation coefficient (− 1 to 1). Red line: negative correlation, black line: positive correlation. Wider line indicated higher correlation value (p value < 0.05)

Furthermore, the mRNA expression of the hub genes (IRF9 was not available) was externally validated in GSE13861 (Additional file 3: Figure S1). Consistently, CD44, OAS3, ISG15, STAT1, and WNT5A were significantly upregulated whereas KIT was significantly downregulated in tumor compared to normal in GSE13861. Moreover, BMP4 was significantly upregulated in tumor in GSE13861. OASL, EGR1, and BMP2 were significantly downregulated in tumor in GSE13861 (Additional file 3: Figure S1). Moreover, the mRNA expression of all the hub genes in specific clinic stages had been analyzed. In fact, only CD44 (p = 0.0146), VIM (p = 1.07e−05) and KIT (0.00759) exhibited significant stage-specific expression (Additional file 4: Figure S2).

Mechanism of hub genes correlations associated with trastuzumab resistance

To further elucidate the underlying mechanism between the DEGs, the STAD of TCGA data was employed based on GEPIA platform. Of note, 87.8% (65/74) gene-gene correlations were positive. What is more, OAS1, 3, and CDH1 featured high degrees and strong correlations with other hub genes. Additionally, VIM was negatively correlated with CCND1, HER2, and CDH1, respectively. KIT was negatively correlated with HER2, ISG15, and OAS1, respectively (Fig. 8b). Meanwhile, to investigate the potential roles of the hub genes in other target therapies, GSE19043 and GSE95414 were retrieved for external investigation (Additional file 5: Table S3). In GSE19043, none of the hub genes exhibited differential expression between gefitinib group and control, whereas in GSE95414, only six of the hub genes, including VIM, BMP2, CD44, OAS3, KIT, and WNT5A, showed slight fold change values > 1 between T-DM1-resistant cell lines and control (Additional file 6: Figure S3). In summary, the hub genes identified in this study may not be directly involved in gefitinib (EGFR inhibition, GSE95414) and T-DM1 (GSE19043) (Additional file 6: Figure S3).

Discussion

Although the overall mortality and morbidity of GC has been declining over the decades around the globe, it is one of the most common causes for cancer-related deaths. Postoperative recurrence remains high even with curable resection and combinational chemotherapy [7,8,9]. Trastuzumab, the only approved treatment for GC with HER2 overexpress, had contributed to the encouraging results in GC clinical trials [13, 14]. However, secondary resistance of trastuzumab remained one of the major challenges in treatment courses. Therefore, identification of potential mechanisms and key genes underlying the acquired trastuzumab resistance could distinguish the sensitive subsets and improve overall benefits.

Generally, individual gene rarely dictate either systematic biochemical physiological actions or sophisticated multilevel network interactions. Up to now, genomic data had been stored in large matrix and processed by well-established bioinformatics pipelines for the ultimate conclusive visualization.

This study provided a systematic bioinformatics analysis of the gene expression profile, GSE77346, containing four trastuzumab-resistant cell lines and one sensitive cell line. Pathways in cancer and ECM-receptor interaction were the most significantly enriched for all DEGs. CD44, STAT1, EGR1, VIM, KIT, and FYN were associated with favorable OS while HER2, CDH1, OAS1, OAS3, ISG15, BMP4, CCND1, and WNT5A were associated with poor OS.

Mechanistically, OAS1, OAS3, and CDH1 featured highest degrees among the hub genes, diverse from the nodes (CD44, HER2, and CDH1) with highest degrees in PPI networks.

OAS1 and OAS3, which encode the key enzymes, 2′, 5′-oligoadenylate synthetase (2′5′AS), are involved in viral genome degradation and inhibits protein synthesis [52, 53]. As classic interferon target genes, OAS1 and OAS3 differ in cellular compartment, conformation, and biological functions [54]. Previously, OAS1 and OAS3 had been participated in apoptosis process [55]. Until now, only OAS3 had been associated with the HPV persistence and progression of cervical cancer [56]. No specific study unveiled the association between OAS1 and OAS3 and GC. This is the first in silico study suggesting the involvement of OAS1and OAS3 in trastuzumab-resistant GC.

CD44, a key cancer stem cell (CSC) marker, was downregulated in trastuzumab-resistant breast cancer and associated with the trastuzumab resistance in GC. [57]. Previously, high expression of CD44 correlated with downregulated HER2 in breast cancer cell lines [58]. SiRNA CD44 led to reduced internalization of trastuzumab, highlighting the involvement of endocytosis and membrane trafficking [58]. Furthermore, Bao et al. revealed that CD44 could directly bind to HER2 and increase invasiveness both in vivo and vitro [59]. Consistently, this study highlighted CD44 as the top hub gene in PPI networks of trastuzumab-resistant GC; however, the correlation between CD44 and HER2 associated with trastuzumab resistance in GC required further validation.

Noteworthy, eight of the 20 hub genes (WNT5A, BMP4, BMP2, CCND1, HER2, CDH1, KIT, STAT1) associated with trastuzumab resistance were commonly enriched in the pathways in cancer (KEGG hsa05200). Thus, the acquired resistance of trastuzumab in GC at least could be partially attributed by the progression of GC itself, if not all. Moreover, the potential impact of the mutations and fusion of the genes in the pathway in cancer on the trastuzumab resistance in GC remains largely unsolved.

In addition, for PPI networks, both degree and betweenness centrality were included for proper evaluation of hub genes. Generally, centrality is not generally equivalent to connectivity. As a local quantity, connectivity does not fully elucidate the importance of certain node in PPI networks. Thus, both connectivity and betweenness centrality were incorporated for a good measurement of hub genes in PPI networks [60].

Remarkably, ion channels, one of the major transmembrane complexes that regulate the communication between the extracellular matrix and intracellular environments, can influence the growth and invasiveness of cancer cells by altered expression or biological activities [61, 62]. In fact, ion channels could be novel molecular targets [62]. Fujimoto et al. indicated that the inhibition of ANO1, a Ca2 + -activated Cl- channel overexpressed in HER2-positive breast cancer, could lead to the transcriptional repression of HER2 in breast cancer cells with resistance to trastuzumab [63]. Another Ca2 + -permeable channel, transient receptor potential canonical 6 (TRPC6), exhibited a vital role in tumor growth, differentiation, and apoptosis with promising pharmaceutic target values [64, 65].

Recently, Huang et al. published a result focusing on the trastuzumab-resistant role of COL4A1 in GC [66]. Validation of COL4A1 in GSE77346 was one of the key steps in their study. However, GSE77346 remained far from fully explored with respect to trastuzumab resistance. In fact, new agents to be discovered against HER2 and other signaling pathways open the way to the improvement of trastuzumab therapy [67].

In breast cancer, trastuzumab remains one of the intensively studied drugs. It has been recommended as combination treatments in breast cancer [67]. In fact, mining the relationships between HER2 signaling pathway and other signaling pathways as well as the potential mechanisms provides greater insights for rational combination therapy. Currently, targets such as mTOR, PI3K, IGF-1R, Akt, HSP90, and VEGF exhibited significant clinical interests in HER2-positive breast cancer [67]. However, insightful evidences to define, refine, and optimize the use of trastuzumab in gastric cancer patients with HER2-positive remain largely lacked. Therefore, this study contributed to the understanding of trastuzumab resistance and the prognostic values of hub genes and opened the way for future research in combination therapy in gastric cancer.

Noteworthy, this was the first in silico study focusing on the bioinformatics analysis of trastuzumab resistance in GC, predicting the key genes and pathways associated with trastuzumab resistance. In addition, this study also investigated the prognostic values of key genes. However, no disease-free survival (DFS) or progression-free survival (PFS) was collected. Further clinical and experimental validation of the study findings was required.

Conclusion

This bioinformatics analysis identified key genes and pathways as potential targets and predictors associated with trastuzumab resistance GC and further opened the way to the improvement of trastuzumab therapy in GC.

Abbreviations

BP:

Biological processes

CC:

Cellular components

CI:

Confidential intervals

CSC:

Cancer stem cell

DAVID:

Database for Annotation, Visualization, and Integrated Discovery

DFS:

Disease-free survival

EGFR:

Epidermal growth factor receptor

FDR:

False discovery rate

GC:

Gastric cancer

GEO:

Gene Expression Omnibus

GFP:

Green fluorescent protein

GO:

Gene ontology

HER2:

Epidermal growth factor receptor 2

HR:

Hazard ratio

KEGG:

Kyoto Encyclopedia of Genes and Genomes

MCODE:

Molecular Complex Detection

MF:

Molecular functions

OS:

Overall survival

PFS:

Progression-free survival

PPI:

The protein-protein interaction

RPMI:

Roswell Park Memorial Institute

STAD:

Stomach Adenocarcinoma

STRING:

Search Tool for the Retrieval of Interacting Genes

TCGA:

The Cancer Genome Atlas

References

  1. 1.

    Siegel RL, Miller KD, Jemal A. Cancer statistics, 2016. CA Cancer J Clin. 2016;66(1):7–30.

  2. 2.

    Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, Jemal A, Yu XQ, He J. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66(2):115–32.

  3. 3.

    Ajani JA, Bentrem DJ, Besh S, D’Amico TA, Das P, Denlinger C, Fakih MG, Fuchs CS, Gerdes H, Glasgow RE, Hayman JA. Gastric cancer, version 2.2013. J Natl Compr Cancer Netw. 2013;11(5):531–46.

  4. 4.

    Guggenheim DE, Shah MA. Gastric cancer epidemiology and risk factors. J Surg Oncol. 2013;107(3):230–6.

  5. 5.

    Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. 2011;61(2):69–90.

  6. 6.

    Japanese Gastric Cancer Association. Japanese gastric cancer treatment guidelines 2014 (ver. 4). Gastric Cancer. 2017;20(1):1–9.

  7. 7.

    Kang WM, Meng QB, Yu JC, Ma ZQ, Li ZT. Factors associated with early recurrence after curative surgery for gastric cancer. World J Gastroenterol: WJG. 2015;21(19):5934.

  8. 8.

    Hartgrink HH, Van de Velde CJ, Putter H, Bonenkamp JJ, Klein Kranenbarg E, Songun I, Welvaart K, Van Krieken JH, Meijer S, Plukker JT, Van Elk PJ. Extended lymph node dissection for gastric cancer: who may benefit? Final results of the randomized Dutch gastric cancer group trial. J Clin Oncol. 2004;22(11):2069–77.

  9. 9.

    Cuschieri A, Weeden S, Fielding J, Bancewicz J, Craven J, Joypaul V, Sydes M. Patient survival after D 1 and D 2 resections for gastric cancer: long-term results of the MRC randomized surgical trial. Br J Cancer. 1999;79(9–10):1522.

  10. 10.

    Wilke H, Muro K, Van Cutsem E, Oh SC, Bodoky G, Shimada Y, Hironaka S, Sugimoto N, Lipatov O, Kim TY, Cunningham D. Ramucirumab plus paclitaxel versus placebo plus paclitaxel in patients with previously treated advanced gastric or gastro-oesophageal junction adenocarcinoma (RAINBOW): a double-blind, randomised phase 3 trial. Lancet Oncol. 2014;15(11):1224–35.

  11. 11.

    Satoh T, Lee KH, Rha SY, Sasaki Y, Park SH, Komatsu Y, Yasui H, Kim TY, Yamaguchi K, Fuse N, Yamada Y. Randomized phase II trial of nimotuzumab plus irinotecan versus irinotecan alone as second-line therapy for patients with advanced gastric cancer. Gastric Cancer. 2015;18(4):824–32.

  12. 12.

    Seidman A, Hudis C, Pierri MK, Shak S, Paton V, Ashby M, Murphy M, Stewart SJ, Keefe D. Cardiac dysfunction in the trastuzumab clinical trials experience. J Clin Oncol. 2002;20(5):1215–21.

  13. 13.

    Bang YJ, Van Cutsem E, Feyereislova A, Chung HC, Shen L, Sawaki A, Lordick F, Ohtsu A, Omuro Y, Satoh T, Aprile G. 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. Lancet. 2010;376(9742):687–97.

  14. 14.

    Lordick F, Janjigian YY. Clinical impact of tumour biology in the management of gastroesophageal cancer. Nat Rev Clin Oncol. 2016;13(6):348.

  15. 15.

    Xu W, Yang Z, Lu N. Molecular targeted therapy for the treatment of gastric cancer. J Exp Clin Cancer Res. 2016;35(1):1.

  16. 16.

    Okines AF, Cunningham D. Trastuzumab in gastric cancer. Eur J Cancer. 2010;46(11):1949–59.

  17. 17.

    Piro G, Carbone C, Cataldo I, Di Nicolantonio F, Giacopuzzi S, Aprile G, Simionato F, Boschi F, Zanotto M, Mina MM, Santoro R. An FGFR3 autocrine loop sustains acquired resistance to trastuzumab in gastric cancer patients. Clin Cancer Res. 2016;22(24):6164–75.

  18. 18.

    Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30(1):207–10.

  19. 19.

    Cho JY, et al. Gene expression signature–based prognostic risk score in gastric cancer. Clin Cancer Res. 2011;17(7):1850-57.

  20. 20.

    Bertotti A, Burbridge MF, Gastaldi S, Galimi F, Torti D, Medico E, Giordano S, Corso S, Rolland-Valognes G, Lockhart BP, Hickman JA. Only a subset of Met-activated pathways are required to sustain oncogene addiction. Sci Signal. 2009;2(100):ra80.

  21. 21.

    Junttila TT, Li G, Parsons K, Phillips GL, Sliwkowski MX. Trastuzumab-DM1 (T-DM1) retains all the mechanisms of action of trastuzumab and efficiently inhibits growth of lapatinib insensitive breast cancer. Breast Cancer Res Treat. 2011;128(2):347–56.

  22. 22.

    Davis S, Meltzer PS. GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor. Bioinformatics. 2007;23(14):1846–7.

  23. 23.

    Pathan M, Keerthikumar S, Ang CS, Gangoda L, Quek CY, Williamson NA, Mouradov D, Sieber OM, Simpson RJ, Salim A, Bacic A. FunRich: an open access standalone functional enrichment and interaction network analysis tool. Proteomics. 2015;15(15):2597–601.

  24. 24.

    Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008;4(1):44.

  25. 25.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25.

  26. 26.

    Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

  27. 27.

    Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, Simonovic M, Roth A, Santos A, Tsafou KP, Kuhn M. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2014;43(D1):D447–52.

  28. 28.

    Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

  29. 29.

    Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4(1):2.

  30. 30.

    Lánczky A, Nagy Á, Bottai G, Munkácsy G, Szabó A, Santarpia L, Győrffy B. miRpower: a web-tool to validate survival-associated miRNAs utilizing expression data from 2178 breast cancer patients. Breast Cancer Res Treat. 2016;160(3):439–46.

  31. 31.

    Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–102.

  32. 32.

    Goldman M, Craft B, Zhu J, Haussler D. The UCSC Xena system for cancer genomics data visualization and interpretation [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017. Washington, DC., Philadelphia: AACR; 2017. Cancer Res 2017;77(13 Suppl):Abstract nr 2584.

  33. 33.

    Subramaniam V, Vincent IR, Gilakjan M, Jothy S. Suppression of human colon cancer tumors in nude mice by siRNA CD44 gene therapy. Exp Mol Pathol. 2007;83(3):332–40.

  34. 34.

    Tan WB, Jiang S, Zhang Y. Quantum-dot based nanoparticles for targeted silencing of HER2/neu gene via RNA interference. Biomaterials. 2007;28(8):1565–71.

  35. 35.

    Herrero-Mendez A, Almeida A, Fernández E, Maestre C, Moncada S, Bolaños JP. The bioenergetic and antioxidant status of neurons is controlled by continuous degradation of a key glycolytic enzyme by APC/C–Cdh1. Nat Cell Biol. 2009;11(6):747.

  36. 36.

    Zhao J, Feng N, Li Z, Wang P, Qi Z, Liang W, Zhou X, Xu X, Liu B. 2ng P, oligoadenylate synthetase 1 (OAS1) inhibits PRRSV replication in Marc-145 cells. Antivir Res. 2016;132:268–73.

  37. 37.

    Bin L, Howell MD, Kim BE, Streib JE, Hall CF, Leung DY. Specificity protein 1 is pivotal in the skin’s antiviral response. J Allergy Clin Immunol. 2011;127(2):430–8.

  38. 38.

    Lin W, Zhu C, Hong J, Zhao L, Jilg N, Fusco DN, Schaefer EA, Brisac C, Liu X, Peng LF, Xu Q. The spliceosome factor SART1 exerts its anti-HCV action through mRNA splicing. J Hepatol. 2015;62(5):1024–32.

  39. 39.

    Zheng S, Zhu D, Lian X, Liu W, Cao R, Chen P. Porcine 2ction-oligoadenylate synthetases inhibit Japanese encephalitis virus replication in vitro. J Med Virol. 2016;88(5):760–8.

  40. 40.

    Chua PK, McCown MF, Rajyaguru S, Kular S, Varma R, Symons J, Chiu SS, Cammack N, Najera I. Modulation of alpha interferon anti-hepatitis C virus activity by ISG15. J Gen Virol. 2009;90(12):2929–39.

  41. 41.

    Xia Y, Paul BY, Sidis Y, Beppu H, Bloch KD, Schneyer AL, Lin HY. Repulsive guidance molecule RGMa alters utilization of bone morphogenetic protein (BMP) type II receptors by BMP2 and BMP4. J Biol Chem. 2007;282(25):18129–40.

  42. 42.

    Lin W, Choe WH, Hiasa Y, Kamegaya Y, Blackard JT, Schmidt EV, Chung RT. Hepatitis C virus expression suppresses interferon signaling by degrading STAT1. Gastroenterol. 2005;128(4):1034–41.

  43. 43.

    Ogishima T, Shiina H, Breault JE, Terashima M, Honda S, Enokida H, Urakami S, Tokizane T, Kawakami T, Ribeiro-Filho LA, Fujime M. Promoter CpG hypomethylation and transcription factor EGR1 hyperactivate heparanase expression in bladder cancer. Oncogene. 2005;24(45):6765.

  44. 44.

    Oridate N, Kim HJ, Xu X, Lotan R. Growth inhibition of head and neck squamous carcinoma cells by small interfering RNAs targeting eIF4E or cyclin D1 alone or combined with cisplatin. Cancer Biol Ther. 2005;4(3):318–23.

  45. 45.

    Walsh N, O’Donovan N, Kennedy S, Henry M, Meleady P, Clynes M, Dowling P. Identification of pancreatic cancer invasion-related proteins by proteomic analysis. Proteome Sci. 2009;7(1):3.

  46. 46.

    Yang L, et al. siRNA-mediated silencing of Wnt5a regulates inflammatory responses in atherosclerosis through the MAPK/NF-kappaB pathways. Int J Mol Med. 2014;34(4):1147–52.

  47. 47.

    Lefevre G, et al. Roles of stem cell factor/c-Kit and effects of Glivec/STI571 in human uveal melanoma cell tumorigenesis. J Biol Chem. 2004;279(30):31769-79.

  48. 48.

    Morrow AN, Schmeisser H, Tsuno T, Zoon KC. A novel role for IFN-stimulated gene factor 3II in IFN-γ signaling and induction of antiviral activity in human cells. J Immunol. 2010:1001359.

  49. 49.

    Toyokawa K, Leite F, Ott TL. Cellular localization and function of the antiviral protein, ovine Mx1 (oMx1): II. The oMx1 protein is a regulator of secretion in an ovine glandular epithelial cell line. Am J Reprod Immunol. 2007;57(1):23–33.

  50. 50.

    Chen S, Charness ME. Ethanol inhibits neuronal differentiation by disrupting activity-dependent neuroprotective protein signaling. Proc Natl Acad Sci. 2008; https://doi.org/10.1073/pnas.0807758105.

  51. 51.

    Arimoto KI, Hishiki T, Kiyonari H, Abe T, Cheng C, Yan M, Fan JB, Futakuchi M, Tsuda H, Murakami Y, Suzuki H. Murine Herc6 plays a critical role in protein ISGylation in vivo and has an ISGylation-independent function in seminal vesicles. J Interf Cytokine Res. 2015;35(5):351–8.

  52. 52.

    Bonnevie-Nielsen V, et al. Variation in antiviral 2',5'- oligoadenylate synthetase (2'5'AS) enzyme activity is controlled by a single-nucleotide polymorphism at a splice-acceptor site in the OAS1 gene. Am J Hum Genet. 2005;76(4):623–33.

  53. 53.

    Bridge AJ, Pebernard S, Ducraux A, Nicoulaz AL, Iggo R. Induction of an interferon response by RNAi vectors in mammalian cells. Nat Genet. 2003;34(3):263.

  54. 54.

    Barkhash AV, Perelygin AA, Babenko VN, Myasnikova NG, Pilipenko PI, Romaschenko AG, Voevoda AG, Brinton MA. Variability in the 2nko A oligoadenylate synthetase gene cluster is associated with human predisposition to tick-borne encephalitis virus-induced disease. J Infect Dis. 2010;202(12):1813–8.

  55. 55.

    Chawla-Sarkar M, Lindner DJ, Liu YF, Williams BR, Sen GC, Silverman RH, Borden EC. Apoptosis and interferons: role of interferon-stimulated genes as mediators of apoptosis. Apoptosis. 2003;8(3):237–49.

  56. 56.

    Wang SS, Gonzalez P, Yu K, Porras C, Li Q, Safaeian M, Rodriguez AC, Sherman ME, Bratti C, Schiffman M, Wacholder S. Common genetic variants and risk for HPV persistence and progression to cervical cancer. PLoS One. 2010;5(1):e8667.

  57. 57.

    Boulbes DR, Chauhan GB, Jin Q, Bartholomeusz C, Esteva FJ. CD44 expression contributes to trastuzumab resistance in HER2-positive breast cancer cells. Breast Cancer Res Treat. 2015;151(3):501–13.

  58. 58.

    Pályi-Krekk Z, Barok M, Isola J, Tammi M, Szöllo J, Nagy P. Hyaluronan-induced masking of ErbB2 and CD44-enhanced trastuzumab internalisation in trastuzumab resistant breast cancer. Eur J Cancer. 2007;43(16):2423–33.

  59. 59.

    Bao W, Fu HJ, Xie QS, Wang L, Zhang R, Guo ZY, Zhao J, Meng YL, Ren XL, Wang T, Li Q. HER2 interacts with CD44 to up-regulate CXCR4 via epigenetic silencing of microRNA-139 in gastric cancer cells. Gastroenterol. 2011;141(6):2076–87.

  60. 60.

    Barthelemy M. Betweenness centrality in large complex networks. Eur Phys J B. 2004;38(2):163–8.

  61. 61.

    Lastraioli E, Iorio J, Arcangeli A. Ion channel expression as promising cancer biomarker. BBA Biomembranes. 2015;1848(10):2685–702.

  62. 62.

    Xia J, Wang H, Li S, Wu Q, Sun L, Huang H, Zeng M. Ion channels or aquaporins as novel molecular targets in gastric cancer. Mol Cancer. 2017;16(1):54.

  63. 63.

    Fujimoto M, Inoue T, Kito H, Niwa S, Suzuki T, Muraki K, Ohya S. Transcriptional repression of HER2 by ANO1 Cl− channel inhibition in human breast cancer cells with resistance to trastuzumab. Biochem Biophys Res Commun. 2017;482(1):188–94.

  64. 64.

    Cai R, Ding X, Zhou K, Shi Y, Ge R, Ren G, Jin Y, Wang Y. Blockade of TRPC6 channels induced G2/M phase arrest and suppressed growth in human gastric cancer cells. Int J Cancer. 2009;125(10):2281–7.

  65. 65.

    Ding M, et al. Pyrazolo [1,5-a] pyrimidine TRPC6 antagonists for the treatment of gastric cancer. Cancer Lett. 2018;432:47-55.

  66. 66.

    Huang R, Gu W, Sun B, Gao L. Identification of COL4A1 as a potential gene conferring trastuzumab resistance in gastric cancer based on bioinformatics analysis. Mol Med Rep. 2018;17(5):6387–96.

  67. 67.

    Arteaga CL, Sliwkowski MX, Osborne CK, Perez EA, Puglisi F, Gianni L. Treatment of HER2-positive breast cancer: current status and future perspectives. Nat Rev Clin Oncol. 2012;9(1):16.

Download references

Acknowledgements

We would like to thank the Shanghai Institute of Digestive Surgery, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine for academic support.

Funding

The study is financially supported by the National Natural Science Foundation of China (NSFC) (81402423, 81572818) and Shanghai Municipal Commission of Health and Family Planning (2017YQ062).

Availability of data and materials

The datasets supporting the conclusion of this article are included within the article.

Human participants and animal rights

This article does not contain any studies with human participants or animals performed by any of the authors.

Author information

CY, LZ, PX, and RP carried out data analysis. CY, ZC, ZH, JS, and MZ drafted the manuscript. JS, MZ, and CY participated in the study design and data collection. All authors read and approved the final manuscript.

Correspondence to Jing Sun or Minhua Zheng.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. 20 hub genes with siRNA synthesizers and sequence [34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52]. (DOCX 19 kb)

Additional file 2:

Table S2. The adj. p value and log fold change (logFC) of the hub genes in TCGA and GSE77346 datasets. (DOCX 19 kb)

Additional file 3:

Figure S1. The mRNA expression of hub genes in GSE13861. (PDF 31 kb)

Additional file 4:

Figure S2. The mRNA expression of hub genes in clinical stages in TCGA. (PDF 2878 kb)

Additional file 5:

Table S3. Targeted drugs with corresponding GEO datasets for gastric cancer. (DOCX 18 kb)

Additional file 6:

Figure S3. The mRNA expression and the fold change of the hub genes. (A) the mRNA expression of hub genes in GSE19043 (Gefitinib); (B) the fold change of hub genes in GSE95414. (PDF 740 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Differentially expressed genes
  • Gene ontology
  • KEGG pathway
  • Gastric cancer
  • Trastuzumab
  • Resistance
  • Protein-protein interaction