Functional annotation of Cotesia congregata bracovirus: identification of viral genes expressed in parasitized host immune tissues

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.

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.

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.

The Complete Genome of Chelonus insularis Reveals Dynamic Arrangement of Genome Components in Parasitoid Wasps That Produce Bracoviruses

Bracoviruses (BVs) are endogenized nudiviruses in parasitoid wasps of the microgastroid complex (family Braconidae). Microgastroid wasps have coopted nudivirus genes to produce replication-defective virions that females use to transfer virulence genes to parasitized hosts. The microgastroid complex further consists of six subfamilies and ∼50,000 species but current understanding of BV gene inventories and organization primarily derives from analysis of two wasp species in the subfamily Microgastrinae (Microplitis demolitor and Cotesia congregata) that produce M. demolitor BV (MdBV) and C. congregata BV (CcBV). Notably, several genomic features of MdBV and CcBV remain conserved since divergence of M. demolitor and C. congregata ∼53 million years ago (MYA). However, it is unknown whether these conserved traits more broadly reflect BV evolution, because no complete genomes exist for any microgastroid wasps outside the Microgastrinae. In this regard, the subfamily Cheloninae is of greatest interest because it diverged earliest from the Microgastrinae (∼85 MYA) after endogenization of the nudivirus ancestor. Here, we present the complete genome of Chelonus insularis, which is an egg-larval parasitoid in the Cheloninae that produces C. insularis BV (CinsBV). We report that the inventory of nudivirus genes in C. insularis is conserved but are dissimilarly organized compared to M. demolitor and C. congregata. Reciprocally, CinsBV proviral segments share organizational features with MdBV and CcBV but virulence gene inventories exhibit almost no overlap. Altogether, our results point to the functional importance of a conserved inventory of nudivirus genes and a dynamic set of virulence genes for the successful parasitism of hosts. Our results also suggest organizational features previously identified in MdBV and CcBV are likely not essential for BV virion formation. IMPORTANCE Bracoviruses are a remarkable example of virus endogenization, because large sets of genes from a nudivirus ancestor continue to produce virions that thousands of wasp species rely upon to parasitize hosts. Understanding how these genes interact and have been coopted by wasps for novel functions is of broad interest in the study of virus evolution. This work characterizes bracovirus genome components in the parasitoid wasp Chelonus insularis, which together with existing wasp genomes captures a large portion of the diversity among wasp species that produce bracoviruses. Results provide new information about how bracovirus genome components are organized in different wasps while also providing additional insights on key features required for function.

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.

Bracoviruses recruit host integrases for their integration into caterpillar's genome

Some DNA viruses infect host animals usually by integrating their DNAs into the host genome. However, the mechanisms for integration remain largely unknown. Here, we find that Cotesia vestalis bracovirus (CvBV), a polydnavirus of the parasitic wasp C. vestalis (Haliday), integrates its DNA circles into host Plutella xylostella (L.) genome by two distinct strategies, conservatively and randomly, through high-throughput sequencing analysis. We confirmed that the conservatively integrating circles contain an essential "8+5" nucleotides motif which is required for integration. Then we find CvBV circles are integrated into the caterpillar's genome in three temporal patterns, the early, mid and late stage-integration. We further identify that three CvBV-encoded integrases are responsible for some, but not all of the virus circle integrations, indeed they mainly participate in the processes of early stage-integration. Strikingly, we find two P. xylostella retroviral integrases (PxIN1 and PxIN2) are highly induced upon wasp parasitism, and PxIN1 is crucial for integration of some other early-integrated CvBV circles, such as CvBV_04, CvBV_12 and CvBV_24, while PxIN2 is important for integration of a late-integrated CvBV circle, CvBV_21. Our data uncover a novel mechanism in which CvBV integrates into the infected host genome, not only by utilizing its own integrases, but also by recruiting host enzymes. These findings will strongly deepen our understanding of how bracoviruses regulate and integrate into their hosts.

A Behavior-Manipulating Virus Relative as a Source of Adaptive Genes for Drosophila Parasitoids

Some species of parasitic wasps have domesticated viral machineries to deliver immunosuppressive factors to their hosts. Up to now, all described cases fall into the Ichneumonoidea superfamily, which only represents around 10% of hymenoptera diversity, raising the question of whether such domestication occurred outside this clade. Furthermore, the biology of the ancestral donor viruses is completely unknown. Since the 1980s, we know that Drosophila parasitoids belonging to the Leptopilina genus, which diverged from the Ichneumonoidea superfamily 225 Ma, do produce immunosuppressive virus-like structure in their reproductive apparatus. However, the viral origin of these structures has been the subject of debate. In this article, we provide genomic and experimental evidence that those structures do derive from an ancestral virus endogenization event. Interestingly, its close relatives induce a behavior manipulation in present-day wasps. Thus, we conclude that virus domestication is more prevalent than previously thought and that behavior manipulation may have been instrumental in the birth of such associations.

Influence of parasitoid-associated viral symbionts on plant-insect interactions and biological control

Insect parasitoids have evolved symbiotic interactions with several viruses and thousands of parasitoid species have established mutualistic associations with polydnaviruses (PDVs). While PDVs have often been described as virulence factors allowing development of immature parasitoids inside their herbivore hosts, there is increasing awareness that PDVs can affect plant-insect interactions. We review recent literature showing that PDVs alter not only host physiology, but also feeding patterns and composition of herbivore's oral secretions. In turn PDV-induced changes in herbivore phenotype affect plant responses to herbivory with consequences ranging from differential expression of plant defense-related genes to wider ecological effects across multiple trophic levels. In this opinion paper we also highlight important missing gaps to fully understand the role of PDVs and other parasitoid-associated viral symbionts in a plant-insect interaction perspective. Because PDVs negatively impact performance and survival of herbivore pests, we conclude arguing that PDV genomes offer potential opportunities for biological control.

Genome-Wide Profiling of Diadegma semiclausum Ichnovirus Integration in Parasitized Plutella xylostella Hemocytes Identifies Host Integration Motifs and Insertion Sites

Polydnaviruses (PDVs), classified into two genera, bracoviruses (BVs) and ichnoviruses (IVs), are large, double-stranded DNA viruses, which are beneficial symbionts of parasitoid wasps. PDVs do not replicate in their infected lepidopteran hosts. BV circles have been demonstrated to be integrated into host genomic DNA after natural parasitization. However, the integrations of IV circles in vivo remain largely unknown. Here, we analyzed the integration of Diadegma semiclausum ichnovirus (DsIV) in the genomic DNA of parasitized Plutella xylostella hemocytes. We found that DsIV circles are present in host hemocytes with non-integrated and integrated forms. Moreover, DsIV integrates its DNA circles into the host genome by two distinct strategies, conservatively, and randomly. We also found that four conserved-broken circles share similar motifs containing two reverse complementary repeats at their breaking sites, which were host integration motifs (HIMs). We also predicted HIMs of eight circles from other ichnoviruses, indicating that a HIM-mediated specific mechanism was conserved in IV integrations. Investigation of DsIV circle insertion sites of the host genome revealed the enrichment of microhomologies between the host genome and the DsIV circles at integration breakpoints. These findings will deepen our understanding of the infections of PDVs, especially IVs.

Detecting bracoviral orthologs distribution in five tsetse fly species and the housefly genomes

Objective

Mutualism between endogenous viruses and eukaryotes is still poorly understood. Several endogenous double-stranded polydnaviruses, bracoviruses, homologous to those present in parasitic braconid wasp genomes were detected in the tsetse fly (Glossina morsitans morsitans). This is peculiar since tsetse flies do not share a reproductive lifestyle similar to wasps, but deliver fully developed larvae that pupate within minutes of exiting their mothers. The objective of this study is to investigate genomic distribution of bracoviral sequences in five tsetse fly species and the housefly, and examine its value as a potential vector control strategy target point. We use comparative genomics to determine the presence, distribution across Glossina species genomes, and evolutionary relationships of bracoviruses of five tsetse fly species and the housefly.

Results

We report on homologous bracoviruses in multiple Dipteran genomes. Phylogenetic reconstruction using within-species concatenated bracoviral orthologs shows great congruence with previously reconstructed insect species phylogenies. Our findings suggest that bracoviruses present in Diptera originate from a single integration event of the viral genome that occurred in an ancestor insect before the evolutionary radiation of different insect orders.

Role of the DNA-Binding Protein pA104R in ASFV Genome Packaging and as a Novel Target for Vaccine and Drug Development

The recent incursions of African swine fever (ASF), a severe, highly contagious, transboundary viral disease that affects members of the Suidae family, in Europe and China have had a catastrophic impact on trade and pig production, with serious implications for global food security. Despite efforts made over past decades, there is no vaccine or treatment available for preventing and controlling the ASF virus (ASFV) infection, and there is an urgent need to develop novel strategies. Genome condensation and packaging are essential processes in the life cycle of viruses. The involvement of viral DNA-binding proteins in the regulation of virulence genes, transcription, DNA replication, and repair make them significant targets. pA104R is a highly conserved HU/IHF-like DNA-packaging protein identified in the ASFV nucleoid that appears to be profoundly involved in the spatial organization and packaging of the ASFV genome. Here, we briefly review the components of the ASFV packaging machinery, the structure, function, and phylogeny of pA104R, and its potential as a target for vaccine and drug development.

Quantitative trait loci involved in the reproductive success of a parasitoid wasp

Dissecting the genetic basis of intraspecific variations in life history traits is essential to understand their evolution, notably for potential biocontrol agents. Such variations are observed in the endoparasitoid Cotesia typhae (Hymenoptera: Braconidae), specialized on the pest Sesamia nonagrioides (Lepidoptera: Noctuidae). Previously, we identified two strains of C. typhae that differed significantly for life history traits on an allopatric host population. To investigate the genetic basis underlying these phenotypic differences, we used a quantitative trait locus (QTL) approach based on restriction site-associated DNA markers. The characteristic of C. typhae reproduction allowed us generating sisters sharing almost the same genetic content, named clonal sibship. Crosses between individuals from the two strains were performed to generate F2 and F8 recombinant CSS. The genotypes of 181 clonal sibships were determined as well as the phenotypes of the corresponding 4,000 females. Informative markers were then used to build a high-quality genetic map. These 465 markers spanned a total length of 1,300 cM and were organized in 10 linkage groups which corresponded to the number of C. typhae chromosomes. Three QTLs were detected for parasitism success and two for offspring number, while none were identified for sex ratio. The QTLs explained, respectively, 27.7% and 24.5% of the phenotypic variation observed. The gene content of the genomic intervals was investigated based on the genome of C. congregata and revealed 67 interesting candidates, as potentially involved in the studied traits, including components of the venom and of the symbiotic virus (bracovirus) shown to be necessary for parasitism success in related wasps.

Regulation of the Larval Transcriptome of Diatraea saccharalis (Lepidoptera: Crambidae) by Maternal and Other Factors of the Parasitoid Cotesia flavipes (Hymenoptera: Braconidae)

Koinobiont endoparasitoid wasps regulate the host's physiology to their own benefit during their growth and development, using maternal, immature and/or derived-tissue weaponry. The tools used to subdue the wasps' hosts interfere directly with host transcription activity. The broad range of host tissues and pathways affected impedes our overall understanding of the host-regulation process during parasitoid development. Next-generation sequencing and de novo transcriptomes are helpful approaches to broad questions, including in non-model organisms. In the present study, we used Illumina sequencing to assemble a de novo reference transcriptome of the sugarcane borer Diatraea saccharalis, to investigate the regulation of host gene expression by the larval endoparasitoid Cotesia flavipes. We obtained 174,809,358 reads and assembled 144,116 transcripts, of which 44,325 were putatively identified as lepidopteran genes and represented a substantial number of pathways that are well described in other lepidopteran species. Comparative transcriptome analyses of unparasitized versus parasitized larvae identified 1,432 transcripts of D. saccharalis that were up-regulated under parasitization by C. flavipes, while 1,027 transcripts were down-regulated. Comparison of the transcriptomes of unparasitized and pseudoparasitized D. saccharalis larvae led to the identification of 1,253 up-regulated transcripts and 972 down-regulated transcripts in the pseudoparasitized larvae. Analysis of the differentially expressed transcripts showed that C. flavipes regulated several pathways, including the Ca+2 transduction signaling pathway, glycolysis/gluconeogenesis, chitin metabolism, and hormone biosynthesis and degradation, as well as the immune system, allowing us to identify key target genes involved in the metabolism and development of D. saccharalis.

Gall Wasp Transcriptomes Unravel Potential Effectors Involved in Molecular Dialogues With Oak and Rose

To gain insight into wasp factors that might be involved in the initial induction of galls on woody plants, we performed high throughput (454) transcriptome analysis of ovaries and venom glands of two cynipid gall wasps, Biorhiza pallida and Diplolepis rosae, inducing galls on oak and rose, respectively. De novo assembled and annotated contigs were compared to sequences from phylogenetically related parasitoid wasps. The relative expression levels of contigs were estimated to identify the most expressed gene sequences in each tissue. We identify for the first time a set of maternally expressed gall wasp proteins potentially involved in the interaction with the plant. Some genes highly expressed in venom glands and ovaries may act to suppress early plant defense signaling. We also identify gall wasp cellulases that could be involved in observed local lysis of plant tissue following oviposition, and which may have been acquired from bacteria by horizontal gene transfer. We find no evidence of virus-related gene expression, in contrast to many non-cynipid parasitoid wasps. By exploring gall wasp effectors, this study is a first step toward understanding the molecular mechanisms underlying cynipid gall induction in woody plants, and the recent sequencing of oak and rose genomes will enable study of plant responses to these factors.

An Ovarian Protein Involved in Passive Avoidance of an Endoparasitoid To Evade Its Host Immune Response

Through a combination of transcriptomic and proteomic analyses, we identified 817 secreted ovarian proteins from an endoparasitoid wasp, Cotesia chilonis, of which five proteins are probably involved in passive evasion. The results of an encapsulation assay revealed that one of these passive evasion-associated proteins (Crp32B), a homologue of a 32-kDa protein (Crp32) from C. rubecula, could protect resin beads from being encapsulated by host hemocytes in a dose-dependent manner. Crp32B is transcribed in ovarian cells, nurse cells, follicular cells, and oocytes, and the protein is located throughout the ovary and on the egg surface. Moreover, Crp32B has antigenic similarity to several host components. These results indicate that C. chilonis may use molecular mimicry as a mechanism to avoid host cellular immune response.

Evolution of an insect immune barrier through horizontal gene transfer mediated by a parasitic wasp

Genome sequencing data have recently demonstrated that eukaryote evolution has been remarkably influenced by the acquisition of a large number of genes by horizontal gene transfer (HGT) across different kingdoms. However, in depth-studies on the physiological traits conferred by these accidental DNA acquisitions are largely lacking. Here we elucidate the functional role of Sl gasmin, a gene of a symbiotic virus of a parasitic wasp that has been transferred to an ancestor of the moth species Spodoptera littoralis and domesticated. This gene is highly expressed in circulating immune cells (haemocytes) of larval stages, where its transcription is rapidly boosted by injection of microorganisms into the body cavity. RNAi silencing of Sl gasmin generates a phenotype characterized by a precocious suppression of phagocytic activity by haemocytes, which is rescued when these immune cells are incubated in plasma samples of control larvae, containing high levels of the encoded protein. Proteomic analysis demonstrates that the protein Sl gasmin is released by haemocytes into the haemolymph, where it opsonizes the invading bacteria to promote their phagocytosis, both in vitro and in vivo. Our results show that important physiological traits do not necessarily originate from evolution of pre-existing genes, but can be acquired by HGT events, through unique pathways of symbiotic evolution. These findings indicate that insects can paradoxically acquire selective advantages with the help of their natural enemies.

Evolutionary relationships of courtship songs in the parasitic wasp genus, Cotesia (Hymenoptera: Braconidae)

Acoustic signals play an important role in premating isolation based on sexual selection within many taxa. Many male parasitic wasps produce characteristic courtship songs used by females in mate selection. In Cotesia (Hymenoptera: Braconidae: Microgastrinae), courtship songs are generated by wing fanning with repetitive pulses in stereotypical patterns. Our objectives were to sample the diversity of courtship songs within Cotesia and to identify e underlying patterns of differentiation. We compared songs among 12 of ca. 80 Cotesia species in North America, including ten species that have not been recorded previously. For Cotesia congregata, we compared songs of wasps originating from six different host-foodplant sources, two of which are considered incipient species. Songs of emergent males from wild caterpillar hosts in five different families were recorded, and pattern, frequency, and duration of song elements analyzed. Principal component analysis converted the seven elements characterized into four uncorrelated components used in a hierarchical cluster analysis and grouped species by similarity of song structure. Species songs varied significantly in duration of repeating pulse and buzz elements and/or in fundamental frequency. Cluster analysis resolved similar species groups in agreement with the most recent molecular phylogeny for Cotesia spp., indicating the potential for using courtship songs as a predictor of genetic relatedness. Courtship song analysis may aid in identifying closely related cryptic species that overlap spatially, and provide insight into the evolution of this highly diverse and agriculturally important taxon.

The Domestication of a Large DNA Virus by the Wasp Venturia canescens Involves Targeted Genome Reduction through Pseudogenization

Polydnaviruses (PDVs) are compelling examples of viral domestication, in which wasps express a large set of genes originating from a chromosomally integrated virus to produce particles necessary for their reproductive success. Parasitoid wasps generally use PDVs as a virulence gene delivery system allowing the protection of their progeny in the body of parasitized host. However, in the wasp Venturia canescens an independent viral domestication process led to an alternative strategy as the wasp incorporates virulence proteins in viral liposomes named virus-like particles (VLPs), instead of DNA molecules. Proteomic analysis of purified VLPs and transcriptome sequencing revealed the loss of some viral functions. In particular, the genes coding for capsid components are no longer expressed, which explains why VLPs do not incorporate DNA. Here a thorough examination of V. canescens genome revealed the presence of the pseudogenes corresponding to most of the genes involved in lost functions. This strongly suggests that an accumulation of mutations that leads to gene specific pseudogenization precedes the loss of viral genes observed during virus domestication. No evidence was found for block loss of collinear genes, although extensive gene order reshuffling of the viral genome was identified from comparisons between endogenous and exogenous viruses. These results provide the first insights on the early stages of large DNA virus domestication implicating massive genome reduction through gene-specific pseudogenization, a process which differs from the large deletions described for bacterial endosymbionts.

Cotesia congregata Bracovirus Circles Encoding PTP and Ankyrin Genes Integrate into the DNA of Parasitized Manduca sexta Hemocytes

Polydnaviruses (PDVs) are essential for the parasitism success of tens of thousands of species of parasitoid wasps. PDVs are present in wasp genomes as proviruses, which serve as the template for the production of double-stranded circular viral DNA carrying virulence genes that are injected into lepidopteran hosts. PDV circles do not contain genes coding for particle production, thereby impeding viral replication in caterpillar hosts during parasitism. Here, we investigated the fate of PDV circles of Cotesia congregata bracovirus during parasitism of the tobacco hornworm, Manduca sexta, by the wasp Cotesia congregata Sequences sharing similarities with host integration motifs (HIMs) of Microplitis demolitor bracovirus (MdBV) circles involved in integration into DNA could be identified in 12 CcBV circles, which encode PTP and VANK gene families involved in host immune disruption. A PCR approach performed on a subset of these circles indicated that they persisted in parasitized M. sexta hemocytes as linear forms, possibly integrated in host DNA. Furthermore, by using a primer extension capture method based on these HIMs and high-throughput sequencing, we could show that 8 out of 9 circles tested were integrated in M. sexta hemocyte genomic DNA and that integration had occurred specifically using the HIM, indicating that an HIM-mediated specific mechanism was involved in their integration. Investigation of BV circle insertion sites at the genome scale revealed that certain genomic regions appeared to be enriched in BV insertions, but no specific M. sexta target site could be identified.IMPORTANCE The identification of a specific and efficient integration mechanism shared by several bracovirus species opens the question of its role in braconid parasitoid wasp parasitism success. Indeed, results obtained here show massive integration of bracovirus DNA in somatic immune cells at each parasitism event of a caterpillar host. Given that bracoviruses do not replicate in infected cells, integration of viral sequences in host DNA might allow the production of PTP and VANK virulence proteins within newly dividing cells of caterpillar hosts that continue to develop during parasitism. Furthermore, this integration process could serve as a basis to understand how PDVs mediate the recently identified gene flux between parasitoid wasps and Lepidoptera and the frequency of these horizontal transfer events in nature.

Persistent expression of Cotesia plutellae bracovirus genes in parasitized host, Plutella xylostella

Cotesia plutellae (= vestalis) bracovirus (CpBV) is symbiotic to an endoparasitoid wasp, C. plutellae, and plays crucial roles in parasitism against the diamondback moth, Plutella xylostella. CpBV virion genome consists of 35 circular DNAs encoding 157 putative open reading frames (ORFs). This study re-annotated 157 ORFs with update genome database and analyzed their gene expressions at early and late parasitic stages. Re-annotation has established 15 different viral gene families, to which 83 ORFs are assigned with remaining 74 hypothetical genes. Among 157 ORFs, 147 genes were expressed at early or late parasitic stages, among which 141 genes were expressed in both parasitic stages, indicating persistent nature of gene expression. Relative frequencies of different viral circles present in the ovarian lumen did not explain the expression variation of the viral ORFs. Furthermore, expression level of each viral gene was varied during parasitism along with host development. Highly up-regulated CpBV genes at early parasitic stage included BEN (BANP, E₅R and NAC₁), ELP (EP1-like protein), IkB (inhibitor kB), P494 (protein 494 kDa) family genes, while those at late stage were mostly hypothetical genes. Along with the viral gene expression, 362 host genes exhibited more than two fold changes in expression levels at early parasitic stage compared to nonparasitized host. At late stage, more number (1,858) of host genes was regulated. These results suggest that persistent expression of most CpBV genes may be necessary to regulate host physiological processes during C. plutellae parasitism.

Impact of Lateral Transfers on the Genomes of Lepidoptera

Transfer of DNA sequences between species regardless of their evolutionary distance is very common in bacteria, but evidence that horizontal gene transfer (HGT) also occurs in multicellular organisms has been accumulating in the past few years. The actual extent of this phenomenon is underestimated due to frequent sequence filtering of "alien" DNA before genome assembly. However, recent studies based on genome sequencing have revealed, and experimentally verified, the presence of foreign DNA sequences in the genetic material of several species of Lepidoptera. Large DNA viruses, such as baculoviruses and the symbiotic viruses of parasitic wasps (bracoviruses), have the potential to mediate these transfers in Lepidoptera. In particular, using ultra-deep sequencing, newly integrated transposons have been identified within baculovirus genomes. Bacterial genes have also been acquired by genomes of Lepidoptera, as in other insects and nematodes. In addition, insertions of bracovirus sequences were present in the genomes of certain moth and butterfly lineages, that were likely corresponding to rearrangements of ancient integrations. The viral genes present in these sequences, sometimes of hymenopteran origin, have been co-opted by lepidopteran species to confer some protection against pathogens.

The recurrent domestication of viruses: major evolutionary transitions in parasitic wasps

Several lineages of endoparasitoid wasps, which develop inside the body of other insects, have domesticated viruses, used as delivery tools of essential virulence factors for the successful development of their progeny. Virus domestications are major evolutionary transitions in highly diverse parasitoid wasps. Much progress has recently been made to characterize the nature of these ancestrally captured endogenous viruses that have evolved within the wasp genomes. Virus domestication from different viral families occurred at least three times in parasitoid wasps. This evolutionary convergence led to different strategies. Polydnaviruses (PDVs) are viral gene transfer agents and virus-like particles of the wasp Venturia canescens deliver proteins. Here, we take the standpoint of parasitoid wasps to review current knowledge on virus domestications by different parasitoid lineages. Then, based on genomic data from parasitoid wasps, PDVs and exogenous viruses, we discuss the different evolutionary steps required to transform viruses into vehicles for the delivery of the virulence molecules that we observe today. Finally, we discuss how endoparasitoid wasps manipulate host physiology and ensure parasitism success, to highlight the possible advantages of viral domestication as compared with other virulence strategies.

Translational Control of Host Gene Expression by a Cys-Motif Protein Encoded in a Bracovirus

Translational control is a strategy that various viruses use to manipulate their hosts to suppress acute antiviral response. Polydnaviruses, a group of insect double-stranded DNA viruses symbiotic to some endoparasitoid wasps, are divided into two genera: ichnovirus (IV) and bracovirus (BV). In IV, some Cys-motif genes are known as host translation-inhibitory factors (HTIF). The genome of endoparasitoid wasp Cotesia plutellae contains a Cys-motif gene (Cp-TSP13) homologous to an HTIF known as teratocyte-secretory protein 14 (TSP14) of Microplitis croceipes. Cp-TSP13 consists of 129 amino acid residues with a predicted molecular weight of 13.987 kDa and pI value of 7.928. Genomic DNA region encoding its open reading frame has three introns. Cp-TSP13 possesses six conserved cysteine residues as other Cys-motif genes functioning as HTIF. Cp-TSP13 was expressed in Plutella xylostella larvae parasitized by C. plutellae. C. plutellae bracovirus (CpBV) was purified and injected into non-parasitized P. xylostella that expressed Cp-TSP13. Cp-TSP13 was cloned into a eukaryotic expression vector and used to infect Sf9 cells to transiently express Cp-TSP13. The synthesized Cp-TSP13 protein was detected in culture broth. An overlaying experiment showed that the purified Cp-TSP13 entered hemocytes. It was localized in the cytosol. Recombinant Cp-TSP13 significantly inhibited protein synthesis of secretory proteins when it was added to in vitro cultured fat body. In addition, the recombinant Cp-TSP13 directly inhibited the translation of fat body mRNAs in in vitro translation assay using rabbit reticulocyte lysate. Moreover, the recombinant Cp-TSP13 significantly suppressed cellular immune responses by inhibiting hemocyte-spreading behavior. It also exhibited significant insecticidal activities by both injection and feeding routes. These results indicate that Cp-TSP13 is a viral HTIF.

Characterization of a venom gland-associated rhabdovirus in the parasitoid wasp Diachasmimorpha longicaudata

Parasitoid wasps reproduce by laying their eggs on or inside of a host insect, which triggers a defense response in the host insect that kills the developing wasp. To counteract the host's lethal response, some parasitoid wasps are associated with symbiotic viruses that alter host metabolism and development to promote successful development of the wasp embryo. These symbiotic viruses display a number of characteristics that differ from those of pathogenic viruses, but are poorly understood with the exception of one group, the polydnaviruses. Here, we characterize the genome of a non-polydnavirus associated with parasitoid wasps, Diachasmimorpha longicaudata rhabdovirus (DlRhV), and assess its role as a potential mutualistic virus. Our results show that the DlRhV genome contains six open reading frames (ORFs). Three ORFs show sequence homology to known viral genes and one ORF encodes a previously identified protein, called parasitism-specific protein 24 (PSP24), that has been hypothesized to play a role in promoting successful parasitism by D. longicaudata. We constructed a phylogeny that shows that DlRhV is most closely related to other insect-infecting rhabdoviruses. Finally, we report that DlRhV infection does not occur in all populations of D. longicaudata, and is not required for successful parasitism.

Glyceraldehyde-3-phosphate dehydrogenase is a mediator of hemocyte-spreading behavior and molecular target of immunosuppressive factor CrV1

Cellular immunity is accompanied by hemocyte-spreading behavior, which undergoes cytoskeletal rearrangement. Polydnaviral factor CpBV-CrV1 can inhibit the hemocyte-spreading behavior and suppress host immune response of Plutella xylostella. However, host target molecule of CpBV-CrV1 that inhibits the hemocyte behavior has not been identified yet. This study used a pull-down approach to identify the target molecule of CpBV-CrV1. A protein bound to CpBV-CrV1 was co-precipitated and identified to be glyceraldehyde-3-phosphate dehydrogenase (GAPDH) by LC-MS/MS analysis. RNA interference (RNAi) specific to GAPDH of P. xylostella was found to be able to inhibit the hemocyte-spreading behavior, while RNAi treatments with other glycolytic genes had no effect on the spreading behavior. An addition of recombinant CpBV-CrV1 on hemocyte monolayer interrupted the association between GAPDH and α-tubulin in the cytoplasm. Overlay of mutant proteins (Y492A or Y501A with tyrosine to alanine at putative GAPDH-binding site) of CpBV-CrV1 on hemocyte monolayer revealed that they could enter hemocytes unlike a mutant in the N-terminal coiled-coil domain. However, they failed to inhibit the hemocyte-spreading behavior without any binding affinity to GAPDH. These results suggest that GAPDH plays a critical role in hemocyte-spreading behavior during immune challenge as a molecular target of viral factor CpBV-CrV1.

Recurrent Domestication by Lepidoptera of Genes from Their Parasites Mediated by Bracoviruses

Bracoviruses are symbiotic viruses associated with tens of thousands of species of parasitic wasps that develop within the body of lepidopteran hosts and that collectively parasitize caterpillars of virtually every lepidopteran species. Viral particles are produced in the wasp ovaries and injected into host larvae with the wasp eggs. Once in the host body, the viral DNA circles enclosed in the particles integrate into lepidopteran host cell DNA. Here we show that bracovirus DNA sequences have been inserted repeatedly into lepidopteran genomes, indicating this viral DNA can also enter germline cells. The original mode of Horizontal Gene Transfer (HGT) unveiled here is based on the integrative properties of an endogenous virus that has evolved as a gene transfer agent within parasitic wasp genomes for ≈100 million years. Among the bracovirus genes thus transferred, a phylogenetic analysis indicated that those encoding C-type-lectins most likely originated from the wasp gene set, showing that a bracovirus-mediated gene flux exists between the 2 insect orders Hymenoptera and Lepidoptera. Furthermore, the acquisition of bracovirus sequences that can be expressed by Lepidoptera has resulted in the domestication of several genes that could result in adaptive advantages for the host. Indeed, functional analyses suggest that two of the acquired genes could have a protective role against a common pathogen in the field, baculovirus. From these results, we hypothesize that bracovirus-mediated HGT has played an important role in the evolutionary arms race between Lepidoptera and their pathogens.

Eukaryotes are generally thought to evolve mainly through the modification of existing genetic information. However, evidence of horizontal gene transfer (HGT) in eukaryotes-the accidental acquisition of a novel gene from another species, allowing acquisition of novel traits—is now recognized as an important factor in their evolution. We show here that in several lineages, lepidopteran genomes have acquired genes from a bracovirus that is symbiotically used by parasitic wasps to inhibit caterpillar host immune defences. Integration of parts of the viral genome into host caterpillar DNA strongly suggests that integration can sporadically occur in the germline, leading to the production of lepidopteran lineages that harbor bracovirus sequences. Moreover, some of the transferred bracovirus genes reported here originate from the wasp genome, demonstrating that a gene flux exists between the two insect orders Hymenoptera and Lepidoptera that diverged ≈300 MYA. As bracovirus gene organisation has evolved to allow expression in Lepidoptera, these transferred genes can be readily domesticated. Additionally, we present functional analyses suggesting that some of the acquired genes confer to caterpillars a protection toward baculovirus, a very common pathogen in the field. This phenomenon may have implications for understanding how caterpillars acquire resistance against baculoviruses used in biological control.

Transcriptomic response of Manduca sexta immune tissues to parasitization by the bracovirus associated wasp Cotesia congregata

During oviposition, Cotesia congregata parasitoid wasps inject into their host, Manduca sexta, some biological factors such as venom, ovarian fluid and a symbiotic polydnavirus (PDV) named Cotesia congregata bracovirus (CcBV). During parasitism, complex interactions occur between wasp-derived factors and host targets that lead to important modifications in host physiology. In particular, the immune response leading to wasp egg encapsulation is inhibited allowing wasp survival. To date, the regulation of host genes during the interaction had only been studied for a limited number of genes. In this study, we analysed the global impact of parasitism on host gene regulation 24 h post oviposition by high throughput 454 transcriptomic analyses of two tissues known to be involved in the host immune response (hemocytes and fat body). To identify specific effects of parasitism on host transcription at this time point, transcriptomes were obtained from non-treated and parasitized larvae, and also from larvae injected with heat-killed bacteria and double stimulated larvae that were parasitized prior to bacterial challenge. Results showed that, immune challenge by bacteria leads to induction of certain antimicrobial peptide (AMP) genes in M. sexta larvae whether they were parasitized or not prior to bacterial challenge. These results show that at 24 h post oviposition pathways leading to expression of AMP genes are not all inactivated suggesting wasps are in an antiseptic environment. In contrast, at this time point genes involved in phenoloxidase activation and cellular immune responses were globally down-regulated after parasitism in accordance with the observed inhibition of wasp egg encapsulation.

Polydnaviruses: From discovery to current insights

The International Committee on Taxonomy of Viruses (ICTV) recognized the Polydnaviridae in 1991 as a virus family associated with insects called parasitoid wasps. Polydnaviruses (PDVs) have historically received limited attention but advances in recent years have elevated interest because their unusual biology sheds interesting light on the question of what viruses are and how they function. Here, we present a succinct history of the PDV literature. We begin with the findings that first led ICTV to recognize the Polydnaviridae. We then discuss what subsequent studies revealed and how these findings have shaped views of PDV evolution.

Origin and evolution of symbiotic viruses associated with parasitoid wasps

The Polydnaviridae (PDV), including the Bracovirus (BV) and Ichnovirus (IV) genera, originated from the integration of viruses in the genomes of two parasitoid wasp lineages. In a remarkable example of convergent evolution BVs evolved from the domestication of a nudivirus, while IVs originate from a different ancestral virus belonging to a new virus entity. In both cases the ancestor genomes have been maintained in wasp genomes as endogenous viral elements involved in production of particles containing DNA encoding virulence genes that are injected into lepidopteran hosts. However many PDV virulence genes appear to be of eukaryotic origin, and expansion and diversification of these genes have led to the production of novel PDVs in different wasp species that promote survival of offspring in particular hosts.