BioMed Research International
Volume 2020 (2020), Article ID 7658230, 12 pages
https://doi.org/10.1155/2020/7658230

Research Article

Identification of Key Genes in Gastric Cancer by Bioinformatics Analysis

Xinyu Chong,¹ Rui Peng,² Yan Sun,¹ Luyu Zhang,¹ and Zheng Zhang¹

¹Department of Molecular Medicine and Cancer Research Center, Chongqing Medical University, Chongqing, China
²Department of Bioinformatics, Chongqing Medical University, Chongqing, China

Correspondence should be addressed to Zheng Zhang; zhangzheng92@163.com

Received 16 May 2020; Revised 12 August 2020; Accepted 1 September 2020; Published 21 September 2020

Academic Editor: Costantino Del Gaudio

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

Abstract.

Gastric cancer (GC) is one of the most common malignancies of the digestive system with few genetic markers for its early detection and prevention. In this study, differentially expressed genes (DEGs) were analyzed using GEO2R from GSE54129 and GSE13911 of the Gene Expression Omnibus (GEO). Then, gene enrichment analysis, protein-protein interaction (PPI) network construction, and topological analysis were performed on the DEGs by the Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, STRING, and Cytoscape. Finally, we performed survival analysis of key genes through the Kaplan-Meier plotter. A total of 1034 DEGs were identified in GC. GO and KEGG results showed that DEGs mainly enriched in plasma membrane, cell adhesion, and PI3K-Akt signaling pathway. Subsequently, the PPI network with 44 nodes and 333 edges was constructed, and 18 candidate genes in the network were focused on by centrality analysis and module analysis. Furthermore, data showed that high expressions of fibronectin 1(FN1), the tissue inhibitor of metalloproteinases 1 (TIMP1), secreted phosphoprotein 1 (SPP1), apolipoprotein E (APOE), and versican (VCAN) were related to poor overall survivals in GC patients. In summary, this study suggests that FN1, TIMP1, SPP1, APOE, and VCAN may act as the key genes in GC.

1. Introduction

Gastric cancer (GC) is one of the malignant tumors threatening human health, and it is the fifth most common cancer and the third leading cause of cancer death in the world [1, 2]. The cause of the onset of it includes diet, lifestyle, genetic, and infection of helicobacter pylori[3–5]. At present, the treatment of gastric cancer includes surgery [6], radiotherapy [7], neoadjuvant chemotherapy [8], and immunotherapy [9]. The survival rate for early gastric cancer patients can reach 90%. However, it is difficult to detect and diagnose the disease during the early period, resulting in a significant decrease in survival after diagnosis [10]. Therefore, it is vitally important to explore potential diagnostic and prognostic biomarkers and therapeutic targets of early GC.

With the advance of the human genome project, cancer has been studied at the genetic level. Gene chips can be used to identify genes that cause early cancer. It has the characteristics of high flux, high sensitivity, and low cost. It is widely used in disease diagnosis and drug screening [11, 12]. At present, DNA microarrays were used to identify potential biomarkers that affect the development of diseases in studies [13]. Gene microarray technology plays an important role in elucidating cell cycle-related genes and evaluating their expression in gastric cancer [14]. However, the pathogenesis of GC remains unclear.

In this study, we aimed to obtain the key genes between GC patients and normal controls. We downloaded the gene expression profiles of GSE54129 and GSE13911 and identified 1034 differentially expressed genes (DEGs) in GC. Moreover, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were carried out for DEGs, and GC-related protein-protein interaction (PPI) network was constructed. Furthermore, 18 candidate genes displayed high centrality values and located at the 1st module, which were found by the centrality analysis and module analysis on the basis of the GC-related PPI network. Data also showed that high expressions of fibronectin 1 (FN1), the tissue inhibitor of metalloproteinases 1 (TIMP1), secreted phosphoprotein 1 (SPP1), apolipoprotein E (APOE), and versican (VCAN) were related to a poor overall survival in gastric cancer patients. These key genes could be used as potential therapeutic targets and biomarkers for gastric cancer at early period.

2. Materials and Methods

2.1. Microarray Data

The Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/geo) is a public genome database [15]. In this study, two gene expression profiles GSE54129 and GSE13911 were downloaded from the GEO. The criteria for both gene expression profiles were (a) the samples included two groups of GC tumors and normal tissues; (b) the sample size of each gene expression profile was greater than 60; (c) they were recently updated in the last two years (2019-2020); and (d) they were derived from the same platform: GPL570 [hg-u133_plus_2] Affymetrix Human Genome U133 Plus 2.0 Array. The GSE54129 data profile contained 111 tumor tissues and 21 normal tissues. The GSE13911 data profile contained 38 tumor tissues and 31 normal tissues.

2.2. Identification of DEGs

GEO2R (https://www.ncbi.nlm.nih.gov/geo/geo2r/) online analysis software was used to analyze the differentially expressed genes in tumor and normal samples of GC. The intersection of DEGs was obtained by Venny 2.1.0 (https://bioinfogp.cnb.csic.es/tools/venny/) [16]. The adj. value was obtained by the Benjamini and Hochberg method to control the probability of false positives. Adj. value < 0.01 and were used as the cut-off criteria [17].

2.3. GO and KEGG Enrichment Analysis of DEGs

The DAVID (https://david.ncifcrf.gov/) is a database that provides systematic and comprehensive annotation information of biological functions for genes and proteins [18]. GO is an international standard classification system of gene function [19]. It includes biological process (BP), cellular component (CC), and molecular function (MF). KEGG is a database for systematic analysis of gene function and genomic information [20]. We used DAVID to carry out GO functional annotation and KEGG pathway enrichment analysis on the DEGs. were considered statistically significant [21].

2.4. PPI Network Construction

The STRING (http://string-db.org) is a database for searching between known proteins and predicting the interactions between proteins [22]. We used it to build PPI network for DEGs. The combined was considered as the cut-off value [23].

2.5. Module Analysis of PPI Network

The Cytoscape is a software for visual networks [24]. Functional modules in the network were identified by using the plug-in MCODE of Cytoscape [25]. The selection criteria were as follows: , , , and . Submodules were sorted by score value. The higher the score was, the stronger the protein correlation in the module was.

2.6. Centrality Analysis of PPI Network

Centrality analysis includes analyzing the degree, betweenness, and eigenvector of network nodes. Cytoscape plug-in CytoNCA was used to calculate the values of degree, betweenness, and eigenvector to predict the key genes [26]. Degree centrality is a measure of the importance of a single node, and it describes the number of edges connecting nodes [27]. Betweenness centrality is the shortest path through which a particular node is analyzed [28]. Eigenvector centrality takes into account the degree of itself and the degree of its neighbors [29]. The distribution characteristics of degree, betweenness, and eigenvector were determined by density analysis. The correlation between degree and betweenness, between degree and eigenvector, and between betweenness and eigenvector was calculated by taking the top 5% of the three parameters. And it was visualized by R language.

2.7. Survival Analysis of Hub Genes

The Kaplan-Meier plotter (KM plotter, http://kmplot.com/analysis/) is an online website survival analysis [30]. There are 21 types of cancer, including breast cancer, ovarian cancer, lung cancer, and gastric cancer. The online analysis database includes GEO, EGA, and TCGA. The KM plotter was used to reflect the survival expression of hub genes. In the KM plotter, we selected patient grouping criteria to automatically select the optimal cutoff, overall survival, and complete follow-up time for systematic and complete survival analysis. was set as the cut-off criterion [31].

3. Results

3.1. Identification of DEGs

To explore the role of systems biology in the pathogenesis of GC, we analyzed two chip data of GSE54129 and GSE13911 by GEO2R. There were 3878 differentially expressed genes in GSE54129 and 3061 differentially expressed genes in GSE13911. Then, venny2.1.0 was used to obtain the intersection of the DEGs of the two chips. The results showed that 1034 differentially expressed genes appeared on both chips (Figure 1(a)). These DEGs included 403 upregulated genes and 631 downregulated genes between tumor and nontumor samples (Figure 1(b)). These data provided basic data for further analysis.

(a)

(b)

Figure 1: Screening and identification of differentially expressed genes. (a) Venn diagram showed the differentially expressed genes of adj. value < 0.01 and . (b) Red points meaned upregulated genes screened on the basis and value < 0.01. Blue points meaned downregulated genes screened on the basis and value < 0.01.

3.2. GO and KEGG Enrichment Analysis of DEGs

In order to better understand the biological function of DEGs, we conducted GO function and KEGG enrichment analysis by DAVID. GO results showed that DEGs significantly enriched in extracellular matrix organization, collagen catabolic process and cell adhesion of BP, extracellular space, extracellular region and extracellular exosome of CC, extracellular matrix structural constituent, and heparin binding and integrin binding of MF (Figure 2(a)). Moreover, KEGG analysis showed that the DEGs were enriched in ECM-receptor interaction, PI3K-Akt signaling pathway and focal adhesion, and so on (Figure 2(b)).

(a)

(b)

Figure 2: GO and KEGG enrichment analysis of the PPI network. (a) Top 25 significantly enriched gene ontology terms, including three groups (biological process, cellular component, and molecular function), . (b) The KEGG pathway analysis of DEGs in GC, .

3.3. Construction of GC-Related PPI Network

To study the molecular mechanism of gastric cancer from a systematic perspective, PPI network was constructed to explore the relationship between proteins. PPI network was constructed by STRING for DEGs with a confidence level of >0.4. The result of network analysis showed that PPI enrichment value < 1.0-16. There were 865 nodes and 4483 edges in the visualization network using the Cyctoscape.

3.4. Module Analysis of the PPI Network

In order to explore more closely related genes in the complex PPI network, we conducted module analysis of the network by MCODE. The result showed that there were 27 modules in PPI network. We found that the first module was the most densely interacted region in PPI network, with a score of 15.488. The module was located at the center of the entire network, including 44 nodes and 333 edges (Figure 3). The above results suggested that the protein association in the first-rank module may be the strongest and the most significant part of the whole network.

Figure 3: Module analysis of PPI networks obtained through Cytoscape’s plug-in MCODE. The most prominent module in the PPI network included 44 nodes and 333 edges.

3.5. Centrality Analysis of PPI Networks

To analyze the key genes in the complex PPI network, we used the centrality analysis to analyze them. First, we analyzed the subcases of the three parameters by their density. The results showed that degree, betweenness, and eigenvector were power-law distributions (Figure 4). Then, we took the top 5% of the three parameters and analyze their correlation. The results showed that the correlation coefficient between degree and betweenness was 0.793, the correlation coefficient between degree and eigenvector was 0.920, and the correlation coefficient between betweenness and eigenvector was 0.620 (Figures 5(a)–5(c)). The results indicated that the degree, betweenness, and eigenvector are positively correlated to each other and have significant correlations. Finally, we studied the top 5% of the genes of each parameter and obtained 18 genes with degree, betweenness, and eigenvector by taking the intersection (Figure 5(d), Table 1). Combined these with the results of module analysis, FN1, TIMP1, SPP1, apolipoprotein B (APOB), APOE, VCAN, and complement C3 (C3) were focused on because these seven genes with high centrality values were located in the first-rank module.

(a)

(b)

(c)

Figure 4: Centrality analysis of PPI networks obtained through Cytoscape’s plug-in CytoNCA. (a) A density diagram of degree centrality. (b) A density diagram of betweenness centrality. (c) A density diagram of eigenvector centrality.

(a)

(b)

(c)

(d)

Figure 5: Correlation analysis of the top 5% of molecules of each centrality (degree, betweenness, and eigenvector). (a) The correlation coefficient between degree and betweenness was 0.793. (b) The correlation coefficient between degree and eigenvector was 0.920. (c) The correlation coefficient between betweenness and eigenvector was 0.620. (d) Venny2.1.0 was used to obtain the intersection of the top 5% of genes of each centrality (degree, betweenness, and eigenvector). The results showed that the 18 keys were further studied because of their high degree, betweenness, and eigenvector values.

Table 1: Top 5% of candidate genes in the centrality analysis.

Gene
Degree
Betweenness
Eigenvector

FN1
127
77922.220
0.253395320

MMP9
93
45189.664
0.198073490

CXCL8
85
33816.297
0.146214620

CD44
84
41183.110
0.172110660

MYC
82
48818.060
0.120246940

CXCL12
70
19844.107
0.157128860

TIMP1
65
8912.090
0.174682420

PTGS2
64
26884.875
0.114611500

SPP1
62
9094.989
0.165916760

APOB
54
28650.291
0.087770930

VCAN
52
9312.416
0.134104030

ICAM1
52
10434.576
0.122722970

CXCL1
52
10233.371
0.106985554

APOE
51
17188.280
0.110692450

STAT1
51
19835.434
0.081156254

KRAS
51
27196.115
0.077731330

BGN
48
11646.911
0.141241250

C3
42
9121.711
0.088115714

3.6. Survival Analysis of Hub Genes

Survival analysis of seven candidate genes was further studied using the KM plotter. The results showed that FN1, TIMP1, SPP1, APOE, and VCAN were related to OS in gastric cancer patients () (Figure 6).

(a)

(b)

(c)

(d)

(e)

Figure 6: Survival analysis of key genes by the KM plotter in gastric cancer. (a) Gastric cancer patients with high expression of FN1 had poor prognosis. (b) Patients of gastric cancer with high expression of TIMP1 had poor prognosis. (c) Patients of gastric cancer with high expression of SPP1 had poor prognosis. (d) Patients of gastric cancer with high expression of APOE had poor prognosis. (e) Patients of gastric cancer with high expression of VCAN had poor prognosis ().

4. Discussion

Microarray technology is a product of the gradual implementation of the human genome project and the rapid development and application of molecular biology. With the rapid development of gene microarray technology, people can quickly measure the expression levels of thousands of genes simultaneously [32, 33]. DNA microarray has a high-throughput speed, the characteristics of high sensitivity, and so on. It is widely used to study gene expression in a variety of organisms, including yeast, plants, and humans [34]. DNA microarrays can effectively identify disease-related genes and use them as biomarkers for diagnosis and treatment [35], including hepatocellular cancer [36], renal cell cancer [37], and colorectal cancer [38]. It has been reported that the methylation status of CpG island can be detected by microarray method, which is helpful for cancer diagnosis and detection of recurrence [39]. And studies have shown that through the microarray analysis between test specimens of gastric cancer and adjacent nontumor specimens, circRNA expression changes reveal the circRNA potential role in gastric cancer [40]. The above studies suggest that we can provide useful information for elucidating the development mechanism of gastric cancer and searching for new therapeutic targets and biomarkers through microarray technology.

In this study, we screened a total of 1034 differentially expressed genes from GSE54129 and GSE13911 gene expression profiles, among which 403 genes were upregulated, and 631 genes were downregulated. The KEGG pathway enrichment analysis revealed that the DEGs were mainly in the pathways in cancer and PI3K-Akt signaling pathway. Previous studies had shown that gastric cancer cell proliferation can be promoted by activating PI3K/Akt signaling pathway [41].

To explore the pathogenesis of gastric cancer, we constructed PPI network for systematic analysis. In the STRING database, we set the minimum interaction score with parameters >0.400 to obtain the PPI network of protein interactions. This setup avoided noise and incomplete data affecting the PPI network. MCODE discovers dense regions in PPI networks based on connection data. This function is not affected by the false-positive effect of high-throughput technology. We selected the parameters , , , and to find out all the modules that interact in the network, so as to further explore the most powerful proteins in the network. Through the module analysis of the network, we found that there were 27 modules in the network. The first-rank module with a score of 15.488 was the most closely related module in the whole network. Previous studies had shown that the related genes can be screened out more accurately by modular analysis, including cervical cancer [42], glioblastoma multiforme [43], and squamous cell cancer of head and neck [44]. These indicated that modular analysis played an important role in screening molecular markers. The genes in the module with higher score were the key genes that affected the occurrence of disease.

In order to further analyze the whole PPI network, the centrality analysis was used to explain the importance of the nodes in the network and the influence of the nodes on the network. We obtained 18 genes with high central values. Moreover, seven of the 18 genes were located in first-order module. Previous studies had shown that central nodes connected more protein-protein interactions, and central nodes also had more information for path enrichment analysis, which had a notable effect in the whole network [45, 46]. The results suggested that these genes may play significant roles in gastric cancer.

The study of cancer survival analysis plays an important role in the evaluation of cancer prevention measures [47]. Our results showed that the abnormal expression of FN1, TIMP1, SPP1, APOE, and VCAN influenced the prognosis of patients with gastric cancer. The high expression of them was associated with the worse OS in gastric cancer. FN1 is a member of the glycoprotein family and is widely expressed in many cell types. It played an important role in cell adhesion, growth, migration, and differentiation [48]. Studies had shown that inhibition of the expression of FN1 can inhibit the invasion and migration of gastric cancer cells [49, 50]. TIMP1 belongs to the TIMP gene family. This gene family encodes proteins that are natural inhibitors of matrix metalloproteinases (MMPS) [51]. TIMP1 can promote the aggregation of tumor-associated fibroblasts in the body and promote the proliferation and migration of cancer cells, and has antiapoptotic function [52, 53]. Studies had shown that TIMP1 was overexpressed and promoted cell proliferation in patients with gastric cancer through the NF-κB-dependent mechanism [54]. SPP1, also known as osteopontin (OPN), is an acidic glycoprotein that secretes several functions. OPN participated in the epithelial mesenchymal metastasis (EMT) pathway and played an important role in tumor metastasis [55]. Previous studies had shown that the high expression of OPN was closely related to the occurrence of gastric cancer [56]. APOE acts as the member of the family of low density lipoprotein (LDL) receptor ligands and interacts with the cell membrane receptor, which involves in cholesterol and other lipid transport [57]. APOE activated the PI3K/AKT/mTOR signaling pathway and played an important regulatory role in angiogenesis, tumor cell growth, and metastasis [58, 59]. A study had shown that tumor-associated macrophages (TAMs) interacted with gastric cancer cells through APOE to promote the migration of gastric cancer cells [60]. However, there are few reports of VCAN in gastric cancer. VCAN is a chondroitin sulfate proteoglycan. The interaction between extracellular matrix and cell surface proteins promoted cell growth, proliferation, and differentiation [61]. A study had shown that upregulation of VCAN promoted the migration and invasion of ovarian cancer cells by activating the NF-κB signaling pathway [62]. The growth of renal cell cancer cells can be inhibited by the activation of TNF signaling pathway through the silencing gene VCAN [63]. In addition, VCAN stimulated mesothelioma growth by weakening the antitumor activity of macrophages [64]. The above conclusions indicated that in GC, VCAN may affect the disease progression of patients through these paths, which were worth further study. The above results indicated that these five genes were closely related to the prognosis of patients with gastric cancer and can be used as a biomarker for GC.

5. Conclusion

In this study, 1034 differentially expressed genes were identified. On based of these genes, GO and KEGG results showed they were mainly in plasma membrane, cell adhesion, and PI3K-Akt signaling pathway. Moreover, 18 topological key genes of the 1st-rank module were focused on. Furthermore, five of them (FN1, TIMP1, SPP1, APOE, and VCAN) were found to be related to gastric cancer. Therefore, it provides new research directions for the detection and treatment of gastric cancer. However, their involvement in the molecular mechanisms of disease needs further clinical studies.

Data Availability

The GSE54129 and GSE13911 data used to support the findings of this study are included within the article.

Disclosure

Xinyu Chong and Rui Peng are co-first authors.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (81570747).

References

P. Rawla and A. Barsouk, “Epidemiology of gastric cancer: global trends, risk factors and prevention,” Gastroenterology Review, vol. 14, no. 1, pp. 26–38, 2019.

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

P. Karimi, F. Islami, S. Anandasabapathy, N. D. Freedman, and F. Kamangar, “Gastric cancer: descriptive epidemiology, risk factors, screening, and prevention,” Cancer Epidemiology, Biomarkers & Prevention, vol. 23, no. 5, pp. 700–713, 2014.

A. R. Yusefi, K. Bagheri Lankarani, P. Bastani, M. Radinmanesh, and Z. Kavosi, “Risk factors for gastric cancer: a systematic review,” Asian Pacific Journal of Cancer Prevention, vol. 19, no. 3, pp. 591–603, 2018.

J. Poorolajal, L. Moradi, Y. Mohammadi, Z. Cheraghi, and F. Gohari-Ensaf, “Risk factors for stomach cancer: a systematic review and meta-analysis,” Epidemiol Health, vol. 42, article e2020004, 2020.

G. Herrera-Almario and V. E. Strong, “Minimally invasive gastric surgery,” Annals of Surgical Oncology, vol. 23, no. 12, pp. 3792–3797, 2016.

J. Tey, Y. Y. Soon, W. Y. Koh et al., “Palliative radiotherapy for gastric cancer: a systematic review and meta-analysis,” Oncotarget, vol. 8, no. 15, pp. 25797–25805, 2017.

O. Zhandossov, G. Kaussova, and A. Koten, “Combined treatment for gastric cancer: immunological approach,” The Turkish Journal of Gastroenterology, vol. 29, no. 2, pp. 151–156, 2018.

E. Niccolai, A. Taddei, D. Prisco, and A. Amedei, “Gastric cancer and the epoch of immunotherapy approaches,” World Journal of Gastroenterology, vol. 21, no. 19, pp. 5778–5793, 2015.

Z. Tan, “Recent advances in the surgical treatment of advanced gastric cancer: a review,” Medical Science Monitor, vol. 25, pp. 3537–3541, 2019.

S. P. Pfeifer, “From next-generation resequencing reads to a high-quality variant data set,” Heredity (Edinb), vol. 118, no. 2, pp. 111–124, 2017.

K. Owzar, W. T. Barry, S. H. Jung, I. Sohn, and S. L. George, “Statistical challenges in preprocessing in microarray experiments in cancer,” Clinical Cancer Research, vol. 14, no. 19, pp. 5959–5966, 2008.

J. Man, X. Zhang, H. Dong et al., “Screening and identification of key biomarkers in lung squamous cell carcinoma by bioinformatics analysis,” Oncology Letters, vol. 18, no. 5, pp. 5185–5196, 2019.

G. D’Angelo, “Microarray analysis in gastric cancer: a review,” World Journal of Gastroenterology, vol. 20, no. 34, p. 11972, 2014.

E. Clough and T. Barrett, “The gene expression omnibus database,” Methods in Molecular Biology, vol. 1418, pp. 93–110, 2016.

H. B. Hu, H. Y. Jie, and X. X. Zheng, “Three circulating LncRNA predict early progress of esophageal squamous cell carcinoma,” Cellular Physiology and Biochemistry, vol. 40, no. 1-2, pp. 117–125, 2016.

L. Li, Q. Lei, S. Zhang, L. Kong, and B. Qin, “Screening and identification of key biomarkers in hepatocellular carcinoma: evidence from bioinformatic analysis,” Oncology Reports, vol. 38, no. 5, pp. 2607–2618, 2017.

D. W. Huang, B. T. Sherman, Q. Tan et al., “DAVID bioinformatics resources: expanded annotation database and novel algorithms to better extract biology from large gene lists,” Nucleic Acids Research, vol. 35, suppl_2, pp. W169–W175, 2007.

P. D. Thomas, “The gene ontology and the meaning of biological function,” Methods in Molecular Biology, vol. 1446, pp. 15–24, 2017.

M. Kanehisa and S. Goto, “KEGG: Kyoto Encyclopedia of Genes and Genomes,” Nucleic Acids Research, vol. 28, no. 1, pp. 27–30, 2000.

W. D. Jiang and Z. H. Ye, “Integrated analysis of a competing endogenous RNA network in renal cell carcinoma using bioinformatics tools,” Bioscience Reports, vol. 39, no. 7, 2019.

D. Szklarczyk, A. L. Gable, D. Lyon et al., “STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets,” Nucleic Acids Research, vol. 47, no. D1, pp. D607–d613, 2019.

L. Wang, Z. Li, X. Song, L. Liu, G. Su, and Y. Cui, “Bioinformatic analysis of genes and microRNAs associated with atrioventricular septal defect in down syndrome patients,” International Heart Journal, vol. 57, no. 4, pp. 490–495, 2016.

M. Kohl, S. Wiese, and B. Warscheid, “Cytoscape: software for visualization and analysis of biological networks,” Methods in Molecular Biology, vol. 696, pp. 291–303, 2011.

G. D. Bader and C. W. Hogue, “An automated method for finding molecular complexes in large protein interaction networks,” BMC Bioinformatics, vol. 4, no. 1, p. 2, 2003.

Y. Tang, M. Li, J. Wang, Y. Pan, and F. X. Wu, “CytoNCA: a cytoscape plugin for centrality analysis and evaluation of protein interaction networks,” Biosystems, vol. 127, pp. 67–72, 2015.

X. Luo, L. Guo, X. J. Dai et al., “Abnormal intrinsic functional hubs in alcohol dependence: evidence from a voxelwise degree centrality analysis,” Neuropsychiatric Disease and Treatment, vol. Volume 13, pp. 2011–2020, 2017.

Y. Li, W. Li, Y. Tan, F. Liu, Y. Cao, and K. Y. Lee, “Hierarchical decomposition for betweenness centrality measure of complex networks,” Scientific Reports, vol. 7, no. 1, 2017.

C. F. A. Negre, U. N. Morzan, H. P. Hendrickson et al., “Eigenvector centrality for characterization of protein allosteric pathways,” Proceedings of the National Academy of Sciences of the United States of America, vol. 115, no. 52, pp. E12201–e12208, 2018.

J. H. Pan, H. Zhou, L. Cooper et al., “LAYN is a prognostic biomarker and correlated with immune infiltrates in gastric and colon cancers,” Frontiers in Immunology, vol. 10, p. 6, 2019.

G. Liu, J. Zhao, B. Pan, G. Ma, and L. Liu, “UBE2C overexpression in melanoma and its essential role in G2/M transition,” Journal of Cancer, vol. 10, no. 10, pp. 2176–2184, 2019.

A. Brazma and J. Vilo, “Gene expression data analysis,” FEBS Letters, vol. 480, no. 1, pp. 17–24, 2000.

R. Bumgarner, “Overview of DNA microarrays: types, applications, and their future,” Current Protocols in Molecular Biology, vol. 22, 2013.

M. Kurella, L. L. Hsiao, T. Yoshida et al., “DNA microarray analysis of complex biologic processes,” J Am Soc Nephrol, vol. 12, no. 5, pp. 1072–1078, 2001.

S. M. Yoo, J. H. Choi, S. Y. Lee, and N. C. Yoo, “Applications of DNA microarray in disease diagnostics,” Journal of Microbiology and Biotechnology, vol. 19, no. 7, pp. 635–646, 2009.

T. L. He, K. L. Zheng, G. Li, B. Song, and Y. J. Zhang, “Identification of typical miRNAs and target genes in hepatocellular carcinoma by DNA microarray technique,” European Review for Medical and Pharmacological Sciences, vol. 18, no. 1, pp. 108–116, 2014.

J.-Y. Feng, X.-W. Diao, M.-Q. Fan et al., “Screening of feature genes of the renal cell carcinoma with DNA microarray,” European Review for Medical and Pharmacological Sciences, vol. 17, no. 22, pp. 2994–3001, 2013.

P. Kwiatkowski, P. Wierzbicki, A. Kmieć, and J. Godlewski, “DNA microarray-based gene expression profiling in diagnosis, assessing prognosis and predicting response to therapy in colorectal cancer,” Postȩpy Higieny i Medycyny Doświadczalnej (Online), vol. 66, pp. 330–338, 2012.

P. Hou, J. Y. Shen, M. J. Ji, N. Y. He, and Z. H. Lu, “Microarray-based method for detecting methylation changes of p16 (Ink4a) gene 5'-CpG islands in gastric carcinomas,” World Journal of Gastroenterology, vol. 10, no. 24, pp. 3553–3558, 2004.

W. Gu, Y. Sun, X. Zheng et al., “Identification of gastric cancer-related circular RNA through microarray analysis and bioinformatics analysis,” BioMed Research International, vol. 2018, Article ID 2381680, 9 pages, 2018.

C. Zhao, T. Tao, L. Yang et al., “Loss of PDZK1 expression activates PI3K/AKT signaling via PTEN phosphorylation in gastric cancer,” Cancer Letters, vol. 453, pp. 107–121, 2019.

L. Xia, X. Su, J. Shen et al., “ANLN functions as a key candidate gene in cervical cancer as determined by integrated bioinformatic analysis,” Cancer Management and Research, vol. 10, pp. 663–670, 2018.

Q. Yang, R. Wang, B. Wei et al., “Candidate biomarkers and molecular mechanism investigation for glioblastoma multiforme utilizing WGCNA,” BioMed Research International, vol. 2018, Article ID 4246703, 10 pages, 2018.

B. Yang, Z. Chen, Y. Huang, G. Han, and W. Li, “Identification of potential biomarkers and analysis of prognostic values in head and neck squamous cell carcinoma by bioinformatics analysis,” Oncotargets and Therapy, vol. Volume 10, pp. 2315–2321, 2017.

D. A.-Y. O. N. G. ZHUANG, L. I. JIANG, Q. I. N. G.-Q. I. N. G. HE, P. E. N. G. ZHOU, and T. A. O. YUE, “Identification of hub subnetwork based on topological features of genes in breast cancer,” International Journal of Molecular Medicine, vol. 35, no. 3, pp. 664–674, 2015.

M. Abedi and Y. Gheisari, “Nodes with high centrality in protein interaction networks are responsible for driving signaling pathways in diabetic nephropathy,” Peer J, vol. 3, p. e1284, 2015.

H. Imai, K. Kaira, and K. Minato, “Clinical significance of post-progression survival in lung cancer,” Thorac Cancer, vol. 8, no. 5, pp. 379–386, 2017.

B. Li, W. Shen, H. Peng et al., “Fibronectin 1 promotes melanoma proliferation and metastasis by inhibiting apoptosis and regulating EMT,” Oncotargets and Therapy, vol. Volume 12, pp. 3207–3221, 2019.

T. P. Xu, M. D. Huang, R. Xia et al., “Decreased expression of the long non-coding RNA FENDRR is associated with poor prognosis in gastric cancer and FENDRR regulates gastric cancer cell metastasis by affecting fibronectin 1 expression,” Journal of Hematology & Oncology, vol. 7, no. 1, p. 63, 2014.

H. B. Chen and H. T. Zheng, “MicroRNA-200c represses migration and invasion of gastric cancer SGC-7901 cells by inhibiting expression of fibronectin 1,” European Review for Medical and Pharmacological Sciences, vol. 21, no. 8, pp. 1753–1758, 2017.

G. Cheng, X. Fan, M. Hao, J. Wang, X. Zhou, and X. Sun, “Higher levels of TIMP-1 expression are associated with a poor prognosis in triple-negative breast cancer,” Molecular Cancer, vol. 15, no. 1, p. 30, 2016.

J. Zhang, T. Wu, S. Zhan et al., “TIMP-1 and CD82, a promising combined evaluation marker for PDAC,” Oncotarget, vol. 8, no. 4, pp. 6496–6512, 2017.

Y. Gong, E. Scott, R. Lu, Y. Xu, W. K. Oh, and Q. Yu, “TIMP-1 promotes accumulation of cancer associated fibroblasts and cancer progression,” PLoS One, vol. 8, no. 10, article e77366, 2013.

O. M. Omar, M. Soutto, N. S. Bhat et al., “TFF1 antagonizes TIMP-1 mediated proliferative functions in gastric cancer,” Molecular Carcinogenesis, vol. 57, no. 11, pp. 1577–1587, 2018.

R. Jia, Y. Liang, R. Chen et al., “Osteopontin facilitates tumor metastasis by regulating epithelial-mesenchymal plasticity,” Cell Death Dis, vol. 7, no. 12, p. e2564, 2016.

T. Chen, L. Sun, C. He, Y. Gong, Q. Xu, and Y. Yuan, “Serum OPN expression for identification of gastric cancer and atrophic gastritis and its influencing factors,” PLoS One, vol. 9, no. 12, article e114005, 2014.

C. E. Yu, E. Cudaback, J. Foraker et al., “Epigenetic signature and enhancer activity of the human APOE gene,” Human Molecular Genetics, vol. 22, no. 24, pp. 5036–5047, 2013.

Z. Zhao, S. Zou, X. Guan et al., “Apolipoprotein E overexpression is associated with tumor progression and poor survival in colorectal cancer,” Frontiers in Genetics, vol. 9, p. 650, 2018.

T. Vogel, N. H. Guo, R. Guy et al., “Apolipoprotein E: a potent inhibitor of endothelial and tumor cell proliferation,” Journal of Cellular Biochemistry, vol. 54, no. 3, pp. 299–308, 1994.

P. Zheng, Q. Luo, W. Wang et al., “Tumor-associated macrophages-derived exosomes promote the migration of gastric cancer cells by transfer of functional Apolipoprotein E,” Cell Death & Disease, vol. 9, no. 4, p. 434, 2018.

K. Asano, C. M. Nelson, S. Nandadasa et al., “Stromal versican regulates tumor growth by promoting angiogenesis,” Scientific Reports, vol. 7, no. 1, p. 17225, 2017.

T. L. Yeung, C. S. Leung, K. K. Wong et al., “TGF-β modulates ovarian cancer invasion by upregulating CAF-derived versican in the tumor microenvironment,” Cancer Research, vol. 73, no. 16, pp. 5016–5028, 2013.

Y. Mitsui, H. Shiina, T. Kato et al., “Versican promotes tumor progression, metastasis and predicts poor prognosis in renal carcinoma,” Molecular Cancer Research, vol. 15, no. 7, pp. 884–895, 2017.

A. G. Pappas, S. Magkouta, I. S. Pateras et al., “Versican modulates tumor-associated macrophage properties to stimulate mesothelioma growth,” Oncoimmunology, vol. 8, no. 2, article e1537427, 2018.