Control of RNA viruses in mosquito cells through the acquisition of vDNA and endogenous viral elements

Molecular antiviral responses, immune priming and inheritance in insects

Viral diseases transmitted by insects to plants cause severe agricultural damage and arboviruses transmitted to humans cause severe disease outbreaks. The interaction between viruses and the insect defences is complex and has evolved into acting-counteracting molecular interplays. Viruses depict complex molecular mechanisms to ensure invasion, replication and exit the insect host cell, to invade other cells. On the other hand, insect cells use molecular strategies to recognize, halt replication and eliminate the invaders. In turn, virus counteract with evasive strategies. The main antiviral defence mechanism RNA interference (RNAi) recognizes and degrades viral RNA, thereby inhibiting viral replication. These in conjunction with other canonical immune pathways, Toll, IMD, JAK/STAT and Akt-ERK developed mainly to combat bacteria, fungi and protozoa, along with mechanisms to eliminate infected cells like apoptosis and phagocytosis comprise a multifactorial system. Insects exposed to an attenuated or sublethal viral infection could respond with faster and enhanced immune responses to the same pathogen (priming), which is like immunological memory in vertebrates. Several mechanisms have been proposed to explain priming, including endoreplication, epigenetic gene modifications by DNA methylation and histone acetylation. Priming could be inherited by the offspring (transgenerational immune priming, TGIP). However, the precise molecular mechanisms underlying TGIP remain to be elucidated. This article reviews the molecular mechanisms employed by insects to combat viral infections, discusses the current information and the outstanding research questions in the area.

Phasmavirus-derived genome sequences and endogenous viral element identified in the small hive beetle, Aethina tumida Murray

The small hive beetle (SHB), Aethina tumida Murray is an invasive pest of the honey bee. This beetle feeds not only on bee resources within the hive such as honey and pollen, but also on bee brood and dead bees. The impact of this beetle's intimate parasitic association with the honey bee on virus transmission is poorly understood. We aimed to characterize the virome of SHB to identify SHB viruses with potential for use in biological control of this pest. We characterized the virome of SHB by sequencing the transcriptomes and small RNAs of SHB collected from multiple geographical regions: Adult and larval SHB were collected from midwestern- (Illinois, Ohio) and southern- (Florida, Texas) states of the USA, and from South Africa. Small RNAs were sequenced for adult beetles from Florida and Ohio, for larvae from Florida, and for an SHB-derived cell line (BCIRL-AtumEN-1129). Assembled transcripts were annotated by BLASTx. In field-caught adult beetles and adults and larvae from South Africa, the near-complete sequences for all three genomic segments of a putative novel phasmavirus (order: Elliovirales, formerly Bunyavirales) were identified. In addition, transcripts from a partial glycoprotein sequence from a different phasmavirus integrated into the genome of SHB were detected in all samples, including the SHB-derived cell line. Apparent PIWI-interacting RNAs derived from the integrated glycoprotein sequence were also detected. Whether the putative extant phasmavirus replicates in SHB remains to be determined.

EEfinder, a general purpose tool for identification of bacterial and viral endogenized elements in eukaryotic genomes

Horizontal gene transfer is a phenomenon of genetic material transmission between species with no parental relationship. It has been characterized among several major branches of life, including among prokaryotes, viruses and eukaryotes. The characterization of endogenous elements derived from viruses or bacteria provides a snapshot of past host-pathogen interactions and coevolution as well as reference information to remove false positive results from metagenomic studies. Currently there is a lack of general purpose standardized tools for endogenous elements screening which limits reproducibility and hinder comparative analysis between studies. Here we describe EEfinder, a new general purpose tool for identification and classification of endogenous elements derived from viruses or bacteria found in eukaryotic genomes. The tool was developed to include six common steps performed in this type of analysis: data cleaning, similarity search through sequence alignment, filtering candidate elements, taxonomy assignment, merging of truncated elements and flanks extraction. We evaluated the sensitivity of EEfinder to identify endogenous elements through comparative analysis using data from the literature and showed that EEfinder automatically detected 97 % of the EVEs compared to published results obtained by manual curation and detected an almost exact full integration of a Wolbachia genome described using wet-lab experiments. Therefore, EEfinder can effectively and systematically identify endogenous elements with bacterial/viral origin integrated in eukaryotic genomes. EEfinder is publicly available on https://github.com/WallauBioinfo/EEfinder.

Evidence of endogenous non-retroviral RNA virus sequences into the genome and transcriptome of the malaria vector Anopheles darlingi

The characterization of non-retroviral integrated RNA virus sequences (NIRVS) in mosquitoes has emerged as a significant area of research that could yield insight into virus-host interactions. This study aimed to characterize NIRVS in the Anopheles darlingi reference genome and identify putative transcribed NIRVS in field-collected mosquitoes from Colombia. The An. darlingi reference genome was analyzed to identify and characterize NIRVS by conducting a BLAST query with all the virus sequences previously identified in arthropods available in the NCBI-virus repository. In addition, An. darlingi field-collected mosquitoes were examined for NIRVS using a metatranscriptomic approach. As a result, 44 NIRVS were identified in the An. darlingi genome, constituting integrations of negative single-stranded RNA viruses (ssRNA-) from the families Rhabdoviridae, Chuviridae and Phasmaviridae, and integrations of double-stranded RNA viruses (dsRNA) from the families Partitiviridae and Sedoreoviridae. These NIRVS were not randomly distributed but clustered in specific regions of the genome enriched with BEL/Pao and Ty3/Gypsy long terminal repeat elements. Furthermore, putative NIRVS-like sequences were present in the transcriptomic data from all the Colombian An. darlingi natural populations. This study is significant as it represents the first identification of NIRVS in the most important malaria vector of the Neotropics. The findings help in understanding the intricate relationship between the mosquito and its virome, and the regulation of viruses' mechanisms in the Anopheles genus.

Endogenous viral elements: insights into data availability and accessibility

Endogenous viral elements (EVEs) are remnants of viral genetic material endogenized into the host genome. They have, in the last decades, attracted attention for their role as potential contributors to pathogenesis, drivers of selective advantage for the host, and genomic remnants of ancient viruses. EVEs have a nuanced and complex influence on both host health and evolution, and can offer insights on the deep evolutionary history of viruses. As an emerging field of research, several factors limit a comprehensive understanding of EVEs: they are currently underestimated and periodically overlooked in studies of the host genome, transcriptome, and virome. The absence of standardized guidelines for ensuring EVE-related data availability and accessibility following the FAIR ('findable, accessible, interoperable, and reusable') principles obstructs our ability to gather and connect information. Here, we discuss challenges to the availability and accessibility of EVE-related data and propose potential solutions. We identified the biological and research focus imbalance between different types of EVEs, and their overall biological complexity as genomic loci with viral ancestry, as potential challenges that can be addressed with the development of a user-oriented identification tool. In addition, reports of EVE identification are scattered between different subfields under different keywords, and EVE sequences and associated data are not properly gathered in databases. While developing an open and dedicated database might be ideal, targeted improvements of generalist databases might provide a pragmatic solution to EVE data and metadata accessibility. The implementation of these solutions, as well as the collective effort by the EVE scientific community in discussing and setting guidelines, is now drastically needed to lead the development of EVE research and offer insights into host-virus interactions and their evolutionary history.

Pan-flavivirus analysis reveals sfRNA-independent, 3' UTR-biased siRNA production from an insect-specific flavivirus

RNA interference (RNAi) plays an essential role in mosquito antiviral immunity, but it is not known whether viral small interfering RNA (siRNA) profiles differ between mosquito-borne and mosquito-specific viruses. A pan-Orthoflavivirus analysis in Aedes albopictus cells revealed that viral siRNAs were evenly distributed across the viral genome of most representatives of the Flavivirus genus. In contrast, siRNA production was biased toward the 3' untranslated region (UTR) of the genomes of classical insect-specific flaviviruses (cISF), which was most pronounced for Kamiti River virus (KRV), a virus with a unique, 1.2 kb long 3' UTR. KRV-derived siRNAs were produced in high quantities and almost exclusively mapped to the 3' UTR. We mapped the 5' end of KRV subgenomic flavivirus RNAs (sfRNAs), products of the 5'-3' exoribonuclease XRN1/Pacman stalling on secondary RNA structures in the 3' UTR of the viral genome. We found that KRV produces high copy numbers of a long, 1,017 nt sfRNA1 and a short, 421 nt sfRNA2, corresponding to two predicted XRN1-resistant elements. Expression of both sfRNA1 and sfRNA2 was reduced in Pacman-deficient Aedes albopictus cells; however, this did not correlate with a shift in viral siRNA profiles. We suggest that cISFs, particularly KRV, developed a unique mechanism to produce high amounts of siRNAs as a decoy for the antiviral RNAi response in an sfRNA-independent manner.IMPORTANCEThe Flavivirus genus contains diverse mosquito viruses ranging from insect-specific viruses circulating exclusively in mosquito populations to mosquito-borne viruses that cause disease in humans and animals. Studying the mechanisms of virus replication and antiviral immunity in mosquitoes is important to understand arbovirus transmission and may inform the development of disease control strategies. In insects, RNA interference (RNAi) provides broad antiviral activity and constitutes a major immune response against viruses. Comparing diverse members of the Flavivirus genus, we found that all flaviviruses are targeted by RNAi. However, the insect-specific Kamiti River virus was unique in that small interfering RNAs are highly skewed toward its uniquely long 3' untranslated region. These results suggest that mosquito-specific viruses have evolved unique mechanisms for genome replication and immune evasion.

Endogenous nege-like viral elements in arthropod genomes reveal virus-host coevolution and ancient history of two plant virus families

Negevirus is a recently proposed taxon of arthropod-infecting virus, which is associated with plant viruses of two families (Virgaviridae and Kitaviridae). Nevertheless, the evolutionary history of negevirus-host and its relationship with plant viruses remain poorly understood. Endogenous nege-like viral elements (ENVEs) are ancient nege-like viral sequences integrated into the arthropod genomes, which can serve as the molecular fossil records of previous viral infection. In this study, 292 ENVEs were identified in 150 published arthropod genomes, revealing the evolutionary history of nege-like viruses and two related plant virus families. We discovered three novel and eight strains of nege-like viruses in 11 aphid species. Further analysis indicated that 10 ENVEs were detected in six aphid genomes, and they were divided into four types (ENVE1-ENVE4). Orthologous integration and phylogenetic analyses revealed that nege-like viruses had a history of infection of over 60 My and coexisted with aphid ancestors throughout the Cenozoic Era. Moreover, two nege-like viral proteins (CP and SP24) were highly homologous to those of plant viruses in the families Virgaviridae and Kitaviridae. CP- and SP24-derived ENVEs were widely integrated into numerous arthropod genomes. These results demonstrate that nege-like viruses have a long-term coexistence with arthropod hosts and plant viruses of the two families, Virgaviridae and Kitaviridae, which may have evolved from the nege-like virus ancestor through horizontal virus transfer events. These findings broaden our perspective on the history of viral infection in arthropods and the origins of plant viruses.

IMPORTANCE

Although negevirus is phylogenetically related to plant virus, the evolutionary history of negevirus-host and its relationship with plant virus remain largely unknown. In this study, we used endogenous nege-like viral elements (ENVEs) as the molecular fossil records to investigate the history of nege-like viral infection in arthropod hosts and the evolution of two related plant virus families (Virgaviridae and Kitaviridae). Our results showed the infection of nege-like viruses for over 60 My during the arthropod evolution. ENVEs highly homologous to viral sequences in Virgaviridae and Kitaviridae were present in a wide range of arthropod genomes but were absent in plant genomes, indicating that plant viruses in these two families possibly evolved from the nege-like virus ancestor through cross-species horizontal virus transmission. Our findings provide a new perspective on the virus-host coevolution and the origins of plant viruses.

Rich diversity of RNA viruses in the biological control agent, Orius laevigatus

Orius laevigatus (Hemiptera, Anthocoridae) is a generalist predator extensively used for the biocontrol of diverse agricultural pests. Previous studies on O. laevigatus have focused on the improvement of insect genetic traits, but little is known about its association with microbes, especially viruses that may influence its production and efficacy. More than 280 RNA viruses have been described in other Hemiptera insects, in line with the continuous discovery of insect-specific viruses (ISVs) boosted by next-generation sequencing. In this study, we characterized the repertoire of RNA viruses associated with O. laevigatus. Its virome comprises 27 RNA viruses, classified within fourteen viral families, of which twenty-three viruses are specific to O. laevigatus and four are likely associated with fungal microbiota. The analysis of viral abundance in five O. laevigatus populations confirmed the presence of simultaneous viral infections and highlighted the ubiquitous presence and high abundance of one solinvivirus and three totiviruses. Moreover, we identified 24 non-retroviral endogenous viral elements (nrEVEs) in the genome of O. laevigatus, suggesting a long-term relationship between the host and its virome. Although no symptoms were described in the insect populations under study, the high diversity of viral species and the high abundance of certain RNA viruses identified indicate that RNA viruses may be significant for the applicability and efficacy of O. laevigatus in biocontrol programs.

Apolipoprotein D3 and LOX product play a role in immune-priming of a lepidopteran insect, Spodoptera exigua

Immune-priming occurs in insects after a prior pathogen exposure. However, its underlying mechanism in insects remains elusive. In the present work, immune-priming was detected in a lepidopteran insect, Spodoptera exigua. Specifically, a prior infection with a heat-killed pathogenic bacterium, Escherichia coli, led to increased survival upon the second infection of different pathogens. Plasma collected from larvae with the prior infection possessed the immune-priming factor(s) that significantly up-regulated cellular and humoral immune responses of naïve larvae. Our study also finds that variations in the timing of plasma collection for priming larvae resulted in distinct impacts on both cellular and humoral responses. However, when the active plasma exhibiting the immune-priming was heat-treated, it lost this priming activity, therefore suggesting that protein factor(s) play a role in this immune-priming. An immunofluorescence assay showed that the hemocytes collected from the immune-primed larvae highly reacted to a polyclonal antibody specific to a vertebrate lipocalin, apolipoprotein D (ApoD). Among 27 ApoD genes (Se-ApoD1 ∼ Se-ApoD27) of S. exigua, Se-ApoD3 was found to be highly induced during the immune-priming, in which it was shown to be expressed in hemocytes and fat body from a fluorescence in situ hybridization analysis. RNA interference of Se-ApoD3 expression significantly impaired the immune-priming of S. exigua larvae. Moreover, the inhibition of eicosanoid biosynthesis suppressed the immune-priming, in which treatment with a lipoxygenase (LOX) inhibitor-and not treatment with a cyclooxygenase inhibitor-suppressed immune-priming. Further, an addition of LOX product such as lipoxin A4 or lipoxin B4 significantly rescued the lost immune-priming activity. Taken together, these results suggest that a complex of ApoD3 and LOX product mediates the immune-priming activity of S. exigua.

Innate and germline immune memory: specificity and heritability of the ancient immune mechanisms for adaptation and survival

The immune memory is one of the defensive strategies developed by both unicellular and multicellular organisms for ensuring their integrity and functionality. While the immune memory of the vertebrate adaptive immune system (based on somatic recombination) is antigen-specific, encompassing the generation of memory T and B cells that only recognize/react to a specific antigen epitope, the capacity of vertebrate innate cells to remember past events is a mostly non-specific mechanism of adaptation. This "innate memory" can be considered as germline-encoded because its effector tools (such as innate receptors) do not need somatic recombination for being active. Also, in several organisms the memory-related information is integrated in the genome of germline cells and can be transmitted to the progeny for several generations, but it can also be erased depending on the environmental conditions. Overall, depending on the organism, its environment and its living habits, innate immune memory appears to be a mechanism for achieving better protection and survival against repeated exposure to microbes/stressful agents present in the same environment or occurring in the same anatomical district, able to adapt to changes in the environmental cues. The anatomical and functional complexity of the organism and its lifespan drive the generation of different immune memory mechanisms, for optimal adaptation to changes in the living/environmental conditions. The concept of innate immunity being non-specific needs to be revisited, as a wealth of evidence suggests a significant degree of specificity both in the primary immune reaction and in the ensuing memory-like responses. This is clearly evident in invertebrate metazoans, in which distinct scenarios can be observed, with both non-specific (immune enhancement) or specific (immune priming) memory-like responses. In the case of mammals, there is evidence that some degree of specificity can be attained in different situations, for instance as organ-specific protection rather than microorganism-specific reaction. Thus, depending on the challenges and conditions, innate memory can be non-specific or specific, can be integrated in the germline and transmitted to the progeny or be short-lived, thereby representing an exceptionally plastic mechanism of defensive adaptation for ensuring individual and species survival.

Transposon and Transgene Tribulations in Mosquitoes: A Perspective of piRNA Proportions

Mosquitoes, like Drosophila, are dipterans, the order of "true flies" characterized by a single set of two wings. Drosophila are prime model organisms for biomedical research, while mosquito researchers struggle to establish robust molecular biology in these that are arguably the most dangerous vectors of human pathogens. Both insects utilize the RNA interference (RNAi) pathway to generate small RNAs to silence transposons and viruses, yet details are emerging that several RNAi features are unique to each insect family, such as how culicine mosquitoes have evolved extreme genomic feature differences connected to their unique RNAi features. A major technical difference in the molecular genetic studies of these insects is that generating stable transgenic animals are routine in Drosophila but still variable in stability in mosquitoes, despite genomic DNA-editing advances. By comparing and contrasting the differences in the RNAi pathways of Drosophila and mosquitoes, in this review we propose a hypothesis that transgene DNAs are possibly more intensely targeted by mosquito RNAi pathways and chromatin regulatory pathways than in Drosophila. We review the latest findings on mosquito RNAi pathways, which are still much less well understood than in Drosophila, and we speculate that deeper study into how mosquitoes modulate transposons and viruses with Piwi-interacting RNAs (piRNAs) will yield clues to improving transgene DNA expression stability in transgenic mosquitoes.

In vitro characterization of cell-fusing agent virus DNA forms in Aedes aegypti mosquitoes

How non-retroviral endogenous viral elements (EVEs) are established is a long-standing question. Viral DNA (vDNA) forms of RNA viruses are likely to be EVE precursors. Cell-fusing agent virus (CFAV) is a major insect-specific virus (ISV) in the Aedes aegypti mosquitoes and one of the few existing non-retroviral RNA viruses found as EVEs. We characterized CFAV-derived vDNA in the cell line to understand the mechanism of why current viruses are rarely endogenized. vDNA production was affected by cell culture media independent of CFAV replication. vDNAs that correspond to different regions covering the entire viral genome were detected, implying multiple initiation sites exist. A considerable proportion of vDNAs corresponded to ssDNA. Higher vDNA copies were detected in the cytoplasm than the nucleus. Our findings provide valuable insights into the intracellular characteristics of ISV-derived vDNAs, which will aid in understanding the underlying mechanisms of non-retroviral EVE formation.

RNA interference in insects: the link between antiviral defense and pest control

RNA interference (RNAi) is a form of gene silencing triggered by double-stranded RNA (dsRNA) that operates in all eukaryotic cells. RNAi has been widely investigated in insects to determine the underlying molecular mechanism, to investigate its role in systemic antiviral defense, and to develop strategies for pest control. When insect cells are infected by viruses, viral dsRNA signatures trigger a local RNAi response to block viral replication and generate virus-derived DNA that confers systemic immunity. RNAi-based insect pest control involves the application of exogenous dsRNA targeting genes essential for insect development or survival, but the efficacy of this approach has limited potency in many pests through a combination of rapid dsRNA degradation, inefficient dsRNA uptake/processing, and ineffective RNAi machinery. This could be addressed by dsRNA screening and evaluation, focusing on dsRNA design and off-target management, as well as dsRNA production and delivery. This review summarizes recent progress to determine the role of RNAi in antiviral defense and as a pest control strategy in insects, addressing gaps between our fundamental understanding of the RNAi mechanism and the exploitation of RNAi-based pest control strategies.

Mosquito transgenerational antiviral immunity is mediated by vertical transfer of virus DNA sequences and RNAi

Mosquitoes are important vectors for transmission of many viruses of public and veterinary health concern. These viruses most commonly have an RNA genome and infect mosquitoes for life. The principal mosquito antiviral response is the RNAi system which destroys virus RNA. Here, we confirm an earlier study that Aedes aegypti mosquitoes infected with positive-stranded RNA arboviruses can transmit specific immunity to their offspring. We show that this trans-generational immunity requires replication of virus RNA and reverse transcription of vRNA to vDNA in the infected parents and intergenerational transfer of vDNA. This vDNA is both genome-integrated and episomal. The episomal vDNA sequences are flanked by retrotransposon long-terminal repeats, predominantly Copia-like. Integrated vDNA sequences are propagated along several generations but specific immunity is effective only for a few generations and correlates with the presence of vRNA and episomal vDNA. This understanding raises new possibilities for the control of important mosquito-borne virus diseases.

Evolutionary genomics of three agricultural pest moths reveals rapid evolution of host adaptation and immune-related genes

BACKGROUND

Understanding the genotype of pest species provides an important baseline for designing integrated pest management (IPM) strategies. Recently developed long-read sequence technologies make it possible to compare genomic features of nonmodel pest species to disclose the evolutionary path underlying the pest species profiles. Here we sequenced and assembled genomes for 3 agricultural pest gelechiid moths: Phthorimaea absoluta (tomato leafminer), Keiferia lycopersicella (tomato pinworm), and Scrobipalpa atriplicella (goosefoot groundling moth). We also compared genomes of tomato leafminer and tomato pinworm with published genomes of Phthorimaea operculella and Pectinophora gossypiella to investigate the gene family evolution related to the pest species profiles.

RESULTS

We found that the 3 solanaceous feeding species, P. absoluta, K. lycopersicella, and P. operculella, are clustered together. Gene family evolution analyses with the 4 species show clear gene family expansions on host plant-associated genes for the 3 solanaceous feeding species. These genes are involved in host compound sensing (e.g., gustatory receptors), detoxification (e.g., ABC transporter C family, cytochrome P450, glucose-methanol-choline oxidoreductase, insect cuticle proteins, and UDP-glucuronosyl), and digestion (e.g., serine proteases and peptidase family S1). A gene ontology enrichment analysis of rapid evolving genes also suggests enriched functions in host sensing and immunity.

CONCLUSIONS

Our results of family evolution analyses indicate that host plant adaptation and pathogen defense could be important drivers in species diversification among gelechiid moths.

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.

Intrinsic factors driving mosquito vector competence and viral evolution: a review

Mosquitoes are responsible for the transmission of numerous viruses of global health significance. The term "vector competence" describes the intrinsic ability of an arthropod vector to transmit an infectious agent. Prior to transmission, the mosquito itself presents a complex and hostile environment through which a virus must transit to ensure propagation and transmission to the next host. Viruses imbibed in an infectious blood meal must pass in and out of the mosquito midgut, traffic through the body cavity or hemocoel, invade the salivary glands, and be expelled with the saliva when the vector takes a subsequent blood meal. Viruses encounter physical, cellular, microbial, and immunological barriers, which are influenced by the genetic background of the mosquito vector as well as environmental conditions. Collectively, these factors place significant selective pressure on the virus that impact its evolution and transmission. Here, we provide an overview of the current state of the field in understanding the mosquito-specific factors that underpin vector competence and how each of these mechanisms may influence virus evolution.

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.

The piRNA pathway is required for nucleopolyhedrovirus replication in Lepidoptera

The Piwi-interacting RNA (piRNA) pathway has been shown to be involved in the antiviral defense against RNA viruses, especially in mosquitoes, but its universality has been questioned. Here, we used the Bombyx mori nucleopolyhedrovirus (BmNPV) -infected silkworm as a model to explore the effects of the key factors of piRNA pathway, BmAgo3 and Siwi, on replication of a large DNA virus (belonging to the family of Baculoviridae). We demonstrated that BmAgo3 and Siwi could promote the replication of BmNPV through both overexpression and knockdown experiments in BmN cell lines and silkworm larvae. In addition, we also studied the effect of PIWI-class genes on Autographa californica nucleopolyhedrovirus (AcMNPV) replication in the Spodoptera frugiperda cell line Sf9. By knocking down the expression of PIWI-class genes in Sf9, we found that Piwi-like-1 and Piwi-like-2-3 could inhibit AcMNPV replication, while Piwi-like-4-5 promoted virus replication. Our study provides compelling evidence that the piRNA pathway affects host infection by exogenous viruses in Lepidoptera. Also, our results reflect the diversity of the roles of PIWI-class genes in virus infection of the host across species. This study is the first to explore the interaction of PIWI-class proteins with DNA viruses, providing new insights into the functional roles of the piRNA pathway.

The Aedes aegypti RNA interference response against Zika virus in the context of co-infection with dengue and chikungunya viruses

Since its detection in 2015 in Brazil, Zika virus (ZIKV) has remained in the spotlight of international public health and research as an emerging arboviral pathogen. In addition to single infection, ZIKV may occur in co-infection with dengue (DENV) and chikungunya (CHIKV) viruses, with whom ZIKV shares geographic distribution and the mosquito Aedes aegypti as a vector. The main mosquito immune response against arboviruses is RNA interference (RNAi). It is unknown whether or not the dynamics of the RNAi response differ between single arboviral infections and co-infections. In this study, we investigated the interaction of ZIKV and DENV, as well as ZIKV and CHIKV co-infections with the RNAi response in Ae. aegypti. Using small RNA sequencing, we found that the efficiency of small RNA production against ZIKV -a hallmark of antiviral RNAi-was mostly similar when comparing single and co-infections with either DENV or CHIKV. Silencing of key antiviral RNAi proteins, showed no change in effect on ZIKV replication when the cell is co-infected with ZIKV and DENV or CHIKV. Interestingly, we observed a negative effect on ZIKV replication during CHIKV co-infection in the context of Ago2-knockout cells, though his effect was absent during DENV co-infection. Overall, this study provides evidence that ZIKV single or co-infections with CHIKV or DENV are equally controlled by RNAi responses. Thus, Ae. aegypti mosquitoes and derived cells support co-infections of ZIKV with either CHIKV or DENV to a similar level than single infections, as long as the RNAi response is functional.

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.

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.

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.

RNA interference is essential to modulating the pathogenesis of mosquito-borne viruses in the yellow fever mosquito Aedes aegypti

While it has long been known that the transmission of mosquito-borne viruses depends on the establishment of persistent and nonlethal infections in the invertebrate host, specific roles for the insects' antiviral immune pathways in modulating the pathogenesis of viral infections is the subject of speculation and debate. Here, we show that a loss-of-function mutation in the Aedes aegypti Dicer-2 (Dcr-2) gene renders the insect acutely susceptible to a disease phenotype upon infection with pathogens in multiple virus families associated with important human diseases. Additional interrogation of the disease phenotype demonstrated that the virus-induced pathology is controlled through a canonical RNA interference (RNAi) pathway, which functions as a resistance mechanism. These results suggest comparatively modest contributions of proposed tolerance mechanisms to the fitness of A. aegypti infected with these pathogens. Similarly, the production of virus-derived piwi-interacting RNAs (vpiRNAs) was not sufficient to prevent the pathology associated with viral infections in Dcr-2 null mutants, also suggesting a less critical, or potentially secondary, role for vpiRNAs in antiviral immunity. These findings have important implications for understanding the ecological and evolutionary interactions occurring between A. aegypti and the pathogens they transmit to human and animal hosts.

Positive-strand RNA viruses-a Keystone Symposia report

Positive-strand RNA viruses have been the cause of several recent outbreaks and epidemics, including the Zika virus epidemic in 2015, the SARS outbreak in 2003, and the ongoing SARS-CoV-2 pandemic. On June 18-22, 2022, researchers focusing on positive-strand RNA viruses met for the Keystone Symposium "Positive-Strand RNA Viruses" to share the latest research in molecular and cell biology, virology, immunology, vaccinology, and antiviral drug development. This report presents concise summaries of the scientific discussions at the symposium.

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.

Antiviral RNAi Response in Culex quinquefasciatus-Derived HSU Cells

Culex spp. mosquitoes are important vectors of viruses, such as West Nile virus, Eastern equine encephalitis virus and Rift valley fever virus. However, their interactions with innate antiviral immunity, especially RNA interference (RNAi), are not well known. Most research on RNAi pathways in mosquitoes is focused on the tropical vector mosquito Aedes aegypti. Here, we investigated the production of arbovirus-specific small RNAs in Cx. quinquefasciatus-derived HSU cells. Furthermore, by silencing RNAi-related proteins, we investigated the antiviral role of these proteins for two different arboviruses: Semliki Forest virus (SFV) and Bunyamwera orthobunyavirus (BUNV). Our results showed an expansion of Ago2 and Piwi6 in Cx. quinquefasciatus compared to Ae. aegypti. While silencing Ago2a and Ago2b increased BUNV replication, only Ago2b showed antiviral activity against SFV. Our results suggest differences in the function of Cx. quinquefasciatus and Ae. aegypti RNAi proteins and highlight the virus-specific function of these proteins in Cx. quinquefasciatus.

Origins and diversification of animal innate immune responses against viral infections

Immune systems are of pivotal importance to any living organism on Earth, as they protect the organism against deleterious effects of viral infections. Though the current knowledge about these systems is still biased towards the immune response in vertebrates, some studies have focused on the identification and characterization of components of invertebrate antiviral immune systems. Two classic model organisms, the insect Drosophila melanogaster and the nematode Caenorhabditis elegans, were instrumental in the discovery of several important components of the innate immune system, such as the Toll-like receptors and the RNA interference pathway. However, these two model organisms provide only a limited view of the evolutionary history of the immune system, as they both are ecdysozoan protostomes. Recent functional studies in non-classic models such as unicellular holozoans (for example, choanoflagellates), lophotrochozoans (for example, oysters) and cnidarians (for example, sea anemones) have added crucial information for understanding the evolution of antiviral systems, as they revealed unexpected ancestral complexity. This Review aims to summarize this information and present the ancestral nature of the antiviral immune response in animals. We also discuss lineage-specific adaptations and future perspectives for the comparative study of the innate immune system that are essential for understanding its evolution.

Culex Mosquito Piwi4 Is Antiviral against Two Negative-Sense RNA Viruses

Culex spp. mosquitoes transmit several pathogens concerning public health, including West Nile virus and Saint Louis encephalitis virus. Understanding the antiviral immune system of Culex spp. mosquitoes is important for reducing the transmission of these viruses. Mosquitoes rely on RNA interference (RNAi) to control viral replication. While the siRNA pathway in mosquitoes is heavily studied, less is known about the piRNA pathway. The piRNA pathway in mosquitoes has recently been connected to mosquito antiviral immunity. In Aedes aegypti, Piwi4 has been implicated in antiviral responses. The antiviral role of the piRNA pathway in Culex spp. mosquitoes is understudied compared to Ae. aegypti. Here, we aimed to identify the role of PIWI genes and piRNAs in Culex quinquefasciatus and Culex tarsalis cells during virus infection. We examined the effect of PIWI gene silencing on virus replication of two arboviruses and three insect-specific viruses in Cx. quinquefasciatus derived cells (Hsu) and Cx. tarsalis derived (CT) cells. We show that Piwi4 is antiviral against the La Crosse orthobunyavirus (LACV) in Hsu and CT cells, and the insect-specific rhabdovirus Merida virus (MERDV) in Hsu cells. None of the silenced PIWI genes impacted replication of the two flaviviruses Usutu virus (USUV) and Calbertado virus, or the phasivirus Phasi-Charoen-like virus. We further used small RNA sequencing to determine that LACV-derived piRNAs, but not USUV-derived piRNAs were generated in Hsu cells and that PIWI gene silencing resulted in a small reduction in vpiRNAs. Finally, we determined that LACV-derived DNA was produced in Hsu cells during infection, but whether this viral DNA is required for vpiRNA production remains unclear. Overall, we expanded our knowledge on the piRNA pathway and how it relates to the antiviral response in Culex spp mosquitoes.

Virome and nrEVEome diversity of Aedes albopictus mosquitoes from La Reunion Island and China

Background

Aedes albopictus is a public health threat for its worldwide spread and ability to transmit arboviruses. Understanding mechanisms of mosquito immunity can provide new tools to control arbovirus spread. The genomes of Aedes mosquitoes contain hundreds of nonretroviral endogenous viral elements (nrEVEs), which are enriched in piRNA clusters and produce piRNAs, with the potential to target cognate viruses. Recently, one nrEVE was shown to limit cognate viral infection through nrEVE-derived piRNAs. These findings suggest that nrEVEs constitute an archive of past viral infection and that the landscape of viral integrations may be variable across populations depending on their viral exposure.

Methods

We used bioinformatics and molecular approaches to identify known and novel (i.e. absent in the reference genome) viral integrations in the genome of wild collected Aedes albopictus mosquitoes and characterize their virome.

Results

We showed that the landscape of viral integrations is dynamic with seven novel viral integrations being characterized, but does not correlate with the virome, which includes both viral species known and unknown to infect mosquitoes. However, the small RNA coverage profile of nrEVEs and the viral genomic contigs we identified confirmed an interaction among these elements and the piRNA and siRNA pathways in mosquitoes.

Conclusions

Mosquitoes nrEVEs have been recently described as a new form of heritable, sequence-specific mechanism of antiviral immunity. Our results contribute to understanding the dynamic distribution of nrEVEs in the genomes of wild Ae. albopictus and their interaction with mosquito viruses.

Understanding pathogen survival and transmission by arthropod vectors to prevent human disease

Many endemic poverty-associated diseases, such as malaria and leishmaniasis, are transmitted by arthropod vectors. Pathogens must interact with specific molecules in the vector gut, the microbiota, and the vector immune system to survive and be transmitted. The vertebrate host, in turn, is infected when the pathogen and vector-derived factors, such as salivary proteins, are delivered into the skin by a vector bite. Here, we review recent progress in our understanding of the biology of pathogen transmission from the human to the vector and back, from the vector to the host. We also highlight recent advances in the biology of vector-borne disease transmission, which have translated into additional strategies to prevent human disease by either reducing vector populations or by disrupting their ability to transmit pathogens.

Optimized In Vitro CRISPR/Cas9 Gene Editing Tool in the West Nile Virus Mosquito Vector, Culex quinquefasciatus

Culex quinquefasciatus mosquitoes are a globally widespread vector of multiple human and animal pathogens, including West Nile virus, Saint Louis encephalitis virus, and lymphatic filariasis. Since the introduction of West Nile virus to the United States in 1999, a cumulative 52,532 cases have been reported to the CDC, including 25,849 (49.2%) neuroinvasive cases and 2456 (5%) deaths. Viral infections elicit immune responses in their mosquito vectors, including the RNA interference (RNAi) pathway considered to be the cornerstone antiviral response in insects. To investigate mosquito host genes involved in pathogen interactions, CRISPR/Cas9-mediated gene-editing can be used for functional studies of mosquito-derived cell lines. Yet, the tools available for the study of Cx. quinquefasciatus-derived (Hsu) cell lines remain largely underdeveloped compared to other mosquito species. In this study, we constructed and characterized a Culex-optimized CRISPR/Cas9 plasmid for use in Hsu cell cultures. By comparing it to the original Drosophila melanogaster CRISPR/Cas9 plasmid, we showed that the Culex-optimized plasmid demonstrated highly efficient editing of the genomic loci of the RNAi proteins Dicer-2 and PIWI4 in Hsu cells. These new tools support our ability to investigate gene targets involved in mosquito antiviral response, and thus the future development of gene-based vector control strategies.

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.

Shrimp genome sequence contains independent clusters of ancient and current Endogenous Viral Elements (EVE) of the parvovirus IHHNV

Background

Shrimp have the ability to accommodate viruses in long term, persistent infections without signs of disease. Endogenous viral elements (EVE) play a role in this process probably via production of negative-sense Piwi-interacting RNA (piRNA)-like fragments. These bind with Piwi proteins to dampen viral replication via the RNA interference (RNAi) pathway. We searched a genome sequence (GenBank record JABERT000000000) of the giant tiger shrimp (Penaeus monodon for the presence of EVE related to a shrimp parvovirus originally named infectious hypodermal and hematopoietic necrosis virus (IHHNV).

Results

The shrimp genome sequence contained three piRNA-like gene clusters containing scrambled IHHNV EVE. Two clusters were located distant from one another in pseudochromosome 35 (PC35). Both PC35 clusters contained multiple sequences with high homology (99%) to GenBank records DQ228358 and EU675312 that were both called “non-infectious IHHNV Type A” (IHHNV-A) when originally discovered. However, our results and those from a recent Australian P. monodon genome assembly indicate that the relevant GenBank records for IHHNV-A are sequence-assembly artifacts derived from scrambled and fragmental IHHNV-EVE. Although the EVE in the two PC35 clusters showed high homology only to IHHNV-A, the clusters were separate and distinct with respect to the arrangement (i.e., order and reading direction) and proportional content of the IHHNV-A GenBank records. We conjecture that these 2 clusters may constitute independent allele-like clusters on a pair of homologous chromosomes. The third EVE cluster was found in pseudochromosome 7 (PC7). It contained EVE with high homology (99%) only to GenBank record AF218266 with the potential to protect shrimp against current types of infectious IHHNV. One disadvantage was that some EVE in PC7 can give false positive PCR test results for infectious IHHNV.

Conclusions

Our results suggested the possibility of viral-type specificity in EVE clusters. Specificity is important because whole EVE clusters for one viral type would be transmitted to offspring as collective hereditary units. This would be advantageous if one or more of the EVE within the cluster were protective against the disease caused by the cognate virus. It would also facilitate gene editing for removal of non-protective EVE clusters or for transfer of protective EVE clusters to genetically improve existing shrimp breeding stocks that might lack them.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08802-3.

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.

Invading viral DNA triggers dsRNA synthesis by RNA polymerase II to activate antiviral RNA interference in Drosophila

dsRNA sensing triggers antiviral responses against RNA and DNA viruses in diverse eukaryotes. In Drosophila, Invertebrate iridescent virus 6 (IIV-6), a large DNA virus, triggers production of small interfering RNAs (siRNAs) by the dsRNA sensor Dicer-2. Here, we show that host RNA polymerase II (RNAPII) bidirectionally transcribes specific AT-rich regions of the IIV-6 DNA genome to generate dsRNA. Both replicative and naked IIV-6 genomes trigger production of dsRNA in Drosophila cells, implying direct sensing of invading DNA. Loquacious-PD, a Dicer-2 co-factor essential for the biogenesis of endogenous siRNAs, is dispensable for processing of IIV-6-derived dsRNAs, which suggests that they are distinct. Consistent with this finding, inhibition of the RNAPII co-factor P-TEFb affects the synthesis of endogenous, but not virus-derived, dsRNA. Altogether, our results suggest that a non-canonical RNAPII complex recognizes invading viral DNA to synthesize virus-derived dsRNA, which activates the antiviral siRNA pathway in Drosophila.

BmCPV-Derived Circular DNA vcDNA-S7 Mediated by Bombyx mori Reverse Transcriptase (RT) Regulates BmCPV Infection

Circular DNAs derived from single-stranded RNA viruses play important roles in counteracting viral infection. However, whether double-stranded RNA viruses generate functional circular DNAs is still unknown. Using circDNA sequencing, divergent PCR, DNA in situ hybridization and rolling circular amplification, we presently confirmed that in silkworm, Bombyx mori cytoplasmic polyhedrosis virus (BmCPV), a double-stranded RNA virus belonging to cypovirus, is prone to produce a BmCPV-derived circular DNA termed as vcDNA-S7. We have also found that vcDNA-S7 formation is mediated by endogenous reverse transcriptase (RT), and the proliferation of BmCPV can be inhibited by vcDNA-S7 in vitro and in vivo. Moreover, we have discovered that the silkworm RNAi immune pathway is activated by vcDNA-S7, while viral small interfering RNAs (vsiRNAs) derived from transcribed RNA by vcDNA-S7 can be detected by small RNA deep sequencing. These results suggest that BmCPV-derived vcDNA-S7, mediated by RT, can serve as a template for the biogenesis of antiviral siRNAs, which may lead to the repression of BmCPV infection. To our knowledge, this is the first demonstration that a circular DNA, produced by double stranded RNA viruses, is capable of regulating virus infection.

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.

Evolution and Diversity of Inherited Viruses in the Nearctic Phantom Midge, Chaoborus americanus

Inherited mutualists, parasites, and commensals occupy one of the most intimate ecological niches available to invertebrate-associated microbes. How this transmission environment influences microbial evolution is increasingly understood for inherited bacterial symbionts, but in viruses, research on the prevalence of vertical transmission and its effects on viral lineages is still maturing. The evolutionary stability of this strategy remains difficult to assess, although phylogenetic evidence of frequent host shifts and selective sweeps have been interpreted as strategies favoring parasite persistence. In this study, we describe and investigate a natural insect system in which species-wide sweeps have been restricted by the isolation of host populations. Previous work identified evidence of pronounced mitochondrial genetic structure among North American populations of the phantom midge, Chaoborus americanus. Here we take advantage of the geographical isolation in this species to investigate the diversity and persistence of its inherited virome. We identify eight novel RNA viruses from six families and use small RNA sequencing in reproductive tissues to provide evidence of vertical transmission. We report region-specific virus strains that mirror the continental phylogeography of the host, demonstrating that members of the inherited virome have independently persisted in parallel host lineages since they last shared a common ancestor in the Mid Pleistocene. We find that the small interfering RNA pathway, a frontline of antiviral defense in insects, targets members of this inherited virome. Finally, our results suggest that the Piwi-mediated RNA silencing pathway is unlikely to function as a general antiviral defense in Chaoborus, in contrast to its role in some mosquitoes. However, we also report that the PIWI-interacting RNA pathway generates abundant piRNAs from endogenous viral elements closely related to actively infecting inherited viruses, potentially helping to explain idiosyncratic patterns of virus-specific Piwi targeting in this insect.

The Neglected Virome of Triatomine Insects

The Triatominae subfamily (Reduviidae) harbors some hematophagous insect species that have been firmly connected to the transmission of Trypanosoma cruzi, the causative agent of Chagas disease. Triatomines not only host and transmit trypanosomatids, but also coexist with a variety of symbiotic microorganisms that generally reside in the insect’s intestinal flora. The microbiome has profound effects on the physiology, immunity, fitness and survival of animals and plants. The interaction between triatomines and bacteria has been investigated to some extent and has revealed important bacteria symbionts. In contrast, the range of viral species that can infect triatomine insects is almost completely unknown. In some cases, genomic and metatranscriptomic approaches have uncovered sequences related to possible viral genomes, but, to date, only eight positive single-strand RNA viruses, namely Triatoma virus and Rhodnius prolixus viruses 1 - 7 have been investigated in more detail. Here, we review the literature available on triatomine viruses and the viruses-insect host relationship. The lack of broader metagenomic and metatranscriptomic studies in these medically relevant insects underscores the importance of expanding our knowledge of the triatomine virome both for surveillance purposes as well as to possibly harness their potential for insect vector population control strategies.

BmNPV p35 Reduces the Accumulation of Virus-Derived siRNAs and Hinders the Function of siRNAs to Facilitate Viral Infection

Antiviral immunity involves various mechanisms and responses, including the RNA interference (RNAi) pathway. During long-term coevolution, viruses have gained the ability to evade this defense by encoding viral suppressors of RNAi (VSRs). It was reported that p35 of baculovirus can inhibit cellular small interference RNA (siRNA) pathway; however, the molecular mechanisms underlying p35 as a VSR remain largely unclear. Here, we showed that p35 of Bombyx mori nucleopolyhedrovirus (BmNPV) reduces the accumulation of virus-derived siRNAs (vsiRNAs) mapped to a particular region in the viral genome, leading to an increased expression of the essential genes in this region, and revealed that p35 disrupts the function of siRNAs by preventing them from loading into Argonaute-2 (Ago2). This repressive effect on the cellular siRNA pathway enhances the replication of BmNPV. Thus, our findings illustrate for the first time the inhibitory mechanism of a baculovirus VSR and how this effect influences viral infection.

Antiviral RNAi Response against the Insect-Specific Agua Salud Alphavirus

Arboviruses transmitted by mosquitoes are responsible for the death of millions of people each year. In addition to arboviruses, many insect-specific viruses (ISVs) have been discovered in mosquitoes in the last decade. ISVs, in contrast to arboviruses transmitted by mosquitoes to vertebrates, cannot replicate in vertebrate cells even when they are evolutionarily closely related to arboviruses. The alphavirus genus includes many arboviruses, although only a few ISVs have been discovered from this genus so far. Here, we investigate the interactions of a recently isolated insect-specific alphavirus, Agua Salud alphavirus (ASALV), with its mosquito host. RNA interference (RNAi) is one of the essential antiviral responses against arboviruses, although there is little knowledge on the interactions of RNAi with ISVs. Through the knockdown of transcripts of the different key RNAi pathway (small interfering RNA [siRNA], microRNA [miRNA], and P-element-induced wimpy testis [PIWI]-interacting RNA [piRNA]) proteins, we show the antiviral role of Ago2 (siRNA), Ago1 (miRNA), and Piwi4 proteins against ASALV in Aedes aegypti-derived cells. ASALV replication was increased in Dicer2 and Ago2 knockout cells, confirming the antiviral role of the siRNA pathway. In infected cells, mainly ASALV-specific siRNAs are produced, while piRNA-like small RNAs, with the characteristic nucleotide bias resulting from ping-pong amplification, are produced only in Dicer2 knockout cells. Taken together, ASALV interactions with the mosquito RNAi response differ from those of arthropod-borne alphaviruses in some aspects, although they also share some commonalities. Further research is needed to understand whether the identified differences can be generalized to other insect-specific alphaviruses. IMPORTANCE Mosquitoes are efficient vectors for many arboviruses that cause emergent infectious diseases in humans. Many insect-specific viruses (ISVs) that can infect mosquitoes but cannot infect vertebrates have been discovered in the last decade. ISVs have attracted great attention due to their potential use in mosquito or arbovirus control, by either decreasing mosquito fitness or restricting arbovirus replication and transmission to humans. However, ISV-mosquito interactions are not well understood. RNA interference (RNAi) is the most important innate immune response against many arboviruses, while it is unknown if it is antiviral against ISVs. Here, we investigate in detail the antiviral effect of the RNAi response in mosquitoes against an ISV for the first time. Using a recently isolated insect-specific alphavirus, we show that the regulation of virus replication was different from that for arthropod-borne alphaviruses despite some similarities. The differences in mosquito-virus interactions could drive the different transmission modes, which could eventually drive the evolution of arboviruses. Hence, an understanding of mosquito-ISV interactions can shed light on the ecology and evolution of both ISVs and the medically important 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.

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.

Somatic piRNAs and Transposons are Differentially Expressed Coincident with Skeletal Muscle Atrophy and Programmed Cell Death

PIWI-interacting RNAs (piRNAs) are small single-stranded RNAs that can repress transposon expression via epigenetic silencing and transcript degradation. They have been identified predominantly in the ovary and testis, where they serve essential roles in transposon silencing in order to protect the integrity of the genome in the germline. The potential expression of piRNAs in somatic cells has been controversial. In the present study we demonstrate the expression of piRNAs derived from both genic and transposon RNAs in the intersegmental muscles (ISMs) from the tobacco hawkmoth Manduca sexta. These piRNAs are abundantly expressed, ∼27 nt long, map antisense to transposons, are oxidation resistant, exhibit a 5’ uridine bias, and amplify via the canonical ping-pong pathway. An RNA-seq analysis demonstrated that 19 piRNA pathway genes are expressed in the ISMs and are developmentally regulated. The abundance of piRNAs does not change when the muscles initiate developmentally-regulated atrophy, but are repressed coincident with the commitment of the muscles undergo programmed cell death at the end of metamorphosis. This change in piRNA expression is correlated with the repression of several retrotransposons and the induction of specific DNA transposons. The developmentally-regulated changes in the expression of piRNAs, piRNA pathway genes, and transposons are all regulated by 20-hydroxyecdysone, the steroid hormone that controls the timing of ISM death. Taken together, these data provide compelling evidence for the existence of piRNA in somatic tissues and suggest that they may play roles in developmental processes such as programmed cell death.

Extending and Running the Mosquito Small RNA Genomics Resource Pipeline

The Mosquito Small RNA Genomics (MSRG) resource is a repository of analyses on the small RNA transcriptomes of mosquito cell cultures and somatic and gonadal tissues. This resource allows for comparing the regulation dynamics of small RNAs generated from transposons and viruses across mosquito species. This chapter covers the procedures to set up the MSRG resource pipeline as a new installation by detailing the necessary collection of genome reference and annotation files and lists of microRNAs (miRNAs) hairpin sequences, transposon repeats consensus sequences, and virus genome sequences. Proper execution of the MSRG resource pipeline yields outputs amenable to biologists to further analyze with desktop and spreadsheet software to gain insights into the balance between arthropod endogenous small RNA populations and the proportions of virus-derived small RNAs that include Piwi-interacting RNAs (piRNAs) and endogenous small interfering RNAs (siRNAs).