Comparative transcriptomic analyzes of human lung epithelial cells infected with wild-type SARS-CoV-2 and its variant with a 12-bp missing in the E gene.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a novel coronavirus that caused a global outbreak of coronavirus disease 2019 (COVID-19) pandemic. To elucidate the mechanism of SARS-CoV-2 replication and immunogenicity, we performed a comparative transcriptome profile of mRNA and long non-coding RNAs (lncRNAs) in human lung epithelial cells infected with the SARS-CoV-2 wild-type strain (8X) and the variant with a 12-bp deletion in the E gene (F8). In total, 3,966 differentially expressed genes (DEGs) and 110 differentially expressed lncRNA (DE-lncRNA) candidates were identified. Of these, 94 DEGs and 32 DE-lncRNAs were found between samples infected with F8 and 8X. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyzes revealed that pathways such as the TNF signaling pathway and viral protein interaction with cytokine and cytokine receptor were involved. Furthermore, we constructed a lncRNA-protein-coding gene co-expression interaction network. The KEGG analysis of the co-expressed genes showed that these differentially expressed lncRNAs were enriched in pathways related to the immune response, which might explain the different replication and immunogenicity properties of the 8X and F8 strains. These results provide a useful resource for studying the pathogenesis of SARS-CoV-2 variants.