Highly divergent and diverse viral community infecting sylvatic mosquitoes from Northeast Brazil.

Mosquitoes can transmit several pathogenic viruses to humans, but their natural viral community is also composed of a myriad of other viruses such as insect-specific viruses (ISVs) and those that infect symbiotic microorganisms. Besides a growing number of studies investigating the mosquito virome, the majority are focused on few urban species, and relatively little is known about the virome of sylvatic mosquitoes, particularly in high biodiverse biomes such as the Brazilian biomes. Here, we characterized the RNA virome of 10 sylvatic mosquito species from Atlantic forest remains at a sylvatic-urban interface in Northeast Brazil employing a metatranscriptomic approach. A total of 16 viral families were detected. The phylogenetic reconstructions of 14 viral families revealed that the majority of the sequences are putative ISVs. The phylogenetic positioning and, in most cases, the association with a high RNA-dependent RNA polymerase amino acid divergence from other known viruses suggests that the viruses characterized here represent at least 34 new viral species. Therefore, the sylvatic mosquito viral community is predominantly composed of highly divergent viruses highlighting the limited knowledge we still have about the natural virome of mosquitoes in general. Moreover, we found that none of the viruses recovered were shared between the species investigated, and only one showed high identity to a virus detected in a mosquito sampled in Peru, South America. These findings add further in-depth understanding about the interactions and coevolution between mosquitoes and viruses in natural environments. IMPORTANCE: Mosquitoes are medically important insects as they transmit pathogenic viruses to humans and animals during blood feeding. However, their natural microbiota is also composed of a diverse set of viruses that cause no harm to the insect and other hosts, such as insect-specific viruses. In this study, we characterized the RNA virome of sylvatic mosquitoes from Northeast Brazil using unbiased metatranscriptomic sequencing and in-depth bioinformatic approaches. Our analysis revealed that these mosquitoes species harbor a diverse set of highly divergent viruses, and the majority comprises new viral species. Our findings revealed many new virus lineages characterized for the first time broadening our understanding about the natural interaction between mosquitoes and viruses. Finally, it also provided several complete genomes that warrant further assessment for mosquito and vertebrate host pathogenicity and their potential interference with pathogenic arboviruses.

Molecular Epidemiology of Western Equine Encephalitis Virus, South America, 2023-2024.

Western equine encephalitis virus (WEEV) is a mosquitoborne virus that reemerged in December 2023 in Argentina and Uruguay, causing a major outbreak. We investigated the outbreak using epidemiologic, entomological, and genomic analyses, focusing on WEEV circulation near the Argentina‒Uruguay border in Rio Grande do Sul state, Brazil. During November 2023‒April 2024, the outbreak in Argentina and Uruguay resulted in 217 human cases, 12 of which were fatal, and 2,548 equine cases. We determined cases on the basis of laboratory and clinical epidemiologic criteria. We characterized 3 fatal equine cases caused by a novel WEEV lineage identified through a nearly complete coding sequence analysis, which we propose as lineage C. Our findings highlight the importance of continued surveillance and equine vaccination to control future WEEV outbreaks in South America.

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.

Infection of South American coatis (Nasua nasua) with highly pathogenic avian influenza H5N1 virus displaying mammalian adaptive mutations.

Deadly outbreaks among poultry, wild birds, and carnivorous mammals by the highly pathogenic H5N1 virus of the clade 2.3.4.4b have been reported in South America. The increasing virus incidence in various mammal species poses a severe zoonotic and pandemic threat. In Uruguay, the clade 2.3.4.4b viruses were first detected in February 2023, affecting wild birds and backyard poultry. Three months after the first reported case in Uruguay, the disease affected a population of 23 coatis (Nasua) in an ecological park. Most animals became infected, likely directly or indirectly from wild birds in the park, and experienced sudden death. Five animals from the colony survived, and four of them developed antibodies. The genomes of the H5N1 strains infecting coatis belonged to the B3.2 genotype of the clade 2.3.4.4b. Genomes from coatis were closely associated with those infecting backyard poultry, but transmission likely occurred through wild birds. Notable, two genomes have a 627K substitution in the RNA polymerase PB2 subunit, a hallmark amino acid linked to mammalian adaptation. Our findings support the ability of the avian influenza virus of the 2.3.4.4b clade to infect and transmit among terrestrial mammals with high pathogenicity and undergo rapid adaptive changes. It also highlights the coatis' ability to develop immunity and naturally clear the infection.

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.