Unfolding the evolutionary story of polydnaviruses

Cellular dynamics of host - parasitoid interactions: Insights from the encapsulation process in a partially resistant host

Cotesia typhae is an eastern African endoparasitoid braconid wasp that targets the larval stage of the lepidopteran stem borer, Sesamia nonagrioides, a maize crop pest in Europe. The French host population is partially resistant to the Makindu strain of the wasp, allowing its development in only 40% of the cases. Resistant larvae can encapsulate the parasitoid and survive the infection. This interaction provides a very interesting frame for investigating the impact of parasitism on host cellular resistance. We characterized the parasitoid ovolarval development in a permissive host and studied the encapsulation process in a resistant host by dissection and histological sectioning compared to that of inert chromatography beads. We measured the total hemocyte count in parasitized and bead-injected larvae over time to monitor the magnitude of the immune reaction. Our results show that parasitism of resistant hosts delayed encapsulation but did not affect immune abilities towards inert beads. Moreover, while bead injection increased total hemocyte count, it remained constant in resistant and permissive larvae. We conclude that while Cotesia spp virulence factors are known to impair the host immune system, our results suggest that passive evasion could also occur.

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.

Insights into the venom protein components of the egg parasitoid Anastatus japonicus (Hymenoptera: Eupelmidae)

BACKGROUND

Parasitoid venom is composed of a complex mixture of various active substances with different biological functions and is injected in the host during the parasitoid oviposition. Anastatus japonicus (Hymenoptera: Eupelmidae) is an egg parasite of Tessaratoma papillosa (Hemiptera: Tessaratomidae). Although the venom of this egg parasitoid plays an important role in the parasitic process, relatively little work has been done to address the mechanism.

RESULTS

In the present study, proteomic analysis was performed to identify the proteins that play an important role in the parasitic process of A. japonicus. A total of 2084 proteins were identified, including 81 putative venom proteins, most of which were identified as Hexamerin, Chitinase 2, Calreticulin, Heat shock protein 83-like, Serine protease, Arginine kinase, Phosphoserine aminotransferase and Actin protein. Together the before (Be) and after (Af) parasitization venom contains 1628 proteins, including 212 DEPs with 181 and 31 significantly up-regulated and down-regulated respectively. In addition, 10 differentially expressed proteins (DEPs) with fold change ≥8.71 were subjected to RT-qPCR to validate the proteomic data. The differential expression analysis revealed that nine proteins were specifically present in the pre-parasitic venom, whereas 26 proteins were specific to the post-parasitic treatments. Results of RT-qPCR analysis showed high expression of the selected DEPs which further validated our proteomics data.

CONCLUSION

These new proteomic data greatly enrich our current knowledge about key venom proteins associated with parasitic process in A. japonicus and contribute to better understanding of the parasitic mechanisms leading to the development of new biological control strategies. © 2020 Society of Chemical Industry.

The Complete Genome of an Endogenous Nimavirus (Nimav-1_LVa) From the Pacific Whiteleg Shrimp Penaeus (Litopenaeus) Vannamei

White spot syndrome virus (WSSV), the lone virus of the genus Whispovirus under the family Nimaviridae, is one of the most devastating viruses affecting the shrimp farming industry. Knowledge about this virus, in particular, its evolution history, has been limited, partly due to its large genome and the lack of other closely related free-living viruses for comparative studies. In this study, we reconstructed a full-length endogenous nimavirus consensus genome, Nimav-1_LVa (279,905 bp), in the genome sequence of Penaeus (Litopenaeus) vannamei breed Kehai No. 1 (ASM378908v1). This endogenous virus seemed to insert exclusively into the telomeric pentanucleotide microsatellite (TAACC/GGTTA)n. It encoded 117 putative genes, with some containing introns, such as g012 (inhibitor of apoptosis, IAP), g046 (crustacean hyperglycemic hormone, CHH), g155 (innexin), g158 (Bax inhibitor 1 like). More than a dozen Nimav-1_LVa genes are involved in the pathogen-host interactions. We hypothesized that g046, g155, g158, and g227 (semaphorin 1A like) were recruited host genes for their roles in immune regulation. Sequence analysis indicated that a total of 43 WSSV genes belonged to the ancestral/core nimavirus gene set, including four genes reported in this study: wsv112 (dUTPase), wsv206, wsv226, and wsv308 (nucleocapsid protein). The availability of the Nimav-1_LVa sequence would help understand the genetic diversity, epidemiology, evolution, and virulence of WSSV.

Effects of ovarian fluid, venom and egg surface characteristics of Tetrastichus brontispae (Hymenoptera: Eulophidae) on the immune response of Octodonta nipae (Coleoptera: Chrysomelidae)

Although the importance of parasitoids as biocontrol agents has long been recognized, systematic studies of the physiological mechanisms are scarce, especially in those parasitoids that are able to successfully invade their hosts by activating host immune responses. This study explored this phenomenon by investigating the effects of ovarian fluid, venom and egg surface characteristics of Tetrastichus brontispae (Hymenoptera: Eulophidae) on host immunity. The results showed that the injection of venom alone induced higher phenoloxidase activity, while a mixture of ovarian plus venom fluids provoked higher granulocyte and plasmatocyte spreading ratios, highlighting the role that egg surface characteristics may play in successful parasitism. After thorough investigation, the presence of a hemomucin homologue was documented on the egg surface (which was named Tetrastichus brontispae adipocyte plasma membrane associated protein-like, TbAPMAP-like), while the absence of polydnaviruses, fibrous layers and virus-like filaments was confirmed. The higher encapsulation index of eggs incubated with TbAPMAP-like polyclonal antibody demonstrated the protection of the protein against encapsulation. These results contribute to our understanding of the mechanisms used by endoparasitoids to evade encapsulation during the early parasitism stage while enriching our knowledge of local active regulatory mechanisms. It is likely that this is the first study to determine the egg protective properties of TbAPMAP-like in host-parasite systems.

Baculoviral p94 homologs encoded in Cotesia plutellae bracovirus suppress both immunity and development of the diamondback moth, Plutellae xylostella

Polydnaviruses (PDVs) are a group of insect DNA viruses, which exhibit a mutual symbiotic relationship with their specific host wasps. Moreover, most encapsidated genes identified so far in PDVs share homologies with insect-originated genes, but not with virus-originated genes. In the meantime, PDVs associated with 2 wasp genera Cotesia and Glytapanteles encode some genes presumably originated from other viruses. Cotesia plutellae bracovirus (CpBV) encodes 4 genes homologous to baculoviral p94: CpBV-E94k1, CpBV-E94k2, CpBV-E94k3, and CpBV-E94k4. This study was conducted to predict the origin of CpBV-E94ks by comparing their sequences with those of baculoviral orthologs and to determine the physiological functions by their transient expressions in nonparasitized larvae and subsequent specific RNA interference. Our phylogenetic analysis indicated that CpBV-E94ks were clustered with other E94ks originated from different PDVs and shared high similarity with betabaculoviral p94s. These 4 CpBV genes were expressed during most developmental stages of the larvae of Plutella xylostella parasitized by C. plutellae. Expression of these 4 E94ks was mainly detected in hemocytes and fat body. Subsequent functional analysis by in vivo transient expression showed that all 4 viral genes significantly inhibited both host immune and developmental processes. These results suggest that CpBV-E94ks share an origin with betabaculoviral p94s and play parasitic roles in suppressing host immune and developmental processes.

Polydnaviruses: Nature's Genetic Engineers

Virus-host associations are usually viewed as parasitic, but several studies in recent years have reported examples of viruses that benefit host organisms. The Polydnaviridae are of particular interest because these viruses are all obligate mutualists of insects called parasitoid wasps. Parasitoids develop during their immature stages by feeding inside the body of other insects, which serve as their hosts. Polydnaviruses are vertically transmitted as proviruses through the germ line of wasps but also function as gene delivery vectors that wasps rely upon to genetically manipulate the hosts they parasitize. Here we review the evolutionary origin of polydnaviruses, the organization and function of their genomes, and some of their roles in parasitism.

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

Bracoviruses (BVs) from the Polydnaviridae family are symbiotic viruses used as biological weapons by parasitoid wasps to manipulate lepidopteran host physiology and induce parasitism success. BV particles are produced by wasp ovaries and injected along with the eggs into the caterpillar host body, where viral gene expression is necessary for wasp development. Recent sequencing of the proviral genome of Cotesia congregata BV (CcBV) identified 222 predicted virulence genes present on 35 proviral segments integrated into the wasp genome. To date, the expressions of only a few selected candidate virulence genes have been studied in the caterpillar host, and we lacked a global vision of viral gene expression. In this study, a large-scale transcriptomic analysis by 454 sequencing of two immune tissues (fat body and hemocytes) of parasitized Manduca sexta caterpillar hosts allowed the detection of expression of 88 CcBV genes expressed 24 h after the onset of parasitism. We linked the expression profiles of these genes to several factors, showing that different regulatory mechanisms control viral gene expression in the host. These factors include the presence of signal peptides in encoded proteins, diversification of promoter regions, and, more surprisingly, gene position on the proviral genome. Indeed, most genes for which expression could be detected are localized in particular proviral regions globally producing higher numbers of circles. Moreover, this polydnavirus (PDV) transcriptomic analysis also reveals that a majority of CcBV genes possess at least one intron and an arthropod transcription start site, consistent with an insect origin of these virulence genes.

IMPORTANCE

Bracoviruses (BVs) are symbiotic polydnaviruses used by parasitoid wasps to manipulate lepidopteran host physiology, ensuring wasp offspring survival. To date, the expressions of only a few selected candidate BV virulence genes have been studied in caterpillar hosts. We performed a large-scale analysis of BV gene expression in two immune tissues of Manduca sexta caterpillars parasitized by Cotesia congregata wasps. Genes for which expression could be detected corresponded to genes localized in particular regions of the viral genome globally producing higher numbers of circles. Our study thus brings an original global vision of viral gene expression and paves the way to the determination of the regulatory mechanisms enabling the expression of BV genes in targeted organisms, such as major insect pests. In addition, we identify sequence features suggesting that most BV virulence genes were acquired from insect genomes.

Functional endogenous viral elements in the genome of the parasitoid wasp Cotesia congregata: insights into the evolutionary dynamics of bracoviruses

Bracoviruses represent the most complex endogenous viral elements (EVEs) described to date. Nudiviral genes have been hosted within parasitoid wasp genomes since approximately 100 Ma. They play a crucial role in the wasp life cycle as they produce bracovirus particles, which are injected into parasitized lepidopteran hosts during wasp oviposition. Bracovirus particles encapsidate multiple dsDNA circles encoding virulence genes. Their expression in parasitized caterpillars is essential for wasp parasitism success. Here, we report on the genomic organization of the proviral segments (i.e. master sequences used to produce the encapsidated dsDNA circles) present in the Cotesia congregata parasitoid wasp genome. The provirus is composed of a macrolocus, comprising two-thirds of the proviral segments and of seven dispersed loci, each containing one to three segments. Comparative genomic analyses with closely related species gave insights into the evolutionary dynamics of bracovirus genomes. Conserved synteny in the different wasp genomes showed the orthology of the proviral macrolocus across different species. The nudiviral gene odv-e66-like1 is conserved within the macrolocus, suggesting an ancient co-localization of the nudiviral genome and bracovirus proviral segments. By contrast, the evolution of proviral segments within the macrolocus has involved a series of lineage-specific duplications.

The bracovirus genome of the parasitoid wasp Cotesia congregata is amplified within 13 replication units, including sequences not packaged in the particles

The relationship between parasitoid wasps and polydnaviruses constitutes one of the few known mutualisms between viruses and eukaryotes. Viral particles are injected with the wasp eggs into parasitized larvae, and the viral genes thus introduced are used to manipulate lepidopteran host physiology. The genome packaged in the particles is composed of 35 double-stranded DNA (dsDNA) circles produced in wasp ovaries by amplification of viral sequences from proviral segments integrated in tandem arrays in the wasp genome. These segments and their flanking regions within the genome of the wasp Cotesia congregata were recently isolated, allowing extensive mapping of amplified sequences. The bracovirus DNAs packaged in the particles were found to be amplified within more than 12 replication units. Strikingly, the nudiviral cluster, the genes of which encode particle structural components, was also amplified, although not encapsidated. Amplification of bracoviral sequences was shown to involve successive head-to-head and tail-to-tail concatemers, which was not expected given the nudiviral origin of bracoviruses.

Adaptive selection on bracovirus genomes drives the specialization of Cotesia parasitoid wasps

The geographic mosaic of coevolution predicts parasite virulence should be locally adapted to the host community. Cotesia parasitoid wasps adapt to local lepidopteran species possibly through their symbiotic bracovirus. The virus, essential for the parasitism success, is at the heart of the complex coevolutionary relationship linking the wasps and their hosts. The large segmented genome contained in the virus particles encodes virulence genes involved in host immune and developmental suppression. Coevolutionary arms race should result in the positive selection of particular beneficial alleles. To understand the global role of bracoviruses in the local adaptation or specialization of parasitoid wasps to their hosts, we studied the molecular evolution of four bracoviruses associated with wasps of the genus Cotesia, including C congregata, C vestalis and new data and annotation on two ecologically differentiated populations of C sesamie, Kitale and Mombasa. Paired orthologs analyses revealed more genes under positive selection when comparing the two C sesamiae bracoviruses belonging to the same species, and more genes under strong evolutionary constraint between species. Furthermore branch-site evolutionary models showed that 17 genes, out of the 54 currently available shared by the four bracoviruses, harboured sites under positive selection including: the histone H4-like, a C-type lectin, two ep1-like, ep2, a viral ankyrin, CrV1, a ben-domain, a Serine-rich, and eight unknown genes. Lastly the phylogenetic analyses of the histone, ep2 and CrV1 genes in different African C sesamiae populations showed that each gene described differently the individual relationships. In particular we found recombination had happened between the ep2 and CrV1 genes, which are localized 37.5 kb apart on the wasp chromosomes. Involved in multidirectional coevolutionary interactions, C sesamiae wasps rely on different bracovirus mediated molecular pathways to overcome local host resistance.

Characterization of the polydnaviral ‘T. rostrale virus’ (TrV) gene family: TrV1 expression inhibits in vitro cell proliferation

Tranosema rostrale ichnovirus (TrIV) is a polydnavirus (PDV) transmitted by the endoparasitic wasp T. rostrale to its host Choristoneura fumiferana during oviposition. PDV genes are expressed in infected caterpillars, causing physiological disturbances that promote the survival of the developing endoparasite. The previously sequenced genome of TrIV contains ~86 genes organized in multigene families and distributed on multiple segments of circular dsDNA. Among these, the ‘T. rostrale virus’ (TrV) family comprises seven genes that are absent in other PDV genomes examined to date and whose function(s) remain(s) unknown. Here, we initiated a functional analysis of the TrV family using qPCR, transfection and RNAi approaches. TrV family genes were weakly expressed in wasp ovaries, but some displayed high transcript abundance in parasitized caterpillars. Whilst TrV1 was the most highly transcribed TrV gene in infected caterpillars, transcript levels for TrV5 and TrV6 were nearly undetectable, indicating that they may be pseudogenes. Temporal and tissue-specific patterns of transcript abundance were similar for all expressed TrV family genes, indicative of an apparent lack of difference in function or tissue specificity. Infection of Cf-203 and Sf-21 insect cells with TrIV led to a dose-dependent inhibition of cell proliferation with no sign of apoptosis. Whilst similar inhibition was observed following transfection of cells with a cloned genome segment carrying the TrV1 gene, RNA interference targeting TrV1 largely restored cell growth in TrIV-infected cells, indicating that TrV1 expression was responsible for the observed inhibition. We suggest that TrV genes may contribute to host developmental disruption by interfering with host-cell proliferation during parasitism.

Conservation of a triple-helix-forming RNA stability element in noncoding and genomic RNAs of diverse viruses

Abundant expression of the long noncoding (lnc) PAN (polyadenylated nuclear) RNA by the human oncogenic gammaherpesvirus Kaposi's sarcoma-associated herpesvirus (KSHV) depends on a cis-element called the expression and nuclear retention element (ENE). The ENE upregulates PAN RNA by inhibiting its rapid nuclear decay through triple-helix formation with the poly(A) tail. Using structure-based bioinformatics, we identified six ENE-like elements in evolutionarily diverse viral genomes. Five are in double-stranded DNA viruses, including mammalian herpesviruses, insect polydnaviruses, and a protist mimivirus. One is in an insect picorna-like positive-strand RNA virus, suggesting that the ENE can counteract cytoplasmic as well as nuclear RNA decay pathways. Functionality of four of the ENEs was demonstrated by increased accumulation of an intronless polyadenylated reporter transcript in human cells. Identification of these ENEs enabled the discovery of PAN RNA homologs in two additional gammaherpesviruses, RRV and EHV2. Our findings demonstrate that searching for structural elements can lead to rapid identification of lncRNAs.

Differential inhibition of BmRelish1-dependent transcription in lepidopteran cells by bracovirus ankyrin-repeat proteins

In the tripartite parasitization system of the lepidopteran host Manduca sexta, the endoparasitoid wasp Cotesia congregata and its endosymbiotic virus, C. congregata Bracovirus (CcBV), the expression of viral proteins is necessary for successful parasitization. Here we have examined the in vitro effects of six members of the ankyrin-repeat protein family (Ank) of CcBV, which are thought to interfere with the host's induced innate immune responses, on the transcriptional activity of a heterologous lepidopteran Rel/NFκB transcription factor, Relish1 of Bombyx mori. Using as transcriptional activator BmRelish1-d2 (R1d2), a constitutively active mutant of the major regulator of the Imd pathway, BmRelish1, in conjunction with a reporter gene controlled by a B. mori antimicrobial peptide gene promoter, we have found that 5 of the 6 examined Anks suppress R1d2-dependent transcriptional activity to various degrees. Immunofluorescence studies have also revealed that while some of the Ank proteins have a rather strict cytoplasmic localization, others are detected both in the cytoplasm and the nucleus of the expressing cells and that colocalization with R1d2 occurs exclusively in the nucleus. Thus, our results suggest that functional and spatial differences among the various CcBV Ank family members may be responsible for the observed differential inhibition of R1d2 activity.

Paleozoic origin of insect large dsDNA viruses

To understand how extant viruses interact with their hosts, we need a historical framework of their evolutionary association. Akin to retrovirus or hepadnavirus viral fossils present in eukaryotic genomes, bracoviruses are integrated in braconid wasp genomes and are transmitted by Mendelian inheritance. However, unlike viral genomic fossils, they have retained functional machineries homologous to those of large dsDNA viruses pathogenic to arthropods. Using a phylogenomic approach, we resolved the relationships between bracoviruses and their closest free relatives: baculoviruses and nudiviruses. The phylogeny showed that bracoviruses are nested within the nudivirus clade. Bracoviruses establish a bridge between the virus and animal worlds. Their inclusion in a virus phylogeny allowed us to relate free viruses to fossils. The ages of the wasps were used to calibrate the virus phylogeny. Bayesian analyses revealed that insect dsDNA viruses first evolved at ∼310 Mya in the Paleozoic Era during the Carboniferous Period with the first insects. Furthermore the virus diversification time frame during the Mesozoic Era appears linked to the diversification of insect orders; baculoviruses that infect larvae evolved at the same period as holometabolous insects. These results imply ancient coevolution by resource tracking between several insect dsDNA virus families and their hosts, dating back to 310 Mya.

Transcriptomic profiling of Microplitis demolitor bracovirus reveals host, tissue and stage-specific patterns of activity

The polydnaviruses (PDVs) are a family of DNA viruses that are symbiotically associated with parasitoid wasps. The transcription of particular genes or gene-family members have been reported for several PDVs, but no studies have characterized the spatio-temporal patterns of expression for the entire complement of predicted genes in the encapsidated genome of any PDV isolate. The braconid wasp Microplitis demolitor carries the PDV Microplitis demolitor bracovirus (MdBV) and parasitizes larval stage Pseudoplusia (Chrysodeixis) includens. The encapsidated genome consists of 15 genomic segments with 51 predicted ORFs encoding proteins ≥100 aa. A majority of these ORFs form four multimember gene families (ptp, ank, glc and egf) while the remaining ORFs consist of single copy (orph) genes. Here we used RT-PCR and quantitative real-time PCR methods to profile the encapsidated transcriptome of MdBV in P. includens and M. demolitor. Our results indicate that most predicted genes are expressed in P. includens. Spatial patterns of expression in P. includens differed among genes, but temporal patterns of expression were generally similar, with transcript abundance progressively declining between 24 and 120 h. A subset of ptp, ank and orph genes were also expressed in adult female but not male M. demolitor. Only one encapsidated gene (ank-H4) was expressed in all life stages of M. demolitor, albeit at much lower levels than in P. includens. However, another encapsidated gene (orph-B1) was expressed in adult M. demolitor at similar levels to those detected in P. includens.

The encapsidated genome of Microplitis demolitor bracovirus integrates into the host Pseudoplusia includens

Polydnaviruses (PDVs) are symbionts of parasitoid wasps that function as gene delivery vehicles in the insects (hosts) that the wasps parasitize. PDVs persist in wasps as integrated proviruses but are packaged as circularized and segmented double-stranded DNAs into the virions that wasps inject into hosts. In contrast, little is known about how PDV genomic DNAs persist in host cells. Microplitis demolitor carries Microplitis demolitor bracovirus (MdBV) and parasitizes the host Pseudoplusia includens. MdBV infects primarily host hemocytes and also infects a hemocyte-derived cell line from P. includens called CiE1 cells. Here we report that all 15 genomic segments of the MdBV encapsidated genome exhibited long-term persistence in CiE1 cells. Most MdBV genes expressed in hemocytes were persistently expressed in CiE1 cells, including members of the glc gene family whose products transformed CiE1 cells into a suspension culture. PCR-based integration assays combined with cloning and sequencing of host-virus junctions confirmed that genomic segments J and C persisted in CiE1 cells by integration. These genomic DNAs also rapidly integrated into parasitized P. includens. Sequence analysis of wasp-viral junction clones showed that the integration of proviral segments in M. demolitor was associated with a wasp excision/integration motif (WIM) known from other bracoviruses. However, integration into host cells occurred in association with a previously unknown domain that we named the host integration motif (HIM). The presence of HIMs in most MdBV genomic DNAs suggests that the integration of each genomic segment into host cells occurs through a shared mechanism.

Functional interactions between polydnavirus and host cellular innexins

Polydnaviruses are double-stranded DNA viruses associated with some subfamilies of ichneumonoid parasitoid wasps. Polydnavirus virions are delivered during wasp parasitization of a host, and virus gene expression in the host induces alterations of host physiology. Infection of susceptible host caterpillars by the polydnavirus Campoletis sonorensis ichnovirus (CsIV) leads to expression of virus genes, resulting in immune and developmental disruptions. CsIV carries four homologues of insect gap junction genes (innexins) termed vinnexins, which are expressed in multiple tissues of infected caterpillars. Previously, we demonstrated that two of these, VinnexinD and VinnexinG, form functional gap junctions in paired Xenopus oocytes. Here we show that VinnexinQ1 and VinnexinQ2, likewise, form junctions in this heterologous system. Moreover, we demonstrate that the vinnexins interact differentially with the Innexin2 orthologue of an ichnovirus host, Spodoptera frugiperda. Cell pairs coexpressing a vinnexin and Innexin2 or pairs in which one cell expresses a vinnexin and the neighboring cell Innexin2 assemble functional junctions with properties that differ from those of junctions composed of Innexin2 alone. These data suggest that altered gap junctional intercellular communication may underlie certain cellular pathologies associated with ichnovirus infection of caterpillar hosts.

Pervasive horizontal transfer of rolling-circle transposons among animals

Horizontal transfer (HT) of genes is known to be an important mechanism of genetic innovation, especially in prokaryotes. The impact of HT of transposable elements (TEs), however, has only recently begun to receive widespread attention and may be significant due to their mutagenic potential, inherent mobility, and abundance. Helitrons, also known as rolling-circle transposons, are a distinctive subclass of TE with a unique transposition mechanism. Here, we describe the first evidence for the repeated HT of four different families of Helitrons in an unprecedented array of organisms, including mammals, reptiles, fish, invertebrates, and insect viruses. The Helitrons present in these species have a patchy distribution and are closely related (80–98% sequence identity), despite the deep divergence times among hosts. Multiple lines of evidence indicate the extreme conservation of sequence identity is not due to selection, including the highly fragmented nature of the Helitrons identified and the lack of any signatures of selection at the nucleotide level. The presence of horizontally transferred Helitrons in insect viruses, in particular, suggests that this may represent a potential mechanism of transfer in some taxa. Unlike genes, Helitrons that have horizontally transferred into new host genomes can amplify, in some cases reaching up to several hundred copies and representing a substantial fraction of the genome. Because Helitrons are known to frequently capture and amplify gene fragments, HT of this unique group of DNA transposons could lead to horizontal gene transfer and incur dramatic shifts in the trajectory of genome evolution.

New dimensions of the virus world discovered through metagenomics

Metagenomic analysis of viruses suggests novel patterns of evolution, changes the existing ideas of the composition of the virus world and reveals novel groups of viruses and virus-like agents. The gene composition of the marine DNA virome is dramatically different from that of known bacteriophages. The virome is dominated by rare genes, many of which might be contained within virus-like entities such as gene transfer agents. Analysis of marine metagenomes thought to consist mostly of bacterial genes revealed a variety of sequences homologous to conserved genes of eukaryotic nucleocytoplasmic large DNA viruses, resulting in the discovery of diverse members of previously undersampled groups and suggesting the existence of new classes of virus-like agents. Unexpectedly, metagenomics of marine RNA viruses showed that representatives of only one superfamily of eukaryotic viruses, the picorna-like viruses, dominate the RNA virome.

The natural genetic engineering of polydnaviruses

The polydnaviruses represent an unusual example of a highly evolved symbiosis between some parasitic wasps, DNA containing particles or viruses, and lepidopteran larval hosts of the wasp. The viruses can no longer replicate independently, as genes that encode viral structural proteins are restricted to the wasp genome and are not encapsidated. Interestingly, the DNA that is encapsidated is more similar in terms of gene identity and gene density to eukaryotic genomes than viral genomes. We compare and relate this unusual example of natural genetic engineering to the well-known system of viral lysogeny. The similarities in the two systems may prove useful in understanding the replication strategy and genomic organization of polydnaviruses and provide some insight into how this unusual virus system may have evolved.

Polydnavirus hidden face: the genes producing virus particles of parasitic wasps

Very few obligatory relationships involve viruses to the remarkable exception of polydnaviruses (PDVs) associated with tens of thousands species of parasitic wasps that develop within the body of lepidopteran larvae. PDV particles, injected along with parasite eggs into the host body, act by manipulating host immune defences, development and physiology, thereby enabling wasp larvae to survive in a potentially harmful environment. Particle production does not occur in infected tissues of parasitized caterpillars, but is restricted to specialized cells of the wasp ovaries. Moreover, the genome enclosed in the particles encodes almost no viral structural protein, but mostly factors used to manipulate the physiology of the parasitized host. We recently unravelled the viral nature of PDVs associated with braconid wasps by characterizing a large set of nudivirus genes residing permanently in the wasp chromosome(s). Many of these genes encode structural components of the bracovirus particles and their expression pattern correlates with particle production. They constitute a viral machinery comprising a large number of core genes shared by nudiviruses and baculoviruses. Thus bracoviruses do not appear to be nudiviruses remnants, but instead complex nudiviral devices carrying DNA for the delivery of virulence genes into lepidopteran hosts. This highlights the fact that viruses should no longer be exclusively considered obligatory parasites, and that in certain cases they are obligatory symbionts.

N-terminal tail of a viral histone H4 encoded in Cotesia plutellae bracovirus is essential to suppress gene expression of host histone H4

An endoparasitoid wasp, Cotesia plutellae, possesses a symbiotic bracovirus (CpBV), which facilitates parasitism of a specific host, such as larvae of the diamondback moth, Plutella xylostella. A viral histone H4 (CpBV-H4) has been found in the CpBV genome and its gene product plays a role in impairing the host insect cellular immune response. Based on its high similarity to histone H4 of P. xylostella apart from its extended N-terminal tail, it has been suspected to alter host gene expression. Histone subunits were purified from parasitized P. xylostella larvae and found to contain both host and viral H4s, confirming a previous report of a possible epigenetic mode of action. Moreover, this study showed that the host H4 levels in the parasitized larvae clearly decreased during the parasitization period, whereas CpBV-H4 levels maintained a significant level without significant changes. To understand the decrease of host H4 levels, transcription levels of host H4 were monitored by quantitative reverse-transcriptase PCR (RT-PCR) and showed a significant decrease in parasitized P. xylostella larvae, whereas no significant change of the mRNA level was detected in nonparasitized larvae. This transcriptional control of host H4 expression was also observed by inducing transient expression of CpBV-H4 in nonparasitized P. xylostella. Moreover, co-injection of CpBV-H4 and its specific double-stranded RNA recovered the host H4 expression level. To identify a functional domain of CpBV-H4 involved in the transcriptional control, the extended N-terminal tail of CpBV-H4 was removed by preparing a truncated viral H4 construct in an expression vector by deleting the N-terminal tail of 38 amino acid residues and inducing its expression in nonparasitized P. xylostella larvae. The truncated CpBV-H4 clearly lost its inhibitory effects on host H4 transcription. Moreover, the presence of CpBV-H4 affects the spreading of host haemocytes by an epigenetic effect, which is at least partly restored in larvae expressing the truncated version of CpBV-H4. This study suggests that the viral H4 encoded in CpBV can alter host gene expression with its extended N-terminal tail.

Comparative genomics of mutualistic viruses of Glyptapanteles parasitic wasps

Comparative genome analysis of two endosymbiotic polydnaviruses from Glyptapanteles parasitic wasps reveals new insights into the evolutionary arms race between host and parasite.

Background

Polydnaviruses, double-stranded DNA viruses with segmented genomes, have evolved as obligate endosymbionts of parasitoid wasps. Virus particles are replication deficient and produced by female wasps from proviral sequences integrated into the wasp genome. These particles are co-injected with eggs into caterpillar hosts, where viral gene expression facilitates parasitoid survival and, thereby, survival of proviral DNA. Here we characterize and compare the encapsidated viral genome sequences of bracoviruses in the family Polydnaviridae associated with Glyptapanteles gypsy moth parasitoids, along with near complete proviral sequences from which both viral genomes are derived.

Results

The encapsidated Glyptapanteles indiensis and Glyptapanteles flavicoxis bracoviral genomes, each composed of 29 different size segments, total approximately 517 and 594 kbp, respectively. They are generated from a minimum of seven distinct loci in the wasp genome. Annotation of these sequences revealed numerous novel features for polydnaviruses, including insect-like sugar transporter genes and transposable elements. Evolutionary analyses suggest that positive selection is widespread among bracoviral genes.

Conclusions

The structure and organization of G. indiensis and G. flavicoxis bracovirus proviral segments as multiple loci containing one to many viral segments, flanked and separated by wasp gene-encoding DNA, is confirmed. Rapid evolution of bracovirus genes supports the hypothesis of bracovirus genes in an 'arms race' between bracovirus and caterpillar. Phylogenetic analyses of the bracoviral genes encoding sugar transporters provides the first robust evidence of a wasp origin for some polydnavirus genes. We hypothesize transposable elements, such as those described here, could facilitate transfer of genes between proviral segments and host DNA.

The Big Bang of picorna-like virus evolution antedates the radiation of eukaryotic supergroups

The recent discovery of RNA viruses in diverse unicellular eukaryotes and developments in evolutionary genomics have provided the means for addressing the origin of eukaryotic RNA viruses. The phylogenetic analyses of RNA polymerases and helicases presented in this Analysis article reveal close evolutionary relationships between RNA viruses infecting hosts from the Chromalveolate and Excavate supergroups and distinct families of picorna-like viruses of plants and animals. Thus, diversification of picorna-like viruses probably occurred in a 'Big Bang' concomitant with key events of eukaryogenesis. The origins of the conserved genes of picorna-like viruses are traced to likely ancestors including bacterial group II retroelements, the family of HtrA proteases and DNA bacteriophages.

Molecular evidence for the evolution of ichnoviruses from ascoviruses by symbiogenesis

BACKGROUND

Female endoparasitic ichneumonid wasps inject virus-like particles into their caterpillar hosts to suppress immunity. These particles are classified as ichnovirus virions and resemble ascovirus virions, which are also transmitted by parasitic wasps and attack caterpillars. Ascoviruses replicate DNA and produce virions. Polydnavirus DNA consists of wasp DNA replicated by the wasp from its genome, which also directs particle synthesis. Structural similarities between ascovirus and ichnovirus particles and the biology of their transmission suggest that ichnoviruses evolved from ascoviruses, although molecular evidence for this hypothesis is lacking.

RESULTS

Here we show that a family of unique pox-D5 NTPase proteins in the Glypta fumiferanae ichnovirus are related to three Diadromus pulchellus ascovirus proteins encoded by ORFs 90, 91 and 93. A new alignment technique also shows that two proteins from a related ichnovirus are orthologs of other ascovirus virion proteins.

CONCLUSION

Our results provide molecular evidence supporting the origin of ichnoviruses from ascoviruses by lateral transfer of ascoviral genes into ichneumonid wasp genomes, perhaps the first example of symbiogenesis between large DNA viruses and eukaryotic organisms. We also discuss the limits of this evidence through complementary studies, which revealed that passive lateral transfer of viral genes among polydnaviral, bacterial, and wasp genomes may have occurred repeatedly through an intimate coupling of both recombination and replication of viral genomes during evolution. The impact of passive lateral transfers on evolutionary relationships between polydnaviruses and viruses with large double-stranded genomes is considered in the context of the theory of symbiogenesis.

Protein tyrosine phosphatase-H2 from a polydnavirus induces apoptosis of insect cells

The family Polydnaviridae is a large group of immunosuppressive insect viruses that are symbiotically associated with parasitoid wasps. The polydnavirus Microplitis demolitor bracovirus (MdBV) causes several alterations that disable the cellular and humoral immune defences of host insects, including apoptosis of the primary phagocytic population of circulating immune cells (haemocytes), called granulocytes. Here, we show that MdBV infection causes granulocytes in the lepidopteran Spodoptera frugiperda to apoptose. An expression screen conducted in the S. frugiperda 21 cell line identified the MdBV gene ptp-H2 as an apoptosis inducer, as indicated by cell fragmentation, annexin V binding, mitochondrial membrane depolarization and caspase activation. PTP-H2 is a classical protein tyrosine phosphatase that has been shown previously to function as an inhibitor of phagocytosis. PTP-H2-mediated death of Sf-21 cells was blocked by the pan-caspase inhibitor benzyloxycarbonyl-Val-Ala-(O-methyl) Asp-fluoromethylketone (Z-VAD-FMK), but cells maintained in this inhibitor still exhibited a suppressed phagocytic response. Mutagenesis experiments indicated that the essential catalytic cysteine residue required for the phosphatase activity of PTP-H2 was required for apoptotic activity in Sf-21 cells. Loss of adhesion was insufficient to stimulate apoptosis of Sf-21 cells. PTP-H2 expression, however, did significantly reduce proliferation of Sf-21 cells, which could contribute to the apoptotic activity of this viral gene. Overall, our results indicate that specific genes expressed by MdBV induce apoptosis of certain insect cells and that this activity contributes to immunosuppression of hosts.

Proteomes of the aphid Macrosiphum euphorbiae in its resistance and susceptibility responses to differently compatible parasitoids

Host insects are either susceptible or resistant to parasitoids, where resistant hosts express immunity factors and compatible parasitoids express virulence factors that may reveal the manipulation of susceptible hosts. Using proteomics we compared responses of the same host, the aphid Macrosiphum euphorbiae, challenged by a well-adapted parasitoid Aphidius nigripes or by a less adapted relative, Aphidius ervi. The host was found to be equally acceptable to both parasitoids, but while A. nigripes normally developed and killed hosts (high susceptibility), development of the incompatible A. ervi was arrested at the primary egg stage (high resistance). Two-dimensional gels at two stages of parasitism revealed divergence in patterns of protein regulation of the M. euphorbiae host, responding to A. ervi or A. nigripes, with the greatest number of protein modulations in the host resistance response. In A. ervi-resistant hosts, proPO was strongly up-regulated, as were also three cuticle proteins, suggesting a PO basis and exoskeleton reinforcement as early and late responses of M. euphorbiae to the risk of parasitism. Resistance also correlated with up-regulation of antioxidative, energy-related, cytoskeleton and heat shock proteins. In A. nigripes-susceptible hosts, various proteins implicated in host and bacterial symbiont metabolism were significantly altered, suggesting complex host nutritional modulation. Over-expression of energy-related proteins also increased when A. nigripes established and developed. Aphid proteomes of compatible and incompatible Aphidius parasitism provide an integrative basis for consolidating our knowledge of host-parasitoid interactions.

Bracovirus gene products are highly divergent from insect proteins

Recently, several polydnavirus (PDV) genomes have been completely sequenced. The dsDNA circles enclosed in virus particles and injected by wasps into caterpillars appear to mainly encode virulence factors potentially involved in altering host immunity and/or development, thereby allowing the survival of the parasitoid larvae within the host tissues. Parasitoid wasps generally inject virulence factors produced in the venom gland. As PDV genomes are inherited vertically by wasps through a proviral form, wasp virulence genes may have been transferred to this chromosomal form, leading to their incorporation into virus particles. Indeed, many gene products from Cotesia congregata bracovirus (CcBV), such as PTPs, IkappaB-like, and cystatins, contain protein domains conserved in metazoans. Surprisingly however, CcBV virulence gene products are not more closely related to insect proteins than to human proteins. To determine whether the distance between CcBV and insect proteins is a specific feature of BV proteins or simply reflects a general high divergence of parasitoid wasp products, which might be due to parasitic lifestyle, we have analyzed the sequences of wasp genes obtained from a cDNA library. Wasp sequences having a high similarity with Apis mellifera genes involved in a variety of biological functions could be identified indicating that the high level of divergence observed for BV products is a hallmark of these viral proteins. We discuss how this divergence might be explained in the context of the current hypotheses on the origin and evolution of wasp-bracovirus associations.

Comparison of three methods of parasitoid polydnavirus genomic DNA isolation to facilitate polydnavirus genomic sequencing

A major long-term goal of polydnavirus (PDV) genome research is to identify novel virally encoded molecules that may serve as biopesticides to target insect pests that threaten agriculture and human health. As PDV viral replication in cell culture in vitro has not yet been achieved, several thousands of wasps must be dissected to yield enough viral DNA from the adult ovaries to carry out PDV genomic sequencing. This study compares three methods of PDV genomic DNA isolation for the PDV of Cotesia flavipes, which parasitizes the sugarcane borer, Diatraea saccharalis, preparatory to sequencing the C. flavipes bracovirus genome. Two of these protocols incorporate phenol-chloroform DNA extraction steps in the procedure and the third protocol uses a modified Qiagen DNA kit method to extract viral DNA. The latter method proved significantly less time-consuming and more cost-effective. Efforts are currently underway to bioengineer insect pathogenic viruses with PDV genes, so that their gene products will enhance baculovirus virulence for agricultural insect pests, either via suppression of the immune system of the host or by PDV-mediated induction of its developmental arrest. Sequencing a growing number of complete PDV genomes will enhance those efforts, which will be facilitated by the study reported here.

Transient expression of an EP1-like gene encoded in Cotesia plutellae bracovirus suppresses the hemocyte population in the diamondback moth, Plutella xylostella

A genome project has been launched and aims to sequence total genome of Cotesia plutellae bracovirus (CpBV). This on-going research has identified seven EP1-like (ELP) genes in the CpBV genome. A group of ELP genes has been speculated as an immunosuppressant encoded in Cotesia-associated bracoviruses. This study analyzed gene expression of these seven CpBV-ELPs in the parasitized diamondback moth, Plutella xylostella. Of these, six CpBV-ELPs were expressed in P. xylostella parasitized by C. plutellae. However, their expression levels varied in different tissues and parasitization stages. Especially, CpBV-ELP1 showed a persistent and ubiquitous expression pattern in both reverse transcriptase-polymerase chain reaction (RT-PCR) and immunofluorescence assays. When nonparasitized P. xylostella was transfected with a recombinant CpBV-ELP1 in a eukaryotic expression vector, CpBV-ELP1 was expressed for at least 3 days and the proteins were detectable in the cytoplasm of hemocytes. The transfected larvae showed significant reduction in total hemocyte numbers, compared with larvae injected with the cloning vector alone. Co-transfection with double-stranded RNA could knock down the expression of CpBV-ELP1 and prevented the reduction of the hemocyte population. This study demonstrates that CpBV-ELP1 plays a physiological role in suppressing host immune response presumably by its hemolytic activity during C. plutellae parasitization.

Transcriptional analysis of polydnaviral genes in the course of parasitization reveals segment-specific patterns

Polydnaviruses are symbiotic viruses of endoparasitic wasps, which are formed in their ovary and injected along with the eggs into the host. They manipulate the host in a way to allow successful parasitoid development. A hallmark of polydnaviruses is their segmented genome consisting of several circles of double-stranded DNA. We are studying the solitary egg-larval parasitoid Chelonus inanitus (Braconidae) parasitizing Spodoptera littoralis (Noctuidae). The polydnavirus of Chelonus inanitus (CiV) protects the parasitoid larva from encapsulation by the host's immune system, slightly modifies host nutritional physiology, and induces a developmental arrest of the host in the prepupal stage. Here we present data on newly identified CiV genes and their expression patterns in the course of parasitization. None of these genes has similarity to other genes and so far no gene families could be found. A rough estimation of transcript quantities revealed that even the most highly expressed CiV genes reach maximal values, which are 250 times lower than actin. This indicates that the CiV-induced alterations of the host are brought about by a concerted action of low levels of transcripts. In an overview, we show the expression patterns of all CiV genes analysed up to now; they indicate that several genes with similar expression patterns (early, persistent, intermediate, or late) are grouped together on the same segment. This is the first observation of this type. It suggests that one function of the segmentation of the polydnavirus genome may be the grouping together of genes, which are regulated in a similar manner.

Shared and species-specific features among ichnovirus genomes

During egg-laying, some endoparasitic wasps transmit a polydnavirus to their caterpillar host, causing physiological disturbances that benefit the wasp larva. Members of the two recognized polydnavirus taxa, ichnovirus (IV) and bracovirus (BV), have large, segmented, dsDNA genomes containing virulence genes expanded into families. A recent comparison of IV and BV genomes revealed taxon-specific features, but the IV database consisted primarily of the genome sequence of a single species, the Campoletis sonorensis IV (CsIV). Here we describe analyses of two additional IV genomes, the Hyposoter fugitivus IV (HfIV) and the Tranosema rostrale IV (TrIV), which we compare to the sequence previously reported for CsIV. The three IV genomes share several features including a low coding density, a strong A+T bias, similar estimated aggregate genome sizes ( approximately 250 kb) and the presence of nested genome segments. In addition, all three IV genomes contain members of six conserved gene families: repeat element, cysteine motif, viral innexin, viral ankyrin, N-family, and a newly defined putative family, the polar-residue-rich proteins. The three genomes, however, differ in their degree of segmentation, in within-family gene frequency and in the presence, in TrIV, of a unique gene family (TrV). These interspecific variations may reflect differences in parasite/host biology, including virus-induced pathologies in the latter.

Genomic and morphological features of a banchine polydnavirus: comparison with bracoviruses and ichnoviruses

Many ichneumonid and braconid endoparasitoids inject a polydnavirus (PDV) into their caterpillar hosts during oviposition. The viral entities carried by wasps of these families are referred to as "ichnoviruses" (IVs) and "bracoviruses" (BVs), respectively. All IV genomes characterized to date are found in wasps of the subfamily Campopleginae; consequently, little is known about PDVs found in wasps of the subfamily Banchinae, the only other ichneumonid taxon thus far shown to carry these viruses. Here we report on the genome sequence and virion morphology of a PDV carried by the banchine parasitoid Glypta fumiferanae. With an aggregate genome size of approximately 290 kb and 105 genome segments, this virus displays a degree of genome segmentation far greater than that reported for BVs or IVs. The size range of its genome segments is also lower than those in the latter two groups. As reported for other PDVs, the predicted open reading frames of this virus cluster into gene families, including the protein tyrosine phosphatase (PTP) and viral ankyrin (ank) families, but phylogenetic analysis indicates that ank genes of the G. fumiferanae virus are not embedded within the IV lineage, while its PTPs and those of BVs form distinct clusters. The banchine PDV genome also encodes a novel family of NTPase-like proteins displaying a pox-D5 domain. The unique genomic features of the first banchine virus examined, along with the morphological singularities of its virions (IV-like nucleocapsids, but enveloped in groups like some of the BVs), suggest that they could have an origin distinct from those of IVs and BVs.

Microplitis demolitor bracovirus genome segments vary in abundance and are individually packaged in virions

Polydnaviruses (PDVs) are distinguished by their unique association with parasitoid wasps and their segmented, double-stranded (ds) DNA genomes that are non-equimolar in abundance. Relatively little is actually known, however, about genome packaging or segment abundance of these viruses. Here, we conducted electron microscopy (EM) and real-time polymerase chain reaction (PCR) studies to characterize packaging and segment abundance of Microplitis demolitor bracovirus (MdBV). Like other PDVs, MdBV replicates in the ovaries of females where virions accumulate to form a suspension called calyx fluid. Wasps then inject a quantity of calyx fluid when ovipositing into hosts. The MdBV genome consists of 15 segments that range from 3.6 (segment A) to 34.3 kb (segment O). EM analysis indicated that MdBV virions contain a single nucleocapsid that encapsidates one circular DNA of variable size. We developed a semi-quantitative real-time PCR assay using SYBR Green I. This assay indicated that five (J, O, H, N and B) segments of the MdBV genome accounted for more than 60% of the viral DNAs in calyx fluid. Estimates of relative segment abundance using our real-time PCR assay were also very similar to DNA size distributions determined from micrographs. Analysis of parasitized Pseudoplusia includens larvae indicated that copy number of MdBV segments C, B and J varied between hosts but their relative abundance within a host was virtually identical to their abundance in calyx fluid. Among-tissue assays indicated that each viral segment was most abundant in hemocytes and least abundant in salivary glands. However, the relative abundance of each segment to one another was similar in all tissues. We also found no clear relationship between MdBV segment and transcript abundance in hemocytes and fat body.

The ancient Virus World and evolution of cells

Background

Recent advances in genomics of viruses and cellular life forms have greatly stimulated interest in the origins and evolution of viruses and, for the first time, offer an opportunity for a data-driven exploration of the deepest roots of viruses. Here we briefly review the current views of virus evolution and propose a new, coherent scenario that appears to be best compatible with comparative-genomic data and is naturally linked to models of cellular evolution that, from independent considerations, seem to be the most parsimonious among the existing ones.

Results

Several genes coding for key proteins involved in viral replication and morphogenesis as well as the major capsid protein of icosahedral virions are shared by many groups of RNA and DNA viruses but are missing in cellular life forms. On the basis of this key observation and the data on extensive genetic exchange between diverse viruses, we propose the concept of the ancient virus world. The virus world is construed as a distinct contingent of viral genes that continuously retained its identity throughout the entire history of life. Under this concept, the principal lineages of viruses and related selfish agents emerged from the primordial pool of primitive genetic elements, the ancestors of both cellular and viral genes. Thus, notwithstanding the numerous gene exchanges and acquisitions attributed to later stages of evolution, most, if not all, modern viruses and other selfish agents are inferred to descend from elements that belonged to the primordial genetic pool. In this pool, RNA viruses would evolve first, followed by retroid elements, and DNA viruses. The Virus World concept is predicated on a model of early evolution whereby emergence of substantial genetic diversity antedates the advent of full-fledged cells, allowing for extensive gene mixing at this early stage of evolution. We outline a scenario of the origin of the main classes of viruses in conjunction with a specific model of precellular evolution under which the primordial gene pool dwelled in a network of inorganic compartments. Somewhat paradoxically, under this scenario, we surmise that selfish genetic elements ancestral to viruses evolved prior to typical cells, to become intracellular parasites once bacteria and archaea arrived at the scene. Selection against excessively aggressive parasites that would kill off the host ensembles of genetic elements would lead to early evolution of temperate virus-like agents and primitive defense mechanisms, possibly, based on the RNA interference principle. The emergence of the eukaryotic cell is construed as the second melting pot of virus evolution from which the major groups of eukaryotic viruses originated as a result of extensive recombination of genes from various bacteriophages, archaeal viruses, plasmids, and the evolving eukaryotic genomes. Again, this vision is predicated on a specific model of the emergence of eukaryotic cell under which archaeo-bacterial symbiosis was the starting point of eukaryogenesis, a scenario that appears to be best compatible with the data.

Conclusion

The existence of several genes that are central to virus replication and structure, are shared by a broad variety of viruses but are missing from cellular genomes (virus hallmark genes) suggests the model of an ancient virus world, a flow of virus-specific genes that went uninterrupted from the precellular stage of life's evolution to this day. This concept is tightly linked to two key conjectures on evolution of cells: existence of a complex, precellular, compartmentalized but extensively mixing and recombining pool of genes, and origin of the eukaryotic cell by archaeo-bacterial fusion. The virus world concept and these models of major transitions in the evolution of cells provide complementary pieces of an emerging coherent picture of life's history.

Reviewers

W. Ford Doolittle, J. Peter Gogarten, and Arcady Mushegian.

Evolution of complexity in the viral world: the dawn of a new vision

Recent sequencing of the genomes of numerous large viruses provide for unprecedented opportunities to study the emergence and evolution of complexity in the virus world. This special issue of Virus Research explores trends in the evolution of complex genomes in most major classes of viruses.