Decomposing the sources of SARS-CoV-2 fitness variation in the United States (preprint)

ABSTRACT

The fitness of a pathogen is composite phenotype determined by many different factors influencing growth rates both within and between hosts. Determining what factors shape fitness at the host population-level is especially challenging because both intrinsic factors like pathogen genetics and extrinsic factors such as host behaviour influence between-host transmission potential. These challenges have been highlighted by controversy surrounding the population-level fitness effects of mutations in the SARS-CoV-2 genome and their relative importance when compared against non-genetic factors shaping transmission dynamics. Building upon phylodynamic birth-death models, we develop a new framework to learn how hundreds of genetic and non-genetic factors have shaped the fitness of SARS-CoV-2. We estimate the fitness effects of all amino acid variants and several structural variants that have circulated in the United States between February and September 2020 from viral phylogenies. We also estimate how much fitness variation among pathogen lineages is attributable to genetic versus non-genetic factors such as spatial heterogeneity in transmission rates. Up to September 2020, most fitness variation between lineages can be explained by background spatial heterogeneity in transmission rates across geographic regions. Furthermore, no genetic variant including the Spike D614G mutation has had a significant effect on population-level fitness. Instead, the rapid increase in the frequency of the Spike D614G can be explained by the variant having a spatial transmission advantage due to first establishing in regions with higher transmission rates during the earliest stages of the pandemic.