Integrated Analysis of Bulk RNA-Seq and Single-Cell RNA-Seq Unravels the Influences of SARS-CoV-2 Infections to Cancer Patients

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a highly contagious and pathogenic coronavirus that emerged in late 2019 and caused a pandemic of respiratory illness termed as coronavirus disease 2019 (COVID-19). Cancer patients are more susceptible to SARS-CoV-2 infection. The treatment of cancer patients infected with SARS-CoV-2 is more complicated, and the patients are at risk of poor prognosis compared to other populations. Patients infected with SARS-CoV-2 are prone to rapid development of acute respiratory distress syndrome (ARDS) of which pulmonary fibrosis (PF) is considered a sequelae. Both ARDS and PF are factors that contribute to poor prognosis in COVID-19 patients. However, the molecular mechanisms among COVID-19, ARDS and PF in COVID-19 patients with cancer are not well-understood. In this study, the common differentially expressed genes (DEGs) between COVID-19 patients with and without cancer were identified. Based on the common DEGs, a series of analyses were performed, including Gene Ontology (GO) and pathway analysis, protein-protein interaction (PPI) network construction and hub gene extraction, transcription factor (TF)-DEG regulatory network construction, TF-DEG-miRNA coregulatory network construction and drug molecule identification. The candidate drug molecules (e.g., Tamibarotene CTD 00002527) obtained by this study might be helpful for effective therapeutic targets in COVID-19 patients with cancer. In addition, the common DEGs among ARDS, PF and COVID-19 patients with and without cancer are TNFSF10 and IFITM2. These two genes may serve as potential therapeutic targets in the treatment of COVID-19 patients with cancer. Changes in the expression levels of TNFSF10 and IFITM2 in CD14+/CD16+ monocytes may affect the immune response of COVID-19 patients. Specifically, changes in the expression level of TNFSF10 in monocytes can be considered as an immune signature in COVID-19 patients with hematologic cancer. Targeting N-6-methyladenosine (m6A) pathways (e.g., METTL3/SERPINA1 axis) to restrict SARS-CoV-2 reproduction has therapeutic potential for COVID-19 patients.