Abstract
Coronavirus disease 2019 (COVID-19) is acutely infectious pneumonia. Currently, the specific causes and treatment targets of COVID-19 are still unclear. Herein, comprehensive bioinformatics methods were employed to analyze the hub genes in COVID-19 and tried to reveal its potential mechanisms. First of all, 34 groups of COVID-19 lung tissues and 17 other diseases' lung tissues were selected from the GSE151764 gene expression profile for research. According to the analysis of the DEGs (differentially expressed genes) in the samples using the limma software package, 84 upregulated DEGs and 46 downregulated DEGs were obtained. Later, by the Database for Annotation, Visualization, and Integrated Discovery (DAVID), they were enriched in the Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. It was found that the upregulated DEGs were enriched in the type I interferon signaling pathway, AGE-RAGE signaling pathway in diabetic complications, coronavirus disease, etc. Downregulated DEGs were in cellular response to cytokine stimulus, IL-17 signaling pathway, FoxO signaling pathway, etc. Then, based on GSEA, the enrichment of the gene set in the sample was analyzed in the GO terms, and the gene set was enriched in the positive regulation of myeloid leukocyte cytokine production involved in immune response, programmed necrotic cell death, translesion synthesis, necroptotic process, and condensed nuclear chromosome. Finally, with the help of STRING tools, the PPI (protein-protein interaction) network diagrams of DEGs were constructed. With degree ≥13 as the cutoff degree, 3 upregulated hub genes (ISG15, FN1, and HLA-G) and 4 downregulated hub genes (FOXP3, CXCR4, MMP9, and CD69) were screened out for high degree. All these findings will help us to understand the potential molecular mechanisms of COVID-19, which is also of great significance for its diagnosis and prevention.