Blood transcriptome analysis revealed the crosstalk between COVID-19 and HIV

BACKGROUND: The severe coronavirus disease 2019 (COVID-19) is an infectious disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which has resulted in the most devastating pandemic in modern history. Human immunodeficiency virus (HIV) destroys immune system cells and weakens the body’s ability to resist daily infections and diseases. Furthermore, HIV-infected individuals had double COVID-19 mortality risk and experienced worse COVID-related outcomes. However, the existing research still lacks the understanding of the molecular mechanism underlying crosstalk between COVID-19 and HIV. The aim of our work was to illustrate blood transcriptome crosstalk between COVID-19 and HIV and to provide potential drugs that might be useful for the treatment of HIV-infected COVID-19 patients. METHODS: COVID-19 datasets (GSE171110 and GSE152418) were downloaded from Gene Expression Omnibus (GEO) database, including 54 whole-blood samples and 33 peripheral blood mononuclear cells samples, respectively. HIV dataset (GSE37250) was also obtained from GEO database, containing 537 whole-blood samples. Next, the “Deseq2” package was used to identify differentially expressed genes (DEGs) between COVID-19 datasets (GSE171110 and GSE152418) and the “limma” package was utilized to identify DEGs between HIV dataset (GSE37250). By intersecting these two DEG sets, we generated common DEGs for further analysis, containing Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Ontology (GO) functional enrichment analysis, protein-protein interaction (PPI) analysis, transcription factor (TF) candidate identification, microRNAs (miRNAs) candidate identification and drug candidate identification. RESULTS: In this study, a total of 3213 DEGs were identified from the merged COVID-19 dataset (GSE171110 and GSE152418), and 1718 DEGs were obtained from GSE37250 dataset. Then, we identified 394 common DEGs from the intersection of the DEGs in COVID-19 and HIV datasets. GO and KEGG enrichment analysis indicated that common DEGs were mainly gathered in chromosome-related and cell cycle-related signal pathways. Top ten hub genes (CCNA2, CCNB1, CDC20, TOP2A, AURKB, PLK1, BUB1B, KIF11, DLGAP5, RRM2) were ranked according to their scores, which were screened out using degree algorithm on the basis of common DEGs. Moreover, top ten drug candidates (LUCANTHONE, Dasatinib, etoposide, Enterolactone, troglitazone, testosterone, estradiol, calcitriol, resveratrol, tetradioxin) ranked by their P values were screened out, which maybe be beneficial for the treatment of HIV-infected COVID-19 patients. CONCLUSION: In this study, we provide potential molecular targets, signaling pathways, small molecular compounds, and promising biomarkers that contribute to worse COVID-19 prognosis in patients with HIV, which might contribute to precise diagnosis and treatment for HIV-infected COVID-19 patients.