Analysis of 6.4 million SARS-CoV-2 genomes identifies mutations associated with fitness

Repeated emergence of SARS-CoV-2 variants with increased fitness underscores the value of rapid detection and characterization of new lineages. We have developed PyR₀, a hierarchical Bayesian multinomial logistic regression model that infers the relative prevalence of all viral lineages across geographic regions, detects lineages increasing in prevalence, and identifies mutations relevant to fitness. Applying PyR₀ to all publicly available SARS-CoV-2 genomes, we identify numerous substitutions that increase fitness, including previously identified spike mutations and many non-spike mutations within the nucleocapsid and nonstructural proteins. PyR₀ forecasts the growth of new lineages from their mutational profile ranks the fitness of lineages as new sequences become available, and prioritizes mutations of biological and public health concerns for functional characterization.