Genetic variant of IRAK2 in the toll-like receptor signaling pathway and survival of non-small cell lung cancer

Both IRAK3 and IRAK1 Activate the MyD88-TRAF6 Pathway in Zebrafish

IL-1R-associated kinases (IRAKs) are signal transducers of the TLR/IL-1R-MyD88-TRAF6 pathways. Vertebrates possess two IRAK lineages, IRAK1/2/3 and IRAK4. In mammals, IRAK4/IRAK1 and IRAK4/IRAK2 are pathway enhancers, whereas IRAK3 is a repressor. However, in bony fish, IRAK2 is absent, and it remains elusive how fish IRAK1/3/4 functionally differ from their mammalian counterparts. In this study, we explored this using the zebrafish model. First, we showed that in human 293T cells, zebrafish IRAK1 and IRAK4 were components of the Myddosome (MyD88-IRAK4-IRAK1) complex, with IRAK1 serving as a potent pathway enhancer. Then, we discovered two zebrafish IRAK3 variants: one (IRAK3a) contains an N-terminal Death domain, a middle pseudokinase domain, and a C-terminal TRAF6-binding domain, whereas the other (IRAK3b) lost both the kinase and TRAF6-binding domains. This truncation of IRAK3 variants could be a conserved phenomenon in fish, because it is also observed in trout and grass carp. We proceeded to show that zebrafish IRAK3a acts as a pathway enhancer by binding with MyD88 and TRAF6, but its activity is milder than IRAK1, possibly because it has no kinase activity. Zebrafish IRAK3b, however, plays a sheer negative role, apparently because of its lack of kinase and TRAF6-binding domains. Moreover, zebrafish IRAK3a/3b inhibit the activity of IRAK1/4, not by interacting with IRAK1/4 but possibly by competing for MyD88 and TRAF6. Finally, we have verified the essential activities of zebrafish IRAK1/3a/3b/4 in zebrafish cells and embryos. In summary, to our knowledge, our findings provide new insights into the molecular functions of fish IRAKs and the evolution of the IRAK functional modes in vertebrates.

IRAK2, an Immune and Radiation-Response Gene, Correlates with Advanced Disease Features but Predicts Higher Post-Irradiation Local Control in Non-Metastatic and Resected Oral Cancer Patients

Gene Ontology (GO) analysis can provide a comprehensive function analysis for investigating genes, allowing us to identify the potential biological roles of genes. The present study conducted GO analysis to explore the biological function of IRAK2 and performed a case analysis to define its clinical role in disease progression and mediating tumor response to RT. Methods: We performed a GO enrichment analysis on the RNA-seq data to validate radiation-induced gene expression. A total of 172 I-IVB specimens from oral squamous cell carcinoma patients were collected for clinical analysis, from which IRAK2 expression was analyzed by immunohistochemistry. This was a retrospective study conducted between IRAK2 expression and the outcomes of oral squamous cell carcinoma patients after radiotherapy treatment. We conducted Gene Ontology (GO) analysis to explore the biological function of IRAK2 and performed a case analysis to define its clinical role in mediating tumor response to radiotherapy. GO enrichment analysis to validate radiation-induced gene expression was performed. Clinically, 172 stage I-IVB resected oral cancer patients were used to validate IRAK2 expression in predicting clinical outcomes. GO enrichment analysis showed that IRAK2 is involved in 10 of the 14 most enriched GO categories for post-irradiation biological processes, focusing on stress response and immune modulation. Clinically, high IRAK2 expression was correlated with adverse disease features, including pT3-4 status (p = 0.01), advanced overall stage (p = 0.02), and positive bone invasion (p = 0.01). In patients who underwent radiotherapy, the IRAK2-high group was associated with reduced post-irradiation local recurrence (p = 0.025) compared to the IRAK2-low group. IRAK2 plays a crucial role in the radiation-induced response. Patients with high IRAK2 expression demonstrated more advanced disease features but predicted higher post-irradiation local control in a clinical setting. These findings support IRAK2 as a potential predictive biomarker for radiotherapy response in non-metastatic and resected oral cancer patients.

Integrative weighted molecular network construction from transcriptomics and genome wide association data to identify shared genetic biomarkers for COPD and lung cancer

Chronic obstructive pulmonary disease (COPD) is a multifactorial progressive airflow obstruction in the lungs, accounting for high morbidity and mortality across the world. This study aims to identify potential COPD blood-based biomarkers by analyzing the dysregulated gene expression patterns in blood and lung tissues with the help of robust computational approaches. The microarray gene expression datasets from blood (136 COPD and 6 controls) and lung tissues (16 COPD and 19 controls) were analyzed to detect shared differentially expressed genes (DEGs). Then these DEGs were used to construct COPD protein network-clusters and functionally enrich them against gene ontology annotation terms. The hub genes in the COPD network clusters were then queried in GWAS catalog and in several cancer expression databases to explore their pathogenic roles in lung cancers. The comparison of blood and lung tissue datasets revealed 63 shared DEGs. Of these DEGs, 12 COPD hub gene-network clusters (SREK1, TMEM67, IRAK2, MECOM, ASB4, C1QTNF2, CDC42BPA, DPF3, DET1, CCDC74B, KHK, and DDX3Y) connected to dysregulations of protein degradation, inflammatory cytokine production, airway remodeling, and immune cell activity were prioritized with the help of protein interactome and functional enrichment analysis. Interestingly, IRAK2 and MECOM hub genes from these COPD network clusters are known for their involvement in different pulmonary diseases. Additional COPD hub genes like SREK1, TMEM67, CDC42BPA, DPF3, and ASB4 were identified as prognostic markers in lung cancer, which is reported in 1% of COPD patients. This study identified 12 gene network- clusters as potential blood based genetic biomarkers for COPD diagnosis and prognosis.

Identification and Validation of lncRNA-SNHG17 in Lung Adenocarcinoma: A Novel Prognostic and Diagnostic Indicator

Background

Lung cancer has the highest death rate among cancers globally. Accumulating evidence has indicated that cancer-related inflammation plays an important role in the initiation and progression of lung cancer. However, the prognosis, immunological role, and associated regulation axis of inflammatory response-related gene (IRRGs) in non-small-cell lung cancer (NSCLC) remains unclear.

Methods

In this study, we perform comprehensive bioinformatics analysis and constructed a prognostic inflammatory response-related gene (IRRGs) and related competing endogenous RNA (ceRNA) network. We also utilized the Pearson’s correlation analysis to determine the correlation between IRRGs expression and tumor mutational burden (TMB), microsatellite instability (MSI), tumor-immune infiltration, and the drug sensitivity in NSCLC. Growth curve and Transwell assay used to verify the function of SNHG17 on NSCLC progression.

Results

First, we found that IRRGs were significantly upregulated in lung cancer, and its high expression was correlated with poor prognosis; high expression of IRRGs was significantly correlated with the tumor stage and poor prognosis in lung cancer patients. Moreover, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment indicated that these IRRGs are mainly involved in the inflammatory and immune response-related signaling pathway in the progression of NSCLC. We utilized 10 prognostic-related genes to construct a prognostic IRRGs model that could predict the overall survival of lung adenocarcinoma (LUAD) patients possessing high specificity and accuracy. Our evidence demonstrated that IRRGs expression was significantly correlated with the TMB, MSI, immune-cell infiltration, and diverse cancer-related drug sensitivity. Finally, we identified the upstream regulatory axis of IRRGs in NSCLC, namely, lncRNA MIR503HG/SNHG17/miR-330-3p/regulatory axis. Finally, knockdown of SNHG17 expression inhibited lung adenocarcinoma (LUAD) cell proliferation and migration. Our findings confirmed that SNHG17 is a novel oncogenic lncRNA and may be a biomarker for the prognosis and diagnosis of LUAD.

Conclusion

DNA hypomethylation/lncRNA MIR503HG/SNHG17/microRNA-330-3p/regulatory axis may be a valuable biomarker for prognosis and is significantly correlated with immune cell infiltration in lung cancer.

IRAK2-NF-κB signaling promotes glycolysis-dependent tumor growth in pancreatic cancer

BACKGROUND

Metabolic reprogramming has emerged as a core hallmark of cancer, and cancer metabolism has long been equated with aerobic glycolysis. Moreover, hypoxia and the hypovascular tumor microenvironment (TME) are major hallmarks of pancreatic ductal adenocarcinoma (PDAC), in which glycolysis is imperative for tumor cell survival and proliferation. Here, we explored the impact of interleukin 1 receptor-associated kinase 2 (IRAK2) on the biological behavior of PDAC and investigated the underlying mechanism.

METHODS

The expression pattern and clinical relevance of IRAK2 was determined in GEO, TCGA and Ren Ji datasets. Loss-of-function and gain-of-function studies were employed to investigate the cellular functions of IRAK2 in vitro and in vivo. Gene set enrichment analysis, Seahorse metabolic analysis, immunohistochemistry and Western blot were applied to reveal the underlying molecular mechanisms.

RESULTS

We found that IRAK2 is highly expressed in PDAC patient samples and is related to a poor prognosis. IRAK2 knockdown led to a significant impairment of PDAC cell proliferation via an aberrant Warburg effect. Opposite results were obtained after exogenous IRAK2 overexpression. Mechanistically, we found that IRAK2 is critical for sustaining the activation of transcription factors such as those of the nuclear factor-κB (NF-κB) family, which have increasingly been recognized as crucial players in many steps of cancer initiation and progression. Treatment with maslinic acid (MA), a NF-κB inhibitor, markedly attenuated the aberrant oncological behavior of PDAC cells caused by IRAK2 overexpression.

CONCLUSIONS

Our data reveal a role of IRAK2 in PDAC metabolic reprogramming. In addition, we obtained novel insights into how immune-related pathways affect PDAC progression and suggest that targeting IRAK2 may serve as a novel therapeutic approach for PDAC.

Cysteine cathepsin C: a novel potential biomarker for the diagnosis and prognosis of glioma

Background

Cysteine cathepsin C encoded by the CTSC gene is an important member of the cysteine cathepsin family that plays a key role regulation of many types of tumors. However, whether CTSC is involved in the pathological process of glioma has not yet been reported. We comprehensively analyzed data from multiple databases and for the first time revealed a role and specific mechanism of action of CTSC in glioma, identifying it as a novel and efficient biomarker for the diagnosis and treatment of this brain tumor.

Methods

The expression of CTSC in glioma and its relationship with clinical characteristics and prognosis of patients with glioma were analyzed at different levels by using clinical sample information from several databases. CTSC expression levels in glioma and normal brain tissues, as well as in glioma cells and normal brain cells, was validated by real-time quantitative polymerase chain reaction (RT-qPCR). Gene set enrichment analysis (GSEA) was used to reveal the signaling pathways that CTSC may participate in. The connectivity map was used to reveal small molecules that may inhibit CTSC expression in glioma, and the putative effect of these compounds was verified by RT-qPCR.

Results

Our analyses showed that the expression of CTSC in glioma was higher than that in non-cancerous cells. GSEA showed that CTSC expression may regulate the malignant development of glioma through Toll-like receptor signaling pathways, pathways in cancer, and extracellular matrix receptor interaction signaling pathways. And we proved piperlongumine and scopoletin could inhibit CTSC expression in glioma cells.

Conclusions

CTSC may serve as an efficient molecular target for the diagnosis and therapy of glioma, thereby improving the poor prognosis of patients with glioma.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02417-6.

IRAK1 and IRAK4 as emerging therapeutic targets in hematologic malignancies

PURPOSE OF REVIEW

Cell intrinsic and extrinsic perturbations to inflammatory signaling pathways are a hallmark of development and progression of hematologic malignancies. The interleukin 1 receptor-associated kinases (IRAKs) are a family of related signaling intermediates (IRAK1, IRAK2, IRAK3, IRAK4) that operate at the nexus of multiple inflammatory pathways implicated in the hematologic malignancies. In this review, we explicate the oncogenic role of these kinases and review recent therapeutic advances in the dawning era of IRAK-targeted therapy.

RECENT FINDINGS

Emerging evidence places IRAK signaling at the confluence of adaptive resistance and oncogenesis in the hematologic malignancies and solid tissue tumors. Preclinical investigations nominate the IRAK kinases as targetable molecular dependencies in diverse cancers.

SUMMARY

IRAK-targeted therapies that have matriculated to early phase trials are yielding promising preliminary results. However, studies of IRAK kinase signaling continue to defy conventional signaling models and raise questions as to the design of optimal treatment strategies. Efforts to refine IRAK signaling mechanisms in the malignant context will inspire deliberate IRAK-targeted drug development and informed combination therapy.

Potentially functional variants of HBEGF and ITPR3 in GnRH signaling pathway genes predict survival of non-small cell lung cancer patients

The gonadotropin-releasing hormone (GnRH) signaling pathway controls reproductive functions and cancer growth and progression. However, few studies investigated roles of genetic variants of GnRH pathway genes in survival of patients with non-small cell lung cancer (NSCLC). Therefore, we first evaluated associations between 22,528 single-nucleotide polymorphisms (SNPs) in 101 GnRH pathway genes and survival of 1185 NSCLC patients using a dataset from Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial. We found 572 SNPs to be significantly associated with overall survival (OS) of NSCLC (P ≤ 0.05, Bayesian false discovery probability ≤0.80). We then validated these SNPs in another dataset with 984 NSCLC patients from the Harvard Lung Cancer Susceptibility Study. Finally, two independent SNPs (HBEGF rs4150236G>A and ITPR3 rs116454384C>T) remained significantly associated with NSCLC OS in the combined analysis with hazards ratios of 0.84 (95% confidence interval = 0.76-0.92, P = 0.0003) and 0.85 (0.78-0.94, 0.0012), respectively; their genetic score (the number of protective genotypes) was associated with a better OS and disease-specific survival (P_trend = 0.0002 and 0.0001, respectively). Further expression quantitative trail loci analysis showed a significant correlation between ITPR3 rs116454384 T allele and an increased mRNA expression level in both whole blood and normal lung tissue, and high ITPR3 mRNA expression levels in tumors were associated with a better survival of NSCLC patients. Because ITPR3 mutations were rare in tumors, ITPR3 rs116454384C>T likely had an effect on cancer progression by regulating the gene expression. Therefore, genetic variants of HBEGF rs4150236G>A and ITPR3 rs116454384C>T may be predictors for NSCLC survival, but HBEGF rs4150236G>A functional relevance remains to be determined.

IRAK2, an IL1R/TLR Immune Mediator, Enhances Radiosensitivity via Modulating Caspase 8/3-Mediated Apoptosis in Oral Squamous Cell Carcinoma

Predicting and overcoming radioresistance are crucial in radiation oncology, including in managing oral squamous cell carcinoma (OSCC). First, we used RNA-sequence to compare expression profiles of parent OML1 and radioresistant OML1-R OSCC cells in order to select candidate genes responsible for radiation sensitivity. We identified IRAK2, a key immune mediator of the IL-1R/TLR signaling, as a potential target in investigating radiosensitivity. In four OSCC cell lines, we observed that intrinsically low IRAK2 expression demonstrated a radioresistant phenotype (i.e., OML1-R and SCC4), and vice versa (i.e., OML1 and SCC25). Next, we overexpressed IRAK2 in low IRAK2-expression OSCC cells and knocked it down in high IRAK2-expression cells to examine changes of irradiation response. After ionizing radiation (IR) exposure, IRAK2 overexpression enhanced the radiosensitivity of radioresistant cells and synergistically suppressed OSCC cell growth both in vitro and in vivo, and vice versa. We found that IRAK2 overexpression restored and enhanced radiosensitivity by enhancing IR-induced cell killing via caspase-8/3-dependent apoptosis. OSCC patients with high IRAK2 expression had better post-irradiation local control than those with low expression (i.e., 87.4% vs. 60.0% at five years, P = 0.055), showing that IRAK2 expression was associated with post-radiation recurrence. Multivariate analysis confirmed high IRAK2 expression as an independent predictor for local control (HR, 0.11; 95% CI, 0.016 – 0.760; P = 0.025). In conclusion, IRAK2 enhances radiosensitivity, via modulating caspase 8/3-medicated apoptosis, potentially playing double roles as a predictive biomarker and a novel therapeutic target in OSCC.

Genetic variants in the human leukocyte antigen region and survival of Chinese patients with non-small cell lung carcinoma

Human leukocyte antigen (HLA) is highly polymorphic, driving antigen presentation, complement cascade and leukocyte maturation against cancer cells. Therefore, we extracted genotyping data in the HLA region from an ongoing Chinese genome-wide association study of non-small cell lung cancer (NSCLC). Using deep sequencing data of 10 689 healthy Han Chinese, we imputed for untyped genetic variants in the HLA region, followed by a two-stage survival analysis of 1531 NSCLC patients. In the discovery stage of 758 patients, we identified 301 out of 15 138 single-nucleotide polymorphisms to be independently associated with overall survival [P < 0.05 and Bayesian false-discovery probability < 0.8]. In further validation of another 773 patients, we confirmed chromosome 6p21, rs241424 (located at intron 3 of TAP2) and rs6457642 as two independent survival predictors. In the combined analysis of 1531 NSCLC patients, rs241424 G>A and rs6457642 C>T were associated with a hazards ratio of 1.26 [95% confidence interval (CI) = 1.14-1.40 and P = 4.04 × 10-6] and 0.76 (95% CI = 0.66-0.87 and P = 1.16 × 10-4), respectively. The analysis of publically available ChIP-sequencing and Hi-C data found that the rs241424 locus was involved in potential cis-regulatory element by a long-range interaction with the HLA-DQA1 promoter. Additional expression quantitative trait loci analysis showed that the rs241424 G>A change decreased HLA-DQA1 mRNA expression. Furthermore, expression levels of HLA-DQA1 were lower in lung cancer tissues than in adjacent normal tissues, and the lower expression was associated with a worse prognosis for patients with lung adenocarcinoma. Collectively, HLA genetic variants may modulate OS of NSCLC patients, possibly via a mechanism of long-range promoter interaction regulating HLA-DQA1 expression.

Targeting Toll like Receptors in Cancer: Role of TLR Natural and Synthetic Modulators

Background:

Toll like receptors (TLRs) are a group of transmembrane receptors belonging to the broad class pattern recognition receptors (PRR), involved in recognition of Pathogen Associated Molecular Patterns (PAMPs) thereby inducing an immune response. Apart from these exogenous PAMPs, numerous endogenous PAMPs are also ligands for various TLRs thereby activating the TLR dependent immune response, subsequently leading to the onset of an inflammatory response. Prolonged activation of TLR by these endogenous PAMPs leads to chronic inflammatory insults to the body and which in turn alters the proliferative patterns of the cells, which ultimately leads to the development of cancer.

Objectives:

The present review aims to provide a detailed outline of the differential roles of various TLRs in cancer and the possible use of them as a therapeutic target.

Methods:

Data were collected from PubMed/Sciencedirect/Web of Science database and sorted; the latest literature on TLRs was incorporated in the review.

Results:

Among the different TLRs, few are reported to be anti-neoplastic, which controls the cell growth and multiplication in response to the endogenous signals. On the contrary, numerous studies have reported the procarcinogenic potentials of TLRs. Hence, TLRs have emerged as a potential target for the prevention and treatment of various types of cancers. Several molecules, such as monoclonal antibodies, small molecule inhibitors and natural products have shown promising anticancer potential by effectively modulating the TLR signalling.

Conclusion:

Toll-like receptors play vital roles in the process of carcinogenesis, hence TLR targeting is a promising approach for cancer prevention.

Novel genetic variants of SYK and ITGA1 related lymphangiogenesis signaling pathway predict non-small cell lung cancer survival

Although lymphangiogenesis is a vital step in lung cancer metastasis, the association between lymphangiogenesis and non-small cell lung cancer (NSCLC) survival remains unclear. Since single-nucleotide polymorphisms (SNPs) have been reported to predict NSCLC survival, we investigated associations between SNPs in lymphangiogenesis-related pathway genes and NSCLC survival in a discovery genotyping dataset of 1,185 patients from the Prostate, Lung, Colorectal and Ovarian (PLCO) Cancer Screening Trial and validated the findings in another genotyping dataset of 984 patients from the Harvard Lung Cancer Susceptibility Study. We evaluated associations between 34,509 genetic variants (3252 genotyped and 31,257 imputed) in 247 genes involved in lymphangiogenesis-related pathway and NSCLC survival. After validation, we finally identified two independent SNPs (SYK rs11787670 A>G and ITGA1 rs67715745 T>C) to be significantly associated with NSCLC overall survival (OS), with adjusted hazards ratios of 0.77 and 0.83 (95% confidence interval =0.66-0.90, P=7.20×10^-4) and 0.84 (95% confidence interval =0.75-0.92, P=3.50×10^-4), respectively. Moreover, an increasing number of combined protective alleles of these two SNPs was significantly associated with an improved NSCLC OS and disease-specific survival (DSS) in the PLCO dataset (P _trend=0.011 and 0.006, respectively). Furthermore, the addition of these protective alleles to the prediction model for the 5-year survival increased the time-dependent area under the curve both from 87% to 87.67% for OS (P=0.029) and from 88.54% to 89.06% for DSS (P=0.022). Subsequent expression quantitative trait loci (eQTL) functional analysis revealed that the rs11787670 G allele was significantly associated with an elevated SYK mRNA expression in normal tissues. Additional analyses suggested a suppressor role for both SYK and ITGA1 in NSCLC survival. Collectively, these findings indicated that SYK rs11787670 A>G and ITGA1 rs67715745 T>C may be independent prognostic factors for NSCLC survival once further validated.

lncRNA NEAT1 mediates sepsis progression by regulating Irak2 via sponging miR-370-3p

Sepsis is a life-threatening condition and often associated with multiple organ failure. Nuclear-enriched abundant transcript 1 (NEAT1), a member of the long non-coding RNAs (lncRNAs), was reported to be involved in the regulation of sepsis progression. However, its precise regulatory mechanism needs to be further explored. In this study, the cell-counting kit-8 assay was used to check cell viability. The quantitative real-time polymerase chain reaction (qRT-PCR) was employed to detect the expression levels of NEAT1, miR-370-3p and Interleukin 1 receptor associated kinase 2 (Irak2). Flow cytometry assay and ELISA were used to check cell apoptosis and the concentrations of inflammatory cytokines, respectively. The starBase was used to predict binding sites between miR-370-3p and NEAT1 or Irak2 and the dual-luciferase reporter assay was performed to verify the interaction. The protein level of Irak2 in samples was measured by western blot. The high concentration of lipopolysaccharide (LPS) led to the high death ratio of RAW 264.7 and HL-1 cells. NEAT1 and Irak2 were upregulated in sepsis tissues and LPS-induced RAW 264.7 and HL-1 cells, opposite to the expression of miR-370-3p. In addition, knockdown of NEAT1 promoted viability, suppressed apoptosis and reduced the expression of inflammatory cytokines in LPS-induced RAW 264.7 and HL-1 cells. Moreover, we found that miR-370-3p interacted with NEAT1 and targeted the 3'UTR of Irak2. Further research indicated that downregulation of miR-370-3p or upregulation of Irak2 rescued NEAT1 silencing-mediated inhibitory effect on sepsis progression. Knockdown of NEAT1 hampered sepsis progression by downregulating Irak2 via interacting with miR-370-3p in LPS-induced RAW 264.7 and HL-1 cells.

Identification of Candidate Biomarkers Correlated With the Pathogenesis and Prognosis of Non-small Cell Lung Cancer via Integrated Bioinformatics Analysis

Background and Objective: Non-small cell lung cancer (NSCLC) accounts for 80-85% of all patients with lung cancer and 5-year relative overall survival (OS) rate is less than 20%, so that identifying novel diagnostic and prognostic biomarkers is urgently demanded. The present study attempted to identify potential key genes associated with the pathogenesis and prognosis of NSCLC. Methods: Four GEO datasets (GSE18842, GSE19804, GSE43458, and GSE62113) were obtained from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) between NSCLC samples and normal ones were analyzed using limma package, and RobustRankAggreg (RRA) package was used to conduct gene integration. Moreover, Search Tool for the Retrieval of Interacting Genes database (STRING), Cytoscape, and Molecular Complex Detection (MCODE) were utilized to establish protein-protein interaction (PPI) network of these DEGs. Furthermore, functional enrichment and pathway enrichment analyses for DEGs were performed by Funrich and OmicShare. While the expressions and prognostic values of top genes were carried out through Gene Expression Profiling Interactive Analysis (GEPIA) and Kaplan Meier-plotter (KM) online dataset. Results: A total of 249 DEGs (113 upregulated and 136 downregulated) were identified after gene integration. Moreover, the PPI network was established with 166 nodes and 1784 protein pairs. Topoisomerase II alpha (TOP2A), a top gene and hub node with higher node degrees in module 1, was significantly enriched in mitotic cell cycle pathway. In addition, Interleukin-6 (IL-6) was enriched in amb2 integrin signaling pathway. The mitotic cell cycle was the most significant pathway in module 1 with the highest P-value. Besides, five hub genes with high degree of connectivity were selected, including TOP2A, CCNB1, CCNA2, UBE2C, and KIF20A, and they were all correlated with worse OS in NSCLC. Conclusion: The results showed that TOP2A, CCNB1, CCNA2, UBE2C, KIF20A, and IL-6 may be potential key genes, while the mitotic cell cycle pathway may be a potential pathway contribute to progression in NSCLC. Further, it could be used as a new biomarker for diagnosis and to direct the synthesis medicine of NSCLC.