Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico

  1. Hugo G Castelán-Sánchez
  2. Luis Delaye
  3. Rhys PD Inward
  4. Simon Dellicour
  5. Bernardo Gutierrez
  6. Natalia Martinez de la Vina
  7. Celia Boukadida
  8. Oliver Pybus
  9. Guillermo de Anda Jáuregui
  10. Plinio Guzmán
  11. Marisol Flores-Garrido
  12. Óscar Fontanelli
  13. Maribel Hernández Rosales
  14. Amilcar Meneses
  15. Gabriela Olmedo-Alvarez
  16. Alfredo Heriberto Herrera-Estrella
  17. Alejandro Sánchez-Flores
  18. José Esteban Muñoz-Medina
  19. Andreu Comas-García
  20. Bruno Gómez-Gil
  21. Selene Zárate
  22. Blanca Taboada
  23. Susana López
  24. Carlos F Arias
  25. Moritz U G Kraemer
  26. Antonio Lazcano
  27. Marina Escalera Zamudio  Is a corresponding author
  1. Consejo Nacional de Ciencia y Tecnología, Mexico
  2. CINVESTAV-Unidad Irapuato, Mexico
  3. University of Oxford, United Kingdom
  4. KU Leuven, Belgium
  5. Instituto Nacional de Enfermedades Respiratorias, Mexico
  6. Instituto Nacional de Medicina Genómica, Mexico
  7. Astronomer LTD, Mexico
  8. Universidad Nacional Autónoma de México, Mexico
  9. CINVESTAV-IPN, Mexico
  10. Instituto Mexicano del Seguro Social, Mexico
  11. Universidad Autónoma de San Luis Potosí, Mexico
  12. Unidad Regional Mazatlán en Acuicultura y Manejo Ambiental, Mexico
  13. Universidad Autónoma de la Ciudad de México, Mexico
  14. Universidad Nacional Autónoma de Méxic, Mexico

Abstract

Over 200 different SARS-CoV-2 lineages have been observed in Mexico by November 2021. To investigate lineage replacement dynamics, we applied a phylodynamic approach and explored the evolutionary trajectories of five dominant lineages that circulated during the first year of local transmission. For most lineages, peaks in sampling frequencies coincided with different epidemiological waves of infection in Mexico. Lineages B.1.1.222 and B.1.1.519 exhibited similar dynamics, constituting clades that likely originated in Mexico and persisted for >12 months. Lineages B.1.1.7, P.1 and B.1.617.2 also displayed similar dynamics, characterized by multiple introduction events leading to a few successful extended local transmission chains that persisted for several months. For the largest B.1.617.2 clades, we further explored viral lineage movements across Mexico. Many clades were located within the south region of the country, suggesting that this area played a key role in the spread of SARS-CoV-2 in Mexico.

Data availability

Virus genome IDs and GISAID accession numbers for the sequences used in each dataset are provided in the Supplementary file 1 file. All genomic and epidemiological data supporting the findings of this study is publicly available from GISAID/GenBank, from the Ministry Of Health Mexico102, and/or from the 'Our World in Data' coronavirus pandemic web portal 29. For the GISAID data used, the corresponding acknowledgement table is available on the 'GISAID Data Acknowledgement Locator' under the EPI_SET_20220405qd and EPI_SET_20220215at keys 49. Our bioinformatic pipeline implementing a migration data and phylogenetically-informed sequence subsampling approach is publicly available at https://github.com/rhysinward/Mexico_subsampling.

Article and author information

Author details

  1. Hugo G Castelán-Sánchez

    Consejo Nacional de Ciencia y Tecnología, Mexico City, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  2. Luis Delaye

    Departamento de Ingeniería Genética, CINVESTAV-Unidad Irapuato, Guanajuato, Mexico
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4193-2720
  3. Rhys PD Inward

    Department of Biology, University of Oxford, Oxford, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0016-661X
  4. Simon Dellicour

    Department of Microbiology, Immunology and Transplantation, KU Leuven, Leuven, Belgium
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-9558-1052
  5. Bernardo Gutierrez

    Department of Biology, University of Oxford, Oxford, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
  6. Natalia Martinez de la Vina

    Department of Biology, University of Oxford, Oxford, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
  7. Celia Boukadida

    Centro de Investigación en Enfermedades Infecciosas, Instituto Nacional de Enfermedades Respiratorias, Mexico City, Mexico
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1744-0083
  8. Oliver Pybus

    Department of Biology, University of Oxford, Oxford, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
  9. Guillermo de Anda Jáuregui

    Instituto Nacional de Medicina Genómica, Mexico City, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  10. Plinio Guzmán

    Astronomer LTD, Mexico City, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  11. Marisol Flores-Garrido

    Escuela Nacional de Estudios Superiores, Universidad Nacional Autónoma de México, Morelia, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  12. Óscar Fontanelli

    Departamento de Ingeniería Genética, CINVESTAV-Unidad Irapuato, Guanajuato, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  13. Maribel Hernández Rosales

    Departamento de Ingeniería Genética, CINVESTAV-Unidad Irapuato, Guanajuato, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  14. Amilcar Meneses

    Departamento de Ciencias de la Computación, CINVESTAV-IPN, Mexico City, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  15. Gabriela Olmedo-Alvarez

    Departamento de Ingeniería Genética, CINVESTAV-Unidad Irapuato, Guanajuato, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  16. Alfredo Heriberto Herrera-Estrella

    Laboratorio de expresión génica y desarrollo en hongos, CINVESTAV-Unidad Irapuato, Irapuato, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  17. Alejandro Sánchez-Flores

    Instituto de Biotecnología, Universidad Nacional Autónoma de México, Mexico City, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  18. José Esteban Muñoz-Medina

    Coordinación de Calidad de Insumos y Laboratorios Especializados, Instituto Mexicano del Seguro Social, Mexico City, Mexico
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1289-4457
  19. Andreu Comas-García

    Facultad de Medicina y Centro de Investigación en Ciencias de la Salud y Biomedicina, Universidad Autónoma de San Luis Potosí, San Luis Potosí, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  20. Bruno Gómez-Gil

    Centro de Investigación en Alimentación y Desarrollo-CIAD, Unidad Regional Mazatlán en Acuicultura y Manejo Ambiental, Sinaloa, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  21. Selene Zárate

    Posgrado en Ciencias Genómicas, Universidad Autónoma de la Ciudad de México, Mexico City, Mexico
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1034-204X
  22. Blanca Taboada

    Departamento de Genética del Desarrollo y Fisiología Molecular, Universidad Nacional Autónoma de México, Mexico City, Mexico
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1896-5962
  23. Susana López

    Departamento de Genética del Desarrollo y Fisiología Molecular, Universidad Nacional Autónoma de México, Mexico City, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  24. Carlos F Arias

    Departamento de Genética del Desarrollo y Fisiología Molecular, Universidad Nacional Autónoma de México, Cuernavaca, Mexico
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3130-4501
  25. Moritz U G Kraemer

    Department of Biology, University of Oxford, Oxford, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
  26. Antonio Lazcano

    Facultad de Ciencias, Universidad Nacional Autónoma de Méxic, Mexico City, Mexico
    Competing interests
    The authors declare that no competing interests exist.
  27. Marina Escalera Zamudio

    Department of Biology, University of Oxford, Oxford, United Kingdom
    For correspondence
    marina.escalerazamudio@zoo.ox.ac.uk
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-4773-2773

Funding

FNRS (F.4515.22)

  • Simon Dellicour

UNAM (DGAPA-PAPIIT (IN214421)

  • Antonio Lazcano

UNAM (DGAPA-PAPIME (PE204921))

  • Antonio Lazcano

Research Foundation Flanders (G098321N)

  • Simon Dellicour

European Horizon 2020 project MOOD (874850)

  • Simon Dellicour

Leverhulme Trust (ECF-2019-542)

  • Marina Escalera Zamudio

European Horizon 2020 project MOOD (874850)

  • Oliver Pybus

European Horizon 2020 project MOOD (874850)

  • Moritz U G Kraemer

CONACyT Vigilancia Genómica del Virus SARS-CoV-2 en México-2022"" (PP-F003)

  • Carlos F Arias

Ministry of Education, Science, Technology and Innovation of Mexico City (057)

  • Carlos F Arias

AHF Global Public Health Institute at the University of Miami (Genomic surveillance for SARS-CoV-2 variants in Mexico"")

  • Carlos F Arias

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Copyright

© 2023, Castelán-Sánchez et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,191
    views
  • 223
    downloads
  • 7
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Hugo G Castelán-Sánchez
  2. Luis Delaye
  3. Rhys PD Inward
  4. Simon Dellicour
  5. Bernardo Gutierrez
  6. Natalia Martinez de la Vina
  7. Celia Boukadida
  8. Oliver Pybus
  9. Guillermo de Anda Jáuregui
  10. Plinio Guzmán
  11. Marisol Flores-Garrido
  12. Óscar Fontanelli
  13. Maribel Hernández Rosales
  14. Amilcar Meneses
  15. Gabriela Olmedo-Alvarez
  16. Alfredo Heriberto Herrera-Estrella
  17. Alejandro Sánchez-Flores
  18. José Esteban Muñoz-Medina
  19. Andreu Comas-García
  20. Bruno Gómez-Gil
  21. Selene Zárate
  22. Blanca Taboada
  23. Susana López
  24. Carlos F Arias
  25. Moritz U G Kraemer
  26. Antonio Lazcano
  27. Marina Escalera Zamudio
(2023)
Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico
eLife 12:e82069.
https://doi.org/10.7554/eLife.82069

Share this article

https://doi.org/10.7554/eLife.82069

Further reading

    1. Epidemiology and Global Health
    2. Microbiology and Infectious Disease
    Gillian AM Tarr, Linda Chui ... Tim A McAllister
    Research Article

    Several areas of the world suffer a notably high incidence of Shiga toxin-producing Escherichia coli. To assess the impact of persistent cross-species transmission systems on the epidemiology of E. coli O157:H7 in Alberta, Canada, we sequenced and assembled E. coli O157:H7 isolates originating from collocated cattle and human populations, 2007–2015. We constructed a timed phylogeny using BEAST2 using a structured coalescent model. We then extended the tree with human isolates through 2019 to assess the long-term disease impact of locally persistent lineages. During 2007–2015, we estimated that 88.5% of human lineages arose from cattle lineages. We identified 11 persistent lineages local to Alberta, which were associated with 38.0% (95% CI 29.3%, 47.3%) of human isolates. During the later period, six locally persistent lineages continued to be associated with human illness, including 74.7% (95% CI 68.3%, 80.3%) of reported cases in 2018 and 2019. Our study identified multiple locally evolving lineages transmitted between cattle and humans persistently associated with E. coli O157:H7 illnesses for up to 13 y. Locally persistent lineages may be a principal cause of the high incidence of E. coli O157:H7 in locations such as Alberta and provide opportunities for focused control efforts.

    1. Epidemiology and Global Health
    Riccardo Spott, Mathias W Pletz ... Christian Brandt
    Research Article

    Given the rapid cross-country spread of SARS-CoV-2 and the resulting difficulty in tracking lineage spread, we investigated the potential of combining mobile service data and fine-granular metadata (such as postal codes and genomic data) to advance integrated genomic surveillance of the pandemic in the federal state of Thuringia, Germany. We sequenced over 6500 SARS-CoV-2 Alpha genomes (B.1.1.7) across 7 months within Thuringia while collecting patients’ isolation dates and postal codes. Our dataset is complemented by over 66,000 publicly available German Alpha genomes and mobile service data for Thuringia. We identified the existence and spread of nine persistent mutation variants within the Alpha lineage, seven of which formed separate phylogenetic clusters with different spreading patterns in Thuringia. The remaining two are subclusters. Mobile service data can indicate these clusters’ spread and highlight a potential sampling bias, especially of low-prevalence variants. Thereby, mobile service data can be used either retrospectively to assess surveillance coverage and efficiency from already collected data or to actively guide part of a surveillance sampling process to districts where these variants are expected to emerge. The latter concept was successfully implemented as a proof-of-concept for a mobility-guided sampling strategy in response to the surveillance of Omicron sublineage BQ.1.1. The combination of mobile service data and SARS-CoV-2 surveillance by genome sequencing is a valuable tool for more targeted and responsive surveillance.