Identification of receptors and factors associated with human coronaviruses in the oral cavity using single-cell RNA sequencing

Abstract

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) ravaged the world, and Coronavirus Disease 2019 (COVID-19) exhibited highly prevalent oral symptoms that had significantly impacted the lives of affected patients. However, the involvement of four human coronavirus (HCoVs), namely SARS-CoV-2, SARS-CoV, MERS-CoV, and HCoV-229E, in oral cavity infections remained poorly understood. We integrated single-cell RNA sequencing (scRNA-seq) data of seven human oral tissues through consistent normalization procedure, including minor salivary gland (MSG), parotid gland (PG), tongue, gingiva, buccal, periodontium and pulp. The Seurat, scDblFinder, Harmony, SingleR, Ucell and scCancer packages were comprehensively used for analysis. We identified specific cell clusters and generated expression profiles of SARS-CoV-2 and coronavirus-associated receptors and factors (SCARFs) in seven oral regions, providing direction for predicting the tropism of four HCoVs for oral tissues, as well as for dental clinical treatment. Based on our analysis, it appears that various SCARFs, including ACE2, ASGR1, KREMEN1, DPP4, ANPEP, CD209, CLEC4G/M, TMPRSS family proteins (including TMPRSS2, TMPRSS4, and TMPRSS11A), and FURIN, are expressed at low levels in the oral cavity. Conversely, BSG, CTSB, and CTSL exhibit enrichment in oral tissues. Our study also demonstrates widespread expression of restriction factors, particularly IFITM1-3 and LY6E, in oral cells. Additionally, some replication, assembly, and trafficking factors appear to exhibit broad oral tissues expression patterns. Overall, the oral cavity could potentially serve as a high-risk site for SARS-CoV-2 infection, while displaying a comparatively lower degree of susceptibility towards other HCoVs (including SARS-CoV, MERS-CoV and HCoV-229E). Specifically, MSG, tongue, and gingiva represent potential sites of vulnerability for four HCoVs infection, with the MSG exhibiting a particularly high susceptibility. However, the expression patterns of SCARFs in other oral sites demonstrate relatively intricate and may only be specifically associated with SARS-CoV-2 infection. Our study sheds light on the mechanisms of HCoVs infection in the oral cavity as well as gains insight into the characteristics and distribution of possible HCoVs target cells in oral tissues, providing potential therapeutic targets for HCoVs infection in the oral cavity.