Genomic architecture of endogenous ichnoviruses reveals distinct evolutionary pathways leading to virus domestication in parasitic wasps

Identification of a viral gene essential for the genome replication of a domesticated endogenous virus in ichneumonid parasitoid wasps

Thousands of endoparasitoid wasp species in the families Braconidae and Ichneumonidae harbor "domesticated endogenous viruses" (DEVs) in their genomes. This study focuses on ichneumonid DEVs, named ichnoviruses (IVs). Large quantities of DNA-containing IV virions are produced in ovary calyx cells during the pupal and adult stages of female wasps. Females parasitize host insects by injecting eggs and virions into the body cavity. After injection, virions rapidly infect host cells which is followed by expression of IV genes that promote the successful development of wasp offspring. IV genomes consist of two components: proviral segment loci that serve as templates for circular dsDNAs that are packaged into capsids, and genes from an ancestral virus that produce virions. In this study, we generated a chromosome-scale genome assembly for Hyposoter didymator that harbors H. didymator ichnovirus (HdIV). We identified a total of 67 HdIV loci that are amplified in calyx cells during the wasp pupal stage. We then focused on an HdIV gene, U16, which is transcribed in calyx cells during the initial stages of replication. Sequence analysis indicated that U16 contains a conserved domain in primases from select other viruses. Knockdown of U16 by RNA interference inhibited virion morphogenesis in calyx cells. Genome-wide analysis indicated U16 knockdown also inhibited amplification of HdIV loci in calyx cells. Altogether, our results identified several previously unknown HdIV loci, demonstrated that all HdIV loci are amplified in calyx cells during the pupal stage, and showed that U16 is required for amplification and virion morphogenesis.

The Viromes of Six Ecosystem Service Provider Parasitoid Wasps

Parasitoid wasps are fundamental insects for the biological control of agricultural pests. Despite the importance of wasps as natural enemies for more sustainable and healthy agriculture, the factors that could impact their species richness, abundance, and fitness, such as viral diseases, remain almost unexplored. Parasitoid wasps have been studied with regard to the endogenization of viral elements and the transmission of endogenous viral proteins that facilitate parasitism. However, circulating viruses are poorly characterized. Here, RNA viromes of six parasitoid wasp species are studied using public libraries of next-generation sequencing through an integrative bioinformatics pipeline. Our analyses led to the identification of 18 viruses classified into 10 families (Iflaviridae, Endornaviridae, Mitoviridae, Partitiviridae, Virgaviridae, Rhabdoviridae, Chuviridae, Orthomyxoviridae, Xinmoviridae, and Narnaviridae) and into the Bunyavirales order. Of these, 16 elements were described for the first time. We also found a known virus previously identified on a wasp prey which suggests viral transmission between the insects. Altogether, our results highlight the importance of virus surveillance in wasps as its service disruption can affect ecology, agriculture and pest management, impacting the economy and threatening human food security.

Developments in the classification and nomenclature of arthropod-infecting large DNA viruses that contain pif genes

Viruses of four families of arthropod-specific, large dsDNA viruses (the nuclear arthropod large DNA viruses, or NALDVs) possess homologs of genes encoding conserved components involved in the baculovirus primary infection mechanism. The presence of such homologs encoding per os infectivity factors (pif genes), along with their absence from other viruses and the occurrence of other shared characteristics, suggests a common origin for the viruses of these families. Therefore, the class Naldaviricetes was recently established, accommodating these four families. In addition, within this class, the ICTV approved the creation of the order Lefavirales for three of these families, whose members carry homologs of the baculovirus genes that code for components of the viral RNA polymerase, which is responsible for late gene expression. We further established a system for the binomial naming of all virus species in the order Lefavirales, in accordance with a decision by the ICTV in 2019 to move towards a standardized nomenclature for all virus species. The binomial species names for members of the order Lefavirales consist of the name of the genus to which the species belongs (e.g., Alphabaculovirus), followed by a single epithet that refers to the host species from which the virus was originally isolated. The common names of viruses and the abbreviations thereof will not change, as the format of virus names lies outside the remit of the ICTV.

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.

Massive Somatic and Germline Chromosomal Integrations of Polydnaviruses in Lepidopterans

Increasing numbers of horizontal transfer (HT) of genes and transposable elements are reported in insects. Yet the mechanisms underlying these transfers remain unknown. Here we first quantify and characterize the patterns of chromosomal integration of the polydnavirus (PDV) encoded by the Campopleginae Hyposoter didymator parasitoid wasp (HdIV) in somatic cells of parasitized fall armyworm (Spodoptera frugiperda). PDVs are domesticated viruses injected by wasps together with their eggs into their hosts in order to facilitate the development of wasp larvae. We found that six HdIV DNA circles integrate into the genome of host somatic cells. Each host haploid genome suffers between 23 and 40 integration events (IEs) on average 72 h post-parasitism. Almost all IEs are mediated by DNA double-strand breaks occurring in the host integration motif (HIM) of HdIV circles. We show that despite their independent evolutionary origins, PDV from both Campopleginae and Braconidae wasps use remarkably similar mechanisms for chromosomal integration. Next, our similarity search performed on 775 genomes reveals that PDVs of both Campopleginae and Braconidae wasps have recurrently colonized the germline of dozens of lepidopteran species through the same mechanisms they use to integrate into somatic host chromosomes during parasitism. We found evidence of HIM-mediated HT of PDV DNA circles in no less than 124 species belonging to 15 lepidopteran families. Thus, this mechanism underlies a major route of HT of genetic material from wasps to lepidopterans with likely important consequences on lepidopterans.

Genome-wide identification of lncRNAs associated with viral infection in Spodoptera frugiperda

The role of lncRNAs in immune defence has been demonstrated in many multicellular and unicellular organisms. However, investigation of the identification and characterization of long non-coding RNAs (lncRNAs) involved in the insect immune response is still limited. In this study, we used RNA sequencing (RNA-seq) to investigate the expression profiles of lncRNAs and mRNAs in the fall armyworm Spodoptera frugiperda in response to virus infection. To assess the tissue- and virus-specificity of lncRNAs, we analysed and compared their expression profiles in haemocytes and fat body of larvae infected with two entomopathogenic viruses with different lifestyles, i.e. the polydnavirus HdIV (Hyposoter didymator IchnoVirus) and the densovirus JcDV (Junonia coenia densovirus). We identified 1883 candidate lncRNAs, of which 529 showed differential expression following viral infection. Expression profiles differed considerably between samples, indicating that many differentially expressed (DE) lncRNAs showed virus- and tissue-specific expression patterns. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and target prediction analyses indicated that DE-LncRNAs were mainly enriched in metabolic process, DNA replication and repair, immune response, metabolism of insect hormone and cell adhesion. In addition, we identified three DE-lncRNAs potentially acting as microRNA host genes, suggesting that they participate in gene regulation by producing miRNAs in response to virus infection. This study provides a catalogue of lncRNAs expressed in two important immune tissues and potential insight into their roles in the antiviral defence in S. frugiperda. The results may help future in-depth functional studies to better understand the biological function of lncRNAs in interaction between viruses and the fall armyworm.

Genome, host genome integration, and gene expression in Diadegma fenestrale ichnovirus from the perspective of coevolutionary hosts

Polydnaviruses (PDVs) exhibit species-specific mutualistic relationships with endoparasitoid wasps. PDVs can be categorized into bracoviruses and ichnoviruses, which have independent evolutionary origins. In our previous study, we identified an ichnovirus of the endoparasitoid Diadegma fenestrale and named it DfIV. Here, DfIV virions from the ovarian calyx of gravid female wasps were characterized. DfIV virion particles were ellipsoidal (246.5 nm × 109.0 nm) with a double-layered envelope. Next-generation sequencing of the DfIV genome revealed 62 non-overlapping circular DNA segments (A1-A5, B1-B9, C1-C15, D1-D23, E1-E7, and F1-F3); the aggregate genome size was approximately 240 kb, and the GC content (43%) was similar to that of other IVs (41%-43%). A total of 123 open reading frames were predicted and included typical IV gene families such as repeat element protein (41 members), cysteine motif (10 members), vankyrin (9 members), polar residue-rich protein (7 members), vinnexin (6 members), and N gene (3 members). Neuromodulin N (2 members) was found to be unique to DfIV, along with 45 hypothetical genes. Among the 62 segments, 54 showed high (76%-98%) sequence similarities to the genome of Diadegma semiclausum ichnovirus (DsIV). Three segments, namely, D22, E3, and F2, contained lepidopteran host genome integration motifs with homologous regions of about 36-46 bp between them (Diadegma fenestrale ichnovirus, DfIV and lepidopteran host, Plutella xylostella). Most of the DfIV genes were expressed in the hymenopteran host and some in the lepidopteran host (P. xylostella), parasitized by D. fenestrale. Five segments (A4, C3, C15, D5, and E4) were differentially expressed at different developmental stages of the parasitized P. xylostella, and two segments (C15 and D14) were highly expressed in the ovaries of D. fenestrale. Comparative analysis between DfIV and DsIV revealed that the genomes differed in the number of segments, composition of sequences, and internal sequence homologies.

[Plants and animals biological functions obtained from viruses]

Viruses can provide new biological functions to plants and animals. Some viruses persisting at low levels in plants might confer resistance to stress and parasites. In animals, more numerous examples of genes originating from viruses and used by different organisms have been described. For examples these genes might contribute to protect from new infections, or to ensure communication between neurons or to enable placenta development. In parasitic wasps, a complex viral machinery has been conserved as an endogenous virus dispersed in the wasp genome, which produces virions. These virions infect the parasitized host resulting in the production of virulence factors that inhibit defense mechanisms against the parasite. Different organisms have used the same viral functions repeatedly during animal evolution.

Chelonus inanitus bracovirus encodes lineage-specific proteins and truncated immune IκB-like factors

Bracoviruses and ichnoviruses are endogenous viruses of parasitic wasps that produce particles containing virulence genes expressed in host tissues and necessary for parasitism success. In the case of bracoviruses the particles are produced by conserved genes of nudiviral origin integrated permanently in the wasp genome, whereas the virulence genes can strikingly differ depending on the wasp lineage. To date most data obtained on bracoviruses concerned species from the braconid subfamily of Microgastrinae. To gain a broader view on the diversity of virulence genes we sequenced the genome packaged in the particles of Chelonus inanitus bracovirus (CiBV) produced by a wasp belonging to a different subfamily: the Cheloninae. These are egg-larval parasitoids, which means that they oviposit into the host egg and the wasp larvae then develop within the larval stages of the host. We found that most of CiBV virulence genes belong to families that are specific to Cheloninae. As other bracoviruses and ichnoviruses however, CiBV encode v-ank genes encoding truncated versions of the immune cactus/IκB factor, which suggests these proteins might play a key role in host–parasite interactions involving domesticated endogenous viruses. We found that the structures of CiBV V-ANKs are different from those previously reported. Phylogenetic analysis supports the hypothesis that they may originate from a cactus/IκB immune gene from the wasp genome acquired by the bracovirus. However, their evolutionary history is different from that shared by other V-ANKs, whose common origin probably reflects horizontal gene transfer events of virus sequences between braconid and ichneumonid wasps.

Many evolutionary roads led to virus domestication in ichneumonoid parasitoid wasps

The investigation of endogenous viral elements (EVEs) has historically focused on only a few lineages of parasitoid wasps, with negative results consistently underreported. Recent studies show that multiple viral lineages were integrated in at least seven instances in Ichneumonoidea and may be much more widespread than previously thought. Increasingly affordable genomic and bioinformatic approaches have made it feasible to search for viral sequences within wasp genomes, opening an extremely promising research avenue. Advances in wasp phylogenetics have shed light on the evolutionary history of EVE integration, although many questions remain. Phylogenetic proximity can be used as a guide to facilitate targeted screening, to estimate the number and age of integration events and to identify taxa involved in major host switches.

Identifying bracovirus and ichnovirus genes involved in virion morphogenesis

Bracoviruses (BVs) and ichnoviruses (IVs) evolved from different endogenized viruses but through convergence have been coopted by parasitoids in the families Braconidae and Ichneumonidae for similar functions in parasitizing hosts. Experimentally studying the role of endogenized viral genes in virion morphogenesis remains a key challenge in the study of BVs and IVs. Here we summarize how multiomics, electron microscopy, and RNA interference (RNAi) methods have provided new insights about BV and IV gene function.

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.

Bracoviruses, ichnoviruses, and virus-like particles from parasitoid wasps retain many features of their virus ancestors

Animal genomes commonly contain genes or sequences that have been acquired from different types of viruses. The vast majority of these endogenous virus elements (EVEs) are inactive or consist of only a small number of components that show no evidence of cooption for new functions or interaction. Unlike most EVEs, bracoviruses (BVs), ichnoviruses (IVs) and virus-like particles (VLPs) in parasitoid wasps have evolved through retention and interaction of many genes from virus ancestors. Here, we discuss current understanding of BV, IV and VLP evolution along with associated implications for what constitutes a virus. We suggest that BVs and IVs are domesticated endogenous viruses (DEVs) that differ in several important ways from other known EVEs.

Genome-Wide Patterns of Bracovirus Chromosomal Integration into Multiple Host Tissues during Parasitism

Bracoviruses are domesticated viruses found in parasitic wasp genomes. They are composed of genes of nudiviral origin that are involved in particle production and proviral segments containing virulence genes that are necessary for parasitism success. During particle production, proviral segments are amplified and individually packaged as DNA circles in nucleocapsids. These particles are injected by parasitic wasps into host larvae together with their eggs. Bracovirus circles of two wasp species were reported to undergo chromosomal integration in parasitized host hemocytes, through a conserved sequence named the host integration motif (HIM). Here, we used bulk Illumina sequencing to survey integrations of Cotesia typhae bracovirus circles in the DNA of its host, the maize corn borer (Sesamia nonagrioides), 7 days after parasitism. First, assembly and annotation of a high-quality genome for C. typhae enabled us to characterize 27 proviral segments clustered in proviral loci. Using these data, we characterized large numbers of chromosomal integrations (from 12 to 85 events per host haploid genome) for all 16 bracovirus circles containing a HIM. Integrations were found in four S. nonagrioides tissues and in the body of a caterpillar in which parasitism had failed. The 12 remaining circles do not integrate but are maintained at high levels in host tissues. Surprisingly, we found that HIM-mediated chromosomal integration in the wasp germ line has occurred accidentally at least six times during evolution. Overall, our study furthers our understanding of wasp-host genome interactions and supports HIM-mediated chromosomal integration as a possible mechanism of horizontal transfer from wasps to their hosts. IMPORTANCE Bracoviruses are endogenous domesticated viruses of parasitoid wasps that are injected together with wasp eggs into wasp host larvae during parasitism. Several studies have shown that some DNA circles packaged into bracovirus particles become integrated into host somatic genomes during parasitism, but the phenomenon has never been studied using nontargeted approaches. Here, we use bulk Illumina sequencing to systematically characterize and quantify bracovirus circle integrations that occur in four tissues of the Mediterranean corn borer (Sesamia nonagrioides) during parasitism by the Cotesia typhae wasp. Our analysis reveals that all circles containing a HIM integrate at substantial levels (from 12 to 85 integrations per host cell, in total) in all tissues, while other circles do not integrate. In addition to shedding new light on wasp-bracovirus-host interactions, our study supports HIM-mediated chromosomal integration of bracovirus as a possible source of wasp-to-host horizontal transfer, with long-term evolutionary consequences.

Proteo-Trancriptomic Analyses Reveal a Large Expansion of Metalloprotease-Like Proteins in Atypical Venom Vesicles of the Wasp Meteorus pulchricornis (Braconidae)

Meteorus pulchricornis (Ichneumonoidea, Braconidae) is an endoparasitoid wasp of lepidopteran caterpillars. Its parasitic success relies on vesicles (named M. pulchricornis Virus-Like Particles or MpVLPs) that are synthesized in the venom gland and injected into the parasitoid host along with the venom during oviposition. In order to define the content and understand the biogenesis of these atypical vesicles, we performed a transcriptome analysis of the venom gland and a proteomic analysis of the venom and purified MpVLPs. About half of the MpVLPs and soluble venom proteins identified were unknown and no similarity with any known viral sequence was found. However, MpVLPs contained a large number of proteins labelled as metalloproteinases while the most abundant protein family in the soluble venom was that of proteins containing the Domain of Unknown Function DUF-4803. The high number of these proteins identified suggests that a large expansion of these two protein families occurred in M. pulchricornis. Therefore, although the exact mechanism of MpVLPs formation remains to be elucidated, these vesicles appear to be “metalloproteinase bombs” that may have several physiological roles in the host including modifying the functions of its immune cells. The role of DUF4803 proteins, also present in the venom of other braconids, remains to be clarified.

Plant-phenotypic changes induced by parasitoid ichnoviruses enhance the performance of both unparasitized and parasitized caterpillars

There is increasing awareness that interactions between plants and insects can be mediated by microbial symbionts. Nonetheless, evidence showing that symbionts associated with organisms beyond the second trophic level affect plant‐insect interactions are restricted to a few cases belonging to parasitoid‐associated bracoviruses. Insect parasitoids harbour a wide array of symbionts which, like bracoviruses, can be injected into their herbivorous hosts to manipulate their physiology and behaviour. Yet, the function of these symbionts in plant‐based trophic webs remains largely overlooked. Here, we provide the first evidence of a parasitoid‐associated symbiont belonging to the group of ichnoviruses which affects the strength of plant‐insect interactions. A comparative proteomic analysis shows that, upon parasitoid injection of calyx fluid containing ichnovirus particles, the composition of salivary glands of caterpillars changes both qualitatively (presence of two viral‐encoded proteins) and quantitatively (abundance of several caterpillar‐resident enzymes, including elicitors such as glucose oxidase). In turn, plant phenotypic changes triggered by the altered composition of caterpillar oral secretions affect the performance of herbivores. Ichnovirus manipulation of plant responses to herbivory leads to benefits for their parasitoid partners in terms of reduced developmental time within the parasitized caterpillar. Interestingly, plant‐mediated ichnovirus‐induced effects also enhance the performances of unparasitized herbivores which in natural conditions may feed alongside parasitized ones. We discuss these findings in the context of ecological costs imposed to the plant by the viral symbiont of the parasitoid. Our results provide intriguing novel findings about the role played by carnivore‐associated symbionts on plant‐insect‐parasitoid systems and underline the importance of placing mutualistic associations in an ecological perspective.

Chromosomal scale assembly of parasitic wasp genome reveals symbiotic virus colonization

Endogenous viruses form an important proportion of eukaryote genomes and a source of novel functions. How large DNA viruses integrated into a genome evolve when they confer a benefit to their host, however, remains unknown. Bracoviruses are essential for the parasitism success of parasitoid wasps, into whose genomes they integrated ~103 million years ago. Here we show, from the assembly of a parasitoid wasp genome at a chromosomal scale, that bracovirus genes colonized all ten chromosomes of Cotesia congregata. Most form clusters of genes involved in particle production or parasitism success. Genomic comparison with another wasp, Microplitis demolitor, revealed that these clusters were already established ~53 mya and thus belong to remarkably stable genomic structures, the architectures of which are evolutionary constrained. Transcriptomic analyses highlight temporal synchronization of viral gene expression without resulting in immune gene induction, suggesting that no conflicts remain between ancient symbiotic partners when benefits to them converge.

The Unconventional Viruses of Ichneumonid Parasitoid Wasps

To ensure their own immature development as parasites, ichneumonid parasitoid wasps use endogenous viruses that they acquired through ancient events of viral genome integration. Thousands of species from the campoplegine and banchine wasp subfamilies rely, for their survival, on their association with these viruses, hijacked from a yet undetermined viral taxon. Here, we give an update of recent findings on the nature of the viral genes retained from the progenitor viruses and how they are organized in the wasp genome.

Genomic architecture of endogenous ichnoviruses reveals distinct evolutionary pathways leading to virus domestication in parasitic wasps

BACKGROUND

Polydnaviruses (PDVs) are mutualistic endogenous viruses inoculated by some lineages of parasitoid wasps into their hosts, where they facilitate successful wasp development. PDVs include the ichnoviruses and bracoviruses that originate from independent viral acquisitions in ichneumonid and braconid wasps respectively. PDV genomes are fully incorporated into the wasp genomes and consist of (1) genes involved in viral particle production, which derive from the viral ancestor and are not encapsidated, and (2) proviral segments harboring virulence genes, which are packaged into the viral particle. To help elucidating the mechanisms that have facilitated viral domestication in ichneumonid wasps, we analyzed the structure of the viral insertions by sequencing the whole genome of two ichnovirus-carrying wasp species, Hyposoter didymator and Campoletis sonorensis.

RESULTS

Assemblies with long scaffold sizes allowed us to unravel the organization of the endogenous ichnovirus and revealed considerable dispersion of the viral loci within the wasp genomes. Proviral segments contained species-specific sets of genes and occupied distinct genomic locations in the two ichneumonid wasps. In contrast, viral machinery genes were organized in clusters showing highly conserved gene content and order, with some loci located in collinear wasp genomic regions. This genomic architecture clearly differs from the organization of PDVs in braconid wasps, in which proviral segments are clustered and viral machinery elements are more dispersed.

CONCLUSIONS

The contrasting structures of the two types of ichnovirus genomic elements are consistent with their different functions: proviral segments are vehicles for virulence proteins expected to adapt according to different host defense systems, whereas the genes involved in virus particle production in the wasp are likely more stable and may reflect ancestral viral architecture. The distinct genomic architectures seen in ichnoviruses versus bracoviruses reveal different evolutionary trajectories that have led to virus domestication in the two wasp lineages.