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.

Metagenomic analysis reveals differences in the co-occurrence and abundance of viral species in SARS-CoV-2 patients with different severity of disease.

BACKGROUND: SARS-CoV-2 infections have a wide spectrum of clinical manifestations whose causes are not completely understood. Some human conditions predispose to severe outcome, like old age or the presence of comorbidities, but many other facets, including coinfections with other viruses, remain poorly characterized. METHODS: In this study, the eukaryotic fraction of the respiratory virome of 120 COVID-19 patients was characterized through whole metagenomic sequencing. RESULTS: Genetic material from respiratory viruses was detected in 25% of all samples, whereas human viruses other than SARS-CoV-2 were found in 80% of them. Samples from hospitalized and deceased patients presented a higher prevalence of different viruses when compared to ambulatory individuals. Small circular DNA viruses from the Anneloviridae (Torque teno midi virus 8, TTV-like mini virus 19 and 26) and Cycloviridae families (Human associated cyclovirus 10), Human betaherpesvirus 6, were found to be significantly more abundant in samples from deceased and hospitalized patients compared to samples from ambulatory individuals. Similarly, Rotavirus A, Measles morbillivirus and Alphapapilomavirus 10 were significantly more prevalent in deceased patients compared to hospitalized and ambulatory individuals. CONCLUSIONS: Results show the suitability of using metagenomics to characterize a broader peripheric virological landscape of the eukaryotic virome in SARS-CoV-2 infected patients with distinct disease outcomes. Identified prevalent viruses in hospitalized and deceased patients may prove important for the targeted exploration of coinfections that may impact prognosis.

Genomic Analysis of Early SARS-CoV-2 Variants Introduced in Mexico.

The coronavirus disease 2019 (COVID-19) pandemic has affected most countries in the world. Studying the evolution and transmission patterns in different countries is crucial to enabling implementation of effective strategies for disease control and prevention. In this work, we present the full genome sequence for 17 SARS-CoV-2 isolates corresponding to the earliest sampled cases in Mexico. Global and local phylogenomics, coupled with mutational analysis, consistently revealed that these viral sequences are distributed within 2 known lineages, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineage A/G, containing mostly sequences from North America, and lineage B/S, containing mainly sequences from Europe. Based on the exposure history of the cases and on the phylogenomic analysis, we characterized 14 independent introduction events. Additionally, three cases with no travel history were identified. We found evidence that two of these cases represented local transmission cases occurring in Mexico during mid-March 2020, denoting the earliest events described for the country. Within this local transmission cluster, we also identified an H49Y amino acid change in the Spike protein. This mutation represents a homoplasy occurring independently through time and space and may function as a molecular marker to follow any further spread of these viral variants throughout the country. Our results provide a general picture of the SARS-CoV-2 variants introduced at the beginning of the outbreak in Mexico, setting the foundation for future surveillance efforts.IMPORTANCE Understanding the introduction, spread, and establishment of SARS-CoV-2 within distinct human populations as well as the evolution of the pandemics is crucial to implement effective control strategies. In this work, we report that the initial virus strains introduced in Mexico came from Europe and the United States and that the virus was circulating locally in the country as early as mid-March. We also found evidence for early local transmission of strains with a H49Y mutation in the Spike protein, which could be further used as a molecular marker to follow viral spread within the country and the region.

Emergence and spread of the potential variant of interest (VOI) B.1.1.519 of SARS-CoV-2 predominantly present in Mexico.

SARS-CoV-2 variants emerged in late 2020, and at least three variants of concern (B.1.1.7, B.1.351, and P1) have been reported by WHO. These variants have several substitutions in the spike protein that affect receptor binding; they exhibit increased transmissibility and may be associated with reduced vaccine effectiveness. In the present work, we report the identification of a potential variant of interest, harboring the mutations T478K, P681H, and T732A in the spike protein, within the newly named lineage B.1.1.519, that rapidly outcompeted the preexisting variants in Mexico and has been the dominant virus in the country during the first trimester of 2021.

Operon-mapper: a web server for precise operon identification in bacterial and archaeal genomes.

Summary: Operon-mapper is a web server that accurately, easily and directly predicts the operons of any bacterial or archaeal genome sequence. The operon predictions are based on the intergenic distance of neighboring genes as well as the functional relationships of their protein-coding products. To this end, Operon-mapper finds all the ORFs within a given nucleotide sequence, along with their genomic coordinates, orthology groups and functional relationships. We believe that Operon-mapper, due to its accuracy, simplicity and speed, as well as the relevant information that it generates, will be a useful tool for annotating and characterizing genomic sequences. Availability and implementation: http://biocomputo.ibt.unam.mx/operon_mapper/.

A novel endogenous betaretrovirus in the common vampire bat (Desmodus rotundus) suggests multiple independent infection and cross-species transmission events.

The Desmodus rotundus endogenous betaretrovirus (DrERV) is fixed in the vampire bat D. rotundus population and in other phyllostomid bats but is not present in all species from this family. DrERV is not phylogenetically related to Old World bat betaretroviruses but to betaretroviruses from rodents and New World primates, suggesting recent cross-species transmission. A recent integration age estimation of the provirus in some taxa indicates that an exogenous counterpart might have been in recent circulation.

Rhinovirus is an important pathogen in upper and lower respiratory tract infections in Mexican children.

BACKGROUND: Most of the studies characterizing the incidence of rhinovirus (RV) have been carried out in hospitalized children and in developed countries. In those studies, RV-C has been associated with more severe respiratory tract infections than RV species A and B. In this study we determined the frequency and diversity of RV strains associated with upper and lower respiratory tract infections (URTI, LRTI) in Mexico, and describe the clinical characteristics of the illness associated with different RV species. METHODS: A prospective surveillance of 526 and 250 children with URTI and LRTI was carried out. Respiratory samples were analyzed by RT-PCR for viruses. The 5' untranslated region of the RV genome was amplified and sequenced. RESULTS: In the case of URTI, 17.5% were positive for RV, while this virus was found in 24.8% of LRTI. The RV species was determined in 73 children with URTI: 61.6% were RV-A, 37% RV-C and, 1.4% RV-B; and in 43 children with LRTI: 51.2% were RV-A, 41.8% RV-C, and 7% RV-B. No significant differences in clinical characteristics were found in patients with RV-A or RV-C infections. A high genetic diversity of RV strains was found in both URTI and LRTI. CONCLUSIONS: Both RV-A and RV-C species were frequently found in hospitalized as well as in outpatient children. This study underlines the high prevalence and genetic diversity of RV strains in Mexico and the potential severity of disease associated with RV-A and RV-C infections.

ProOpDB: Prokaryotic Operon DataBase.

The Prokaryotic Operon DataBase (ProOpDB, http://operons.ibt.unam.mx/OperonPredictor) constitutes one of the most precise and complete repositories of operon predictions now available. Using our novel and highly accurate operon identification algorithm, we have predicted the operon structures of more than 1200 prokaryotic genomes. ProOpDB offers diverse alternatives by which a set of operon predictions can be retrieved including: (i) organism name, (ii) metabolic pathways, as defined by the KEGG database, (iii) gene orthology, as defined by the COG database, (iv) conserved protein domains, as defined by the Pfam database, (v) reference gene and (vi) reference operon, among others. In order to limit the operon output to non-redundant organisms, ProOpDB offers an efficient method to select the most representative organisms based on a precompiled phylogenetic distances matrix. In addition, the ProOpDB operon predictions are used directly as the input data of our Gene Context Tool to visualize their genomic context and retrieve the sequence of their corresponding 5' regulatory regions, as well as the nucleotide or amino acid sequences of their genes.

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.

RegulonDB version 7.0: transcriptional regulation of Escherichia coli K-12 integrated within genetic sensory response units (Gensor Units).

RegulonDB (http://regulondb.ccg.unam.mx/) is the primary reference database of the best-known regulatory network of any free-living organism, that of Escherichia coli K-12. The major conceptual change since 3 years ago is an expanded biological context so that transcriptional regulation is now part of a unit that initiates with the signal and continues with the signal transduction to the core of regulation, modifying expression of the affected target genes responsible for the response. We call these genetic sensory response units, or Gensor Units. We have initiated their high-level curation, with graphic maps and superreactions with links to other databases. Additional connectivity uses expandable submaps. RegulonDB has summaries for every transcription factor (TF) and TF-binding sites with internal symmetry. Several DNA-binding motifs and their sizes have been redefined and relocated. In addition to data from the literature, we have incorporated our own information on transcription start sites (TSSs) and transcriptional units (TUs), obtained by using high-throughput whole-genome sequencing technologies. A new portable drawing tool for genomic features is also now available, as well as new ways to download the data, including web services, files for several relational database manager systems and text files including BioPAX format.

Evaluation of tools for taxonomic classification of viruses.

Viruses are the most abundant infectious agents on earth, and they infect living organisms such as bacteria, plants and animals, among others. They play an important role in the balance of different ecosystems by modulating microbial populations. In humans, they are responsible for some common diseases and may cause severe illnesses. Viral metagenomic studies have become essential and offer the possibility to understand and extend the knowledge of virus diversity and functionality. For these approaches, an essential step is the classification of viral sequences. In this work, 11 taxonomic classification tools were compared by analysing their performances, in terms of sensitivity and precision, to classify reads at the species and family levels using the same (viral and nonviral) datasets and evaluation metrics, as well as their processing times and memory requirements. The results showed that factors such as richness (numbers of viral species in samples), taxonomic level in the classification and read length influence tool performance. High values of viral richness in samples decreased the performances of most tools. Additionally, the classifications were better at higher taxonomic levels, such as families, compared to lower taxonomic levels, such as species, and were more evident in short reads. The results also indicated that BLAST and Kraken2 were the best tools for classifying all types of reads, while FastViromeExplorer and VirusFinder were only good when used for long reads and Centrifuge, DIAMOND, and One Codex when used for short reads. Regarding nonviral datasets (human and bacterial), all tools correctly classified them as nonviral.

The fecal and oropharyngeal eukaryotic viromes of healthy infants during the first year of life are personal.

Using a metagenomic sequencing approach, we described and compared the diversity and dynamics of the oropharyngeal and fecal eukaryotic virome of nine asymptomatic children in a semi-rural community setting located in the State of Morelos, Mexico. Ninety oropharyngeal swabs and 97 fecal samples were collected starting 2 weeks after birth and monthly thereafter until 12 months of age. In both niches, more than 95% of the total sequence reads were represented by viruses that replicate either in humans or in plants. Regarding human viruses, three families were most abundant and frequent in the oropharynx: Herpesviridae, Picornaviridae, and Reoviridae; in fecal samples, four virus families predominated: Caliciviridae, Picornaviridae, Reoviridae, and Anelloviridae. Both niches showed a high abundance of plant viruses of the family Virgaviridae. Differences in the frequency and abundance of sequence reads and diversity of virus species were observed in both niches and throughout the year of study, with some viruses already present in the first months of life. Our results suggest that the children's virome is dynamic and likely shaped by the environment, feeding, and age. Moreover, composition analysis suggests that the virome composition is mostly individual. Whether this constant exposition to different viruses has a long-term impact on children's health or development remains to be studied.

Longitudinal anti-SARS-CoV-2 antibody immune response in acute and convalescent patients.

Despite global efforts to assess the early response and persistence of SARS-CoV-2 antibodies in patients infected with or recovered from COVID-19, our understanding of the factors affecting its dynamics remains limited. This work aimed to evaluate the early and convalescent immunity of outpatients infected with SARS-CoV-2 and to determine the factors that affect the dynamics and persistence of the IgM and IgG antibody response. Seropositivity of volunteers from Mexico City and the State of Mexico, Mexico, was evaluated by ELISA using the recombinant receptor-binding domain (RBD) of the SARS-CoV-2 Spike protein for 90 days, at different time points (1, 15, 45, 60, and 90 days) after molecular diagnosis (RT-qPCR). Gender, age range, body mass index (BMI), comorbidities, and clinical spectrum of disease were analyzed to determine associations with the dynamics of anti-SARS-CoV-2 antibodies. On 90 days post-infection, individuals with moderate and asymptomatic disease presented the lowest levels of IgM, while for IgG, at the same time, the highest levels occurred with mild and moderate disease. The IgM and IgG levels were related to the clinical spectrum of disease, BMI, and the presence/absence of comorbidities through regression trees. The results suggest that the dynamics of anti-SARS-CoV-2 IgM and IgG antibodies in outpatients could be influenced by the clinical spectrum of the disease. In addition, the persistence of antibodies against SARS-CoV-2 could be related to the clinical spectrum of the disease, BMI, and the presence/absence of comorbidities.

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.

Omicron-BA.1 Dispersion Rates in Mexico Varied According to the Regional Epidemic Patterns and the Diversity of Local Delta Subvariants.

PURPOSE: The Omicron subvariant BA.1 of SARS-CoV-2 was first detected in November 2021 and quickly spread worldwide, displacing the Delta variant. In this work, a characterization of the spread of this variant in Mexico is presented. METHODS: The time to fixation of BA.1, the diversity of Delta sublineages, the population density, and the level of virus circulation during the inter-wave interval were determined to analyze differences in BA.1 spread. RESULTS: BA.1 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 Mexico. Unlike previous variants, BA.1 did not exhibit a geographically distinct circulation pattern. However, a regional difference in the speed of the replacement of the Delta variant was observed. CONCLUSIONS: 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, in addition to population density. Nonetheless, if there is a significant difference in the fitness of the variants, or if the time allowed for the competition is sufficiently long, it seems the fitter virus will eventually become dominant, as observed in the eventual dominance of the BA.1.x variant in Mexico.

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

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.

High accuracy operon prediction method based on STRING database scores.

We present a simple and highly accurate computational method for operon prediction, based on intergenic distances and functional relationships between the protein products of contiguous genes, as defined by STRING database (Jensen,L.J., Kuhn,M., Stark,M., Chaffron,S., Creevey,C., Muller,J., Doerks,T., Julien,P., Roth,A., Simonovic,M. et al. (2009) STRING 8-a global view on proteins and their functional interactions in 630 organisms. Nucleic Acids Res., 37, D412-D416). These two parameters were used to train a neural network on a subset of experimentally characterized Escherichia coli and Bacillus subtilis operons. Our predictive model was successfully tested on the set of experimentally defined operons in E. coli and B. subtilis, with accuracies of 94.6 and 93.3%, respectively. As far as we know, these are the highest accuracies ever obtained for predicting bacterial operons. Furthermore, in order to evaluate the predictable accuracy of our model when using an organism's data set for the training procedure, and a different organism's data set for testing, we repeated the E. coli operon prediction analysis using a neural network trained with B. subtilis data, and a B. subtilis analysis using a neural network trained with E. coli data. Even for these cases, the accuracies reached with our method were outstandingly high, 91.5 and 93%, respectively. These results show the potential use of our method for accurately predicting the operons of any other organism. Our operon predictions for fully-sequenced genomes are available at http://operons.ibt.unam.mx/OperonPredictor/.

High Prevalence and Diversity of Caliciviruses in a Community Setting Determined by a Metagenomic Approach.

We recently carried out a metagenomic study to determine the fecal virome of infants during their first year of life in a semirural community in Mexico. A total of 97 stool samples from nine children were collected starting 2 weeks after birth and monthly thereafter until 12 months of age. In this work, we describe the prevalence and incidence of caliciviruses in this birth cohort. We found that 54 (56%) and 24 (25%) of the samples were positive for norovirus and sapovirus sequence reads detected by next-generation sequencing, respectively. Potential infections were arbitrarily considered when at least 20% of the complete virus genome was determined. Considering only these samples, there were 3 cases per child/year for norovirus and 0.33 cases per child/year for sapovirus. All nine children had sequence reads related to norovirus in at least 2 and up to 10 samples, and 8 children excreted sapovirus sequence reads in 1 and up to 5 samples during the study. The virus in 35 samples could be genotyped. The results showed a high diversity of both norovirus (GI.3[P13], GI.5, GII.4, GII.4[P16], GII.7[P7], and GII.17[P17]) and sapovirus (GI.1, GI.7, and GII.4) in the community. Of interest, despite the frequent detection of caliciviruses in the stools, all children remained asymptomatic during the study. Our results clearly show that metagenomic studies in stools may reveal a detailed picture of the prevalence and diversity of gastrointestinal viruses in the human gut during the first year of life. IMPORTANCE Human caliciviruses are important etiological agents of acute gastroenteritis in children under 5 years of age. Several studies have characterized their association with childhood diarrhea and their presence in nondiarrheal stool samples. In this work, we used a next-generation sequencing approach to determine, in a longitudinal study, the fecal virome of infants during their first year of life. Using this method, we found that caliciviruses can be detected significantly more frequently than previously reported, providing a more detailed picture of the prevalence and genetic diversity of these viruses in the human gut during early life.

Genomic Characterization of SARS-CoV-2 Isolated from Patients with Distinct Disease Outcomes in Mexico.

The coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), has shown a wide spectrum of clinical manifestations ranging from asymptomatic infections to severe disease and death. Pre-existing medical conditions and age have been mainly linked to the development of severe disease; however, the potential association of viral genetic characteristics with different clinical conditions remains unclear. SARS-CoV-2 variants with increased transmissibility were detected early in the pandemics, and several variants with potential relevance for public health are currently circulating around the world. In this study, we characterized 57 complete SARS-CoV-2 genomes during the exponential growth phase of the early epidemiological curve in Mexico, in April 2020. Patients were categorized under distinct disease severity outcomes: mild disease or ambulatory care, severe disease or hospitalized, and deceased. To reduce bias related to risk factors, the patients were less than 60 years old and with no diagnosed comorbidities A trait-association phylogenomic approach was used to explore genotype-phenotype associations, represented by the co-occurrence of mutations, disease severity outcome categories, and clusters of Mexican sequences. Phylogenetic results revealed a higher genomic diversity compared to the initial viruses detected during the early stage of the local epidemic. We identified a total of 90 single nucleotide variants compared to the Wuhan-Hu-1 genome, including 54 nonsynonymous mutations. We did not find evidence for the co-occurrence of mutations associated with specific disease outcomes. Therefore, in the group of patients studied, disease severity was likely mainly driven by the host genetic background and other demographic factors. IMPORTANCE The genetic association of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) with different clinical conditions remains unclear and needs further investigation. In this study, we characterized 57 complete SARS-CoV-2 genomes from patients in Mexico with distinct disease severity outcomes: mild disease or ambulatory care, severe disease or hospitalized, and deceased. To reduce bias related to risk factors the patients were less than 60 years old and with no diagnosed comorbidities. We did not find evidence for the co-occurrence of mutations associated with specific disease outcomes. Therefore, in the group of patients studied, disease severity was likely mainly driven by the host genetic background and other demographic factors.

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.