Phylogenomics of Ichneumonoidea (Hymenoptera) and implications for evolution of mode of parasitism and viral endogenization

Parasitoid wasps (Ichneumonoidea) collected from faba bean fields, Kharga Oasis, New Valley, Egypt, with new records and the description of a new species

A faunistic inventory of ichneumonoid wasps collected from faba bean fields (Vicia faba L.) in Kharga Oasis, New Valley, Egypt is presented. Fourteen species of Braconidae in 11 genera and five subfamilies (Aphidiinae, Braconinae, Euphorinae, Microgastrinae, and Opiinae); and six ichneumonid species in five genera and three subfamilies (Campopleginae, Diplazontinae, and Ichneumoninae) were collected and identified. Among those, three genera Barichneumon Thomson, and Virgichneumon Heinrich (Ichneumoninae, Ichneumonidae), and Centistes Haliday (Euphorinae, Braconidae), and five species: Cotesia icipe Fernndez & Fiaboe, 2017 (Microgastrinae) (Braconidae), Diadegma trochanteratum (Thomson, 1887), and Sinophorus tibialis Sanborne, 1984 (Campopleginae), Barichneumon bilunulatus (Gravenhorst, 1829), and Virgichneumon callicerus (Gravenhorst, 1820) (Ichneumoninae, Ichneumonidae) are new records to the Egyptian fauna. A new species, Trioxys khargaiensis Gadallah & Edmardash (Aphidiinae: Braconidae) is described and illustrated. A faunistic list comprising diagnostic characters and distribution in the Middle East and North Africa is provided. Illustrations of species characters are also provided.

The PWWP domain and the evolution of unique DNA methylation toolkits in Hymenoptera

DNMT3 in Hymenoptera has a unique duplication of the essential PWWP domain. Using GST-tagged PWWP fusion proteins and histone arrays we show that these domains have gained new properties and represent the first case of PWWP domains binding to H3K27 chromatin modifications, including H3K27me3, a key modification that is important during development. Phylogenetic analyses of 107 genomes indicate that the duplicated PWWP domains separated into two sister clades, and their distinct binding capacities are supported by 3D modeling. Other features of this unique DNA methylation system include variable copies, losses, and duplications of DNMT1 and DNMT3, and combinatorial generations of DNMT3 isoforms including variants missing the catalytic domain. Some of these losses and duplications of are found only in parasitic wasps. We discuss our findings in the context of the crosstalk between DNA methylation and histone methylation, and the expanded potential of epigenomic modifications in Hymenoptera to drive evolutionary novelties.

A revised terminology for male genitalia in Hymenoptera (Insecta), with a special emphasis on Ichneumonoidea

Applying consistent terminology for morphological traits across different taxa is a highly pertinent task in the study of morphology and evolution. Different terminologies for the same traits can generate bias in phylogeny and prevent correct homology assessments. This situation is exacerbated in the male genitalia of Hymenoptera, and specifically in Ichneumonoidea, in which the terminology is not standardized and has not been fully aligned with the rest of Hymenoptera. In the current contribution, we review the terms used to describe the skeletal features of the male genitalia in Hymenoptera, and provide a list of authors associated with previously used terminology. We propose a unified terminology for the male genitalia that can be utilized across the order and a list of recommended terms. Further, we review and discuss the genital musculature for the superfamily Ichneumonoidea based on previous literature and novel observations and align the terms used for muscles across the literature.

Chromosome-Scale Genome Assembly of the Solitary Parasitoid Wasp Microplitis manilae Ashmead, 1904 (Braconidae: Microgastrinae)

Parasitoid wasps are invaluable natural enemies extensively used to control coleopteran, dipteran, and lepidopteran pests in agriculture and forestry owing to their killing and reproductive actions on hosts. The important larval endoparasitoid wasp Microplitis manilae, which belongs to the Microgastrinae subfamily, parasitizes the larval stages of Spodoptera spp., such as Spodoptera litura and Spodoptera frugiperda. The absence of a genomic resource for M. manilae has impeded studies on chemosensory- and detoxification-related genes. This study presents a chromosome-level genome assembly of M. manilae with a genome size of 293.18 Mb, which includes 222 contigs (N50 size, 7.58 Mb) and 134 scaffolds (N50 size, 27.33 Mb). A major proportion of the genome (284.76 Mb; 97.13%) was anchored to 11 pseudochromosomes with a single-copy BUSCO score of 98.4%. Furthermore, 14,316 protein-coding genes, 165.14 Mb (57.99%) repetitive elements, and 871 noncoding RNAs were annotated and identified. Additionally, a manual annotation of 399 genes associated with chemosensation and 168 genes involved in detoxification was conducted. This study provides a valuable and high-quality genomic resource to facilitate further functional genomics research on parasitoid wasps.

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.

Biodiversity of hymenopteran parasitoids

Parasitoid wasps are the most successful group of insect parasitoids, comprising more than half the known diversity of Hymenoptera and probably most of the unknown diversity. This lifestyle has enabled them to be used as pest control agents conferring substantial economic benefits to global agriculture. Major lineages of parasitoid wasps include Ichneumonoidea, Ceraphronoidea, Proctotrupomorpha, and a number of aculeate families. The parasitoid lifestyle arose only once among basal Hymenoptera, in the common ancestor of the Orussidae and Apocrita some 200+ Ma ago. The ancestral parasitoid wasp was probably an idiobiont on wood-living beetle larvae. From this comparatively simple biology, Hymenoptera radiated into an incredible diversity of hosts and parasitoid lifestyles, including hyperparasitoidism, kleptoparasitoidism, egg parasitoidism, and polyembryony, in several instances co-opting viruses to subdue their hosts. Many lineages evolved beyond the parasitoid niche, becoming secondarily herbivorous or predatory nest provisioners and eventually giving rise to most instances of insect societies.

Viriforms-A New Category of Classifiable Virus-Derived Genetic Elements

The International Committee on Taxonomy of Viruses (ICTV) recently accepted viriforms as a new polyphyletic category of classifiable virus-derived genetic elements, juxtaposed to the polyphyletic virus, viroid, and satellite nucleic acid categories. Viriforms are endogenized former viruses that have been exapted by their cellular hosts to fulfill functions important for the host's life cycle. While morphologically resembling virions, particles made by viriforms do not package the viriform genomes but instead transport host genetic material. Known viriforms are highly diverse: members of family Polydnaviriformidae (former Polydnaviridae) have thus far been found exclusively in the genomes of braconid and ichneumonid parasitoid wasps, whereas the completely unrelated gene transfer agents (GTAs) are widely distributed among prokaryotes. In addition, recent discoveries likely extend viriforms to mammalian genomes. Here, we briefly outline the properties of these viriform groups and the first accepted and proposed ICTV frameworks for viriform classification.

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.

Conserved Viral Transcription Plays a Key Role in Virus-Like Particle Production of the Parasitoid Wasp Venturia canescens

Nudiviruses are large double-stranded DNA viruses related to baculoviruses known to be endogenized in the genomes of certain parasitic wasp species. These wasp-virus associations allow the production of viral particles or virus-like particles that ensure wasp parasitism success within lepidopteran hosts. Venturia canescens is an ichneumonid wasp belonging to the Campopleginae subfamily that has endogenized nudivirus genes belonging to the Alphanudivirus genus to produce "virus-like particles" (Venturia canescens virus-like particles [VcVLPs]), which package proteic virulence factors. The main aim of this study was to determine whether alphanudivirus gene functions have been conserved following endogenization. The expression dynamics of alphanudivirus genes was monitored by a high throughput transcriptional approach, and the functional role of lef-4 and lef-8 genes predicted to encode viral RNA polymerase components was investigated by RNA interference. As described for baculovirus infections and for endogenized nudivirus genes in braconid wasp species producing bracoviruses, a transcriptional cascade involving early and late expressed alphanudivirus genes could be observed. The expression of lef-4 and lef-8 was also shown to be required for the expression of alphanudivirus late genes allowing correct particle formation. Together with previous literature, the results show that endogenization of nudiviruses in parasitoid wasps has repeatedly led to the conservation of the viral RNA polymerase function, allowing the production of viruses or viral-like particles that differ in composition but enable wasp parasitic success. IMPORTANCE This study shows that endogenization of a nudivirus genome in a Campopleginae parasitoid wasp has led to the conservation, as for endogenized nudiviruses in braconid parasitoid wasps, of the viral RNA polymerase function, required for the transcription of genes encoding viral particles involved in wasp parasitism success. We also showed for the first time that RNA interference (RNAi) can be successfully used to downregulate gene expression in this species, a model in behavioral ecology. This opens the opportunity to investigate the function of genes involved in other traits important for parasitism success, such as reproductive strategies and host choice. Fundamental data acquired on gene function in Venturia canescens are likely to be transferable to other parasitoid wasp species used in biological control programs. This study also renders possible the investigation of other nudivirus gene functions, for which little data are available.

Mitochondrial phylogenomics and mitogenome organization in the parasitoid wasp family Braconidae (Hymenoptera: Ichneumonoidea)

Background

Mitochondrial (mt) nucleotide sequence data has been by far the most common tool employed to investigate evolutionary relationships. While often considered to be more useful for shallow evolutionary scales, mt genomes have been increasingly shown also to contain valuable phylogenetic information about deep relationships. Further, mt genome organization provides another important source of phylogenetic information and gene reorganizations which are known to be relatively frequent within the insect order Hymenoptera. Here we used a dense taxon sampling comprising 148 mt genomes (132 newly generated) collectively representing members of most of the currently recognised subfamilies of the parasitoid wasp family Braconidae, which is one of the largest radiations of hymenopterans. We employed this data to investigate the evolutionary relationships within the family and to assess the phylogenetic informativeness of previously known and newly discovered mt gene rearrangements.

Results

Most subfamilial relationships and their composition obtained were similar to those recovered in a previous phylogenomic study, such as the restoration of Trachypetinae and the recognition of Apozyginae and Proteropinae as valid braconid subfamilies. We confirmed and detected phylogenetic signal in previously known as well as novel mt gene rearrangements, including mt rearrangements within the cyclostome subfamilies Doryctinae and Rogadinae.

Conclusions

Our results showed that both the mt genome DNA sequence data and gene organization contain valuable phylogenetic signal to elucidate the evolution within Braconidae at different taxonomic levels. This study serves as a basis for further investigation of mt gene rearrangements at different taxonomic scales within the family.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01983-1.

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.

A chromosome-scale genome and transcriptomic analysis of the endangered tropical tree Vatica mangachapoi (Dipterocarpaceae)

Vatica mangachapoi is a tropical tree species native to Southeast Asia. It has long been valued as a timber species because the wood resists decay, but it is now considered vulnerable to extinction due to habitat loss and overexploitation. Here, we present the first chromosome-level genome assembly of V. mangachapoi that we created by combining data from PacBio long read sequencing with Hi-C proximity ligation and Illumina short-read sequencing. The assembled genome was 456.21 Mb, containing 11 chromosome and a BUSCO score of 93.4%. From the newly assembled genome, 46,811 protein-coding genes were predicted. Repetitive DNA accounted for 53% of the genome. Phylogenomic and gene family analyses showed that V. mangachapoi diverged from a common ancestor of Gossypium raimondii 70 million years ago. Transcriptome analyses found 227 genes that were differentially expressed in the leaves of plants grown in normal soil relative to plants grown in dry, coastal, sandy soil. For these genes, we identified three significantly enriched with GO terms: responses to organonitrogen compounds, chitin-triggered immunity, and wound response. This genome provides an important comparative benchmark not only for future conservation work on V. mangachapoi but also for phylogenomics work on Dipterocarpaceae.

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.

Mitochondrial Genomes Yield Insights into the Basal Lineages of Ichneumonid Wasps (Hymenoptera: Ichneumonidae)

We obtained four mitochondrial genomes of Odontocolon sp., Xorides funiuensis, Euceros kiushuensis and Euceros serricornis, which represent the first two representatives from subfamily Eucerotinae and Xoridinae (Ichneumonidae), respectively. All of the 4 newly sequenced mitochondrial genomes contain 13 protein-coding genes (PCGs) and most 24 RNA genes. Furthermore, they all have novel tRNA rearrangement patterns comparing with published mitogenomes in Ichneumonidae. For the tRNA cluster trnI-trnQ-trnM, X. funiuensis is shuffled as trnM-trnI-trnQ with trnQ inversed, while Odontocolon sp. with a remote translocation of trnK, shuffling as trnI-trnM-trnQ. E. kiushuensis and E. serricornis shared the same cluster trnQ-trnY-trnW-trnC. Finally, we reconstructed the phylogenetic relationships among the sequenced subfamilies of Ichneumonidae based on nucleotides and amino acids sequences of 13 PCGs in mitochondrial genomes, and the results of both the maximum likelihood and Bayesian inference analyses highly support that Eucerotinae is the basal ichneumonid lineage rather than Xoridinae.

Comparative Mitochondrial Genomics of 104 Darwin Wasps (Hymenoptera: Ichneumonidae) and Its Implication for Phylogeny

Simple Summary

Nearly a hundred mitochondrial genomes of ichneumonid wasps are newly reported. Comparative mitogenomics of 104 mitochondrial genomes representing 33 subfamilies of Ichneumonidae, as well as its implications for phylogeny, were studied. We found that the mitochondrial genomes of ichneumonid wasps were highly conserved in their base composition and had low evolutionary rates, but were diverse in gene order. There are 38 types of gene rearrangement events in 104 ichneumonid mitochondrial genomes, of which 30 novel rearrangement types (R3-6, R8-R10, R12-R15, R17-R18, R20-R35 and R38) and a hot spot rearrangement around R1, with a shuffled tRNA cluster trnW-trnY-trnC and trnI-trnQ-trnM, were detected. The relationships among these subfamilies are firstly discussed based on mitochondrial genomes at a large scale. We suggest five subfamily groupings of Ichneumonidae: Brachycyrtiformes, Ichneumoniformes, Ophioniformes, Pimpliformes and Xoridiformes. Two formerly unplaced subfamilies, Eucerotinae and Microleptinae, were placed in Brachycyrtiformes and Ichneumoniformes, respectively.

Abstract

Ichneumonidae is one of the largest families of insects with a mega-diversity of specialized morphological and biological characteristics. We newly sequenced 92 mitochondrial genomes of ichneumonid wasps and found that they have a conserved base composition and a lower evolutionary rate than that of other families of parasitic Hymenoptera. There are 38 types of gene order in the ichneumonid mitochondrial genome, with 30 novel types identified in 104 ichneumonids. We also found that the rearrangement events occur more frequently in Ophioniformes than in Ichneumoniformes and Pimpliformes. Furthermore, the higher Ophioniformes and their relative lineages shared the transposition of trnL2 to trnI-trnQ-trnM tRNA cluster. We confirmed five higher-level groupings of Ichneumonidae: Brachycyrtiformes, Ichneumoniformes, Ophioniformes, Pimpliformes and Xoridiformes. Two formerly unplaced subfamilies, Eucerotinae and Microleptinae, were placed in Brachycyrtiformes and Ichneumoniformes, respectively. The results will improve our understanding of the diversity and evolution of Ichneumonidae.

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.

Including fossils in phylogeny: a glimpse into the evolution of the superfamily Evanioidea (Hymenoptera: Apocrita) under tip-dating and the fossilized birth–death process

Using a fossilized birth–death model, a new phylogeny of the superfamily Evanioidea (including ensign wasps, nightshade wasps and hatchet wasps) is proposed, with estimates of divergence times for its constitutive families and for corroborating the monophyly of Evanioidea. Additionally, our Bayesian analyses demonstrate the monophyly of †Anomopterellidae, †Othniodellithidae, †Andreneliidae, Aulacidae, Gasteruptiida and Evaniidae, whereas †Praeaulacidae and †Baissidae appear to be paraphyletic. Vectevania vetula and Hyptiogastrites electrinus are transferred to Aulacidae. We estimate the divergence time of Evanioidea to be in the Late Triassic (~203 Mya). Additionally, three new othniodellithid wasps are described and figured from mid-Cretaceous Burmese amber as the new genus Keratodellitha, with three new species: Keratodellitha anubis sp. nov., Keratodellitha basilisci sp. nov. and Keratodellitha kirin sp. nov. We also document a temporal shift in relative species richness between Ichneumonoidea and Evanioidea.

The Presence of Ancient Core Genes Reveals Endogenization from Diverse Viral Ancestors in Parasitoid Wasps

The Ichneumonoidea (Ichneumonidae and Braconidae) is an incredibly diverse superfamily of parasitoid wasps that includes species that produce virus-like entities in their reproductive tracts to promote successful parasitism of host insects. Research on these entities has traditionally focused upon two viral genera Bracovirus (in Braconidae) and Ichnovirus (in Ichneumonidae). These viruses are produced using genes known collectively as endogenous viral elements (EVEs) that represent historical, now heritable viral integration events in wasp genomes. Here, new genome sequence assemblies for 11 species and 6 publicly available genomes from the Ichneumonoidea were screened with the goal of identifying novel EVEs and characterizing the breadth of species in lineages with known EVEs. Exhaustive similarity searches combined with the identification of ancient core genes revealed sequences from both known and novel EVEs. One species harbored a novel, independently derived EVE related to a divergent large double-stranded DNA (dsDNA) virus that manipulates behavior in other hymenopteran species. Although bracovirus or ichnovirus EVEs were identified as expected in three species, the absence of ichnoviruses in several species suggests that they are independently derived and present in two younger, less widespread lineages than previously thought. Overall, this study presents a novel bioinformatic approach for EVE discovery in genomes and shows that three divergent virus families (nudiviruses, the ancestors of ichnoviruses, and Leptopilina boulardi Filamentous Virus-like viruses) are recurrently acquired as EVEs in parasitoid wasps. Virus acquisition in the parasitoid wasps is a common process that has occurred in many more than two lineages from a diverse range of arthropod-infecting dsDNA viruses.