Genomic epidemiology of circulating SARS-CoV-2 variants during first two years of the pandemic in Colombia

The emergence of highly transmissible SARS-CoV-2 variants has led to surges in cases and the need for global genomic surveillance. While some variants rapidly spread worldwide, other variants only persist nationally. There is a need for more fine-scale analysis to understand transmission dynamics at a country scale. For instance, the Mu variant of interest, also known as lineage B.1.621, was first detected in Colombia and was responsible for a large local wave but only a few sporadic cases elsewhere. To provide a better understanding of the epidemiology of SARS-Cov-2 variants in Colombia, we used 14,049 complete SARS-CoV-2 genomes from the 32 states of Colombia, and performed Bayesian phylodynamic analyses to estimate the time of variants introduction, their respective effective reproductive number, and effective population size, and the impact of disease control measures. We detected a total of 188 SARS-CoV-2 Pango lineages circulating in Colombia since the start of the pandemic. We showed that the effective reproduction number oscillated drastically throughout the first two years of the pandemic, with Mu showing the highest transmissibility (Re and growth rate estimation). Our results reinforce that genomic surveillance programs are essential for countries to make evidence-driven interventions towards the emergence and circulation of novel SARS-CoV-2 variants.

[1]  M. Suchard,et al.  Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil , 2021, Science.

[2]  L. Luo,et al.  Transmission and containment of the SARS-CoV-2 Delta variant of concern in Guangzhou, China: A population-based study , 2022, PLoS neglected tropical diseases.

[3]  A. Iafrate,et al.  Multiple SARS-CoV-2 variants escape neutralization by vaccine-induced humoral immunity , 2021, Cell.

[4]  Yuelong Shu,et al.  GISAID: Global initiative on sharing all influenza data – from vision to reality , 2017, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[5]  R. Bouckaert,et al.  Real-Time Genomics for Tracking Severe Acute Respiratory Syndrome Coronavirus 2 Border Incursions after Virus Elimination, New Zealand , 2021, Emerging infectious diseases.

[6]  An efficient coalescent epoch model for Bayesian phylogenetic inference , 2021, bioRxiv.

[7]  L. Poon,et al.  Genomic epidemiology of SARS-CoV-2 under an elimination strategy in Hong Kong , 2022, Nature communications.

[8]  R. Neher,et al.  Nextclade: clade assignment, mutation calling and quality control for viral genomes , 2021, J. Open Source Softw..

[9]  J. Lopez Bernal,et al.  Effectiveness of Covid-19 Vaccines against the B.1.617.2 (Delta) Variant. Reply. , 2021, The New England journal of medicine.

[10]  M. V. Van Kerkhove,et al.  Increased transmissibility and global spread of SARS-CoV-2 variants of concern as at June 2021 , 2021, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[11]  H. Kishino,et al.  Dating of the human-ape splitting by a molecular clock of mitochondrial DNA , 2005, Journal of Molecular Evolution.

[12]  D. A. Jackson,et al.  Evaluating the Effects of SARS-CoV-2 Spike Mutation D614G on Transmissibility and Pathogenicity , 2020, Cell.

[13]  T. Choudhri,et al.  Policy Interventions, Social Distancing, and SARS-CoV-2 Transmission in the United States: A Retrospective State-level Analysis , 2021, The American Journal of the Medical Sciences.

[14]  Md. Rabiul Islam,et al.  New Coronavirus Variants are Creating More Challenges to Global Healthcare System: A Brief Report on the Current Knowledge , 2022, Clinical pathology.

[15]  O. Pybus,et al.  Bayesian coalescent inference of past population dynamics from molecular sequences. , 2005, Molecular biology and evolution.

[16]  N. Ajami,et al.  Genomic epidemiology of SARS-CoV-2 in Colombia , 2020, medRxiv.

[17]  S. Bhatt,et al.  Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe , 2020, Nature.

[18]  Gerardo del Cerro Santamaría,et al.  Characterization of the emerging B.1.621 variant of interest of SARS-CoV-2 , 2021, Infection, Genetics and Evolution.

[19]  Edward C. Holmes,et al.  A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology , 2020, Nature Microbiology.

[20]  M. Kraemer,et al.  Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa , 2021, Nature.

[21]  Victor M Corman,et al.  Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR , 2020, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[22]  Christiane,et al.  World Medical Association Declaration of Helsinki: ethical principles for medical research involving human subjects. , 2004, Journal international de bioethique = International journal of bioethics.

[23]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[24]  S. Eckstein Ethical principles for medical research involving human subjects. , 2001, European journal of emergency medicine : official journal of the European Society for Emergency Medicine.

[25]  K. V. Parag,et al.  Establishment & lineage dynamics of the SARS-CoV-2 epidemic in the UK , 2020, medRxiv.

[26]  A. Paniz‐Mondolfi,et al.  Impact of SARS‐CoV‐2 Mu variant on vaccine effectiveness: A comparative genomics study at the peak of the third wave in Bogota, Colombia , 2022, Journal of medical virology.

[27]  Guangchuang Yu,et al.  Using ggtree to Visualize Data on Tree‐Like Structures , 2020, Current protocols in bioinformatics.

[28]  M. Ajelli,et al.  Modeling transmission of SARS-CoV-2 Omicron in China , 2022, Nature Medicine.

[29]  N. Alvis-Guzmán,et al.  Severe Acute Respiratory Syndrome Coronavirus 2 Seroprevalence Among Adults in a Tropical City of the Caribbean Area, Colombia: Are We Much Closer to Herd Immunity Than Developed Countries? , 2020, Open forum infectious diseases.

[30]  A. Ko,et al.  Insights into the limited global spread of the immune evasive SARS-CoV-2 variant Mu , 2022, medRxiv.

[31]  N. G. Davies,et al.  Increased mortality in community-tested cases of SARS-CoV-2 lineage B.1.1.7 , 2021, Nature.

[32]  Khalil Abudahab,et al.  Microreact: visualizing and sharing data for genomic epidemiology and phylogeography , 2016, Microbial genomics.

[33]  Elizabeth Weirich,et al.  Guidelines for safe work practices in human and animal medical diagnostic laboratories. Recommendations of a CDC-convened, Biosafety Blue Ribbon Panel. , 2012, MMWR supplements.

[34]  Korbinian Strimmer,et al.  APE: Analyses of Phylogenetics and Evolution in R language , 2004, Bioinform..

[35]  A. Kaneda,et al.  Neutralization of the SARS-CoV-2 Mu Variant by Convalescent and Vaccine Serum , 2021, The New England journal of medicine.

[36]  N. Linton,et al.  Serial interval of novel coronavirus (COVID-19) infections , 2020, International Journal of Infectious Diseases.

[37]  D. Paredes-Sabja,et al.  Evolution and Epidemic Spread of SARS-CoV-2 in Colombia: A Year into the Pandemic , 2021, Vaccines.

[38]  S. Bonhoeffer,et al.  Birth–death skyline plot reveals temporal changes of epidemic spread in HIV and hepatitis C virus (HCV) , 2012, Proceedings of the National Academy of Sciences.

[39]  Colin Simpson,et al.  Genomic epidemiology reveals transmission patterns and dynamics of SARS-CoV-2 in Aotearoa New Zealand , 2020, Nature Communications.

[40]  O. Pybus,et al.  Assignment of epidemiological lineages in an emerging pandemic using the pangolin tool , 2021, Virus evolution.

[41]  R. Myers,et al.  Effectiveness of Covid-19 Vaccines against the B.1.617.2 (Delta) Variant , 2021, The New England journal of medicine.

[42]  M. Mulligan,et al.  High-titer neutralization of Mu and C.1.2 SARS-CoV-2 variants by vaccine-elicited antibodies of previously infected individuals , 2021, Cell Reports.

[43]  I. Jaja,et al.  Social distancing: how religion, culture and burial ceremony undermine the effort to curb COVID-19 in South Africa , 2020, Emerging microbes & infections.

[44]  G. Wallau,et al.  Evolutionary Dynamics and Dissemination Pattern of the SARS-CoV-2 Lineage B.1.1.33 During the Early Pandemic Phase in Brazil , 2021, Frontiers in Microbiology.

[45]  Trevor Bedford,et al.  Nextstrain: real-time tracking of pathogen evolution , 2017, bioRxiv.

[46]  G. Wallau,et al.  COVID-19 in Amazonas, Brazil, was driven by the persistence of endemic lineages and P.1 emergence , 2021, Nature Medicine.

[47]  Colin Simpson,et al.  Phylodynamics reveals the role of human travel and contact tracing in controlling COVID-19 in four island nations , 2020, medRxiv.

[48]  Andrew Rambaut,et al.  Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen) , 2016, Virus evolution.

[49]  Y. Liu,et al.  The reproductive number of the Delta variant of SARS-CoV-2 is far higher compared to the ancestral SARS-CoV-2 virus , 2021, Journal of travel medicine.

[50]  I. Majsterek,et al.  The Emerging Concern and Interest SARS-CoV-2 Variants , 2021, Pathogens.

[51]  A. Londei,et al.  Ranking the effectiveness of worldwide COVID-19 government interventions , 2020, Nature Human Behaviour.

[52]  Kazutaka Katoh,et al.  Multiple alignment of DNA sequences with MAFFT. , 2009, Methods in molecular biology.

[53]  R. Bouckaert,et al.  bModelTest: Bayesian phylogenetic site model averaging and model comparison , 2015, BMC Evolutionary Biology.

[54]  M. Lipsitch,et al.  Population impact of SARS-CoV-2 variants with enhanced transmissibility and/or partial immune escape , 2021, Cell.

[55]  Ryan R. Wick,et al.  Performance of neural network basecalling tools for Oxford Nanopore sequencing , 2019, Genome Biology.

[56]  M. Roser,et al.  Coronavirus Pandemic (COVID-19) , 2020 .

[57]  S. Mattar,et al.  Effectiveness of the CoronaVac® vaccine in a region of the Colombian Amazon, was herd immunity achieved? , 2022, Tropical Diseases, Travel Medicine and Vaccines.

[58]  Z. Cucunubá,et al.  COVID-19 spread, detection, and dynamics in Bogota, Colombia , 2021, Nature Communications.

[59]  H. van Bakel,et al.  The arrival and spread of SARS‐CoV‐2 in Colombia , 2020, Journal of medical virology.

[60]  M. Pallen,et al.  SARS-CoV-2 Variants of Interest and Concern naming scheme conducive for global discourse , 2021, Nature Microbiology.

[61]  Alexei J. Drummond,et al.  Bayesian Phylogeography Finds Its Roots , 2009, PLoS Comput. Biol..

[62]  J. Quick nCoV-2019 sequencing protocol v3 (LoCost) , 2020 .