Regional epidemic dynamics and Delta variant diversity resulted in varying rates of spread of Omicron-BA.1 in Mexico
Authors: Zárate Selene , Taboada Blanca , Rosales-Rivera Mauricio , García-López Rodrigo , Muñoz-Medina José Esteban , Sanchez-Flores Alejandro , Herrera-Estrella Alfredo , Gómez-Gil Bruno , Mojica Nelly Selem , Salas-Lais Angel Gustavo , Vazquez-Perez Joel Armando , Cabrera-Gaytán David Alejandro , Fernandes-Matano Larissa , Uribe-Noguez Luis Antonio , Chale-Dzul Juan Bautista , Meza Brenda Irasema Maldonado , Mejía-Nepomuceno Fidencio , Pérez-Padilla Rogelio , Gutiérrez-Ríos Rosa María , Loza Antonio , López Susana , Arias Carlos F. 
Published: October 2022

The Omicron subvariant BA.1 of SARS-CoV-2 was first detected in November 2021 and quickly spread worldwide, displacing the Delta variant. In Mexico, this subvariant began spreading during the first week of December 2021 and became dominant in the next three weeks, causing the fourth COVID-19 epidemiological surge in the country. Unlike previous SARS-CoV-2 variants, BA.1 did not acquire local substitutions nor exhibited a geographically distinct circulation pattern in Mexico. However, a regional difference in the speed of the replacement of the Delta variant was observed, as some northern states showed persistence of Delta lineages well into February 2022. Mexican states were divided into four regions (North, Central North, Central South, and Southeast) based on the lineage circulation before the dominance of BA.1 to study possible causes for this difference. For each region, the time to fixation of BA.1, the diversity of Delta sublineages in the weeks preceding BA.1 entry, the population density, and the level of virus circulation during the inter-wave interval were determined. An association between a faster Omicron spread and lower Delta diversity, as well as fewer COVID-19 cases during the Delta-BA.1.x inter-wave period, was observed. For example, the North region exhibited the slowest spread but had the highest diversity of Delta sublineages and the greatest number of inter-wave cases relative to the maximum amount of the virus circulating in the region, whereas the Southeast region showed the opposite. Viral diversity and the relative abundance of the virus in a particular area around the time of the introduction of a new lineage seem to have influenced the spread dynamics. Nonetheless, if there is a significant difference in the fitness of the variants or the time allowed for the competition is sufficient, it seems the fitter virus will eventually become dominant, as observed in the eventual dominance of the BA.1.x variant in Mexico.

Impact statement

The surveillance of lineage circulation of SARS-CoV-2 has helped identify variants that have a transmission advantage and are of concern to public health and to track the virus dispersion accurately. However, many factors contributing to differences in lineage spread dynamics beyond the acquisition of specific mutations remain poorly understood. In this work, a description of BA.1 entry and dispersion within Mexico is presented, and which factors potentially affected the spread rates of the Omicron variant BA.1 among geographical regions in the country are analyzed, underlining the importance of population density, the proportion of active cases, and viral lineage diversity and identity before the entry of BA.1.

Data summary

This work was carried out using data shared through the GISAID initiative. All sequences and metadate are available through GISAID with the accession EPI_SET_220927gw, accession numbers and metadata are also reported in the supplemental material of this article. Epidemiological data was obtained though the Secretaría de Salud website (https://www.gob.mx/salud/documentos/datos-abiertos-152127),