Mathematical modeling of SARS-CoV-2 variant substitutions in European countries: transmission dynamics and epidemiological insights

Abstract

Background: Countries across Europe have faced similar evolutions of SARS-CoV-2 variants of concern, including the Alpha, Delta, and Omicron variants. Materials and methods: We used data from GISAID and applied a robust, automated mathematical substitution model to study the dynamics of COVID-19 variants in Europe over a period of more than 2 years, from late 2020 to early 2023. This model identifies variant substitution patterns and distinguishes between residual and dominant behavior. We used weekly sequencing data from 19 European countries to estimate the increase in transmissibility (Formula presented.) between consecutive SARS-CoV-2 variants. In addition, we focused on large countries with separate regional outbreaks and complex scenarios of multiple competing variants. Results: Our model accurately reproduced the observed substitution patterns between the Alpha, Delta, and Omicron major variants. We estimated the daily variant prevalence and calculated (Formula presented.) between variants, revealing that: ((Formula presented.)) (Formula presented.) increased progressively from the Alpha to the Omicron variant; ((Formula presented.)) (Formula presented.) showed a high degree of variability within Omicron variants; ((Formula presented.)) a higher (Formula presented.) was associated with a later emergence of the variant within a country; ((Formula presented.)) a higher degree of immunization of the population against previous variants was associated with a higher (Formula presented.) for the Delta variant; ((Formula presented.)) larger countries exhibited smaller (Formula presented.) suggesting regionally diverse outbreaks within the same country; and finally ((Formula presented.)) the model reliably captures the dynamics of competing variants, even in complex scenarios. Conclusion: The use of mathematical models allows for precise and reliable estimation of daily cases of each variant. By quantifying (Formula presented.) we have tracked the spread of the different variants across Europe, highlighting a robust increase in transmissibility trend from Alpha to Omicron. Additionally, we have shown that the geographical characteristics of a country, as well as the timing of new variant entrances, can explain some of the observed differences in variant substitution dynamics across countries.