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,207
    views
  • 226
    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
    Bo Zheng, Bronner P Gonçalves ... Caoyi Xue
    Research Article

    Background:

    In many settings, a large fraction of the population has both been vaccinated against and infected by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Hence, quantifying the protection provided by post-infection vaccination has become critical for policy. We aimed to estimate the protective effect against SARS-CoV-2 reinfection of an additional vaccine dose after an initial Omicron variant infection.

    Methods:

    We report a retrospective, population-based cohort study performed in Shanghai, China, using electronic databases with information on SARS-CoV-2 infections and vaccination history. We compared reinfection incidence by post-infection vaccination status in individuals initially infected during the April–May 2022 Omicron variant surge in Shanghai and who had been vaccinated before that period. Cox models were fit to estimate adjusted hazard ratios (aHRs).

    Results:

    275,896 individuals were diagnosed with real-time polymerase chain reaction-confirmed SARS-CoV-2 infection in April–May 2022; 199,312/275,896 were included in analyses on the effect of a post-infection vaccine dose. Post-infection vaccination provided protection against reinfection (aHR 0.82; 95% confidence interval 0.79–0.85). For patients who had received one, two, or three vaccine doses before their first infection, hazard ratios for the post-infection vaccination effect were 0.84 (0.76–0.93), 0.87 (0.83–0.90), and 0.96 (0.74–1.23), respectively. Post-infection vaccination within 30 and 90 days before the second Omicron wave provided different degrees of protection (in aHR): 0.51 (0.44–0.58) and 0.67 (0.61–0.74), respectively. Moreover, for all vaccine types, but to different extents, a post-infection dose given to individuals who were fully vaccinated before first infection was protective.

    Conclusions:

    In previously vaccinated and infected individuals, an additional vaccine dose provided protection against Omicron variant reinfection. These observations will inform future policy decisions on COVID-19 vaccination in China and other countries.

    Funding:

    This study was funded the Key Discipline Program of Pudong New Area Health System (PWZxk2022-25), the Development and Application of Intelligent Epidemic Surveillance and AI Analysis System (21002411400), the Shanghai Public Health System Construction (GWVI-11.2-XD08), the Shanghai Health Commission Key Disciplines (GWVI-11.1-02), the Shanghai Health Commission Clinical Research Program (20214Y0020), the Shanghai Natural Science Foundation (22ZR1414600), and the Shanghai Young Health Talents Program (2022YQ076).

    1. Epidemiology and Global Health
    Marina Padilha, Victor Nahuel Keller ... Gilberto Kac
    Research Article Updated

    Background:

    The role of circulating metabolites on child development is understudied. We investigated associations between children’s serum metabolome and early childhood development (ECD).

    Methods:

    Untargeted metabolomics was performed on serum samples of 5004 children aged 6–59 months, a subset of participants from the Brazilian National Survey on Child Nutrition (ENANI-2019). ECD was assessed using the Survey of Well-being of Young Children’s milestones questionnaire. The graded response model was used to estimate developmental age. Developmental quotient (DQ) was calculated as the developmental age divided by chronological age. Partial least square regression selected metabolites with a variable importance projection ≥1. The interaction between significant metabolites and the child’s age was tested.

    Results:

    Twenty-eight top-ranked metabolites were included in linear regression models adjusted for the child’s nutritional status, diet quality, and infant age. Cresol sulfate (β=–0.07; adjusted-p <0.001), hippuric acid (β=–0.06; adjusted-p <0.001), phenylacetylglutamine (β=–0.06; adjusted-p <0.001), and trimethylamine-N-oxide (β=–0.05; adjusted-p=0.002) showed inverse associations with DQ. We observed opposite directions in the association of DQ for creatinine (for children aged –1 SD: β=–0.05; pP=0.01;+1 SD: β=0.05; p=0.02) and methylhistidine (–1 SD: β = - 0.04; p=0.04;+1 SD: β=0.04; p=0.03).

    Conclusions:

    Serum biomarkers, including dietary and microbial-derived metabolites involved in the gut-brain axis, may potentially be used to track children at risk for developmental delays.

    Funding:

    Supported by the Brazilian Ministry of Health and the Brazilian National Research Council.