Bioinformatic analysis of pivotal genes associated with septic shock

We aimed to identify important genes associated with septic shock and then explore the possibly significant mechanisms of this disease. We downloaded GSE26440 expression data of samples from 98 children with septic shock and 32 normal controls from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) in samples from patients with septic shock were analyzed in comparison with those in samples from normal controls using a limma package. Functional enrichment analysis for DEGs was performed using DAVID, and a protein–protein interaction (PPI) network was constructed. Upstream transcription factors for DEGs were predicted using the CHIPBase database, and a transcriptional regulation network was constructed. A total of 383 significantly DEGs, including 141 downregulated and 242 upregulated genes, were obtained in the sepsis shock group compared with the normal group. The top five nodes in the PPI network were lysine (K)-specific demethylase 6B (KDM6B), histone deacetylase 2 (HDAC2), V-Myc avian myelocytomatosis viral oncogene homolog (MYC), heat-shock protein 90 kDa alpha (cytosolic), class B member 1 (HSP90AB1), and poly (A)-binding protein, cytoplasmic 1 (PABPC1). Nuclear factor kappa-light-chain-enhancer of activated B cells (NFkB) was the transcription factor targeted by most genes, and it regulated the expression of KDM6B, HDAC2, MYC, HSP90AB1, and PABPC1. In conclusion, KDM6B, HDAC2, MYC, HSP90AB1, and PABPC1 may play important roles in the development of septic shock. Furthermore, NFκB may be involved in septic shock by regulating the expression of KDM6B, HDAC2, MYC, HSP90AB1, and PABPC1.