Identification of the feature genes involved in cytokine release syndrome in COVID-19.

Abstract

Screening of feature genes involved in cytokine release syndrome (CRS) from the coronavirus disease 19 (COVID-19). The data sets related to COVID-19 were retrieved using Gene Expression Omnibus (GEO) database, the differentially expressed genes (DEGs) related to CRS were analyzed with R software and Venn diagram, and the biological processes and signaling pathways involved in DEGs were analyzed with GO and KEGG enrichment. Core genes were screened using Betweenness and MCC algorithms. GSE164805 and GSE171110 dataset were used to verify the expression level of core genes. Immunoinfiltration analysis was performed by ssGSEA algorithm in the GSVA package. The DrugBank database was used to analyze the feature genes for potential therapeutic drugs. This study obtained 6950 DEGs, of which 971 corresponded with CRS disease genes (common genes). GO and KEGG enrichment showed that multiple biological processes and signaling pathways associated with common genes were closely related to the inflammatory response. Furthermore, the analysis revealed that transcription factors that regulate these common genes are also involved in inflammatory response. Betweenness and MCC algorithms were used for common gene screening, yielding seven key genes. GSE164805 and GSE171110 dataset validation revealed significant differences between the COVID-19 and normal controls in four core genes (feature genes), namely IL6R, TLR4, TLR2, and IFNG. The upregulated IL6R, TLR4, and TLR2 genes were mainly involved in the Toll-like receptor signaling pathway of the inflammatory pathway, while the downregulated IFNG genes primarily participated in the necroptosis and JAK-STAT signaling pathways. Moreover, immune infiltration analysis indicated that higher expression of these genes was associated with immune cell infiltration that mediates inflammatory response. In addition, potential therapeutic drugs for these four feature genes were identified via the DrugBank database. IL6R, TLR4, TLR2, and IFNG may be potential pathogenic genes and therapeutic targets for the CRS associated with COVID-19.