New Insights into Avian Infectious Bronchitis Virus in Colombia from Whole-Genome Analysis.

Infectious Bronchitis (IB) is a respiratory disease caused by a highly variable Gammacoronavirus, which generates a negative impact on poultry health worldwide. GI-11 and GI-16 lineages have been identified in South America based on Infectious Bronchitis virus (IBV) partial S1 sequences. However, full genome sequence information is limited. In this study we report, for the first time, the whole-genome sequence of IBV from Colombia. Seven IBV isolates obtained during 2012 and 2013 from farms with respiratory disease compatible with IB were selected and the complete genome sequence was obtained by NGS. According to S1 sequence phylogenetic analysis, six isolates belong to lineage GI-1 and one to lineage GVI-1. When whole genome was analyzed, five isolates were related to the vaccine strain Ma5 2016 and two showed mosaic genomes. Results from complete S1 sequence analysis provides further support for the hypothesis that GVI-1, considered a geographically confined lineage in Asia, could have originated in Colombia. Complete genome information reported in this research allow a deeper understanding of the phylogenetic evolution of variants and the recombination events between strains that are circulating worldwide, contributing to the knowledge of coronavirus in Latin America and the world.

Characterization of strains of Anaplasma marginale from clinical cases in bovine using major surface protein 1a in Uruguay.

The major surface protein 1a (MSP1a) gene has been used to characterize Anaplasma marginale genetic diversity. This pathogen causes significant productivity and economic losses to the cattle industry. The objective of the present study was to report the first characterization of A. marginale genetic diversity in Uruguay based on MSP1a genotypes and their putative relationship with Rhipicephalus microplus. This cross-sectional study was conducted between 2016 and 2020. The study included whole blood samples from clinical cases of bovine anaplasmosis obtained from 30 outbreaks located in six Uruguay territorial departments. Diagnosis was performed using Giemsa-stained smears and confirmed by nested Polymerase Chance Reaction (nPCR) targeting the A. marginale major surface protein 5 gene. The genetic diversity of A. marginale strains was characterized by analyzing the microsatellite and tandem repeats of MSP1a. Based on the microsatellite structure, four genotypes were identified. Genotype E was the most prevalent. Analysis of MSP1a tandem repeats showed 28 different strains from the combination of 31 repeats, with τ-10-15 and α-ß-ß-ß-Γ being the most common. Repeats Γ, ß, α, and γ were associated with the absence of R. microplus with statistical significance (p < 0.05). Molecular observations showed that 46.7% of the strains identified in our samples lacked the ability to bind to tick cells; therefore, they were probably transmitted by other vectors. Strain genetic diversity provides valuable information for understanding the epidemiological behavior of A. marginale and could contribute to the development of effective vaccines for the control of this disease.

Draft genome assembly for the colombian freshwater bocachico fish, Prochilodus magdalenae.

We report the first draft genome assembly for Prochilodus magdalenae, the leading representative species of the Prochilodontidae family in Colombia. This 1.2-Gb assembly, with a GC content of 42.0% and a repetitive content of around 31.0%, is in the range of previously reported characid species genomes. Annotation identified 34,725 nuclear genes, and BUSCO completeness value was 94.9%. Gene ontology and primary metabolic pathway annotations indicate similar gene profiles for P. magdalenae and the closest species with annotated genomes: blind cave fish (Astyanax mexicanus) and red piranha (Pygocentrus nattereri). A comparative analysis showed similar genome traits to other characid species. The fully sequenced and annotated mitochondrial genome reproduces the taxonomic classification of P. magdalenae and confirms the low mitochondrial genetic divergence inside the Prochilodus genus. Phylogenomic analysis, using nuclear single-copy orthologous genes, also confirmed the evolutionary position of the species. This genome assembly provides a high-resolution genetic resource for sustainable P. magdalenae management in Colombia and, as the first genome assembly for the Prochilodontidae family, will contribute to fish genomics throughout South America.

Emergence and Spread of a B.1.1.28-Derived P.6 Lineage with Q675H and Q677H Spike Mutations in Uruguay.

Uruguay controlled the viral dissemination during the first nine months of the SARS-CoV-2 pandemic. Unfortunately, towards the end of 2020, the number of daily new cases exponentially increased. Herein, we analyzed the country-wide genetic diversity of SARS-CoV-2 between November 2020 and April 2021. We identified that the most prevalent viral variant during the first epidemic wave in Uruguay (December 2020-February 2021) was a B.1.1.28 sublineage carrying Spike mutations Q675H + Q677H, now designated as P.6, followed by lineages P.2 and P.7. P.6 probably arose around November 2020, in Montevideo, Uruguay's capital department, and rapidly spread to other departments, with evidence of further local transmission clusters; it also spread sporadically to the USA and Spain. The more efficient dissemination of lineage P.6 with respect to P.2 and P.7 and the presence of mutations (Q675H and Q677H) in the proximity of the key cleavage site at the S1/S2 boundary suggest that P.6 may be more transmissible than other lineages co-circulating in Uruguay. Although P.6 was replaced by the variant of concern (VOC) P.1 as the predominant lineage in Uruguay since April 2021, the monitoring of the concurrent emergence of Q675H + Q677H in VOCs should be of worldwide interest.

First evidence of enterovirus A71 and echovirus 30 in Uruguay and genetic relationship with strains circulating in the South American region.

Human enteroviruses (EVs) comprise more than 100 types of coxsackievirus, echovirus, poliovirus and numbered enteroviruses, which are mainly transmitted by the faecal-oral route leading to diverse diseases such as aseptic meningitis, encephalitis, and acute flaccid paralysis, among others. Since enteroviruses are excreted in faeces, wastewater-based epidemiology approaches are useful to describe EV diversity in a community. In Uruguay, knowledge about enteroviruses is extremely limited. This study assessed the diversity of enteroviruses through Illumina next-generation sequencing of VP1-amplicons obtained by RT-PCR directly applied to viral concentrates of 84 wastewater samples collected in Uruguay during 2011-2012 and 2017-2018. Fifty out of the 84 samples were positive for enteroviruses. There were detected 27 different types belonging to Enterovirus A species (CVA2-A6, A10, A16, EV-A71, A90), Enterovirus B species (CVA9, B1-B5, E1, E6, E11, E14, E21, E30) and Enterovirus C species (CVA1, A13, A19, A22, A24, EV-C99). Enterovirus A71 (EV-A71) and echovirus 30 (E30) strains were studied more in depth through phylogenetic analysis, together with some strains previously detected by us in Argentina. Results unveiled that EV-A71 sub-genogroup C2 circulates in both countries at least since 2011-2012, and that the C1-like emerging variant recently entered in Argentina. We also confirmed the circulation of echovirus 30 genotypes E and F in Argentina, and reported the detection of genotype E in Uruguay. To the best of our knowledge this is the first report of the EV-A71 C1-like emerging variant in South-America, and the first report of EV-A71 and E30 in Uruguay.

Recurrent Dissemination of SARS-CoV-2 Through the Uruguayan-Brazilian Border.

Uruguay is one of the few countries in the Americas that successfully contained the coronavirus disease 19 (COVID-19) epidemic during the first half of 2020. Nevertheless, the intensive human mobility across the dry border with Brazil is a major challenge for public health authorities. We aimed to investigate the origin of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains detected in Uruguayan localities bordering Brazil as well as to measure the viral flux across this ∼1,100 km uninterrupted dry frontier. Using complete SARS-CoV-2 genomes from the Uruguayan-Brazilian bordering region and phylogeographic analyses, we inferred the virus dissemination frequency between Brazil and Uruguay and characterized local outbreak dynamics during the first months (May-July) of the pandemic. Phylogenetic analyses revealed multiple introductions of SARS-CoV-2 Brazilian lineages B.1.1.28 and B.1.1.33 into Uruguayan localities at the bordering region. The most probable sources of viral strains introduced to Uruguay were the Southeast Brazilian region and the state of Rio Grande do Sul. Some of the viral strains introduced in Uruguayan border localities between early May and mid-July were able to locally spread and originated the first outbreaks detected outside the metropolitan region. The viral lineages responsible for Uruguayan urban outbreaks were defined by a set of between four and 11 mutations (synonymous and non-synonymous) with respect to the ancestral B.1.1.28 and B.1.1.33 viruses that arose in Brazil, supporting the notion of a rapid genetic differentiation between SARS-CoV-2 subpopulations spreading in South America. Although Uruguayan borders have remained essentially closed to non-Uruguayan citizens, the inevitable flow of people across the dry border with Brazil allowed the repeated entry of the virus into Uruguay and the subsequent emergence of local outbreaks in Uruguayan border localities. Implementation of coordinated bi-national surveillance systems is crucial to achieve an efficient control of the SARS-CoV-2 spread across this kind of highly permeable borderland regions around the world.

Environmental Surveillance through Next-Generation Sequencing to Unveil the Diversity of Human Enteroviruses beyond the Reported Clinical Cases.

The knowledge about circulation of Human Enteroviruses (EVs) obtained through medical diagnosis in Argentina is scarce. Wastewater samples monthly collected in Córdoba, Argentina during 2011-2012, and then in 2017-2018 were retrospectively studied to assess the diversity of EVs in the community. Partial VP1 gene was amplified by PCR from wastewater concentrates, and amplicons were subject of next-generation sequencing and genetic analyses. There were 41 EVs detected, from which ~50% had not been previously reported in Argentina. Most of the characterized EVs (60%) were detected at both sampling periods, with similar values of intratype nucleotide diversity. Exceptions were enterovirus A71, coxsackievirus B4, echovirus 14, and echovirus 30, which diversified in 2017-2018. There was a predominance of types from EV-C in 2017-2018, evidencing a common circulation of these types throughout the year in the community. Interestingly, high genetic similarity was evidenced among environmental strains of echovirus 30 circulating in 2011-2012 and co-temporal isolates obtained from patients suffering aseptic meningitis in different locations of Argentina. This study provides an updated insight about EVs circulating in an important region of South America, and suggests a valuable role of wastewater-based epidemiology in predicting outbreaks before the onset of cases in the community.

Tracking the onset date of the community spread of SARS-CoV-2 in western countries.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) rapidly spread around the world during 2020, but the precise time in which the virus began to spread locally is difficult to trace for most countries. Here, we estimate the probable onset date of the community spread of SARS-CoV-2 for heavily affected countries from Western Europe and the Americas on the basis of the cumulative number of deaths reported during the early stage of the epidemic. Our results support that SARS-CoV-2 probably started to spread locally in all western countries analysed between mid-January and mid-February 2020, thus long before community transmission was officially recognised and control measures were implemented.

Coxsackievirus A10 causing hand-foot-and-mouth disease in Uruguay.

We report for the first time in South America an HFMD case associated with Coxsackievirus A10. The viral strain belongs to a lineage involved in important European outbreaks and probably entered Uruguay after 2017 with a Greek origin. These findings call for strengthening the regional surveillance of HFMD.

Evolutionary Dynamics and Dissemination Pattern of the SARS-CoV-2 Lineage B.1.1.33 During the Early Pandemic Phase in Brazil.

A previous study demonstrates that most of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) Brazilian strains fell in three local clades that were introduced from Europe around late February 2020. Here we investigated in more detail the origin of the major and most widely disseminated SARS-CoV-2 Brazilian lineage B.1.1.33. We recovered 190 whole viral genomes collected from 13 Brazilian states from February 29 to April 31, 2020 and combined them with other B.1.1 genomes collected globally. Our genomic survey confirms that lineage B.1.1.33 is responsible for a variable fraction of the community viral transmissions in Brazilian states, ranging from 2% of all SARS-CoV-2 genomes from Pernambuco to 80% of those from Rio de Janeiro. We detected a moderate prevalence (5-18%) of lineage B.1.1.33 in some South American countries and a very low prevalence (<1%) in North America, Europe, and Oceania. Our study reveals that lineage B.1.1.33 evolved from an ancestral clade, here designated B.1.1.33-like, that carries one of the two B.1.1.33 synapomorphic mutations. The B.1.1.33-like lineage may have been introduced from Europe or arose in Brazil in early February 2020 and a few weeks later gave origin to the lineage B.1.1.33. These SARS-CoV-2 lineages probably circulated during February 2020 and reached all Brazilian regions and multiple countries around the world by mid-March, before the implementation of air travel restrictions in Brazil. Our phylodynamic analysis also indicates that public health interventions were partially effective to control the expansion of lineage B.1.1.33 in Rio de Janeiro because its median effective reproductive number (R e ) was drastically reduced by about 66% during March 2020, but failed to bring it to below one. Continuous genomic surveillance of lineage B.1.1.33 might provide valuable information about epidemic dynamics and the effectiveness of public health interventions in some Brazilian states.

Tracking the onset date of the community spread of SARS-CoV-2 in western countries

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) rapidly spread around the world during 2020, but the precise time in which the virus began to spread locally is difficult to trace for most countries. Here, we estimate the probable onset date of the community spread of SARS-CoV-2 for heavily affected countries from Western Europe and the Americas on the basis of the cumulative number of deaths reported during the early stage of the epidemic. Our results support that SARS-CoV-2 probably started to spread locally in all western countries analysed between mid-January and mid-February 2020, thus long before community transmission was officially recognised and control measures were implemented.

Description of a human Bocavirus recombinant strain in the Americas.

Human bocaviruses (HBoV) are mainly associated with respiratory and gastroenteric infections. These viruses belong to the family Parvoviridae, genus Bocaparvovirus and are classified in four subtypes (HBoV1-4). Recombination and point mutation have been described as basis of parvovirus evolution. In this study three viral sequences were obtained from positives HBoV sewage samples collected in two Uruguayan cities and were characterised by different methods as recombinant strains. This recombination event was localised in the 5' end of VP1 gene and the parental strains belonged to subtypes 3 and 4. These three Uruguayan strains are identical at the nucleotide sequences in the analysed genome region of the virus. As far as we known, this study represents the first detection of HBoV recombinants strains in the Americas.

Description of a human Bocavirus recombinant strain in the Americas

Human bocaviruses (HBoV) are mainly associated with respiratory and gastroenteric infections. These viruses belong to the family Parvoviridae, genus Bocaparvovirus and are classified in four subtypes (HBoV1-4). Recombination and point mutation have been described as basis of parvovirus evolution. In this study three viral sequences were obtained from positives HBoV sewage samples collected in two Uruguayan cities and were characterised by different methods as recombinant strains. This recombination event was localised in the 5' end of VP1 gene and the parental strains belonged to subtypes 3 and 4. These three Uruguayan strains are identical at the nucleotide sequences in the analysed genome region of the virus. As far as we known, this study represents the first detection of HBoV recombinants strains in the Americas.

The HIV-1 Subtype B Epidemic in French Guiana and Suriname Is Driven by Ongoing Transmissions of Pandemic and Non-pandemic Lineages.

The HIV-1 subtype B epidemic in French Guiana and Suriname is characterized by the co-circulation of the globally disseminated "BPANDEMIC" lineage and of non-pandemic subtype B lineages of Caribbean origin (BCAR). To reconstruct the spatiotemporal pattern of spread of those viral lineages circulating in these two countries, a total of 361 HIV-1 subtype B pol sequences recovered from treatment-naive adult patients from French Guiana and Suriname between 2006 and 2012 were combined with BPANDEMIC and BCAR reference sequences. Major Guianese/Surinamese BPANDEMIC and BCAR lineages were identified by Maximum Likelihood phylogenetic analysis and the spatiotemporal and demographic parameters estimated using a Bayesian coalescent-based method. We detected four BCAR and three BPANDEMIC transmission chains of large size that together comprise most pandemic and non-pandemic subtype B sequences from French Guiana (≥52%) and Suriname (≥70%) here analyzed. These major lineages were probably introduced into French Guiana and Suriname from the Caribbean (BCAR) and North/South America (BPANDEMIC) between the middle 1970s and the late 1980s and spread among populations from both countries with roughly comparable demographic growth rates. We detected a significant trend for higher viral loads and higher proportion of homosexual/bisexual men among subjects infected with BPANDEMIC relative to BCAR strains in French Guiana. These results show that the HIV subtype B epidemic in French Guiana and Suriname has been driven by multiple active BCAR and BPANDEMIC transmission chains that arose since the middle 1970s onward and operate in both countries simultaneously. Although no significant differences in the epidemic potential of major BCAR and BPANDEMIC lineages were observed, relevant associations between the infecting subtype B lineage and epidemiological and clinical characteristics were detected in French Guiana.

Inferring population dynamics of HIV-1 subtype C epidemics in Eastern Africa and Southern Brazil applying different Bayesian phylodynamics approaches.

The subtype C Eastern Africa clade (CEA), a particularly successful HIV-1 subtype C lineage, has seeded several sub-epidemics in Eastern African countries and Southern Brazil during the 1960s and 1970s. Here, we characterized the past population dynamics of the major CEA sub-epidemics in Eastern Africa and Brazil by using Bayesian phylodynamic approaches based on coalescent and birth-death models. All phylodynamic models support similar epidemic dynamics and exponential growth rates until roughly the mid-1980s for all the CEA sub-epidemics. Divergent growth patterns, however, were supported afterwards. The Bayesian skygrid coalescent model (BSKG) and the birth-death skyline model (BDSKY) supported longer exponential growth phases than the Bayesian skyline coalescent model (BSKL). The BDSKY model uncovers patterns of a recent decline for the CEA sub-epidemics in Burundi/Rwanda and Tanzania (Re < 1) and a recent growth for Southern Brazil (Re > 1); whereas coalescent models infer an epidemic stabilization. To the contrary, the BSKG model captured a decline of Ethiopian CEA sub-epidemic between the mid-1990s and mid-2000s that was not uncovered by the BDSKY model. These results underscore that the joint use of different phylodynamic approaches may yield complementary insights into the past HIV population dynamics.

Genomic and structural features of the yellow fever virus from the 2016-2017 Brazilian outbreak.

Southeastern Brazil has been suffering a rapid expansion of a severe sylvatic yellow fever virus (YFV) outbreak since late 2016, which has reached one of the most populated zones in Brazil and South America, heretofore a yellow fever-free zone for more than 70 years. In the current study, we describe the complete genome of 12 YFV samples from mosquitoes, humans and non-human primates from the Brazilian 2017 epidemic. All of the YFV sequences belong to the modern lineage (sub-lineage 1E) of South American genotype I, having been circulating for several months prior to the December 2016 detection. Our data confirm that viral strains associated with the most severe YF epidemic in South America in the last 70 years display unique amino acid substitutions that are mainly located in highly conserved positions in non-structural proteins. Our data also corroborate that YFV has spread southward into Rio de Janeiro state following two main sylvatic dispersion routes that converged at the border of the great metropolitan area comprising nearly 12 million unvaccinated inhabitants. Our original results can help public health authorities to guide the surveillance, prophylaxis and control measures required to face such a severe epidemiological problem. Finally, it will also inspire other workers to further investigate the epidemiological and biological significance of the amino acid polymorphisms detected in the Brazilian 2017 YFV strains.

Tracing the origin of the NS1 A188V substitution responsible for recent enhancement of Zika virus Asian genotype infectivity

A recent study showed that infectivity of Zika virus (ZIKV) Asian genotype was enhanced by an alanine-to-valine amino acid substitution at residue 188 of the NS1 protein, but the precise time and location of origin of this mutation were not formally estimated. Here, we applied a Bayesian coalescent-based framework to estimate the age and location of the ancestral viral strain carrying the A188V substitution. Our results support that the ancestral ZIKV strain carrying the A188V substitution arose in Southeastern Asia at the early 2000s and circulated in that region for some time (5-10 years) before being disseminated to Southern Pacific islands and the Americas.

Tracing the origin of the NS1 A188V substitution responsible for recent enhancement of Zika virus Asian genotype infectivity.

A recent study showed that infectivity of Zika virus (ZIKV) Asian genotype was enhanced by an alanine-to-valine amino acid substitution at residue 188 of the NS1 protein, but the precise time and location of origin of this mutation were not formally estimated. Here, we applied a Bayesian coalescent-based framework to estimate the age and location of the ancestral viral strain carrying the A188V substitution. Our results support that the ancestral ZIKV strain carrying the A188V substitution arose in Southeastern Asia at the early 2000s and circulated in that region for some time (5-10 years) before being disseminated to Southern Pacific islands and the Americas.

Phylodynamics of Yellow Fever Virus in the Americas: new insights into the origin of the 2017 Brazilian outbreak.

Yellow fever virus (YFV) strains circulating in the Americas belong to two distinct genotypes (I and II) that have diversified into several concurrent enzootic lineages. Since 1999, YFV genotype I has spread outside endemic regions and its recent (2017) reemergence in non-endemic Southeastern Brazilian states fuels one of the largest epizootic of jungle Yellow Fever registered in the country. To better understand this phenomenon, we reconstructed the phylodynamics of YFV American genotypes using sequences from nine countries sampled along 60 years, including strains from Brazilian 2017 outbreak. Our analyses reveals that YFV genotypes I and II follow roughly similar evolutionary and demographic dynamics until the early 1990s, when a dramatic change in the diversification process of the genotype I occurred associated with the emergence and dissemination of a new lineage (here called modern). Trinidad and Tobago was the most likely source of the YFV modern-lineage that spread to Brazil and Venezuela around the late 1980s, where it replaced all lineages previously circulating. The modern-lineage caused all major YFV outbreaks detected in non-endemic South American regions since 2000, including the 2017 Brazilian outbreak, and its dissemination was coupled to the accumulation of several amino acid substitutions particularly within non-structural viral proteins.

Detection and molecular characterization of emergent GII.P17/GII.17 Norovirus in Brazil, 2015.

A newly GII.17 Kawazaki_2014 variant strain was detected recently in Brazil. Phylogenetic analysis reveals at least four independent introduction events of this lineage into this country that took place throughout 2014, coinciding with FIFA World Cup in Brazil, 2014, and Hong Kong has been identified as the most likely source of introduction. This variant emerged in Asia causing outbreaks and replacing prevalent GII.4. Emergence of GII.P17/GII.17 variant emphasizes the need for active laboratory surveillance for NoV including molecular epidemiology and studies on virus evolution.