During the first year of the SARS-CoV-2 pandemic in Mexico, more than two million people were infected. In this study, we analyzed full genome sequences from 27 February 2020 to 28 February 2021 to characterize the geographical and temporal distribution of SARS-CoV-2 lineages and identify the most common circulating lineages during this period. We defined six different geographical regions with particular dynamics of lineage circulation. The Northeast and Northwest regions were the ones that exhibited the highest lineage diversity, while the Central south and South/Southeast regions presented less diversity with predominance of a certain lineage. Additionally, by late February 2021, lineage B.1.1.519 represented more than 89% of all circulating lineages in the country.

The COVID-19 outbreak has caused over three million deaths worldwide. Understanding the pathology of the disease and the factors that drive severe and fatal clinical outcomes is of special relevance. Studying the role of the respiratory microbiota in COVID-19 is especially important as the respiratory microbiota is known to interact with the host immune system, contributing to clinical outcomes in chronic and acute respiratory diseases. Here, we characterized the microbiota in the respiratory tract of patients with mild, severe, or fatal COVID-19, and compared it to healthy controls and patients with non-COVID-19-pneumonia. We comparatively studied the microbial composition, diversity, and microbiota structure between the study groups and correlated the results with clinical data. We found differences in the microbial composition for COVID-19 patients, healthy controls, and non-COVID-19 pneumonia controls. In particular, we detected a high number of potentially opportunistic pathogens associated with severe and fatal levels of the disease. Also, we found higher levels of dysbiosis in the respiratory microbiota of patients with COVID-19 compared to the healthy controls. In addition, we detected differences in diversity structure between the microbiota of patients with mild, severe, and fatal COVID-19, as well as the presence of specific bacteria that correlated with clinical variables associated with increased risk of mortality. In summary, our results demonstrate that increased dysbiosis of the respiratory tract microbiota in patients with COVID-19 along with a continuous loss of microbial complexity structure found in mild to fatal COVID-19 cases may potentially alter clinical outcomes in patients. Taken together, our findings identify the respiratory microbiota as a factor potentially associated with the severity of COVID-19.

The COVID-19 pandemic has affected most countries in the world. Studying the evolution and transmission patterns in different countries is crucial to implement 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, the SARS-CoV-2 lineage A/G, containing mostly sequences from North America, and the lineage B/S containing mainly sequences from Europe. Based on the exposure history of the cases and on the phylogenomic analysis, we characterized fourteen independent introduction events. Additionally, three cases with no travel history were identified. We found evidence that two of these cases represent 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 the H49Y amino acid change in the Spike protein. This mutation is a homoplasy occurring independently through time and space, and may function as a molecular marker to follow on any further spread of these viral variants throughout the country. Our results depict the 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 is crucial to implement effective control strategies as well as the evolution of the pandemics. In this work, we describe that the initial virus strains introduced in Mexico came from Europe and the United States and the virus was circulating locally in the country as early as mid-March. We also found evidence for early local transmission of strains having the mutation H49Y in the Spike protein, that could be further used as a molecular marker to follow viral spread within the country and the region.

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.

The emergence and rapid increase of the B.1.1.7 (Alpha) lineage of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) throughout the world were due to its increased transmissibility. However, it did not displace cocirculating lineages in most of Mexico, particularly B.1.1.519, which dominated the country from February to May 2021.