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 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.

BackgroundPandemic type A (H1N1) influenza arose in early 2009, probably in Mexico and the United States, and reappeared in North America in September for seven more months. An amino acid substitution in the hemagglutinin (HA), D222G, has been reported in a significant proportion of patients with a severe and fatal outcome. We studied the prevalence of HA222 substitutions in patients in Mexico during the second wave and its association with clinical outcome and pathogenicity in a mouse model.MethodsThe nucleotide sequences of hemagglutinin (HA) from viruses collected from 77 patients were determined including 50 severe and fatal cases and 27 ambulatory cases. Deep sequencing was done on 5 samples from severe or fatal cases in order to determine the quasispecies proportion. Weight loss and mortality due to infection with cultured influenza viruses were analyzed in a mouse model.ResultsViruses from 14 out of 50 hospitalized patients (28%) had a non aspartic acid residue at the HA 222 position (nD222), while all 27 ambulatory patients had D222 (p = 0.0014). G222 was detected as sole species or in coexistence with N222 and D222 in 12 patients with severe disease including 8 who died. N222 in coexistence with D222 was detected in 1 patient who died and co-occurrence of A222 and V222, together with D222, was detected in another patient who died. The patients with a nD222 residue had higher mortality (71.4%), compared to the group with only D222 (22.2%, p = 0.0008). Four of the 14 viruses from hospitalized patients were cultured and intranasally infected into mice. Two viruses with G222 were lethal while a third virus, with G222, caused only mild illness in mice similar to the fourth virus that contained D222.ConclusionsWe confirm the elevated incidence of HA222 (H1N1)pdm09 variants in severe disease and mortality. Both clinical and mouse infection data support the idea that nD222 mutations contribute to increased severity of disease but additional determinants in disease outcome may be present.

Brazil experienced a large dengue virus (DENV) epidemic in 2019, highlighting a continuous struggle with effective control and public health preparedness. Using Oxford Nanopore sequencing, we led field and classroom initiatives for the monitoring of DENV in Brazil, generating 227 novel genome sequences of DENV1-2 from 85 municipalities (2015–2019). This equated to an over 50% increase in the number of DENV genomes from Brazil available in public databases. Using both phylogenetic and epidemiological models we retrospectively reconstructed the recent transmission history of DENV1-2. Phylogenetic analysis revealed complex patterns of transmission, with both lineage co-circulation and replacement. We identified two lineages within the DENV2 BR-4 clade, for which we estimated the effective reproduction number and pattern of seasonality. Overall, the surveillance outputs and training initiative described here serve as a proof-of-concept for the utility of real-time portable sequencing for research and local capacity building in the genomic surveillance of emerging viruses.

To assess the possible circulation of Zika virus (ZIKV) prior to the first documented case in Mexico, we reanalyzed the stored samples from the states of Veracruz and Yucatán, which were originally collected to test for dengue (DENV) and chikungunya (CHIKV) but were negative for these viruses despite the symptomatology. The samples were originally collected between the 30 and 46 epidemiological weeks (EW) when the ZIKV was not yet declared as a Public Health Emergency of International Concern (PHEIC). From the total 4016 negative samples, a total of one hundred samples, 50 from Veracruz (CHIK -DENV -) and 50 from Yucatán (4 CHIK -DENV -and 46 CHIK -or DENV -), were tested for Zika virus by using RT-PCR. Results showed that in Veracruz and Yucatán, 20 % (10/50) and 70 % (35/50) were, respectively, ZIKV positive, indicating unequivocally the presence of ZIKV at least since July 2015. We also tested non-confirmed suspect measles cases from early 2015 for ZIKV by RT-PCR. Remarkably in 11 Mexican states, 86 % (18/21) were positive with the earlier symptoms onset as early as May 2015. Finally, RT-PCR analyses on RNA extracted from Aedes aegypti mosquitoes captured from January to March 2015 showed the presence of ZIKV, strongly suggesting that the vector was already carrying the virus at the start of 2015.

The mosquito-borne chikungunya virus, an alphavirus of the Togaviridae family, is responsible for acute polyarthralgia epidemics. Here, we report the complete genome sequences of two chikungunya virus strains, InDRE04 and InDRE51, identified in the Mexican states of Jalisco and Chiapas in 2014. Phylogenetic analysis showed that both strains belong to the Asian genotype.

BackgroundHuman enterovirus D68 (EV‐D68) recently caused an increase in mild‐to‐severe pediatric respiratory cases in North America and some European countries. Even though few of these children presented with acute paralytic disease, direct causal relationship cannot yet be assumed.ObjectivesThe purposes of this report were to describe the clinical findings of an outbreak of EV‐D68 infection in Mexico City and identify the genetic relationship with previously reported strains.Patients/MethodsBetween September and December 2014, 126 nasopharyngeal samples (NPS) of hospitalized children <15 years of age with ARI were tested for the presence of respiratory viruses using a multiplex RT‐qPCR and EV‐D68‐specific RT‐qPCR. Clinical, epidemiological, and demographic data were collected and associated with symptomatology and viral infections. Phylogenetic analyses were performed using VP1 region.ResultsEnterovirus/rhinovirus infection was detected in 40 patients (31·7%), of which 24 patients were EV‐D68‐positive. EV‐D68 infection prevailed over September and October 2014 and was associated with neutrophilia and lymphopenia, and patients were more likely to develop hypoxemia. Phylogenetic analyses showed that Mexican EV‐D68 belongs to the new B1 clade.ConclusionsThis is the first EV‐D68 outbreak described in Mexico and occurred few weeks after the United States reported similar infections. Although EV‐D68 belongs to new B1 clade, no neurological affection was observed.