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.

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.

Mosquito genomes are frequently invaded by transposable elements through horizontal transfer.

Transposable elements (TEs) are mobile genetic elements that parasitize basically all eukaryotic species genomes. Due to their complexity, an in-depth TE characterization is only available for a handful of model organisms. In the present study, we performed a de novo and homology-based characterization of TEs in the genomes of 24 mosquito species and investigated their mode of inheritance. More than 40% of the genome of Aedes aegypti, Aedes albopictus, and Culex quinquefasciatus is composed of TEs, while it varied substantially among Anopheles species (0.13%-19.55%). Class I TEs are the most abundant among mosquitoes and at least 24 TE superfamilies were found. Interestingly, TEs have been extensively exchanged by horizontal transfer (172 TE families of 16 different superfamilies) among mosquitoes in the last 30 million years. Horizontally transferred TEs represents around 7% of the genome in Aedes species and a small fraction in Anopheles genomes. Most of these horizontally transferred TEs are from the three ubiquitous LTR superfamilies: Gypsy, Bel-Pao and Copia. Searching more than 32,000 genomes, we also uncovered transfers between mosquitoes and two different Phyla-Cnidaria and Nematoda-and two subphyla-Chelicerata and Crustacea, identifying a vector, the worm Wuchereria bancrofti, that enabled the horizontal spread of a Tc1-mariner element among various Anopheles species. These data also allowed us to reconstruct the horizontal transfer network of this TE involving more than 40 species. In summary, our results suggest that TEs are frequently exchanged by horizontal transfers among mosquitoes, influencing mosquito's genome size and variability.

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.

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.

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.

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.

Mitogenome of Coprophanaeus ensifer and phylogenetic analysis of the Scarabaeidae family (Coleoptera).

Several studies about the phylogenetic relationships of the Scarabaeinae subfamily (Coleoptera: Scarabaeidae) have been performed, but some phylogenetic uncertainties persist including the relationship and monophyly of different tribes and some genera. The aim of this study was to characterize the mitogenome of Coprophanaeus ensifer in order to establish its position within the Scarabaeidae family and to contribute to the resolution of some phylogenetic uncertainties. The mitogenome was sequenced on the Illumina HiSeq 4000, assembled using the Mitobim software and annotated in MITOS WebServer. The phylogenetic trees were reconstructed by Bayesian inference. The C. ensifer mitogenome is a molecule of 14,964 bp that contains the number and organization of the genes similar to those of most Coleoptera species. Phylogenetic reconstruction suggests monophyly of the tribe Phanaeini and supports the hypothesis that Coprini is a sister group of Phanaeini. The results also revealed the position of the tribe Oniticellini which is grouped with Onthophagini and Onitini. The geographic distribution of these species that form the most ancestral clade suggests with Scarabaeinae originated in Africa.

Polymorphisms in GSTE2 is associated with temephos resistance in Aedes aegypti.

The glutathione S-transferases (GSTs) are enzymes involved in several distinct biological processes. In insects, the GSTs, especially delta and epsilon classes, play a key role in the metabolism of xenobiotics used to control insect populations. Here, we investigated its potential role in temephos resistance, examining the GSTE2 gene from susceptible (RecL) and resistant (RecR) strains of the mosquito Aedes aegypti, vector for several pathogenic arboviruses. Total GST enzymatic activity and the GSTE2 gene expression profile were evaluated, with the GSTE2 cDNA and genomic loci sequenced from both strains. Recombinant GSTE2 and mutants were produced in a heterologous expression system and assayed for enzyme kinetic parameters. These proteins also had their 3D structure predicted through molecular modeling. Our results showed that RecR has a profile of total GST enzymatic activity higher than RecL, with the expression of the GSTE2 gene in resistant larvae increasing six folds. Four exclusive RecR mutations were observed (L111S, I150V, E178A and A198E), which were absent in the laboratory susceptible strains. The enzymatic activity of the recombinant GSTE2 showed different kinetic parameters, with the GSTE2 RecR showing an enhanced ability to metabolize its substrate. The I150V mutation was shown to induce significant changes in catalytic parameters and a 3D modeling of GSTE2 mapped two of the RecR changes (L111S and I150V) near the enzyme's catalytic pocket, also implying an impact on its catalytic activity. Our results reinforce a potential role for GSTE2 in the metabolic resistance phenotype while contributing to the understanding of the molecular basis for the resistance mechanism.

VHICA, a New Method to Discriminate between Vertical and Horizontal Transposon Transfer: Application to the Mariner Family within Drosophila.

Transposable elements (TEs) are genomic repeated sequences that display complex evolutionary patterns. They are usually inherited vertically, but can occasionally be transmitted between sexually independent species, through so-called horizontal transposon transfers (HTTs). Recurrent HTTs are supposed to be essential in life cycle of TEs, which are otherwise destined for eventual decay. HTTs also impact the host genome evolution. However, the extent of HTTs in eukaryotes is largely unknown, due to the lack of efficient, statistically supported methods that can be applied to multiple species sequence data sets. Here, we developed a new automated method available as a R package "vhica" that discriminates whether a given TE family was vertically or horizontally transferred, and potentially infers donor and receptor species. The method is well suited for TE sequences extracted from complete genomes, and applicable to multiple TEs and species at the same time. We first validated our method using Drosophila TE families with well-known evolutionary histories, displaying both HTTs and vertical transmission. We then tested 26 different lineages of mariner elements recently characterized in 20 Drosophila genomes, and found HTTs in 24 of them. Furthermore, several independent HTT events could often be detected within the same mariner lineage. The VHICA (Vertical and Horizontal Inheritance Consistence Analysis) method thus appears as a valuable tool to analyze the evolutionary history of TEs across a large range of species.

Phylogenetic positioning of the Antarctic alga Prasiola crispa (Trebouxiophyceae) using organellar genomes and their structural analysis.

Antarctica is one of the most difficult habitats for sustaining life on earth; organisms that live there have developed different strategies for survival. Among these organisms is the green alga Prasiola crispa, belonging to the class Trebouxiophyceae. The literature on P. crispa taxonomy is scarce, and many gaps in the evolutionary relationship with its closest relatives remain. The goal of this study was to analyze the evolutionary relationships between P. crispa and other green algae using plastid and mitochondrial genomes. In addition, we analyzed the synteny conservation of these genomes of P. crispa with those of closely related species. Based on the plastid genome, P. crispa grouped with Prasiolopsis sp. SAG 84.81, another Trebouxiophyceaen species from the Prasiola clade. Based on the mitochondrial genome analysis, P. crispa grouped with other Trebouxiophyceaen species but had a basal position. The structure of the P. crispa chloroplast genome had low synteny with Prasiolopsis sp. SAG 84.81, despite some conserved gene blocks. The same was observed in the mitochondrial genome compared with Coccomyxa subellipsoidea C-169. We were able to establish the phylogenetic position of P. crispa with other species of Trebouxiophyceae using its genomes. In addition, we described the plasticity of these genomes using a structural analysis. The plastid and mitochondrial genomes of P. crispa will be useful for further genetic studies, phylogenetic analysis and resource protection of P. crispa as well as for further phylogenetic analysis of Trebouxiophyceaen green algae.

HTT-DB: horizontally transferred transposable elements database.

MOTIVATION: Horizontal transfer of transposable (HTT) elements among eukaryotes was discovered in the mid-1980s. As then, >300 new cases have been described. New findings about HTT are revealing the evolutionary impact of this phenomenon on host genomes. In order to provide an up to date, interactive and expandable database for such events, we developed the HTT-DB database. RESULTS: HTT-DB allows easy access to most of HTT cases reported along with rich information about each case. Moreover, it allows the user to generate tables and graphs based on searches using Transposable elements and/or host species classification and export them in several formats. AVAILABILITY AND IMPLEMENTATION: This database is freely available on the web at http://lpa.saogabriel.unipampa.edu.br:8080/httdatabase. HTT-DB was developed based on Java and MySQL with all major browsers supported. Tools and software packages used are free for personal or non-profit projects. CONTACT: bdotto82@gmail.com or gabriel.wallau@gmail.com.

An evaluation of the ecological relationship between Drosophila species and their parasitoid wasps as an opportunity for horizontal transposon transfer.

Evidences of horizontal transfer, the exchange of genetic material between reproductively isolated species, have accumulated over the last decades, including for multicellular eukaryotic organisms. However, the mechanisms and ecological relationships that promote such phenomenon is still poorly known. Host-parasite interaction is one type of relationship usually pointed in the literature that could potentially increase the probability of the horizontal transfer between species, because the species involved in such relationships are generally in close contact. Transposable elements, which are well-known genomic parasites, are DNA entities that tend to be involved in horizontal transfer due to their ability to mobilize between different genomic locations. Using Drosophila species and their parasitoid wasps as a host-parasite model, we evaluated the hypothesis that horizontal transposon transfers (HTTs) are more frequent in this set of species than in species that do not exhibit a close ecological and phylogenetic relationship. For this purpose, we sequenced two sets of species using a metagenomic and single-species genomic sampling approach through next-generation DNA sequencing. The first set was composed of five generalist Drosophila (D. maculifrons, D. bandeirantorum, D. polymorpha, D. mercatorum and D. willistoni) species and their associated parasitoid wasps, whereas the second set was composed of D. incompta, which is a flower specialist species, and its parasitoid wasp. We did not find strong evidence of HTT in the two sets of Drosophila and wasp parasites. However, at least five cases of HTT were observed between the generalist and specialist Drosophila species. Moreover, we detected an HT event involving a Wolbachia lineage between generalist and specialist species, indicating that these endosymbiotic bacteria could play a role as HTT vectors. In summary, our results do not support the hypothesis of prevalent HTT between species with a host-parasite relationship, at least for the studied wasp-Drosophila pairs. Moreover, it suggests that other mechanisms or parasites are involved in promoting HTT between Drosophila species as the Wolbachia endosymbiotic bacteria.

Genomic landscape and evolutionary dynamics of mariner transposable elements within the Drosophila genus.

BACKGROUND: The mariner family of transposable elements is one of the most widespread in the Metazoa. It is subdivided into several subfamilies that do not mirror the phylogeny of these species, suggesting an ancient diversification. Previous hybridization and PCR studies allowed a partial survey of mariner diversity in the Metazoa. In this work, we used a comparative genomics approach to access the genus-wide diversity and evolution of mariner transposable elements in twenty Drosophila sequenced genomes. RESULTS: We identified 36 different mariner lineages belonging to six distinct subfamilies, including a subfamily not described previously. Wide variation in lineage abundance and copy number were observed among species and among mariner lineages, suggesting continuous turn-over. Most mariner lineages are inactive and contain a high proportion of damaged copies. We showed that, in addition to substitutions that rapidly inactivate copies, internal deletion is a major mechanism contributing to element decay and the generation of non-autonomous sublineages. Hence, 23% of copies correspond to several Miniature Inverted-repeat Transposable Elements (MITE) sublineages, the first ever described in Drosophila for mariner. In the most successful MITEs, internal deletion is often associated with internal rearrangement, which sheds light on the process of MITE origin. The estimation of the transposition rates over time revealed that all lineages followed a similar progression consisting of a rapid amplification burst followed by a rapid decrease in transposition. We detected some instances of multiple or ongoing transposition bursts. Different amplification times were observed for mariner lineages shared by different species, a finding best explained by either horizontal transmission or a reactivation process. Different lineages within one species have also amplified at different times, corresponding to successive invasions. Finally, we detected a preference for insertion into short TA-rich regions, which appears to be specific to some subfamilies. CONCLUSIONS: This analysis is the first comprehensive survey of this family of transposable elements at a genus scale. It provides precise measures of the different evolutionary processes that were hypothesized previously for this family based on PCR data analysis. mariner lineages were observed at almost all "life cycle" stages: recent amplification, subsequent decay and potential (re)-invasion or invasion of genomes.

Highly pathogenic avian influenza H5N1 virus infections in pinnipeds and seabirds in Uruguay: Implications for bird-mammal transmission in South America.

The highly pathogenic avian influenza viruses of clade 2.3.4.4b have caused unprecedented deaths in South American wild birds, poultry, and marine mammals. In September 2023, pinnipeds and seabirds appeared dead on the Uruguayan Atlantic coast. Sixteen influenza virus strains were characterized by real-time reverse transcription PCR and genome sequencing in samples from sea lions (Otaria flavescens), fur seals (Arctocephalus australis), and terns (Sterna hirundinacea). Phylogenetic and ancestral reconstruction analysis showed that these strains have pinnipeds most likely as the ancestral host, representing a recent introduction of clade 2.3.4.4b in Uruguay. The Uruguayan and closely related strains from Peru (sea lions) and Chile (sea lions and a human case) carry mammalian adaptative residues 591K and 701N in the viral polymerase basic protein 2 (PB2). Our findings suggest that clade 2.3.4.4b strains in South America may have spread from mammals to mammals and seabirds, revealing a new transmission route.

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.

A new lineage nomenclature to aid genomic surveillance of dengue virus.

Dengue virus (DENV) is currently causing epidemics of unprecedented scope in endemic settings and expanding to new geographical areas. It is therefore critical to track this virus using genomic surveillance. However, the complex patterns of viral genomic diversity make it challenging to use the existing genotype classification system. Here we propose adding two sub-genotypic levels of virus classification, named major and minor lineages. These lineages have high thresholds for phylogenetic distance and clade size, rendering them stable between phylogenetic studies. We present an assignment tool to show that the proposed lineages are useful for regional, national and sub-national discussions of relevant DENV diversity. Moreover, the proposed lineages are robust to classification using partial genome sequences. We provide a standardized neutral descriptor of DENV diversity with which we can identify and track lineages of potential epidemiological and/or clinical importance. Information about our lineage system, including methods to assign lineages to sequence data and propose new lineages, can be found at: dengue-lineages.org.

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.

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.