RT Journal Article SR Electronic T1 Genome-wide association study between SARS-CoV-2 single nucleotide polymorphisms and virus copies during infections JF medRxiv FD Cold Spring Harbor Laboratory Press SP 2024.03.02.24303655 DO 10.1101/2024.03.02.24303655 A1 Li, Ke A1 Chaguza, Chrispin A1 Stamp, Julian A1 Chew, Yi Ting A1 Chen, Nicholas F.G. A1 Ferguson, David A1 Pandya, Sameer A1 Kerantzas, Nick A1 Schulz, Wade A1 Yale SARS-CoV-2 Genomic Surveillance Initiative A1 Hahn, Anne M. A1 Ogbunugafor, C Brandon A1 Pitzer, Virginia E. A1 Crawford, Lorin A1 Weinberger, Daniel M. A1 Grubaugh, Nathan D. YR 2024 UL http://medrxiv.org/content/early/2024/07/22/2024.03.02.24303655.abstract AB 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.Competing Interest StatementNDG is a paid consultant for BioNTech, DMW has received consulting fees from Pfizer, Merck, and GSK, unrelated to this manuscript, and has been PI on research grants from Pfizer and Merck to Yale, unrelated to this manuscript.Funding StatementThis project is supported by the Centers for Disease Control and Prevention (CDC) Broad Agency Announcement Contract 75D30122C14697 (NDG). This work does not necessarily represent the views of the CDC.Author DeclarationsI confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.YesThe details of the IRB/oversight body that provided approval or exemption for the research described are given below:The Institutional Review Board from the Yale University Human Research Protection Program determined that the RT-qPCR testing and sequencing of de-identified remnant COVID-19 clinical samples obtained from clinical partners conducted in this study is not research involving human subjects (IRB Protocol ID: 2000028599).I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.YesI understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).YesI have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.YesData and code used in this study are publicly available on Github: https://github.com/grubaughlab/2024_paper_GWAS. All genome sequences used for the GWAS analysis and a subset of the associated metadata (accession number, virus name, collection date, originating lab and submitting lab, and the list of authors) in this dataset are published in GISAIDs EpiCoV database: https://doi.org/10.55876/gis8.240219fh. The de-identified and coded clinical metadata associated with the sequenced samples are available upon request with IRB approval.