SARS-CoV-2 infections are highly correlated with the overexpression of pro-inflammatory cytokines in what is known as a cytokine storm, leading to high fatality rates. Such infections are accompanied by SIRS, ARDS, and sepsis, suggesting a potential link between the three phenotypes. Currently, little is known about the transcriptional similarity between these conditions. Herein, weighted gene co-expression network analysis (WGCNA) clustering was applied to RNA-seq datasets (GSE147902, GSE66890, GSE74224, GSE177477) to identify modules of highly co-expressed and correlated genes, cross referenced with dataset GSE160163, across the samples. To assess the transcriptome similarities between the conditions, module preservation analysis was performed and functional enrichment was analyzed in DAVID webserver. The hub genes of significantly preserved modules were identified, classified into upregulated or downregulated, and used to screen candidate drugs using Connectivity Map (CMap) to identify repurposed drugs. Results show that several immune pathways (chemokine signaling, NOD-like signaling, and Th1 and Th2 cell differentiation) are conserved across the four diseases. Hub genes screened using intramodular connectivity show significant relevance with the pathogenesis of cytokine storms. Transcriptomic-driven drug repurposing identified seven candidate drugs (SB-202190, eicosatetraenoic-acid, loratadine, TPCA-1, pinocembrin, mepacrine, and CAY-10470) that targeted several immune-related processes. These identified drugs warrant further study into their efficacy for treating cytokine storms, and in vitro and in vivo experiments are recommended to confirm the findings of this study.
Keywords: bioinformatics; computational biology; cytokine storm; drug repurposing; gene expression; genetics; systems biology.