Chromosomal scale assembly of parasitic wasp genome reveals symbiotic virus colonization

Chromosome-level genome assembly of the bethylid ectoparasitoid wasp Sclerodermus sp. 'alternatusi'

The Bethylidae are the most diverse of Hymenoptera chrysidoid families. As external parasitoids, the bethylids have been widely adopted as biocontrol agents to control insect pests worldwide. Thus far, the genomic information of the family Bethylidae has not been reported yet. In this study, we crystallized into a high-quality chromosome-level genome of ant-like bethylid wasps Sclerodermus sp. 'alternatusi' (Hymenoptera: Bethylidae) using PacBio sequencing as well as Hi-C technology. The assembled S. alternatusi genome was 162.30 Mb in size with a contig N50 size of 3.83 Mb and scaffold N50 size of 11.10 Mb. Totally, 92.85% assembled sequences anchored to 15 pseudo-chromosomes. A total of 10,204 protein-coding genes were annotated, and 23.01 Mb repetitive sequences occupying 14.17% of genome were pinpointed. The BUSCO results showed that 97.9% of the complete core Insecta genes were identified in the genome, while 97.1% in the gene sets. The high-quality genome of S. alternatusi will not only provide valuable genomic information, but also show insights into parasitoid wasp evolution and bio-control application in future studies.

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.

A chromosome scale assembly of the parasitoid wasp Venturia canescens provides insight into the process of virus domestication

The parasitoid wasp Venturia canescens is an important biological control agent of stored products moth pests and serves as a model to study the function and evolution of domesticated endogenous viruses (DEVs). The DEVs discovered in V. canescens are known as virus-like particles (VcVLPs), which are produced using nudivirus-derived components and incorporate wasp-derived virulence proteins instead of packaged nucleic acids. Previous studies of virus-derived components in the V. canescens genome identified 53 nudivirus-like genes organized in six gene clusters and several viral pseudogenes, but how VcVLP genes are organized among wasp chromosomes following their integration in the ancestral wasp genome is largely unknown. Here, we present a chromosomal scale genome of V. canescens consisting of 11 chromosomes and 56 unplaced small scaffolds. The genome size is 290.8 Mbp with a N50 scaffold size of 24.99 Mbp. A high-quality gene set including 11,831 protein-coding genes were produced using RNA-Seq data as well as publicly available peptide sequences from related Hymenoptera. A manual annotation of genes of viral origin produced 61 intact and 19 pseudogenized nudivirus-derived genes. The genome assembly revealed that two previously identified clusters were joined into a single cluster and a total of 5 gene clusters comprising of 60 intact nudivirus-derived genes were located in three chromosomes. In contrast, pseudogenes are dispersed among 8 chromosomes with only 4 pseudogenes associated with nudivirus gene clusters. The architecture of genes encoding VcVLP components suggests it originates from a recent virus acquisition and there is a link between the processes of dispersal and pseudogenization. This high-quality genome assembly and annotation represents the first chromosome-scale assembly for parasitoid wasps associated with VLPs, and is publicly available in the National Center for Biotechnology Information Genome and RefSeq databases, providing a valuable resource for future studies of DEVs in parasitoid wasps.

Chromosome-level genome assembly of Microplitis manilae Ashmead, 1904 (Hymenoptera: Braconidae)

Microplitis manilae Ashmead (Hymenoptera: Braconidae) is an important parasitoid of agricultural pests in lepidopteran species. So far, two extant genome assembles from the genus Microplitis are fragmented. Here, we offered a high-quality genome assembly of M. manilae at the chromosome level with high accuracy and contiguity, assembled by ONT long-read, MGI-SEQ short-read, and Hi-C sequencing methods. The final assembled genome size was 282.85 Mb, with 268.17 Mb assigned to 11 pseudochromosomes. The scaffold N50 length was 25.23 Mb, and the complete BUSCO score was 98.61%. The genome contained 152.37 Mb of repetitive elements, representing 53.87% of the total genome size. We predicted 15,689 protein-coding genes, of which 13,580 genes were annotated functionally. Gene family evolution investigations of M. manilae revealed 615 expanded and 635 contracted gene families. The high-quality genome of M. manilae reported in this paper will be a useful genomic resource for research on parasitoid wasps in the future.

Genome assembly of the ectoparasitoid wasp Theocolax elegans

The ectoparasitoid wasp Theocolax elegans is a cosmopolitan and generalist pteromalid parasitoid of several major storage insect pests, and can effectively suppress a host population in warehouses. However, little molecular information about this wasp is currently available. In this study, we assembled the genome of T. elegans using PacBio long-read sequencing, Illumina sequencing, and Hi-C methods. The genome assembly is 662.73 Mb in length with contig and scaffold N50 values of 1.15 Mb and 88.8 Mb, respectively. The genome contains 56.4% repeat sequences and 23,212 protein-coding genes were annotated. Phylogenomic analyses revealed that T. elegans diverged from the lineage leading to subfamily Pteromalinae (Nasonia vitripennis and Pteromalus puparum) approximately 110.5 million years ago. We identified 130 significantly expanded gene families, 34 contracted families, 248 fast-evolving genes, and 365 positively selected genes in T. elegans. Additionally, 260 olfactory receptors and 285 venom proteins were identified. This genome assembly provides valuable genetic bases for future investigations on evolution, molecular biology and application of T. elegans.

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.

[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.

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

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

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.

Viral histones: pickpocket's prize or primordial progenitor?

The common histones H2A, H2B, H3, and H4 are the characteristic components of eukaryotic nucleosomes, which function to wrap DNA and compact the genome as well as to regulate access to DNA for transcription and replication in all eukaryotes. In the past two decades, histones have also been found to be encoded in some DNA viruses, where their functions and properties are largely unknown, though recently histones from two related viruses have been shown to form nucleosome-like structures in vitro. Viral histones can be highly similar to eukaryotic histones in primary sequence, suggesting they have been recently picked up from eukaryotic hosts, or they can be radically divergent in primary sequence and may occur as conjoined histone doublets, triplets, or quadruplets, suggesting ancient origins prior to the divergence of modern eukaryotes. Here, we review what is known of viral histones and discuss their possible origins and functions. We consider how the viral life cycle may affect their properties and histories, and reflect on the possible roles of viruses in the origin of the nucleus of modern eukaryotic cells.

Genome of the parasitoid wasp Cotesia chilonis sheds light on amino acid resource exploitation

Background

A fundamental feature of parasitism is the nutritional exploitation of host organisms by their parasites. Parasitoid wasps lay eggs on arthropod hosts, exploiting them for nutrition to support larval development by using diverse effectors aimed at regulating host metabolism. However, the genetic components and molecular mechanisms at the basis of such exploitation, especially the utilization of host amino acid resources, remain largely unknown. To address this question, here, we present a chromosome-level genome assembly of the parasitoid wasp Cotesia chilonis and reconstruct its amino acid biosynthetic pathway.

Results

Analyses of the amino acid synthetic pathway indicate that C. chilonis lost the ability to synthesize ten amino acids, which was confirmed by feeding experiments with amino acid-depleted media. Of the ten pathways, nine are known to have been lost in the common ancestor of animals. We find that the ability to synthesize arginine was also lost in C. chilonis because of the absence of two key genes in the arginine synthesis pathway. Further analyses of the genomes of 72 arthropods species show that the loss of arginine synthesis is common in arthropods. Metabolomic analyses by UPLC-MS/MS reveal that the temporal concentrations of arginine, serine, tyrosine, and alanine are significantly higher in host (Chilo suppressalis) hemolymph at 3 days after parasitism, whereas the temporal levels of 5-hydroxylysine, glutamic acid, methionine, and lysine are significantly lower. We sequence the transcriptomes of a parasitized host and non-parasitized control. Differential gene expression analyses using these transcriptomes indicate that parasitoid wasps inhibit amino acid utilization and activate protein degradation in the host, likely resulting in the increase of amino acid content in host hemolymph.

Conclusions

We sequenced the genome of a parasitoid wasp, C. chilonis, and revealed the features of trait loss in amino acid biosynthesis. Our work provides new insights into amino acid exploitation by parasitoid wasps, and this knowledge can specifically be used to design parasitoid artificial diets that potentially benefit mass rearing of parasitoids for pest control.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01313-3.

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.

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.

OUP accepted manuscript

Chelonus formosanus Sonan 1932 (Hymenoptera: Braconidae) is a wasp capable of parasitizing a variety of lepidopteran pests at the “egg-larval” stage which distributes throughout Taiwan, Guangdong, Zhejiang, and Hainan provinces of China. This wasp has been successfully used to control pests such as Spodoptera litura Fabricius, 1775, Spodoptera frugiperda (JE Smith, 1797), Spodoptera exigua (Hübner, 1808), and Helicoverpa armigera (Hübner, 1808). So far, there is only one genome assembled from the Chelonus genus [Chelonus insularis (Cresson, 1865)] and it is fragmented with 455 scaffolds. Here, we report a chromosome-level genome assembly of C. formosanus, which was sequenced using PacBio, Illumina, and Hi-C technologies. The long reads were 35.4 Gb (∼150× coverage) with an average length of 15.23 kb. The size of the genome assembly was 139.59 Mb. More than 99.46% of the assembled sequences were anchored to seven pseudochromosomes (138.84 Mb). The Benchmarking University Single-Copy Orthologs (BUSCO) assessment results showed 99.0% of the 1,367 genes (insect_odb10 database) were completely present. We annotated 11,242 protein-coding genes including 98.6% of BUSCO complete genes that were recovered. Nearly one-fourth of the genome assembly (22.25%) was annotated as repetitive sequences and 324 noncoding RNAs were predicted. There were 58 gene families found with significant expansion including allelopathic families (odorant receptors and ionotropic receptors), which may play a crucial role in efficiently locating a wide range of hosts. This high-quality genome assembly and annotation could provide a highly valuable resource of parasitic wasp for the biological control of Lepidoptera pest.

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.

Paleovirology of the DNA viruses of eukaryotes

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

Quantitative trait locus analysis of parasitoid counteradaptation to symbiont-conferred resistance

Insect hosts and parasitoids are engaged in an intense struggle of antagonistic coevolution. Infection with heritable bacterial endosymbionts can substantially increase the resistance of aphids to parasitoid wasps, which exerts selection on parasitoids to overcome this symbiont-conferred protection (counteradaptation). Experimental evolution in the laboratory has produced counteradapted populations of the parasitoid wasp Lysiphlebus fabarum. These populations can parasitize black bean aphids (Aphis fabae) protected by the bacterial endosymbiont Hamiltonella defensa, which confers high resistance against L. fabarum. We used two experimentally evolved parasitoid populations to study the genetic architecture of the counteradaptation to symbiont-conferred resistance by QTL analysis. With simple crossing experiments, we showed that the counteradaptation is a recessive trait depending on the maternal genotype. Based on these results, we designed a customized crossing scheme to genotype a mapping population phenotyped for the ability to parasitize Hamiltonella-protected aphids. Using 1835 SNP markers obtained by ddRAD sequencing, we constructed a high-density linkage map consisting of six linkage groups (LGs) with an overall length of 828.3 cM and an average marker spacing of 0.45 cM. We identified a single QTL associated with the counteradaptation to Hamiltonella in L. fabarum on linkage group 2. Out of 120 genes located in this QTL, several genes encoding putative venoms may represent candidates for counteradaptation, as parasitoid wasps inject venoms into their hosts during oviposition.

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.

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.

A chromosome-level genome assembly of the parasitoid wasp, Cotesia glomerata (Hymenoptera: Braconidae)

Hymenopterans make up about 20% of all animal species, but most are poorly known and lack high-quality genomic resources. One group of important, yet under-studied hymenopterans, are parasitoid wasps in the family Braconidae. Among this under-studied group are braconid wasps in the genus Cotesia; a clade of ~1,000 species routinely used in studies of physiology, ecology, biological control, and genetics. However, our ability to understand these organisms has been hindered by a lack of genomic resources. We helped bridge this gap by generating a high-quality genome assembly for the parasitoid wasp, Cotesia glomerata (Braconidae; Microgastrinae). We generated this assembly using multiple sequencing technologies, including Oxford Nanopore, whole-genome shotgun sequencing, and 3-D chromatin contact information (Hi-C). Our assembly is one of the most contiguous, complete, and publicly available hymenopteran genomes, represented by 3,355 scaffolds with a scaffold N50 of ~28Mb and a BUSCO score of ~99%. Given the genome sizes found in closely related species, our genome assembly was ~50% larger than expected, which was apparently induced by runaway amplification of three types of repetitive elements: simple repeats, Long Terminal Repeats (LTRs), and Long Interspersed Nuclear Elements (LINEs). This assembly is another step forward for genomics across this hyper-diverse, yet understudied, order of insects. The assembled genomic data and metadata files are publicly available via Figshare (https://doi.org/10.6084/m9.figshare.13010549).