Dominance of Three Sublineages of the SARS-CoV-2 Delta Variant in Mexico

Sequence signatures within the genome of SARS-CoV-2 can be used to predict host source

We conducted an in silico analysis to better understand the potential factors impacting host adaptation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in white-tailed deer, humans, and mink due to the strong evidence of sustained transmission within these hosts. Classification models trained on single nucleotide and amino acid differences between samples effectively identified white-tailed deer-, human-, and mink-derived SARS-CoV-2. For example, the balanced accuracy score of Extremely Randomized Trees classifiers was 0.984 ± 0.006. Eighty-eight commonly identified predictive mutations are found at sites under strong positive and negative selective pressure. A large fraction of sites under selection (86.9%) or identified by machine learning (87.1%) are found in genes other than the spike. Some locations encoded by these gene regions are predicted to be B- and T-cell epitopes or are implicated in modulating the immune response suggesting that host adaptation may involve the evasion of the host immune system, modulation of the class-I major-histocompatibility complex, and the diminished recognition of immune epitopes by CD8+ T cells. Our selection and machine learning analysis also identified that silent mutations, such as C7303T and C9430T, play an important role in discriminating deer-derived samples across multiple clades. Finally, our investigation into the origin of the B.1.641 lineage from white-tailed deer in Canada discovered an additional human sequence from Michigan related to the B.1.641 lineage sampled near the emergence of this lineage. These findings demonstrate that machine-learning approaches can be used in combination with evolutionary genomics to identify factors possibly involved in the cross-species transmission of viruses and the emergence of novel viral lineages.IMPORTANCESevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a highly transmissible virus capable of infecting and establishing itself in human and wildlife populations, such as white-tailed deer. This fact highlights the importance of developing novel ways to identify genetic factors that contribute to its spread and adaptation to new host species. This is especially important since these populations can serve as reservoirs that potentially facilitate the re-introduction of new variants into human populations. In this study, we apply machine learning and phylogenetic methods to uncover biomarkers of SARS-CoV-2 adaptation in mink and white-tailed deer. We find evidence demonstrating that both non-synonymous and silent mutations can be used to differentiate animal-derived sequences from human-derived ones and each other. This evidence also suggests that host adaptation involves the evasion of the immune system and the suppression of antigen presentation. Finally, the methods developed here are general and can be used to investigate host adaptation in viruses other than SARS-CoV-2.

Genomic epidemiology of SARS-CoV-2 δ sublineages of the second wave of 2021 in Antioquia, Colombia

Introduction. During the development of the SARS-CoV-2 pandemic in Antioquia, we experienced epidemiological peaks related to the α, ɣ, β, ƛ, and δ variants. δ had the highest incidence and prevalence. This lineage is of concern due to its clinical manifestations and epidemiological characteristics. A total of 253 δ sublineages have been reported in the PANGOLIN database. The sublineage identification through genomic analysis has made it possible to trace their evolution and propagation. Objective. To characterize the genetic diversity of the different SARS-CoV-2 δ sublineages in Antioquia and to describe its prevalence. Materials and methods. We collected sociodemographic information from 2,675 samples, and obtained 1,115 genomes from the GISAID database between July 12th, 2021, and January 18th, 2022. From the analyzed genomes, 515 were selected because of their high coverage values (>90%) to perform phylogenetic analysis and to infer allele frequencies of mutations of interest. Results. We characterized 24 sublineages. The most prevalent was AY.25. Mutations of interest as L452R, P681R, and P681H were identified in this sublineage, comprising a frequency close to 0.99. Conclusions. This study identified that the AY.25 sublineage has a transmission advantage compared to the other δ sublineages. This attribute may be related to the presence of the L452R and P681R mutations associated in other studies with higher evasion of the immune system and less efficacy of drugs against SARS-CoV-2.

Exploration of low-frequency allelic variants of SARS-CoV-2 genomes reveals coinfections in Mexico occurred during periods of VOCs turnover

A total of 14 973 alleles in 29 661 sequenced samples collected between March 2021 and January 2023 by the Mexican Consortium for Genomic Surveillance (CoViGen-Mex) and collaborators were used to construct a thorough map of mutations of the Mexican SARS-CoV-2 genomic landscape containing Intra-Patient Minor Allelic Variants (IPMAVs), which are low-frequency alleles not ordinarily present in a genomic consensus sequence. This additional information proved critical in identifying putative coinfecting variants included alongside the most common variants, B.1.1.222, B.1.1.519, and variants of concern (VOCs) Alpha, Gamma, Delta, and Omicron. A total of 379 coinfection events were recorded in the dataset (a rate of 1.28 %), resulting in the first such catalogue in Mexico. The most common putative coinfections occurred during the spread of Delta or after the introduction of Omicron BA.2 and its descendants. Coinfections occurred constantly during periods of variant turnover when more than one variant shared the same niche and high infection rate was observed, which was dependent on the local variants and time. Coinfections might occur at a higher frequency than customarily reported, but they are often ignored as only the consensus sequence is reported for lineage identification.

Immune response to SARS-CoV-2 variants after immunization with different vaccines in Mexico

There is limited information on the antibody responses against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in subjects from developing countries with populations having a high incidence of co-morbidities. Here, we analysed the immunogenicity of homologous schemes using the ChAdOx1-S, Sputnik V, or BNT162b2 vaccines and the effect of a booster dose with ChAdOx1-S in middle-aged adults who were seropositive or seronegative to the SARS-CoV-2 spike protein before vaccination. The study was conducted post-vaccination with a follow-up of 4 months for antibody titre using enzyme-linked immunosorbent assay (ELISA) and pseudovirus (PV) neutralization assays (PNAs). All three vaccines elicited a superior IgG anti-receptor-binding domain (RBD) and neutralization response against the Alpha and Delta variants when administered to individuals with a previous infection by SARS-CoV-2. The booster dose spiked the neutralization activity among individuals with and without a prior SARS-CoV-2 infection. The ChAdOx1-S vaccine induced weaker antibody responses in infection-naive subjects. A follow-up of 4 months post-vaccination showed a drop in antibody titre, with about 20% of the infection-naive and 100% of SARS-CoV-2 pre-exposed participants with detectable neutralization capacity against Alpha pseudovirus (Alpha-PV) and Delta PV (Delta-PV). Our observations support the use of different vaccines in a country with high seroprevalence at the vaccination time.

SARS-CoV-2 Omicron variants BA.4 and BA.5 dominated the fifth COVID-19 epidemiological wave in Mexico

In Mexico, the BA.4 and BA.5 Omicron variants dominated the fifth epidemic wave (summer 2022), superseding BA.2, which had circulated during the inter-wave period. The present study uses genome sequencing and statistical and phylogenetic analyses to examine these variants' abundance, distribution, and genetic diversity in Mexico from April to August 2022. Over 35 % of the sequenced genomes in this period corresponded to the BA.2 variant, 8 % to the BA.4 and 56 % to the BA.5 variant. Multiple subvariants were identified, but the most abundant, BA.2.9, BA.2.12.1, BA.5.1, BA.5.2, BA.5.2.1 and BA.4.1, circulated across the entire country, not forming geographical clusters. Contrastingly, other subvariants exhibited a geographically restricted distribution, most notably in the Southeast region, which showed a distinct subvariant dynamic. This study supports previous results showing that this region may be a significant entry point and contributed to introducing and evolving novel variants in Mexico. Furthermore, a differential distribution was observed for certain subvariants among specific States through time, which may have contributed to the overall increased diversity observed during this wave compared to the previous ones. This study highlights the importance of sustaining genomic surveillance to identify novel variants that may impact public health.

SARS-CoV-2 BW lineage, a fast-growing Omicron variant from southeast Mexico bearing relevant escape mutations

PURPOSE

The swift expansion of the BW.1 SARS-CoV-2 variant coincided with a rapid increase of COVID-19 cases occurring in Southeast Mexico in October, 2022, which marked the start of Mexico's sixth epidemiological wave. In Yucatan, up to 92% (58 of 73) of weekly sequenced genomes between epidemiological week 42 and 47 were identified as either BW.1 or its descendant, BW.1.1 in the region, during the last trimester of 2022. In the current study, a comprehensive genomic comparison was carried out to characterize the evolutionary history of the BW lineage, identifying its origins and its most important mutations.

METHODS

An alignment of all the genomes of the BW lineage and its parental BA.5.6.2 variant was carried out to identify their mutations. A phylogenetic and ancestral sequence reconstruction analysis with geographical inference, as well as a longitudinal analysis of point mutations, were performed to trace back their origin and contrast them with key RBD mutations in variant BQ.1, one of the fastest-growing lineages to date.

RESULTS

Our ancestral reconstruction analysis portrayed Mexico as the most probable origin of the BW.1 and BW.1.1 variants. Two synonymous substitutions, T7666C and C14599T, support their Mexican origin, whereas other two mutations are specific to BW.1: S:N460K and ORF1a:V627I. Two additional substitutions and a deletion are found in its descending subvariant, BW.1.1. Mutations found in the receptor binding domain, S:K444T, S:L452R, S:N460K, and S:F486V in BW.1 have been reported to be relevant for immune escape and are also key mutations in the BQ.1 lineage.

CONCLUSIONS

BW.1 appears to have arisen in the Yucatan Peninsula in Southeast Mexico sometime around July 2022 during the fifth COVID-19 wave. Its rapid growth may be in part explained by the relevant escape mutations also found in BQ.1.

Clinical outcomes and phylogenetic analysis in reflection with three predominant clades of SARS-CoV-2 variants

BACKGROUND

The pandemic of coronavirus disease 2019 (COVID-19) has a broad spectrum of clinical manifestations. The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) undergoes continuous evolution, resulting in the emergence of several variants. Each variant has a different severity and mortality rate.

MATERIALS AND METHODS

In this study, 1174 COVID-19 patients were studied for mortality and severity over three SARS-CoV-2 predominating variant periods in 2021 and 2022 in Sulaimani Province, Iraq. In each period, a representative, variant virus was subjected to phylogenetic and molecular and clinical analysis.

RESULTS

Phylogenetic analysis revealed three SARS-CoV-2 variants, belonging to: Delta B.1.617.2, Omicron BA.1.17.2, and Omicron BA.5.6. The Delta variants showed more severe symptoms and a lower PCR-Ct value than Omicron variants regardless of gender, and only 4.3% of the cases were asymptomatic. The mortality rate was lower with Omicron (.5% for BA.5.2 and 1.3% for BA.1.17.2) compared with Delta variants (2.5%). The higher mortality rate with Delta variants was in males (2.84%), while that with Omicron BA1.17.2 and BA.5.2 was in females, 1.05% and .0%, respectively. Age group (≥70) years had the highest mortality rate; however, it was (.0%) in the age group (30-49) years with Omicron variants, compared with (.96%) in Delta variants.

CONCLUSIONS

There has been a surge in COVID-19 infection in the city due to the predominant lineages of SARS-CoV-2, B.1.617, Omicron BA.1.17.2 and Omicron BA.5.6, respectively. A higher PCR-Ct value and severity of the Delta variant over Omicron BA.1.17.2 and/or BA.5.2 variants were significantly correlated with a higher death rate in the same order.

Two-year follow-up of the COVID-19 pandemic in Mexico

Background

After the initial outbreak in China (December 2019), the World Health Organization declared COVID-19 a pandemic on March 11th, 2020. This paper aims to describe the first 2 years of the pandemic in Mexico.

Design and methods

This is a population-based longitudinal study. We analyzed data from the national COVID-19 registry to describe the evolution of the pandemic in terms of the number of confirmed cases, hospitalizations, deaths and reported symptoms in relation to health policies and circulating variants. We also carried out logistic regression to investigate the major risk factors for disease severity.

Results

From March 2020 to March 2022, the coronavirus disease 2019 (COVID-19) pandemic in Mexico underwent four epidemic waves. Out of 5,702,143 confirmed cases, 680,063 were hospitalized (11.9%), and 324,436 (5.7%) died. Even if there was no difference in susceptibility by gender, males had a higher risk of death (CFP: 7.3 vs. 4.2%) and hospital admission risk (HP: 14.4 vs. 9.5%). Severity increased with age. With respect to younger ages (0-17 years), the 60+ years or older group reached adjusted odds ratios of 9.63 in the case of admission and 53.05 (95% CI: 27.94-118.62) in the case of death. The presence of any comorbidity more than doubled the odds ratio, with hypertension-diabetes as the riskiest combination. While the wave peaks increased over time, the odds ratios for developing severe disease (waves 2, 3, and 4 to wave 1) decreased to 0.15 (95% CI: 0.12-0.18) in the fourth wave.

Conclusion

The health policy promoted by the Mexican government decreased hospitalizations and deaths, particularly among older adults with the highest risk of admission and death. Comorbidities augment the risk of developing severe illness, which is shown to rise by double in the Mexican population, particularly for those reported with hypertension-diabetes. Factors such as the decrease in the severity of the SARS-CoV2 variants, changes in symptomatology, and advances in the management of patients, vaccination, and treatments influenced the decrease in mortality and hospitalizations.

Prevalence of symptoms, comorbidities, and reinfections in individuals infected with Wild-Type SARS-CoV-2, Delta, or Omicron variants: a comparative study in western Mexico

Introduction

The variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been classified into variants of interest (VOIs) or concern (VOCs) to prioritize global monitoring and research on variants with potential risks to public health. The SARS-CoV-2 high-rate mutation can directly impact the clinical disease progression, epidemiological behavior, immune evasion, vaccine efficacy, and transmission rates. Therefore, epidemiological surveillance is crucial for controlling the COVID-19 pandemic. In the present study, we aimed to describe the prevalence of wild-type (WT) SARS-CoV-2 and Delta and Omicron variants in Jalisco State, Mexico, from 2021 to 2022, and evaluate the possible association of these variants with clinical manifestations of COVID-19.

Methods

Four thousand and ninety-eight patients diagnosed with COVID-19 by real-time PCR (COVIFLU, Genes2Life, Mexico) from nasopharyngeal samples from January 2021 to January 2022 were included. Variant identification was performed by the RT-qPCR Master Mut Kit (Genes2Life, Mexico). A study population follow-up was performed to identify patients who had experienced reinfection after being vaccinated.

Results and Discussion

Samples were grouped into variants according to the identified mutations: 46.3% were Omicron, 27.9% were Delta, and 25.8% were WT. The proportions of dry cough, fatigue, headache, muscle pain, conjunctivitis, fast breathing, diarrhea, anosmia, and dysgeusia were significantly different among the abovementioned groups (p < 0.001). Anosmia and dysgeusia were mainly found in WT-infected patients, while rhinorrhea and sore throat were more prevalent in patients infected with the Omicron variant. For the reinfection follow-up, 836 patients answered, from which 85 cases of reinfection were identified (9.6%); Omicron was the VOC that caused all reported reinfection cases. In this study, we demonstrate that the Omicron variant caused the biggest outbreak in Jalisco during the pandemic from late December 2021 to mid-February 2022 but with a less severe form than the one demonstrated by Delta and WT. The co-analysis of mutations and clinical outcomes is a public health strategy with the potential to infer mutations or variants that could increase disease severity and even be an indicator of long-term sequelae of COVID-19.