Identification of Key Biomarkers and Pathways in Small-Cell Lung Cancer Using Biological Analysis

Liu, Huanqing; Li, Tingting; Ye, Xunda; Lyu, Jun

doi:https://doi.org/10.1155/2021/5953386

BioMed Research International

On this page

Abstract Background Materials and Methods Results Discussion Conclusion Data Availability Consent Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 5953386 | https://doi.org/10.1155/2021/5953386

Identification of Key Biomarkers and Pathways in Small-Cell Lung Cancer Using Biological Analysis

Huanqing Liu,¹Tingting Li,²Xunda Ye,³and Jun Lyu¹

Academic Editor: Pei Jiang

Received09 Jun 2021

Accepted25 Sept 2021

Published19 Oct 2021

Abstract

Background. Small-cell lung cancer (SCLC) is a major cause of carcinoma-related deaths worldwide. The aim of this study was to identify the key biomarkers and pathways in SCLC using biological analysis. Methods. Key genes involved in the development of SCLC were identified by downloading three datasets from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were screened using the GEO2R online analyzer; for the functional annotation and pathway enrichment analysis of genes, Funrich software was used. Construction of protein-to-protein interaction (PPI) networks was accomplished using the Search Tool for the Retrieval of Interacting Genes (STRING), and network visualization and module identification were performed using Cytoscape. Results. A total of 268 DEGs were ultimately obtained. The enriched functions and pathways of the upregulated DEGs included cell cycle, mitotic, and DNA replication, and the downregulated DEGs were enriched in epithelial-to-mesenchymal transition, serotonin degradation, and noradrenaline. Analysis of significant modules demonstrated that the upregulated genes are primarily concentrated in functions related to cell cycle and DNA replication. Kaplan-Meier analysis of hub genes revealed that they may promote the carcinogenesis and progression of SCLC. The result of ONCOMINE demonstrated that these 10 hub genes were significantly overexpressed in SCLC compared with normal samples. Conclusion. Identification of the molecular functions and signaling pathways of participating DEGs can deepen the current understanding of the molecular mechanisms of SCLC. The knowledge gained from this work may contribute to the development of treatment options and improve the prognosis of SCLC in the future.

1. Background

Lung cancer is a major cause of carcinoma-related deaths worldwide, and approximately 2.21 million new cases of this disease are estimated in 2020 according to the latest WHO data on lung cancer. Gene mutations and cell environment changes may affect the formation, growth, and metastasis of tumors [1]. Small-cell lung cancer (SCLC) is the main histological form of pulmonary carcinoma. Conventional treatment methods include chemotherapy, radiotherapy, and surgery. Chemotherapy is the most important treatment method for SCLC, but issues such as high drug resistance and recurrence rates limit its effectiveness. Because most patients are usually at advanced stages of the disease at the time of diagnosis, SCLC is often characterized by low survival rates and poor quality of life. In fact, SCLC has a 5-year survival rate of <6% and high mortality; moreover, it is highly invasive and prone to early hematogenesis and lymphatic metastasis [2]. To date, no molecular targeted drugs have yet been shown to significantly prolong patient survival [3]. More importantly, the molecular mechanisms underlying the occurrence, development, invasion, and metastasis of SCLC remain poorly understood. Thus, finding methods to obtain a prognosis of SCLC and identify potential biomarkers for targeted therapy is of great importance to improve the clinical efficacy of lung cancer.

Advances in gene expression profile chip technology over the last few decades have established a strong foundation for the overall exploration of differentially expressed genes (DEGs) in lung cancer and their biological functions. Studies have shown that the expression profiles of cancerous tissues differ from those of neighboring noncancerous tissues [4]. Given theoretical facts on cancer and medical database as well as data mining technology in big data era, we rationally speculated that DEGs may influence the occurrence and development of many diseases, including malignancies. Although some RNA-sequencing (RNA-seq) are insensitive to ribonuclease because of their unique structure, different genes can exist in tissues and serum and possibly even function as biomarkers for cancer. Microarray and bioinformatics technologies show broad applications in disease research, especially for identifying DEGs, mRNA, and miRNA, as well as their elucidation of molecular mechanisms [5–7]. Udhaya et al. used bioinformatics methods to identify four DEGs and associated pathways in systemic lupus erythematosus [8]. Many non-small-cell lung cancer (NSCLC) gene expression profile studies have been conducted using microarray technology, and numbers of NSCLC-related DEGs have been identified [9]. Mao et al. used bioinformatics methods to compare the differential expressions of mRNA and microRNA between SCLC tissues and normal lung tissues to explore the pathogenesis and potential molecular markers of SCLC [10]. Due to the complexity of the biological characteristics of SCLC, the key biomarkers and specific targets for prognosis of SCLC remain unclear. Therefore, it is necessary to explore more genetic information and screen out potential or promising biomarkers for prognosis of SCLC. Chen et al. identified 8 DEGs served as new biomarkers for prognosis of SCLC [11]. However, due to the lack of SCLC-related gene chip data, we need to integrate more chip database for research and analysis of this disease. In the current study, three gene expression profiles, i.e., GSE30219, GSE99316, and GSE149507, were downloaded from public databases to identify DEGs in SCLC and their related pathways. Moreover, Kaplan-Meier analysis was used to explore the relationship between prognosis and hub gene expression level, and the Oncomine database was used to explore hub genes' expressions between SCLC tissues and normal tissues. The obtained data indicate that the identified DEGs may be used as key biomarkers of SCLC. The related pathways also offer insights into the pathogenesis of the disease.

2. Materials and Methods

2.1. Microarray Data

Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/) was adopted to perform gene chip screening. The target chip access criteria are as follows: (1) clinical SCLC specimens of the patients were excluded from cell lines and animal experiments. (2) The selected chip should contain SCLC and normal tissue samples. (3) Only mRNA chips that have been standardized were employed. According to the screening criteria, several datasets were screened out. We selected the datasets with large numbers of samples as the research object. Three gene profiles, i.e., GSE30219, GSE99316, and GSE149507, were then obtained from the GEO database. The GSE30219 dataset contained 21 SCLC and 14 normal tissue samples, the GSE99316 dataset consisted of 23 SCLC and 42 normal tissue samples, while the GSE149507 dataset consisted of 18 SCLC and 18 normal tissue samples.

2.2. Identifications of DEGs

We used the GEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r) statistical tool to recalculate and evaluate genes expressed differently in human SCLC tissues and adjacent noncancerous lung tissues.

Benjamini and Hochberg (error detection rate) and -test methods were used in conjunction with GEO2R to calculate FDR and values, respectively, to identify DEGs. The DEGs were confirmed by the following criteria: and . If , the gene expression is considered upregulated; if , the gene is considered downregulated.

2.3. Functional Enrichment Analysis of DEGs

Funrich (http://www.funrich.org/), a stand-alone software tool used mainly for functional enrichment and interaction network analysis of genes and proteins, can help users load a customized database against which functional enrichment analysis can be carried out. In our study, molecular function, biological process, cell composition, and biological pathway was applied to analyze the DEGs.

2.4. Construction of PPI Networks and Module Analysis

Search Tool for the Retrieval of Interacting Genes version 11.0 (STRING; http://string-db.org) is an online tool that could be used to identify the interactions of proteins and obtain insights into the mechanisms of certain diseases. In the present study, an overall was set as the cut-off point. Cytoscape version 3.7.1 software was used to visualize the Protein–proteininteractions (PPI) networks. Molecular Complex Detection version 1.6.1 (MCODE), a Cytoscape plug-in, can confirm areas with dense connections to select a statistically significant model. Key modules within the PPI networks were identified using MCODE (MCODE , degree cut-, node score cut-, max , ). Enrichment analysis of the DEGs in this module was subsequently conducted using Funrich.

2.5. Selection and Analysis of Hub Genes

The top 10 genes in significant modules were selected as hub genes. The cBioPortal online platform was employed to analyze the gene networks obtained, as well as the relationships among coexpressed genes. Biological Networks Gene Oncology version 3.0.3, another Cytoscape plug-in, was used to assess the performance and visualize the results of the bioprocess analysis for the hub genes. The UCSC Cancer Genomics Browser (http://genome-cancer.ucsc.edu) was also utilized to achieve the hierarchical clustering of these genes. We used Kaplan-Meier analysis to explore the relationship between prognosis and hub gene expression level and draw the corresponding survival curves. The ONCOMINE database was used to analyze the expression of the hub gene between tumor and normal tissues in clinical SCLC and the specific expression at each stage by histogram.

2.6. Statistical Analysis

values and FDR were calculated using GEO2R’s built-in -test, Benjamini and Hochberg (false detection rate), and other methods to determine the DEGs of SCLC patients and controls.

3. Results

3.1. Identification of DEGs

According to the inclusion criteria, a total of two mRNA microarray datasets that met the requirements were screened out, namely, GSE30219, GSE99316, and GSE149507. A total of 2659 DEGs (1534 in GSE30219, 1639 in GSE99316, and 614 in GSE149507) (Table S1) were identified from the three datasets. GEO2R analysis showed that 268 DEGs were expressed in the three datasets; Venn diagram analysis of these genes is shown in Figure 1, including 192 upregulated and 76 downregulated genes in both SCLC and normal tissues (Table S2); the list of DEGs were shown in Table 1.

(a)

(b)

(c)

(d)

3.2. Functional Enrichment Analyses of DEGs

The results of GO (including three categories: biological processes (BP), cellular components (CC), and molecular functions (MF)) and biological pathway analyses of these DEGs are shown in Figure 2 (Table S3, S4). The upregulated DEGs related to BP included those for chromosome segregation, regulation of nucleobase, nucleoside, nucleotide and nucleicacid metabolism, cell cycle, and spindle assembly. Most of the upregulated DEGs obtained were related to MF, including motor activity and protein binding. The upregulated DEGs related to CC were mainly located in kinetochore, nucleus, microtubuls, nucleoplasm, chromosome, centrometric region, and condensed chromosome kinetochore.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 2

Functional annotation and path enrichment analysis. (a) BP enriched upregulated DEGs. (b) BP enriched downregulated DEGs. (c) MF enriched upregulated DEGs. (d) MF enriched downregulated DEGs. (e) CC enriched upregulated DEGs. (f) CC enriched downregulated DEGs. The axis is a detailed term for functional annotation and path enrichment. The axis represents the percentage of genes, . The red line is , which is the reference value of statistical analysis truncation standard. The yellow line represents .

While the analysis of BP in downregulated DEGs were most but not significantly concentrated in immune response, cell communication, and cellular defense response, most of the downregulated DEGs obtained were related to MF, including extracellular matrix structual constituent, water channel activity, and catalytic activity (there was no statistical significance). The downregulated DEGs related to CC were mainly located in extracellular and cell surface.

Furthermore, biological pathway analysis was shown in Figure 3; the upregulated DEGs were mainly enriched in cell cycle, mitotic, DNA replication, mitotic M-M/G1 phases, M phase, mitotic prometaphase, and mitotic G1-G1/S phases. The downregulated DEGs were mainly enriched in epithelial-to-mesenchymal transition, serotonin degradation, noradrenaline and adrenaline degradation, FOXA transcription factor networks, and FOXA1 transcription factor network.

(a)

(b)

3.3. Construction of PPI Networks and Significant Module Analysis

Understanding interactions between DEGs that may be related to the development of SCLC is necessary to explore the underlying mechanisms in SCLC. PPI network was detected by STRING, as shown in Figure 4. A total of 224 nodes and 5302 edges, with each node representing a protein (gene) and each edge representing an interaction relationship were obtained. Genes with the most significant modules were screened by MCODE with (Figure 4), and results indicated that the identified genes were all upregulated genes (Figure 5). The results of biological pathway analysis of significant module DEGs were shown in Figure 6 (Table S5); results showed that DEGs in significant module were mainly enriched in cell cycle, mitotic, and cell replication.

3.4. Hub Gene Selection and Analysis

Identification of 10 hub genes () was conducted, and the results are shown in Table 2. The hierarchical clustering results in Figure 7 reveal that the identified hub genes could distinguish SCLC samples from normal ones.

3.5. Prognostic Analysis of Hub Genes

Kaplan-Meier analysis was conducted to analyze the relation between overall survival and the 10 hub genes identified earlier and predict the association of these genes with lung cancer prognosis. The median survival time of the group showing high expression of CHEK1, DTL, KIF14, MCM4, CENPU, NEK2 CDC20, KIF4A, and NCAPG2 was significantly shorter than that of the group demonstrating low expression of these genes, and the difference between groups was statistically significant () as shown in Figure 8.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

3.6. Hub Gene Analysis

The ONCOMINE database could be utilized to determine the expression of hub genes in normal and SCLC tissues. The results in Figure 9 show that the hub genes (except KIF4A) were significantly overexpressed in the SCLC tissues of different studies.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

4. Discussion

SCLC, a subtype of lung cancer, is the sixth most common cause of all cancer-related deaths worldwide and has significant clinical and pathological features, including early metastasis and poor prognosis. The survival time of patients in advanced stages of the disease is less than 1 year, and the 2-year survival rate is only approximately 5% [12]. The high drug resistance and recurrence rates of SCLC are mainly attributed to the high mutation rate of genes involved in this malignancy and genomic instability. Studies have shown that the mutation frequencies of P53 and RB1 in SCLC are 85% and 57%, respectively. Thus, mutations of P53 and RB1 indicate poor prognosis [13]. The molecular mechanisms of SCLC remain unclear, and identification of potential key genes that can serve as biomarkers is an urgent undertaking. Bioinformatics may be used to explore gene-level changes in SCLC and identify potential biomarkers.

In our study, three datasets, GSE30219, GSE99316, and GSE149507, were applied to screen for DEGs between SCLC and normal adjacent tissues, and a total of 268 DEGs were obtained. GO and KEGG enrichment analyses were performed on these 268 DEGs, and results indicated that the upregulated DEGs were mainly enriched in cell cycle, mitotic, DNA replication, mitotic M-M/G1 phases, M phase, mitotic prometaphase, and mitotic G1-G1/S phases and the downregulated DEGs were mainly enriched in epithelial-to-mesenchymal transition, serotonin degradation, noradrenaline and adrenaline degradation, FOXA transcription factor networks, and FOXA1 transcription factor network. Mutations in tumor cells often result in changes in the cell cycle leading to unrestricted growth compared with that of normal cells. Zhang et al. found that cell cycle inhibitors could be used in SCLC to interfere with the cell cycle, induce DNA replication stress and genomic instability, and trigger immune response signal [14, 15]. These studies showed consistency with the finding that dysregulation of the cell cycle promotes tumorigenesis and progression.

In the analysis of significant module genes, it was found that all 112 genes in cluster 1 were upregulated, and we selected the top 10 genes as hub genes. Genes with a were selected as hub genes. CDC20, CENPU, CHEK1, DTL, KIF4A, KIF14, MCM4, NCAPG2, NEK2, and FOXM1 were all located in core nodes in the PPI network, which means these 10 genes may be critical therapeutic targets for SCLC.

CDC20, a homolog of Saccharomyces cerevisiae cellular division cycle 20 protein, serves as an activator of the anaphase-promoting complex, which performs an essential function in governing cell cycle progression for cell division [16]. CDC20 is highly overexpressed in NSCLC patients [17], and downregulation of CDC20 expression can slow down the growth and colony formation rate of lung cancer cells [18]. According to GO analysis results, CDC20 is mainly involved in the biological process of cell cycle; there is evidence that CDC20 regulates the cell cycle progression of cell division by targeting several key degradation substrates [18]. Therefore, our results are consistent with previous studies, suggesting that CDC20 may be a key factor in the onset and progression of SCLC.

Centromeric protein U (CENPU), also known as myelodysplastic/myeloid leukemia factor 1 Interaction protein (MLF1IP) [19], is an important component of spindle recovery after injury. CENPU is reported to be abnormally high expressed in various human tumor tissues and is involved in tumor progression, such as prostate cancer, breast cancer, bladder cancer, and ovarian cancer, and its overexpression has been shown to predict poor prognosis [20–22]. A study has demonstrated that CENPU regulated the proliferation and migration of lung adenocarcinoma cells through the PI3K/AKT pathway [23]; however, there are no data on the carcinogenic effect and clinical significance of CENPU in SCLC. Unfortunately, the biological process related to CENPU was not screened in our study (the result is unknown). In pathway analysis, CENPU was found to be involved in the PLK pathway in addition to the pathways related to cell cycle, DNA replication, and mitosis, which has not been verified by basic experiments.

CHEK1, an evolutionarily conserved serine/threonine kinase, has been shown to regulate cell cycle checkpoints, coordinating cellular activity involved in DNA repair, and cell cycle arrest [24]. A growing number of studies have found that CHEK1 is highly expressed in multiple cancer species [25–27] and is considered a potential target for cancer treatment. At present, there is little evidence on the relationship between CHEK1 and SCLC. Gali-Muhtasib et al. [27] found that CHEK1 was significantly overexpressed in SCLC compared with NSCLC samples and inhibited CHEK1 or ATR could induce genotoxic stress and apoptosis. Through our analysis, it was found that the biological process CHEK1 was mainly involved in was the cell cycle, and its related pathways included ATR and ATM signaling pathways in addition to the cell cycle [28, 29], of which is consistent with the reported CHEK signaling pathway. It is worth noting that CHEK1-mediated inactivation of Cyclin B is implicated in our analysis, which requires further verification.

Until now, DTL was thought to be involved in the regulation of cell cycle and DNA replication to influence tumor progression [30]. Notably, previous studies have shown that DTL expression is elevated in many cancers [31, 32], and the abnormal expression of DTL is also associated with poor prognosis [33]. However, systematic studies of the function of DTL in tumors remain to be evaluated.

KIF4A is one of the members of the kinesin superfamily. KIF4A is involved in a variety of cellular activities and plays a critical role in biological processes such as mitotic spindle formation and DNA damage repair [34]. Our analysis shows that KIF4A is involved in biological processes of regulation of nucleobase, nucleoside, nucleotide, and nucleic acid metabolism, which is consistent with the previous study [35]. A large number of studies have confirmed that KIF4A is overexpressed in colorectal cancer, liver cancer, and lung cancer [36–38] and is an independent prognostic risk factor [38]. It can also be seen from the analysis results of the GEO database that KIF4A is highly expressed in SCLC tissues, and its high expression is related to poor prognosis of patients, which proves that KIF4A can be used as a potential target for the treatment of SCLC.

Like KIF4A, KIF14 is a member of the driver superfamily. KIF14 is widely believed to play a role in tumorigenesis. The overexpression of KIF14 may lead to rapid and error-prone mitosis [39] and is involved in the progression of a variety of malignancies, such as retinoblastoma and gastric cancer [40, 41], providing evidence that KIF14 may be an oncogene in the progression of a variety of cancers, while in the study of lung adenocarcinoma, it was found that KIF14 was underexpressed in 30% of cancer tissue samples [42], and the decreased expression of KIF14 was significantly correlated with the overall survival rate of lung cancer patients [43]. However, it can be seen from our analysis that KIF14 is highly expressed in SCLC cancer tissues, reminding us that we can distinguish lung adenocarcinoma from SCLC by the level of KIF14 expression. KIF14 is mainly involved in cell growth and/or maintenance biological processes. However, it is a pity that the main signaling pathway of KIF14’s major involvement in SCLC process has not been found, and further research needs to be carried out by future researchers.

MCM4 is a microchromosomal maintenance (MCM) protein complex, which is involved in cell cycle and cell replication. In cancer cells, abnormal expression of members of the MCM family has been reported in a wide range of cancers, and knockout of these genes can inhibit malignant phenotypes of cancer cells. Recent studies by Sanada et al. [44] have shown that siRNA-mediated MCM4 knockdown attenuates the invasiveness of lung adenocarcinoma cells. In addition, their team also confirmed that siRNA-mediated MCM4 knockdown enhanced the sensitivity of SCLC cells to cisplatin [45], suggesting that MCM4 could be used as a therapeutic target for SCLC.

The non-SMC lectin II complex subunit G2 (NCAPG2), a component of the lectin II complex, interacts with PLK1 to regulate correct chromosome separation. Similar to our results, NCAPG2 is highly expressed in a variety of cancers, such as liver cancer and NSCLC [46, 47]. According to our analysis, its overexpression leads to a short survival time, which is consistent with the results of previous studies in other cancer type [47]. At the same time, we also confirmed that NACPG2 affects tumor process by participating in cell growth, and Meng et al. also proposed the conclusion that NCAPG2 can affect cell proliferation in their study [46]. These findings confirm that NCAPG2 is both an oncogene of SCLC and a biomarker that predicts poor prognosis in patients.

Never in mitosis- (NIMA-) related kinase 2 (NEK2) is a member of the serine/threonine kinase family, mainly involved in regulating the cell cycle progression and microtubule organization and stabilization [47]. Recent reports have consistently identified high expression of NEK2 in various cancer types, including breast cancer, cervical cancer, liver cancer, and lung cancer [48–52]. In addition, NEK2 overexpression was enhanced in advanced lung adenocarcinoma, suggesting a role for NEK2 in tumor progression [53]. In terms of overall survival, patients with high NEK2 expression in NSCLC and its subtypes of lung adenocarcinoma have a poor prognosis [54, 55]. However, there are few related studies on NEK in SCLC, and further clinical studies and basic researches are needed. Therefore, the expression of NEK2 may be helpful in determining tumor progression and disease prognosis.

Forkhead Box M1 (FOXM1) is a member of the Forkhead family of proteins and is involved in cell cycle regulation. FOXM1 expression was low in quiescent cells but elevated in most tumors, including liver cancer and gastric cancer [56, 57]. Moreover, Hu et al. [58] found that FOXM1 and KIF4A proteins were upregulated in clinical liver cancer tissue samples, which was positively correlated with poor prognosis of patients with hepatocellular carcinoma. Currently, there are few studies on the biological function and clinical significance of FOXM1 in SCLC. Notably, Liang et al. [59] found that FOXM1 knockout inhibited SCLC formation in mouse models through increased levels of neuroendocrine markers Ascl1 and Cgrp and decreased levels of Yap1. In addition, this study also confirmed that SCLC with high FOXM1 expression was significantly associated with reduced clinical stage, extracthoracic metastasis, and OS with shorter progression-free survival. These evidences support the application of FOXM1 as a prognostic biomarker and potential molecular target for SCLC.

According to the Oncomine analysis results, the 10 hub genes identified in this work were able to distinguish SCLC samples from normal ones, thus suggesting their potential use as diagnostic biomarkers. We used Kaplan-Meier analysis to obtain the relations between these hub genes and lung cancer prognosis and found that the median survival time of the group with high expression of CDC20, CENPU, CHEK1, DTL, KIF14, MCM4, NCAPG2, NEK2, and FOXM1 is significantly shorter than that of the group with low expression of these genes. This finding suggests that the overexpression of these genes could predict the poor prognosis of patients with SCLC.

At present, several scholars have studied the GEO database of SCLC and obtained a certain number of DEGs. Liao et al. [60] screened out 5 highly expressed hub genes (NDC80, BUB1B, PLK1, CDC20, and MAD2L1) from 4 datasets (GSE60052, GSE43346, GSE15240, and GSE6044), and the cell cycle pathway was considered to be the main pathway for the diagnosis and treatment of SCLC of these five hub genes. Mao et al. [9] studied two databases, namely, GSE6044 and GSE19945, and finally identified 32 miRNAs and 32 regulated genes by using the bioinformatics platform “miRNAWalk.” It was suggested that bioinformatics analysis may contribute to a better understanding of the roles of DEGs, DEM, and miRNA genes in cell proliferation and signal transduction, and their related hub genes can be used as biomarkers for diagnosis and prognosis of SCLC, as well as potential drug targets. Wen et al. [61] studied two databases (GSE11969 and GSE6044) and finally confirmed 10 hub genes (TOP2A, PCNA, RFC4, CHEK1, TYMS, MCM2, CDC20, CDKN3, MCM3, and CDC6). At the same time, the signaling pathways of these 10 hub genes were also analyzed to provide molecular targets and diagnostic markers for the treatment and early diagnosis of SCLC. In the latest study, Chen et al. [11] studied 3 databases, namely, GSE40275, GSE99316, and GSE6052, and used GEO2R analysis tool to filter DEGs and Funrich for functional annotation; then, 8 hub genes (CDC20, BUB1, TOP2A, RRM2, CCNA2, UBE2C, MAD2L1, and BUB1B) were identified by the PPI network, module analysis, and mRNA expression level verification of hub genes in the ONCOMINE database. RT-qPCR was also used in clinical samples to verify that these hub genes may become prognostic markers or therapeutic targets for SCLC.

At present, the sample size of SCLC mRNA tissue in GEO database is relatively limited. GSE6044 and GSE60052 have been studied for several times. We need to analyze more datasets in the database (with a large enough sample size) to expand the number of DEGs. Therefore, new datasets GSE30219, GSE9316, and GSE149507 were selected for analysis in this study. Our innovation lies in the prognostic analysis, namely, the survival curve analysis, of the newly screened 10 hub genes, which improves the clinical application value of hub genes. But there are limitations to our study. Firstly, the data used in this study were all from public databases, but the quality of the data could not be evaluated. Secondly, the sample size of relevant data is relatively small. Third, it has not been validated in clinical samples. Therefore, a lot of valuable biological information may be ignored in our research. Finally, all 10 hub genes are overexpressed in SCLC, but the related mechanisms have not been fully clarified. Therefore, more molecular evidence is needed. In addition, current SCLC studies lack prognostic data of these hub genes, such as survival curves, which limits the clinical application value of hub genes. In this paper, the expression levels of 8 pivotal genes were analyzed. Whether these hub genes can be used as biomarkers or therapeutic targets for SCLC needs further study.

5. Conclusion

In summary, we aimed to find DEGs associated with the carcinogenesis and progression of SCLC. The DEGs we obtained revealed a significant function in the occurrence and metastasis of SCLC. This work provides new molecular targets at the genetic level, as well as new insights into precision SCLC treatment. Further experiments are necessary to verify the results.

Data Availability

The datasets analyzed during current study are available from the corresponding author upon reasonable request.

The patients involved in the database have obtained ethical approval.

Conflicts of Interest

The authors have no financial or commercial conflicts of interest about our study.

Authors’ Contributions

Huanqing Liu performed statistical analysis and performed the data analysis, and Tingting Li was responsible for the quality control of data and performed literature research. Xunda Ye performed data interpretation. Jun Lyu contributed to the study concept and study design. All authors contributed to writing of the manuscript and approved the final version.

Acknowledgments

We thank all colleagues involved in the study for their contributions. We acknowledge GEO database for providing their platforms and contributors for uploading their meaningful datasets.

Supplementary Materials

Table S1: gene expression microarray datasets for SCLC in GEO datasets, including GSE30219, GSE99316, and GSE149507. Table S2: information for the DEGs identified from the GEO datasets (, ). Table S3: enrichment analysis of the upregulated DEGs. Table S4: enrichment analysis of the downregulated DEGs. Table S5: enrichment analysis of the DEGs of significant module. (Supplementary Materials)

References

C. Zappa and S. A. Mousa, “Non-small cell lung cancer: current treatment and future advances,” Translational Lung Cancer Research, vol. 5, no. 3, pp. 288–300, 2016.
View at: Publisher Site | Google Scholar
D. Xie, R. Marks, M. Zhang et al., “Nomograms predict overall survival for patients with small-cell lung cancer incorporating pretreatment peripheral blood markers,” Journal of Thoracic Oncology, vol. 10, no. 8, pp. 1213–1220, 2015.
View at: Publisher Site | Google Scholar
S. N. Waqar and D. Morgensztern, “Treatment advances in small cell lung cancer (SCLC),” Pharmacology & Therapeutics, vol. 180, pp. 16–23, 2017.
View at: Publisher Site | Google Scholar
N. B. Prasad, H. Somervell, R. P. Tufano et al., “Identification of genes differentially expressed in benign versus malignant thyroid tumors,” Clinical Cancer Research, vol. 14, no. 11, pp. 3327–3337, 2008.
View at: Publisher Site | Google Scholar
S. Udhaya Kumar, D. Thirumal Kumar, R. Bithia et al., “Analysis of differentially expressed genes and molecular pathways in familial hypercholesterolemia involved in atherosclerosis: a systematic and bioinformatics approach,” Frontiers in Genetics, vol. 11, p. 734, 2020.
View at: Publisher Site | Google Scholar
J. Wan, S. Jiang, Y. Jiang et al., “Data mining and expression analysis of differential lncRNA ADAMTS9-AS1 in prostate cancer,” Frontiers in Genetics, vol. 10, p. 1377, 2020.
View at: Publisher Site | Google Scholar
D. Fu, B. Zhang, L. Yang, S. Huang, and W. Xin, “Development of an immune-related risk signature for predicting prognosis in lung squamous cell carcinoma,” Frontiers in Genetics, vol. 11, p. 978, 2020.
View at: Publisher Site | Google Scholar
S. Udhaya Kumar, D. Thirumal Kumar, R. Siva et al., “Dysregulation of signaling pathways due to differentially expressed genes from the B-cell transcriptomes of systemic lupus erythematosus patients - a bioinformatics approach,” Frontiers in Bioengineering and Biotechnology, vol. 30, no. 8, p. 276, 2020.
View at: Publisher Site | Google Scholar
K. A. Mitchell, A. Zingone, L. Toulabi, J. Boeckelman, and B. M. Ryan, “Comparative transcriptome profiling reveals coding and noncoding RNA differences in NSCLC from African Americans and European Americans,” Clinical Cancer Research, vol. 23, no. 23, pp. 7412–7425, 2017.
View at: Publisher Site | Google Scholar
Y. Mao, P. Xue, L. Li et al., “Bioinformatics analysis of mRNA and miRNA microarray to identify the key miRNA-gene pairs in small-cell lung cancer,” Molecular Medicine Reports, vol. 20, no. 3, pp. 2199–2208, 2019.
View at: Publisher Site | Google Scholar
X. Chen, L. Wang, X. Su, S. Y. Luo, X. Tang, and Y. Huang, “Identification of potential target genes and crucial pathways in small cell lung cancer based on bioinformatic strategy and human samples,” PLoS One, vol. 15, no. 11, p. e0242194, 2020.
View at: Publisher Site | Google Scholar
W. D. Travis, “Update on small cell carcinoma and its differentiation from squamous cell carcinoma and other non-small cell carcinomas,” Modern Pathology, vol. 25, Suppl 1, pp. S18–S30, 2012.
View at: Publisher Site | Google Scholar
V. Sundaresan, V. T. Lin, F. Liang et al., “Significantly mutated genes and regulatory pathways in SCLC--a meta-analysis,” Cancer Genetics, vol. 216-217, pp. 20–28, 2017.
View at: Publisher Site | Google Scholar
L. S. Penna, J. A. P. Henriques, and D. Bonatto, “Anti-mitotic agents: are they emerging molecules for cancer treatment?” Pharmacology & Therapeutics, vol. 173, pp. 67–82, 2017.
View at: Publisher Site | Google Scholar
H. Zhang, C. L. Christensen, R. Dries et al., “CDK7 inhibition potentiates genome instability triggering anti-tumor immunity in small cell lung cancer,” Cancer Cell, vol. 37, no. 1, pp. 37–54.e9, 2020.
View at: Publisher Site | Google Scholar
L. Wang, J. Zhang, L. Wan, X. Zhou, Z. Wang, and W. Wei, “Targeting Cdc20 as a novel cancer therapeutic strategy,” Pharmacology & Therapeutics, vol. 151, pp. 141–151, 2015.
View at: Publisher Site | Google Scholar
T. Kato, Y. Daigo, M. Aragaki, K. Ishikawa, M. Sato, and M. Kaji, “Overexpression of CDC20 predicts poor prognosis in primary non-small cell lung cancer patients,” Journal of Surgical Oncology, vol. 106, no. 4, pp. 423–430, 2012.
View at: Publisher Site | Google Scholar
T. Kidokoro, C. Tanikawa, Y. Furukawa, T. Katagiri, Y. Nakamura, and K. Matsuda, “CDC20, a potential cancer therapeutic target, is negatively regulated by p53,” Oncogene, vol. 27, no. 11, pp. 1562–1571, 2008.
View at: Publisher Site | Google Scholar
S. H. Hanissian, U. Akbar, B. Teng et al., “cDNA cloning and characterization of a novel gene encoding the MLF1-interacting protein MLF1IP,” Oncogene, vol. 23, no. 20, pp. 3700–3707, 2004.
View at: Publisher Site | Google Scholar
L. Zhang, G. Ji, Y. Shao et al., “MLF1 interacting protein: a potential gene therapy target for human prostate cancer?” Medical Oncology, vol. 32, no. 2, p. 454, 2015.
View at: Publisher Site | Google Scholar
D. P. Huang and R. C. Luo, “MLF1IP is correlated with progression and prognosis in luminal breast cancer,” Biochemical and Biophysical Research Communications, vol. 477, no. 4, pp. 923–926, 2016.
View at: Publisher Site | Google Scholar
S. Wang, B. Liu, J. Zhang et al., “Centromere protein U is a potential target for gene therapy of human bladder cancer,” Oncology Reports, vol. 38, no. 2, pp. 735–744, 2017.
View at: Publisher Site | Google Scholar
J. Li, Z. G. Wang, L. B. Pang, R. H. Zhang, and Y. Y. Wang, “Reduced CENPU expression inhibits lung adenocarcinoma cell proliferation and migration through PI3K/AKT signaling,” Bioscience, Biotechnology, and Biochemistry, vol. 83, no. 6, pp. 1077–1084, 2019.
View at: Publisher Site | Google Scholar
K. A. Cole, J. Huggins, M. Laquaglia et al., “RNAi screen of the protein kinome identifies checkpoint kinase 1 (CHK1) as a therapeutic target in neuroblastoma,” PNAS, vol. 108, no. 8, pp. 3336–3341, 2011.
View at: Publisher Site | Google Scholar
H. O. Ebili, V. O. Iyawe, K. R. Adeleke et al., “Checkpoint kinase 1 expression predicts poor prognosis in Nigerian breast cancer patients,” Molecular Diagnosis & Therapy, vol. 22, no. 1, pp. 79–90, 2018.
View at: Publisher Site | Google Scholar
A. Alcaraz-Sanabria, C. Nieto-Jiménez, V. Corrales-Sánchez et al., “Synthetic lethality interaction between Aurora kinases and CHEK1 inhibitors in ovarian cancer,” Molecular Cancer Therapeutics, vol. 16, no. 11, pp. 2552–2562, 2017.
View at: Publisher Site | Google Scholar
H. Gali-Muhtasib, D. Kuester, C. Mawrin et al., “Thymoquinone triggers inactivation of the stress response pathway sensor CHEK1 and contributes to apoptosis in colorectal cancer cells,” Cancer Research, vol. 68, no. 14, pp. 5609–5618, 2008.
View at: Publisher Site | Google Scholar
F. Doerr, J. George, A. Schmitt et al., “Targeting a non-oncogene addiction to the ATR/CHK1 axis for the treatment of small cell lung cancer,” Scientific Reports, vol. 7, no. 1, p. 15511, 2017.
View at: Publisher Site | Google Scholar
H. Lee, Y. Kim, J. H. Jeong, J. H. Ryu, and W. Y. Kim, “ATM/CHK/p53 pathway dependent chemopreventive and therapeutic activity on lung cancer by pterostilbene,” PLoS One, vol. 11, no. 9, article e0162335, 2016.
View at: Publisher Site | Google Scholar
T. J. Slenn, B. Morris, C. G. Havens, R. M. Freeman Jr., T. S. Takahashi, and J. C. Walter, “Thymine DNA glycosylase is a CRL4Cdt2 substrate,” The Journal of Biological Chemistry, vol. 289, no. 33, pp. 23043–23055, 2014.
View at: Publisher Site | Google Scholar
H. Kobayashi, S. Komatsu, D. Ichikawa et al., “Overexpression of denticleless E3 ubiquitin protein ligase homolog (DTL) is related to poor outcome in gastric carcinoma,” Oncotarget, vol. 6, no. 34, pp. 36615–36624, 2015.
View at: Publisher Site | Google Scholar
Y. C. Chen, I. S. Chen, G. J. Huang et al., “Targeting DTL induces cell cycle arrest and senescence and suppresses cell growth and colony formation through TPX2 inhibition in human hepatocellular carcinoma cells,” Oncotargets and Therapy, vol. 11, pp. 1601–1616, 2018.
View at: Publisher Site | Google Scholar
J. Perez-Peña, V. Corrales-Sánchez, E. Amir, A. Pandiella, and A. Ocana, “Ubiquitin-conjugating enzyme E2T (UBE2T) and denticleless protein homolog (DTL) are linked to poor outcome in breast and lung cancers,” Scientific Reports, vol. 7, no. 1, p. 17530, 2017.
View at: Publisher Site | Google Scholar
G. Wu, L. Zhou, L. Khidr et al., “A novel role of the chromokinesin Kif4A in DNA damage response,” Cell Cycle, vol. 7, no. 13, pp. 2013–2020, 2008.
View at: Publisher Site | Google Scholar
Y. M. Lee and W. Kim, “Association of human kinesin superfamily protein member 4 with BRCA2-associated factor 35,” The Biochemical Journal, vol. 374, Part 2, pp. 497–503, 2003.
View at: Publisher Site | Google Scholar
P. F. Hou, T. Jiang, F. Chen et al., “KIF4A facilitates cell proliferation via induction of p21-mediated cell cycle progression and promotes metastasis in colorectal cancer,” Cell Death & Disease, vol. 9, no. 5, p. 477, 2018.
View at: Publisher Site | Google Scholar
C. L. Zhu, D. Z. Cheng, F. Liu et al., “Hepatitis B virus upregulates the expression of kinesin family member 4A,” Molecular Medicine Reports, vol. 12, no. 3, pp. 3503–3507, 2015.
View at: Publisher Site | Google Scholar
M. Taniwaki, A. Takano, N. Ishikawa et al., “Activation of KIF4A as a prognostic biomarker and therapeutic target for lung cancer,” Clinical Cancer Research, vol. 13, 22 Part 1, pp. 6624–6631, 2007.
View at: Publisher Site | Google Scholar
T. W. Corson, A. Huang, M. S. Tsao, and B. L. Gallie, “KIF14 is a candidate oncogene in the 1q minimal region of genomic gain in multiple cancers,” Oncogene, vol. 24, no. 30, pp. 4741–4753, 2005.
View at: Publisher Site | Google Scholar
J. Madhavan, K. Coral, K. Mallikarjuna et al., “High expression of KIF14 in retinoblastoma: association with older age at diagnosis,” Investigative Ophthalmology & Visual Science, vol. 48, no. 11, pp. 4901–4906, 2007.
View at: Publisher Site | Google Scholar
X. Chen, D. Zhang, F. Jiang et al., “Prognostic prediction using a stemness index-related signature in a cohort of gastric cancer,” Frontiers in Molecular Biosciences, vol. 7, p. 570702, 2020.
View at: Publisher Site | Google Scholar
P. F. Hung, T. M. Hong, Y. C. Hsu et al., “The motor protein KIF14 inhibits tumor growth and cancer metastasis in lung adenocarcinoma,” PLoS One, vol. 8, no. 4, article e61664, 2013.
View at: Publisher Site | Google Scholar
T. W. Corson, C. Q. Zhu, S. K. Lau, F. A. Shepherd, M. S. Tsao, and B. L. Gallie, “KIF14 messenger RNA expression is independently prognostic for outcome in lung cancer,” Clinical Cancer Research, vol. 13, no. 11, pp. 3229–3234, 2007.
View at: Publisher Site | Google Scholar
H. Sanada, N. Seki, K. Mizuno et al., “Involvement of dual strands of miR-143 (miR-143-5p and miR-143-3p) and their target oncogenes in the molecular pathogenesis of lung adenocarcinoma,” International Journal of Molecular Sciences, vol. 20, no. 18, p. 4482, 2019.
View at: Publisher Site | Google Scholar
S. Misono, K. Mizuno, T. Suetsugu et al., “Molecular signature of small cell lung cancer after treatment failure: the MCM complex as therapeutic target,” Cancers, vol. 13, no. 6, p. 1187, 2021.
View at: Publisher Site | Google Scholar
F. Meng, S. Zhang, R. Song et al., “NCAPG2 overexpression promotes hepatocellular carcinoma proliferation and metastasis through activating the STAT3 and NF-κB/miR-188-3p pathways,” eBioMedicine, vol. 44, pp. 237–249, 2019.
View at: Publisher Site | Google Scholar
P. Zhan, G. M. Xi, B. Zhang et al., “NCAPG2 promotes tumour proliferation by regulating G2/M phase and associates with poor prognosis in lung adenocarcinoma,” Journal of Cellular and Molecular Medicine, vol. 21, no. 4, pp. 665–676, 2017.
View at: Publisher Site | Google Scholar
Y. Jeong, J. Lee, K. Kim, J. C. Yoo, and K. Rhee, “Characterization of NIP2/centrobin, a novel substrate of Nek 2, and its potential role in microtubule stabilization,” Journal of Cell Science, vol. 120, pp. 2106–2116, 2007.
View at: Publisher Site | Google Scholar
M. Marina and H. I. Saavedra, “Nek2 and Plk4: prognostic markers, drivers of breast tumorigenesis and drug resistance,” Frontiers in Bioscience (Landmark Ed)., vol. 19, no. 2, pp. 352–365, 2014.
View at: Publisher Site | Google Scholar
C. Naro, F. Barbagallo, P. Chieffi, C. F. Bourgeois, M. P. Paronetto, and C. Sette, “The centrosomal kinase NEK2 is a novel splicing factor kinase involved in cell survival,” Nucleic Acids Research, vol. 42, no. 5, pp. 3218–3227, 2014.
View at: Publisher Site | Google Scholar
G. Y. Wubetu, Y. Morine, H. Teraoku et al., “High NEK2 expression is a predictor of tumor recurrence in hepatocellular carcinoma patients after hepatectomy,” Anticancer Research, vol. 36, no. 2, pp. 757–762, 2016.
View at: Google Scholar
X. Zhong, X. Guan, W. Liu, and L. Zhang, “Aberrant expression of NEK2 and its clinical significance in non-small cell lung cancer,” Oncology Letters, vol. 8, no. 4, pp. 1470–1476, 2014.
View at: Publisher Site | Google Scholar
G. Bidkhori, Z. Narimani, S. Hosseini Ashtiani, A. Moeini, A. Nowzari-Dalini, and A. Masoudi-Nejad, “Reconstruction of an integrated genome-scale co-expression network reveals key modules involved in lung adenocarcinoma,” PLoS One, vol. 8, no. 7, article e67552, 2013.
View at: Publisher Site | Google Scholar
X. Zhong, X. Guan, Q. Dong, S. Yang, W. Liu, and L. Zhang, “Examining Nek2 as a better proliferation marker in non-small cell lung cancer prognosis,” Tumour Biology, vol. 35, no. 7, pp. 7155–7162, 2014.
View at: Publisher Site | Google Scholar
W. Zhou, Y. Yang, J. Xia et al., “NEK2 induces drug resistance mainly through activation of efflux drug pumps and is associated with poor prognosis in myeloma and other cancers,” Cancer Cell, vol. 23, no. 1, pp. 48–62, 2013.
View at: Publisher Site | Google Scholar
S. M. E. Weiler, F. Pinna, T. Wolf et al., “Induction of chromosome instability by activation of yes-associated protein and forkhead box M1 in liver cancer,” Gastroenterology, vol. 152, no. 8, pp. 2037–2051.e22, 2017.
View at: Publisher Site | Google Scholar
P. Francica, L. Nisa, D. M. Aebersold et al., “Depletion of FOXM1 via MET targeting underlies establishment of a DNA damage-induced senescence program in gastric cancer,” Clinical Cancer Research, vol. 22, no. 21, pp. 5322–5336, 2016.
View at: Publisher Site | Google Scholar
G. Hu, Z. Yan, C. Zhang et al., “FOXM1 promotes hepatocellular carcinoma progression by regulating KIF4A expression,” Journal of Experimental & Clinical Cancer Research, vol. 38, no. 1, p. 188, 2019.
View at: Publisher Site | Google Scholar
S. K. Liang, C. C. Hsu, H. L. Song et al., “FOXM1 is required for small cell lung cancer tumorigenesis and associated with poor clinical prognosis,” Oncogene, vol. 40, no. 30, pp. 4847–4858, 2021.
View at: Publisher Site | Google Scholar
Y. Liao, G. Yin, X. Wang, P. Zhong, X. Fan, and C. Huang, “Identification of candidate genes associated with the pathogenesis of small cell lung cancer via integrated bioinformatics analysis,” Oncology Letters, vol. 18, no. 4, pp. 3723–3733, 2019.
View at: Publisher Site | Google Scholar
P. Wen, T. Chidanguro, Z. Shi et al., “Identification of candidate biomarkers and pathways associated with SCLC by bioinformatics analysis,” Molecular Medicine Reports, vol. 18, no. 2, pp. 1538–1550, 2018.
View at: Publisher Site | Google Scholar

Copyright

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

PDF Download Citation

Download other formats

Order printed copies

Views

832

Downloads

1093

Citations