Chenyu Sun1†
Yue Chen2†
Na Hyun Kim1
Scott Lowe3Shaodi Ma4
Zhen Zhou5Rachel Bentley3
Yi-Sheng Chen6Margarita Whitaker Tuason7
Wenchao Gu8Chandur Bhan1John Pocholo Whitaker Tuason1Pratikshya Thapa1
Ce Cheng9,10
Qin Zhou11
Yanzhe Zhu12*- 1AMITA Health Saint Joseph Hospital Chicago, Chicago, IL, United States
- 2Department of Clinical Medicine, School of the First Clinical Medicine, Anhui Medical University, Hefei, China
- 3College of Osteopathic Medicine, Kansas City University, Kansas City, MO, United States
- 4Department of Epidemiology and Health Statistics, School of Public Health, Anhui Medical University, Hefei, China
- 5Menzies Institute for Medical Research, University of Tasmania, Hobart, TAS, Australia
- 6Department of Orthopedics, Shanghai General Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai Jiao Tong University, Shanghai, China
- 7Faculty of Medicine and Surgery, University of Santo Thomas, Metro Manila, Philippines
- 8Department of Diagnostic Radiology and Nuclear Medicine, Gunma University Graduate School of Medicine, Maebashi, Japan
- 9The University of Arizona College of Medicine, Tucson, AZ, United States
- 10Banner-University Medical Center South, Tucson, AZ, United States
- 11Mayo Clinic, Rochester, MN, United States
- 12Department of Oncology, The First Affiliated Hospital of Anhui Medical University, Hefei, China
Background: Gastric cancer (GC) is a common cancer with high mortality. This study aimed to identify its differentially expressed genes (DEGs) using bioinformatics methods.
Methods: DEGs were screened from four GEO (Gene Expression Omnibus) gene expression profiles. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed. A protein–protein interaction (PPI) network was constructed. Expression and prognosis were assessed. Meta-analysis was conducted to further validate prognosis. The receiver operating characteristic curve (ROC) was analyzed to identify diagnostic markers, and a nomogram was developed. Exploration of drugs and immune cell infiltration analysis were conducted.
Results: Nine up-regulated and three down-regulated hub genes were identified, with close relations to gastric functions, extracellular activities, and structures. Overexpressed Collagen Type VIII Alpha 1 Chain (COL8A1), Collagen Type X Alpha 1 Chain (COL10A1), Collagen Triple Helix Repeat Containing 1 (CTHRC1), and Fibroblast Activation Protein (FAP) correlated with poor prognosis. The area under the curve (AUC) of ADAM Metallopeptidase With Thrombospondin Type 1 Motif 2 (ADAMTS2), COL10A1, Collagen Type XI Alpha 1 Chain (COL11A1), and CTHRC1 was >0.9. A nomogram model based on CTHRC1 was developed. Infiltration of macrophages, neutrophils, and dendritic cells positively correlated with COL8A1, COL10A1, CTHRC1, and FAP. Meta-analysis confirmed poor prognosis of overexpressed CTHRC1.
Conclusion: ADAMTS2, COL10A1, COL11A1, and CTHRC1 have diagnostic values in GC. COL8A1, COL10A1, CTHRC1, and FAP correlated with worse prognosis, showing prognostic and therapeutic values. The immune cell infiltration needs further investigations.
Background
As one of the top five common malignancies, gastric cancer (GC) is the fourth leading cause of cancer-associated death around the globe, according to Global cancer statistics 2020. It accounted for more than one million new cases in 2020 and is responsible for one in every 13 deaths globally (Sung et al., 2021). The incidence and mortality of GC vary among different populations and geographical locations. GC is the leading cause of cancer death in several South Central Asian countries. Incidence rates are highest in Eastern Asia and Eastern Europe, while the incidence rates are generally low in Northern America, Northern Europe, and African regions (Sung et al., 2021). These differences could be attributed to the various environmental risk factors, such as the different prevalence of Helicobacter pylori (H. pylori) infections, alcohol consumptions, tobacco smoking, consumption of preserved salty food and processed meat, ingestion of grilled or barbecued meat and fish, and viral infection (Hooi et al., 2017; Research. WCRFAIfC, 2018; Palrasu et al., 2021; Sung et al., 2021). In addition to these environmental factors, genetic factors were also thought to affect the carcinogenesis of GC, as less than 5% of H. pylori infected hosts will develop GC, and evidence of genetic alterations, such as aberrantly expressed activation-induced cytidine deaminase (AID), has emerged (Nagata et al., 2014; Nakanishi et al., 2021). Significant progress in the diagnosis and treatment of GC has been made, such as development of novel human epidermal growth factor receptor 2 (HER2)-targeted drugs for GC (Zhu et al., 2021), the development of minimally invasive surgery, and endoscopic mucosal resection (EMR) and endoscopic submucosal dissection (ESD) (Roh et al., 2020; Wang et al., 2021). Moreover, several biomarkers have been found for the target therapy of GC, such as programmed death 1 (PD-1), HER2, and MNNG HOS transforming gene (MET) (Choi et al., 2022). PD-1 is an inhibitory checkpoint receptor protein expressed on cytotoxic T cells and other immune cells (Pardoll, 2012). Some tumor cells express high levels of PD-L1 to evade from immune system, as PD-1/PD-L1 interaction induces cytotoxic T cell inactivation and downregulation of immune responses (Schreiber et al., 2011), and PD-L1 expression was proposed to be a potential biomarker of response to pembrolizumab (Ghidini et al., 2021). HER2 overexpression is particularly important in GC, as targeted therapy trastuzumab has been widely used to treat HER2+ GC (Bang et al., 2010). MET activation triggers a downstream cascade of phosphoinositide 3-kinases (PI3K) and Rat sarcoma virus (RAS) signaling and regulates cell survival and proliferation (Zhang et al., 2018), Thus, the over-activation of MET plays a critical role in cancer development and is frequently identified in various types of tumors, including GC (Ariyawutyakorn et al., 2016). However, the prognosis of GC remains quite unsatisfactory due to its low early diagnosis rate, with a 5-year overall survival (OS) of less than 40% (Tan, 2019). Therefore, exploration of novel biomarkers that are sensitive and specific for early diagnosis, as well as predictors of prognosis and response to potential targeted treatment, is pivotal in the management of GC.
With the development of next-generation sequencing (NGS) and other techniques, the availability of information related to these potential biomarkers and knowledge of the relevant gene expressions available for various tumors have increased significantly (Behjati and Tarpey, 2013; Levy and Myers, 2016; Giunchi et al., 2021; Shirdarreh et al., 2021). As a result, the mechanisms of various cancers and other diseases have become more widely studied based on bioinformatic analysis, a field combining molecular biology and information technology. Bioinformatics methods, such as data-mining, are now commonly used to explore the carcinogenesis at the molecular level, and to explore biomarkers for potential diagnostic markers, prognostic predictors, and therapeutic targets (Li et al., 2018a; Shen et al., 2019; Yang et al., 2019). As several microarray profiling studies have been performed in GC, this study integrating publicly available data of some of these existing studies, to search for the differentially expressed genes (DEGs) and ultimately biomarkers that could show potential diagnostic values, predict prognosis, and those that might become therapeutic targets.
Methods
Data Collection
Four gene expression profiles [GSE13911 (D'Errico et al., 2009; Wang et al., 2012a), GSE1982626, GSE54129, and GSE79973 (He et al., 2016; Jin et al., 2017)] were downloaded from the GEO database. The patient’s data from GEO datasets were obtained if the pathological sample/biopsy of the gastric cancer in the experiment group or normal gastric tissues in control group were used. GSE13911 dataset included 38 cancer tissues and 31 non-cancerous tissues, whilst GSE19826 dataset included 12 cancer tissues and 12 non-cancer tissues, GSE54129 dataset included 111 cancer tissues and 21 non-cancer tissues, and GSE79973 dataset included 10 cancer tissues and 10 non-cancer tissues. All datasets were based on GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array.
Identification of DEGs
The linear models for microarray data (LIMMA) package (Ritchie et al., 2015) based on R software was utilized to screen up-regulated DEGs within the adj. p < 0.01 and Log2FC > 2, and down-regulated DEGs within the adj. p < 0.01 and Log2FC < -2 between samples in cancer group and non-cancer group. To identify overlapping DEGs, a Venn diagram was constructed using the bioinformatics & evolutionary genomics website (http://bioinformatics.psb.ugent.be/webtools/Venn/). Volcano plot of the four datasets was drawn by Hiplot (https://hiplot.org).
Protein–Protein Interaction Network Construction
PPI was conducted by using the Search Tool For the Retrieval of Interacting Genes (STRING) database (Szklarczyk et al., 2019), and Cytoscape software (Shannon et al., 2003). The network was constructed based on setting the medium confidence as >0.4 in the STRING database, and then imported the network into Cytoscape software for further analysis. Hub genes were selected based on plugin Cytohubba to identify hub genes through 12 algorithms (Chin et al., 2014; Zhang et al., 2019a). The top 15 genes ranked by score from each algorithm were extracted and mutual genes that were overlapped in all 12 algorithms were selected as hub genes. PPI network of the hub genes was visualized by GeneMANIA (http://www.genemania.org) (Vlasblom et al., 2015).
Validation of Hub Genes: Survival Analysis, Expression Analysis, and Receiver Operating Characteristic Curve Analysis
Survival analysis and direct tumor/normal differential expression were conducted for the selected hub genes via data obtained from The Cancer Genome Atlas (TCGA) dataset (https://portal.gdc.cancer.gov/). Log2 transformed FPKM (fragments per kilobase exon-model per million reads mapped) were used. Images of immunohistochemical (IHC) staining for the protein expressed by up-regulated genes were obtained from the Human Protein Atlas (HPA) (http://www.proteinatlas.org/) to evaluate their expressions in GC (Pontén et al., 2011; Zhu et al., 2022). The Kaplan–Meier (KM) survival analysis with log-rank test was also used to compare the OS difference between the high expression and low expression group. KM curves, with p-values and hazard ratio (HR) with 95% confidence interval (CI), generated by log-rank tests and univariate Cox proportional hazards regression were performed using R software version v3.6.3 (The R Foundation for Statistical Computing, 2020) with “survminer,” and “survival” packages. Expression analysis was performed by using Wilcoxon rank sum test, and visualized by “ggplot2” package of R software. p < 0.05 was considered as statistically significant. Then, genes with significant worse overall survival (OS) were also verified by Kaplan–Meier Plotter (Szász et al., 2016).
The best discriminate cut-off point of overexpressed DEGs between the high and low expression groups were assessed by the receiver operating characteristic (ROC) curve and area under the curve (AUC) values, based on data obtained from TCGA. Log2 transformed FPKM were used. R software with “pROC” and “ggplot2” packages were used.
Nomogram Development
A predictive model was established to predict the mortality risk based on the overexpressed hub genes with worst outcomes and other potential predictors (Iasonos et al., 2008; Balachandran et al., 2015; Liu et al., 2018a). A nomogram was developed based on the results of multivariate Cox proportional hazards analysis through “rms” and “survival” R packages. Data were obtained from TCGA. The nomogram provided a graphical representation of the factors to calculate the risk of mortality at 1, 3,and 5-year time points for an individual patient by the points associated with each risk factor. C-index was also calculated to assess the discriminatory performance of the model (Liu et al., 2018a; Kramer and Zimmerman, 2007; Pencina and D'Agostino, 2004).
Exploration of Potential Drugs That Are Interacted With Hub Genes That Were Associated With Poor Prognosis in LC
We explored potential drugs that are interacted with hub genes linked to poor prognosis explored by using RNAactDrug (http://bio-bigdata.hrbmu.edu.cn/RNAactDrug/index.jsp), which is a comprehensive database for exploring associations between drug sensitivity and RNA molecules at expression level and other molecular levels from integrated analysis of three large-scale pharmacogenomic databases (GDSC, CellMiner and CCLE) (Dong et al., 2020).
Go Enrichment and KEGG Pathway Analysis
The overlapping up-regulated and down-regulated DEGs were analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) in the Database for Annotation, Visualization, and Integrated Discovery (DAVID) database (https://david.ncifcrf.gov/summary.jsp) (Huang et al., 2009a; Huang et al., 2009b). GO enrichment analysis predicted the function based on biological processes (BP), cellular components (CC), and molecular functions (MF), while KEGG analysis determined the related pathways of hub genes and their associated interactors. The results of GO and KEGG analyses were visualized by the bioinformatics online tool (http://www.bioinformatics.com.cn) (Wang et al., 2020a; Zhang et al., 2021).
Immune Cell Infiltration of the Hub Genes With Worse Prognosis in GC
The infiltration of different immune cells and their clinical impact were assessed via Tumor Immune Estimation Resource (TIMER) (https://cistrome.shinyapps.io/timer/), an online tool for comprehensive investigation of molecular characterization of tumor-immune interactions based on 10,897 tumors from 32 cancer types (Li et al., 2016; Li et al., 2017). Hub genes that were associated with poor prognosis were entered into the “Gene module” to generate plots for analyzing the correlation between their expressions and immune infiltration level in GC. Positive correlation was considered for the cuff value of Cor >0.2 and p < 0.05 (Liu et al., 2021; Zhong et al., 2021).
Meta-Analysis to Verify the Results of Survival Analysis of the Hub Genes
Meta-analysis was conducted to verify the results of the survival analysis of the hub genes associated with poor prognosis. Electronic databases including China National Knowledge Infrastructure database (CNKI), Web of science, and PubMed were searched to find eligible articles for conducting a meta-analysis to explore and verify survival analyses of the hub genes that were associated with worse prognosis in gastric cancer, namely, Collagen Triple Helix Repeat Containing 1 (CTHRC1). The search strategy included the following: (CTHRC1 OR collagen triple helix repeat containing 1 OR COL8A1 OR collagen type VIII alpha 1 chain OR COL10A1 OR collagen type X alpha 1 chain OR FAP OR fibroblast activation protein) AND (gastric cancer OR stomach adenocarcinoma OR gastric adenocarcinoma OR stomach cancer OR STAD). For CNKI database, the corresponding Chinese expression was used. The meta-analysis was conducted according to the Preferred Reporting Items declared by the Systematic Review and Meta-Analysis (PRISMA) (Page et al., 2021). Finally, the HR estimate with 95% CI was calculated using the effect values extracted from the incorporated articles and the corresponding result of our survival analysis. Q test and I2 statistics were used to evaluate the extent of heterogeneity across the studies (DerSimonian and Laird, 1986). If significant heterogeneity (I2 statistic >50% or Q test <0.1) was observed, then a random-effects model was used, otherwise a fixed-effects model was applied for combined HRs (Borenstein et al., 2010). Sensitivity analysis was performed by switching between fixed and random effects models, for testing the stability of the study results (Hernandez et al., 2020; Dan Song and Zhang, 2021). All statistical analyses were performed using STATA software (version 15.0).
Results
Identification of DEGs
As shown in Figure 1A, there were 63 up-regulated DEGs from GSE79973, 31 DEGS from GSE19826, 41 DEGs from GSE19311, and 303 DEGs from GES54129. Among them, there were 23 up-regulated overlapping DEGs. As shown in Figure 1B, there were 155 down-regulated DEGs from GSE79973, 74 DEGS from GSE19826, 62 DEGs from GSE19311, and 225 DEGs from GES54129. Among them, 77 down-regulated overlapping DEGs were found. Volcano plot of gene expression profile in non-cancer compared to cancer groups of the four GSE datasets is shown in Figure 1C.
FIGURE 1. (A) Venn diagram of overlapped up-regulated differentially expressed genes (DEGs); (B) Venn diagram of overlapped down-regulated DEGs; (C) Volcano plot of gene expression profiles in non-cancer compared to cancer groups.
PPI Network Construction and Identification of Hub Genes
Based on the STRING database and Cytoscape software, PPI networks of up-regulated overlapping DEGs (22 nodes and 46 edges) and down-regulated DEGs (71 nodes and 65 edges) were constructed. Up-regulated hub genes included Collagen Type I Alpha 2 Chain (COL1A2), thrombospondin 2 (THBS2), Collagen Type XI Alpha 1 Chain (COL11A1), Collagen Type VIII Alpha 1 Chain (COL8A1), Collagen Type X Alpha 1 Chain (COL10A1), ADAM Metallopeptidase With Thrombospondin Type 1 Motif 2 (ADAMTS2), CTHRC1, fibroblast activation protein (FAP), and WNT1-inducible-signaling pathway protein 1 (WISP1), and down-regulated hub genes included Trefoil Factor 2 (TFF2), gastric intrinsic factor (GIF), and Cytochrome P450 family 2 subfamily C member 9 (CYP2C9). PPI network of these hub genes were visualized via GeneMANIA as shown in Figure 2.
FIGURE 2. PPI network of different expressed hub genes and their interactors visualized by GeneMANIA.
Validation of Hub Gene Expression Analysis and Survival Analysis
The expression levels of the up-regulated hub genes between cancer tissues and non-cancer tissues were all statically significant. Expressions of ADAMTS2, COL1A2, COL8A1, COL10A1, COL11A1, CTHRC1, FAP, THBS2, and WISP1 were higher in cancer samples. However, among down-regulated hub genes, there was no statistically significant difference between the expression level of CYP2C9 in cancer samples and normal gastric samples, while the expressions of TFF2 and GIF were statistically significantly lower in cancer samples than normal gastric tissues (Figure 3A). Survival analysis found that higher expressed COL8A1, COL10A1, CTHRC1, and FAP were associated with poor OS. (Figure 3B; Table 1) In addition, the OS of COL8A1, COL10A1, CTHRC1, and FAP verified by Kaplan–Meier plotter firmed their association with worse OS except for FAP. (Figure 4) Immune staining of ADAMTS2, COL1A2, COL8A1, CTHRC1, FAP, THBS2, and WISP1 proteins were shown in Figure 5, however, COL10A1 and COL11A1 were not available from HPA.
FIGURE 3. (A) Expression levels of the hub genes between cancer samples and normal gastric tissues (B) Survival analysis of the hub genes.
TABLE 1. Overall survival analysis of Collagen Type VIII Alpha 1 Chain (COL8A1), Collagen Type X Alpha 1 Chain (COL10A1), Collagen Triple Helix Repeat Containing 1 (CTHRC1), and Fibroblast Activation Protein (FAP).
FIGURE 4. Survival analysis of Collagen Type VIII Alpha 1 Chain (COL8A1), Collagen Type X Alpha 1 Chain (COL10A1), Collagen Triple Helix Repeat Containing 1 (CTHRC1), and Fibroblast Activation Protein (FAP) by Kaplan–Meier Plotter.
FIGURE 5. Representative images of immunohistochemical staining of normal gastric tissue (panel a) and gastric cancer tissues (panel b) based on the Human Protein Atlas: (A) ADAM Metallopeptidase With Thrombospondin Type 1 Motif 2 (ADAMTS2); (B) Collagen Type I Alpha 2 Chain (COL1A2); (C) Collagen Type VIII Alpha 1 Chain (COL8A1); (D) CTHRC1; (E) FAP; (F) Thrombospondin 2 (THBS2); (G) WNT1-inducible-signaling pathway protein 1 (WISP1).
ROC Analysis
Further analysis of ROC curve demonstrated that AUC values of ADAMTS2, COL10A1, COL11A1, and CTHRC1 were 0.937 (95% CI: 0.910–0.963), 0.973 (95% CI: 0.959–0.988), 0.934 (95% CI: 0.906–0.962), and 0.966 (95% CI: 0.949–0.983), respectively, (Figure 6). Cutoff values of ADAMTS2, COL10A1, COL11A1, and CTHRC1 were 1.512, 0.382, 0.194, and 2.410.
FIGURE 6. Receiver operating characteristics curve (ROC) curve analysis for the hub genes in gastric cancer.
Nomogram Development
A nomogram model incorporating CTHRC1, the hub gene with the worst OS, and other predictors (age, gender, reflux history, Barrett’s esophagus, H. Pylori infection, pathologic stage, histologic grade, resident tumor) is shown in Figure 7. The C-index of the nomogram was 0.709 (95% CI, 0.678–0.740).
FIGURE 7. Developed overexpressed CTHRC1 nomogram. Note: The nomogram was developed in the cohort, with age, gender, reflux history, Barret’s esophagus, Helicobacter pylori (H. Pylori) infection, pathologic stage, histologic grade, resident tumor. (C-index: 0.709, 95% CI, 0.678–0.740).
Drug Interactions With Hub Genes That Were Associated With Poor Prognosis in LC
As overexpressed COL8A1, COL10A1, CTHRC1, and FAP were found to be associated with worse OS, the top ten drugs or compounds that demonstrated the strongest association with these genes were identified, based on the significance of Spearman correlation. (Table 2)
TABLE 2. Top ten drugs associated with COL8A1, COL10A1, and CTHRC1 at expression level.
Go Enrichment and KEGG Pathway Analysis
Go enrichment and KEGG pathway analysis of the up-regulated and down-regulated overlapping DEGs was conducted by using DAVID and visualized by using bioinformatics online tool (www.bioinformatics.com.cn). Digestion, collagen catabolic process, and extracellular matrix organization were the top three biological processes that were associated with the up-regulated and down-regulated overlapping hub genes (Figure 8A). The extracellular space, extracellular region, and proteinaceous extracellular matrix were the top three major cellular components of these hub genes (Figure 8B). As for molecular function, extracellular matrix structural constituent, inward rectifier potassium channel activity, and extracellular matrix binding were the top three functions (Figure 8C). In regard to KEGG pathways, gastric acid secretion, chemical carcinogenesis, and extracellular matrix (ECM)-receptor interaction were the top three pathways involved (Figure 8D).
FIGURE 8. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis (DAVID). GO enrichment analysis of target genes based on (A) biological process, (B) cellular component, and (C) molecular function. (D) KEGG pathway enrichment analysis of target genes.
Immune Cell Infiltration of the Hub Genes With Worse Prognosis in GC
The TIMER database was utilized to investigate the association between COL8A1, COL10A1, CTHRC1, and FAP, and immune cell infiltration, as immune cell levels correlate with the proliferation and progression of cancer cells (Figure 9). The infiltrations of macrophages, neutrophils, and dendritic cells positively correlated with COL8A1, COL10A1, CTHRC1, and FAP. In addition, the Cox proportional hazard model showed that macrophage (p = 0.001) and CTHRC1 (p = 0.021) were significantly associated with adverse clinical outcomes in GC patients (Table 3).
FIGURE 9. Correlations between COL8A1, COL10A1, CTHRC1, and FAP and immune cell infiltration in gastric cancer (TIMER).
TABLE 3. The Cox proportional hazard model of COL8A1, COL10A1, CTHRC1, and FAP and six tumor-infiltrating immune cells in gastric cancer (TIMER).
Meta-Analysis
To verify the results of survival analysis of the hub genes that were associated with worse OS, a meta-analysis was performed. Despite a comprehensive literature search was performed, only three articles investigating CTHRC1 were initially selected for full-text review, while no eligible articles on other hub genes were identified for full-text review. CTHRC1 was associated with the highest HR compared with other hub genes based on our survival analysis, and was shown to have increased HR based on immune cell infiltration analysis; therefore, we pooled the result of the survival analysis of our bioinformatic study with three original articles retrieved from databases (Gu et al., 2014; Wang, 2016a; Dan Song and Zhang, 2021). However, Wang’s article (Wang, 2016b) was excluded due to the incorrect data. The pooled result showed that positive expression of CTHRC1 was associated with poor prognosis of gastric cancer patients. (HR: 1.93, 95% CI: 1.32–2.82, I2 = 66.1%). (Figure 10) Sensitivity analysis by changing random-effect model to fixed-effect model did not change the result significantly.
FIGURE 10. Forrest plot of meta-analysis for CTHRC1.
Discussion
As one of the most common malignant tumors and top leading causes of cancer-related death, the GC is induced by a variety of factors (Nagata et al., 2014; Hooi et al., 2017; Research. WCRFAIfC, 2018; Nakanishi et al., 2021; Palrasu et al., 2021; Sung et al., 2021). In this study, 9 up-regulated hub genes (COL1A2, THBS2, COL11A1, COL8A1, COL10A1, ADAMTS2, CTHRC1, FAP, and WISP1) and 3 down-regulated hub genes (TFF2, GIF, and CYP2C9) were identified through integrating 4 gene expression profiles to screen for DEGs and perform PPI network analysis. Previous studies have already identified a number of potential hub genes in GC by using similar research methods. For example, a bioinformatics analysis by Liu et al. (Liu et al., 2018b) found nine hub genes, including DNA Topoisomerase II Alpha (TOP2A), collagen type I alpha 1 chain (COL1A1), COL1A2, NDC80 Kinetochore Complex Component (NDC80), Collagen Type III Alpha 1 Chain (COL3A1), Cyclin Dependent Kinase Inhibitor 3 (CDKN3), Centrosomal Protein 55 (CEP55), Targeting protein for Xklp2 (TPX2), and TIMP Metallopeptidase Inhibitor 1 (TIMP1), that might be associated with the pathogenesis of GC. In addition, the results of the article by Wang et al. (Wang et al., 2020b) identified that Fibronectin 1 (FN1), COL1A1, Inhibin beta-A (INHBA), and cystatin SN (CST1) might be potential biomarkers and therapeutic targets for GC patients. Moreover, Lu et al. (Lu et al., 2021) found 5 key genes, namely, Hyaluronan-Mediated Motility Receptor (HMMR), Cyclin B1(CCNB1), C-X-C motif chemokine ligand 8 (CXCL8), mitotic arrest deficient 2 like 1 (MAD2L1), and Cyclin A2 (CCNA2), in GC patients with poor prognosis using the datasets from GEO database. Additionally, Liu et al. suggested that the expression levels of (ATPase H+/K + Transporting Subunit Alpha (ATP4A), carbonic anhydrase 9 (CA9), Fibrinogen Alpha Chain (FGA), Aldehyde Dehydrogenase 1 Family Member A1 (ALDH1A1), and Ghrelin And Obestatin Prepropeptide (GHRL) were reduced, whereas those of TIMP1, Secreted Phosphoprotein 1 (SPP1), CXCL8, Thy-1 Cell Surface Antigen (THY1), and COL1A1 were increased in GC. Another study demonstrated that COL1A1, Collagen Type V Alpha 2 Chain (COL5A2), Prolyl 4-Hydroxylase Subunit Alpha 3 (P4HA3), and Secreted Protein Acidic And Cysteine Rich (SPARC) showed vital values in prognosis and diagnosis of GC (Niu et al., 2022). Furthermore, another study conducted by Dalkilic to analyze GC transcriptomic data revealed that Secreted Frizzled Related Protein 2 (SFRP2), Early Growth Response 1 (EGR1), Chitinase 3 Like 1 (CHI3L1), COL8A1, Nuclear Enriched Abundant Transcript 1 (NEAT1), INHBA, CXCL8, and Myosin Light Chain 9 (MYL9) were highly expressed, while expression of Gastrin (GAST), GIF, Gastrokine 1 (GKN1), Gastrokine 2 (GKN2), Secretoglobin Family 2A Member 1 (SCGB2A1), and HRAS-like suppressor 2 (HRASLS2) were downregulated (Dalkilic, 2020). Although a variation of hub genes were found in previous studies, our study added additional evidence of potential hub genes in GC that could be served as important biomarkers.
As expected, GO analysis found that all the hub genes were related to the digestion pathway for the biological process. Collagen catabolic process and extracellular matrix were also involved significantly in the biological process. For cellular components, the extracellular space, and extracellular regions, as well as proteinaceous extracellular matrix, were mostly involved. In regards to molecular functions, extracellular matrix structural constituent, inward rectifier potassium channel activity, and extracellular matrix binding were greatly involved. KEGG analysis identified that the gastric acid secretion pathway was significantly involved, which was unsurprising given that interference with gastric acid secretion may have caused damage to gastric mucosa. In addition, chemical carcinogenesis and ECM-receptor interaction were also found to be involved. Most of the pathways found in GO and KEGG analyses were closely related to either gastric functions or extracellular activities and structures. A summary of these biological functions and pathways that the hub genes were mostly involved in are shown in Table 4.
TABLE 4. Biological functions and pathways that the hub genes are mostly involved in.
Further validation of these hub genes found that all up-regulated hub genes were expressed higher in GC tissues, while among the down-regulated hub genes, the expression level of CYP2C9 was not significantly lower in GC samples. Further survival analysis of the up-regulated hub genes found that COL8A1, COL10A1, CTHRC1, and FAP were associated with poor OS. COL8A1 was known to encode the alpha 1 chain of collagen, type VIII, and thus may modulate migration, proliferation, and adherence of various cells (Zhao et al., 2009a). COL8A1 was also proposed to promote the migration of certain cancer cells by mediating ECM-receptor interaction (Peng et al., 2020). Zhao et al. showed that knockdown of COL8A1 induced the inhibition of hepatocellular carcinoma growth and invasion (Zhao et al., 2009a). Previous studies also found that COL8A1 was associated with poor prognosis of GC, consistent with our findings (Wu et al., 2020a; Chen et al., 2020; Wang et al., 2020c; Zhou et al., 2020). COL10A1, which encodes Collagen type X alpha 1, was found to be overexpressed in different types of cancer, such as esophageal cancer and breast cancer, and promoted the malignant progression by upregulating the expression of Prolyl 4-hydroxylase beta polypeptide (P4HB) (Song et al., 2021; Yang et al., 2021). As a member of the collagen family, COL10A1could activate ECM remodeling and the epithelial-to-mesenchymal transition (EMT), vascular endothelial growth factor receptor 3 (VEGFR3) and Wnt signaling pathway, and its aberrant expression might affect the development of cancer (Song et al., 2021). A recent study identified COL10A1 as a potential inducer of EMT (Li et al., 2018b). Silencing of COL10A1 was found to induce inhibition of cell proliferation, migration, and invasion in GC (Li et al., 2020). Consistent with previous studies, we also found that COL10A1 was associated with adverse outcomes in GC (Chen et al., 2021). Collagen triple helix repeat containing-1 (CTHRC1) was known as a cancer-related protein, and overexpression of CTHRC1 was believed to be involved in tumorigenesis, proliferation, invasion, and metastasis in various gastrointestinal malignancies, including gastric cancer, as well as other non-gastrointestinal malignant tumors (Wang et al., 2012b; Mei et al., 2020). A recent in vitro study found that CTHRC1 increased CXCR4 expression through upregulating hypoxia-inducible factor 1-alpha (HIF-1α) expression, leading to the promotion of cell migration and invasion in GC (Ding et al., 2020). Another study found that repression of CTHRC1 protein activity could inhibit cell proliferation, migration, and invasion in GC (Yu et al., 2015). In addition, a previous clinical study also concluded that higher expression of CTHRC1 was associated with worse prognosis (Gu et al., 2014). Moreover, the results of our meta-analysis further confirmed that CTHRC1 overexpression was associated with worse OS. As for FAP, it is one of the active members of the S9b protease family and known to promote EMT of oral squamous cell carcinoma (OSCC) (Wu et al., 2020b; Huang et al., 2021). However, comprehensive research of FAP on GC is limited. Our finding of higher expression of FAP in GC and its association with worse OS suggests that more studies are needed to investigate the mechanism and effects of FAP in GC. It also has to be noted that the verification by Kaplan–Meier plotter only found marginal association with worse OS, and thus more original studies on its prognostic role in GC are also needed.
As CTHRC1 was found to have the highest HR among other up-regulated hub genes, a nomogram was created based on CTHRC1, age, gender, reflux history, Barrett’s esophagus, H. Pylori infection, pathologic stage, histologic grade, and resident tumor. This showed a relatively high accuracy of prediction given its C-index above 0.7. Thus, CTHRCI demonstrated its utility as not only a diagnostic biomarker, but also as one of the predictors in the nomogram developed in this study.
Further analysis of the up-regulated hub genes also identified the potential diagnostic markers found that ADAMTS2, COL10A1, COL11A1, and CTHRC1 as potential diagnostic markers given that their AUC value and 95% CI were all above 0.9. There were a few studies on the prognostic value of ADAMTS2 in GC, but the evidence of its diagnostic value on GC is limited (Jiang et al., 2019; Liang et al., 2020), while previous study showed that high levels of COL10A1 in plasma could provide diagnostic value in GC with AUC of 0.9171 (p = 0.0002) (Necula et al., 2020). As for COL11A1, a previous study suggested it could be used for differentiating the malignant lesions from premalignant tissues in stomach cancer based on 42 tissues samples (Zhao et al., 2009b), and our result further supports their conclusions though more studies are still warranted. Regarding CTHRC1, although its prognostic predicting value was more extensively studied, its diagnostic value has not been fully investigated and our results provide new evidence of its potential use in diagnosis of GC.
To further explore the potential drugs that may target the hub genes that are associated with worse prognosis, analysis of drug interactions was performed. For COL8A1, the compound or drug that showed the highest interaction was not investigated previously on GC. The second strongest interactor was lovastatin, an agonist of Src homology-2 domain-containing protein tyrosine phosphatase-2 (SHP2) that was found to significantly enhance the efficacy of chemotherapy in colon cancer (Wei et al., 2021). Other in vitro studies found that lovastatin inhibited gastric cancer cells (Cheng-Qian et al., 2014; Zhang et al., 2019b). As for COL10A1, although the top interacting compound, 2-amino-4-(2-hydroxy-4-methylphenyl)-5-phenylpyrimidine, was not previously investigated, the second drug on the list was sendanin. This compound was shown to inhibit the cancer cell lines according to previous in vitro study (Kim et al., 1994). In regards to CTHRC1, the first four compounds on the list were not well studied on cancer, while PF2341066, an inhibitor of anaplastic lymphoma kinase and c-Met later named as crizotinib, was previously shown to have antitumor activity of PF-2341066 in experimental models of anaplastic large-cell lymphoma. Its clinical use is increasingly reported in gastric cancer patients (Christensen et al., 2007; Sabree et al., 2018; Hou et al., 2019; Aparicio et al., 2021). Regarding FAP, although the first few drugs on the list were not well studied, kinetin riboside was found to inhibit colon cancer cells or even stimulated apoptosis (Cheong et al., 2009; Dudzik et al., 2011; Rajabi et al., 2012). Despite a number of studies investigating the compounds or drugs found in our analysis, more research studying on their effects on GC patients is still needed.
Further analysis on the relationship between COL8A1, COL10A1, CTHRC1, and FAP, and immune cell infiltrations in GC found that macrophages, neutrophils, and dendritic cells are positively correlated with these genes. Tumor-infiltrating macrophages were known to play a vital role in tumorigenesis by promoting tumor growth, migration, and invasion, as well as suppression of anti-tumor activity and progression (Van Overmeire et al., 2014; Rihawi et al., 2021). Moreover, tumor-associated macrophages (TAMs) were proposed to have the ability of significantly interfering with treatment response to chemotherapy, immune checkpoint inhibitors (ICIs), antiangiogenic drugs, and even radiotherapy or other treatment methods, leading to failure of treatment (Mantovani and Allavena, 2015; Qiu et al., 2018; Rihawi et al., 2021). Previous studies based on gastric cancer samples found increased neutrophil infiltration in GC, which is consistent with our finding (Kim et al., 2017). It was also shown that neutrophils activated by granulocyte-macrophage colony-stimulating factor (GM-CSF) could express CD54 and B7-H4 that are associated with reduced overall survival of GC patients following surgery (Shan et al., 2021). Tumor-associated neutrophils (TANs) were proposed to induce lymphangiogenesis and angiogenesis, and it was shown that local infiltration of certain types of TANs may play a role in the metastasis in GC (Hiramatsu et al., 2018). Infiltration of dendritic cells in gastric cancer and lymph nodes of GC patients are well known (Tsujitani et al., 1992; Tsujitani et al., 1995), and its infiltrating level in GC was also found to be closely correlated with macrophage infiltration (Xiang et al., 2020). One study found that higher level of dendritic cell infiltration was associated with longer OS (Higgins and Thompson, 2002; Hu et al., 2014). Given that our findings of the roles of the tumor-infiltrating immune cells and the selected up-regulated hub genes in the tumor microenvironment, further investigation and more comprehensive studies on the associations of tumor-infiltrating immune cells and these genes in GC are needed.
This bioinformatic study also has some limitations: First, all data were retrieved from online databases; therefore, the results need to be validated with other cohorts and experiments. Second, as this study mainly aimed to explore the potential clinical values of selected hub genes in the diagnosis and therapy of GC, the details of their mechanisms were not comprehensively explored, especially FAP with very limited number of previous studies of its effect on GC. Third, the TIMER database was mainly based on the TCGA database; therefore, the results need to be verified in the future with other cohorts and experiments.
Conclusion
In conclusion, among the nine up-regulated hub genes (COL1A2, THBS2, COL11A1, COL8A1, COL10A1, ADAMTS2, CTHRC1, FAP, and WISP1) and three down-regulated hub genes (TFF2, GIF, and CYP2C9), ADAMTS2, COL10A1, COL11A1, and CTHRC1 have demonstrated potential as diagnostic markers. COL8A1, COL10A1, CTHRC1, and FAP are associated with worse prognosis and could be potential prognostic biomarkers and therapeutic targets. The infiltration of macrophages, neutrophils, and dendritic cells are positively correlated with these COL8A1, COL10A1, CTHRC1, and FAP, suggesting the need for exploration of their roles in the tumor microenvironment of GC. Among these hub genes, CTHRC1 was found to have the highest prognostic accuracy and associated with worst prognosis as compared with other hub genes, and could be the next research hot spot. A nomogram based on CTHRC1 and other clinical predictors might be useful in clinical decision-making.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
YZ and CS designed the research study. CS, YC, and SM selected and collected the data. CS, YC, ZZ, SM, and Y-SC analyzed the data. CS, YC, NK, and SL wrote the manuscript. Y-SC, MWT, SM, WG, CB, JPWT, PT, CC, and QZ provided critical opinions and revised the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by The Training Program Foundation of Youth Scholars by the First Affiliated Hospital of Anhui Medical University (Grant Number 2020kj24).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Aparicio, T., Cozic, N., de la Fouchardière, C., Meriaux, E., Plaza, J., Mineur, L., et al. (2021). The Activity of Crizotinib in Chemo-Refractory MET-Amplified Esophageal and Gastric Adenocarcinomas: Results from the AcSé-Crizotinib Program. Targ. Oncol. 16 (3), 381–388. doi:10.1007/s11523-021-00811-8
Ariyawutyakorn, W., Saichaemchan, S., and Varella-Garcia, M. (2016). Understanding and Targeting MET Signaling in Solid Tumors - Are We There yet? J. Cancer 7, 633–649. doi:10.7150/jca.12663
Balachandran, V. P., Gonen, M., Smith, J. J., and DeMatteo, R. P. (2015). Nomograms in Oncology: More Than Meets the Eye. Lancet Oncol. 16 (4), e173–e180. doi:10.1016/s1470-2045(14)71116-7
Bang, Y.-J., Van Cutsem, E., Feyereislova, A., Chung, H. C., Shen, L., Sawaki, A., et al. (2010). 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 376, 687–697. doi:10.1016/s0140-6736(10)61121-x
Behjati, S., and Tarpey, P. S. (2013). What Is Next Generation Sequencing? Arch. Dis. Child. Educ. Pract. Ed. 98 (6), 236–238. doi:10.1136/archdischild-2013-304340
Borenstein, M., Hedges, L. V., Higgins, J. P. T., and Rothstein, H. R. (2010). A Basic Introduction to Fixed-Effect and Random-Effects Models for Meta-Analysis. Res. Synth. Method 1 (2), 97–111. doi:10.1002/jrsm.12
Chen, Q., Hu, L., and Chen, K. (2020). Construction of a Nomogram Based on a Hypoxia-Related lncRNA Signature to Improve the Prediction of Gastric Cancer Prognosis. Front. Genet. 11, 570325. doi:10.3389/fgene.2020.570325
Chen, S., Wei, Y., Liu, H., Gong, Y., Zhou, Y., Yang, H., et al. (2021). Analysis of Collagen Type X Alpha 1 (COL10A1) Expression and Prognostic Significance in Gastric Cancer Based on Bioinformatics. Bioengineered 12 (1), 127–137. doi:10.1080/21655979.2020.1864912
Cheng-Qian, Y., Xin-Jing, W., Wei-Xinbing, X., Zhuang-Lei, G., Hong-Peng, Z., Songde, X., et al. (2014). Lovastatin Inhibited the Growth of Gastric Cancer Cells. Hepatogastroenterology 61 (129), 1–4.
Cheong, J., Goh, D., Wan Hong Yong, J., Ngin Tan, S., and Shi Ong, E. (2009). Inhibitory Effect of Kinetin Riboside in Human Heptamoa, HepG2. Mol. Biosyst. 5 (1), 91–98. doi:10.1039/b712807j
Chin, C. H., Chen, S. H., Wu, H. H., Ho, C. W., Ko, M. T., and Lin, C. Y. (2014). cytoHubba: Identifying Hub Objects and Sub-networks from Complex Interactome. BMC Syst. Biol. 8 (Suppl. 4), S11. doi:10.1186/1752-0509-8-S4-S11
Choi, S., Park, S., Kim, H., Kang, S. Y., Ahn, S., and Kim, K.-M. (2022). Gastric Cancer: Mechanisms, Biomarkers, and Therapeutic Approaches. Biomedicines 10 (3), 543. doi:10.3390/biomedicines10030543
Christensen, J. G., Zou, H. Y., Arango, M. E., Li, Q., Lee, J. H., McDonnell, S. R., et al. (2007). Cytoreductive Antitumor Activity of PF-2341066, a Novel Inhibitor of Anaplastic Lymphoma Kinase and C-Met, in Experimental Models of Anaplastic Large-Cell Lymphoma. Mol. Cancer Ther. 6 (12 Pt 1), 3314–3322. doi:10.1158/1535-7163.MCT-07-0365
D'Errico, M., de Rinaldis, E., Blasi, M. F., Viti, V., Falchetti, M., Calcagnile, A., et al. (2009). Genome-wide Expression Profile of Sporadic Gastric Cancers with Microsatellite Instability. Eur. J. Cancer 45 (3), 461–469. doi:10.1016/j.ejca.2008.10.032
Dalkilic, S. (2020). Analysis of Gastric Cancer Transcriptomic Data by Bioinformatics Tools and Detection of Candidate Diagnostic Biomarker Genes. Prog. Nutr. 22 (2), e2020003.
Dan Song, Y. W., and Zhang, Jie (2021). Expression of BRSK2, ZFX and CTHRC1 Proteins in Gastric Cancer Tissues and its Relationship with Pathological Parameters and Prognosis [in Chineses]. Med. Pharm. J. Chin. People's Liberation Army 33 (02), 18–23.
DerSimonian, R., and Laird, N. (1986). Meta-analysis in Clinical Trials. Control. Clin. Trials 7 (3), 177–188. doi:10.1016/0197-2456(86)90046-2
Ding, X., Huang, R., Zhong, Y., Cui, N., Wang, Y., Weng, J., et al. (2020). CTHRC1 Promotes Gastric Cancer Metastasis via HIF-1α/CXCR4 Signaling Pathway. Biomed. Pharmacother. 123, 109742. doi:10.1016/j.biopha.2019.109742
Dong, Q., Li, F., Xu, Y., Xiao, J., Xu, Y., Shang, D., et al. (2020). RNAactDrug: a Comprehensive Database of RNAs Associated with Drug Sensitivity from Multi-Omics Data. Brief. Bioinform 21 (6), 2167–2174. doi:10.1093/bib/bbz142
Dudzik, P., Dulińska-Litewka, J., Wyszko, E., Jędrychowska, P., Opałka, M., Barciszewski, J., et al. (2011). Effects of Kinetin Riboside on Proliferation and Proapoptotic Activities in Human Normal and Cancer Cell Lines. J. Cell. Biochem. 112 (8), 2115–2124. doi:10.1002/jcb.23132
Ghidini, M., Petrillo, A., Botticelli, A., Trapani, D., Parisi, A., La Salvia, A., et al. (2021). How to Best Exploit Immunotherapeutics in Advanced Gastric Cancer: Between Biomarkers and Novel Cell-Based Approaches. J. Clin. Med. 10 (7). doi:10.3390/jcm10071412
Giunchi, F., Franceschini, T., and Fiorentino, M. (2021). A Narrative Review of Individualized Treatments of Genitourinary Tumors: Is the Future Brighter with Molecular Evaluations? Transl. Androl. Urol. 10 (3), 1553–1561. doi:10.21037/tau-20-1185
Gu, L., Liu, L., Zhong, L., Bai, Y., Sui, H., Wei, X., et al. (2014). Cthrc1 Overexpression Is an Independent Prognostic Marker in Gastric Cancer. Hum. Pathol. 45 (5), 1031–1038. doi:10.1016/j.humpath.2013.12.020
He, J., Jin, Y., Chen, Y., Yao, H.-B., Xia, Y.-J., Ma, Y.-Y., et al. (2016). Downregulation of ALDOB Is Associated with Poor Prognosis of Patients with Gastric Cancer. Ott 9, 6099–6109. doi:10.2147/ott.s110203
Hernandez, A. V., Marti, K. M., and Roman, Y. M. (2020). Meta-Analysis. Chest 158 (1S), S97–S102. doi:10.1016/j.chest.2020.03.003
Higgins, J. P. T., and Thompson, S. G. (2002). Quantifying Heterogeneity in a Meta-Analysis. Stat. Med. 21 (11), 1539–1558. doi:10.1002/sim.1186
Hiramatsu, S., Tanaka, H., Nishimura, J., Sakimura, C., Tamura, T., Toyokawa, T., et al. (2018). Neutrophils in Primary Gastric Tumors Are Correlated with Neutrophil Infiltration in Tumor-Draining Lymph Nodes and the Systemic Inflammatory Response. BMC Immunol. 19 (1), 13. doi:10.1186/s12865-018-0251-2
Hooi, J. K. Y., Lai, W. Y., Ng, W. K., Suen, M. M. Y., Underwood, F. E., Tanyingoh, D., et al. (2017). Global Prevalence of Helicobacter pylori Infection: Systematic Review and Meta-Analysis. Gastroenterology 153 (2), 420–429. doi:10.1053/j.gastro.2017.04.022
Hou, G., Song, B. B., and Song, B. (2019). Gastric Cancer Patient with C-MET Amplification Treated with Crizotinib after Failed Multi-Line Treatment: A Case Report and Literature Review. Math. Biosci. Eng. 16 (5), 5923–5930. doi:10.3934/mbe.2019296
Hu, M., Li, K., Maskey, N., Xu, Z., Peng, C., Wang, B., et al. (2014). Decreased Intratumoral Foxp3 Tregs and Increased Dendritic Cell Density by Neoadjuvant Chemotherapy Associated with Favorable Prognosis in Advanced Gastric Cancer. Int. J. Clin. Exp. Pathol. 7 (8), 4685–4694.
Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2009). Bioinformatics Enrichment Tools: Paths toward the Comprehensive Functional Analysis of Large Gene Lists. Nucleic Acids Res. 37 (1), 1–13. doi:10.1093/nar/gkn923
Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2009). Systematic and Integrative Analysis of Large Gene Lists Using DAVID Bioinformatics Resources. Nat. Protoc. 4 (1), 44–57. doi:10.1038/nprot.2008.211
Huang, J. C., Emran, A. A., Endaya, J. M., McCaughan, G. W., Gorrell, M. D., and Zhang, H. E. (2021). DPP9: Comprehensive In Silico Analyses of Loss of Function Gene Variants and Associated Gene Expression Signatures in Human Hepatocellular Carcinoma. Cancers (Basel) 13 (7). doi:10.3390/cancers13071637
Iasonos, A., Schrag, D., Raj, G. V., and Panageas, K. S. (2008). How to Build and Interpret a Nomogram for Cancer Prognosis. Jco 26 (8), 1364–1370. doi:10.1200/jco.2007.12.9791
Jiang, C., Zhou, Y., Huang, Y., Wang, Y., Wang, W., and Kuai, X. (2019). Overexpression of ADAMTS-2 in Tumor Cells and Stroma Is Predictive of Poor Clinical Prognosis in Gastric Cancer. Hum. Pathol. 84, 44–51. doi:10.1016/j.humpath.2018.08.030
Jin, Y., He, J., Du, J., Zhang, R. X., Yao, H. B., and Shao, Q. S. (2017). Overexpression of HS6ST2 Is Associated with Poor Prognosis in Patients with Gastric Cancer. Oncol. Lett. 14 (5), 6191–6197. doi:10.3892/ol.2017.6944
Kim, H. M., Oh, G. T., Han, S. B., Hong, D. H., Hwang, B. Y., Kim, Y. H., et al. (1994). Comparative Studies of Adriamycin and 28-Deacetyl Sendanin Onin Vitro Growth Inhibition of Human Cancer Cell Lines. Arch. Pharm. Res. 17 (2), 100–103. doi:10.1007/bf02974231
Kim, J.-W., Lee, H. S., Nam, K. H., Ahn, S., Kim, J. W., Ahn, S.-H., et al. (2017). PIK3CA Mutations Are Associated with Increased Tumor Aggressiveness and Akt Activation in Gastric Cancer. Oncotarget 8 (53), 90948–90958. doi:10.18632/oncotarget.18770
Kramer, A. A., and Zimmerman, J. E. (2007). Assessing the Calibration of Mortality Benchmarks in Critical Care: The Hosmer-Lemeshow Test Revisited*. Crit. Care Med. 35 (9), 2052–2056. doi:10.1097/01.ccm.0000275267.64078.b0
Levy, S. E., and Myers, R. M. (2016). Advancements in Next-Generation Sequencing. Annu. Rev. Genom. Hum. Genet. 17, 95–115. doi:10.1146/annurev-genom-083115-022413
Li, B., Severson, E., Pignon, J.-C., Zhao, H., Li, T., Novak, J., et al. (2016). Comprehensive Analyses of Tumor Immunity: Implications for Cancer Immunotherapy. Genome Biol. 17 (1), 174. doi:10.1186/s13059-016-1028-7
Li, H. H., Wang, J. D., Wang, W., Wang, H. F., and Lv, J. Q. (2020). Effect of miR-26a-5p on Gastric Cancer Cell Proliferation, Migration and Invasion by Targeting COL10A1. Eur. Rev. Med. Pharmacol. Sci. 24 (3), 1186–1194. doi:10.26355/eurrev_202002_20170
Li, S., Yang, R., Sun, X., Miao, S., Lu, T., Wang, Y., et al. (2018). Identification of SPP1 as a Promising Biomarker to Predict Clinical Outcome of Lung Adenocarcinoma Individuals. Gene 679, 398–404. doi:10.1016/j.gene.2018.09.030
Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 77 (21), e108–e110. doi:10.1158/0008-5472.can-17-0307
Li, T., Huang, H., Shi, G., Zhao, L., Li, T., Zhang, Z., et al. (2018). TGF-β1-SOX9 Axis-Inducible COL10A1 Promotes Invasion and Metastasis in Gastric Cancer via Epithelial-To-Mesenchymal Transition. Cell. Death Dis. 9 (9), 849. doi:10.1038/s41419-018-0877-2
Liang, L., Zhu, J. H., Chen, G., Qin, X. G., and Chen, J. Q. (2020). Prognostic Values for the mRNA Expression of the ADAMTS Family of Genes in Gastric Cancer. J. Oncol. 2020, 9431560. doi:10.1155/2020/9431560
Liu, J., Lichtenberg, T., Hoadley, K. A., Poisson, L. M., Lazar, A. J., Cherniack, A. D., et al. (2018). An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics. Cell. 173 (2), 400–e11. doi:10.1016/j.cell.2018.02.052
Liu, S., Geng, R., Lin, E., Zhao, P., and Chen, Y. (2021). ERBB1/2/3 Expression, Prognosis, and Immune Infiltration in Cutaneous Melanoma. Front. Genet. 12, 602160. doi:10.3389/fgene.2021.602160
Liu, X., Wu, J., Zhang, D., Bing, Z., Tian, J., Ni, M., et al. (2018). Identification of Potential Key Genes Associated with the Pathogenesis and Prognosis of Gastric Cancer Based on Integrated Bioinformatics Analysis. Front. Genet. 9, 265. doi:10.3389/fgene.2018.00265
Lu, X.-Q., Zhang, J.-Q., Zhang, S.-X., Qiao, J., Qiu, M.-T., Liu, X.-R., et al. (2021). Identification of Novel Hub Genes Associated with Gastric Cancer Using Integrated Bioinformatics Analysis. BMC Cancer 21 (1), 697. doi:10.1186/s12885-021-08358-7
Mantovani, A., and Allavena, P. (2015). The Interaction of Anticancer Therapies with Tumor-Associated Macrophages. J. Exp. Med. 212 (4), 435–445. doi:10.1084/jem.20150295
Mei, D., Zhu, Y., Zhang, L., and Wei, W. (2020). The Role of CTHRC1 in Regulation of Multiple Signaling and Tumor Progression and Metastasis. Mediat. Inflamm. 2020, 9578701. doi:10.1155/2020/9578701
Nagata, N., Akiyama, J., Marusawa, H., Shimbo, T., Liu, Y., Igari, T., et al. (2014). Enhanced Expression of Activation-Induced Cytidine Deaminase in Human Gastric Mucosa Infected by Helicobacter pylori and its Decrease Following Eradication. J. Gastroenterol. 49 (3), 427–435. doi:10.1007/s00535-013-0808-z
Nakanishi, R., Shimizu, T., Kumagai, K., Takai, A., and Marusawa, H. (2021). Genetic Pathogenesis of Inflammation-Associated Cancers in Digestive Organs. Pathogens 10 (4). doi:10.3390/pathogens10040453
Necula, L., Matei, L., Dragu, D., Pitica, I., Neagu, A. I., Bleotu, C., et al. (2020). High Plasma Levels of COL10A1 Are Associated with Advanced Tumor Stage in Gastric Cancer Patients. Wjg 26 (22), 3024–3033. doi:10.3748/wjg.v26.i22.3024
Niu, X., Ren, L., Hu, A., Zhang, S., and Qi, H. (2022). Identification of Potential Diagnostic and Prognostic Biomarkers for Gastric Cancer Based on Bioinformatic Analysis. Front. Genet. 13, 862105. doi:10.3389/fgene.2022.862105
Page, M. J., McKenzie, J. E., Bossuyt, P. M., Boutron, I., Hoffmann, T. C., Mulrow, C. D., et al. (2021). The PRISMA 2020 Statement: An Updated Guideline for Reporting Systematic Reviews. J. Clin. Epidemiol. doi:10.1186/s13643-021-01626-4
Palrasu, M., Zaika, E., El-Rifai, W., Que, J., and Zaika, A. I. (2021). Role of Bacterial and Viral Pathogens in Gastric Carcinogenesis. Cancers (Basel) 13 (8). doi:10.3390/cancers13081878
Pardoll, D. M. (2012). The Blockade of Immune Checkpoints in Cancer Immunotherapy. Nat. Rev. Cancer 12, 252–264. doi:10.1038/nrc3239
Pencina, M. J., and D'Agostino, R. B. (2004). OverallC as a Measure of Discrimination in Survival Analysis: Model Specific Population Value and Confidence Interval Estimation. Stat. Med. 23 (13), 2109–2123. doi:10.1002/sim.1802
Peng, W., Li, J.-D., Zeng, J.-J., Zou, X.-P., Tang, D., Tang, W., et al. (2020). Clinical Value and Potential Mechanisms of COL8A1 Upregulation in Breast Cancer: a Comprehensive Analysis. Cancer Cell. Int. 20, 392. doi:10.1186/s12935-020-01465-8
Pontén, F., Schwenk, J. M., Asplund, A., and Edqvist, P.-H. D. (2011). The Human Protein Atlas as a Proteomic Resource for Biomarker Discovery. J. Intern Med. 270 (5), 428–446. doi:10.1111/j.1365-2796.2011.02427.x
Qiu, S.-Q., Waaijer, S. J. H., Zwager, M. C., de Vries, E. G. E., van der Vegt, B., and Schröder, C. P. (2018). Tumor-associated Macrophages in Breast Cancer: Innocent Bystander or Important Player? Cancer Treat. Rev. 70, 178–189. doi:10.1016/j.ctrv.2018.08.010
Rajabi, M., Gorincioi, E., and Santaniello, E. (2012). Antiproliferative Activity of Kinetin Riboside on HCT-15 Colon Cancer Cell Line. Nucleosides, Nucleotides Nucleic Acids 31 (6), 474–481. doi:10.1080/15257770.2012.681825
Research. WCRFAIfC, (2018). “The Continuous Update Project Expert Report,” in Diet, Nutrition, Physical Activity and Cancer: Colorectal Cancer.
Rihawi, K., Ricci, A. D., Rizzo, A., Brocchi, S., Marasco, G., Pastore, L. V., et al. (2021). Tumor-Associated Macrophages and Inflammatory Microenvironment in Gastric Cancer: Novel Translational Implications. Int. J. Mol. Sci. 22 (8). doi:10.3390/ijms22083805
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res. 43 (7), e47. doi:10.1093/nar/gkv007
Roh, C. K., Choi, S., Seo, W. J., Cho, M., Son, T., Kim, H.-I., et al. (2020). Indocyanine Green Fluorescence Lymphography during Gastrectomy after Initial Endoscopic Submucosal Dissection for Early Gastric Cancer. Br. J. Surg. 107 (6), 712–719. doi:10.1002/bjs.11438
Sabree, S., Berg, D., and Sato, M. (2018). Treatment of a Pediatric Patient with MET -amplified Signet Ring Cell Adenocarcinoma of the Stomach with Crizotinib. Pediatr. Blood Cancer 65 (6), e26984. doi:10.1002/pbc.26984
Schreiber, R. D., Old, L. J., and SmythCancer immunoediting, M. J. (2011). Cancer Immunoediting: Integrating Immunity's Roles in Cancer Suppression and Promotion. Science 331, 1565–1570. doi:10.1126/science.1203486
Shan, Z. G., Yan, Z. B., Peng, L. S., Cheng, P., Teng, Y. S., Mao, F. Y., et al. (2021). Granulocyte-Macrophage Colony-Stimulating Factor-Activated Neutrophils Express B7-H4 that Correlates with Gastric Cancer Progression and Poor Patient Survival. J. Immunol. Res. 2021, 6613247. doi:10.1155/2021/6613247
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13 (11), 2498–2504. doi:10.1101/gr.1239303
Shen, Y., Liu, J., Zhang, L., Dong, S., Zhang, J., Liu, Y., et al. (2019). Identification of Potential Biomarkers and Survival Analysis for Head and Neck Squamous Cell Carcinoma Using Bioinformatics Strategy: A Study Based on TCGA and GEO Datasets. Biomed. Res. Int. 2019, 7376034. doi:10.1155/2019/7376034
Shirdarreh, M., Aziza, O., Pezo, R. C., Jerzak, K. J., and Warner, E. (2021). Patients' and Oncologists' Knowledge and Expectations Regarding Tumor Multigene Next-Generation Sequencing: A Narrative Review. Oncologist. doi:10.1002/onco.13783
Song, Y., Wang, X., Wang, F., Peng, X., Li, P., Liu, S., et al. (2021). Identification of Four Genes and Biological Characteristics of Esophageal Squamous Cell Carcinoma by Integrated Bioinformatics Analysis. Cancer Cell. Int. 21 (1), 123. doi:10.1186/s12935-021-01814-1
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Szász, A. M., Lánczky, A., Nagy, Á., Förster, S., Hark, K., Green, J. E., et al. (2016). Cross-validation of Survival Associated Biomarkers in Gastric Cancer Using Transcriptomic Data of 1,065 Patients. Oncotarget 7 (31), 49322–49333. doi:10.18632/oncotarget.10337
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING V11: Protein-Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-wide Experimental Datasets. Nucleic Acids Res. 47 (D1), D607–D613. doi:10.1093/nar/gky1131
Tan, Z. (2019). Recent Advances in the Surgical Treatment of Advanced Gastric Cancer: A Review. Med. Sci. Monit. 25, 3537–3541. doi:10.12659/msm.916475
Tsujitani, S., Kakeji, Y., Orita, H., Watanabe, A., Kohnoe, S., Baba, H., et al. (1992). Postoperative Adjuvant Immunochemotherapy and Infiltration of Dendritic Cells for Patients with Advanced Gastric Cancer. Anticancer Res. 12 (3), 645–648.
Tsujitani, S., Oka, A., Kondo, A., Ikeguchi, M., Maeta, M., and Kaibara, N. (1995). Infiltration of Dendritic Cells into Regional Lymph Nodes in Gastric Cancer. Cancer 75 (6 Suppl. l), 1478–1483. doi:10.1002/1097-0142(19950315)75:6+<1478::aid-cncr2820751516>3.0.co;2-o
Van Overmeire, E., Laoui, D., Keirsse, J., Van Ginderachter, J. A., and Sarukhan, A. (2014). Mechanisms Driving Macrophage Diversity and Specialization in Distinct Tumor Microenvironments and Parallelisms with Other Tissues. Front. Immunol. 5, 127. doi:10.3389/fimmu.2014.00127
Vlasblom, J., Zuberi, K., Rodriguez, H., Arnold, R., Gagarinova, A., Deineko, V., et al. (2015). Novel Function Discovery with GeneMANIA: a New Integrated Resource for Gene Function Prediction in Escherichia coli. Bioinformatics 31 (3), 306–310. doi:10.1093/bioinformatics/btu671
Wang, F., Xue, Q., Xu, D., Jiang, Y., Tang, C., and Liu, X. (2020). Identifying the Hub Gene in Gastric Cancer by Bioinformatics Analysis and In Vitro Experiments. Cell. Cycle 19 (11), 1326–1337. doi:10.1080/15384101.2020.1749789
Wang, H. (2016). Expression of CTHRC1, VEGF-C in Gastric Carcinoma. Jingzhou: Yangtze University. [in Chinese] MSc diss.
Wang, P., Wang, Y.-C., Chen, X.-Y., Shen, Z.-Y., Cao, H., Zhang, Y.-J., et al. (2012). CTHRC1is Upregulated by Promoter Demethylation and Transforming Growth Factor-Β1 and May Be Associated with Metastasis in Human Gastric Cancer. Cancer Sci. 103 (7), 1327–1333. doi:10.1111/j.1349-7006.2012.02292.x
Wang, Q., Wen, Y.-G., Li, D.-P., Xia, J., Zhou, C.-Z., Yan, D.-W., et al. (2012). Upregulated INHBA Expression Is Associated with Poor Survival in Gastric Cancer. Med. Oncol. 29 (1), 77–83. doi:10.1007/s12032-010-9766-y
Wang, T., Feng, Y., Wang, H., Huo, G., Cai, Y., Wang, L., et al. (2020). The Mechanisms of Sijunzi Decoction in the Treatment of Chronic Gastritis Revealed by Network Pharmacology. Evid. Based Complement. Altern. Med. 2020, 8850259. doi:10.1155/2020/8850259
Wang, W., He, Y., Zhao, Q., Zhao, X., and Li, Z. (2020). Identification of Potential Key Genes in Gastric Cancer Using Bioinformatics Analysis. Biomed. Rep. 12 (4), 178–192. doi:10.3892/br.2020.1281
Wang, Y. (2016). Expression of CTHRC1 in Digestive System Neoplasms and its Prognostic Value in Gastric Carcinoma. Wuhu: MSc diss, WanNan Medical College. [in Chinese].
Wang, Y., Zhang, L., Yang, Y., Lu, S., and Chen, H. (2021). Progress of Gastric Cancer Surgery in the Era of Precision Medicine. Int. J. Biol. Sci. 17 (4), 1041–1049. doi:10.7150/ijbs.56735
Wei, B., Xu, L., Guo, W., Wang, Y., Wu, J., Li, X., et al. (2021). SHP2-mediated Inhibition of DNA Repair Contributes to cGAS-STING Activation and Chemotherapeutic Sensitivity in Colon Cancer. Cancer Res. doi:10.1158/0008-5472.can-20-3738
Wu, K.-Z., Xu, X.-H., Zhan, C.-P., Li, J., and Jiang, J.-L. (2020). Identification of a Nine-Gene Prognostic Signature for Gastric Carcinoma Using Integrated Bioinformatics Analyses. Wjgo 12 (9), 975–991. doi:10.4251/wjgo.v12.i9.975
Wu, Q.-q., Zhao, M., Huang, G.-z., Zheng, Z.-n., Chen, Y., Zeng, W.-s., et al. (2020). Fibroblast Activation Protein (FAP) Overexpression Induces Epithelial-Mesenchymal Transition (EMT) in Oral Squamous Cell Carcinoma by Down-Regulating Dipeptidyl Peptidase 9 (DPP9). Ott 13, 2599–2611. doi:10.2147/ott.s243417
Xiang, A., Lin, X., Xu, L., Chen, H., Guo, J., and Zhou, F. (2020). PCOLCE Is Potent Prognostic Biomarker and Associates with Immune Infiltration in Gastric Cancer. Front. Mol. Biosci. 7, 544895. doi:10.3389/fmolb.2020.544895
Yang, H., Wu, J., Zhang, J., Yang, Z., Jin, W., Li, Y., et al. (2019). Integrated Bioinformatics Analysis of Key Genes Involved in Progress of Colon Cancer. Mol. Genet. Genomic Med. 7 (4), e00588. doi:10.1002/mgg3.588
Yang, W., Wu, X., and Zhou, F. (2021). Collagen Type X Alpha 1 (COL10A1) Contributes to Cell Proliferation, Migration, and Invasion by Targeting Prolyl 4-Hydroxylase Beta Polypeptide (P4HB) in Breast Cancer. Med. Sci. Monit. 27, e928919.
Yu, J., Feng, J., Zhi, X., Tang, J., Li, Z., Xu, Y., et al. (2015). Let-7b Inhibits Cell Proliferation, Migration, and Invasion through Targeting Cthrc1 in Gastric Cancer. Tumor Biol. 36 (5), 3221–3229. doi:10.1007/s13277-014-2950-5
Zhang, H., Guo, L., Zhang, Z., Sun, Y., Kang, H., Song, C., et al. (2019). Co-Expression Network Analysis Identified Gene Signatures in Osteosarcoma as a Predictive Tool for Lung Metastasis and Survival. J. Cancer 10 (16), 3706–3716. doi:10.7150/jca.32092
Zhang, L., Kang, W., Lu, X., Ma, S., Dong, L., and Zou, B. (2019). Weighted Gene Co-expression Network Analysis and Connectivity Map Identifies Lovastatin as a Treatment Option of Gastric Cancer by Inhibiting HDAC2. Gene 681, 15–25. doi:10.1016/j.gene.2018.09.040
Zhang, M.-M., Wang, D., Lu, F., Zhao, R., Ye, X., He, L., et al. (2021). Identification of the Active Substances and Mechanisms of Ginger for the Treatment of Colon Cancer Based on Network Pharmacology and Molecular Docking. BioData Min. 14 (1), 1. doi:10.1186/s13040-020-00232-9
Zhang, Y., Xia, M., Jin, K., Wang, S., Wei, H., Fan, C., et al. (2018). Function of the C-Met Receptor Tyrosine Kinase in Carcinogenesis and Associated Therapeutic Opportunities. Mol. Cancer. 17, 45. doi:10.1186/s12943-018-0796-y
Zhao, Y., Jia, L., Mao, X., Xu, H., Wang, B., and Liu, Y. (2009). siRNA-Targeted COL8A1 Inhibits Proliferation, Reduces Invasion and Enhances Sensitivity to D-Limonence Treatment in Hepatocarcinoma Cells. IUBMB Life 61 (1), 74–79. doi:10.1002/iub.151
Zhao, Y., Zhou, T., Li, A., Yao, H., He, F., Wang, L., et al. (2009). A Potential Role of Collagens Expression in Distinguishing between Premalignant and Malignant Lesions in Stomach. Anat. Rec. 292 (5), 692–700. doi:10.1002/ar.20874
Zhong, L., Yang, Z., Lei, D., Li, L., Song, S., Cao, D., et al. (2021). Bromodomain 4 Is a Potent Prognostic Marker Associated with Immune Cell Infiltration in Breast Cancer. Basic Clin. Pharmacol. Toxicol. 128 (1), 169–182. doi:10.1111/bcpt.13481
Zhou, J., Song, Y., Gan, W., Liu, L., Chen, G., Chen, Z., et al. (2020). Upregulation of COL8A1 Indicates Poor Prognosis across Human Cancer Types and Promotes the Proliferation of Gastric Cancer Cells. Oncol. Lett. 20 (4), 34. doi:10.3892/ol.2020.11895
Zhu, H., Hu, X., Feng, S., Jian, Z., Xu, X., Gu, L., et al. (2022). The Hypoxia-Related Gene COL5A1 Is a Prognostic and Immunological Biomarker for Multiple Human Tumors. Oxid. Med. Cell. Longev. 2022, 6419695. doi:10.1155/2022/6419695
Keywords: gastric cancer, hub genes, bioinformatics analysis, meta-analysis, ROC (receiver operating curve), prognostic value, immune infiltration
Citation: Sun C, Chen Y, Kim NH, Lowe S, Ma S, Zhou Z, Bentley R, Chen Y-S, Tuason MW, Gu W, Bhan C, Tuason JPW, Thapa P, Cheng C, Zhou Q and Zhu Y (2022) Identification and Verification of Potential Biomarkers in Gastric Cancer By Integrated Bioinformatic Analysis. Front. Genet. 13:911740. doi: 10.3389/fgene.2022.911740
Received: 12 April 2022; Accepted: 08 June 2022;
Published: 15 July 2022.
Edited by:
Qiu-Ning Liu, Yancheng Teachers University, ChinaReviewed by:
Tung Nguyen-Thanh, Hue University of Medicine and Pharmacy, VietnamShihai Liu, The Affiliated Hospital of Qingdao University, China
Copyright © 2022 Sun, Chen, Kim, Lowe, Ma, Zhou, Bentley, Chen, Tuason, Gu, Bhan, Tuason, Thapa, Cheng, Zhou and Zhu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yanzhe Zhu, ZHJ6aHV5YW56aGVAMTI2LmNvbQ==
†These authors have contributed equally to this work