Endogenous Viral Elements Are Widespread in Arthropod Genomes and Commonly Give Rise to PIWI-Interacting RNAs

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.

piRNA associates with immune diseases

PIWI-interacting RNA (piRNA) is the most abundant small non-coding RNA in animal cells, typically 26-31 nucleotides in length and it binds with PIWI proteins, a subfamily of Argonaute proteins. Initially discovered in germ cells, piRNA is well known for its role in silencing transposons and maintaining genome integrity. However, piRNA is also present in somatic cells as well as in extracellular vesicles and exosomes. While piRNA has been extensively studied in various diseases, particular cancer, its function in immune diseases remains unclear. In this review, we summarize current research on piRNA in immune diseases. We first introduce the basic characteristics, biogenesis and functions of piRNA. Then, we review the association of piRNA with different types of immune diseases, including autoimmune diseases, immunodeficiency diseases, infectious diseases, and other immune-related diseases. piRNA is considered a promising biomarker for diseases, highlighting the need for further research into its potential mechanisms in disease pathogenesis.

Genetic diversity of Flaviviridae and Rhabdoviridae EVEs in Aedes aegypti and Aedes albopictus on Hainan Island and the Leizhou Peninsula, China

BACKGROUND

Hainan Island and the Leizhou Peninsula, the southernmost part of mainland China, are areas where Aedes aegypti and Ae. albopictus are sympatric and are also high-incidence areas of dengue outbreaks in China. Many studies have suggested that Aedes endogenous viral components (EVEs) are enriched in piRNA clusters which can silence incoming viral genomes. Investigation the EVEs present in the piRNA clusters associated with viral infection of Aedes mosquitoes in these regions may provide a theoretical basis for novel transmission-blocking vector control strategies.

METHODS

In this study, specific primers for endogenous Flaviviridae elements (EFVEs) and endogenous Rhabdoviridae elements (ERVEs) were used to detect the distribution of Zika virus infection associated EVEs in the genomes of individuals of the two Aedes mosquitoes. Genetic diversity of EVEs with a high detection rate was also analyzed.

RESULTS

The results showed that many EVEs associated with Zika virus infection were detected in both Aedes species, with the detection rates were 47.68% to 100% in Ae. aegypti and 36.15% to 92.31% in sympatric Ae. albopictus populations. EVEs detection rates in another 17 Ae. albopictus populations ranged from 29.39% to 89.85%. Genetic diversity analyses of the four EVEs (AaFlavi53, AaRha61, AaRha91 and AaRha100) of Ae. aegypti showed that each had high haplotype diversity and low nucleotide diversity. The number of haplotypes in AaFlavi53 was 8, with the dominant haplotype being Hap_1 and the other 7 haplotypes being further mutated from Hap_1 in a lineage direction. In contrast, the haplotype diversity of the other three ERVEs (AaRha61, AaRha91 and AaRha100) was more diverse and richer, with the haplotype numbers were 9, 15 and 19 respectively. In addition, these EVEs all showed inconsistent patterns of both population differentiation and dispersal compared to neutral evolutionary genes such as the Mitochondrial COI gene.

CONCLUSION

The EFVEs and ERVEs tested were present at high frequencies in the field Aedes mosquito populations. The haplotype diversity of the EFVE AaFlavi53 was relatively lower and the three ERVEs (AaRha61, AaRha91, AaRha100) were higher. None of the four EVEs could be indicative of the genetic diversity of the Ae. aegypti population. This study provided theoretical support for the use of EVEs to block arbovirus transmission, but further research is needed into the mechanisms by which these EVEs are antiviral to Aedes mosquitoes.

Metatranscriptomic investigation of single Ixodes pacificus ticks reveals diverse microbes, viruses, and novel mRNA-like endogenous viral elements

Ticks are increasingly important vectors of human and agricultural diseases. While many studies have focused on tick-borne bacteria, far less is known about tick-associated viruses and their roles in public health or tick physiology. To address this, we investigated patterns of bacterial and viral communities across two field populations of western black-legged ticks (Ixodes pacificus). Through metatranscriptomic analysis of 100 individual ticks, we quantified taxon prevalence, abundance, and co-occurrence with other members of the tick microbiome. In addition to commonly found tick-associated microbes, we assembled 11 novel RNA virus genomes from Rhabdoviridae, Chuviridae, Picornaviridae, Phenuiviridae, Reoviridae, Solemovidiae, Narnaviridae and two highly divergent RNA virus genomes lacking sequence similarity to any known viral families. We experimentally verified the presence of these in I. pacificus ticks across several life stages. We also unexpectedly identified numerous virus-like transcripts that are likely encoded by tick genomic DNA, and which are distinct from known endogenous viral element-mediated immunity pathways in invertebrates. Taken together, our work reveals that I. pacificus ticks carry a greater diversity of viruses than previously appreciated, in some cases resulting in evolutionarily acquired virus-like transcripts. Our findings highlight how pervasive and intimate tick-virus interactions are, with major implications for both the fundamental biology and vectorial capacity of I. pacificus ticks.

IMPORTANCE

Ticks are increasingly important vectors of disease, particularly in the United States where expanding tick ranges and intrusion into previously wild areas has resulted in increasing human exposure to ticks. Emerging human pathogens have been identified in ticks at an increasing rate, and yet little is known about the full community of microbes circulating in various tick species, a crucial first step to understanding how they interact with each and their tick host, as well as their ability to cause disease in humans. We investigated the bacterial and viral communities of the Western blacklegged tick in California and found 11 previously uncharacterized viruses circulating in this population.

Unveiling the hidden viromes across the animal tree of life: insights from a taxonomic classification pipeline applied to invertebrates of 31 metazoan phyla

Invertebrates constitute the majority of animal species on Earth, including most disease-causing agents or vectors, with more diverse viromes when compared to vertebrates. Recent advancements in high-throughput sequencing have significantly expanded our understanding of invertebrate viruses, yet this knowledge remains biased toward a few well-studied animal lineages. In this study, we analyze invertebrate DNA and RNA viromes for 31 phyla using 417 publicly available RNA-Seq data sets from diverse environments in the marine-terrestrial and marine-freshwater gradients. This study aims to (i) estimate virome compositions at the family level for the first time across the animal tree of life, including the first exploration of the virome in several phyla, (ii) quantify the diversity of invertebrate viromes and characterize the structure of invertebrate-virus infection networks, and (iii) investigate host phylum and habitat influence on virome differences. Results showed that a set of few viral families of eukaryotes, comprising Retroviridae, Flaviviridae, and several families of giant DNA viruses, were ubiquitous and highly abundant. Nevertheless, some differences emerged between phyla, revealing for instance a less diverse virome in Ctenophora compared to the other animal phyla. Compositional analysis of the viromes showed that the host phylum explained over five times more variance in composition than its habitat. Moreover, significant similarities were observed between the viromes of some phylogenetically related phyla, which could highlight the influence of co-evolution in shaping invertebrate viromes.IMPORTANCEThis study significantly enhances our understanding of the global animal virome by characterizing the viromes of previously unexamined invertebrate lineages from a large number of animal phyla. It showcases the great diversity of viromes within each phylum and investigates the role of habitat shaping animal viral communities. Furthermore, our research identifies dominant virus families in invertebrates and distinguishes phyla with analogous viromes. This study sets the road toward a deeper understanding of the virome across the animal tree of life.

Rhabdoviral Endogenous Sequences Identified in the Leishmaniasis Vector Lutzomyia longipalpis Are Widespread in Sandflies from South America

Sandflies are known vectors of leishmaniasis. In the Old World, sandflies are also vectors of viruses while little is known about the capacity of New World insects to transmit viruses to humans. Here, we relate the identification of RNA sequences with homology to rhabdovirus nucleocapsids (NcPs) genes, initially in the Lutzomyia longipalpis LL5 cell lineage, named NcP1.1 and NcP2. The Rhabdoviridae family never retrotranscribes its RNA genome to DNA. The sequences here described were identified in cDNA and DNA from LL-5 cells and in adult insects indicating that they are transcribed endogenous viral elements (EVEs). The presence of NcP1.1 and NcP2 in the L. longipalpis genome was confirmed in silico. In addition to showing the genomic location of NcP1.1 and NcP2, we identified another rhabdoviral insertion named NcP1.2. Analysis of small RNA molecules derived from these sequences showed that NcP1.1 and NcP1.2 present a profile consistent with elements targeted by primary piRNAs, while NcP2 was restricted to the degradation profile. The presence of NcP1.1 and NcP2 was investigated in sandfly populations from South America and the Old World. These EVEs are shared by different sandfly populations in South America while none of the Old World species studied presented the insertions.

A tale of caution: How endogenous viral elements affect virus discovery in transcriptomic data

Large-scale metagenomic and -transcriptomic studies have revolutionized our understanding of viral diversity and abundance. In contrast, endogenous viral elements (EVEs), remnants of viral sequences integrated into host genomes, have received limited attention in the context of virus discovery, especially in RNA-Seq data. EVEs resemble their original viruses, a challenge that makes distinguishing between active infections and integrated remnants difficult, affecting virus classification and biases downstream analyses. Here, we systematically assess the effects of EVEs on a prototypical virus discovery pipeline, evaluate their impact on data integrity and classification accuracy, and provide some recommendations for better practices. We examined EVEs and exogenous viral sequences linked to Orthomyxoviridae, a diverse family of negative-sense segmented RNA viruses, in 13 genomic and 538 transcriptomic datasets of Culicinae mosquitoes. Our analysis revealed a substantial number of viral sequences in transcriptomic datasets. However, a significant portion appeared not to be exogenous viruses but transcripts derived from EVEs. Distinguishing between transcribed EVEs and exogenous virus sequences was especially difficult in samples with low viral abundance. For example, three transcribed EVEs showed full-length segments, devoid of frameshift and nonsense mutations, exhibiting sufficient mean read depths that qualify them as exogenous virus hits. Mapping reads on a host genome containing EVEs before assembly somewhat alleviated the EVE burden, but it led to a drastic reduction of viral hits and reduced quality of assemblies, especially in regions of the viral genome relatively similar to EVEs. Our study highlights that our knowledge of the genetic diversity of viruses can be altered by the underestimated presence of EVEs in transcriptomic datasets, leading to false positives and altered or missing sequence information. Thus, recognizing and addressing the influence of EVEs in virus discovery pipelines will be key in enhancing our ability to capture the full spectrum of viral diversity.

Non-retroviral Endogenous Viral Elements in Tephritid Fruit Flies Reveal Former Viral Infections Not Related to Known Circulating Viruses

A wide variety of insect-specific non-retroviral RNA viruses specifically infect insects. During viral infection, fragments of viral sequences can integrate into the host genomes creating non-retroviral endogenous viral elements (nrEVEs). Although the exact function of nrEVEs is so far unknown, some studies suggest that nrEVEs may interfere with virus replication by producing PIWI-interacting RNAs (piRNAs) that recognize and degrade viral RNAs through sequence complementarity. In this article, we identified the nrEVEs repertoire of ten species within the dipteran family Tephritidae (true fruit flies), which are considered a major threat to agriculture worldwide. Our results suggest that each of these species contains nrEVEs, although in limited numbers, and that nrEVE integration may have occurred both before and after speciation. Furthermore, the majority of nrEVEs originated from viruses with negative single-stranded RNA genomes and represent structural viral functions. Notably, these nrEVEs exhibit low similarity to currently known circulating viruses. To explore the potential role of nrEVEs, we investigated their transcription pattern and the production of piRNAs in different tissues of Ceratitis capitata. We successfully identified piRNAs that are complementary to the sequence of one nrEVE in C. capitata, thereby highlighting a potential link between nrEVEs and the piRNA pathway. Overall, our results provide valuable insights into the comparative landscape of nrEVEs in true fruit flies, contributing to the understanding of the intimate relation between fruit flies and their past and present viral pathogens.

Co-option of a non-retroviral endogenous viral element in planthoppers

Non-retroviral endogenous viral elements (nrEVEs) are widely dispersed throughout the genomes of eukaryotes. Although nrEVEs are known to be involved in host antiviral immunity, it remains an open question whether they can be domesticated as functional proteins to serve cellular innovations in arthropods. In this study, we found that endogenous toti-like viral elements (ToEVEs) are ubiquitously integrated into the genomes of three planthopper species, with highly variable distributions and polymorphism levels in planthopper populations. Three ToEVEs display exon‒intron structures and active transcription, suggesting that they might have been domesticated by planthoppers. CRISPR/Cas9 experiments revealed that one ToEVE in Nilaparvata lugens, NlToEVE14, has been co-opted by its host and plays essential roles in planthopper development and fecundity. Large-scale analysis of ToEVEs in arthropod genomes indicated that the number of arthropod nrEVEs is currently underestimated and that they may contribute to the functional diversity of arthropod genes.

Untangling an insect's virome from its endogenous viral elements

BACKGROUND

Insects are an important reservoir of viral biodiversity, but the vast majority of viruses associated with insects have not been discovered. Recent studies have employed high-throughput RNA sequencing, which has led to rapid advances in our understanding of insect viral diversity. However, insect genomes frequently contain transcribed endogenous viral elements (EVEs) with significant homology to exogenous viruses, complicating the use of RNAseq for viral discovery.

METHODS

In this study, we used a multi-pronged sequencing approach to study the virome of an important agricultural pest and prolific vector of plant pathogens, the potato aphid Macrosiphum euphorbiae. We first used rRNA-depleted RNAseq to characterize the microbes found in individual insects. We then used PCR screening to measure the frequency of two heritable viruses in a local aphid population. Lastly, we generated a quality draft genome assembly for M. euphorbiae using Illumina-corrected Nanopore sequencing to identify transcriptionally active EVEs in the host genome.

RESULTS

We found reads from two insect-specific viruses (a Flavivirus and an Ambidensovirus) in our RNAseq data, as well as a parasitoid virus (Bracovirus), a plant pathogenic virus (Tombusvirus), and two phages (Acinetobacter and APSE). However, our genome assembly showed that part of the 'virome' of this insect can be attributed to EVEs in the host genome.

CONCLUSION

Our work shows that EVEs have led to the misidentification of aphid viruses from RNAseq data, and we argue that this is a widespread challenge for the study of viral diversity in insects.

Mosquito community composition shapes virus prevalence patterns along anthropogenic disturbance gradients

Previously unknown pathogens often emerge from primary ecosystems, but there is little knowledge on the mechanisms of emergence. Most studies analyzing the influence of land-use change on pathogen emergence focus on a single host-pathogen system and often observe contradictory effects. Here, we studied virus diversity and prevalence patterns in natural and disturbed ecosystems using a multi-host and multi-taxa approach. Mosquitoes sampled along a disturbance gradient in Côte d'Ivoire were tested by generic RT-PCR assays established for all major arbovirus and insect-specific virus taxa including novel viruses previously discovered in these samples based on cell culture isolates enabling an unbiased and comprehensive approach. The taxonomic composition of detected viruses was characterized and viral infection rates according to habitat and host were analyzed. We detected 331 viral sequences pertaining to 34 novel and 15 previously identified viruses of the families Flavi-, Rhabdo-, Reo-, Toga-, Mesoni- and Iflaviridae and the order Bunyavirales. Highest host and virus diversity was observed in pristine and intermediately disturbed habitats. The majority of the 49 viruses was detected with low prevalence. However, nine viruses were found frequently across different habitats of which five viruses increased in prevalence towards disturbed habitats, in congruence with the dilution effect hypothesis. These viruses were mainly associated with one specific mosquito species (Culex nebulosus), which increased in relative abundance from pristine (3%) to disturbed habitats (38%). Interestingly, the observed increased prevalence of these five viruses in disturbed habitats was not caused by higher host infection rates but by increased host abundance, an effect tentatively named abundance effect. Our data show that host species composition is critical for virus abundance. Environmental changes that lead to an uneven host community composition and to more individuals of a single species are a key driver of virus emergence.

Endoparasitoid lifestyle promotes endogenization and domestication of dsDNA viruses

The accidental endogenization of viral elements within eukaryotic genomes can occasionally provide significant evolutionary benefits, giving rise to their long-term retention, that is, to viral domestication. For instance, in some endoparasitoid wasps (whose immature stages develop inside their hosts), the membrane-fusion property of double-stranded DNA viruses have been repeatedly domesticated following ancestral endogenizations. The endogenized genes provide female wasps with a delivery tool to inject virulence factors that are essential to the developmental success of their offspring. Because all known cases of viral domestication involve endoparasitic wasps, we hypothesized that this lifestyle, relying on a close interaction between individuals, may have promoted the endogenization and domestication of viruses. By analyzing the composition of 124 Hymenoptera genomes, spread over the diversity of this clade and including free-living, ecto, and endoparasitoid species, we tested this hypothesis. Our analysis first revealed that double-stranded DNA viruses, in comparison with other viral genomic structures (ssDNA, dsRNA, ssRNA), are more often endogenized and domesticated (that is, retained by selection) than expected from their estimated abundance in insect viral communities. Second, our analysis indicates that the rate at which dsDNA viruses are endogenized is higher in endoparasitoids than in ectoparasitoids or free-living hymenopterans, which also translates into more frequent events of domestication. Hence, these results are consistent with the hypothesis that the endoparasitoid lifestyle has facilitated the endogenization of dsDNA viruses, in turn, increasing the opportunities of domestications that now play a central role in the biology of many endoparasitoid lineages.

Endogenous viral elements reveal associations between a non-retroviral RNA virus and symbiotic dinoflagellate genomes

Endogenous viral elements (EVEs) offer insight into the evolutionary histories and hosts of contemporary viruses. This study leveraged DNA metagenomics and genomics to detect and infer the host of a non-retroviral dinoflagellate-infecting +ssRNA virus (dinoRNAV) common in coral reefs. As part of the Tara Pacific Expedition, this study surveyed 269 newly sequenced cnidarians and their resident symbiotic dinoflagellates (Symbiodiniaceae), associated metabarcodes, and publicly available metagenomes, revealing 178 dinoRNAV EVEs, predominantly among hydrocoral-dinoflagellate metagenomes. Putative associations between Symbiodiniaceae and dinoRNAV EVEs were corroborated by the characterization of dinoRNAV-like sequences in 17 of 18 scaffold-scale and one chromosome-scale dinoflagellate genome assembly, flanked by characteristically cellular sequences and in proximity to retroelements, suggesting potential mechanisms of integration. EVEs were not detected in dinoflagellate-free (aposymbiotic) cnidarian genome assemblies, including stony corals, hydrocorals, jellyfish, or seawater. The pervasive nature of dinoRNAV EVEs within dinoflagellate genomes (especially Symbiodinium), as well as their inconsistent within-genome distribution and fragmented nature, suggest ancestral or recurrent integration of this virus with variable conservation. Broadly, these findings illustrate how +ssRNA viruses may obscure their genomes as members of nested symbioses, with implications for host evolution, exaptation, and immunity in the context of reef health and disease.

Divergent RNA viruses infecting sea lice, major ectoparasites of fish

Sea lice, the major ectoparasites of fish, have significant economic impacts on wild and farmed finfish, and have been implicated in the decline of wild salmon populations. As blood-feeding arthropods, sea lice may also be reservoirs for viruses infecting fish. However, except for two groups of negative-strand RNA viruses within the order Mononegavirales, nothing is known about viruses of sea lice. Here, we used transcriptomic data from three key species of sea lice (Lepeophtheirus salmonis, Caligus clemensi, and Caligus rogercresseyi) to identify 32 previously unknown RNA viruses. The viruses encompassed all the existing phyla of RNA viruses, with many placed in deeply branching lineages that likely represent new families and genera. Importantly, the presence of canonical virus-derived small interfering RNAs (viRNAs) indicates that most of these viruses infect sea lice, even though in some cases their closest classified relatives are only known to infect plants or fungi. We also identified both viRNAs and PIWI-interacting RNAs (piRNAs) from sequences of a bunya-like and two qin-like viruses in C. rogercresseyi. Our analyses showed that most of the viruses found in C. rogercresseyi occurred in multiple life stages, spanning from planktonic to parasitic stages. Phylogenetic analysis revealed that many of the viruses infecting sea lice were closely related to those that infect a wide array of eukaryotes with which arthropods associate, including fungi and parasitic tapeworms, implying that over evolutionary time there has been cross-phylum and cross-kingdom switching of viruses between arthropods and other eukaryotes. Overall, this study greatly expands our view of virus diversity in crustaceans, identifies viruses that infect and replicate in sea lice, and provides evidence that over evolutionary time, viruses have switched between arthropods and eukaryotic hosts in other phyla and kingdoms.

Chromatin profiling identifies transcriptional readthrough as a conserved mechanism for piRNA biogenesis in mosquitoes

The piRNA pathway in mosquitoes differs substantially from other model organisms, with an expanded PIWI gene family and functions in antiviral defense. Here, we define core piRNA clusters as genomic loci that show ubiquitous piRNA expression in both somatic and germline tissues. These core piRNA clusters are enriched for non-retroviral endogenous viral elements (nrEVEs) in antisense orientation and depend on key biogenesis factors, Veneno, Tejas, Yb, and Shutdown. Combined transcriptome and chromatin state analyses identify transcriptional readthrough as a conserved mechanism for cluster-derived piRNA biogenesis in the vector mosquitoes Aedes aegypti, Aedes albopictus, Culex quinquefasciatus, and Anopheles gambiae. Comparative analyses between the two Aedes species suggest that piRNA clusters function as traps for nrEVEs, allowing adaptation to environmental challenges such as virus infection. Our systematic transcriptome and chromatin state analyses lay the foundation for studies of gene regulation, genome evolution, and piRNA function in these important vector species.

What Are the Functional Roles of Piwi Proteins and piRNAs in Insects?

Research on Piwi proteins and piRNAs in insects has focused on three experimental models: oogenesis and spermatogenesis in Drosophila melanogaster, the antiviral response in Aedes mosquitoes and the molecular analysis of primary and secondary piRNA biogenesis in Bombyx mori-derived BmN4 cells. Significant unique and complementary information has been acquired and has led to a greater appreciation of the complexity of piRNA biogenesis and Piwi protein function. Studies performed in other insect species are emerging and promise to add to the current state of the art on the roles of piRNAs and Piwi proteins. Although the primary role of the piRNA pathway is genome defense against transposons, particularly in the germline, recent findings also indicate an expansion of its functions. In this review, an extensive overview is presented of the knowledge of the piRNA pathway that so far has accumulated in insects. Following a presentation of the three major models, data from other insects were also discussed. Finally, the mechanisms for the expansion of the function of the piRNA pathway from transposon control to gene regulation were considered.

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.

Integrated Jingmenvirus Polymerase Gene in Ixodes ricinus Genome

Members of the jingmenviruses group have been found in arthropods and mammals on all continents except Australia and Antarctica. Two viruses of this group were isolated from patients with fever after a tick bite. Using a nested RT-PCR assay targeting a jingmenvirus polymerase gene fragment, we screened ticks collected in seven regions of Russia and found that the abundant jingmenvirus-positive were of Ixodes ricinus species, with the prevalence ranging from 19.8% to 34.3%. In all cases, DNase/RNase treatment suggested that the detected molecule was DNA and subsequent next generation sequencing (NGS) proved that the viral polymerase gene was integrated in the I. ricinus genome. The copy number of the integrated polymerase gene was quantified by qPCR relative to the ITS2 gene and estimated as 1.32 copies per cell. At least three different genetic variants of the integrated polymerase gene were found in the territory of Russia. Phylogenetic analysis of the integrated jingmenvirus polymerase gene showed the highest similarity with the sequence of the correspondent gene obtained in Serbia from I. ricinus.

Mycovirus Diversity and Evolution Revealed/Inferred from Recent Studies

High-throughput virome analyses with various fungi, from cultured or uncultured sources, have led to the discovery of diverse viruses with unique genome structures and even neo-lifestyles. Examples in the former category include splipalmiviruses and ambiviruses. Splipalmiviruses, related to yeast narnaviruses, have multiple positive-sense (+) single-stranded (ss) RNA genomic segments that separately encode the RNA-dependent RNA polymerase motifs, the hallmark of RNA viruses (members of the kingdom Orthornavirae). Ambiviruses appear to have an undivided ssRNA genome of 3∼5 kb with two large open reading frames (ORFs) separated by intergenic regions. Another narna-like virus group has two fully overlapping ORFs on both strands of a genomic segment that span more than 90% of the genome size. New virus lifestyles exhibited by mycoviruses include the yado-kari/yado-nushi nature characterized by the partnership between the (+)ssRNA yadokarivirus and an unrelated dsRNA virus (donor of the capsid for the former) and the hadaka nature of capsidless 10-11 segmented (+)ssRNA accessible by RNase in infected mycelial homogenates. Furthermore, dsRNA polymycoviruses with phylogenetic affinity to (+)ssRNA animal caliciviruses have been shown to be infectious as dsRNA-protein complexes or deproteinized naked dsRNA. Many previous phylogenetic gaps have been filled by recently discovered fungal and other viruses, which haveprovided interesting evolutionary insights. Phylogenetic analyses and the discovery of natural and experimental cross-kingdom infections suggest that horizontal virus transfer may have occurred and continue to occur between fungi and other kingdoms.

PIWI Proteins Play an Antiviral Role in Lepidopteran Cell Lines

Insect antiviral immunity primarily relies on RNAi mechanisms. While a key role of small interfering (si)RNAs and AGO proteins has been well established in this regard, the situation for PIWI proteins and PIWI-interacting (pi)RNAs is not as clear. In the present study, we investigate whether PIWI proteins and viral piRNAs are involved in the immunity against single-stranded RNA viruses in lepidopteran cells, where two PIWIs are identified (Siwi and Ago3). Via loss- and gain-of-function studies in Bombyx mori BmN4 cells and in Trichoplusia ni High Five cells, we demonstrated an antiviral role of Siwi and Ago3. However, small RNA analysis suggests that viral piRNAs can be absent in these lepidopteran cells. Together with the current literature, our results support a functional diversification of PIWI proteins in insects.

High-quality, chromosome-scale genome assemblies: comparisons of three Diaphorina citri (Asian citrus psyllid) geographic populations

The Asian citrus psyllid, Diaphorina citri, is the insect vector of the causal agent of huanglongbing (HLB), a devastating bacterial disease of commercial citrus. Presently, few genomic resources exist for D. citri. In this study, we utilized PacBio HiFi and chromatin confirmation contact (Hi-C) sequencing to sequence, assemble, and compare three high-quality, chromosome-scale genome assemblies of D. citri collected from California, Taiwan, and Uruguay. Our assemblies had final sizes of 282.67 Mb (California), 282.89 Mb (Taiwan), and 266.67 Mb (Uruguay) assembled into 13 pseudomolecules—a reduction in assembly size of 41–45% compared with previous assemblies which we validated using flow cytometry. We identified the X chromosome in D. citri and annotated each assembly for repetitive elements, protein-coding genes, transfer RNAs, ribosomal RNAs, piwi-interacting RNA clusters, and endogenous viral elements. Between 19,083 and 20,357 protein-coding genes were predicted. Repetitive DNA accounts for 36.87–38.26% of each assembly. Comparative analyses and mitochondrial haplotype networks suggest that Taiwan and Uruguay D. citri are more closely related, while California D. citri are closely related to Florida D. citri. These high-quality, chromosome-scale assemblies provide new genomic resources to researchers to further D. citri and HLB research.

First Evidence of Past and Present Interactions between Viruses and the Black Soldier Fly, Hermetia illucens

Black soldier flies (BSFs, Hermetia illucens) are becoming a prominent research model encouraged by the insect as food and feed and waste bioconversion industries. Insect mass-rearing facilities are at risk from the spread of viruses, but so far, none have been described in BSFs. To fill this knowledge gap, a bioinformatic approach was undertaken to discover viruses specifically associated with BSFs. First, BSF genomes were screened for the presence of endogenous viral elements (EVEs). This led to the discovery and mapping of seven orthologous EVEs integrated into three BSF genomes originating from five viral families. Secondly, a virus discovery pipeline was used to screen BSF transcriptomes. This led to detecting a new exogenous totivirus that we named hermetia illucens totivirus 1 (HiTV1). Phylogenetic analyses showed this virus belongs to a clade of insect-specific totiviruses and is closely related to the largest EVE located on chromosome 1 of the BSF genome. Lastly, this EVE was found to express a small transcript in some BSFs infected by HiTV1. Altogether, this data mining study showed that far from being unscathed from viruses, BSFs bear traces of past interactions with several viral families and of present interactions with the exogenous HiTV1.

Transposable element repression using piRNAs, and its relevance to endogenous viral elements (EVEs) and immunity in insects

The piRNA system controls transposable element (TE) mobility by transcriptional gene silencing and post-transcriptional gene silencing. Dispersed in insect genomes, piRNA clusters contain TE copies, from which they produce piRNAs (specific small RNAs). These piRNAs can both target the nascent transcripts produced by active TE copies and directly repress them by heterochromatinization. They can also target mature transcripts and cleave them following amplification by the so-called 'ping-pong' loop mechanism. Moreover, piRNA clusters contain endogenous viral elements (EVEs), from which they produce piRNAs. The current idea is that these piRNAs could participate in the antiviral response against exogenous viral infection. In this review, we show that among insects, to date, this antiviral response by the piRNA system appears mainly restricted to mosquitoes, but this could be due to the focus of most studies on arboviruses.

The diversity of endogenous viral elements in insects

We provide an overview of the currently known diversity of viral sequences integrated into insect genomes. Such endogenous viral elements (EVE) have so far been annotated in at least eight insect orders and can be assigned to at least three families of large double-stranded (ds) DNA viruses, at least 22 families of RNA viruses, and three families of single-stranded DNA viruses. The study of these EVE has already produced important insights into insect-virus interactions, including the discovery of a new form of adaptive antiviral immunity. Insect EVE diversity will continue to increase as new insect genomes and exogenous viruses are sequenced, which will continue to make paleovirology a vibrant research field in this group of animals in the years to come.

Innate immune memory in invertebrates: Concept and potential mechanisms

Invertebrates are the protagonists of a recent paradigm shift because they now show that vertebrates are not the only group with immune memory. This review discusses the concept of immune priming, its characteristics, and differences with trained immunity and immune enhancement. We include an update of the current status of immune priming within generations in different groups of invertebrates which now include work in 5 Phyla: Ctenophora, Cnidaria, Mollusca, Nematoda, and Arthropoda. Clearly, few Phyla have been studied. We also resume and discuss the effector mechanism related to immune memory, including integrating viral elements into the genome, endoreplication, and epigenetics. The roles of other elements are incorporated, such as hemocytes, immune pathways, and metabolisms. We conclude that taking care of the experimental procedure will discern if results provide or do not support the invertebrates' immune memory and that regarding mechanisms, indeed, there are no studies on the immune memory mechanisms, this is how specificity is reached, and how and where the immune memory is stored and how is recall upon subsequent encounters. Finally, we discuss the possibility of having more than one mechanism working in different groups of invertebrates depending on the environmental conditions.

Virus and endogenous viral element-derived small non-coding RNAs and their roles in insect-virus interaction

RNA interference pathways mediated by different types of small non-coding RNAs (siRNAs, miRNAs and piRNAs) are conserved biological responses to exotic stresses, including viral infection. Aside from the well-established siRNA pathway, the miRNA pathway and the piRNA pathway process viral sequences, exogenously or endogenously, into miRNAs and piRNAs, respectively. During the host-virus interaction, viral sequences, including both coding and non-coding sequences, can be integrated as endogenous viral elements (EVEs) and thereby become present within the germline of a non-viral organism. In recent years, significant progress has been made in characterizing the biogenesis and function of viruses and EVEs associated with snRNAs. Overall, the siRNA pathway acts as the primarily antiviral defense against a wide range of exogenous viruses; the miRNA pathways associated with viruses or EVEs function in antiviral response and host gene regulation; EVE derived piRNAs with a ping-pong signature have the potential to limit cognate viral infection.

RNA virus EVEs in insect genomes

Insects are infected by a diverse set of RNA viruses that are more broadly distinguished by their ability to infect single or multiple host species. During replication into the host cell, partial or complete double strand DNA derived from the viral genome may be integrated into their host genomes giving origin to endogenous viral elements (EVEs). EVEs from RNA viruses have been identified in a variety of insect genomes showing different evolutionary trajectories: from highly degraded viral genomic remains to partial and complete viral coding regions. Limited functional knowledge exists about RNA EVEs impact on hosts and circulating viruses, but exciting results are emerging showing a complex arms race interplay that influences the evolutionary trajectory of these interacting entities.

Endogenous viral elements in mosquito genomes: current knowledge and outstanding questions

Integrations from non-retroviral RNA viruses (nrEVEs) have been identified across several taxa, including mosquitoes. Amongst all Culicinae species, the viral vectors Aedes aegypti and Aedes albopictus stand out for their high number of nrEVEs. In addition, Aedes nrEVEs are enriched in piRNA clusters and generate piRNAs that can silence incoming viral genomes. As such, nrEVEs represent a new form of inherited antiviral immunity. To propel this discovery into novel transmission-blocking vector control strategies, a deeper understanding of nrEVE biology and evolution is essential because differences in the landscape of nrEVEs have been identified in wild-caught mosquitoes, the piRNA profile of nrEVEs is not homogeneous and nrEVEs outside piRNA clusters exist and are expressed at the mRNA level. Here we summarise current knowledge on nrEVEs in mosquitoes and we point out the many unanswered questions and potentials of these genomic elements.

Computational Methods for the Discovery and Annotation of Viral Integrations

The transfer of genetic material between viruses and eukaryotic cells is pervasive. Somatic integrations of DNA viruses and retroviruses have been linked to persistent viral infection and genotoxic effects. Integrations into germline cells, referred to as Endogenous Viral Elements (EVEs), can be co-opted for host functions. Besides DNA viruses and retroviruses, EVEs can also derive from nonretroviral RNA viruses, which have often been observed in piRNA clusters. Here, we describe a bioinformatic framework to annotate EVEs in a genome assembly, study their widespread occurrence and polymorphism and identify sample-specific viral integrations using whole genome sequencing data.

Vector Specificity of Arbovirus Transmission

More than 25% of human infectious diseases are vector-borne diseases (VBDs). These diseases, caused by pathogens shared between animals and humans, are a growing threat to global health with more than 2.5 million annual deaths. Mosquitoes and ticks are the main vectors of arboviruses including flaviviruses, which greatly affect humans. However, all tick or mosquito species are not able to transmit all viruses, suggesting important molecular mechanisms regulating viral infection, dissemination, and transmission by vectors. Despite the large distribution of arthropods (mosquitoes and ticks) and arboviruses, only a few pairings of arthropods (family, genus, and population) and viruses (family, genus, and genotype) successfully transmit. Here, we review the factors that might limit pathogen transmission: internal (vector genetics, immune responses, microbiome including insect-specific viruses, and coinfections) and external, either biotic (adult and larvae nutrition) or abiotic (temperature, chemicals, and altitude). This review will demonstrate the dynamic nature and complexity of virus–vector interactions to help in designing appropriate practices in surveillance and prevention to reduce VBD threats.

Exosomes/microvesicles target SARS-CoV-2 via innate and RNA-induced immunity with PIWI-piRNA system

Murine neural stem cell exosomes/microvesicles can work to reduce SARS-CoV-2, an effect that can be adaptively enhanced via viral RNA fragment stimulation, which requires the PIWI-piRNA system.

Murine neural stem cells (NSCs) were recently shown to release piRNA-containing exosomes/microvesicles (Ex/Mv) for exerting antiviral immunity, but it remains unknown if these Ex/Mv could target SARS-CoV-2 and whether the PIWI-piRNA system is important for these antiviral actions. Here, using in vitro infection models, we show that hypothalamic NSCs (htNSCs) Ex/Mv provided an innate immunity protection against SARS-CoV-2. Importantly, enhanced antiviral actions were achieved by using induced Ex/Mv that were derived from induced htNSCs through twice being exposed to several RNA fragments of SARS-CoV-2 genome, a process that was designed not to involve protein translation of these RNA fragments. The increased antiviral effects of these induced Ex/Mv were associated with increased expression of piRNA species some of which could predictably target SARS-CoV-2 genome. Knockout of piRNA-interacting protein PIWIL2 in htNSCs led to reductions in both innate and induced antiviral effects of Ex/Mv in targeting SARS-CoV-2. Taken together, this study demonstrates a case suggesting Ex/Mv from certain cell types have innate and adaptive immunity against SARS-CoV-2, and the PIWI-piRNA system is important for these antiviral actions.

Persistent Spodoptera frugiperda rhabdovirus infection in Sf9 cells is not restricted by Wolbachia wMelPop-CLA and wAlbB strains and is targeted by the RNAi machinery

The endosymbiotic bacterium Wolbachia pipientis confers RNA virus refractoriness in Drosophila and Aedes mosquitoes. Questions remain about the Wolbachia-virus restriction phenotype and how extensive this phenomenon may be within other arthropods. Here, we generated two Spodoptera frugiperda cell lines stably transinfected with two strains of Wolbachia, wAlbB and wMelPop-CLA. Despite the high density of Wolbachia in stably infected Sf9 cells, RT-PCR indicated the presence of the negative-sense RNA virus Spodoptera frugiperda rhabdovirus (SfRV) in Wolbachia-infected and uninfected cell lines. No differences in the replication of SfRV between Sf9 and Wolbachia-infected cells was found. RNA-Seq analysis of the parental Sf9 cells supported SfRV's presence in these cells with abundant 20 nt virus-derived small RNAs indicating active replication of SfRV in these cells. Overall, this study supports a growing body of evidence that Wolbachia does not restrict negative-sense RNA viruses and generates an in vitro model to examine Lepidoptera-Wolbachia virus interactions.

Paleovirology of the DNA viruses of eukaryotes

Paleovirology is the study of ancient viruses and how they have coevolved with their hosts. An increasingly detailed understanding of the diversity, origins, and evolution of the DNA viruses of eukaryotes has been obtained through the lens of paleovirology in recent years. Members of multiple viral families have been found integrated in the genomes of eukaryotes, providing a rich fossil record to study. These elements have extended our knowledge of exogenous viral diversity, host ranges, and the timing of viral evolution, and are revealing the existence of entire new families of eukaryotic integrating dsDNA viruses and transposons. Future work in paleovirology will continue to provide insights into antiviral immunity, viral diversity, and potential applications, and reveal other secrets of the viral world.

RNAi-based immunity in insects against baculoviruses and the strategies of baculoviruses involved in siRNA and miRNA pathways to weaken the defense

Protection against viral infection in hosts concerns diverse cellular and molecular mechanisms, among which RNA interference (RNAi) response is a vital one. Small interfering RNAs (siRNAs), microRNAs (miRNAs) and PIWI interacting RNAs (piRNAs) are primary categories of small RNAs involved in RNAi response, playing significant roles in restraining viral invasion. However, during a long-term coevolution, viruses have gained the ability to evade, avoid, or suppress antiviral immunity to ensure efficient replication and transmission. Baculoviruses are enveloped, insect-pathogenic viruses with double-stranded circular DNA genomes, which encode suppressors of siRNA pathway and miRNAs targeting immune-related genes to mask the antiviral activity of their hosts. This review summarized recent findings for the RNAi-based antiviral immunity in insects as well as the strategies that baculoviruses exploit to break the shield of host siRNA pathway, and hijack cellular miRNAs or encode their own miRNAs that regulate both viral and cellular gene expression to create a favorable environment for viral infection.

Uncovering the Worldwide Diversity and Evolution of the Virome of the Mosquitoes Aedes aegypti and Aedes albopictus

Aedes aegypti, the yellow fever mosquito, and Aedes albopictus, the Asian tiger mosquito, are the most significant vectors of dengue, Zika, and Chikungunya viruses globally. Studies examining host factors that control arbovirus transmission demonstrate that insect-specific viruses (ISVs) can modulate mosquitoes’ susceptibility to arbovirus infection in both in vivo and in vitro co-infection models. While research is ongoing to implicate individual ISVs as proviral or antiviral factors, we have a limited understanding of the composition and diversity of the Aedes virome. To address this gap, we used a meta-analysis approach to uncover virome diversity by analysing ~3000 available RNA sequencing libraries representing a worldwide geographic range for both mosquitoes. We identified ten novel viruses and previously characterised viruses, including mononegaviruses, orthomyxoviruses, negeviruses, and a novel bi-segmented negev-like group. Phylogenetic analysis suggests close relatedness to mosquito viruses implying likely insect host range except for one arbovirus, the multi-segmented Jingmen tick virus (Flaviviridae) in an Italian colony of Ae. albopictus. Individual mosquito transcriptomes revealed remarkable inter-host variation of ISVs within individuals from the same colony and heterogeneity between different laboratory strains. Additionally, we identified striking virus diversity in Wolbachia infected Aedes cell lines. This study expands our understanding of the virome of these important vectors. It provides a resource for further assessing the ecology, evolution, and interaction of ISVs with their mosquito hosts and the arboviruses they transmit.

Nudivirus Remnants in the Genomes of Arthropods

Endogenous viral elements (EVEs), derived from all major types of viruses, have been discovered in many eukaryotic genomes, representing "fossil records" of past viral infections. The endogenization of nudiviruses has been reported in several insects, leading to the question of whether genomic integration is a common phenomenon for these viruses. In this study, genomic assemblies of insects and other arthropods were analyzed to identify endogenous sequences related to Nudiviridae. A total of 359 nudivirus-like genes were identified in 43 species belonging to different groups; however, none of these genes were detected in the known hosts of nudiviruses. A large proportion of the putative EVEs identified in this study encode intact open reading frames or are transcribed as mRNAs, suggesting that they result from recent endogenization of nudiviruses. Phylogenetic analyses of the identified EVEs and inspections of their flanking regions indicated that integration of nudiviruses has occurred recurrently during the evolution of arthropods. This is the first report of a comprehensive screening for nudivirus-derived EVEs in arthropod genomes. The results of this study demonstrated that a large variety of arthropods, especially hemipteran and hymenopteran insects, have previously been or are still infected by nudiviruses. These findings have greatly extended the host range of Nudiviridae and provide new insights into viral diversity, evolution, and host-virus interactions.

Flavivirus integrations in Aedes aegypti are limited and highly conserved across samples from different geographic regions unlike integrations in Aedes albopictus

Mosquitoes of the genus Aedes are the main vectors of many viruses, e.g. dengue and Zika, which affect millions of people each year and for which there are limited treatment options. Understanding how Aedes mosquitoes tolerate high viral loads may lead to better disease control strategies. Elucidating endogenous viral elements (EVEs) within vector genomes may give exploitable biological insights. Previous studies have reported the presence of a large number of EVEs in Aedes genomes. Here we investigated if flavivirus EVEs are conserved across populations and different Aedes species by using ~ 500 whole genome sequence libraries from Aedes aegypti and Aedes albopictus, sourced from colonies and field mosquitoes across continents. We found that nearly all flavivirus EVEs in the Ae. aegypti reference genome originate from four separate putative viral integration events, and that they are highly conserved across geographically diverse samples. By contrast, flavivirus EVEs in the Ae. albopictus reference genome originate from up to nine distinct integration events and show low levels of conservation, even within samples from narrow geographical ranges. Our analysis suggests that flaviviruses integrated as long sequences and were subsequently fragmented and shuffled by transposable elements. Given that EVEs of Ae. aegypti and Ae. albopictus belong to different phylogenetic clades and have very differing levels of conservation, they may have different evolutionary origins and potentially different functional roles.

Novel RNA Viruses from the Transcriptome of Pheromone Glands in the Pink Bollworm Moth, Pectinophora gossypiella

Simple Summary

The pink bollworm, Pectinophora gossypiella (Lepidoptera: Gelechiidae), is a major pest of cotton. In this study, we analyzed the mRNA from pheromone glands of two populations in Israel. We found several virus sequences that were the same in these populations. We identified these viruses based on high-throughput sequencing data and analysis of the assembled transcripts. Through analysis of the sequences, we identified several unique viral sequences representing possible novel viral species. Two of the viral sequences were found in relatively high abundance in pheromone glands. One of the virus sequences was also found through analysis of previous transcriptome sequencing data from the midgut of pink bollworm larvae. This is the first report of these unique viral sequences found in the pink bollworm, and these viruses could be developed to help control this pest around the world, but more research is needed to determine their utility as biological control agents.

Abstract

In this study, we analyzed the transcriptome obtained from the pheromone gland isolated from two Israeli populations of the pink bollworm Pectinophora gossypiella to identify viral sequences. The lab population and the field samples carried the same viral sequences. We discovered four novel viruses: two positive-sense single-stranded RNA viruses, Pectinophora gossypiella virus 1 (PecgV1, a virus of Iflaviridae) and Pectinophora gossypiella virus 4 (PecgV4, unclassified), and two negative-sense single-stranded RNA viruses, Pectinophora gossypiella virus 2 (PecgV2, a virus of Phasmaviridae) and Pectinophora gossypiella virus 3 (PecgV3, a virus of Phenuiviridae). In addition, sequences derived from two negative-sense single-stranded RNA viruses that belong to Mononegavirales were found in the data. Analysis of previous transcriptome sequencing data derived from the midgut of pink bollworm larvae of a USA population only identified PecgV1, but no other viruses. High viral sequence coverages of PecgV1 and PecgV4 were observed in both field and lab populations. This is the first report of viral sequences discovered from the pink bollworm. Results from this investigation suggest that the pink bollworm harbors multiple viruses. Further investigation of the viral pathogens may help to develop novel pest management strategies for control of the pink bollworm.

Persistence of Ambigrammatic Narnaviruses Requires Translation of the Reverse Open Reading Frame

Narnaviruses are RNA viruses detected in diverse fungi, plants, protists, arthropods, and nematodes. Though initially described as simple single-gene nonsegmented viruses encoding RNA-dependent RNA polymerase (RdRp), a subset of narnaviruses referred to as "ambigrammatic" harbor a unique genomic configuration consisting of overlapping open reading frames (ORFs) encoded on opposite strands. Phylogenetic analysis supports selection to maintain this unusual genome organization, but functional investigations are lacking. Here, we establish the mosquito-infecting Culex narnavirus 1 (CxNV1) as a model to investigate the functional role of overlapping ORFs in narnavirus replication. In CxNV1, a reverse ORF without homology to known proteins covers nearly the entire 3.2-kb segment encoding the RdRp. Additionally, two opposing and nearly completely overlapping novel ORFs are found on the second putative CxNV1 segment, the 0.8-kb "Robin" RNA. We developed a system to launch CxNV1 in a naive mosquito cell line and then showed that functional RdRp is required for persistence of both segments, and an intact reverse ORF is required on the RdRp segment for persistence. Mass spectrometry of persistently CxNV1-infected cells provided evidence for translation of this reverse ORF. Finally, ribosome profiling yielded a striking pattern of footprints for all four CxNV1 RNA strands that was distinct from actively translating ribosomes on host mRNA or coinfecting RNA viruses. Taken together, these data raise the possibility that the process of translation itself is important for persistence of ambigrammatic narnaviruses, potentially by protecting viral RNA with ribosomes, thus suggesting a heretofore undescribed viral tactic for replication and transmission. IMPORTANCE Fundamental to our understanding of RNA viruses is a description of which strand(s) of RNA are transmitted as the viral genome relative to which encode the viral proteins. Ambigrammatic narnaviruses break the mold. These viruses, found broadly in fungi, plants, and insects, have the unique feature of two overlapping genes encoded on opposite strands, comprising nearly the full length of the viral genome. Such extensive overlap is not seen in other RNA viruses and comes at the cost of reduced evolutionary flexibility in the sequence. The present study is motivated by investigating the benefits which balance that cost. We show for the first time a functional requirement for the ambigrammatic genome configuration in Culex narnavirus 1, which suggests a model for how translation of both strands might benefit this virus. Our work highlights a new blueprint for viral persistence, distinct from strategies defined by canonical definitions of the coding strand.

Population genomics in the arboviral vector Aedes aegypti reveals the genomic architecture and evolution of endogenous viral elements

Horizontal gene transfer from viruses to eukaryotic cells is a pervasive phenomenon. Somatic viral integrations are linked to persistent viral infection whereas integrations into germline cells are maintained in host genomes by vertical transmission and may be co-opted for host functions. In the arboviral vector Aedes aegypti, an endogenous viral element from a nonretroviral RNA virus (nrEVE) was shown to produce PIWI-interacting RNAs (piRNAs) to limit infection with a cognate virus. Thus, nrEVEs may constitute a heritable, sequence-specific mechanism for antiviral immunity, analogous to piRNA-mediated silencing of transposable elements. Here, we combine population genomics and evolutionary approaches to analyse the genomic architecture of nrEVEs in A. aegypti. We conducted a genome-wide screen for adaptive nrEVEs and searched for novel population-specific nrEVEs in the genomes of 80 individual wild-caught mosquitoes from five geographical populations. We show a dynamic landscape of nrEVEs in mosquito genomes and identified five novel nrEVEs derived from two currently circulating viruses, providing evidence of the environmental-dependent modification of a piRNA cluster. Overall, our results show that virus endogenization events are complex with only a few nrEVEs contributing to adaptive evolution in A. aegypti.

Profile of Small RNAs, vDNA Forms and Viral Integrations in Late Chikungunya Virus Infection of Aedes albopictus Mosquitoes

The Asian tiger mosquito Aedes albopictus is contributing to the (re)-emergence of Chikungunya virus (CHIKV). To gain insights into the molecular underpinning of viral persistence, which renders a mosquito a life-long vector, we coupled small RNA and whole genome sequencing approaches on carcasses and ovaries of mosquitoes sampled 14 days post CHIKV infection and investigated the profile of small RNAs and the presence of vDNA fragments. Since Aedes genomes harbor nonretroviral Endogenous Viral Elements (nrEVEs) which confers tolerance to cognate viral infections in ovaries, we also tested whether nrEVEs are formed after CHIKV infection. We show that while small interfering (si)RNAs are evenly distributed along the full viral genome, PIWI-interacting (pi)RNAs mostly arise from a ~1000 bp window, from which a unique vDNA fragment is identified. CHIKV infection does not result in the formation of new nrEVEs, but piRNAs derived from existing nrEVEs correlate with differential expression of an endogenous transcript. These results demonstrate that all three RNAi pathways contribute to the homeostasis during the late stage of CHIKV infection, but in different ways, ranging from directly targeting the viral sequence to regulating the expression of mosquito transcripts and expand the role of nrEVEs beyond immunity against cognate viruses.

piRNAs as Modulators of Disease Pathogenesis

Advances in understanding disease pathogenesis correlates to modifications in gene expression within different tissues and organ systems. In depth knowledge about the dysregulation of gene expression profiles is fundamental to fully uncover mechanisms in disease development and changes in host homeostasis. The body of knowledge surrounding mammalian regulatory elements, specifically regulators of chromatin structure, transcriptional and translational activation, has considerably surged within the past decade. A set of key regulators whose function still needs to be fully elucidated are small non-coding RNAs (sncRNAs). Due to their broad range of unfolding functions in the regulation of gene expression during transcription and translation, sncRNAs are becoming vital to many cellular processes. Within the past decade, a novel class of sncRNAs called PIWI-interacting RNAs (piRNAs) have been implicated in various diseases, and understanding their complete function is of vital importance. Historically, piRNAs have been shown to be indispensable in germline integrity and stem cell development. Accumulating research evidence continue to reveal the many arms of piRNA function. Although piRNA function and biogenesis has been extensively studied in Drosophila, it is thought that they play similar roles in vertebrate species, including humans. Compounding evidence suggests that piRNAs encompass a wider functional range than small interfering RNAs (siRNAs) and microRNAs (miRNAs), which have been studied more in terms of cellular homeostasis and disease. This review aims to summarize contemporary knowledge regarding biogenesis, and homeostatic function of piRNAs and their emerging roles in the development of pathologies related to cardiomyopathies, cancer, and infectious diseases.

ViR: a tool to solve intrasample variability in the prediction of viral integration sites using whole genome sequencing data

Background

Several bioinformatics pipelines have been developed to detect sequences from viruses that integrate into the human genome because of the health relevance of these integrations, such as in the persistence of viral infection and/or in generating genotoxic effects, often progressing into cancer. Recent genomics and metagenomics analyses have shown that viruses also integrate into the genome of non-model organisms (i.e., arthropods, fish, plants, vertebrates). However, rarely studies of endogenous viral elements (EVEs) in non-model organisms have gone beyond their characterization from reference genome assemblies. In non-model organisms, we lack a thorough understanding of the widespread occurrence of EVEs and their biological relevance, apart from sporadic cases which nevertheless point to significant roles of EVEs in immunity and regulation of expression. The concomitance of repetitive DNA, duplications and/or assembly fragmentations in a genome sequence and intrasample variability in whole-genome sequencing (WGS) data could determine misalignments when mapping data to a genome assembly. This phenomenon hinders our ability to properly identify integration sites.

Results

To fill this gap, we developed ViR, a pipeline which solves the dispersion of reads due to intrasample variability in sequencing data from both single and pooled DNA samples thus ameliorating the detection of integration sites. We tested ViR to work with both in silico and real sequencing data from a non-model organism, the arboviral vector Aedes albopictus. Potential viral integrations predicted by ViR were molecularly validated supporting the accuracy of ViR results.

Conclusion

ViR will open new venues to explore the biology of EVEs, especially in non-model organisms. Importantly, while we generated ViR with the identification of EVEs in mind, its application can be extended to detect any lateral transfer event providing an ad-hoc sequence to interrogate.

Leaning Into the Bite: The piRNA Pathway as an Exemplar for the Genetic Engineering Need in Mosquitoes

The piRNA pathway is a specialized small RNA interference that in mosquitoes is mechanistically distant from analogous biology in the Drosophila model. Current genetic engineering methods, such as targeted genome manipulation, have a high potential to tease out the functional complexity of this intricate molecular pathway. However, progress in utilizing these methods in arthropod vectors has been geared mostly toward the development of new vector control strategies rather than to study cellular functions. Herein we propose that genetic engineering methods will be essential to uncover the full functionality of PIWI/piRNA biology in mosquitoes and that extending the applications of genetic engineering on other aspects of mosquito biology will grant access to a much larger pool of knowledge in disease vectors that is just out of reach. We discuss motivations for and impediments to expanding the utility of genetic engineering to study the underlying biology and disease transmission and describe specific areas where efforts can be placed to achieve the full potential for genetic engineering in basic biology in mosquito vectors. Such efforts will generate a refreshed intellectual source of novel approaches to disease control and strong support for the effective use of approaches currently in development.

The Interplay Between Viruses and RNAi Pathways in Insects

As an overarching immune mechanism, RNA interference (RNAi) displays pathogen specificity and memory via different pathways. The small interfering RNA (siRNA) pathway is the primary antiviral defense mechanism against RNA viruses of insects and plays a lesser role in defense against DNA viruses. Reflecting the pivotal role of the siRNA pathway in virus selection, different virus families have independently evolved unique strategies to counter this host response, including protein-mediated, decoy RNA-based, and microRNA-based strategies. In this review, we outline the interplay between insect viruses and the different pathways of the RNAi antiviral response; describe practical application of these interactions for improved expression systems and for pest and disease management; and highlight research avenues for advancement of the field.

Viral fossils in marsupial genomes: secret cellular guardians

Genomic viral integrations, termed endogenous viral elements (EVEs), are fragments of viruses in host chromosomes that provide information about viral evolution and could even help protect the host from infection. In the present study we examined EVEs in thirteen different Australian marsupial species to identify trends in their integration, commonality and to investigate their possible cellular function. We found that marsupial EVEs are commonly derived from viruses of the Bornaviridae, Filoviridae and Parvoviridae families, and circulated up to 160 million years ago. We also show the EVEs are actively transcribed into both long and short RNA molecules in marsupials, and propose they are involved in a cellular defence mechanism to protect the germline from viral genomic invasion.

Diversity of Sea Star-Associated Densoviruses and Transcribed Endogenous Viral Elements of Densovirus Origin

The primary interest in sea star densoviruses, specifically SSaDV, has been their association with sea star wasting syndrome (SSWS), a disease that has decimated sea star populations across the West Coast of the United States since 2013. The association of SSaDV with SSWS was originally drawn from metagenomic analysis, which was further studied through field surveys using quantitative PCR (qPCR), with the conclusion that it was the most likely viral candidate in the metagenomic data based on its representation in symptomatic sea stars compared to asymptomatic sea stars. We reexamined the original metagenomic data with additional genomic data sets and found that SSaDV was 1 of 10 densoviruses present in the original data set and was no more represented in symptomatic sea stars than in asymptomatic sea stars. Instead, SSaDV appears to be a widespread, generalist virus that exists among a large diversity of densoviruses present in sea star populations.

A viral etiology of sea star wasting syndrome (SSWS) was originally explored with virus-sized material challenge experiments, field surveys, and metagenomics, leading to the conclusion that a densovirus is the predominant DNA virus associated with this syndrome and, thus, the most promising viral candidate pathogen. Single-stranded DNA viruses are, however, highly diverse and pervasive among eukaryotic organisms, which we hypothesize may confound the association between densoviruses and SSWS. To test this hypothesis and assess the association of densoviruses with SSWS, we compiled past metagenomic data with new metagenomic-derived viral genomes from sea stars collected from Antarctica, California, Washington, and Alaska. We used 179 publicly available sea star transcriptomes to complement our approaches for densovirus discovery. Lastly, we focus the study on sea star-associated densovirus (SSaDV), the first sea star densovirus discovered, by documenting its biogeography and putative tissue tropism. Transcriptomes contained only endogenized densovirus elements similar to the NS1 gene, while numerous extant densoviral genomes were recovered from viral metagenomes. SSaDV was associated with nearly all tested species from southern California to Alaska, and in contrast to previous work, we show that SSaDV is one genotype among a high diversity of densoviruses present in sea stars across the West Coast of the United States and globally that are commonly associated with grossly normal (i.e., healthy or asymptomatic) animals. The diversity and ubiquity of these viruses in sea stars confound the original hypothesis that one densovirus is the etiological agent of SSWS.

IMPORTANCE The primary interest in sea star densoviruses, specifically SSaDV, has been their association with sea star wasting syndrome (SSWS), a disease that has decimated sea star populations across the West Coast of the United States since 2013. The association of SSaDV with SSWS was originally drawn from metagenomic analysis, which was further studied through field surveys using quantitative PCR (qPCR), with the conclusion that it was the most likely viral candidate in the metagenomic data based on its representation in symptomatic sea stars compared to asymptomatic sea stars. We reexamined the original metagenomic data with additional genomic data sets and found that SSaDV was 1 of 10 densoviruses present in the original data set and was no more represented in symptomatic sea stars than in asymptomatic sea stars. Instead, SSaDV appears to be a widespread, generalist virus that exists among a large diversity of densoviruses present in sea star populations.

Soybean Thrips (Thysanoptera: Thripidae) Harbor Highly Diverse Populations of Arthropod, Fungal and Plant Viruses

Soybean thrips (Neohydatothrips variabilis) are one of the most efficient vectors of soybean vein necrosis virus, which can cause severe necrotic symptoms in sensitive soybean plants. To determine which other viruses are associated with soybean thrips, the metatranscriptome of soybean thrips, collected by the Midwest Suction Trap Network during 2018, was analyzed. Contigs assembled from the data revealed a remarkable diversity of virus-like sequences. Of the 181 virus-like sequences identified, 155 were novel and associated primarily with taxa of arthropod-infecting viruses, but sequences similar to plant and fungus-infecting viruses were also identified. The novel viruses were predicted to have positive-sense RNA, negative-stranded RNA, double-stranded RNA, and single-stranded DNA genomes. The assembled sequences included 100 contigs that represented at least 95% coverage of a virus genome or genome segment. Sequences represented 12 previously described arthropod viruses including eight viruses reported from Hubei Province in China, and 12 plant virus sequences of which six have been previously described. The presence of diverse populations of plant viruses within soybean thrips suggests they feed on and acquire viruses from multiple host plant species that could be transmitted to soybean. Assessment of the virome of soybean thrips provides, for the first time, information on the diversity of viruses present in thrips.

Endogenous Viral Element-Derived Piwi-Interacting RNAs (piRNAs) Are Not Required for Production of Ping-Pong-Dependent piRNAs from Diaphorina citri Densovirus

Piwi-interacting RNAs (piRNAs) are a class of small RNAs primarily responsible for silencing transposons in the animal germ line. The ping-pong cycle, the posttranscriptional silencing branch of the piRNA pathway, relies on piRNAs produced from endogenous transposon remnants to direct cleavage of transposon RNA via association with Piwi-family Argonaute proteins. In some mosquito species and mosquito-derived cell lines expressing a functionally expanded group of Piwi-family Argonaute proteins, both RNA and DNA viruses are targeted by piRNAs in a manner thought to involve direct processing of exogenous viral RNA into piRNAs. Whether viruses are targeted by piRNAs in nonmosquito species is unknown. Partial integrations of DNA and nonretroviral RNA virus genomes, termed endogenous viral elements (EVEs), are abundant in arthropod genomes and often produce piRNAs that are speculated to target cognate viruses through the ping-pong cycle. Here, we describe a Diaphorina citri densovirus (DcDV)-derived EVE in the genome of Diaphorina citri We found that this EVE gives rise to DcDV-specific primary piRNAs and is unevenly distributed among D. citri populations. Unexpectedly, we found that DcDV is targeted by ping-pong-dependent virus-derived piRNAs (vpiRNAs) in D. citri lacking the DcDV-derived EVE, while four naturally infecting RNA viruses of D. citri are not targeted by vpiRNAs. Furthermore, a recombinant Cricket paralysis virus containing a portion of the DcDV genome corresponding to the DcDV-derived EVE was not targeted by vpiRNAs during infection in D. citri harboring the EVE. These results demonstrate that viruses can be targeted by piRNAs in a nonmosquito species independently of endogenous piRNAs.IMPORTANCE Small RNAs serve as specificity determinants of antiviral responses in insects. Piwi-interacting RNAs (piRNAs) are a class of small RNAs found in animals, and their primary role is to direct antitransposon responses. These responses require endogenous piRNAs complementary to transposon RNA. Additionally, piRNAs have been shown to target RNA and DNA viruses in some mosquito species. In contrast to transposons, targeting of viruses by the piRNA pathway in these mosquito species does not require endogenous piRNAs. Here, we show that piRNAs target a DNA virus, but not RNA viruses, in an agricultural insect pest. We found that targeting of this DNA virus did not require endogenous piRNAs and that endogenous piRNAs did not mediate targeting of an RNA virus with which they shared complementary sequence. Our results highlight differences between mosquitoes and our experimental system and raise the possibility that DNA viruses may be targeted by piRNAs in other species.