Next Article in Journal
Evaluating Environmental Persistence and Disinfection of the Ebola Virus Makona Variant
Next Article in Special Issue
Genome Variability and Gene Content in Chordopoxviruses: Dependence on Microsatellites
Previous Article in Journal
Insect-Specific Flaviviruses: A Systematic Review of Their Discovery, Host Range, Mode of Transmission, Superinfection Exclusion Potential and Genomic Organization
Previous Article in Special Issue
The Evolution of Poxvirus Vaccines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gene Acquisition Convergence between Entomopoxviruses and Baculoviruses

1
Institut de Recherche sur la Biologie de l'Insecte, CNRS UMR 7261, Université François-Rabelais, UFR Sciences et Techniques, 37200 Tours, France
2
Université de Poitiers, UMR CNRS 7267 Ecologie et Biologie des Interactions, Equipe Ecologie Evolution Symbiose, 86073 Poitiers, France
3
Forestry and Forest Products Research Institute, Matsunosato, Tsukuba, Ibaraki 305-8687, Japan
4
Institute of Agriculture, Tokyo University of Agriculture and Technology, Saiwai, Fuchu, Tokyo 183-8509, Japan
5
Great Lakes Forestry Centre, Sault Sainte Marie, ON P6A 2E5, Canada
*
Author to whom correspondence should be addressed.
Present address: Department of Zoology, University of Oxford, South Parks Road, Oxford, OX1 3PS, UK.
Submission received: 17 December 2014 / Revised: 25 March 2015 / Accepted: 8 April 2015 / Published: 13 April 2015
(This article belongs to the Special Issue Poxvirus Evolution)

Abstract

:
Organisms from diverse phylogenetic origins can thrive within the same ecological niches. They might be induced to evolve convergent adaptations in response to a similar landscape of selective pressures. Their genomes should bear the signature of this process. The study of unrelated virus lineages infecting the same host panels guarantees a clear identification of phyletically independent convergent adaptation. Here, we investigate the evolutionary history of genes in the accessory genome shared by unrelated insect large dsDNA viruses: the entomopoxviruses (EPVs, Poxviridae) and the baculoviruses (BVs). EPVs and BVs have overlapping ecological niches and have independently evolved similar infection processes. They are, in theory, subjected to the same selective pressures from their host’s immune responses. Their accessory genomes might, therefore, bear analogous genomic signatures of convergent adaption and could point out key genomic mechanisms of adaptation hitherto undetected in viruses. We uncovered 32 homologous, yet independent acquisitions of genes originating from insect hosts, different eukaryotes, bacteria and viruses. We showed different evolutionary levels of gene acquisition convergence in these viruses, underlining a continuous evolutionary process. We found both recent and ancient gene acquisitions possibly involved to the adaptation to both specific and distantly related hosts. Multidirectional and multipartite gene exchange networks appear to constantly drive exogenous gene assimilations, bringing key adaptive innovations and shaping the life histories of large DNA viruses. This evolutionary process might lead to genome level adaptive convergence.

Graphical Abstract

1. Introduction

All cellular and viral species are engaged in multitrophic interactions involving a wide diversity of beneficial or antagonist partners [1]. It is therefore not uncommon for many parasites, including viruses, of unrelated families to thrive within the same host species [2]. In theory, sympatric parasites subjected to similar ecological landscapes (i.e., host immune responses) may select similar mutations to adapt to their hosts [3,4]. The strongest evidence of adaptation though is the convergent evolution in independent lineages of similar traits and genomic features [5,6].
Involved in antagonistic yet durable interactions with their hosts, viruses are strongly locked in a coevolutionary fitness race [7], requiring rapid and efficient genomic innovations. Large dsDNA viruses with their relatively fluid genomes provide sufficient genomic information to test if adaptation towards the same hosts might have produced similar traces on their genomes, particularly the acquisition of beneficial genes [8]. So far, most studies have focused on the convergent acquisition of specific beneficial genes between phylogenetically related viruses infecting Bacteria [9,10]. However, convergence at the genomic level is rarely demonstrated because phylogenetic independence is seldom achieved.
For pathogens, insects are taxonomically the most diverse group of potential mutlicellular eukaryotic host species [11]. Interest in insect pest pathogens for biocontrol applications [12] generated an abundance of viral genomic data. The Entomopoxvirinae (EPV) a sub-family of the Poxviridae and the family Baculoviridae (BV) are large dsDNA viruses exclusively pathogenic to insects [13]. Similar in size, EPVs and BVs, respectively, encode 250–350 and 100–200 genes on 230–330 and 80–180 kb genomes [13,14]. However, their divergent genomic structure (linear versus circular) and gene content indicate the virus families are phylogenetically unrelated, and probably have different origins with no common ancestors [15]. The EPV and BV core genomes, containing 49 and 37 essential genes [16,17], only overlap for the DNA polymerase, sustaining the replication of their DNA genomes. The remaining genes found in EPV and BV genomes have not been inherited from the origin of either families and form the so-called accessory genomes [18]. EPV and BV share a number of homologous accessory genes [19], including some similar to cellular genes possibly acquired from horizontal gene transfers (HGTs) [20,21,22,23]. Locked in a coevolutionary arms-race, it is presumed that viruses retain only the genes that increase their fitness [8] in a selective trade-off between adaptive novelties and efficient genome replication. This indicates the adaptive role of accessory genomes. EPVs and BVs have largely overlapping ecological niches and can infect the same host species. They might therefore bear analogous signatures of their adaption to these hosts immune response. Furthermore, they have independently evolved a similar infection process. The onset of infections occurs after the ingestion by insect larvae of virus particles, embedded within non-homologous occlusion bodies (OBs), which confer some protection against inactivating environmental factors such as heat, desiccation and UV light [24].
Here we investigated the convergent evolution of the accessory genome content of unrelated large dsDNA viruses based on four recently available EPV genomes [14], infecting the same moths as previously sequenced BVs. This guarantees a clear phyletically independent identification of common genomic features. Furthermore, we estimated the putative phyletic origin of shared accessory genes to determine whether particular gene acquisition processes might be favored. Finally, we discuss the potential role played by mobile elements in gene exchange networks.

2. Materials and Methods

2.1. Viral Orthologous Clustering

An in-house bioinformatic pipeline was written to cluster together homologous EPV and BV proteins and to search for cellular and viral homologs of each EPV/BV clusters. First, we performed a protein secondary structure alignment clustering, using the jackhmmer program of the hmmer3 package [25], in an all-against-all alignment (threshold: e-value = 10−6) of all protein members of Entomopoxvirinae sub-family and Baculoviridae family found in the non-redundant protein database of the NCBI, including the newly described Adoxophyes honmai EPV (AHEV), Choristoneura biennis EPV (CBEV), Choristoneura rosaceana EPV (CREV) and Mythimna separata (MySEV) predicted proteins [14]. Then, with the ClustalOmega program [26], protein multiple alignments were conducted on clusters shared by both virus families, prior to hidden Markov chain (HMM) profile building. By using hmmsearch program of the hmmer3 package, HMM profiles were aligned to the non-redundant protein database and to the insect EST and TSA databases of the NCBI, both of which were initially converted into protein sequences using the metagenomic option of Prodigal program [27]. Finally, for each sequence found, we searched for taxonomic information and selected six categories: (1) other insect large dsDNA (including nudiviruses, ascoviruses, iridoviruses and hytrosaviruses); (2) non-insect large dsDNA viruses; (3) invertebrates; (4) vertebrates; (5) plants’ and (6) bacteria. If needed, we limited the results to the top 20 matches for each category for efficient phylogenetic reconstruction.

2.2. Phylogenetic Analyses

Maximum likelihood (ML) phylogenetic inferences were performed on each HMM profile and its associated matches, after conversion into protein multiple alignments. ML analyses were performed with RAxML 7.4.2 [28], using the substitution model and model parameters Whelan And Goldman model with Gamma distribution (WAG+G) and support for node in ML trees were obtained from 100 bootstrap iterations.

3. Results

3.1. Entomopoxviruses and Baculoviruses Share Accessory Genes

We developed a pipeline, to first cluster all available EPV and BV proteins by secondary structure alignment and to align them against different NCBI databases. Out of 487 viral protein clusters, we found 33 clusters of homologous genes shared by EPVs and BVs (blastp e-value threshold: 1e-6), including 26 with similarities to cellular proteins and 21 with similarities to other families of large dsDNA viruses (Table 1). As expected, we found only one core gene, the DNA polymerase, among these homologous clusters. The remaining genes were from the accessory genomes. Several genes encoded essential cellular enzymes such as: the Cu/Zn superoxide dismutase, the ribonucleotide reductase small subunit, the dUTPase, the ubiquitin and the matrixin metalloproteinase. We identified genes potentially involved in insect immune interactions including two families of apoptosis inhibitors (BV iap3 and iap4 and p35/p49), two kinds of protein tyrosine phosphatase, conotoxins and a phage antirepressor. We found multicopy viral gene families, like the bro (baculovirus repeated ORFs), MTG (myeloid translocation gene) motif, leucine rich, and nla (nuclear localization of G-actin). Furthermore, we detected two genes related to host ecology: the fusolin and the DNA photolyase. The remaining genes had previously been described as unknown BV genes.

3.2. Viruses Infecting the Same Hosts Share Homologous Genes

Four EPVs in our dataset had overlapping host species with BVs, allowing the investigation of gene convergence in sympatry. The cluster “unknown XecnGV orf138” (xc138) was the only one to show substantially higher sequence similarities (Table 1, e-value = 0) between an EPV, MySEV and a BV, Xestia c-nigrum granulovirus (XecnGV) (92% over 650 aa), than between XecnGV and its close relative Pseudaletia unipuncta granulovirus (PsunGV) (71% over 650 aa). Singularly, these three viruses are able to infect the host Mythimna separata [29] as well as a number of other noctuid moths. Pseudogenized forms of xc138 were detected in the genome of Helicoverpa armigera granulovirus (HearGV), which is a sister species of XecnGV, and in the genome of Heliothis virescens ascovirus 3g (Ascoviridae, HvAV-3g). Interestingly both viruses can infect M. separata in the lab though they have only been isolated from H. armigera and Spodoptera exigua, respectively. More distant baculovirus xc138 homologs were also identified in Chrysodeixis chalcites nucleopolyhedrovirus and Orgyia leucostigma nucleopolyhedrovirus, and a second copy in XecnGV.
The acquisition of a xc138 gene copy might confer to MySEV and XecnGV specific virulence towards their hosts, including for example M. separata. In different ecological contexts, as in HearGV, which is more specific to H. armigera, the removal of particular host pressure could lead to gene loss in a fitness trade-off between genome replication efficacy and host adaptation. We first hypothesized that both EPVs and BVs might have independently acquired the xc138 gene from their hosts. However, we could not detect any xc138 cellular homolog either by bioinformatics data mining or by PCR from M. separata, M. unipuncta, M. riparia, Xestia c-nigrum, Trichoplusia ni, Chrysodeixis chalcites or Agrotis segetum genomic DNA. In contrast, we were able to amplify a xc138 fragment in a closely related EPV of MySEV: Heliothis armigera entomopoxvirus (HAEV) (Genbank accession number: KJ683043).
Table 1. Origin of homologous gene clusters found in EPVs and BVs.
Table 1. Origin of homologous gene clusters found in EPVs and BVs.
Viral homologous gene cluster nameBest blastp hit between EPVs and BVs (e-value)Presence in other large dsDNA virusesPutative origins of EPV/BV genesPhylogeny (Figure)
DNA polymerase6 × 10−9X--
bro gene family protein9 × 10−26XBacteriaS1
helicase 21 × 10−67XBacteriaS2
matrixin metalloproteinase1 × 10−9XBacteriaS3
nla gene he65 (AcMNPV orf105)8 × 10−97-BacteriaS4
protein tyrosine phosphatase 2 (ptp-2)2 × 10−33XBacteriaS5
putative phage antirepressor9 × 10−26XBacteriaS6
unknown LdMNPV orf1292 × 10−57-BacteriaS7
chitin binding protein (AcMNPV orf145)1 × 10−17XInsectaS8
Cu/Zn superoxide dismutase4 × 10−56XInsecta2b/S9
inhibitor of apoptosis4 × 10−79XInsectaS10
MTG motif gene family protein7 × 10−10XInsectaS11
ubiquitin6 × 10−38XInsectaS12
unknown AcMNPV orf75 × 10−9-InsectaS13
unknown XecnGV orf1063 × 10−13XInsectaS14
unknown XecnGV orf226 × 10−103XInsectaS15
acetyltransferase7 × 10−13-Insecta (Lepidoptera)S16
protein phosphatase 1, regulary subunit 15A5 × 10−15XInsecta (Lepidoptera)S17
ribonucleotide reductase small subunit homolog8 × 10−134XInsecta (Lepidoptera)2a/S18
unknown XecnGV orf728 × 10−133-Insecta (Lepidoptera)S19
protein tyrosine phosphatase 1 (ptp-1)4 × 10−41XInsecta (Lepidoptera) & EukaryoteS20
DNA photolyase3 × 10−140XBilateria/Insecta (Lepidoptera)S21
dUTPase1 × 10−43XInsecta/UnknownS22
unknown ClanGV orf0859 × 10−17XEukaryoteS23
leucine rich gene family protein4 × 10−27-Unicellular eukaryoteS24
fusolin/spindlin/gp37 (AcMNPV orf64)4 × 10−67-Unicellular eukaryote (Amoebozoa)2c/S25
conotoxin-like protein3 × 10−20-VirusS26
p35/p49 apoptosis inhibitor4 × 10−30-VirusS27
unknown AcMNPV orf188 × 10−11-VirusS28
unknown AdorNPV orf1108 × 10−11XVirusS29
unknown AgseGV orf47 × 10−103-VirusNA
unknown ChocGV orf113 × 10−14-VirusNA
unknown XecnGV orf1380XVirus1a
The comparison of the DNA polymerase phylogeny, representing the virus species tree, with the xc138 phylogeny highlighted two xc138 HGTs have occurred between EPVs and BVs (Figure 1a). Reconciliation analyses by bipartition dissimilarity [30] proposed that xc138 was transferred once between MySEV and the ancestor of XecnGV and HearGV, and a second time between HAEV and PsunGV (Figure 1a,b). The direction of HGTs between these EPVs and BVs is unknown because the presence of a xc138 gene copy in the common ancestors of both EPV and BV lineages cannot be established. However, this suggests unrelated viruses, such as XecnGV and MySEV, sharing the same hosts can exchange specific genes such as xc138. Its subsequent loss in HearGV, suggest xc138 might be involved in viral adaptation to particular ecological context and convergently maintained in the genomes of XecnGV and MySEV.
Figure 1. Acquisition convergence of the gene xc138. (a) Phylogenetic reconciliation between xc138 gene tree (left) and DNA polymerase gene tree, representing the virus species tree (right) shows two horizontal gene transfers between Mythimna separata entomopoxvirus (MySEV) and the Xestia c-nigrum granulovirus (XecnGV)/Helicoverpa armigera granulovirus (HearGV) common ancestor (in orange) and between Heliothis armigera entomopoxvirus (HAEV) and Pseudaletia unipuncta granulovirus (PsunGV) (in blue), both with high confidence supports (100%); (b) Schematic phylogenies of Poxviridae and Baculoviridae families show the most parsimonious evolutionary history of xc138 gene. The black circles denote the acquisition of the xc138 gene by the entomopoxvirus or baculovirus lineages from an unknown origin and later the two recent horizontal gene transfers between entomopoxviruses and baculoviruses.
Figure 1. Acquisition convergence of the gene xc138. (a) Phylogenetic reconciliation between xc138 gene tree (left) and DNA polymerase gene tree, representing the virus species tree (right) shows two horizontal gene transfers between Mythimna separata entomopoxvirus (MySEV) and the Xestia c-nigrum granulovirus (XecnGV)/Helicoverpa armigera granulovirus (HearGV) common ancestor (in orange) and between Heliothis armigera entomopoxvirus (HAEV) and Pseudaletia unipuncta granulovirus (PsunGV) (in blue), both with high confidence supports (100%); (b) Schematic phylogenies of Poxviridae and Baculoviridae families show the most parsimonious evolutionary history of xc138 gene. The black circles denote the acquisition of the xc138 gene by the entomopoxvirus or baculovirus lineages from an unknown origin and later the two recent horizontal gene transfers between entomopoxviruses and baculoviruses.
Viruses 07 01960 g001

3.3. Convergence between Viruses Infecting Different Hosts

At the evolutionary timescale, gene acquisition convergence could point out adaptations linked to essential conserved virus-host interactions. We performed maximum likelihood phylogenetic inferences on each of the 31 remaining protein clusters shared between EPVs and BVs (Table 1, Supporting information Figures S1–S29). Unexpectedly, we found 15 genes had putative insect origins, suggesting they were derived from HGTs from host genomes to the virus genomes. Three genes had putative eukaryotic origin outside of the metazoa and seven putatively derived from bacteria, suggesting viruses could have acquired them from the environment or from co-infecting agents. Seven genes, including xc138, had no cellular homologs, and were, therefore, assigned a putative viral origin.
The ribonucleotide reductase small subunit phylogeny (rr2, Figure 2a and Supporting information Figure S18) shows the paraphyly of the insects with multiple inclusions of lepidopteran virus clades (BVs, EPVs and ascoviruses) in the Lepidoptera clade. This implies independent and relatively recent rr2 HGTs from diverse moth genomes to various large DNA viruses, suggesting convergent adaptive acquisitions of host genes across several unrelated virus families infecting a spectrum of related hosts.
Figure 2. Schematic phylogenies of rr2, sod and fusolin genes. The trees were obtained from maximum likelihood inferences of the multiple amino acid alignments of (a) the rr2; (b) the sod and (c) the fusolin gene products. Support for nodes above branches indicate maximum likelihood nonparametric bootstraps (100 replicates). Bootstraps values between 50 and 75 and upper to 75 are denoted with * and ** respectively.
Figure 2. Schematic phylogenies of rr2, sod and fusolin genes. The trees were obtained from maximum likelihood inferences of the multiple amino acid alignments of (a) the rr2; (b) the sod and (c) the fusolin gene products. Support for nodes above branches indicate maximum likelihood nonparametric bootstraps (100 replicates). Bootstraps values between 50 and 75 and upper to 75 are denoted with * and ** respectively.
Viruses 07 01960 g002
Similarly, the Cu/Zn-superoxide dismutase phylogeny (sod, Figure 2b and Supporting information Figure S9) shows the inclusion of lepidopteran EPV and BV lineages in the insect sod clade, and thus the convergent adaptive acquisitions of insect sod genes by EPVs and BVs. Moreover, within the BV clade, the sod phylogeny reflects the Baculoviridae phylogeny [13]. This implies the BV sod gene was acquired from an ancestral insect sod gene early in baculovirus history, and possibly as early as 90 million years ago [31]. During this time viral sequences have greatly diverged from host sequences, which might explain the low bootstrap support we could recover for the deepest nodes of the tree.
Strikingly both the rr2 and sod phylogeny also reveal higher levels of convergent adaptive gene acquisitions. The inclusions of a dipteran virus clade (Hytrosaviridae family) within the Diptera clade in the rr2 tree (Figure 2a and Supporting information Figure S18) and of the Chordopoxvirus (ChPV) clade within the vertebrate lineage in the sod tree (Figure 2b and Supporting information Figure S9) indicate convergent gene acquisitions can be independent of the type of host. Both rr2 and sod genes have been recurrently acquired from their hosts by unrelated viruses infecting both insects and vertebrates. This is strong evidence of viral adaptive convergence towards essential cellular processes encountered in all metazoans.
The fusolin phylogeny (Figure 2c and Supporting information Figure S25) shows convergent ancestral acquisitions of the fusolin gene occurred early in the evolution of both EPVs and BVs, as it reflects both the Entomopoxvirinae [14] and Baculoviridae phylogenies [13]. The fusolin appears absent in metazoans and its closest cellular relatives are chitin-binding proteins of amoeba and bacteria. This suggests EPVs and BVs have captured the fusolin from an alternate cellular organism than their natural insect hosts.

4. Discussion

Previous studies on genomic convergent evolution highlighted similar mutations and HGTs outcomes for bacteria sharing particular ecological niches [3,4]. Phyletic constraint could have imposed convergent adaptation on the bacteria, due to their shared history as cellular organisms. As viruses have multiple independent origins [32], we studied unrelated virus lineages infecting the same panel of hosts to identify phyletically independent convergent adaptations underpinned by key genomic adaptive mechanisms, as yet undetected in eukaryotic viruses. Although they are phylogenetically unrelated, EPVs and BVs share remarkably similar replication strategies. Both EPVs and BVs produce pseudo-crystalline OBs for per os infection of insect larvae as well as spindle bodies enhancing the initial stages of midgut infection [24]. As this phenotypic convergence should be traceable at the genomic level, we looked for homologous genes between EPVs and BVs and investigated their evolutionary history using phylogenetic analyses.
We identified 32 clusters of homologous genes present in the accessory genomes of both EPVs and BVs (Table 1) signaling substantial convergent gene acquisitions. Strikingly, the xc138 gene of XecnGV is 92% identical to that of MySEV, suggesting a recent HGT between these virus lineages (Figure 1b). Notably, the presence of xc138 might be correlated with virus host range, as both MySEV and XecnGV are able to infect M. separata. Close homologs, recently acquired by other EPVs, BVs and ascoviruses, might be implicated in specific host interactions as all infect closely related noctuid host species. However these five viruses carrying xc138 are all able to infect a larger range of host species than most of EPVs and BVs that do not carry this gene. Thus, this gene might alternatively confer a wider host range. Further functional analyses are however necessary to specify the phenotype linked with this gene.
Besides these recent HGTs, we found gene acquisition convergence anchored deeper in time, early in the history of virus families, as illustrated by the rr2, sod and fusolin phylogenies (Figure 2). In particular, the sod phylogeny (Figure 2b) underlines ancestral gene acquisition convergence between unrelated virus families towards insects, but also within the Poxviridae between two virus lineages of non-overlapping host spectrum (vertebrates and insects). Markedly, phylogenetic analyses on the origin of common EPV and BV accessory genes show HGT from insect hosts to be the most common source of exogenous gene assimilations (Table 1). But viruses appear to also have acquired genes from cellular organisms other than their natural hosts. This implies that eukaryotic large DNA viruses could undergo HGTs all across the tree of life [33].
Despite extensive effort to find cellular homologs, some genes seem exclusively viral (Table 1). Their presence, in independent phylogenetic lineages, shows gene exchanges could occur between viruses. For example the xc138 and fusolin trees (Figure 1 and Figure 2c) showed HGTs from EPVs to BVs. In the case of the fusolin, codon bias studies had suggested the HGT [34]. Furthermore, in EPVs the product of the fusolin are spindle bodies, which can enhance BV per oral infectivity [35]. This implies the immediate adaptive benefit BV could have gained from the fusolin HGT from an EPV. Apparently absent from metazoans, the fusolin closest cellular relatives are chitin-binding proteins of amoeba and bacteria. Thus EPVs might have captured an ancestral form of the fusolin from other cellular organisms than their natural hosts. The window for such a gene capture is necessarily smaller than for a gene from a natural host, reducing the opportunity for recurrent gene transfer into different virus families. Consequently, the opportunity to acquire this gene from a cellular organism might have happened only once thus in only one virus family. Fusolin confers a significant adaptive advantage in per oral infectivity [36], which is the natural route of infection for both EPVs and BVs. Transmission of an ancestral fusolin gene from EPVs to BVs could have occurred afterwards, but still early in the evolution of both groups, allowing the BV lineage to adapt more efficiently to lepidopteran hosts. Even if gene transfers between viruses are possible, we posit that they are relatively rare events in comparison with HGT from hosts, which are the majority in our dataset. The EPV/BV models are particularly appropriate for observing these outstanding exchanges, as they have similar infection processes and mixed infections of the same hosts occur [37,38]. Though replicating in different cell compartments (BVs in the nucleus and EPVs in the cytoplasm), mixed infections of the same cells could support viral gene exchanges, although the molecular mechanisms remain elusive.
Gene acquisition seems to be a widespread phenomenon in large dsDNA viruses. However, it is presumed that viruses retain only the genes that increase their fitness. So during the course of evolutionary history, genes have profusely been captured and lost. Only those genes conferring substantial adaptive benefits could spread within virus populations through natural selection. Previous studies have shown different gene exchanges involving virus, multicellular eukaryotes and unicellular organisms. Gene transfers can occur from viruses to multicellular eukaryote hosts, as in the case of endogenous viral elements [39] and by whole genome integration, as in polydnaviruses [6]. Recurrent HGTs have also been shown between bacteria and multicellular eukaryote hosts [40] and between viruses and bacteria [41]. As viruses interact with other viruses and cellular organisms, within what makes up the host hologenome [42], each interacting partner becomes a potential source of genes. We propose that these hologenomes support networks of multipartite and multidirectional gene exchanges (Figure 3) providing a large reservoir of new potentially adaptive genes to all partners including viruses.
Several families of transposable elements integrate into BV genomes during infection cycles and are thus present at low frequency within BV populations [43]. It is possible that other large DNA viruses, such as EPV could also be the targets of transposable elements. Moreover, the discovery of new mobile elements has recently soared and their implication in genome evolution seems substantial [44,45]. The genomes of these mobile elements can include genes from various origins (eukaryotic, bacterial and viral) that could have been acquired from host genomes as well as from microbial genomes. These mobile elements spreading within the genomes and populations of various organisms, including viruses [46], could be assimilated afterwards, allowing the delivery of horizontally transferred genes. Although research on these elements remains in its infancy, they could possibly facilitate gene exchanges during coinfections between viruses replicating in different cellular compartments, such as EPVs and BVs, and thus participate in the evolution of gene contents within hologenomes.
Figure 3. Hologenome gene exchange networks. Many routes of horizontal gene transfers within the hologenome of host species and its interacting partners can lead to gene acquisition convergence between entomopoxviruses and baculoviruses. The thickness of arrows indicates the occurrence of horizontal gene transfers (bold: very frequent, regular: frequent) and dotted arrows indicate possible gene exchanges.
Figure 3. Hologenome gene exchange networks. Many routes of horizontal gene transfers within the hologenome of host species and its interacting partners can lead to gene acquisition convergence between entomopoxviruses and baculoviruses. The thickness of arrows indicates the occurrence of horizontal gene transfers (bold: very frequent, regular: frequent) and dotted arrows indicate possible gene exchanges.
Viruses 07 01960 g003
Depending on their phylogenetic origin, the genes of large dsDNA virus are broadly classified as belonging to the essential core genome or to the accessory genome. By containing beneficial non-essential genes, the accessory genome is in essence more plastic than the core genome and is derived principally from the assimilation (acquisition and adaptive molecular evolution) of exogenous genes. Viruses probably rely on this accessory gene reservoir to adapt to ecological niches. The signatures of adaptive evolution should be obvious in the accessory genome. Accessory genes have been found to undergo episodes of positive selection [21], duplication and loss events [47] that participate in virus genome dynamics. We find that unrelated large dsDNA viruses infecting closely or distantly related hosts converge towards homologous genomic adaptation. This demonstrates adaptive gene assimilation is a general evolutionary process. Our study based on EPVs and BVs with unbalanced numbers of sequenced genomes (~60 for BVs and only 6 for EPVs), identified 32 homologous genes between EPVs and BVs (Table 1), which could perform analogous adaptive functions. The facts that 19 of these genes were found in other large dsDNA virus families, including the chordopoxviruses, nudiviruses, ascoviruses, iridoviruses and hytrosaviruses further emphasize the recurrence of adaptive convergence and the fundamental influence of ecological constraints on genome evolution.
We observed adaptive HGTs throughout the evolutionary history of insect viruses. Some early HGTs are phylogenetically conserved after assimilation and might have been instrumental in overcoming particularly complex infections barriers, such as host innate immunity. The production of superoxide [48] and induction of apoptosis are innate host mechanisms that control viral infection [49]. The recurrent convergent acquisitions and evolutionary maintenance of the sod and antiapoptotic genes (iap and p35) [50] in viruses are evidence of the evolutionary leaps that follow adaptive gene assimilation. Other genes such as the different ribonucleotide reductase subunits or the dUTPase are periodically acquired and lost, showing their importance for viral replication success might depend more on specific environmental host conditions. The similarity in host ecology also drives the accessory genome convergence, as in the case of fusolin, enhancing per oral infectivity.
Our definition of networks of multipartite and multidirectional gene exchanges supports that the accessory genomes of large dsDNA virus can be considered as collections of gene acquired from cells and other viruses, and which can be expanded via intragenomic duplications. The consequence is an entanglement of the origins of viral genes. This could explain why, so far, the evolution of the first viruses remains mysterious. Closer to us, seeking adaptive gene convergence in viruses adapted to specific host species could bring new developments in virus biocontrol and in vaccines applications. HGTs and convergence in cellular and viral gene acquisition gave major adaptive impulses to large dsDNA virus evolution, fostering evolutionary leaps in virus life histories.

Supplementary Files

Supplementary File 1

Acknowledgments

European Research Council grant 205206 GENOVIR funded J.T. and E.A.H. JSPS KAKENHI Grant numbers 24380035 and 24657064 supported J.T. We thank Juan Pino for providing noctuid moth samples.

Author Contributions

J.T. and E.A.H. conceived the analyses; J.T. analyzed the data; J.T., M.N. and B.A. contributed material; J.T., B.A. and E.A.H. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Eckburg, P.B.; Bik, E.M.; Bernstein, C.N.; Purdom, E.; Dethlefsen, L.; Sargent, M.; Gill, S.R.; Nelson, K.E.; Relman, D.A. Diversity of the human intestinal microbial flora. Science 2005, 308, 1635–1638. [Google Scholar]
  2. Cox, F.E.G. Concomitant infections, parasites and immune responses. Parasitology 2011, 122, S23–S38. [Google Scholar] [CrossRef]
  3. Tenaillon, O.; Rodríguez-Verdugo, A.; Gaut, R.L.; McDonald, P.; Bennett, A.F.; Long, A.D.; Gaut, B.S. The molecular diversity of adaptive convergence. Science 2012, 335, 457–461. [Google Scholar] [CrossRef] [PubMed]
  4. Smillie, C.S.; Smith, M.B.; Friedman, J.; Cordero, O.X.; David, L.A.; Alm, E.J. Ecology drives a global network of gene exchange connecting the human microbiome. Nature 2011, 480, 241–244. [Google Scholar] [CrossRef] [PubMed]
  5. Parker, J.; Tsagkogeorga, G.; Cotton, J.A.; Liu, Y.; Provero, P.; Stupka, E.; Rossiter, S.J. Genome-wide signatures of convergent evolution in echolocating mammals. Nature 2013, 502, 228–231. [Google Scholar] [CrossRef] [PubMed]
  6. Herniou, E.A.; Huguet, E.; Thézé, J.; Bézier, A.; Periquet, G.; Drezen, J.-M. When parasitic wasps hijacked viruses: Genomic and functional evolution of polydnaviruses. Phil. Trans. R. Soc. B 2013, 368, e20130051. [Google Scholar] [CrossRef]
  7. Combes, C. Parasitism: The Ecology and Evolution of Intimate Interactions; University of Chicago Press: Chicago, IL, USA, 2001. [Google Scholar]
  8. Shackelton, L.A.; Holmes, E.C. The evolution of large DNA viruses: Combining genomic information of viruses and their hosts. Trends Microbiol. 2004, 12, 458–465. [Google Scholar] [CrossRef] [PubMed]
  9. Sabehi, G.; Shaulov, L.; Silver, D.H.; Yanai, I.; Harel, A.; Lindell, D. A novel lineage of myoviruses infecting cyanobacteria is widespread in the oceans. Proc. Natl. Acad. Sci. USA 2012, 109, 2037–2042. [Google Scholar] [CrossRef] [PubMed]
  10. Anantharaman, K.; Duhaime, M.B.; Breier, J.A.; Wendt, K.A.; Toner, B.M.; Dick, G.J. Sulfur oxidation genes in diverse deep-sea viruses. Science 2014, 344, 757–760. [Google Scholar] [CrossRef] [PubMed]
  11. May, R.M. How many species are there on earth? Science 1988, 241, 1441. [Google Scholar] [CrossRef] [PubMed]
  12. Cory, J. Use of baculoviruses as biological insecticides. Mol. Biotechnol. 1997, 7, 303–313. [Google Scholar] [CrossRef] [PubMed]
  13. King, A.M.; Lefkowitz, E.; Adams, M.J.; Carstens, E.B. Virus Taxonomy: Ninth Report of the International Committee of Taxonomy of Viruses; Elsevier: Amsterdam, The Netherlands, 2011. [Google Scholar]
  14. Thézé, J.; Takatsuka, J.; Li, Z.; Gallais, J.; Doucet, D.; Arif, B.; Nakai, M.; Herniou, E.A. New insights into the evolution of Entomopoxvirinae from the complete genome sequences of four entomopoxviruses infecting Adoxophyes honmai, Choristoneura biennis, Choristoneura rosaceana, and Mythimna separata. J. Virol. 2013, 87, 7992–8003. [Google Scholar] [CrossRef] [PubMed]
  15. Forterre, P.; Prangishvili, D. The origin of viruses. Res. Microbiol. 2009, 160, 466–472. [Google Scholar] [CrossRef] [PubMed]
  16. Upton, C.; Slack, S.; Hunter, A.L.; Ehlers, A.; Roper, R.L. Poxvirus orthologous clusters: Toward defining the minimum essential poxvirus genome. J. Virol. 2003, 77, 7590–7600. [Google Scholar] [CrossRef] [PubMed]
  17. Garavaglia, M.J.; Miele, S.A.B.; Iserte, J.A.; Belaich, M.N.; Ghiringhelli, P.D. The ac53, ac78, ac101, and ac103 genes are newly discovered core genes in the family Baculoviridae. J. Virol. 2012, 86, 12069–12079. [Google Scholar] [CrossRef] [PubMed]
  18. De Andrade Zanotto, P.M.; Krakauer, D.C. Complete genome viral phylogenies suggests the concerted evolution of regulatory cores and accessory satellites. PLoS ONE 2008, 3, e3500. [Google Scholar] [CrossRef] [PubMed]
  19. Dall, D.; Luque, T.; O’Reilly, D. Insect-virus relationships: Sifting by informatics. Bioessays 2001, 23, 184–193. [Google Scholar] [PubMed]
  20. Hughes, A.L.; Friedman, R. Genome-wide survey for genes horizontally transferred from cellular organisms to baculoviruses. Mol. Biol. Evol. 2003, 20, 979–987. [Google Scholar] [CrossRef] [PubMed]
  21. McLysaght, A.; Baldi, P.F.; Gaut, B.S. Extensive gene gain associated with adaptive evolution of poxviruses. Proc. Natl. Acad. Sci. USA 2003, 100, 15655–15660. [Google Scholar] [CrossRef] [PubMed]
  22. Hughes, A.L.; Friedman, R. Poxvirus genome evolution by gene gain and loss. Mol. Phylogenet. Evol. 2005, 35, 186–195. [Google Scholar] [CrossRef] [PubMed]
  23. Bratke, K.A.; McLysaght, A. Identification of multiple independent horizontal gene transfers into poxviruses using a comparative genomics approach. BMC Evol. Biol. 2008, 8, e67. [Google Scholar] [CrossRef]
  24. Insect Virology; Asgari, S.; Johnson, K.N. (Eds.) Caister Academic Press: Norfolk, UK, 2010.
  25. Eddy, S.R. Profile hidden Markov models. Bioinformatics 1998, 14, 755–763. [Google Scholar] [CrossRef] [PubMed]
  26. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; Lopez, R.; McWilliam, H.; Remmert, M.; Söding, J.; et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol. Syst. Biol. 2011, 7, e539. [Google Scholar] [CrossRef]
  27. Hyatt, D.; Chen, G.-L.; Locascio, P.F.; Land, M.L.; Larimer, F.W.; Hauser, L.J. Prodigal: Prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 2010, 11, 119. [Google Scholar] [CrossRef]
  28. Stamatakis, A. RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 2006, 22, 2688–2690. [Google Scholar] [CrossRef] [PubMed]
  29. Mukawa, S.; Goto, C. In vivo characterization of two granuloviruses in larvae of Mythimna separata (Lepidoptera: Noctuidae). J. Gen. Virol. 2008, 89, 915–921. [Google Scholar] [CrossRef] [PubMed]
  30. Boc, A.; Philippe, H.; Makarenkov, V. Inferring and validating horizontal gene transfer events using bipartition dissimilarity. Syst. Biol. 2010, 59, 195–211. [Google Scholar] [CrossRef] [PubMed]
  31. Thézé, J.; Bézier, A.; Periquet, G.; Drezen, J.-M.; Herniou, E.A. Paleozoic origin of insect large dsDNA viruses. Proc. Natl. Acad. Sci. USA 2011, 108, 15931–15935. [Google Scholar] [CrossRef] [PubMed]
  32. Brüssow, H. The not so universal tree of life or the place of viruses in the living world. Phil. Trans. R. Soc. B 2009, 364, 2263–2274. [Google Scholar] [CrossRef] [PubMed]
  33. Filée, J.; Pouget, N.; Chandler, M. Phylogenetic evidence for extensive lateral acquisition of cellular genes by Nucleocytoplasmic large DNA viruses. BMC Evol. Biol. 2008, 8, e320. [Google Scholar] [CrossRef]
  34. Possee, R.D.; Rohrmann, G.F. Baculovirus genome organization and evolution. In The Baculoviruses (The Viruses); Miller, L.K., Ed.; Springer: Berlin, Germany, 1997; pp. 109–134. [Google Scholar]
  35. Xu, J.; Hukuhara, T. Enhanced infection of a nuclear polyhedrosis virus in larvae of the armyworm, Pseudaletia separata, by a factor in the spheroids of an entomopoxvirus. J. Invertebr. Pathol. 1992, 60, 259–264. [Google Scholar] [CrossRef]
  36. Mitsuhashi, W.; Kawakita, H.; Murakami, R.; Takemoto, Y.; Saiki, T.; Miyamoto, K.; Wada, S. Spindles of an entomopoxvirus facilitate its infection of the host insect by disrupting the peritrophic membrane. J. Virol. 2007, 81, 4235–4243. [Google Scholar] [CrossRef] [PubMed]
  37. Takatsuka, J.; Okuno, S.; Ishii, T.; Nakai, M.; Kunimi, Y. Fitness-related traits of entomopoxviruses isolated from Adoxophyes honmai (Lepidoptera: Tortricidae) at three localities in Japan. J. Invertebr. Pathol. 2010, 105, 121–131. [Google Scholar] [CrossRef] [PubMed]
  38. Nakai, M.; Shiotsuki, T.; Kunimi, Y. An entomopoxvirus and a granulovirus use different mechanisms to prevent pupation of Adoxophyes honmai. Virus Res. 2004, 101, 185–191. [Google Scholar] [CrossRef] [PubMed]
  39. Katzourakis, A.; Gifford, R.J. Endogenous viral elements in animal genomes. PLoS Genet. 2010, 6, e1001191. [Google Scholar] [CrossRef] [PubMed]
  40. Dunning Hotopp, J.C.; Clark, M.E.; Oliveira, D.C.S.G.; Foster, J.M.; Fischer, P.; Muñoz Torres, M.C.; Giebel, J.D.; Kumar, N.; Ishmael, N.; Wang, S.; et al. Widespread lateral gene transfer from intracellular bacteria to multicellular eukaryotes. Science 2007, 317, 1753–1756. [Google Scholar] [CrossRef] [PubMed]
  41. Sharon, I.; Alperovitch, A.; Rohwer, F.; Haynes, M.; Glaser, F.; Atamna-Ismaeel, N.; Pinter, R.Y.; Partensky, F.; Koonin, E.V.; Wolf, Y.I.; et al. Photosystem I gene cassettes are present in marine virus genomes. Nature 2009, 461, 258–262. [Google Scholar] [CrossRef] [PubMed]
  42. Zilber-Rosenberg, I.; Rosenberg, E. Role of microorganisms in the evolution of animals and plants: The hologenome theory of evolution. FEMS Microbiol. Rev. 2008, 32, 723–735. [Google Scholar] [CrossRef] [PubMed]
  43. Gilbert, C.; Chateigner, A.; Ernenwein, L.; Barbe, V.; Bézier, A.; Herniou, E.A.; Cordaux, R. Population genomics supports baculoviruses as vectors of horizontal transfer of insect transposons. Nat. Commun. 2014, 5, 3348. [Google Scholar] [PubMed]
  44. Pace, J.K.; Gilbert, C.; Clark, M.S.; Feschotte, C. Repeated horizontal transfer of a DNA transposon in mammals and other tetrapods. Proc. Natl. Acad. Sci. USA 2008, 105, 17023–17028. [Google Scholar] [CrossRef] [PubMed]
  45. Feschotte, C.; Pritham, E.J. A cornucopia of Helitrons shapes the maize genome. Proc. Natl. Acad. Sci. USA 2009, 106, 19747–19748. [Google Scholar] [CrossRef] [PubMed]
  46. Breitbart, M.; Rohwer, F. Here a virus, there a virus, everywhere the same virus? Trends Microbiol. 2005, 13, 278–284. [Google Scholar] [CrossRef] [PubMed]
  47. Elde, N.C.; Child, S.J.; Eickbush, M.T.; Kitzman, J.O.; Rogers, K.S.; Shendure, J.; Geballe, A.P.; Malik, H.S. Poxviruses deploy genomic accordions to adapt rapidly against host antiviral defenses. Cell 2012, 150, 831–841. [Google Scholar] [CrossRef] [PubMed]
  48. Bergin, D.; Reeves, E.P.; Renwick, J.; Wientjes, F.B.; Kavanagh, K. Superoxide production in Galleria mellonella hemocytes: Identification of proteins homologous to the NADPH oxidase complex of human neutrophils. Infect. Immun. 2005, 73, 4161–4170. [Google Scholar] [CrossRef] [PubMed]
  49. Clem, R.J. Baculoviruses and apoptosis: A diversity of genes and responses. Curr. Drug Targets 2007, 8, 1069–1074. [Google Scholar] [CrossRef] [PubMed]
  50. Hughes, A.L. Evolution of inhibitors of apoptosis in baculoviruses and their insect hosts. Infect. Genet. Evol. 2002, 2, 3–10. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Thézé, J.; Takatsuka, J.; Nakai, M.; Arif, B.; Herniou, E.A. Gene Acquisition Convergence between Entomopoxviruses and Baculoviruses. Viruses 2015, 7, 1960-1974. https://0-doi-org.brum.beds.ac.uk/10.3390/v7041960

AMA Style

Thézé J, Takatsuka J, Nakai M, Arif B, Herniou EA. Gene Acquisition Convergence between Entomopoxviruses and Baculoviruses. Viruses. 2015; 7(4):1960-1974. https://0-doi-org.brum.beds.ac.uk/10.3390/v7041960

Chicago/Turabian Style

Thézé, Julien, Jun Takatsuka, Madoka Nakai, Basil Arif, and Elisabeth A. Herniou. 2015. "Gene Acquisition Convergence between Entomopoxviruses and Baculoviruses" Viruses 7, no. 4: 1960-1974. https://0-doi-org.brum.beds.ac.uk/10.3390/v7041960

Article Metrics

Back to TopTop