Li, Ke, Chaguza, Chrispin, Stamp, Julian, Chew, Yi Ting, Chen, Nicholas F. G., Ferguson, David, Pandya, Sameer, Kerantzas, Nick, Schulz, Wade, Hahn, Anne M., Ogbunugafor, C. Brandon, Pitzer, Virginia E., Crawford, Lorin, Weinberger, Daniel M., and Grubaugh, Nathan D.
Significant variations have been observed in viral copies generated during SARS-CoV-2 infections. However, the factors that impact viral copies and infection dynamics are not fully understood, and may be inherently dependent upon different viral and host factors. Here, we conducted virus whole genome sequencing and measured viral copies using RT-qPCR from 9,902 SARS-CoV-2 infections over a 2-year period to examine the impact of virus genetic variation on changes in viral copies adjusted for host age and vaccination status. Using a genome-wide association study (GWAS) approach, we identified multiple single-nucleotide polymorphisms (SNPs) corresponding to amino acid changes in the SARS-CoV-2 genome associated with variations in viral copies. We further applied a marginal epistasis test to detect interactions among SNPs and identified multiple pairs of substitutions located in the spike gene that have non-linear effects on viral copies. We also analyzed the temporal patterns and found that SNPs associated with increased viral copies were predominantly observed in Delta and Omicron BA.2/BA.4/BA.5/XBB infections, whereas those associated with decreased viral copies were only observed in infections with Omicron BA.1 variants. Our work showcases how GWAS can be a useful tool for probing phenotypes related to SNPs in viral genomes that are worth further exploration. We argue that this approach can be used more broadly across pathogens to characterize emerging variants and monitor therapeutic interventions. Author summary: Our study explores why viral load (copies measured by RT-qPCR) varies during SARS-CoV-2 infections by analyzing viral mutations and measuring viral copies in 9,902 individuals over two years. We aimed to understand how genetic differences in SARS-CoV-2 influence viral copies, considering host age and vaccination status. Using a genome-wide association study (GWAS), we identified several single-nucleotide polymorphisms (SNPs) in the virus linked to variations in viral levels. Notably, interactions between certain SNPs in the spike gene had non-linear effects on viral copies. Our analysis revealed that SNPs associated with higher viral copies were common in Delta and Omicron BA.2/BA.4/BA.5/XBB variants, while those linked to lower levels were mainly found in Omicron BA.1. This research highlights GWAS as a powerful tool for exploring virus genetics and suggests it can be broadly applied to monitor new variants of COVID-19 and other infectious diseases. [ABSTRACT FROM AUTHOR]