Abstract
COVID-19 pandemic poses a severe threat to public health. However, so far, there are no effective drugs for COVID-19. Transcriptomic changes and key genes related to Th2 cells in COVID-19 have not been reported. These genes play an important role in host interactions with SARS-COV-2 and may be used as promising target. We analyzed five COVID-19-associated GEO datasets (GSE157103, GSE152641, GSE171110, GSE152418, and GSE179627) using the xCell algorithm and weighted gene co-expression network analysis (WGCNA). Results showed that 5 closely correlated modular genes to COVID-19 and Th2 cell enrichment levels, including purple, blue, pink, tan and turquoise, were intersected with differentially expressed genes (DEGs) and 648 shared genes were obtained. GO and KEGG pathway enrichment analyses revealed that they were enriched in cell proliferation, differentiation, and immune responses after virus infection. The most significantly enriched pathway involved the regulation of viral life cycle. Three key genes, namely CCNB1, BUB1, and UBE2C, may clarify the pathogenesis of COVID-19 associated with Th2 cells. 11 drug candidates were identified that could down-regulate three key genes using the cMAP database and demonstrated strong drugs binding energies aganist the three keygenes using molecular docking methods. BUB1, CCNB1 and UBE2C were identified key genes for COVID-19 and could be promising therapeutic targets.