Abstract:Due to the emergence of new variants of the SARS-CoV-2 coronavirus, the question of how the viral genomes evolved, leading to the formation of highly infectious strains, becomes particularly important. Three major emergent strains, Alpha, Beta and Delta, characterized by a significant number of missense mutations, provide a natural test field. We accumulated and aligned 4.7 million SARS-CoV-2 genomes from the GISAID database and carried out a comprehensive set of analyses. This collection covers the period until the end of October 2021, i.e., the beginnings of the Omicron variant. First, we explored combinatorial complexity of the genomic variants emerging and their timing, indicating very strong, albeit hidden, selection forces. Our analyses show that the mutations that define variants of concern did not arise gradually but rather co-evolved rapidly, leading to the emergence of the full variant strain. To explore in more detail the evolutionary forces at work, we developed time trajectories of mutations at all 29,903 sites of the SARS-CoV-2 genome, week by week, and stratified them into trends related to (i) point substitutions, (ii) deletions and (iii) non-sequenceable regions. We focused on classifying the genetic forces active at different ranges of the mutational spectrum. We observed the agreement of the lowest-frequency mutation spectrum with the Griffiths-Tavaré theory, under the Infinite Sites Model and neutrality. If we widen the frequency range, we observe the site frequency spectra much more consistently with the Tung-Durrett model assuming clone competition and selection. The coefficients of the fitting model indicate the possibility of selection acting to promote gradual growth slowdown, as observed in the history of the variants of concern. These results add up to a model of genomic evolution, which partly fits into the classical drift barrier ideas. Certain observations, such as mutation "bands" persistent over the epidemic history, suggest contribution of genetic forces different from mutation, drift and selection, including recombination or other genome transformations. In addition, we show that a "toy" mathematical model can qualitatively reproduce how new variants (clones) stem from rare advantageous driver mutations, and then acquire neutral or disadvantageous passenger mutations which gradually reduce their fitness so they can be then outcompeted by new variants due to other driver mutations.

Early evolution of BA.2.86 sheds light on the origins of highly divergent SARS-CoV-2 lineages

Distinct evolution of SARS-CoV-2 Omicron XBB and BA.2.86/JN.1 lineages combining increased fitness and antibody evasion

SARS-CoV-2 evolution in the Omicron era

Evolution of a Distinct SARS-CoV-2 Lineage Identified during an Investigation of a Hospital Outbreak

SARS-CoV-2 evolution from the BA.2.86 to JN.1 variants: unexpected consequences

Fast Evolution of SARS-CoV-2 BA.2.86 to JN.1 under Heavy Immune Pressure

Evolutionary origins of the SARS-CoV-2 sarbecovirus lineage responsible for the COVID-19 pandemic

On the Origin and Continuing Evolution of SARS-CoV-2

Genomic Analysis of SARS-CoV-2 Alpha, Beta and Delta Variants of Concern Uncovers Signatures of Neutral and Non-Neutral Evolution

Intrahost evolution leading to distinct lineages in the upper and lower respiratory tracts during SARS-CoV-2 prolonged infection

SARS-CoV-2 infectivity and antigenic evasion: spotlight on isolated Omicron sub-lineages

Unraveling the genetic evolution of SARS-CoV-2 Recombinants using mutational dynamics across the different lineages

Evolutionary Trajectory for the Emergence of Novel Coronavirus SARS-CoV-2

The origins and molecular evolution of SARS-CoV-2 lineage B. 1.1. 7 in the UK

Tracking SARS-CoV-2 Omicron diverse spike gene mutations identifies multiple inter-variant recombination events

Generation and transmission of interlineage recombinants in the SARS-CoV-2 pandemic

Immune escape and attenuated severity associated with the SARS-CoV-2 BA.2.86/JN.1 lineage

Distinct SARS-CoV-2 Populational Immune Backgrounds Tolerate Divergent RBD Evolutionary Preferences

Molecular and structural insights into SARS-CoV-2 evolution: from BA.2 to XBB subvariants

Bayesian Phylodynamic Inferences on the Temporal Evolution and Global Transmission of SARS-CoV-2