Human phospho-signaling networks of SARS-CoV-2 infection are rewired by population genetic variants

SARS-CoV-2 infection hijacks signaling pathways and induces protein-protein interactions between human and viral proteins. Human genetic variation may impact SARS-CoV-2 infection and COVID-19 pathology; however, the role of genetic variation in these signaling networks remains uncharacterized. We studied human single nucleotide variants (SNVs) affecting phosphorylation sites modulated by SARS-CoV-2 infection, using machine learning to identify amino acid changes altering kinase-bound sequence motifs. We found 2033 infrequent phosphorylation-associated SNVs (pSNVs) that are enriched in sequence motif alterations, potentially reflecting the evolution of signaling networks regulating host defenses. Proteins with pSNVs are involved in viral life cycle processes and host responses, including regulators of RNA splicing and interferon response, as well as glucose homeostasis pathways with potential associations with COVID-19 co-morbidities. Certain pSNVs disrupt CDK and MAPK substrate motifs and replace these with motifs recognized by Tank Binding Kinase 1 (TBK1) involved in innate immune responses, indicating consistent rewiring of infection signaling networks. Our analysis highlights potential genetic factors contributing to the variation of SARS-CoV-2 infection and COVID-19 and suggests leads for mechanistic and translational studies.