Abstract
In the post COVID-19 era, new SARS-CoV-2 variant strains may continue emerging and long COVID is poised to be another public health challenge. Deciphering the molecular susceptibility of receptors to SARS-CoV-2 spike protein is critical for understanding the immune responses in COVID-19 and the rationale of multi-organ injuries. Currently, such systematic exploration remains limited. Here, we conduct multi-omic analysis of protein binding affinities, transcriptomic expressions, and single-cell atlases to characterize the molecular susceptibility of receptors to SARS-CoV-2 spike protein. Initial affinity analysis explains the domination of delta and omicron variants and demonstrates the strongest affinities between BSG (CD147) receptor and most variants. Further transcriptomic data analysis on 4100 experimental samples and single-cell atlases of 1.4 million cells suggest the potential involvement of BSG in multi-organ injuries and long COVID, and explain the high prevalence of COVID-19 in elders as well as the different risks for patients with underlying diseases. Correlation analysis validated moderate associations between BSG and viral RNA abundance in multiple cell types. Moreover, similar patterns were observed in primates and validated in proteomic expressions. Overall, our findings implicate important therapeutic targets for the development of receptor-specific vaccines and drugs for COVID-19.