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.

Massive Somatic and Germline Chromosomal Integrations of Polydnaviruses in Lepidopterans

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

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

RNA interference identifies domesticated viral genes involved in assembly and trafficking of virus-derived particles in ichneumonid wasps

There are many documented examples of viral genes retained in the genomes of multicellular organisms that may in some cases bring new beneficial functions to the receivers. The ability of certain ichneumonid parasitic wasps to produce virus-derived particles, the so-called ichnoviruses (IVs), not only results from the capture and domestication of single viral genes but of almost entire ancestral virus genome(s). Indeed, following integration into wasp chromosomal DNA, the putative and still undetermined IV ancestor(s) evolved into encoding a ‘virulence gene delivery vehicle’ that is now required for successful infestation of wasp hosts. Several putative viral genes, which are clustered in distinct regions of wasp genomes referred to as IVSPERs (Ichnovirus Structural Protein Encoding Regions), have been assumed to be involved in virus-derived particles morphogenesis, but this question has not been previously functionally addressed. In the present study, we have successfully combined RNA interference and transmission electron microscopy to specifically identify IVSPER genes that are responsible for the morphogenesis and trafficking of the virus-derived particles in ovarian cells of the ichneumonid wasp Hyposoter didymator. We suggest that ancestral viral genes retained within the genomes of certain ichneumonid parasitoids possess conserved functions which were domesticated for the purpose of assembling viral vectors for the delivery of virulence genes to parasitized host animals.

Thousands of parasitic wasp from the ichneumonid family rely on virus-derived particles, named Ichnoviruses (Polydnavirus family), to ensure their successful development. The particles are produced in a specialized ovarian tissue of the female wasp named calyx. Virions are assembled in the calyx cell nuclei and stored in the oviduct before being transferred to the parasitoid host upon female wasp oviposition. Genes encoding proteins associated with the particles had been previously identified. These genes are localized in clusters of genes in the wasp genome (named IVSPER for “Ichnovirus structural proteins encoding regions”), they are specifically transcribed in the calyx but not encapsidated. IVSPER genes were thus hypothesized to derive from the integration of a virus, however still undetermined. Indeed, none of the identified genes had similarity to known sequence, making in addition unclear their function in particle production. In this work, we use the RNA interference technology to decipher the function of six IVSPER genes from the ichneumonid wasp Hyposoter didymator. Thanks to this approach, combined with transmission electron microscopy, we show that the studied IVSPER genes are required in different steps of particle morphogenesis and trafficking, and that their functions are those expected of a typical virus.

First extensive characterization of the venom gland from an egg parasitoid: structure, transcriptome and functional role

The venom gland is a ubiquitous organ in Hymenoptera. In insect parasitoids, the venom gland has been shown to have multiple functions including regulation of host immune response, host paralysis, host castration and developmental alteration. However, the role played by the venom gland has been mainly studied in parasitoids developing in larval or pupal hosts while little is known for parasitoids developing in insect eggs. We conducted the first extensive characterization of the venom of the endoparasitoid Ooencyrtus telenomicida (Vassiliev), a species that develops in eggs of the stink bug Nezara viridula (L.). In particular we investigated the structure of the venom apparatus, its functional role and conducted a transcriptomic analysis of the venom gland. We found that injection of O. telenomicida venom induces: 1) a melanized-like process in N. viridula host eggs (host-parasitoid interaction), 2) impairment of the larval development of the competitor Trissolcus basalis (Wollaston) (parasitoid-parasitoid interaction). The O. telenomicida venom gland transcriptome reveals a majority of digestive enzymes (peptidases and glycosylases) and oxidoreductases (laccases) among the most expressed genes. The former enzymes are likely to be involved in degradation of the host resources for the specific benefit of the O. telenomicida offspring. In turn, alteration of host resources caused by these enzymes may negatively affect the larval development of the competitor T. basalis. We hypothesize that the melanization process induced by venom injection could be related to the presence of laccases, which are multicopper oxidases that belong to the phenoloxidases group. This work contributed to a better understanding of the venom in insect parasitoids and allowed to identify candidate genes whose functional role can be investigated in future studies.

Recurrent DNA virus domestication leading to different parasite virulence strategies

Relics of ancient infections are abundant in eukaryote genomes, but little is known about how they evolve when they confer a functional benefit on their host. We show here, for the first time, that the virus-like particles shown to protect Venturia canescens eggs against host immunity are derived from a nudivirus genome incorporated by the parasitic wasp into its own genetic material. Nudivirus hijacking was also at the origin of protective particles from braconid wasps. However, we show here that the viral genes produce "liposomes" that wrap and deliver V. canescens virulence proteins, whereas the particles are used as gene transfer agents in braconid wasps. Our findings indicate that virus domestication has occurred repeatedly during parasitic wasp evolution but with different evolutionary trajectories after endogenization, resulting in different virulence molecule delivery strategies.

Specificities of ichnoviruses associated with campoplegine wasps: genome, genes and role in host-parasitoid interaction

Ichnoviruses (IVs), unique symbiotic viruses carried by ichneumonid campoplegine wasps, derive from integration of a paleo-ichnovirus into an ancestral wasp genome. The modern 'genome' is composed of both regions that are amplified, circularized and encapsidated into viral particles and non-encapsidated viral genomic regions involved in particle morphogenesis. Packaged genomes include multiple circular dsDNAs encoding many genes mostly organized in gene families. Virus particles are assembled in specialized ovarian cells from which they exit into the oviduct lumen; mature virions are injected during oviposition into the insect host. Expression of viral proteins in infected cells correlates with physiological alterations of the host enabling success of parasitism.

Extensive transcription analysis of the Hyposoter didymator Ichnovirus genome in permissive and non-permissive lepidopteran host species

Ichnoviruses are large dsDNA viruses that belong to the Polydnaviridae family. They are specifically associated with endoparasitic wasps of the family Ichneumonidae and essential for host parasitization by these wasps. We sequenced the Hyposoter didymator Ichnovirus (HdIV) encapsidated genome for further analysis of the transcription pattern of the entire set of HdIV genes following the parasitization of four different lepidopteran host species. The HdIV genome was found to consist of at least 50 circular dsDNA molecules, carrying 135 genes, 98 of which formed 18 gene families. The HdIV genome had general features typical of Ichnovirus (IV) genomes and closely resembled that of the IV carried by Hyposoter fugitivus. Subsequent transcriptomic analysis with Illumina technology during the course of Spodoptera frugiperda parasitization led to the identification of a small subset of less than 30 genes with high RPKM values in permissive hosts, consisting with these genes encoding crucial virulence proteins. Comparisons of HdIV expression profiles between host species revealed differences in transcript levels for given HdIV genes between two permissive hosts, S. frugiperda and Pseudoplusia includens. However, we found no evident intrafamily gene-specific transcription pattern consistent with the presence of multigenic families within IV genomes reflecting an ability of the wasps concerned to exploit different host species. Interestingly, in two non-permissive hosts, Mamestra brassiccae and Anticarsia gemmatalis (most of the parasitoid eggs were eliminated by the host cellular immune response), HdIV genes were generally less strongly transcribed than in permissive hosts. This suggests that successful parasitism is dependent on the expression of given HdIV genes exceeding a particular threshold value. These results raise questions about the mecanisms involved in regulating IV gene expression according to the nature of the lepidopteran host species encountered.

Hyposoter didymator uses a combination of passive and active strategies to escape from the Spodoptera frugiperda cellular immune response

An endoparasitic life style is widespread among Hymenoptera, and various different strategies allowing parasitoids to escape from the host encapsulation response have been reported. Species carrying polydnaviruses (PDVs), such as the ichneumonid Hyposoter didymator, generally rely on the viral symbionts to evade host immune responses. In this work, we show that H. didymator eggs can evade encapsulation by the host in the absence of calyx fluid (containing the viral particles), whereas protection of the larvae requires the presence of calyx fluid. This evasion by the eggs depends on proteins associated with the exochorion. This type of local passive strategy has been described for a few species carrying PDVs. Immune evasion by braconid eggs appears to be related to PDVs or proteins synthesized in the oviducts being associated with the egg. We report that in H. didymator, by contrast, proteins already present in the ovarian follicles are responsible for the eggs avoiding encapsulation. Mass spectrometry analysis of the egg surface proteins revealed the presence of host immune-related proteins, including one with similarities with apolipophorin-III, and also the presence of three viral proteins encoded by IVSPERs (Ichnovirus Structural Protein Encoding Regions).

Venom gland extract is not required for successful parasitism in the polydnavirus-associated endoparasitoid Hyposoter didymator (Hym. Ichneumonidae) despite the presence of numerous novel and conserved venom proteins

The venom gland is a conserved organ in Hymenoptera that shows adaptations associated with life-style diversification. Few studies have investigated venom components and function in the highly diverse parasitic wasps and all suggest that the venom regulates host physiology. We explored the venom of the endoparasitoid Hyposoter didymator (Campopleginae), a species with an associated polydnavirus produced in the ovarian tissue. We investigated the effects of the H. didymator venom on two physiological traits of the host Spodoptera frugiperda (Noctuidae): encapsulation response and growth rate. We found that H. didymator venom had no significant effect on host cellular immunity or development, suggesting that it does not contribute to parasitism success. The host physiology seemed to be modified essentially by the ovarian fluid containing the symbiotic polydnaviruses. Proteomic analyses indicated that the H. didymator venom gland produces a large variety of proteins, consistent with the classical hymenopteran venom protein signature, including: reprolysin-like, dipeptidyl peptidase IV, hyaluronidase, arginine kinase or allergen proteins. The venom extracts also contained novel proteins, encoded by venom genes conserved in Campopleginae ichneumonids, and proteins with similarities to active molecules identified in other parasitoid species, such as calreticulin, reprolysin, superoxide dismutase and serpin. However, some of these proteins appear to be produced only in small amounts or to not be secreted. Possibly, in Campopleginae carrying polydnaviruses, the host-modifying activities of venom became redundant following the acquisition of polydnaviruses by the lineage.

Cecropins as a marker of Spodoptera frugiperda immunosuppression during entomopathogenic bacterial challenge

An antimicrobial peptide (AMP) of the cecropin family was isolated by HPLC from plasma of the insect pest, Spodoptera frugiperda. Its molecular mass is 3910.9 Da as determined by mass spectrometry. Thanks to the EST database Spodobase, we were able to describe 13 cDNAs encoding six different cecropins which belong to the sub-families CecA, CecB, CecC and CecD. The purified peptide identified as CecB1 was chemically synthesized (syCecB1). It was shown to be active against Gram-positive and Gram-negative bacteria as well as fungi. Two closely related entomopathogenic bacteria, Xenorhabdus nematophila F1 and Xenorhabdus mauleonii VC01(T) showed different susceptibility to syCecB1. Indeed, X. nematophila was sensitive to syCecB1 whereas X. mauleonii had a minimal inhibitory concentration (MIC) eight times higher. Interestingly, injection of live X. nematophila into insects did not induce the expression of AMPs in hemolymph. This effect was not observed when this bacterium was heat-killed before injection. On the opposite, both live and heat-killed X. mauleonii induced the expression of AMPs in the hemolymph of S. frugiperda. The same phenomenon was observed for another immune-related protein lacking antimicrobial activity. Altogether, our data suggest that Xenorhabdus strains have developed different strategies to supplant the humoral defense mechanisms of S. frugiperda, either by increasing their resistance to AMPs or by preventing their expression during such host-pathogen interaction.

Multigenic families in Ichnovirus: a tissue and host specificity study through expression analysis of vankyrins from Hyposoter didymator Ichnovirus

The viral ankyrin (vankyrin) gene family is represented in all polydnavirus (PDVs) genomes and encodes proteins homologous to I-kappaBs, inhibitors of NF-kappaB transcription factors. The structural similarities led to the hypothesis that vankyrins mimic eukaryotic factors to subvert important physiological pathways in the infected host. Here, we identified nine vankyrin genes in the genome of the Hyposoter didymator Ichnovirus (HdIV). Time-course gene expression experiments indicate that all members are expressed throughout parasitism of Spodoptera frugiperda, as assessed using RNA extracted from whole larvae. To study tissue and/or species specificity transcriptions, the expression of HdIV vankyrin genes were compared between HdIV-injected larvae of S. frugiperda and S. littoralis. The transcriptional profiles were similar in the two species, including the largely predominant expression of Hd27-vank1 in all tissues examined. However, in various insect cell lines, the expression patterns of HdIV vankyrins differed according to species. No clear relationship between vankyrin expression patterns and abundance of vankyrin-bearing genomic segments were found in the lepidopteran cell lines. Moreover, in these cells, the amount of vankyrin-bearing genomic segments differed substantially between cytosol and nuclei of infected cells, implying the existence of an unexpected step regulating the copy number of HdIV segments in cell nuclei. Our in vitro results reveal a host-specific transcriptional profile of vankyrins that may be related to the success of parasitism in different hosts. In Spodoptera hosts, the predominant expression of Hd27-vank1 suggests that this protein might have pleiotropic functions during parasitism of these insect species.

Lepidopteran transcriptome analysis following infection by phylogenetically unrelated polydnaviruses highlights differential and common responses

The Polydnaviridae is a family of double-stranded DNA viruses that are symbionts of parasitoid wasps. The family is currently divided into two genera, the Ichnovirus (IV) and Bracovirus (BV), which are associated with wasps in the families Ichneumonidae and Braconidae, respectively. IVs and BVs have similar immunosuppressive and developmental effects on parasitized hosts but their encapsidated genomes largely encode different genes. To assess whether IV and BV infection has similar or disparate effects on the transcriptome of shared hosts, we characterized the effects of Hyposoter didymator Ichnovirus (HdIV) and Microplitis demolitor Bracovirus (MdBV) on the fat body and hemocyte transcriptome of Spodoptera frugiperda (Lepidoptera: Noctuidae). Our results indicated that HdIV and MdBV infection alters the abundance of a relatively low proportion of S. frugiperda transcripts at 24 h post-infection. A majority of the transcripts affected by infection also differed between MdBV and HdIV. However, we did identify some host transcripts that were similarly affected by both viruses. A majority of these genes were transcribed in the fat body and most belonged to functional classes with roles in immunity, detoxification, or cell structure. Particularly prominent in this suite of transcripts were genes encoding for predicted motor-related and collagen IV-like proteins. Overall, our data suggest that the broadly similar effects that HdIV and MdBV have on host growth and immunity are not due to these viruses inducing profound changes in host gene expression. Given though that IVs and BVs encode few shared genes, the host transcripts that are similarly affected by HdIV and MdBV could indicate convergence by each virus to target a few processes at the level of transcription that are important for successful parasitism of hosts by H. didymator and M. demolitor.

Analysis of virion structural components reveals vestiges of the ancestral ichnovirus genome

Many thousands of endoparasitic wasp species are known to inject polydnavirus (PDV) particles into their caterpillar host during oviposition, causing immune and developmental dysfunctions that benefit the wasp larva. PDVs associated with braconid and ichneumonid wasps, bracoviruses and ichnoviruses respectively, both deliver multiple circular dsDNA molecules to the caterpillar. These molecules contain virulence genes but lack core genes typically involved in particle production. This is not completely unexpected given that no PDV replication takes place in the caterpillar. Particle production is confined to the wasp ovary where viral DNAs are generated from proviral copies maintained within the wasp genome. We recently showed that the genes involved in bracovirus particle production reside within the wasp genome and are related to nudiviruses. In the present work we characterized genes involved in ichnovirus particle production by analyzing the components of purified Hyposoter didymator Ichnovirus particles by LC-MS/MS and studying their organization in the wasp genome. Their products are conserved among ichnovirus-associated wasps and constitute a specific set of proteins in the virosphere. Strikingly, these genes are clustered in specialized regions of the wasp genome which are amplified along with proviral DNA during virus particle replication, but are not packaged in the particles. Clearly our results show that ichnoviruses and bracoviruses particles originated from different viral entities, thus providing an example of convergent evolution where two groups of wasps have independently domesticated viruses to deliver genes into their hosts.

Inhibition of in vitro angiogenesis by platelet factor-4-derived peptides and mechanism of action

In this study, we examined in detail the interaction of platelet factor-4 (PF-4) with fibroblast growth factor-2 (FGF-2) and vascular endothelial growth factor (VEGF) and the effect of PF-4-derived synthetic peptides. We show that a peptide between amino acids 47 and 70 that contains the heparin-binding lysine-rich site inhibits FGF-2 or VEGF function. This is based on the following observations: PF-4 peptide 47-70 inhibited FGF-2 or VEGF binding to endothelial cells; it inhibited FGF-2 or VEGF binding to FGFRs or VEGFRs in heparan sulfate-deficient CHO cells transfected with FGFR1 (CHOFGFR1) or VEGFR2 (CHOmVEGFR2) cDNA; it blocked proliferation or tube formation in three-dimensional angiogenesis assays; and, finally, it competed with the direct association of (125)I-PF-4 with FGF-2 or VEGF, respectively, and inhibited heparin-induced FGF-2 dimerization. A shorter C-terminal peptide (peptide 58-70), which still contained the heparin-binding lysin-rich site, had no effect. Peptide 17-58, which is located in the central part of the molecule, although it does not inhibit FGF-2 or VEGF binding or biologic activity in endothelial cells, inhibited heparin-dependent binding of (125)I-FGF-2 or (125)I-VEGF to CHOmFGFR1 or CHOmVEGFR2 cells, respectively. Shorter peptides (peptides 34-58 and 47-58) did not show any of these effects.