Molecular epidemiology of Western equine encephalitis virus in Brazil, 2023-2024.

During the ongoing western equine encephalitis virus (WEEV) outbreak in South America, we described three fatal cases in horses from Rio Grande do Sul, Brazil. We sequenced WEEV strains and identified a novel lineage causing these cases. Continued surveillance and horse immunization are needed to mitigate the WEEV burden.

Evidence of Zika Virus Reinfection by Genome Diversity and Antibody Response Analysis, Brazil.

We generated 238 Zika virus (ZIKV) genomes from 135 persons in Brazil who had samples collected over 1 year to evaluate virus persistence. Phylogenetic inference clustered the genomes together with previously reported ZIKV strains from northern Brazil, showing that ZIKV has been remained relatively stable over time. Temporal phylogenetic analysis revealed limited within-host diversity among most ZIKV-persistent infected associated samples. However, we detected unusual virus temporal diversity from >5 persons, uncovering the existence of divergent genomes within the same patient. All those patients showed an increase in neutralizing antibody levels, followed by a decline at the convalescent phase of ZIKV infection. Of interest, in 3 of those patients, titers of neutralizing antibodies increased again after 6 months of ZIKV infection, concomitantly with real-time reverse transcription PCR re-positivity, supporting ZIKV reinfection events. Altogether, our findings provide evidence for the existence of ZIKV reinfection events.

Chikungunya virus infection in the southernmost state of Brazil was characterised by self-limited transmission (2017-2019) and a larger 2021 outbreak.

BACKGROUND: Chikungunya is a mosquito-borne virus that has been causing large outbreaks in the Americas since 2014. In Brazil, Asian-Caribbean (AC) and East-Central-South-African (ECSA) genotypes have been detected and lead to large outbreaks in several Brazilian states. In Rio Grande do Sul (RS), the southernmost state of Brazil, the first cases were reported in 2016. OBJECTIVES AND METHODS: We employed genome sequencing and epidemiological investigation to characterise the Chikungunya fever (CHIKF) burden in RS between 2017-2021. FINDINGS: We detected an increasing CHIKF burden linked to travel associated introductions and communitary transmission of distinct lineages of the ECSA genotype during this period. MAIN CONCLUSIONS: Until 2020, CHIKV introductions were most travel associated and transmission was limited. Then, in 2021, the largest outbreak occurred in the state associated with the introduction of a new ECSA lineage. CHIKV outbreaks are likely to occur in the near future due to abundant competent vectors and a susceptible population, exposing more than 11 million inhabitants to an increasing infection risk.

Reemergence of Dengue Virus Serotype 3, Brazil, 2023.

We characterized 3 autochthonous dengue virus serotype 3 cases and 1 imported case from 2 states in the North and South Regions of Brazil, 15 years after Brazil's last outbreak involving this serotype. We also identified a new Asian lineage recently introduced into the Americas, raising concerns about future outbreaks.

SARS-CoV-2 intra-host diversity, antibody response, and disease severity after reinfection by the variant of concern Gamma in Brazil.

The rapid spread of the SARS-CoV-2 Variant of Concern (VOC) Gamma in Amazonas during early 2021 fueled a second large COVID-19 epidemic wave and raised concern about the potential role of reinfections. Very few cases of reinfection associated with the VOC Gamma have been reported to date, and their potential impact on clinical, immunological, and virological parameters remains largely unexplored. Here we describe 25 cases of SARS-CoV-2 reinfection in Brazil. SARS-CoV-2 genomic analysis confirmed that individuals were primo-infected with distinct viral lineages between March and December 2020 (B.1.1, B.1.1.28, B.1.1.33, B.1.195, and P.2) and reinfected with the VOC Gamma between 3 to 12 months after primo-infection. We found a similar mean cycle threshold (Ct) value and limited intra-host viral diversity in both primo-infection and reinfection samples. Sera of 14 patients tested 10-75 days after reinfection displayed detectable neutralizing antibodies (NAb) titers against SARS-CoV-2 variants that circulated before (B.1.*), during (Gamma), and after (Delta and Omicron) the second epidemic wave in Brazil. All individuals had milder or no symptoms after reinfection, and none required hospitalization. These findings demonstrate that individuals reinfected with the VOC Gamma may display relatively high RNA viral loads at the upper respiratory tract after reinfection, thus contributing to onward viral transmissions. Despite this, our study points to a low overall risk of severe Gamma reinfections, supporting that the abrupt increase in hospital admissions and deaths observed in Amazonas and other Brazilian states during the Gamma wave was mostly driven by primary infections. Our findings also indicate that most individuals analyzed developed a high anti-SARS-CoV-2 NAb response after reinfection that may provide some protection against reinfection or disease by different SARS-CoV-2 variants.

Ancient origin of Jingchuvirales derived glycoproteins integrated in arthropod genomes.

Endogenous virus elements (EVEs) are viral-derived sequences integrated into their host genomes. EVEs of the Jingchuvirales order were detected in a wide range of insect genomes covering several distantly related families. Moreover, Jingchuvirales-derived glycoproteins were recently associated by our group with the origin of a putative new retrovirus based on a glycoprotein captured by a mosquito retrotransposon. But, except for mosquitoes, there is a lack of a more detailed understanding of the endogenization mechanism, timing, and frequency per Jingchuvirales viral lineages. Here we screened Jingchuvirales glycoprotein-derived EVEs (Jg-EVEs) in eukaryotic genomes. We found six distinct endogenization events of Jg-EVEs, that belong to two out of five known Jingchuvirales families (Chuviridae and Natareviridae). For seven arthropod families bearing Jg-EVEs there is no register of bona fide circulating chuvirus infection. Hence, our results show that Jingchuvirales viruses infected or still infect these host families. Although we found abundant evidence of LTR-Gypsy retrotransposons fragments associated with the glycoprotein in Hymenoptera and other insect orders, our results show that the widespread distribution of Jingchuvirales glycoproteins in extant Arhtropods is a result of multiple ancient endogenization events and that these virus fossils are being vertically inherited in Arthropods genomes for millions of years.

Comparative epidemic expansion of SARS-CoV-2 variants Delta and Omicron in the Brazilian State of Amazonas.

The SARS-CoV-2 variants of concern (VOCs) Delta and Omicron spread globally during mid and late 2021, respectively. In this study, we compare the dissemination dynamics of these VOCs in the Amazonas state, one of Brazil's most heavily affected regions. We sequenced the virus genome from 4128 patients collected in Amazonas between July 1st, 2021, and January 31st, 2022, and investigated the viral dynamics using a phylodynamic approach. The VOCs Delta and Omicron BA.1 displayed similar patterns of phylogeographic spread but different epidemic dynamics. The replacement of Gamma by Delta was gradual and occurred without an upsurge of COVID-19 cases, while the rise of Omicron BA.1 was extremely fast and fueled a sharp increase in cases. Thus, the dissemination dynamics and population-level impact of new SARS-CoV-2 variants introduced in the Amazonian population after mid-2021, a setting with high levels of acquired immunity, greatly vary according to their viral phenotype.

The Virome of Bats Inhabiting Brazilian Biomes: Knowledge Gaps and Biases towards Zoonotic Viruses.

Bats host a large variety of viruses, including some that may infect other vertebrates and humans. Research on bat-borne viruses attracted significant attention in recent years mainly due to epizootics caused by viruses having bats as hosts. The characterization of the viral communities of bats was then prioritized, but despite increasing efforts, there are large disparities in the geographical ranges covered and the methodologies employed around the world. As a result, large gaps remain in our current understanding of bat viromes and their role in disease emergence. This is particularly true for megadiverse regions in Latin America. This review aims to summarize the current understanding about bat viruses that inhabit Brazilian biomes, one of the most bat species-rich and diverse regions of the globe. Taking into account all known bat-associated viral families studied in Brazilian biomes, we found that almost half of all bat species (86/181 species) were not investigated for viruses at all. Moreover, only a small fraction of viral lineages or families have been studied more in depth, usually employing targeted methods with limited power to characterize a broad virus diversity. Additionally, these studies relied on limited spatiotemporal sampling and small sample sizes. Therefore, our current understanding of bat viral communities in the Brazilian biomes is limited and biased at different levels, limiting zoonotic risk assessments of bat-borne viruses. Considering these limitations, we propose strategies to bridge the existing gaps in the near future. IMPORTANCE Bat-borne viruses have attracted much attention due to zoonotic outbreaks with large consequences to humans. Because of that, virus characterization in bats has been prioritized in tropical regions of the globe. However, bat-virus research in Latin America and particularly in Brazil, which are among the most bat species-rich regions of the world, are highly biased toward zoonotic viruses and known bat reservoir species. These results have direct implication for virus studies in general but also for new zoonotic virus and spillover events characterization. The limited knowledge we currently have about the virome of Brazilian bats drastically limits any broad assessment of zoonotic viruses they carry and calls for coordinated and large-scale studies to fill this crucial knowledge gap.

Real-life evaluation of a rapid antigen test (DPP SARS-CoV-2 Antigen) for COVID-19 diagnosis of primary healthcare patients, in the context of the Omicron-dominant wave in Brazil.

OBJECTIVES: We aimed to investigate the real-life performance of the rapid antigen test in the context of a primary healthcare setting, including symptomatic and asymptomatic individuals that sought diagnosis during an Omicron infection wave. METHODS: We prospectively accessed the performance of the DPP SARS-CoV-2 Antigen test in the context of an Omicron-dominant real-life setting. We evaluated 347 unselected individuals (all-comers) from a public testing centre in Brazil, performing the rapid antigen test diagnosis at point-of-care with fresh samples. The combinatory result from two distinct real-time quantitative PCR (RT-qPCR) methods was employed as a reference and 13 samples with discordant PCR results were excluded. RESULTS: The assessment of the rapid test in 67 PCR-positive and 265 negative samples revealed an overall sensitivity of 80.5% (CI 95% = 69.1%-89.2%), specificity of 99.2% (CI 95% = 97.3%-99.1%) and positive/negative predictive values higher than 95%. However, we observed that the sensitivity was dependent on the viral load (sensitivity in Ct < 31 = 93.7%, CI = 82.8%-98.7%; Ct > 31 = 47.4%, CI = 24.4%-71.1%). The positive samples evaluated in the study were Omicron (BA.1/BA.1.1) by whole-genome sequencing (n = 40) and multiplex RT-qPCR (n = 17). CONCLUSIONS: Altogether, the data obtained from a real-life prospective cohort supports that the rapid antigen test sensitivity for Omicron remains high and underscores the reliability of the test for COVID-19 diagnosis in settings with high disease prevalence and limited PCR testing capability.

Chikungunya virus infection in the southernmost state of Brazil was characterised by self-limited transmission (2017-2019) and a larger 2021 outbreak

BACKGROUND Chikungunya is a mosquito-borne virus that has been causing large outbreaks in the Americas since 2014. In Brazil, Asian-Caribbean (AC) and East-Central-South-African (ECSA) genotypes have been detected and lead to large outbreaks in several Brazilian states. In Rio Grande do Sul (RS), the southernmost state of Brazil, the first cases were reported in 2016. OBJECTIVES AND METHODS We employed genome sequencing and epidemiological investigation to characterise the Chikungunya fever (CHIKF) burden in RS between 2017-2021. FINDINGS We detected an increasing CHIKF burden linked to travel associated introductions and communitary transmission of distinct lineages of the ECSA genotype during this period. MAIN CONCLUSIONS Until 2020, CHIKV introductions were most travel associated and transmission was limited. Then, in 2021, the largest outbreak occurred in the state associated with the introduction of a new ECSA lineage. CHIKV outbreaks are likely to occur in the near future due to abundant competent vectors and a susceptible population, exposing more than 11 million inhabitants to an increasing infection risk.

Multiple introductions and country-wide spread of DENV-2 genotype II (Cosmopolitan) in Brazil.

Dengue virus serotype 2, genotype Cosmopolitan (DENV-2-GII), is one of the most widespread DENV strains globally. In the USA, DENV-2 epidemics have been dominated by DENV-2 genotype Asian-American (DENV-2-GIII), and the first cases of DENV-2-GII were only described in 2019, in Peru, and in 2021 in Brazil. To gain new information about the circulation of DENV-2-GII in Brazil, we sequenced 237 DENV-2 confirmed cases sampled between March 2021 and March 2023 and revealed that DENV-2-GII is already present in all geographic regions of Brazil. The phylogeographic analysis inferred that DENV-2-GII was introduced at least four times in Brazil, between May 2020 and August 2022, generating multiple clades that spread throughout the country with different success. Despite multiple introductions of DENV-2-GII, analysis of the country-wide laboratory surveillance data showed that the Brazilian dengue epidemic in 2022 was dominated by DENV-1 in most states. We hypothesize that massive circulation of DENV-2-GIII in previous years in Brazil might have created a population immune barrier against symptomatic homotypic reinfections by DENV-2-GII, leading to sustained cryptic circulation in asymptomatic cases and localized outbreaks of this new genotype. In summary, our study stresses the importance of arboviral genomic surveillance to close monitoring and better understanding the potential impact of DENV-2-GII in the coming years.

Molecular characterization of a new SARS-CoV-2 recombinant cluster XAG identified in Brazil.

Recombination events have been described in the Coronaviridae family. Since the beginning of the SARS-CoV-2 pandemic, a variable degree of selection pressure has acted upon the virus, generating new strains with increased fitness in terms of viral transmission and antibody scape. Most of the SC2 variants of concern (VOC) detected so far carry a combination of key amino acid changes and indels. Recombination may also reshuffle existing genetic profiles of distinct strains, potentially giving origin to recombinant strains with altered phenotypes. However, co-infection and recombination events are challenging to detect and require in-depth curation of assembled genomes and sequencing reds. Here, we present the molecular characterization of a new SARS-CoV-2 recombinant between BA.1.1 and BA.2.23 Omicron lineages identified in Brazil. We characterized four mutations that had not been previously described in any of the recombinants already identified worldwide and described the likely breaking points. Moreover, through phylogenetic analysis, we showed that the newly named XAG lineage groups in a highly supported monophyletic clade confirmed its common evolutionary history from parental Omicron lineages and other recombinants already described. These observations were only possible thanks to the joint effort of bioinformatics tools auxiliary in genomic surveillance and the manual curation of experienced personnel, demonstrating the importance of genetic, and bioinformatic knowledge in genomics.

Possible Occupational Infection of Healthcare Workers with Monkeypox Virus, Brazil.

We evaluated epidemiologic and molecular characteristics of monkeypox virus (MPXV) infections sampled from 2 healthcare nurses. Five days after collecting samples from an infected patient, the nurses showed typical MPXV manifestations; quantitative PCR and whole-genome sequencing confirmed MPXV infection, most likely transmitted through contact with fomites.

Phylogenetic-based inference reveals distinct transmission dynamics of SARS-CoV-2 lineages Gamma and P.2 in Brazil.

The COVID-19 epidemic in Brazil experienced two major lineage replacements until mid-2021. The first was driven by lineage P.2, in late 2020, and the second by lineage Gamma, in early 2021. To understand how these SARS-CoV-2 lineages spread in Brazil, we analyzed 11,724 genomes collected throughout the country between September 2020 and April 2021. Our findings indicate that lineage P.2 probably emerged in July 2020 in the Rio de Janeiro state and Gamma in November 2020 in the Amazonas state. Both states were the main hubs of viral disseminations to other Brazilian locations. We estimate that Gamma was 1.56-3.06 times more transmissible than P.2 in Rio de Janeiro and that the median effective reproductive number (Re) of Gamma varied according to the geographic context (Re = 1.59-3.55). In summary, our findings support that lineage Gamma was more transmissible and spread faster than P.2 in Brazil.

Unusual SARS-CoV-2 intrahost diversity reveals lineage superinfection.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has infected almost 200 million people worldwide by July 2021 and the pandemic has been characterized by infection waves of viral lineages showing distinct fitness profiles. The simultaneous infection of a single individual by two distinct SARS-CoV-2 lineages may impact COVID-19 disease progression and provides a window of opportunity for viral recombination and the emergence of new lineages with differential phenotype. Several hundred SARS-CoV-2 lineages are currently well phylogenetically defined, but two main factors have precluded major coinfection/codetection and recombination analysis thus far: (i) the low diversity of SARS-CoV-2 lineages during the first year of the pandemic, which limited the identification of lineage defining mutations necessary to distinguish coinfecting/recombining viral lineages; and the (ii) limited availability of raw sequencing data where abundance and distribution of intrasample/intrahost variability can be accessed. Here, we assembled a large sequencing dataset from Brazilian samples covering a period of 18 May 2020 to 30 April 2021 and probed it for unexpected patterns of high intrasample/intrahost variability. This approach enabled us to detect nine cases of SARS-CoV-2 coinfection with well characterized lineage-defining mutations, representing 0.61 % of all samples investigated. In addition, we matched these SARS-CoV-2 coinfections with spatio-temporal epidemiological data confirming its plausibility with the cocirculating lineages at the timeframe investigated. Our data suggests that coinfection with distinct SARS-CoV-2 lineages is a rare phenomenon, although it is certainly a lower bound estimate considering the difficulty to detect coinfections with very similar SARS-CoV-2 lineages and the low number of samples sequenced from the total number of infections.

Spread of Gamma (P.1) Sub-Lineages Carrying Spike Mutations Close to the Furin Cleavage Site and Deletions in the N-Terminal Domain Drives Ongoing Transmission of SARS-CoV-2 in Amazonas, Brazil.

The Amazonas was one of the most heavily affected Brazilian states by the COVID-19 epidemic. Despite a large number of infected people, particularly during the second wave associated with the spread of the Variant of Concern (VOC) Gamma (lineage P.1), SARS-CoV-2 continues to circulate in the Amazonas. To understand how SARS-CoV-2 persisted in a human population with a high immunity barrier, we generated 1,188 SARS-CoV-2 whole-genome sequences from individuals diagnosed in the Amazonas state from 1st January to 6th July 2021, of which 38 were vaccine breakthrough infections. Our study reveals a sharp increase in the relative prevalence of Gamma plus (P.1+) variants, designated Pango Lineages P.1.3 to P.1.6, harboring two types of additional Spike changes: deletions in the N-terminal (NTD) domain (particularly Δ144 or Δ141-144) associated with resistance to anti-NTD neutralizing antibodies or mutations at the S1/S2 junction (N679K or P681H) that probably enhance the binding affinity to the furin cleavage site, as suggested by our molecular dynamics simulations. As lineages P.1.4 (S:N679K) and P.1.6 (S:P681H) expanded (Re > 1) from March to July 2021, the lineage P.1 declined (Re < 1) and the median Ct value of SARS-CoV-2 positive cases in Amazonas significantly decreases. Still, we did not find an increased incidence of P.1+ variants among breakthrough cases of fully vaccinated patients (71%) in comparison to unvaccinated individuals (93%). This evidence supports that the ongoing endemic transmission of SARS-CoV-2 in the Amazonas is driven by the spread of new local Gamma/P.1 sublineages that are more transmissible, although not more efficient to evade vaccine-elicited immunity than the parental VOC. Finally, as SARS-CoV-2 continues to spread in human populations with a declining density of susceptible hosts, the risk of selecting more infectious variants or antibody evasion mutations is expected to increase. IMPORTANCE The continuous evolution of SARS-CoV-2 is an expected phenomenon that will continue to happen due to the high number of cases worldwide. The present study analyzed how a Variant of Concern (VOC) could still circulate in a population hardly affected by two COVID-19 waves and with vaccination in progress. Our results showed that the answer behind that was a new generation of Gamma-like viruses, which emerged locally carrying mutations that made it more transmissible and more capable of spreading, partially evading prior immunity triggered by natural infections or vaccines. With thousands of new cases daily, the current pandemics scenario suggests that SARS-CoV-2 will continue to evolve and efforts to reduce the number of infected subjects, including global equitable access to COVID-19 vaccines, are mandatory. Thus, until the end of pandemics, the SARS-CoV-2 genomic surveillance will be an essential tool to better understand the drivers of the viral evolutionary process.

ViralFlow: A Versatile Automated Workflow for SARS-CoV-2 Genome Assembly, Lineage Assignment, Mutations and Intrahost Variant Detection.

The COVID-19 pandemic is driven by Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2) that emerged in 2019 and quickly spread worldwide. Genomic surveillance has become the gold standard methodology used to monitor and study this fast-spreading virus and its constantly emerging lineages. The current deluge of SARS-CoV-2 genomic data generated worldwide has put additional pressure on the urgent need for streamlined bioinformatics workflows. Here, we describe a workflow developed by our group to process and analyze large-scale SARS-CoV-2 Illumina amplicon sequencing data. This workflow automates all steps of SARS-CoV-2 reference-based genomic analysis: data processing, genome assembly, PANGO lineage assignment, mutation analysis and the screening of intrahost variants. The pipeline is capable of processing a batch of around 100 samples in less than half an hour on a personal laptop or in less than five minutes on a server with 50 threads. The workflow presented here is available through Docker or Singularity images, allowing for implementation on laptops for small-scale analyses or on high processing capacity servers or clusters. Moreover, the low requirements for memory and CPU cores and the standardized results provided by ViralFlow highlight it as a versatile tool for SARS-CoV-2 genomic analysis.

Mosquito long non-coding RNAs are enriched with Transposable Elements.

Long non-coding RNAs (lncRNAs) lack coding capacity and mounting evidence suggests that they have a regulatory role in diverse organisms. Most knowledge about lncRNAs comes from studies on vertebrates, including a structural association between lncRNAs and transposable elements (TEs). TE sequences are genomic parasites found in all branches of life and are particularly active and abundant in insect genomes. Here we investigate the contribution of TEs to lncRNA biogenesis in Aedes albopictus and Culex quinquefasciatus. We found that a large fraction of lncRNA loci co-occurs with TE loci in both species. Around 40% of A. albopictus and 52% of C. quinquefasciatus lncRNAs show some association with TEs. Most of the lncRNA/TE associations are represented by TE-derived sequences that are expressed as one or all exons of lncRNAs, including five lncRNAs that seem to influence immune-related genes involved in antiviral response. The contribution of TEs to lncRNAs also varies among the different types of TEs. The Gypsi superfamily is particularly enriched in lncRNAs sequences. In sum, this study demonstrates that transposable elements substantially contribute to lncRNAs biogenesis in A. albopictus and C. quinquefasciatus and may have an impact on regulatory modulation in these species.