Bioinformatic analysis and preliminary validation of potential therapeutic targets for COVID-19 infection in asthma patients
Authors: Li Yue , Liu Ye , Duo Mengjie , Wu Ruhao , Jiang Tianci , Li Pengfei , Wang Yu , Cheng Zhe 
Published: December 2022

Abstract
Background
Severe acute respiratory syndrome coronavirus 2 causes coronavirus disease 19 (COVID-19). The number of confirmed cases of COVID-19 is also rapidly increasing worldwide, posing a significant challenge to human safety. Asthma is a risk factor for COVID-19, but the underlying molecular mechanisms of the asthma–COVID-19 interaction remain unclear.

Methods
We used transcriptome analysis to discover molecular biomarkers common to asthma and COVID-19. Gene Expression Omnibus database RNA-seq datasets (GSE195599 and GSE196822) were used to identify differentially expressed genes (DEGs) in asthma and COVID-19 patients. After intersecting the differentially expressed mRNAs, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to identify the common pathogenic molecular mechanism. Bioinformatic methods were used to construct protein–protein interaction (PPI) networks and identify key genes from the networks. An online database was used to predict interactions between transcription factors and key genes. The differentially expressed long noncoding RNAs (lncRNAs) in the GSE195599 and GSE196822 datasets were intersected to construct a competing endogenous RNA (ceRNA) regulatory network. Interaction networks were constructed for key genes with RNA-binding proteins (RBPs) and oxidative stress-related proteins. The diagnostic efficacy of key genes in COVID-19 was verified with the GSE171110 dataset. The differential expression of key genes in asthma was verified with the GSE69683 dataset. An asthma cell model was established with interleukins (IL-4, IL-13 and IL-17A) and transfected with siRNA-CXCR1. The role of CXCR1 in asthma development was preliminarily confirmed.

Results
By intersecting the differentially expressed genes for COVID-19 and asthma, 393 common DEGs were obtained. GO and KEGG enrichment analyses of the DEGs showed that they mainly affected inflammation-, cytokine- and immune-related functions and inflammation-related signaling pathways. By analyzing the PPI network, we obtained 10 key genes: TLR4, TLR2, MMP9, EGF, HCK, FCGR2A, SELP, NFKBIA, CXCR1, and SELL. By intersecting the differentially expressed lncRNAs for COVID-19 and asthma, 13 common differentially expressed lncRNAs were obtained. LncRNAs that regulated microRNAs (miRNAs) were mainly concentrated in intercellular signal transduction, apoptosis, immunity and other related functional pathways. The ceRNA network suggested that there were a variety of regulatory miRNAs and lncRNAs upstream of the key genes. The key genes could also bind a variety of RBPs and oxidative stress-related genes. The key genes also had good diagnostic value in the verification set. In the validation set, the expression of key genes was statistically significant in both the COVID-19 group and the asthma group compared with the healthy control group. CXCR1 expression was upregulated in asthma cell models, and interference with CXCR1 expression significantly reduced cell viability.

Conclusions
Key genes may become diagnostic and predictive biomarkers of outcomes in COVID-19 and asthma.