Abstract
Background:
Small cell carcinoma of the esophagus (SCCE) is a rare form of esophageal cancer, which also belongs to the category of neuroendocrine malignant tumors, with a low incidence but high aggressiveness, and a very poor prognosis for patients. Currently, there is a lack of unique staging and treatment guidelines for SCCE. Therefore, a deeper understanding of the therapeutic targets and the mechanisms underlying its occurrence and development is of great importance for early diagnosis, identification of potential therapeutic agents and improvement of the prognosis for patients.
Methods:
Firstly, the dataset of SCCE was downloaded from the GEO database and GEO2R tool was employed for the analysis of differentially expressed genes (DEGs). GO and KEGG analysis of DEGs were carried out by using the Bioinformatics and OmicStudio tools. Then, up- and down-regulated genes were intersected with the oncogenes and the tumor suppressor genes respectively, to obtain the differentially expressed onco/tumor suppressor genes in SCCE. The STRING database was employed to conduct protein-protein interaction (PPI) analysis of differentially expressed onco/tumor suppressor genes, the network was further constructed in Cytoscape, and hub genes of the network were obtained through the Cytohubba plugin. In addition, miRDB, miRwalk, Targetscan, OncomiR, starbase, Lncbase were used to predict miRNAs and lncRNAs that regulate hub genes, the ceRNA network was built based on this. Transcription factor-miRNA co-regulatory network was analyzed in the NetworkAnalyst database and embellished in Cytoscape. Finally, drugs that may target hub genes were searched through the DGIdb and ConnectivityMAP, and docking verification was performed using Schrodinger.
Results:
A total of 820 genes were upregulated and 716 were downregulated, of these, 54 were identified as oncogenes and 85 as tumor suppressor genes. Seven hub genes were identified from the PPI network, which were AURKA, BIRC5, CDK1, EZH2, FOXM1, KLF4 and UBE2C. Furthermore, a total of 38 drugs were searched and filtered in DGIdb and ConnectivityMAP, in the molecular docking results of drugs with hub genes, the docking score of AURKA, CDK1, and EZH2 with multiple drugs were low (<6). In addition, crizotinib with AURKA, lapatinib with CDK1, rucaparib with EZH2, rucaparib with UBE2C were the lowest energy of all molecular docking results.
Conclusion:
AURKA, BIRC5, CDK1, EZH2, FOXM1, KLF4 and UBE2C are the hub genes of SCCE, among them, AURKA, CDK1 and EZH2 may be used as targets of multiple drugs. Crizotinib, lapatinib, and rucaparib can act on the above targets to inhibit the progression of SCCE and play a therapeutic role.
1 Introduction
Small cell carcinoma of the esophagus (SCCE) is a rare histological type of esophageal cancer, which belongs to the category of neuroendocrine malignant tumors, and accounts for 0.5%–5.9% of all esophageal cancers in China (Law et al., 1994; Chen et al., 2011; Li et al., 2010a; Li et al., 2010b). Some existing studies (Chen et al., 2011; Hudson et al., 2007; Ku et al., 2008) have shown that SCCE has rapid progression, high metastasis, and extremely poor prognosis. Analysis of clinical data revealed that the tumor size, invasion depth, and metastasis rate of SCCE are significantly higher than those of esophageal squamous cell carcinoma and esophageal adenocarcinoma of the same period (Cai et al., 2019). Most patients die within 2 years of diagnosis, and their median survival is only 8–13 months (Wang et al., 2018). Clinically, SCCE is mainly diagnosed by imaging examinations and endoscopic pathological tissue biopsy, and its immunohistochemistry can detect endocrine markers (ZHANG et al., 2019). The occurrence and development of SCCE is a multi-factor, multi-stage dynamic process. Its specific mechanism is still unclear, but existing studies have revealed some relevant aspects. PTEN is a common tumor suppressor gene in the human body. It can inhibit the PI3K-AKT signaling pathway to cause cells to arrest in the G1 phase or induce apoptosis of tumor cells. Zhang et al. detected the mutation of genes such as EGFR, KRAS, PIK3CA and PTEN, and found that the incidence of PTEN mutation in Chinese PSCCE patients is higher than that of other esophageal cancer histological subtypes (Zhang et al., 2014). In addition, the functional status of PTEN will also affect the tumor immune microenvironment. In terms of cell origin, SOX2 is a transcription factor that maintains the pluripotency of tumor stem cells, studies have found that in PSCCE, SCLC and embryonic esophageal tissue, SOX2 gene is overexpressed and tumor suppressor gene Rb1 is expressed at a low level, but in poorly differentiated squamous cell carcinoma, the opposite is true, suggesting that PSCCE may be derived from embryonic stem cells. Mutation of the Rb1 gene is one of the early molecular changes in PSCCE (Ishida, 2017). In addition, the occurrence of SCCE is also related to DNA damage repair. Some studies have found that PAK1 is overexpressed in PSCCE and is positively correlated with the DNA damage marker γH2AX, suggesting that PAK1 may be involved in the DNA damage/repair process of PSCCE and promote the invasion of SCCE cells. Moreover, overactivation of PAK1 is closely related to tumor location, lymph node metastasis, and overall survival. The OS of patients with overactivated PAK1 will also be significantly reduced. PAK1 has the potential to become a direct target for the treatment of PSCCE (Gan et al., 2015). Through the study of genomic expression profiles, it was found that mutations in TP53, RB1 and NOTCH family are widely present in PSCCE and SCLC, while mutations in PDE3A, CBLN3 and PTPRM3 genes have never been reported in esophageal cancer. Mutations in these three genes may be specific molecular markers in the development of PSCCE. Changes in Wnt pathway and NOTCH signaling pathway are also important events in the development of SCCE (Wang et al., 2018). In recent years, studies on the tumor immune microenvironment have also found that the immune checkpoint receptor TIGIT and its ligand CD155 of SCCE are widely upregulated in PSCCE and are significantly related to distant metastasis, Ki-67 index, prognosis (Zhao K. et al., 2020).
Due to the low prevalence of SCCE, there is a lack of clinically relevant studies, usually retrospective studies or a few case-report studies, and there is no unique staging and grading system and treatment guidelines for reference. Some clinical data analyses have shown that patients with SCCE have difficulty benefiting from surgical treatment alone, which may be related to the high metastasis rate accompanying their diagnosis (Hudson et al., 2007; Cai et al., 2019). Compared with surgery, radiotherapy and chemotherapy are more meaningful for the whole treatment of SCCE, especially for patients with extensive lesions (Hudson et al., 2007; Xiao et al., 2019), which can significantly prolong the survival of patients. In the treatment plan for patients with extensive lesions, chemotherapy should be used as the basis of the entire treatment (Song et al., 2009; Vos et al., 2011; Ding et al., 2013; Xie et al., 2015), while for patients with limited lesions, radical surgical resection combined with postoperative radiotherapy and chemotherapy are the main treatment tools (Chen et al., 2011; Zhu et al., 2014). SCCE and small cell lung cancer have many similar histological and clinical pathological features, so internationally, its staging and treatment strategies are mostly based on the staging and treatment protocol of small cell lung cancer. The treatment guidelines for small cell lung cancer are generally platinum-based, in combination with other commonly used chemotherapy agents, such as etoposide or irinotecan (Amarasena et al., 2015; Ito et al., 2021). However, chemotherapy drugs are usually accompanied with tumor resistance and their effectiveness is reduced, and there are very few studies on emerging targeted therapy and immunotherapy in SCCE. In terms of treatment, it refers to the treatment experience of small cell lung cancer and neuroendocrine tumors (Cicin et al., 2007). More experimental results are still needed to confirm their effectiveness in SCCE. In summary, exploring the occurrence and development mechanism of SCCE, clarifying its tumor biological characteristics at different stages and its prognostic factors are of great significance for finding new drugs to treat SCCE or overcome chemotherapy resistance and then formulate the optimal comprehensive treatment strategy.
Since SCCE is extremely invasive, progresses rapidly, and has limited treatment options, early diagnosis and early treatment are essential. Many databases now share disease-related microarray data, which can identify gene expression profiles under different disease states. Biomarkers of a specific disease are often expressed differently in different physiological and pathological states, especially in disease states, which are extremely different from normal physiological states. Therefore, finding and evaluating these genes based on expression profiles may be an effective way for diagnosis and monitoring the prognosis of the disease (Huang et al., 2010). In addition, taking corresponding measures against hub genes in SCCE is a powerful method to manage and alleviate disease prognosis. Therefore, this study aims to mine specific markers for SCCE through gene expression data in the database, and to find suitable candidate drugs for hub genes through systems biology methods.
2 Materials and methods
2.1 Acquisition of differentially expressed genes in SCCE
Searching for “Small cell carcinoma of the esophagus” in the Gene Expression Omnibus (GEO) database, only one dataset, GSE111299, was sequenced from the genomes of patients with SCCE, which was eligible to be screened in this study (Liu et al., 2018). The platform tool GEO2R was used for online analysis, the GPL570 platform was selected and the samples were grouped into “SCCE” and “ctrl” according to their source, with three samples each. After online analysis, a complete table of expression of genes was obtained. The conditions for screening differentially expressed genes (DEGs) were: |log2FC| ≥ 2, P. Value < 0.01, removing the blanks and duplicates of “Gene. Symbol”, if there is more than one expression value for a gene, the maximum value is taken as the final value. Finally, the list of DEGs was obtained, the volcano map and heat map were performed using the OmicStudio tools (https://www.omicstudio.cn/tool) (Lyu et al., 2023).
2.2 GO and KEGG enrichment analysis of DEGs in SCCE
In order to elucidate the biological functions and pathways associated with the identified DEGs, a functional enrichment analysis was conducted using the Gene Ontology (GO) database. Additionally, the Kyoto Encyclopedia of Genes and Genomes database (KEGG) database was employed to examine the signaling pathways with significantly enriched differential genes. Enrichment analysis was performed by https://www.bioinformatics.com.cn, an online platform designed for analysis and visualization of data (Tang et al., 2023). P < 0.05 in the enrichment results is considered significant. The GO enrichment analysis diagram was automatically generated and visualized in the bioinformatics platform. The top10 signal pathways were selected based on the P value of the KEGG enrichment analysis results and visualized in the OmicStudio tools.
2.3 Onco and tumor suppressor genes and their PPI network in SCCE
Human oncogenes and tumor suppressor genes (TSGs) were downloaded from oncogene database (https://bioinfo-minzhao.org/ongene/) (Liu et al., 2017) and TUMOR SUPPRESSOR GENE DATABASE (https://bioinfo.uth.edu/TSGene/index.html) (Zhao et al., 2016), respectively. Then VENNY2.1.0 was used to intersect oncogenes with upregulated genes and TSGs with downregulated genes, to gain upregulated oncogenes and downregulated TSGs in SCCE. Next, the PPI network of oncogenes and TSGs in SCCE was built and visualized in a search tool for the retrieval of interacting genes/proteins (STRING) (https://cn.string-db.org/) version 12.0 and Cytoscape (https://cytoscape.org/) version 3.10.0.
2.4 Identification of hub genes from PPI network
In order to further screen hub genes in SCCE, the cytohubba plugin (https://apps.cytoscape.org/apps/cytohubba) was used in the cytoscape and its’ 12 topological analysis methods (MCC, DMNC, MNC, Degree, EPC, BottleNeck, EcCentricity, Closeness, Radiality, Betweenness, Stress, ClusteringCoefficient) were utilized to calculate. The top 10 genes in each method were obtained, and the genes that were in the top 10 in six or more algorithms were selected and considered to be hub genes in SCCE.
2.5 Survival analysis of hub genes
To further verify the obtained hub genes, we used the GEPIA database for gene expression analysis and survival analysis (Tang et al., 2017). Data from TCGA and GETx were used to first explore the expression of hub genes in ESCA, and then explore the association between survival probability and the expression of hub genes. We focused on ESCA patients and evaluated overall survival (OS) and disease free survival (DFS). Plots with a median cutoff of 50% and hazard ratios with 95% confidence intervals were generated, with statistical significance both set at P < 0.05.
2.6 Construction of ceRNA and TF-miRNA co-regulatory network
In order to explore the regulatory network related to hub genes and better clarify the role of hub genes, the upstream miRNAs of seven hub genes were predicted through three databases: miRwalk (http://mirwalk.umm.uni-heidelberg.de/), miRDB (https://mirdb.org/), and Targetscan (https://www.targetscan.org/vert_80/), and the overlapping data of the three were taken; in the OncomiR (https://oncomir.org/oncomir) database, select ESCA as a cancer type to obtain relevant miRNAs. The significant miRNAs that appear in the four databases at the same time were considered to be the miRNAs that regulate hub genes. The filtered miRNAs were entered into starbase (starBase or ENCORI: Decoding the Encyclopedia of RNA Interactomes (rnasysu.com)) and Lncbase (https://dianalab.e-ce.uth.gr/html/diana/web/index.php?r=lncbasev2) for search, the results obtained were paired one by one with miRNA-lncRNA into a set of data, duplicates were removed. After sorting the correlation data of hub genes, miRNAs, and lncRNAs, the ceRNA network visualization analysis was performed in Cytoscape.
NetworkAnalyst (https://www.networkanalyst.ca/) version 3.0 can embed the target gene into the biological network of 15 different databases and extract the genes/proteins, miRNAs, drugs, chemicals or diseases that are most closely related to the uploaded genes (seed genes) (Zhou et al., 2019). Enter the hub genes’ name in the NetworkAnalyst, and select TF-miRNA coregulatory interaction database, the TF-miRNA co-regulatory network related to the hub genes was obtained after analysis, then visualized in Cytoscape.
2.7 Construction of GeneMANIA-based functional association network
To further identify genes related to the functions of hub genes, we used a novel approach from the GeneMANIA database to find additional 10 genes with the strongest associations for each hub gene (Rajadnya et al., 2024). This analysis selected three key parameters, namely, co-expression, gene interactions, and physical interactions, to improve the accuracy of target identification. Subsequently, all newly discovered genes were merged with the initially obtained hub genes to create a GMFA-based SCCE target expansion database network. Finally, GO and KEGG enrichment analysis were performed on these expanded genes.
2.8 Screening for drug candidates
DGIdb (https://www.dgidb.org/) version 5.0.6 is an open source search engine for drug-gene interactions and druggable genomes (Cannon et al., 2024). In order to find potential therapeutic drugs for SCCE, a search was conducted in the DGIdb. Hub genes were used as input to predict the drugs targeting them, anti-tumor drugs approved by the FDA were screened as drug candidates in this study.
Based on information about differences in gene expression, the ConnectivityMAP (https://clue.io/) uses a pattern matching algorithms to infer whether there is a functional correlation between drugs, genes and diseases. Its query application can be used to find positive and negative connections between the gene expression features of interest and all features in ConnectivityMAP (Lamb, 2007). The upregulated oncogenes and downregulated TSGs were input and queried in the Query tool of ConnectivityMAP. In the results, a negative value of “raw_cs” indicates that the compound may be effective against the disease. Therefore, the top 10 compounds were selected as drug candidates based on the sorting and screening of FDA-approved anti-tumor drugs.
2.9 Molecular docking
The structures of drug candidates were sourced from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/). Firstly, imported these drug structures into the Schrodinger software version 2024.01 to establish a database, following hydrogenation, structural optimization, and energy minimization, the resulting ligand molecules were saved for use in molecular docking. The target protein AURKA (PDB ID: 4UYN), BIRC5 (PDB ID: 7LBQ), CDK1 (PDB ID: 6GU2), EZH2 (PDB ID: 5HYN), FOXM1 (PDB ID: 3G73), KLF4 (PDB ID: 6VTX), UBE2C (PDB ID: 4YII) structures were sourced from the RCSB database (https://www.rcsb.org/), then protein structures were processed on the Maestro (version 2024.01) platform, and the proteins were processed using the ProteinPreparationWizard of Schrodinger. Finally, the proteins were minimized and the geometry structures were optimized. The receptor underwent pre-processing, optimization, and minimization (with the OPLS4 force field being utilized for constrained minimization). The LigPrep module’s default settings were followed in preparing the compound structures. The constructed receptor was imported, the protein’s natural ligand was selected as the target’s active site center, and the active site position was generated in the receptor grid for screening in the Glide module. Lastly, the molecular docking was carried out using the recognized standard docking methodology (Standard Precision, SP) (Friesner et al., 2004; Rajeswari et al., 2014; Fazi et al., 2015). The Schrödinger software support was provided by Nanjing University of Chinese Medicine.
3 Results
3.1 Identification and enrichment analysis of differentially expressed genes
The process of this study is shown in Figure 1. In the GSE111299 dataset, a total of 1,536 differentially expressed genes in SCCE (|log2FC| ≥ 2 and P. Value < 0.01) were screened, 820 genes were upregulated and 716 genes were downregulated. The volcano plot shows the distribution of these DEGs (Figure 2A), and the clustering heat map of DEGs is shown in Figure 2B, GSM3020874, GSM3020876, and GSM3020878 are tumor tissues of patients with SCCE, while GSM3020875, GSM3020877, and GSM3020879 are normal tissues of the corresponding patients. It can be seen that the gene expression of tumor and normal tissues is significantly different.
FIGURE 1
FIGURE 2
To gain insight into the biological classification of the identified DEGs, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed. Figure 2C shows the enrichment analysis results of biological processes, which show that DEGs are particularly enriched in aspects such as mitotic nuclear division, organelle fission, chromosome segregation, etc. DEGs also significantly enriched cellular components such as chromosome, centromeric region, spindle, kinetochore. In terms of molecular functions, DEGs are mainly related to microtubule binding, tubulin binding, and DNA helicase activity. etc. Bar length represents the enrichment score of each GO term. Figure 2C shows the top 10 results for biological processes (BP), cellular components (CC), and molecular functions (MF).
In the KEGG pathway analysis, several important pathways were identified, including cell cycle, Chemical carcinogenesis-receptor activation, p53 signaling pathway, and so on. The left side of Figure 2D shows the interaction between representative genes and these 10 KEGG pathways in the form of a Sankey diagram. The size of the bubble chart circle on the right side indicates the number of differentially enriched genes in the pathway, and the color of the circle indicates the enrichment p-value of DEGs in each pathway. The darker the color, the higher the enrichment p-value. The X-axis indicates the enrichment score of each pathway.
3.2 Onco and tumor suppressor genes in SCCE and protein-protein interaction
The intersections between upregulated DEGs and oncogenes and between downregulated DEGs and TSGs were analyzed using Venny 2.1.0. The results are shown in Figures 3A, B, among the 820 upregulated genes and 716 downregulated genes, there were 54 upregulated oncogenes and 85 downregulated TSGs. The protein interactions of upregulated oncogenes and downregulated TSGs showed 136 interaction nodes, 407 edges, an average node degree of 5.99, and an average local clustering coefficient of 0.411 (Figures 3C, D). The protein interaction network was visualized using Cytoscape and sorted by degree value. The larger the circle and the darker the color, the higher degree of the protein. It is in a key position in this interaction network and is the core of the entire network. It can be considered as a key target for SCCE.
FIGURE 3
3.3 Identification of hub genes from PPI network
Further analysis of the PPI network by 12 different topological analysis methods found that seven genes (AURKA, BIRC5, CDK1, EZH2, FOXM1, KLF4, UBE2C) appeared repeatedly at top 10 among at least 50% of the topological analysis methods (Figure 4; Table 1). They are considered to be the hub genes of SCCE, which may have certain value in the diagnosis, treatment and prognosis of SCCE. However, the effectiveness of these hub genes still needs further experimental verification.
FIGURE 4
TABLE 1
| Different topological analysis | Genes | Full name | Function |
|---|---|---|---|
| MCC, MNC, Degree, EPC, BottleNeck, Closeness | AURKA | Aurora Kinase A | The protein encoded by this gene is a cell cycle-regulated kinase that appears to be involved in microtubule formation and/or stabilization at the spindle pole during chromosome segregation. The encoded protein is found at the centrosome in interphase cells and at the spindle poles in mitosis |
| MCC, MNC, Degree, EPC, Closeness, Radiality | BIRC5 | Baculoviral IAP Repeat Containing 5 | This gene is a member of the inhibitor of apoptosis (IAP) gene family, which encode negative regulatory proteins that prevent apoptotic cell death. IAP family members usually contain multiple baculovirus IAP repeat (BIR) domains, but this gene encodes proteins with only a single BIR domain |
| MCC, MNC, Degree, EPC, BottleNeck, EcCentricity, Closeness, Radiality, Betweenness, Stress | CDK1 | Cyclin-dependent kinase 1 | The protein encoded by this gene is a member of the Ser/Thr protein kinase family. This protein is a catalytic subunit of the highly conserved protein kinase complex known as M-phase promoting factor (MPF), which is essential for G2/M phase transitions of eukaryotic cell cycle. Mitotic cyclins stably associate with this protein and function as regulatory subunits |
| MCC, MNC, Degree, EPC, BottleNeck, Closeness, Radiality, Betweenness, Stress | EZH2 | Enhancer Of Zeste 2 Polycomb Repressive Complex 2 Subunit | This gene encodes a member of the Polycomb-group (PcG) family. PcG family members form multimeric protein complexes, which are involved in maintaining the transcriptional repressive state of genes over successive cell generations. This protein associates with the embryonic ectoderm development protein, the VAV1 oncoprotein, and the X-linked nuclear protein |
| MCC, MNC, Degree, EPC, Closeness, Radiality | FOXM1 | Forkhead Box M1 | The protein encoded by this gene is a transcriptional activator involved in cell proliferation. The encoded protein is phosphorylated in M phase and regulates the expression of several cell cycle genes, such as cyclin B1 and cyclin D1 |
| Degree, EPC, Closeness, Radiality, Betweenness, Stress | KLF4 | KLF Transcription Factor 4 | This gene encodes a protein that belongs to the Kruppel family of transcription factors. The encoded zinc finger protein is required for normal development of the barrier function of skin. The encoded protein is thought to control the G1-to-S transition of the cell cycle following DNA damage by mediating the tumor suppressor gene p53 |
| MCC, MNC, Degree, EPC, Closeness, Radiality | UBE2C | Ubiquitin Conjugating Enzyme E2 C | The modification of proteins with ubiquitin is an important cellular mechanism for targeting abnormal or short-lived proteins for degradation. Ubiquitination involves at least three classes of enzymes: ubiquitin-activating enzymes, ubiquitin-conjugating enzymes, and ubiquitin-protein ligases. This gene encodes a member of the E2 ubiquitin-conjugating enzyme family |
Hub genes identified via the different topological analysis and their functional roles.
3.4 Survival analysis of hub genes
Since we lack clinical survival data for SCCE, and considering that SCCE is a subtype of esophageal cancer, we used the data of ESCA in TCGA and GETx for analysis. Using GEPIA, we first analyzed the expression of seven hub genes in ESCA and found that their expression was roughly similar to that in SCCE, and there was no statistically significant difference in the expression of EZH2 and KLF4 (Figure 5A). The overall survival (OS) and disease free survival (DFS) of the hub genes were then analyzed, with a median cutoff value of 50% and a 95% confidence interval for the hazard ratio. In the OS analysis, with the progression of ESCA, high expression of BIRC5, EZH2, FOXM1, and UBE2C promoted patient death, while high expression of KLF4 prolonged patient survival (Figure 5B). In DFS analysis, high expression of BIRC5, EZH2 and UBE2C will reduce the survival rate of patients, on the contrary, high expression of KLF4 will significantly improve the survival rate of patients, with significant statistical significance (p (HR) = 0.0087) (Figure 5C). These results indicate that the expression differences of the hub genes are largely associated with the prognosis of ESCA patients, further suggesting that hub genes have a certain significance for the prognosis of SCCE.
FIGURE 5
3.5 Construction of ceRNA and TF-miRNA co-regulatory network
In order to better understand the hub genes and their associated regulatory mechanisms, relevant miRNAs, lncRNAs (Supplementary Figure S1), and transcription factors were predicted from the database, the ceRNA and TF-mRNA regulatory networks based on hub genes were constructed, as shown in Figure 6, including 15 miRNAs (Table 2) and 27 relevant lncRNAs. These miRNAs, lncRNAs and transcription factors may have certain regulatory capabilities for the expression or function of hub genes, which is worthy of further study.
FIGURE 6
TABLE 2
| miRNA name | Cancer abbreviation | T-Test P-value | t-Test FDR | Upregulated in | Tumor Log2 mean expression | Normal Log2 mean expression |
|---|---|---|---|---|---|---|
| hsa-miR-30c-2-3p | ESCA | 0.0237 | 0.082 | Normal | 4.40 | 5.19 |
| hsa-miR-153-3p | ESCA | 0.029 | 0.0975 | Normal | 0.74 | 2.06 |
| hsa-miR-877-5p | ESCA | 0.0000274 | 0.00129 | Tumor | 1.56 | 0.27 |
| hsa-miR-103a-3p | ESCA | 0.0000615 | 0.00215 | Tumor | 14.53 | 13.93 |
| hsa-miR-1301-3p | ESCA | 0.000128 | 0.00279 | Tumor | 3.11 | 1.58 |
| hsa-miR-29b-1-5p | ESCA | 0.000523 | 0.00618 | Tumor | 2.79 | 1.48 |
| hsa-miR-204-5p | ESCA | 0.00215 | 0.0155 | Normal | 0.63 | 3.94 |
| hsa-miR-330-5p | ESCA | 0.00545 | 0.0289 | Tumor | 4.43 | 3.81 |
| hsa-miR-542-3p | ESCA | 0.0076 | 0.0385 | Tumor | 7.35 | 6.54 |
| hsa-miR-767-5p | ESCA | 0.0107 | 0.0488 | Tumor | 3.17 | 0.00 |
| hsa-miR-423-5p | ESCA | 0.0113 | 0.0506 | Tumor | 6.08 | 5.64 |
| hsa-miR-34a-5p | ESCA | 0.0218 | 0.0787 | Tumor | 7.17 | 6.49 |
| hsa-miR-320b | ESCA | 0.0231 | 0.0816 | Tumor | 2.46 | 1.64 |
| hsa-miR-330-3p | ESCA | 0.0251 | 0.0861 | Tumor | 1.30 | 0.82 |
| hsa-miR-135a-5p | ESCA | 0.0347 | 0.113 | Normal | 0.87 | 2.30 |
Differential expression analysis of the predicted miRNAs in OncomiR.
3.6 Construction of GeneMANIA-based functional association network and enrichment analysis
To further explore genes associated with hub gene functions, a novel analysis method of GeneMANIA was used to identify 10 additional genes for each hub gene (Figure 7A), which were organized into the GMFA extended database (GMFA-ED), which contains 67 genes after removing duplicates. The GMFA method integrates co-expression, genetic interaction, and physical interaction parameters to capture various genes that are strongly associated with hub genes and form a network diagram, as shown in Figure 7B. GO and KEGG enrichment analysis of GMFA-ED found that these genes were significantly enriched in nuclear division, mitotic nuclear division, organelle fission and the regulation of cell cycle in BP, significantly enriched in spindle, ESC/E(Z) complex, chromosome, centromeric region in CC, and significantly enriched in kinase regulator activity, protein kinase regulator activity, protein kinase activator activity, et al. in MF (Figure 7C). The pathway enrichment results showed that, like the previous enrichment, the most significantly enriched was still Cell cycle, but more surprisingly, the “FoxO signaling pathway” and “Small cell lung cancer” were enriched (Figure 7D). Based on the similarities between small cell lung cancer and SCCE, this enrichment result confirms that these genes are likely to play important roles in SCCE.
FIGURE 7
3.7 Searching for drug candidates targeting hub genes
There are many databases containing information on drug action, which can be used as a basis for predicting drugs that may target certain genes. Through conditional screening, it may lead to the discovery of new drugs for the treatment of SCCE. This study used the DGIdb database to search drugs that targeting hub genes. Rucaparib for targeting CDK1, dabrafenib and tazemetostat hydrobromide for targeting EZH2, and vorinostat, omacetaxine mepesuccinate, paclitaxel, epirubicin, doxorubicin, erlotinib, dactinomycin, lapatinib, cisplatin, imatinib, irinotecan hydrochloride, plicamycin, carboplatin, lonafarnib, trastuzumab, cytarabine, arsenic trioxide, flutamide, docetaxel, INFβ-1A, romidepsin, fluorouracil for targeting BIRC5, and that targeted AURKA are paclitaxel, fluorouracil, cisplatin, entrectinib, sorafenib, pazopanib, tamoxifen, and drug targeting KLF4 is hydroxyurea. In addition, the ConnectivityMAP can query potentially effective drugs based on DEGs, upregulated oncogenes and downregulated TSGs in SCCE were imported for query. After filtering, drug candidates were selected including etoposide, valrubicin, pralatrexate, KPT-330 (selinexor), daunorubicin, teniposide, lonafarnib, irinotecan, cladribine, and crizotinib, as shown in Table 3. All of the above drugs were considered as therapeutic drug candidates for SCCE, and their effectiveness can be further verified.
TABLE 3
| Databases | Genes | Drugs |
|---|---|---|
| DGIdb | AURKA | Paclitaxel, Fluorouracil, Cisplatin, Entrectinib, Sorafenib, Pazopanib, Tamoxifen |
| BIRC5 | Vorinostat, Omacetaxine mepesuccinate, Paclitaxel, Epirubicin, Doxorubicin, Erlotinib, Dactinomycin, Lapatinib, Cisplatin, Imatinib, Irinotecan hydrochloride, Plicamycin, Carboplatin, Lonafarnib, Trastuzumab, Cytarabine, Arsenic trioxide, Flutamide, Docetaxel, INFβ-1A, Romidepsin, Fluorouracil | |
| CDK1 | Rucaparib | |
| EZH2 | Dabrafenib, Tazemetostat hydrobromide | |
| KLF4 | Hydroxyurea | |
| ConnectivityMap | 54 upregulated oncogenes and 85 downregulated TSGs | Etoposide, Valrubicin, Pralatrexate, KPT-330 (Selinexor), Daunorubicin, Teniposide, Lonafarnib, Irinotecan, Cladribine, Crizotinib |
Potential drug candidates analyzed by the DGIdb and ConnectivityMap (FDA approved anti-tumor drugs).
3.8 Molecular docking of hub genes with drug candidates
In order to ascertain the binding ability between drug candidates and hub genes, molecular docking experiments were conducted to determine the interaction between drugs and protein residues, including hydrogen bonding, π-π interactions, hydrophobic interactions, and so forth. The docking scores were then referenced to infer whether the drugs exhibited specific active effects.
In this study, 38 clinically approved anti-cancerous drug candidates were molecularly docked with the target proteins AURKA, BIRC5, CDK1, EZH2, FOXM1, KLF4 and UBE2C. The molecular docking results indicated that drugs had strong binding interactions with the target proteins and a high degree of match (docking score less than −6). The top five compounds with the strongest affinity for each target are listed in Table 4, and the interactions between drugs and the highest binding affinity for each target protein (docking score less than −6) are shown in Figure 8. Based on the comprehensive results, AURKA, CDK1 and EZH2 may be used as therapeutic targets for SCCE, and the docking score of the corresponding drugs are all less than −6. In addition, the crizotinib can form hydrogen bond interactions with the GLU-211, ALA-213 and PRO-214 amino acids of the AURKA protein, and can also form strong hydrophobic interactions with ASN-261, LEU-263, LYS-162, etc., the compound forms van der Waals interactions with the surrounding amino acids. Lapatinib can form hydrogen bond interactions with amino acids GLN-132, SER-84 and LYS-20 of CDK1 protein, and hydrophobic interactions with amino acids such as MET-85, ILE-10, and TYR-15, which enhance the binding ability between the two. Rucaparib is inserted deep into the EZH2 protein pocket, fit well with the protein pocket, and form hydrogen bond interactions with the ARG-685, TYR-726 and LEU-666 amino acid and hydrophobic interactions with the surrounding hydrophobic amino acids (PHE-686, TYR-641, PHE-665, LEU-734, etc.). Similarly, rucaparib compounds are also able to form hydrogen bonding (ILE-113, VAL-96, ALA-93), hydrophobic (LEU-118, PRO-86) and other interactions with UBE2C proteins, and these interactions are effective in inducing the formation of stable complexes between small molecules and proteins. The above hub proteins have been fully studied, and their potential binding sites have been illustrated by researchers through molecular docking and molecular dynamics simulations, which overlap with the binding sites in the results of this study. In summary, the drug candidates are able to form different kinds of interactions such as hydrogen bonding, hydrophobic bonding, and conjugation with protein targets, these interactions can effectively help proteins and compounds to form stable complexes.
TABLE 4
| Target | Drug | Docking score | Target | Drug | Docking score |
|---|---|---|---|---|---|
| AURKA | Crizotinib | −8.025 | FOXM1 | Doxorubicin | −5.681 |
| Cladribine | −7.87 | Cladribine | −5.163 | ||
| Rucaparib | −7.845 | Daunorubicin | −5.044 | ||
| Pazopanib | −7.578 | Hydroxyurea | −5.004 | ||
| Valrubicin | −7.35 | Cytarabine | −4.872 | ||
| BIRC5 | Pralatrexate | −5.726 | KLF4 | Cytarabine | −5.722 |
| Cytarabine | −5.586 | Valrubicin | −5.177 | ||
| Cladribine | −5.097 | Pralatrexate | −5.136 | ||
| Valrubicin | −4.882 | Fluorouracil | −4.959 | ||
| Hydroxyurea | −4.859 | Entrectinib | −4.887 | ||
| CDK1 | Lapatinib | −8.505 | UBE2C | Rucaparib | −6.455 |
| Valrubicin | −7.76 | Fluorouracil | −6.336 | ||
| Pazopanib | −6.996 | Cladribine | −5.773 | ||
| Rucaparib | −6.93 | Pralatrexate | −5.091 | ||
| Daunorubicin | −6.891 | Cytarabine | −4.996 | ||
| EZH2 | Rucaparib | −8.936 | |||
| Lonafarnib | −8.536 | ||||
| Dabrafenib | −7.989 | ||||
| Valrubicin | −7.967 | ||||
| Flutamide | −7.449 |
Docking results of drug candidates with hub genes (top five respectively).
FIGURE 8
4 Discussion
SCCE is a rare neuroendocrine carcinoma originating from the esophagus, studies have shown that it has a low incidence, strong invasiveness, and poor prognosis. Precisely because of its low incidence, there are few related studies, and most of them are case study reports. It is difficult to study the occurrence and development mechanism of SCCE, which leads to a single clinical treatment method and a lack of chemotherapy drugs. Therefore, this study aims to mine the transcriptome sequencing data of SCCE by bioinformatics methods, analyze the hub genes of SCCE, and then find potential therapeutic drugs, so as to establish a theoretical basis for elucidating the molecular mechanism of SCCE, exploring diagnostic and prognostic factors, and finding therapeutic drugs.
Over the past few years, with the continuous development of omics technology, researchers have analyzed the gene expression profiles of SCCE and revealed some of the possible molecular mechanisms. Wang et al. constructed a comprehensive genome map of PESC for the first time (Wang et al., 2018), They conducted whole-exome, SNP chip detection and high-depth targeted sequencing analysis on 55 patients, and found three significant but not yet reported and mutated genes PDE3A, PTPRM, CBLN2. Li et al. conducted genome and transcriptome analysis on 65 SCCE patient tissues and found that SCCE has RB1 destruction mediated by multiple mechanisms (Li et al., 2021a). These sequencing analyses have identified potential therapeutic targets for SCCE from different perspectives. In addition, Okumura et al. examined miRNA expression in SCCE tumor tissues and established a cell line (TYUC-1), they found that eight miRNAs were significantly associated with tumors, and validation in the cell line revealed that downregulate MiR-625 can significantly inhibited the migration of TYUC-1, suggesting that miRNAs may also have a role in the diagnosis of SCCE (Okumura, 2015). Furthermore, studies have shown that PAK1 (Gan et al., 2015), Wnt pathway (Chen et al., 2014), tumor immune microenvironment (Zhao Q. et al., 2020), etc. also play an important role in the development of SCCE.
In this study, data from the GEO database were used to conduct a new analysis of the gene expression results of SCCE tissue samples and normal tissue samples. The screening conditions for DEGs were: |log2FC| ≥ 2, P. Value < 0.01. A total of 1,536 DEGs were finally identified, the up- and downregulated genes were intersected with oncogenes and TSGs, respectively, to obtain 54 upregulated oncogenes and 85 downregulated TSGs. The functional effects and pathways associated with these genes were confirmed by GO and KEGG pathway enrichment analysis. The hub genes (AURKA, BIRC5, CDK1, EZH2, FOXM1, KLF4, UBE2C) were identified by constructing a PPI network and topological analysis methods, and other regulatory networks that interact with hub genes were analyzed through multiple databases. Small cell lung cancer and SCCE are both neuroendocrine tumors, gene expression profiling analysis found that the two have very similar mitotic and DNA repair gene expression patterns (Liu et al., 2021). Feng Wang et al.'s genome sequencing analysis and comparison of SCCE showed that SCCE is very similar to esophageal squamous cell carcinoma, even more similar than SCLC (Wang et al., 2018). Therefore, referring to the relevant studies on small cell lung cancer and esophageal squamous cell carcinoma, we may be able to obtain effective information about SCCE. Among the hub genes, CDK1 encodes protein which is a member of the Ser/Thr protein kinase family, it is the catalytic subunit of a highly conserved protein kinase complex called M phase promoting factor (MPF), which plays a pivotal role in the G2/M phase transition of the eukaryotic cell cycle (Fung and Poon, 2005; Wang et al., 2023). In addition, as early as 2005, Donna E Hansel et al. compared the gene expression profiles of normal esophageal epithelial cell lines and esophageal adenocarcinoma cell lines and found that CDK1 was significantly upregulated in esophageal adenocarcinoma and was a marker of dysplasia and a drug potential targets (Hansel et al., 2005). In SCLC, clinical studies by Kexin Han et al. have shown that CDK1 can not only assist in the diagnosis of SCLC, but is also very important in the differential diagnosis of SCLC and NSCLC (Han et al., 2024). The primary function of EZH2 is to catalyze the methylation of the H3 histone, specifically at residue 27 of the histone, resulting in the formation of H3K27Me3. This process is associated with the repression of transcription of specific genes, including tumor suppressor genes. Additionally, it can form complexes with transcription factors or bind directly to the promoters of target genes, thereby regulating gene transcription (Pasini and Di Croce, 2016; Simon and Lange, 2008). The results of Atsushi Yamada et al. showed that the changes of P53 in esophageal squamous cell carcinoma are related to the abnormal expression of EZH2, and the high expression of EZH2 may promote the progression of esophageal squamous cell carcinoma (Yamada et al., 2011). Xie rui et al.'s study showed that microRNA-30d is a tumor suppressor factor in esophageal squamous cell carcinoma, Luciferase reporter gene assay revealed that EZH2 is a direct target gene of microRNA-30d. MicroRNA-30d can inhibit the movement of cancer cells by targeting and inhibiting EZH2 (Xie et al., 2017). In addition, upregulation of EZH2 expression enhances SCLC progression and radioresistance (Zhai et al., 2020),Eric E Gardner et al. found that in patient-derived xenograft models of SCLC, adding EZH2 inhibitors to standard chemotherapy regimens can prevent the emergence of acquired drug resistance and enhance the chemotherapy sensitivity of SLCL (Gardner et al., 2017). More studies have confirmed that a variety of cancers are associated with EZH2 mutations and expression imbalance, and it is an important target for cancer therapy, with 319 targeted drugs currently available (Liu and Yang, 2023). The protein encoded by KLF4 belongs to the Kruppel family of transcription factors, which play a role in the proliferation and differentiation of epithelial cells. Additionally, this protein is involved in regulating the G1 to S transition of the cell cycle following DNA damage, acting as a mediator for the p53 protein (Wei et al., 2006). KLF4 plays a significant role in tumor development; however, its function in specific cancer types is not uniform. For instance, it acts as a tumor suppressor in gastrointestinal and skin squamous cell carcinoma (Leng et al., 2020), but as a tumor promoter in skin melanoma and breast cancer (Hu et al., 2015; Li et al., 2012). Guo Zhang et al. demonstrated through CHIP experiments that KLF4 binds to the promoter of surviving in esophageal squamous cell carcinoma and inhibits its activity, which can induce apoptosis of cancer cells (Zhang et al., 2009), KLF4 expression is also altered in neuroendocrine lung tumors, and its downregulation may be associated with invasiveness (Naranjo Gómez et al., 2014). Results of this research showed that KLF4 was significantly downregulated in SCCE and act as a tumor suppressor. BIRC5 belongs to the inhibitor of apoptosis (IAP) gene family, which encodes negatively regulated proteins that prevent apoptosis (Unruhe et al., 2016). Bioinformatics analysis has revealed that BIRC5 is highly expressed in many cancers, promotes tumor progression and is strongly associated with poor prognosis (Ye et al., 2022). Interestingly, BIRC5 has been shown to be overexpressed in esophageal adenocarcinoma and is one of the autophagy-related genes, which may serve as a biomarker for the diagnosis and prognosis of esophageal adenocarcinoma (Zhu et al., 2020). In SCLC, BIRC5 can also be used as a new therapeutic target. Yang Yunchu et al. knocked down BIRC5, which significantly inhibited cell growth and migration and induced cell apoptosis (Yunchu et al., 2024). The protein encoded by FOXM1 is a transcriptional activator involved in cell proliferation, it undergoes phosphorylation during the M phase of the cell cycle and regulates the expression of multiple cell cycle-related genes (Laoukili et al., 2005). Research by L Gai et al. found that the protein expression of FOXM1 is significantly upregulated in esophageal squamous cell carcinoma, knockdown the FOXM1 can inhibit the migration and invasion ability of cancer cells, and is significantly positively correlated with lymph node metastasis, clinical stage, and tumor infiltration depth (Gai et al., 2016). Epigenomic analysis by Benjamin Ziman et al. found that FOXM1 plays an important role in the development of esophageal adenocarcinoma and regulates anti-tumor immune response in esophageal adenocarcinoma (Ziman et al., 2024). FOXM1 is also upregulated in SCLC and is significantly associated with poor prognosis. Knockdown of FOXM1 induces apoptosis of SCLC cells and enhances chemotherapy sensitivity, and can be used as a prognostic biomarker (Liang et al., 2021). A large number of studies have shown that FOXM1 can promote the progression and metastasis of various cancers, as well as induce chemotherapy resistance, and has become a target for the treatment of various cancers (Khan et al., 2023). UBE2C is a member of the E2 ubiquitin ligase family, the encoded protein is required to disrupt mitotic cell cycle proteins and cell cycle progression. It has been reported to have pro-carcinogenic effects in esophageal, lung, colorectal, ovarian and liver cancer (Bavi et al., 2011; Pagano, 1997). Previous studies have shown that UBE2C is highly expressed in esophageal squamous cell carcinoma and plays different roles in different stages of esophageal squamous cell carcinoma, mainly affecting the biological function of esophageal squamous cell carcinoma through synergistic effects with CDK1, PTTG1 and SKP2 (Li et al., 2021b; Palumbo et al., 2016), its upstream ECRG4 downregulates the expression of UBE2C in ESCC cells through NF-κB signal transduction, and UBE2C is involved in anti-proliferation and pro-apoptosis functions (Li et al., 2018). High expression of UBE2C and sensitivity of tumor cells after silencing UBE2C were also detected in esophageal adenocarcinoma (Lin et al., 2006). Through bioinformatics and identification of patient tissue samples, it was found that UBE2C expression was also upregulated in SCLC and was identified as a hub gene (Chen et al., 2020). In addition, studies have shown that UBE2C is a transcriptional target of FOXM1, and the two are closely related (Nicolau-Neto et al., 2018). The protein encoded by AURKA is a cell cycle-regulated kinase that promotes the formation and/or stabilization of microtubules at both poles of the spindle during chromosome segregation. It may contribute to the development of tumors by participating in the proliferation of cancer cells, epithelial-mesenchymal transition (EMT), apoptosis and cancer stem cell self-renewal (Du et al., 2021). Previous studies have shown that AURKA can cooperate with TPX2 to regulate the progression of esophageal squamous cell carcinoma through EGFR/PI3K/Akt pathways (Du et al., 2023; Du et al., 2020), it can also regulate the expression and phosphorylation level of JAK2 to promote the activity of STAT3, this mechanism plays an important role in gastric cancer and esophageal cancer (Katsha et al., 2014). In addition, knockdown of AURKA can enhance the ferroptosis effect of esophageal squamous cell carcinoma cells (Mi et al., 2024). AURKA has been shown to be an important factor in the development of SCLC, Knockdown of AURKA can induce SCLC cell cycle arrest and apoptosis, thereby inhibiting the progression of SCLC (Lu et al., 2014). Yixiang Li et al. found that some SCLCs are highly sensitive to AURKA inhibitors, but the persistence is poor, the combination of specific AURKA inhibitors and PD-L1 showed a lasting inhibitory effect (Li et al., 2023). In summary, the seven genes found in our study were not only identified as hub genes in this study, but also played an important role in similar esophageal squamous cell carcinoma, esophageal adenocarcinoma, and small cell lung cancer, which also shows that the conclusions drawn in this study have a certain degree of reliability. However, they must play different roles in different types of cancer, and clinical analysis can be performed to further study whether they can be used as differential diagnostic markers for several similar types of cancer. There is no overlap between the hub genes identified in this study and those found in other bioinformatics studies on SCCE, as the selection of hub genes would differ when using different analysis processes.
Using hub genes and DEGs for prediction in the database, a total of 38 drug candidates were screened, including commonly used cancer chemotherapy drugs such as etoposide, cisplatin, fluorouracil, and cladribine. The 38 small molecule compounds were molecularly docked with seven hub genes, and the results showed that most drugs had good binding affinity with AURKA, CDK1 and EZH2 proteins, and had the potential to become inhibitors of the corresponding targets, thereby playing a role in the treatment of SCCE. Hub genes obtained in this study have been previously demonstrated by a large number of literature to be oncogenes or tumor suppressor genes, and they play an important role in various cancers. However, this is the first time that they have been reported in SCCE, and may be related to certain mechanisms of SCCE pathogenesis. Later, more correlation analysis of clinical pathological characteristics may further reveal more mechanisms of SCCE pathogenesis, and can also verify the value of these genes as SCCE diagnosis or prognosis. The anti-tumor drugs predicted based on these genes were also reported for the first time in SCCE. Clinically, there are very limited drugs used to treat SCCE, and most chemotherapy drugs have a certain degree of drug resistance in patients. Therefore, some new drugs are urgently needed, and completely new drugs take a long time and consume a lot of resources from research and development to final clinical application. Our research provides drugs that have not yet been used in SCCE but have been used in other tumors clinically, studies have shown that they have great potential for treating SCCE. This study further explored the specific sites of action between drugs and targets, providing the possibility of studying these drugs in the direction of targeted therapy, which may have surprising discoveries. Preclinical and clinical studies on their safety and efficacy may enable them to be used in clinical practice as soon as possible, alleviating the current clinical dilemma.
However, the dataset used in this study only has three samples, the sample size is small, and due to the rarity of SCCE, it was difficult to experimentally verify the reliability of hub genes. Therefore, whether these drugs are truly effective for SCCE still requires a large number of in vitro and in vivo experiments to confirm. The present study offers a novel approach to investigating the etiology of rare diseases and identifying prospective pharmaceutical agents.
5 Conclusion
In summary, it can be concluded that AURKA, BIRC5, CDK1, EZH2, FOXM1, KLF4 and UBE2C may be the hub genes closely associated with the development of SCCE. In addition, crizotinib with AURKA, lapatinib with CDK1, rucaparib with EZH2 and UBE2C have good binding affinity, suggesting that these drugs may work against the aforementioned targets. CDK1, EZH2, AURKA proteins have good binding ability with several drugs, further suggesting that targeting these three genes to find drugs may be a powerful approach for the treatment of SCCE. However, it is necessary to validate these conclusions through in vitro and in vivo experiments.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
MZ: Conceptualization, Data curation, Methodology, Validation, Visualization, Writing–original draft. YC: Data curation, Software, Writing–original draft. QY: Methodology, Software, Writing–original draft. JL: Methodology, Validation, Visualization, Writing–review and editing. SC: Project administration, Resources, Software, Supervision, Writing–review and editing. LL: Conceptualization, Funding acquisition, Methodology, Supervision, Validation, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was financially supported by the Science and Technology Project of Nanchong City (No. 23JCYJPT0026), the Medical Research Project of Sichuan Medical Association (No. S23088), the Beijing Medical and Health Economy Research Association Research Project (No. BJHE2023-PRP- 006), and the Sichuan Medical and Health Care Promotion Institute Research Project (No. KY2023SJ0244).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbinf.2025.1495052/full#supplementary-material
SUPPLEMENTARY FIGURE S1Prediction of targeted miRNAs and lncRNAs. (A) Venn diagram of the predicted miRNAs via miRDB, miRwalk and Targetscan database. (B) Venn diagram of predicted miRNAs and OncomiR. (C) Venn diagram of the predicted lncRNAs via LncBase and starbase.
References
1
AmarasenaI. U.ChatterjeeS.WaltersJ. A. E.Wood-BakerR.FongK. M. (2015). Platinum versus non-platinum chemotherapy regimens for small cell lung cancer. Cochrane Database Syst. Rev.2015 (8). 10.1002/14651858.CD006849.pub3
2
BaviP.UddinS.AhmedM.JehanZ.BuR.AbubakerJ.et al (2011). Bortezomib stabilizes mitotic cyclins and prevents cell cycle progression via inhibition of Ube2c in colorectal carcinoma. Am. J. Pathol.178 (5), 2109–2120. 10.1016/j.ajpath.2011.01.034
3
CaiW.GeW.YuanY.DingK.TanY.WuD.et al (2019). A 10-year population-based study of the differences between necs and carcinomas of the esophagus in terms of clinicopathology and survival. J. Cancer10 (6), 1520–1527. 10.7150/jca.29483
4
CannonM.StevensonJ.StahlK.BasuR.CoffmanA.KiwalaS.et al (2024). Dgidb 5.0: rebuilding the drug-gene interaction database for precision medicine and drug discovery platforms. Nucleic Acids Res.52 (D1), D1227–D1235. 10.1093/nar/gkad1040
5
ChenS. B.YangJ. S.YangW. P.WengH. R.LiH.LiuD. T.et al (2011). Treatment and prognosis of limited disease primary small cell carcinoma of esophagus. Dis. Esophagus24 (2), 114–119. 10.1111/j.1442-2050.2010.01112.x
6
ChenW. W.WangF.ZhangD. S.LuoH. Y.WangZ. Q.WangF. H.et al (2014). Primary small cell carcinoma of the esophagus: clinicopathological study of 44 cases. BMC Cancer14, 222. Epub 20140325. 10.1186/1471-2407-14-222
7
ChenX.WangL.SuX.LuoS. Y.TangX.HuangY. (2020). Identification of potential target genes and crucial pathways in small cell lung cancer based on bioinformatic strategy and human samples. PLoS One15 (11), e0242194. Epub 20201113. 10.1371/journal.pone.0242194
8
CicinI.KaragolH.UzunogluS.UygunK.UstaU.KocakZ.et al (2007). Extrapulmonary small‐cell carcinoma compared with small‐cell lung carcinoma. Cancer110 (5), 1068–1076. 10.1002/cncr.22887
9
DingJ.JiJ.ZhuW.ZhouK.HanJ.ZhangY.et al (2013). A retrospective study of different treatments of limited-stage small-cell esophageal carcinoma and associated prognostic factor analysis. Dis. Esophagus26 (7), 696–702. 10.1111/dote.12017
10
DuR.HuangC.ChenH.LiuK.XiangP.YaoN.et al (2020). Sdcbp/Mda-9/Syntenin phosphorylation by aurka promotes esophageal squamous cell carcinoma progression through the egfr-pi3k-akt signaling pathway. Oncogene39 (31), 5405–5419. Epub 20200622. 10.1038/s41388-020-1369-2
11
DuR.HuangC.LiuK.LiX.DongZ. (2021). Targeting aurka in cancer: molecular mechanisms and opportunities for cancer therapy. Mol. Cancer20 (1), 15. Epub 20210115. 10.1186/s12943-020-01305-3
12
DuR.LiK.ZhouZ.HuangY.GuoK.ZhangH.et al (2023). Bioinformatics and experimental validation of an aurka/tpx2 Axis as a potential target in esophageal squamous cell carcinoma. Oncol. Rep.49 (6), 116. Epub 20230421. 10.3892/or.2023.8553
13
FaziR.TintoriC.BraiA.BottaL.SelvarajM.GarbelliA.et al (2015). Homology model-based virtual screening for the identification of human helicase Ddx3 inhibitors. J. Chem. Inf. Model55 (11), 2443–2454. Epub 20151106. 10.1021/acs.jcim.5b00419
14
FriesnerR. A.BanksJ. L.MurphyR. B.HalgrenT. A.KlicicJ. J.MainzD. T.et al (2004). Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J. Med. Chem.47 (7), 1739–1749. 10.1021/jm0306430
15
FungT. K.PoonR. Y. (2005). A roller coaster ride with the mitotic cyclins. Semin. Cell Dev. Biol.16 (3), 335–342. 10.1016/j.semcdb.2005.02.014
16
GaiL.MaoG. X.LiuJ.HuangH.WangX.YaoN. H. (2016). Expression and significance of Foxm1 in esophageal squamous cell carcinoma cells in vitro and in vivo. Zhonghua Zhong Liu Za Zhi38 (3), 179–184. 10.3760/cma.j.issn.0253-3766.2016.03.004
17
GanJ.ZhangY.KeX.TanC.RenH.DongH.et al (2015). Dysregulation of Pak1 is associated with DNA damage and is of prognostic importance in primary esophageal small cell carcinoma. Int. J. Mol. Sci.16 (6), 12035–12050. Epub 20150527. 10.3390/ijms160612035
18
GardnerE. E.LokB. H.SchneebergerV. E.DesmeulesP.MilesL. A.ArnoldP. K.et al (2017). Chemosensitive relapse in small cell lung cancer proceeds through an ezh2-slfn11 Axis. Cancer Cell31 (2), 286–299. 10.1016/j.ccell.2017.01.006
19
HanK.ChenY.SunX.WenL.WuY.ChenS.et al (2024). Combining serum Cdk1 with tumor markers for the diagnosis of small cell lung cancer. Clin. Transl. Oncol. Epub 20241013. 10.1007/s12094-024-03722-y
20
HanselD. E.DharaS.HuangR. C.AshfaqR.DeaselM.ShimadaY.et al (2005). Cdc2/Cdk1 expression in esophageal adenocarcinoma and precursor lesions serves as a diagnostic and cancer progression marker and potential novel drug target. Am. J. Surg. Pathol.29 (3), 390–399. 10.1097/00000478-200503000-00014
21
HuD.GurM.ZhouZ.GamperA.HungM. C.FujitaN.et al (2015). Interplay between arginine methylation and ubiquitylation regulates klf4-mediated genome stability and carcinogenesis. Nat. Commun.6, 8419. Epub 20150930. 10.1038/ncomms9419
22
HuangH. C.ZhengS.VanBurenV.ZhaoZ. (2010). Discovering disease-specific biomarker genes for cancer diagnosis and prognosis. Technol. Cancer Res. Treat.9 (3), 219–229. 10.1177/153303461000900301
23
HudsonE.PowellJ.MukherjeeS.CrosbyT. D. L.BrewsterA. E.MaughanT. S.et al (2007). Small cell oesophageal carcinoma: an institutional experience and review of the literature. Br. J. Cancer96 (5), 708–711. 10.1038/sj.bjc.6603611
24
IshidaH.KasajimaA.KameiT.MiuraT.OkaN.YazdaniS.et al (2017). Sox2 and Rb1 in esophageal small-cell carcinoma: their possible involvement in pathogenesis. Mod. Pathol.30 (5), 660–671. Epub 20170120. 10.1038/modpathol.2016.222
25
ItoT.MasuiT.KomotoI.DoiR.OsamuraR. Y.SakuraiA.et al (2021). Jnets clinical practice guidelines for gastroenteropancreatic neuroendocrine neoplasms: diagnosis, treatment, and follow-up: a synopsis. J. Gastroenterology56 (11), 1033–1044. 10.1007/s00535-021-01827-7
26
KatshaA.ArrasJ.SouttoM.BelkhiriA.El-RifaiW. (2014). Aurka regulates jak2-stat3 activity in human gastric and esophageal cancers. Mol. Oncol.8 (8), 1419–1428. Epub 20140603. 10.1016/j.molonc.2014.05.012
27
KhanM. A.KhanP.AhmadA.FatimaM.NasserM. W. (2023). Foxm1: a small fox that makes more tracks for cancer progression and metastasis. Semin. Cancer Biol.92, 1–15. Epub 20230322. 10.1016/j.semcancer.2023.03.007
28
KuG. Y.MinskyB. D.RuschV. W.BainsM.KelsenD. P.IlsonD. H. (2008). Small-cell carcinoma of the esophagus and gastroesophageal junction: review of the memorial sloan-kettering experience. Ann. Oncol.19 (3), 533–537. 10.1093/annonc/mdm476
29
LambJ. (2007). The connectivity map: a new tool for biomedical research. Nat. Rev. Cancer7 (1), 54–60. 10.1038/nrc2044
30
LaoukiliJ.KooistraM. R.BrásA.KauwJ.KerkhovenR. M.MorrisonA.et al (2005). Foxm1 is required for execution of the mitotic programme and chromosome stability. Nat. Cell Biol.7 (2), 126–136. Epub 20050116. 10.1038/ncb1217
31
LawS. Y. K.FokM.LamK. Y.LokeS. L.MaL. T.WongJ. (1994). Small cell carcinoma of the esophagus. Cancer73 (12), 2894–2899. 10.1002/1097-0142(19940615)73:12<2894::Aid-cncr2820731204>3.0.Co;2-m
32
LengZ.LiY.ZhouG.LvX.AiW.LiJ.et al (2020). Krüppel‐like factor 4 regulates stemness and mesenchymal properties of colorectal cancer stem cells through the TGF‐β1/Smad/snail pathway. J. Cell Mol. Med.24 (2), 1866–1877. Epub 20191212. 10.1111/jcmm.14882
33
LiA. F. Y.HsuH. S.HsuC. Y.LiA. C. H.LiW. Y.LiangW. Y.et al (2010b). A 20‐year retrospective study of small‐cell carcinomas in taiwan. J. Surg. Oncol.102 (5), 497–502. 10.1002/jso.21629
34
LiA. F. Y.LiA. C. H.HsuC. Y.LiW. Y.HsuH. S.ChenJ. Y. (2010a). Small cell carcinomas in gastrointestinal tract: immunohistochemical and clinicopathological features. J. Clin. Pathology63 (7), 620–625. 10.1136/jcp.2010.077024
35
LiJ.ZhengH.YuF.YuT.LiuC.HuangS.et al (2012). Deficiency of the kruppel-like factor Klf4 correlates with increased cell proliferation and enhanced skin tumorigenesis. Carcinogenesis33 (6), 1239–1246. 10.1093/carcin/bgs143
36
LiL.LiX.WangW.GaoT.ShiZ. (2018). Ube2c is involved in the functions of Ecrg4 on esophageal squamous cell carcinoma. Biomed. Pharmacother.98, 201–206. Epub 20171227. 10.1016/j.biopha.2017.12.066
37
LiR.PangX. F.HuangZ. G.YangL. H.PengZ. G.MaJ.et al (2021b). Overexpression of Ube2c in esophageal squamous cell carcinoma tissues and molecular analysis. BMC Cancer21 (1), 996. Epub 20210906. 10.1186/s12885-021-08634-6
38
LiR.YangZ.ShaoF.ChengH.WenY.SunS.et al (2021a). Multi-omics profiling of primary small cell carcinoma of the esophagus reveals Rb1 disruption and additional molecular subtypes. Nat. Commun.12 (1), 3785. Epub 20210618. 10.1038/s41467-021-24043-6
39
LiY.MahadevanN. R.DuplaquetL.HongD.DurmazY. T.JonesK. L.et al (2023). Aurora a kinase inhibition induces accumulation of sclc tumor cells in mitosis with restored interferon signaling to increase response to Pd-L1. Cell Rep. Med.4 (11), 101282. 10.1016/j.xcrm.2023.101282
40
LiangS. K.HsuC. C.SongH. L.HuangY. C.KuoC. W.YaoX.et al (2021). Foxm1 is required for small cell lung cancer tumorigenesis and associated with poor clinical prognosis. Oncogene40 (30), 4847–4858. Epub 20210621. 10.1038/s41388-021-01895-2
41
LinJ.RaoofD. A.WangZ.LinM. Y.ThomasD. G.GreensonJ. K.et al (2006). Expression and effect of inhibition of the ubiquitin-conjugating enzyme E2c on esophageal adenocarcinoma. Neoplasia8 (12), 1062–1071. 10.1593/neo.05832
42
LiuD.WenJ.ChenJ.WangB.XuX.ZhangZ.et al (2021). A comparative analysis of the gene expression profiles of small cell esophageal carcinoma, small cell lung cancer, and esophageal adeno/squamous carcinoma. Front. Surg.8, 655159. Epub 20210730. 10.3389/fsurg.2021.655159
43
LiuD.XuX.WenJ.XieL.ZhangJ.ShenY.et al (2018). Integrated genome-wide analysis of gene expression and DNA copy number variations highlights stem cell-related pathways in small cell esophageal carcinoma. Stem Cells Int.2018, 1–8. Epub 20180830. 10.1155/2018/3481783
44
LiuY.SunJ.ZhaoM. (2017). Ongene: a literature-based database for human oncogenes. J. Genet. Genomics44 (2), 119–121. Epub 20161226. 10.1016/j.jgg.2016.12.004
45
LiuY.YangQ. (2023). The roles of Ezh2 in cancer and its inhibitors. Med. Oncol.40 (6), 167. Epub 20230506. 10.1007/s12032-023-02025-6
46
LuY.LiuY.JiangJ.XiZ.ZhongN.ShiS.et al (2014). Knocking down the expression of aurora-a gene inhibits cell proliferation and induces G2/M phase arrest in human small cell lung cancer cells. Oncol. Rep.32 (1), 243–249. Epub 20140516. 10.3892/or.2014.3194
47
LyuF.HanF.GeC.MaoW.ChenL.HuH.et al (2023). Omicstudio: a composable bioinformatics cloud platform with real-time feedback that can generate high-quality graphs for publication. Imeta2 (1), e85. Epub 20230206. 10.1002/imt2.85
48
MiY.ChenL.WangC.MiaoY.SongC.SuJ.et al (2024). Aurka knockdown inhibits esophageal squamous cell carcinoma progression through ferroptosis. Heliyon10 (7), e28365. Epub 20240319. 10.1016/j.heliyon.2024.e28365
49
Naranjo GómezJ. M.BernalJ. F.ArranzP. G.FernándezS. L.RomanJ. J. (2014). Alterations in the expression of P53, Klf4, and P21 in neuroendocrine lung tumors. Arch. Pathol. Lab. Med.138 (7), 936–942. 10.5858/arpa.2013-0119-OA
50
Nicolau-NetoP.PalumboA.De MartinoM.EspositoF.de Almeida SimãoT.FuscoA.et al (2018). Ube2c is a transcriptional target of the cell cycle regulator Foxm1. Genes (Basel)9 (4), 188. Epub 20180329. 10.3390/genes9040188
51
OkumuraT.ShimadaY.OmuraT.HiranoK.NagataT.TsukadaK. (2015). Microrna profiles to predict postoperative prognosis in patients with small cell carcinoma of the esophagus. Anticancer Res.35 (2), 719–727.
52
PaganoM. (1997). Cell cycle regulation by the ubiquitin pathway. Faseb J.11 (13), 1067–1075. 10.1096/fasebj.11.13.9367342
53
PalumboA.Jr.Da CostaN. M.De MartinoM.SepeR.PellecchiaS.de SousaV. P.et al (2016). Ube2c is overexpressed in escc tissues and its abrogation attenuates the malignant phenotype of escc cell lines. Oncotarget7 (40), 65876–65887. 10.18632/oncotarget.11674
54
PasiniD.Di CroceL. (2016). Emerging roles for polycomb proteins in cancer. Curr. Opin. Genet. Dev.36, 50–58. Epub 20160502. 10.1016/j.gde.2016.03.013
55
RajadnyaR.SharmaN.MahajanA.UlheA.PatilR.HegdeM.et al (2024). Novel systems biology experimental pipeline reveals matairesinol's antimetastatic potential in prostate cancer: an integrated approach of network pharmacology, bioinformatics, and experimental validation. Brief. Bioinform25 (5), bbae466. 10.1093/bib/bbae466
56
RajeswariM.SanthiN.BhuvaneswariV. (2014). Pharmacophore and virtual screening of Jak3 inhibitors. Bioinformation10 (3), 157–163. Epub 20140319. 10.6026/97320630010157
57
SimonJ. A.LangeC. A. (2008). Roles of the Ezh2 histone methyltransferase in cancer epigenetics. Mutat. Res.647 (1-2), 21–29. Epub 20080803. 10.1016/j.mrfmmm.2008.07.010
58
SongY.WangL.-H.HeJ.WangJ.-W. (2009). Treatment and prognosis of primary esophageal small cell carcinoma. Ai Zheng28 (3), 303–307.
59
TangD.ChenM.HuangX.ZhangG.ZengL.ZhangG.et al (2023). Srplot: a free online platform for data visualization and graphing. PLoS One18 (11), e0294236. Epub 20231109. 10.1371/journal.pone.0294236
60
TangZ.LiC.KangB.GaoG.LiC.ZhangZ. (2017). Gepia: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res.45 (W1), W98–w102. 10.1093/nar/gkx247
61
UnruheB.SchröderE.WünschD.KnauerS. K. (2016). An old flame never dies: survivin in cancer and cellular senescence. Gerontology62 (2), 173–181. Epub 20150703. 10.1159/000432398
62
VosB.RozemaT.MillerR. C.HendliszA.Van LaethemJ. L.KhanfirK.et al (2011). Small cell carcinoma of the esophagus: a multicentre rare cancer network study. Dis. Esophagus24 (4), 258–264. Epub 20101112. 10.1111/j.1442-2050.2010.01133.x
63
WangF.LiuD. B.ZhaoQ.ChenG.LiuX. M.WangY. N.et al (2018). The genomic landscape of small cell carcinoma of the esophagus. Cell Res.28 (7), 771–774. Epub 20180504. 10.1038/s41422-018-0039-1
64
WangQ.BodeA. M.ZhangT. (2023). Targeting Cdk1 in cancer: mechanisms and implications. NPJ Precis. Oncol.7 (1), 58. Epub 20230613. 10.1038/s41698-023-00407-7
65
WeiD.KanaiM.HuangS.XieK. (2006). Emerging role of Klf4 in human gastrointestinal cancer. Carcinogenesis27 (1), 23–31. Epub 20051011. 10.1093/carcin/bgi243
66
XiaoQ.XiaoH.OuyangS.TangJ.ZhangB.WangH. (2019). Primary small cell carcinoma of the esophagus: comparison between a Chinese cohort and surveillance, epidemiology, and end results (seer) data. Cancer Med.8 (3), 1074–1085. 10.1002/cam4.2001
67
XieM. R.XuS. B.SunX. H.KeL.MeiX. Y.LiuC. Q.et al (2015). Role of surgery in the management and prognosis of limited-stage small cell carcinoma of the esophagus. Dis. Esophagus28 (5), 476–482. 10.1111/dote.12230
68
XieR.WuS. N.GaoC. C.YangX. Z.WangH. G.ZhangJ. L.et al (2017). Microrna-30d inhibits the migration and invasion of human esophageal squamous cell carcinoma cells via the post-transcriptional regulation of enhancer of zeste homolog 2. Oncol. Rep.37 (3), 1682–1690. Epub 20170125. 10.3892/or.2017.5405
69
YamadaA.FujiiS.DaikoH.NishimuraM.ChibaT.OchiaiA. (2011). Aberrant expression of Ezh2 is associated with a poor outcome and P53 alteration in squamous cell carcinoma of the esophagus. Int. J. Oncol.38 (2), 345–353. Epub 20101209. 10.3892/ijo.2010.868
70
YeH. B.MaB. J.MengG. Q.TaoS.WangY.ChenZ.et al (2022). Bioinformatics analysis of Birc5 in human cancers. Ann. Transl. Med.10 (16), 888. 10.21037/atm-22-3496
71
YunchuY.MiyanagaA.MatsudaK.KamioK.SeikeM. (2024). Exploring effective biomarkers and potential immune related gene in small cell lung cancer. Sci. Rep.14 (1), 7604. Epub 20240331. 10.1038/s41598-024-58454-4
72
ZhaiG.LiJ.ZhengJ.AnP.ChenX.WangX.et al (2020). Htert promoter methylation promotes small cell lung cancer progression and radiotherapy resistance. J. Radiat. Res.61 (5), 674–683. 10.1093/jrr/rraa052
73
ZhangG.ZhuH.WangY.YangS.LiuM.ZhangW.et al (2009). Krüppel-like factor 4 represses transcription of the survivin gene in esophageal cancer cell lines. Biol. Chem.390 (5-6), 463–469. 10.1515/bc.2009.060
74
ZhangH. J.ZhouS. L.ZhaoL. F.ShiM. L.JianL. (2019). Clinical analysis of 22 cases with esophageal primary small cell neuroendocrine carcinoma. Chin. J. cancer Prev. Treat.26 (16), 1200–1207.
75
ZhangZ.XiaoH.XieF.ZhangH.ChenC.XiaoH.et al (2014). High-incidence of pten mutations in Chinese patients with primary small cell carcinoma of the esophagus. BMC Cancer14, 19. Epub 20140114. 10.1186/1471-2407-14-19
76
ZhaoK.MaL.FengL.HuangZ.MengX.YuJ. (2020a). Cd155 overexpression correlates with poor prognosis in primary small cell carcinoma of the esophagus. Front. Mol. Biosci.7, 608404. Epub 20210107. 10.3389/fmolb.2020.608404
77
ZhaoM.KimP.MitraR.ZhaoJ.ZhaoZ. (2016). Tsgene 2.0: an updated literature-based knowledgebase for tumor suppressor genes. Nucleic Acids Res.44 (D1), D1023–D1031. Epub 20151120. 10.1093/nar/gkv1268
78
ZhaoQ.ChenY. X.WuQ. N.ZhangC.LiuM.WangY. N.et al (2020b). Systematic analysis of the transcriptome in small-cell carcinoma of the oesophagus reveals its immune microenvironment. Clin. Transl. Immunol.9 (10), e1173. Epub 20201005. 10.1002/cti2.1173
79
ZhouG.SoufanO.EwaldJ.HancockR. E. W.BasuN.XiaJ. (2019). Networkanalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res.47 (W1), W234–W241. 10.1093/nar/gkz240
80
ZhuL.DongL.FengM.YangF.JiangW.HuangZ.et al (2020). Profiles of autophagy-related genes in esophageal adenocarcinoma. BMC Cancer20 (1), 943. Epub 20201001. 10.1186/s12885-020-07416-w
81
ZhuY.QiuB.LiuH.LiQ.XiaoW.HuY.et al (2014). Primary small cell carcinoma of the esophagus: review of 64 cases from a single institution. Dis. Esophagus27 (2), 152–158. 10.1111/dote.12069
82
ZimanB.YangQ.ZhengY.ShethM.NamC.ZhaoH.et al (2024). Epigenomic analyses identify Foxm1 as a key regulator of anti-tumor immune response in esophageal adenocarcinoma. Cell Death Dis.15 (2), 152. Epub 20240219. 10.1038/s41419-024-06488-x
Summary
Keywords
bioinformatics, small cell carcinoma of the esophagus, CDK1, AURKA, DGIdb, ConnectivityMAP, anti-tumor drug
Citation
Zhu M, Chu Y, Yuan Q, Li J, Chen S and Li L (2025) Integrated bioinformatics analysis to explore potential therapeutic targets and drugs for small cell carcinoma of the esophagus. Front. Bioinform. 5:1495052. doi: 10.3389/fbinf.2025.1495052
Received
14 September 2024
Accepted
03 January 2025
Published
28 January 2025
Volume
5 - 2025
Edited by
Minh Nguyen, Bioinformatics Institute (A∗STAR), Singapore
Reviewed by
Rajesh Kumar, National Institutes of Health (NIH), United States
Bhavika Mam, Independent Researcher, Palo Alto, CA, United States
Updates
Copyright
© 2025 Zhu, Chu, Yuan, Li, Chen and Li.
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: Lin Li, linli82@cqu.edu.cn
Disclaimer
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.