Skip to main content

Identification of hub genes and therapeutic drugs in esophageal squamous cell carcinoma based on integrated bioinformatics strategy



Esophageal squamous cell carcinoma (ESCC) is one of leading malignant cancers of gastrointestinal tract worldwide. Until now, the involved mechanisms during the development of ESCC are largely unknown. This study aims to explore the driven-genes and biological pathways in ESCC.


mRNA expression datasets of GSE29001, GSE20347, GSE100942, and GSE38129, containing 63 pairs of ESCC and non-tumor tissues data, were integrated and deeply analyzed. The bioinformatics approaches include identification of differentially expressed genes (DEGs) and hub genes, gene ontology (GO) terms analysis and biological pathway enrichment analysis, construction and analysis of protein–protein interaction (PPI) network, and miRNA–gene network construction. Subsequently, GEPIA2 database and qPCR assay were utilized to validate the expression of hub genes. DGIdb database was performed to search the candidate drugs for ESCC.


Finally, 120 upregulated and 26 downregulated DEGs were identified. The functional enrichment of DEGs in ESCC were mainly correlated with cell cycle, DNA replication, deleted in colorectal cancer (DCC) mediated attractive signaling pathway, and Netrin-1 signaling pathway. The PPI network was constructed using STRING software with 146 nodes and 2392 edges. The most significant three modules in PPI were filtered and analyzed. Totally ten genes were selected and considered as the hub genes and nuclear division cycle 80 (NDC80) was closely related to the survival of ESCC patients. DGIdb database predicted 33 small molecules as the possible drugs for treating ESCC.


In summary, the data may provide new insights into ESCC pathogenesis and treatments. The candidate drugs may improve the efficiency of personalized therapy in future.


Esophageal cancer (EC) ranks seventh in terms of incidence and sixth in cancer deaths worldwide, responsible for about 572,000 new cases and 509,000 deaths last year [1]. Although we have made great progress on the early diagnosis and novel therapy, EC still is one of challengeable diseases in Eastern Asian [1]. Generally, EC includes two most common histologic subtypes: esophageal squamous cell carcinoma (ESCC) and esophageal adenocarcinoma (EAC) [2]. ESCC comprises over 90% of all EC cases [1]. And risk factors, such as smoking and hot drinks, are closely related to the initiation of ESCC [1, 2]. However, the underlying mechanisms of ESCC are not well understood. And due to the lack of specific biomarkers, most ESCC patients are diagnosed at a late stage, leading to particularly poor outcomes of patients [3]. Even worse, some of ESCC patients suffer from tumor recurrence due to the chemotherapy resistance [3]. Therefore, it is of paramount importance to find novel biomarkers and effective targets for ESCC patients.

Recently, gene profile and gene chip have been extensively applied in the field of scientific researches [4, 5]. Gene expression analysis based on these methods can quickly detect the differentially expressed genes (DEGs) that may have a strong influence on cancer progression [6]. However, most of the gene chip or gene profile data have been only deposited in public databases. And re-analyzing these data can be an efficient way to provide the new insights into further studies. So far, many studies have used gene chip or gene profile to identify key genes for ESCC, and numerous DEGs have been detected [7]. Nevertheless, the results may be inconsistent and variable because of the existence of tumor heterogeneity. To date, few reliable biomarkers and therapeutic targets have been identified for ESCC [8]. Thus, it’s urgent to discover new markers and therapeutic targets for ESCC patients.

Many chemotherapeutic drugs have shown activity against ESCC, including docetaxel [9,10,11], cisplatin [10, 11], fluorouracil [9,10,11], and nedaplatin [9]. Moreover, the combinations of these agents are also recommended because of the existence of chemotherapy resistance. A recent study found that concurrent chemoradiotherapy (CCRT) with 5-fluorouracil plus cisplatin were more effective and less toxic than CCRT with the docetaxel plus cisplatin as the first-line treatment for ESCC patients [11]. However, the progression-free survival and overall survival (OS) of ESCC patients remained short, highlighting the importance of developing some molecular drugs.

In the study, four mRNA expression profiles were downloaded (GSE29001 [12], GSE20347 [13], GSE100942 [14], and GSE38129 [15]) from GEO database, from which there are 63 pairs of ESCC and non-tumor tissues data available. Integrated analyses included identifying DEGs using the GEO2R tool, overlapping four datasets using a Venn diagram tool, GO terms analysis, biological pathway enrichment analysis, PPI construction, hub genes identification and verification, miRNA–hub genes network construction, and exploration of the candidate small molecular drugs for ESCC.

Materials and methods

Data collection

ESCC and adjacent normal tissue gene expression profiles of GSE20347 [13], GSE29001 [12], GSE100942 [14], and GSE38129 [15] were downloaded from GEO ( database [16]. The microarray data of GSE29001 was based on GPL571 Platforms (Affymetrix Human Genome U133A 2.0 Array) and included 12 pairs of ESCC and non-tumor tissues (Submission date: May 02, 2011). The GSE20347 data was based on GPL571 Platforms (Affymetrix Human Genome U133A 2.0 Array) and included 17 ESCC tissues and 17 normal tissues (Submission date: Feb 16, 2010). The GSE100942 data was based on GPL570 Platforms (Affymetrix Human Genome U133 Plus 2.0 Array) and included 4 ESCC tissues and 4 non-tumor tissues (Submission date: Jul 07, 2017). The GSE38129 data was based on GPL571 Platforms (Affymetrix Human Genome U133A 2.0 Array) and included 30 pairs of ESCC and non-tumor tissues (Submission date: May 22, 2012). The above datasets met the following criteria: (1) they used tissue samples from human ESCC tissues and paired adjacent or non-tumor tissues; (2) each dataset involved more than eight samples.

DEGs identification

GEO2R ( was used to pick out the DEGs in ESCC tissues and adjacent non-tumor tissues [17]. p < 0.05 and |logFC| > 1 were set as the cut-off criterion to select DEGs for every dataset microarray respectively [7, 17]. Finally, the overlapping DEGs among the four datasets was identified by Venn diagram tool (

Cell culture, RNA extraction and quantitative PCR (qPCR)

Human ESCC cell line EC109 and human esophageal squamous epithelial cell line Het-1A were cultured in RPMI-1640 medium (Gibco) with 10% fetal bovine serum (Gibco) at 37 °C in a humidified atmosphere with 5% CO2. Total RNA was extracted from cells using the E.Z.N.A.™ Total RNA Kit I (OMEGA). PrimeScript™ RT Master Mix (Perfect Real Time) was used for RNA reverse transcription. SYBR Premix Ex Taq (TaKaRa) was employed to conduct qPCR assay. PCR primers were designed and synthesized by TaKaRa (Additional file 1: Table S1). The experiments were performed in three times. GAPDH was used as the internal control.

GO and biological pathway enrichment analysis

GO terms analysis of selected DEGs were performed using the DAVID database (; version: 6.8) [18]. We submitted the DEGs, including 120 upregulated genes and 26 downregulated genes, into DAVID with p < 0.05 as the cut-off criterion. The GO results of significant terms for cellular component (CC), biological process (BP), and molecular function (MF) were ranked by p-value and exhibited as bar charts. The FunRich tool (version: 3.0) was mainly used for analyzing the functional enrichment and interaction networks of genes and proteins [19]. In this study, the FunRich was used to analyze the biological pathways of DEGs. Finally, the top 10 biological pathways of upregulated genes and downregulated genes were presented as bar charts, respectively. p-value < 0.05 was considered as statistically significant.

PPI network construction and analysis

PPI networks are the networks of protein complexes formed as the results of biochemical or electrostatic forces [20]. PPI network is crucial for molecular processes, and abnormal PPI is the basis of many diseases, including tumors [21]. In this study, the Search Tool for the Retrieval of Interacting Genes (STRING) database (; version: 11.0) [20], Cytoscape software (version: 3.6.1) [22], and FunRich were utilized to construct PPI networks. Cytoscape and FunRich tool were applied to present the PPI networks with the cut-off criterion as confidence score ≥ 0.4, maximum number of interactors = 0. The Molecular Complex Detection (MCODE) plug-in of Cytoscape tool was employed to visualize the significant gene modules in ESCC with degree cutoff = 2, node score cutoff = 0.2, k-core = 2, and max. depth = 100. The criteria for selecting the top 3 significant modules were set as follows: MCODE scores ≥ 4 and number of nodes ≥ 4 [23]. FunRich tool was performed to do the functional enrichment for each module. 10 hub genes with high degree of connectivity were selected and mapped into PPI based on STRING following confidence score ≥ 0.4, maximum number of interactors ≤ 5. Furthermore, STRING was used to perform the co-expression analysis of hub genes.

Validation of the hub genes

The GEPIA2 ( is an online database for analyzing gene expression profiles of 9736 tumors and 8587 normal samples from the Cancer Genome Atlas (TCGA) and the genotype-tissue expression (GTEx) projects [24]. Thus, we can validate the expression levels and genes correlations of hub genes in ESCC tissues and normal tissues. The cBio Cancer Genomics Portal (; version: 2.2.0) is an open access tool which provides analysis, visualization, and downloads of cancer genomics datasets of many types of tumors [25]. Complex cancer genomics profiles are accessible from the cBioPortal tool, thus enabling us to compare the genetic alterations of the selected ten hub genes in ESCC.

miRNA–hub gene network

The targeted miRNAs of hub genes were predicted by four established miRNA target prediction databases [miRanda, PITA, PicTar, and TargetScan (version: 3.1)]. The miRNAs predicted by at least two programs were selected as the targeted miRNAs of hub genes. A co-expression network based on correlation analysis of hub genes and miRNAs associated with cancer was constructed by Cytoscape software. In the network, a green circular node represented the miRNA and a red circular node represented the hub gene, their interaction was represented by an arrow. The numbers of arrows in the networks indicated the contribution of one miRNA to the surrounding genes, and the higher the degree, the more central the hub gene was within the network.

Drug-hub gene interaction

The 10 hub genes were also served as the promising targets for searching drugs through the DGIdb (; version: 3.0.2—sha1 ec916b2) [26]. This database contains drug-gene interaction data from 30 disparate sources including ChEMBL, DrugBank, Ensembl, NCBI Entrez, PharmGKB, PubChem, clinical trial databases, and literature in NCBI PubMed. The results of this process were arranged so that each entry was a specific drug-gene interaction associated with its source link [27]. Drugs supported by more than one databases or PubMed references were selected as the potential drugs. The final list only involved the drugs that have been approved by the Food and Drug Administration (FDA). The identified target network was visualized using STITCH (; version: 5.0), a program similar to STRING which also incorporated drug-gene relationships [27, 28].


Identification of DEGs

The four mRNA expression profiles (GSE29001 [12], GSE20347 [13], GSE100942 [14], and GSE38129 [15]), including 63 pairs of ESCC tissues and adjacent normal tissues, were included in this study. Using p < 0.05 and |logFC| > 1 as cut-off criterion [7, 17], we extracted 2594, 1295, 833, and 520 DEGs from the expression profile datasets GSE29001, GSE20347, GSE100942, and GSE38129, respectively. By using Venn diagrams to overlap the DEGs of the four profile datasets, a total of 146 overlapping DEGs were identified (Fig. 1; Table 1), including 120 upregulated genes and 26 downregulated genes. Employing FunRich, we constructed a heatmap of the 120 upregulated and 26 downregulated DEGs using data profile GSE20347 as a reference. Additional file 1: Figure S1 showed the differential distribution of the 146 DEGs.

Fig. 1

Identification of 164 DEGs from the four cohort profile datasets (GSE29001, GSE20347, GSE100942, and GSE38129). DEGs were screened out by GEO2R ( tool, statistically significant DEGs were defined with p < 0.05 and |logFC| > 1 as the cut-off criterion. Venn diagram tool ( was used to identify the overlapping DEGs in the four datasets. Different color areas represented different datasets. The cross areas meant the overlapping DEGs

Table 1 146 DEGs were identified from the four cohort profile datasets, including 120 upregulated genes and 26 downregulated genes in the ESCC tissues compared to normal controls

Enrichment analysis and PPI network

Functional enrichment analysis of DEGs

GO and pathways enrichment were analyzed through multiple databases or software, including DAVID [18], KEGG pathway (; release 89.0) [29], and FunRich software [19] with p < 0.05 as the cut-off criterion.

GO analysis of DEGs classified DEGs into three functional groups: CC, BP, and MF group (Fig. 2). As shown in Fig. 2a, in the CC group, the upregulated genes were enriched in kinetochore, nucleus, microtubule, nucleoplasm, and condensed chromosome kinetochore, whereas the downregulated genes were related to collagen type XIV, mast cell granule, collagen, Z disc, and extracellular matrix (Fig. 2b). In the MF group, the upregulated genes were mainly enriched in chromatin binding, motor activity, protein binding, protein serine/threonine kinase activity, and DNA binding, whereas the downregulated genes were involved in heat shock protein activity, transcription factor binding, and peroxidase activity. As for BP, the upregulated genes were correlated to cell growth and/or maintenance, cell cycle, chromosome segregation, cell communication, and signal transduction. The downregulated genes were significantly connected with muscle contraction. GO term analysis showed that most of the DEGs were enriched in kinetochore, collagen, binding functions, cell cycle, and cell growth. The data were in keeping with the knowledge that abnormality of cell cycle and cell growth regulators was the major cause of tumorigenesis [30]. Moreover, the metabolism of the nuclear components and intercellular substances of cancer cells was different from that in normal cells [31, 32].

Fig. 2

GO analysis and significant enriched GO terms of DEGs in ESCC a, upregulated DEGs; b downregulated DEGs). GO analysis classified the DEGs into 3 groups (cellular component, molecular function, and biological process)

As for biological pathway enrichment, the upregulated genes were enriched in cell cycle, DNA replication, mitotic M–M/G1 phases, M phase, and Polo-like kinase 1 (PLK1) signaling pathway (Fig. 3a). Previous investigations have demonstrated that some cell cycle-related genes in ESCC development could predict the OS of ESCC patients [33]. Recent evidence also indicated that PLK1 signaling pathway could regulate cell cycle [34]. Moreover, the genes involved in DNA replication system might be useful markers to predict tumor progression [35]. The downregulated genes were enriched in DCC mediated pathway, Netrin-1 signaling pathway, Flavin-containing monooxygenases (FMO) oxidizes nucleophiles, noradrenaline and adrenaline degradation, and ethanol degradation II (cytosol) (Fig. 3b). Previously, we have reported that substance P (SP)/NK-1R signaling could promote the growth and metaseries of ESCC, suggesting that interactions between cancers and nervous system were indispensable for understanding the biological mechanisms of tumorigenesis [36]. Here, we found that the downregulated DEGs were involved in two nervous system-related pathways (DCC mediated signaling pathway and Netrin-1 signaling pathway). Progress toward exploring the links between neuronal activity and oncology might provide new insight into cancer biology. Moreover, specific metabolic activities can be directly involved in cancer progression [35]. The downregulated DEGs were related to metabolic activities such as noradrenaline and adrenaline degradation, suggesting that these metabolism-targeted pathways may be valuable for improving the treatment efficiency.

Fig. 3

Significantly enriched biological pathway terms of DEGs in ESCC. a Biological pathway for upregulated DEGs. b Biological pathway for downregulated DEGs. DEGs functional and signaling pathway enrichment were conducted using KEGG pathway ( and FunRich tool

PPI network construction and modules analysis

Using STRING database [20] and Cytoscape software [22], totally 146 DEGs were mapped into the PPI network, including 146 nodes and 2392 edges (Additional file 1: Figure S2). The PPI enrichment p-value was 1.0 × 10−16. Additional file 1: Figure S3 shows the interaction network of the 146 DEGs and their related genes, allowing us to evaluate their biological functions. For example, AURKA and TPX2 belong to the “Role of Ran in mitotic spindle regulation” pathway, and TPX2 knockdown could inhibit the cell proliferation of ESCC cells [37]. The top three significant clusters within PPI network were selected using MCODE plug-in in Cytoscape software (Module 1, MCODE score = 51.778; Module 2, MCODE score = 7; Module 3, MCODE score = 4). We also analyzed the functions of each module (Fig. 4). Pathway enrichment analysis indicated that Module 1 consisted of 55 nodes and 1398 edges (Fig. 4a, b), which were mainly associated with cell cycle, DNA replication, and PLK1 signaling pathway. Module 2 consisted of 7 nodes and 21 edges (Fig. 4c, d), which were mainly associated with platelet adhesion to exposed collagen, epithelial-to-mesenchymal transition, and VEGFR3 signaling in lymphatic endothelium. Module 3 consisted of 4 nodes and 6 edges (Fig. 4e, f), which were associated with regulation of Insulin-like growth factor (IGF) activity by Insulin-like growth factor binding proteins (IGFBPs), EGF receptor (ErbB1) signaling pathway, and class I Phosphoinositide 3-kinase (PI3K) signaling events.

Fig. 4

The top 3 modules from the PPI network. a Module 1; b the top 5 enriched pathways of module 1; c module 2; d the top 5 enriched pathways of module 2; e module 3; f the top 5 enriched pathways of module 3

Using cytoHubba software, ten genes (Cyclin-dependent kinase 1 (CDK1), Cyclin B1 (CCNB1), DNA topoisomerase II alpha (TOP2A), Cyclin B2 (CCNB2), BUB1 mitotic checkpoint serine/threonine kinase (BUB1), Cyclin A2 (CCNA2), Non-SMC condensin I complex subunit G (NCAPG), Aurora kinase B (AURKB), NDC80, and BUB1B) with higher degree of connectivity were identified as hub genes (Table 2). Moreover, the PPI network of ten hub genes was established using Cytoscape software (Fig. 5a). The interaction network of ten hub genes and their related genes was also established by the FunRich (Fig. 5b). The related genes here were defined as these genes connected to the hub genes. The hub genes and the related genes could be enriched in biological pathways according to the enrichment functions of FunRich tool. The gene co-expression analysis of the ten hub genes performed by STRING database showed that these genes might be actively interacted with each other (Fig. 5c). The above findings suggested that these hub genes might play a crucial role in ESCC progression. For example, AURKB was found to increase in the early development stage of ESCC and might influence the initiation of ESCC [38]. Interestingly, two important genes (BUB1 and BUB1B), which affect the chemotherapy of ESCC [39], could interact with AURKB according to the Fig. 5a. The pathways consisted by these three genes might be important for improving the drug sensitivity of early ESCC.

Table 2 Top 10 hub genes with higher degree of connectivity
Fig. 5

PPI network construction and co-expression analysis for the 10 hub genes in ESCC. a A total of 10 hub genes with higher degree of connectivity were selected and filtered into the PPI network complex using the STRING online database. b The PPI network of the 10 hub genes and their related genes, created by the FunRich software. c The co-expression analysis of 10 hub genes using the STRING online database

Genetic information and hub genes expression

Kaplan–Meier-plotter website ( was used to analyze the prognostic information of the ten hub genes. The result showed that the upregulated NDC80 was closely related to the OS of patients with ESCC (Fig. 6a). The deregulation of NDC80 caused by amplification and mutation might lead to poor OS. The remaining 9 genes presented similar trends to that of NDC80, but not statistically significant (Additional file 1: Table S2). Then, we used cBioPortal to enquiry the genetic alterations of the hub genes. And Fig. 6b presented the network constructed by the 10 hub genes and their 40 most frequently altered neighbor genes. Besides, drugs targeting the 10 genes were illustrated. Figure 6b showed that only TOP2A, CDK1, CCNB1, and AURKB were identified as chemotherapy targets currently. We therefore supposed that the other 6 genes (CCNB2, BUB1, BUB1B, CCNA2, NCAPG, and NDC80) might be the novel targets in the future. Figure 6c, d presented the alteration information of the ten genes. The ten hub genes were changed in 24 (25%) of 96 sequenced patients (96 total). NDC80 and BUB1 were changed most often (6% and 5%), these include amplification, mutation and so on.

Fig. 6

a High expression of NDC80 was significantly associated with poor OS in ESCC patients, using a Kaplan–Meier curve and a log-rank test (Kaplan–Meier-plotter website; b The network contained 50 nodes, including 10 hub genes and the 40 most frequently altered neighbor genes. The relationships between 10 hub genes and drugs were also presented. c A visual summary across a set of ESCC (data from esophageal squamous cell carcinoma, TCGA, Provisional) showed the genetic alterations connected with the 10 hub genes which were altered in 24 (25%) of 96 sequenced cases/patients (96 total). d An overview of changes in the 10 hub genes in the genomics datasets of ESCC in TCGA database

We further conducted the expression analysis using the data from GEPIA2 database. The expression levels of the 10 hub genes were significantly different between ESCC and normal tissues (Fig. 7). The expression trends of the 10 genes from GEPIA2 database were in accordance to the data in GEO datasets. Real-Time PCR results revealed that the mRNA expressions of the 9 hub genes (except for CCNA2) were upregulated in EC109 cells (ESCC cell line) as compared to Het-1A cells (esophageal squamous epithelial cell line) (Additional file 1: Figure S5). In addition, the expression levels of the 10 hub genes in ESCC were positively correlated with each other using GEPIA2 database (data not shown).

Fig. 7

Validation of the hub genes in the Cancer Genome Atlas (TCGA) database. Box plots showed the mRNA expressions of the 10 hub genes using data from the TCGA database in GEPIA2 ( The validation results of 10 hub genes were in accordance with the profiles in our study, and their p-values < 0.05

miRNA–hub genes network

To investigate the regulatory relationships of identified hub genes and miRNAs, four miRNA targets prediction databases was used to predicted the targeted miRNAs of hub genes. The miRNAs predicted by at least two databases were selected as the targeted miRNAs of hub genes. The co-expression network based on the correlation analysis between the hub genes and miRNAs was constructed by Cytoscape software (Fig. 8). The numbers of miRNAs and mRNAs in the network were 175 and 10, respectively. In the network, the numbers of arrows in the networks indicated the contribution of one miRNA to the surrounding hub genes, and the higher the degree, the more central the hub gene was within the network. CDK1, TOP2A, and CCNA2 were identified as the three hub genes which were targeted by the most miRNAs. miR-543, miR-495-3p, and miR-590-3p were the top three miRNAs with the most target genes. Previously, Ma et al. demonstrated that miR-219-5p/CCNA2 axis could inhibit the cell proliferation and cell cycle distribution of ESCC cells, highlighting the role of CCNA2 in cell cycle and tumor growth [40]. In addition, miR-543 could facilitate cell mobility and invasion of ESCC by repressing PLA2G4A [41]. Therefore, the miRNA–hub genes network could provide powerful basis for understanding the molecular mechanisms of ESCC.

Fig. 8

miRNA–hub gene interaction network of ESCC. The green circular node represented the miRNA. The red circular node represented the hub gene. The arrow represented the interaction between the miRNAs and hub genes

Drug-gene interaction

Using the 10 hub genes to explore the drug-gene interactions, 33 drugs for possibly treating ESCC were compiled and selected (Table 3). Promising targets of these drugs include CDK1, TOP2A, CCNA2, and AURKB. Among these four genes, CDK1, TOP2A, and AURKB have been currently considered as the drug targets according to the cBioPortal database (Fig. 6b). The final list comprised only the drugs which were approved by FDA, and several drugs have been tested in clinical trials (teniposide, etoposide, paclitaxel, and epirubicin). Additionally, Table 3 showed that most of the drugs (28/33) might target TOP2A in an inhibitory manner. Among the listed drugs, paclitaxel was considered as a potential drug for cancer therapy based on its interaction with TOP2A [42]. Etoposide, another inhibitor of TOP2A, might inhibit the progress of cancer by inducing DNA damaging [43]. Additionally, TOP2A might be an important therapeutic target in etoposide resistant breast cancer [44]. Using STITCH database, we constructed an extended downstream network of TOP2A to investigate the additional effects caused by TOP2A inhibition. Our model showed that TOP2A inhibition might have possible downstream influence on DNA topoisomerase I (TOP1), DNA topoisomerase II beta (TOP2B), Ubiquitin C (UBC), Proliferating cell nuclear antigen (PCNA), Small ubiquitin-like modifier 1 (SUMO1), and SUMO2 (Additional file 1: Figure S4). According to the network, amsacrine, levofloxacin, dexrazoxane, and etoposide might act as key regulators in these processes. To the best of our knowledge, few TOP2A inhibitors have been tested for ESCC treatments. Some of them were even not to be considered as anti-cancer drugs (such as levofloxacin and dexrazoxane). The data might provide new clues for targeted therapy in ESCC patients.

Table 3 Candidate drugs targeting hub genes


Numerous researches have been performed to explore the mechanisms of ESCC during the past years, but the trends in incidence and mortality of ESCC is still increasing worldwide. Compared to the previous studies that only focused on several genes or a single cohort, this study selected 4 high-quality gene profile datasets from different research teams to integratedly explore the driven-genes and biological pathways in ESCC. Finally, we identified 146 DEGs (120 upregulated and 26 downregulated). Biological pathway enrichment analysis showed that DNA replication, cell cycle, DCC mediated signaling pathway, and Netrin-1 signaling pathway might paly crucial roles in the progression of ESCC. The PPI network was constructed with 146 nodes and 2392 edges. We then selected the top 3 significant modules from the PPI network, and these three modules were mainly related to DNA replication, cell cycle, PLK1 signaling pathway, EMT process, and ErbB1 signaling pathway, etc. According to the degree of connectivity, the top 10 genes in PPI network were considered as hub genes and they were verified in TCGA database. NDC80 was clearly associated with the poor outcome of ESCC patients. miRNA–hub gene network revealed the importance of epigenetic regulation in ESCC. Additionally, small molecular drugs found here provided new insights into the targeted therapies of ESCC.

Driven-genes play crucial roles during carcinogenesis and progression, and they usually serve a distinct biological function as a module. Using integrated bioinformatics analysis, we have identified 3 important modules. The first module (Fig. 4a) included 55 nodes, and its biological pathways (Fig. 4b) were correlated to DNA replication, cell cycle, and PLK1 signaling pathway. The aberrant cell cycle was one of marked features of tumor cells [45]. In this study, the genes related to cell cycle and mitotic regulation, such as CDK1, CCNB1, CCNB2, and NDC80, were apparently altered in patients with ESCC (Fig. 5c). Importantly, high level of NDC80 predicted poor OS in patients with ESCC. The alterations of these genes included amplification, missense, and mutation. We supposed that the above genes might be the driven-genes for ESCC development. Moreover, DNA replication stress could not only cause cell cycle abnormalities, but also accumulate genome alterations [46, 47]. In this study, the function of Module 1 was associated with the DNA replication, indicating that dysregulation of DNA replication might function as a promoter of sustained proliferation and genome instability in ESCC [46]. PLK1 was one of the most extensively studied genes in cell cycle regulation [48], and it was highly expressed in various cancers, especially in gastric cancer [49], lung cancer [50], and pancreatic carcinoma [34]. Recently, researchers have investigated the PLK1-based targeted therapy and found that PLK1 was involved in several pathways of drug resistance [51].

Module 2 (Fig. 4c) mainly consisted of collagen family members COL3A1, COL1A1, COL1A2, and COL5A2. Dysregulated expression of collagen family members was the foundation of cancer invasion and migration [52]. Many studies have demonstrated that the ectopic expression of the above genes could be the cause of cancer development, resulting in genetic mutations, epigenetic alterations, and activation of oncogenic signaling pathways or processes (such as EMT, extracellular matrix (ECM) remodeling, VEGFR3 signaling pathway, and Wnt signaling pathway, etc.) [53,54,55,56]. So far, studies have demonstrated the abnormal expression of collagen family members in several cancers [57,58,59]. However, their crucial role in ESCC was rarely mentioned [21]. In the future, in-depth studies on the roles of collagen family members in ESCC might provide new clues for inhibiting ESCC.

Module 3 (Fig. 4e) were closely related to IGF activity regulation, ErbB1 signaling pathway, and class I PI3K signaling pathway. Previously, Imsumran et al. have found that the expression level of IGF-I receptor (IGF-Ir) and IGF-II were related to the metastasis, invasion depth, and recurrence in patients with ESCC [60], indicating the potential values of using IGF members as the biomarkers for the prediction of recurrence and outcomes of ESCC patients. Moreover, miR-375 could inhibit the proliferation and migration abilities of ESCC cells through regulating the activity and expression of IGF1R [61]. ErbB1 was overexpressed and mutated in several tumors, including breast cancer [62]. The downstream signaling modules of ErbB included the PI3K/Akt signaling pathway, the phospholipase C (PLCγ) pathway and Ras/Raf/MEK/ERK1/2 pathway [63]. Fichter et al. have found that ErbB inhibitors could inhibit cell migration of ESCC cells through distinct signaling pathways (ERK1/2, Akt, STAT3, and RhoA), suggesting the powerful clues for developing ErbB targeted therapies. Since PI3K/Akt pathway also played important roles in ESCC cell growth, invasion, and migration [64, 65], we thus supposed that Module 3 was the cluster that regulated the growth and metastasis of ESCC cells.

In the study, ten genes were recognized as the hub genes, and their expression levels were all verified in the TCGA database. A list of 33 drugs with potential therapeutic efficacy against ESCC were identified. Among the 10 hub genes, the potential gene targets of the drugs are CDK1, TOP2A, CCNA2, and AURKB. In Table 3, most of drugs were TOP2A inhibitors. However, only a few of TOP2A inhibitors have been used for ESCC. More studies and clinical trials were needed to identify and explore the effective drugs for ESCC. Still, the study might push valuable insights into the individualized treatment and targeted therapy in ESCC, and the conventional drug was of potentially new use.

There were still several limitations worth mentioning in this study. First of all, we majorly explored the functions and potential roles of the hub genes without deeply analyzing the other DEGs. In the future, in-depth studies considering this field is required. Secondly, we only used TCGA data and qPCR assay to validate the expression levels of hub genes, and further experimental studies were required to demonstrate the above findings. Finally, the clinical information of ESCC patients were not deeply analyzed due to the inaccessible of data. Despite this, our study provided novel findings for ESCC studies. Compared to the single dataset studies, this study might provide more accurate results by using integrated bioinformatics analysis. Moreover, the therapeutic targets and drugs found in this study are promising and novel for personalized therapy. Additionally, we constructed the miRNA–hub gene network which might reveal the importance of epigenetic regulation in ESCC.


Using integrated bioinformatics analysis, the study identified commonly changed 146 DEGs in ESCC, which were enriched in DNA replication, cell cycle, DCC mediated pathway, and Netrin-1 signaling pathway. We also identified 10 hub genes, including CDK1, CCNB1, TOP2A, CCNB2, BUB1, CCNA2, NCAPG, AURKB, NDC80, and BUB1B, that might play important roles in ESCC. The 10 hub genes might function as novel markers and/or targets for the early cancer detection, prognostic judgment, and targeted therapy of ESCC. Additionally, a group of drugs was identified, and they could be potentially utilized for treatment of ESCC patients. And this study provided powerful basis for ESCC studies, and in-depth experimental studies were needed.

Availability of data and materials

The authors declare that the data supporting the findings of this study are available within the article.



esophageal cancer


esophageal squamous cell carcinoma


esophageal adenocarcinoma


differentially expressed genes


gene ontology


Kyoto Encyclopedia of Genes and Genomes


protein–protein interaction


Search Tool for the Retrieval of Interacting Genes


The Cancer Genome Atlas


  1. 1.

    Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: gLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394–424.

    Article  Google Scholar 

  2. 2.

    Siewert JR, Ott K. Are squamous and adenocarcinomas of the esophagus the same disease? Semin Radiat Oncol. 2007;17:38–44.

    Article  Google Scholar 

  3. 3.

    Yang W, Ma J, Zhou W, Zhou X, Cao B, Zhang H, Zhao Q, Fan D, Hong L. Molecular mechanisms and clinical implications of miRNAs in drug resistance of esophageal cancer. Expert Rev Gastroenterol Hepatol. 2017;11:1151–63.

    CAS  Article  Google Scholar 

  4. 4.

    Nik-Zainal S, Davies H, Staaf J, Ramakrishna M, Glodzik D, Zou X, Martincorena I, Alexandrov LB, Martin S, Wedge DC, Van Loo P, Ju YS, Smid M, et al. Landscape of somatic mutations in 560 breast cancer whole-genome sequences. Nature. 2016;534:47–54.

    CAS  Article  Google Scholar 

  5. 5.

    Vogelstein B, Papadopoulos N, Velculescu VE, Zhou S, Diaz LA Jr, Kinzler KW. Cancer genome landscapes. Science. 2013;339:1546–58.

    CAS  Article  Google Scholar 

  6. 6.

    Kulasingam V, Diamandis EP. Strategies for discovering novel cancer biomarkers through utilization of emerging technologies. Nat Clin Pract Oncol. 2008;5:588–99.

    CAS  Article  Google Scholar 

  7. 7.

    Guo Y, Bao Y, Ma M, Yang W. Identification of key candidate genes and pathways in colorectal cancer by integrated bioinformatical analysis. Int J Mol Sci. 2017;18:772.

    Article  Google Scholar 

  8. 8.

    Xu J, Chen Y, Zhang R, He J, Song Y, Wang J, Wang H, Wang L, Zhan Q, Abliz Z. Global metabolomics reveals potential urinary biomarkers of esophageal squamous cell carcinoma for diagnosis and staging. Sci Rep. 2016;6:35010.

    CAS  Article  Google Scholar 

  9. 9.

    Ohnuma H, Sato Y, Hayasaka N, Matsuno T, Fujita C, Sato M, Osuga T, Hirakawa M, Miyanishi K, Sagawa T, Fujikawa K, Ohi M, Okagawa Y, et al. Neoadjuvant chemotherapy with docetaxel, nedaplatin, and fluorouracil for resectable esophageal cancer: a phase II study. Cancer Sci. 2018;109:3554–63.

    CAS  Article  Google Scholar 

  10. 10.

    Okamoto H, Taniyama Y, Sakurai T, Heishi T, Teshima J, Sato C, Maruyama S, Ito K, Onodera Y, Konno-Kumagai T, Ishida H, Kamei T. Definitive chemoradiotherapy with docetaxel, cisplatin, and 5-fluorouracil (DCF-R) for advanced cervical esophageal cancer. Esophagus. 2018;15:281–5.

    Article  Google Scholar 

  11. 11.

    Zhu Y, Zhang W, Li Q, Li Q, Qiu B, Liu H, Liu M, Hu Y. A phase II randomized controlled trial: definitive concurrent chemoradiotherapy with docetaxel plus cisplatin versus 5-fluorouracil plus cisplatin in patients with oesophageal squamous cell carcinoma. J Cancer. 2017;8:3657–66.

    Article  Google Scholar 

  12. 12.

    Yan W, Shih JH, Rodriguez-Canales J, Tangrea MA, Ylaya K, Hipp J, Player A, Hu N, Goldstein AM, Taylor PR, Emmert-Buck MR, Erickson HS. Identification of unique expression signatures and therapeutic targets in esophageal squamous cell carcinoma. BMC Res Notes. 2012;5:73.

    CAS  Article  Google Scholar 

  13. 13.

    Hu N, Clifford RJ, Yang HH, Wang C, Goldstein AM, Ding T, Taylor PR, Lee MP. Genome wide analysis of DNA copy number neutral loss of heterozygosity (CNNLOH) and its relation to gene expression in esophageal squamous cell carcinoma. BMC Genomics. 2010;11:576.

    Article  Google Scholar 

  14. 14.

    Ming XY, Zhang X, Cao TT, Zhang LY, Qi JL, Kam NW, Tang XM, Cui YZ, Zhang BZ, Li Y, Qin YR, Guan XY. RHCG suppresses tumorigenicity and metastasis in esophageal squamous cell carcinoma via inhibiting NF-kappaB signaling and MMP1 expression. Theranostics. 2018;8:185–98.

    CAS  Article  Google Scholar 

  15. 15.

    Hu N, Wang C, Clifford RJ, Yang HH, Su H, Wang L, Wang Y, Xu Y, Tang ZZ, Ding T, Zhang T, Goldstein AM, Giffen C, et al. Integrative genomics analysis of genes with biallelic loss and its relation to the expression of mRNA and micro-RNA in esophageal squamous cell carcinoma. BMC Genomics. 2015;16:732.

    Article  Google Scholar 

  16. 16.

    Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, Marshall KA, Phillippy KH, Sherman PM, Holko M, Yefanov A, Lee H, Zhang N, et al. NCBI GEO: archive for functional genomics data sets-update. Nucleic Acids Res. 2013;41:D991–5.

    CAS  Article  Google Scholar 

  17. 17.

    Lebrec JJ, Huizinga TW, Toes RE, Houwing-Duistermaat JJ, van Houwelingen HC. Integration of gene ontology pathways with North American Rheumatoid Arthritis Consortium genome-wide association data via linear modeling. BMC Proc. 2009;3(Suppl 7):S94.

    Article  Google Scholar 

  18. 18.

    Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4:P3.

    Article  Google Scholar 

  19. 19.

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

    CAS  Article  Google Scholar 

  20. 20.

    Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, Lin J, Minguez P, Bork P, von Mering C, Jensen LJ. STRING v9.1: protein–protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;41:D808–15.

    CAS  Article  Google Scholar 

  21. 21.

    Wong FH, Huang CY, Su LJ, Wu YC, Lin YS, Hsia JY, Tsai HT, Lee SA, Lin CH, Tzeng CH, Chen PM, Chen YJ, Liang SC, et al. Combination of microarray profiling and protein–protein interaction databases delineates the minimal discriminators as a metastasis network for esophageal squamous cell carcinoma. Int J Oncol. 2009;34:117–28.

    CAS  PubMed  Google Scholar 

  22. 22.

    Su G, Morris JH, Demchak B, Bader GD. Biological network exploration with Cytoscape 3. Curr Protoc Bioinform. 2014;47:1–24.

    Article  Google Scholar 

  23. 23.

    Yuan L, Zeng G, Chen L, Wang G, Wang X, Cao X, Lu M, Liu X, Qian G, Xiao Y, Wang X. Identification of key genes and pathways in human clear cell renal cell carcinoma (ccRCC) by co-expression analysis. Int J Biol Sci. 2018;14:266–79.

    CAS  Article  Google Scholar 

  24. 24.

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

    CAS  Article  Google Scholar 

  25. 25.

    Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, Antipin Y, Reva B, Goldberg AP, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2:401–4.

    Article  Google Scholar 

  26. 26.

    Cotto KC, Wagner AH, Feng YY, Kiwala S, Coffman AC, Spies G, Wollam A, Spies NC, Griffith OL, Griffith M. DGIdb 3.0: a redesign and expansion of the drug-gene interaction database. Nucleic Acids Res. 2018;46:D1068–73.

    CAS  Article  Google Scholar 

  27. 27.

    Kirk J, Shah N, Noll B, Stevens CB, Lawler M, Mougeot FB, Mougeot JC. Text mining-based in silico drug discovery in oral mucositis caused by high-dose cancer therapy. Support Care Cancer. 2018;26:2695–705.

    Article  Google Scholar 

  28. 28.

    Szklarczyk D, Santos A, von Mering C, Jensen LJ, Bork P, Kuhn M. STITCH 5: augmenting protein–chemical interaction networks with tissue and affinity data. Nucleic Acids Res. 2016;44:D380–4.

    CAS  Article  Google Scholar 

  29. 29.

    Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45:D353–61.

    CAS  Article  Google Scholar 

  30. 30.

    Dibb M, Ang YS. Targeting the cell cycle in esophageal adenocarcinoma: an adjunct to anticancer treatment. World J Gastroenterol. 2011;17:2063–9.

    CAS  Article  Google Scholar 

  31. 31.

    Bonnans C, Chou J, Werb Z. Remodelling the extracellular matrix in development and disease. Nat Rev Mol Cell Biol. 2014;15:786–801.

    CAS  Article  Google Scholar 

  32. 32.

    Tran TQ, Lowman XH, Kong M. Molecular pathways: metabolic control of histone methylation and gene expression in cancer. Clin Cancer Res. 2017;23:4004–9.

    CAS  Article  Google Scholar 

  33. 33.

    Song Y, Li L, Ou Y, Gao Z, Li E, Li X, Zhang W, Wang J, Xu L, Zhou Y, Ma X, Liu L, Zhao Z, et al. Identification of genomic alterations in oesophageal squamous cell cancer. Nature. 2014;509:91–5.

    CAS  Article  Google Scholar 

  34. 34.

    Li J, Wang R, Schweickert PG, Karki A, Yang Y, Kong Y, Ahmad N, Konieczny SF, Liu X. Plk1 inhibition enhances the efficacy of gemcitabine in human pancreatic cancer. Cell Cycle. 2016;15:711–9.

    Article  Google Scholar 

  35. 35.

    Uehara H, Miyamoto M, Kato K, Cho Y, Kurokawa T, Murakami S, Fukunaga A, Ebihara Y, Kaneko H, Hashimoto H, Murakami Y, Shichinohe T, Kawarada Y, et al. Deficiency of hMLH1 and hMSH2 expression is a poor prognostic factor in esophageal squamous cell carcinoma. J Surg Oncol. 2005;92:109–15.

    CAS  Article  Google Scholar 

  36. 36.

    Dong J, Feng F, Xu G, Zhang H, Hong L, Yang J. Elevated SP/NK-1R in esophageal carcinoma promotes esophageal carcinoma cell proliferation and migration. Gene. 2015;560:205–10.

    CAS  Article  Google Scholar 

  37. 37.

    Hsu PK, Chen HY, Yeh YC, Yen CC, Wu YC, Hsu CP, Hsu WH, Chou TY. TPX2 expression is associated with cell proliferation and patient outcome in esophageal squamous cell carcinoma. J Gastroenterol. 2014;49:1231–40.

    CAS  Article  Google Scholar 

  38. 38.

    Kwon YJ, Lee SJ, Koh JS, Kim SH, Kim YJ, Park JH. Expression patterns of aurora kinase B, heat shock protein 47, and periostin in esophageal squamous cell carcinoma. Oncol Res. 2009;18:141–51.

    CAS  Article  Google Scholar 

  39. 39.

    He W, Chen L, Yuan K, Zhou Q, Peng L, Han Y. Gene set enrichment analysis and meta-analysis to identify six key genes regulating and controlling the prognosis of esophageal squamous cell carcinoma. J Thorac Dis. 2018;10:5714–26.

    Article  Google Scholar 

  40. 40.

    Ma Q. MiR-219-5p suppresses cell proliferation and cell cycle progression in esophageal squamous cell carcinoma by targeting CCNA2. Cell Mol Biol Lett. 2019;24:4.

    Article  Google Scholar 

  41. 41.

    Zhao H, Diao C, Wang X, Xie Y, Liu Y, Gao X, Han J, Li S. MiR-543 promotes migration, invasion and epithelial–mesenchymal transition of esophageal cancer cells by targeting phospholipase A2 group IVA. Cell Physiol Biochem. 2018;48:1595–604.

    CAS  Article  Google Scholar 

  42. 42.

    Srikantan S, Abdelmohsen K, Lee EK, Tominaga K, Subaran SS, Kuwano Y, Kulshrestha R, Panchakshari R, Kim HH, Yang X, Martindale JL, Marasa BS, Kim MM, et al. Translational control of TOP2A influences doxorubicin efficacy. Mol Cell Biol. 2011;31:3790–801.

    CAS  Article  Google Scholar 

  43. 43.

    Atwal M, Lishman EL, Austin CA, Cowell IG. Myeloperoxidase enhances etoposide and mitoxantrone-mediated DNA damage: a target for myeloprotection in cancer chemotherapy. Mol Pharmacol. 2017;91:49–57.

    CAS  Article  Google Scholar 

  44. 44.

    Kaplan E, Gunduz U. Expression analysis of TOP2A, MSH2 and MLH1 genes in MCF7 cells at different levels of etoposide resistance. Biomed Pharmacother. 2012;66:29–35.

    CAS  Article  Google Scholar 

  45. 45.

    Malumbres M, Barbacid M. Cell cycle, CDKs and cancer: a changing paradigm. Nat Rev Cancer. 2009;9:153–66.

    CAS  Article  Google Scholar 

  46. 46.

    Taylor EM, Lindsay HD. DNA replication stress and cancer: cause or cure? Future Oncol. 2016;12:221–37.

    CAS  Article  Google Scholar 

  47. 47.

    Tian H, Gao Z, Li H, Zhang B, Wang G, Zhang Q, Pei D, Zheng J. DNA damage response—a double-edged sword in cancer prevention and cancer therapy. Cancer Lett. 2015;358:8–16.

    CAS  Article  Google Scholar 

  48. 48.

    Lens SM, Voest EE, Medema RH. Shared and separate functions of polo-like kinases and aurora kinases in cancer. Nat Rev Cancer. 2010;10:825–41.

    CAS  Article  Google Scholar 

  49. 49.

    Otsu H, Iimori M, Ando K, Saeki H, Aishima S, Oda Y, Morita M, Matsuo K, Kitao H, Oki E, Maehara Y. Gastric cancer patients with High PLK1 expression and DNA aneuploidy correlate with poor prognosis. Oncology. 2016;91:31–40.

    CAS  Article  Google Scholar 

  50. 50.

    Xu C, Li S, Chen T, Hu H, Ding C, Xu Z, Chen J, Liu Z, Lei Z, Zhang HT, Li C, Zhao J. miR-296-5p suppresses cell viability by directly targeting PLK1 in non-small cell lung cancer. Oncol Rep. 2016;35:497–503.

    CAS  Article  Google Scholar 

  51. 51.

    Gutteridge RE, Ndiaye MA, Liu X, Ahmad N. Plk1 inhibitors in cancer therapy: from laboratory to clinics. Mol Cancer Ther. 2016;15:1427–35.

    CAS  Article  Google Scholar 

  52. 52.

    Wolf K, Alexander S, Schacht V, Coussens LM, von Andrian UH, van Rheenen J, Deryugina E, Friedl P. Collagen-based cell migration models in vitro and in vivo. Semin Cell Dev Biol. 2009;20:931–41.

    CAS  Article  Google Scholar 

  53. 53.

    Jablonska-Trypuc A, Matejczyk M, Rosochacki S. Matrix metalloproteinases (MMPs), the main extracellular matrix (ECM) enzymes in collagen degradation, as a target for anticancer drugs. J Enzyme Inhib Med Chem. 2016;31:177–83.

    CAS  Article  Google Scholar 

  54. 54.

    Jha SK, Rauniyar K, Karpanen T, Leppanen VM, Brouillard P, Vikkula M, Alitalo K, Jeltsch M. Efficient activation of the lymphangiogenic growth factor VEGF-C requires the C-terminal domain of VEGF-C and the N-terminal domain of CCBE1. Sci Rep. 2017;7:4916.

    Article  Google Scholar 

  55. 55.

    Dettman RW, Simon HG. Rebooting the collagen gel: artificial hydrogels for the study of epithelial mesenchymal transformation. Dev Dyn. 2018;247:332–9.

    CAS  Article  Google Scholar 

  56. 56.

    Willis CM, Kluppel M. Chondroitin sulfate-E is a negative regulator of a pro-tumorigenic Wnt/beta-catenin-Collagen 1 axis in breast cancer cells. PLoS ONE. 2014;9:e103966.

    Article  Google Scholar 

  57. 57.

    Hayashi M, Nomoto S, Hishida M, Inokawa Y, Kanda M, Okamura Y, Nishikawa Y, Tanaka C, Kobayashi D, Yamada S, Nakayama G, Fujii T, Sugimoto H, et al. Identification of the collagen type 1 alpha 1 gene (COL1A1) as a candidate survival-related factor associated with hepatocellular carcinoma. BMC Cancer. 2014;14:108.

    Article  Google Scholar 

  58. 58.

    Li J, Ding Y, Li A. Identification of COL1A1 and COL1A2 as candidate prognostic factors in gastric cancer. World J Surg Oncol. 2016;14:297.

    Article  Google Scholar 

  59. 59.

    Yu Y, Liu D, Liu Z, Li S, Ge Y, Sun W, Liu B. The inhibitory effects of COL1A2 on colorectal cancer cell proliferation, migration, and invasion. J Cancer. 2018;9:2953–62.

    Article  Google Scholar 

  60. 60.

    Imsumran A, Adachi Y, Yamamoto H, Li R, Wang Y, Min Y, Piao W, Nosho K, Arimura Y, Shinomura Y, Hosokawa M, Lee CT, Carbone DP, et al. Insulin-like growth factor-I receptor as a marker for prognosis and a therapeutic target in human esophageal squamous cell carcinoma. Carcinogenesis. 2007;28:947–56.

    CAS  Article  Google Scholar 

  61. 61.

    Kong KL, Kwong DL, Chan TH, Law SY, Chen L, Li Y, Qin YR, Guan XY. MicroRNA-375 inhibits tumour growth and metastasis in oesophageal squamous cell carcinoma through repressing insulin-like growth factor 1 receptor. Gut. 2012;61:33–42.

    CAS  Article  Google Scholar 

  62. 62.

    Wang Z. ErbB Receptors and Cancer. Methods Mol Biol. 2017;1652:3–35.

    CAS  Article  Google Scholar 

  63. 63.

    Roskoski R Jr. The ErbB/HER family of protein-tyrosine kinases and cancer. Pharmacol Res. 2014;79:34–74.

    CAS  Article  Google Scholar 

  64. 64.

    Chen J, Lan T, Zhang W, Dong L, Kang N, Fu M, Liu B, Liu K, Zhang C, Hou J, Zhan Q. Dasatinib enhances cisplatin sensitivity in human esophageal squamous cell carcinoma (ESCC) cells via suppression of PI3K/AKT and Stat3 pathways. Arch Biochem Biophys. 2015;575:38–45.

    CAS  Article  Google Scholar 

  65. 65.

    Tian B, Chen X, Zhang H, Li X, Wang J, Han W, Zhang LY, Fu L, Li Y, Nie C, Zhao Y, Tan X, Wang H, et al. Urokinase plasminogen activator secreted by cancer-associated fibroblasts induces tumor progression via PI3K/AKT and ERK signaling in esophageal squamous cell carcinoma. Oncotarget. 2017;8:42300–13.

    PubMed  PubMed Central  Google Scholar 

Download references


We thank the GEO, DAVID, STRING, GEPIA2, and DGIdb databases for providing their platforms and contributors for their valuable data sets.


This study was supported in part by Grant from the National Natural Scientific Foundation of China (No. 81171923), Grant from the Scientific Foundation of Shaanxi Province (No. S2019-YF-ZDCXL-ZDLSF-0134), Grant from the National Clinical Research Center for Digestive Diseases (No. 2015BAI13B07) and Grant from the State Key Laboratory of Cancer Biology (No. CBSKL2014Z13).

Author information




WY, XZ, YH, and LD contributed equally to this work, and should be considered as the co-first authors. WY analyzed the data and drafted the manuscript; XZ, YH, and LD collected and analyzed the data. LH designed and supervised the study; XL, XW, and YZ revised the manuscript; WZ and JL made the diagrams and tables of the article. HZ, QZ, and DF reviewed and revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yu Han or Liu Hong.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Additional file

Additional file 1. Figure. S1:

Expression heatmap of the 146 DEGs. Figure S2: Construction of PPI network for 146 DEGs in ESCC. Figure S3: Construction of PPI network for 146 DEGs and their related genes in ESCC. Figure S4: Targetable TOP2A subnetwork. Figure S5: Relative mRNA expression levels of the 10 hub genes in EC109 cells compare to that in Het-1A. Table S1: The 10 hub genes and corresponding primer sets. Table S2: Prognostic information of the 10 hub genes in ESCC patients.

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yang, W., Zhao, X., Han, Y. et al. Identification of hub genes and therapeutic drugs in esophageal squamous cell carcinoma based on integrated bioinformatics strategy. Cancer Cell Int 19, 142 (2019).

Download citation


  • Esophageal squamous cell carcinoma
  • Bioinformatics
  • Hub genes
  • Cell cycle
  • Differentially expressed genes
  • Drug