Background: Tuberculosis (TB) is a global health problem that brings us numerous difficulties. Diverse genetic factors play a significant role in the progress of TB disease. However, still no key genes for TB susceptibility have been reported. This study aimed to identify the key genes of TB through comprehensive bioinformatics analysis.Methods: The series microarray datasets from the gene expression omnibus (GEO) database were analyzed. We used the online tool GEO2R to filtrate differentially expressed genes (DEGs) between TB and health control. Database for annotation can complete gene ontology function analysis as well as Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. Protein-protein interaction (PPI) networks of DEGs were established by STRING online tool and visualized by Cytoscape software. Molecular Complex Detection can complete the analysis of modules in the PPI networks. Finally, the significant hub genes were confirmed by plug-in Genemania of Cytoscape, and verified by the verification cohort and protein test.Results: There are a total of 143 genes were confirmed as DEGs, containing 48 up-regulated genes and 50 down-regulated genes. The gene ontology and Kyoto Encyclopedia of Genes and Genomes analysis show that upregulated DEGs were associated with cancer and phylogenetic, whereas downregulated DEGs mainly concentrate on inflammatory immunity. PPI networks show that signal transducer and activator of transcription 1 (STAT1), guanylate binding protein 5 (GBP5), 2 0 -5 0 -oligoadenylate synthetase 1 (OAS1), catenin beta 1 (CTNNB1), and guanylate binding protein 1 (GBP1) were identified as significantly different hub genes.
Conclusion:We conclude that these genes, including TAT1, GBP5, OAS1, CTNNB1, GBP1 are a candidate as potential core genes in TB and treatment of TB in the future.Abbreviations: CTNNB1 = catenin beta 1, DEG = differentially expressed gene, GBP1 = guanylate binding protein 1, GBP5 = guanylate binding protein 5, GEO = gene expression omnibus, GO = gene ontology, IFI44 = interferon-induced protein 44 like, KEGG = Kyoto Encyclopedia of Genes and Genomes, MTB = Mycobacterium tuberculosis, OAS1 = 2 0 -5 0 -oligoadenylate synthetase 1, OAS2 = 2 0 -5 0 -oligoadenylate synthetase 2, OAS3 = 2 0 -5 0 -oligoadenylate synthetase 3, PARP9 = poly(ADP-ribose) polymerase family member 9, PPI = protein-protein interaction , SOCS3 = suppressor of cytokine signaling 3, STAT1 = signal transducer and activator of transcription 1, TB = tuberculosis, XAF1 = XIAP associated factor 1.