Ancient coevolution of baculoviruses and their insect hosts

Molecular antiviral responses, immune priming and inheritance in insects

Viral diseases transmitted by insects to plants cause severe agricultural damage and arboviruses transmitted to humans cause severe disease outbreaks. The interaction between viruses and the insect defences is complex and has evolved into acting-counteracting molecular interplays. Viruses depict complex molecular mechanisms to ensure invasion, replication and exit the insect host cell, to invade other cells. On the other hand, insect cells use molecular strategies to recognize, halt replication and eliminate the invaders. In turn, virus counteract with evasive strategies. The main antiviral defence mechanism RNA interference (RNAi) recognizes and degrades viral RNA, thereby inhibiting viral replication. These in conjunction with other canonical immune pathways, Toll, IMD, JAK/STAT and Akt-ERK developed mainly to combat bacteria, fungi and protozoa, along with mechanisms to eliminate infected cells like apoptosis and phagocytosis comprise a multifactorial system. Insects exposed to an attenuated or sublethal viral infection could respond with faster and enhanced immune responses to the same pathogen (priming), which is like immunological memory in vertebrates. Several mechanisms have been proposed to explain priming, including endoreplication, epigenetic gene modifications by DNA methylation and histone acetylation. Priming could be inherited by the offspring (transgenerational immune priming, TGIP). However, the precise molecular mechanisms underlying TGIP remain to be elucidated. This article reviews the molecular mechanisms employed by insects to combat viral infections, discusses the current information and the outstanding research questions in the area.

Metabolites-induced co-evolutionary warfare between plants, viruses, and their associated vectors: So close yet so far away

Agriculture and global food security encounter significant challenges due to viral threats. In the following decades, several molecular studies have focused on discovering biosynthetic pathways of numerous defensive and signaling compounds, as key regulators of plant interactions, either with viruses or their associated vectors. Nevertheless, the complexities of specialized metabolites mediated plant-virus-vector tripartite viewpoint and the identification of their co-evolutionary crossroads toward antiviral defense system, remain elusive. The current study reviews the various roles of plant-specialized metabolites (PSMs) and how plants use these metabolites to defend against viruses. It discusses recent examples of specialized metabolites that have broad-spectrum antiviral properties. Additionally, the study presents the co-evolutionary basis of metabolite-mediated plant-virus-insect interactions as a potential bioinspired approach to combat viral threats. The prospects also show promising metabolic engineering strategies aimed at discovering a wide range of PSMs that are effective in fending off viruses and their related vectors. These advances in understanding the potential role of PSMs in plant-virus interactions not only serve as a cornerstone for developing plant antiviral systems, but also highlight essential principles of biological control.

Leucoma salicis nucleopolyhedrovirus (LesaNPV) genome sequence shed new light on the origin of the Alphabaculovirus orpseudotsugatae species

LesaNPV (Leucoma salicis nucleopolyhedrovirus) is an alphabaculovirus group Ib. Potentially, it can be an eco-friendly agent to control the white satin moth Leucoma salicis population. In this study, we have established the relationship between LesaNPV and other closely related alphabaculoviruses. Environmental samples of late instar of white satin moth collected in Poland infected with baculovirus have been homogenized, polyhedra were purified and subjected to scanning and transmission electron microscopy. Viral DNA was sequenced using the Illumina platform and the whole-genome sequence was established by de novo assembly of paired reads. Genome annotation and phylogenetic analyses were performed with the use of bioinformatics tools. The genome of LesaNPV is 132 549 bp long with 154 ORFs and 54.9% GC content. Whole-genome sequencing revealed deletion of dUTPase as well as ribonucleoside reductases small and large subunits region in LesaNPV genome compared to Dasychira pudibunda nucleopolyhedrovirus (DapuNPV) and Orgyia pseudotsugata multiple nucleopolyhedrovirus (OpMNPV) where this region is complete. Phylogenetic analysis of Baculoviridae family members showed that LesaNPV is less divergent from a common ancestor than closely related species DapuNPV and OpMNPV. This is interesting because their hosts do not occur in the same area. The baculoviruses described in this manuscript are probably isolates of one species and could be assigned to recently denominated species Alphabaculovirus orpseudotsugatae, historically originating from OpMNPV. This finding could have significant implications for the classification and understanding of the phylogeographical spread of baculoviruses.

Virus-host coevolutionary analyses of an Alphabaculovirus with a wide host range

Baculoviruses are highly host specific, and their host range is usually restricted to a single or a few closely related insect species, except for few virus species, e.g. Alphabaculovirus aucalifonicae and Alphabaculovirus mabrassicae. In this study, two new alphabaculovirus isolates were isolated from the larvae of Mamestra brassicae and Mythimna separata, which were named as Mamestra brassicae multiple nucleopolyhedrovirus isolate QD (MbMNPV-QD) and Mythimna separata multiple nucleopolyhedrovirus isolate Hb (MyseMNPV-Hb), respectively. The Kimura two-parameter values based on the concatenated 38 core genes of baculovirus revealed that MbMNPV (isolates QD/CHb1/K1/CTa), MyseMNPV-Hb, Helicoverpa armigera multiple nucleopolyhedrovirus (HearMNPV) and Mamestra configurata nucleopolyhedrovirus B (MacoNPV-B) were different isolates of a same virus species. A phylogenetic tree of baculoviruses and nudiviruses constructed from their 20 homologous gene sequences, and that of their isolated hosts constructed from 13 protein-coding genes of the insect mitochondrial genomes, were used to analyse the coevolution of baculoviruses with their isolated hosts. The results showed that M. brassicae was the most likely ancestral host of these virus isolates, included MbMNPV isolates, MyseMNPV-Hb, HearMNPV, and MacoNPV-B. Therefore, we concluded that these virus isolates belong to the existing virus species - Alphabaculovirus mabrassicae with M. brassicae as their ancestral host.

Complete genome sequence analysis and genome organization of Dasychira pudibunda nucleopolyhedrovirus (DapuNPV-T1) from Turkey

The economically important pale tussock moth Dasychira pudibunda L. (Lepidoptera: Lymantriidae), known as a beech pest in Europe, primarily inflicts damage on beech leaves. In the present study, we aim to reveal the genome characteristics of Dasychira pudibunda nucleopolyhedrovirus (DapuNPV-T1), which was detected for the first time in Turkey and compare it with the reference genome and other baculovirus genomes. The DapuNPV-T1 genome was determined to be a circular, double-stranded DNA molecule with 136,920 bp and a nucleotide distribution of 54.4% G + C. Bioinformatics analysis showed that the DapuNPV-T1 genome contains 163 open reading frames with more than 150 nucleotides. Fifty-four ORFs of unknown function, 6 homologous regions (hrs), 1 AC-rich region, and 3 bro genes (bro-a, bro-b, and bro-c) were determined in the genome sequence. Comparative analysis with other baculovirus strains revealed distinctions in the DapuNPV-T1 genome based on ORF. The gene parity plot and phylogenetic analysis confirmed that DapuNPV-T1 belongs to the alphabaculovirus group Ib. In addition, the DapuNPV-T1 isolate was found to be close to the nucleopolyhedrovirus Poland isolate in Dasychira pudibunda and Orgyia pseudotsugata multiple nucleopolyhedrovirus. With this study, the first genome analysis of DapuNPV from Turkey became the second in the world to enter the literature. Comprehensive information on a wide range of isolates will provide a more detailed overview of baculoviruses and help overcome their shortcomings as biocontrol agents.

Parasitoid Serpins Evolve Novel Functions to Manipulate Host Homeostasis

Parasitoids introduce various virulence factors when parasitism occurs, and some taxa generate teratocytes to manipulate the host immune system and metabolic homeostasis for the survival and development of their progeny. Host-parasitoid interactions are extremely diverse and complex, yet the evolutionary dynamics are still poorly understood. A category of serpin genes, named CvT-serpins, was discovered to be specifically expressed and secreted by the teratocytes of Cotesia vestalis, an endoparasitoid of the diamondback moth Plutella xylostella. Genomic and phylogenetic analysis indicated that the C. vestalis serpin genes are duplicated and most of them are clustered into 1 monophyletic clade. Intense positive selection was detected at the residues around the P1-P1' cleavage sites of the Cv-serpin reactive center loop domain. Functional analyses revealed that, in addition to the conserved function of melanization inhibition (CvT-serpins 1, 16, 18, and 21), CvT-serpins exhibited novel functions, i.e. bacteriostasis (CvT-serpins 3 and 5) and nutrient metabolism regulation (CvT-serpins 8 and 10). When the host-parasitoid system is challenged with foreign bacteria, CvT-serpins act as an immune regulator to reprogram the host immune system through sustained inhibition of host melanization while simultaneously functioning as immune effectors to compensate for this suppression. In addition, we provided evidence that CvT-serpin8 and 10 participate in the regulation of host trehalose and lipid levels by affecting genes involved in these metabolic pathways. These findings illustrate an exquisite tactic by which parasitoids win out in the parasite-host evolutionary arms race by manipulating host immune and nutrition homeostasis via adaptive gene evolution and neofunctionalization.

MicroRNA-Mediated Host Immune Genes Manipulation Benefits AcMNPV Proliferation in Spodoptera frugiperda

Spodoptera frugiperda is a highly destructive migratory pest that threatens various crops globally. Autographa californica multiple nucleopolyhedrovirus (AcMNPV) is an effective biocontrol agent against lepidopteran pests. Here, we explored the molecular mechanisms underlying the immune response to AcMNPV infection in S. frugiperda. RNA-seq and quantitative reverse transcription-polymerase chain reaction (qRT-PCR) analyses identified the Toll, IMD, and apoptosis pathways as primary immune responses. Investigation into AcMNPV-induced apoptosis in the S. frugiperda cell line (Sf9) revealed that the Toll pathway activated the JNK via the TRAF6 (TNF receptor-associated factor 6) adapter. In addition, AcMNPV-induced the differential expression of several host-encoded microRNAs (miRNAs), with significant negative regulatory effects, on S. frugiperda antiviral immune genes. RNAi and miRNA-mimic mediated silencing of these genes resulted in increased AcMNPV proliferation. Our findings reinforce the potential of AcMNPV as a potent biocontrol agent and further our understanding of developing biotechnology-based targeted pest control agents.

Antiviral Treatment Reveals a Cooperative Pathogenicity of Baculovirus and Iflavirus in Spodoptera exigua, a Lepidopteran Insect

NPVThe beet armyworm, Spodoptera exigua, is a serious insect pest infesting various vegetable crops. Two infectious insect viruses, baculovirus and iflavirus, are known to induce epizootics in S. exigua populations. Indeed, some laboratory colonies have appeared to be covertly infected by these viruses. Diagnostic PCR tests detected two different viruses: Spodoptera exigua multiple nucleopolyhedrosis virus (SeMNPV) and iflaviruses (SeIfV1 and SeIfV2). Viral extract from dead larvae of S. exigua could infect Sf9 cells and produce occlusion bodies (OBs). Feeding OBs to asymptomatic larvae of S. exigua caused significant viral disease. Interestingly, both SeIfV1 and SeIfV2 increased their titers at late larval stages. Sterilization of laid eggs with 1% sodium hypochloride significantly reduced SeMNPV titers and increased larval survival rate. Doublestranded RNA (dsRNA) specific to SeIfV1 or SeIfV2 significantly reduced viral titers and increased larval survival rate. To continuously feed dsRNA, a recombinant Escherichia coli HT115 expressing SeIfV1-dsRNA was constructed with an L4440 expression vector. Adding this recombinant E. coli to the artificial diet significantly reduced the SeIfV1 titer and increased larval survival. These results indicate that laboratory colony collapse of S. exigua is induced by multiple viral infections. In addition, either suppression of SeMNPV or SeIfV infection significantly increased larval survival, suggesting a cooperative pathogenicity between baculovirus and iflavirus against S. exigua.

Mosaic genome evolution and phylogenetics of Chrysodeixis includens nucleopolyhedrovirus (ChinNPV) and virulence of seven new isolates from the Brazilian states of Minas Gerais and Mato Grosso

In a comparative analysis of genome sequences from isolates of the baculovirus Chrysodeixis includens nucleopolyhedrovirus (ChinNPV) from Brazil and Guatemala, we identified a subset of isolates possessing chimeric genomes. We identified six distinct phylogenetically incongruous regions (PIRs) dispersed in the genomes, of between 279 and 3345 bp in length. The individual PIRs possessed high sequence similarity among the affected ChinNPV isolates but varied in coverage in some instances. The donor for four of the PIRs implicated in horizontal gene transfer (HGT) was identified as Trichoplusia ni single nucleopolyhedrovirus (TnSNPV), an alphabaculovirus closely related to ChinNPV, or another unknown but closely related virus. BLAST searches of the other two PIRs returned only ChinNPV sequences, but HGT from an unknown donor baculovirus cannot be excluded. Although Chrysodeixis includens and Trichoplusia ni are frequently co-collected from soybean fields in Brazil, pathogenicity data suggest that natural coinfection of C. includens larvae with ChinNPV and TnSNPV is probably uncommon. Additionally, since the chimeric ChinNPV genomes with tracts of TnSNPV sequence were restricted to a single monophyletic lineage of closely related isolates, a model of progressive restoration of the native DNA sequence by recombination with ChinNPV possessing a fully or partially non-chimeric genome is reasonable. However, multiple independent HGT from TnSNPV to ChinNPV during the evolution of these isolates cannot be excluded. Mortality data suggest that the ChinNPV isolates with chimeric genomes are not significantly different in pathogenicity towards C. includens when compared to most other ChinNPV isolates. Exclusion of the PIRs prior to phylogenetic analysis had a large impact on the topology of part of the maximum-likelihood tree, revealing a homogenous clade of three isolates (IB, IC and ID) from Paraná state in Brazil collected in 2006, together with an isolate from Guatemala collected in 1972 (IA), comprising the lineage uniquely affected by HGT from TnSNPV. The other 10 Brazilian ChinNPV isolates from Paraná, Mato Grosso, and Minas Gerais states showed higher variability, where only three isolates from Paraná state formed a monophyletic group correlating with geographical origin.

Genome Analysis of a Novel Clade b Betabaculovirus Isolated from the Legume Pest Matsumuraeses phaseoli (Lepidoptera: Tortricidae)

Matsumuraeses phaseoli is a Lepidopteran pest that primarily feeds on numerous species of cultivated legumes, such as Glycine and Phaseolus. It is widely distributed in northeast Asia. A novel granulovirus, designated as Matsumuraeses phaseoli granulovirus (MaphGV), was isolated from pathogenic M. phaseoli larvae that dwell in rolled leaves of Astragalus membranaceus, a Chinese medicinal herb. In this study, using next-generation sequencing, we report the complete genome of MaphGV. MaphGV genome comprises a double-stranded DNA of 116,875 bp, with 37.18% GC content. It has 128 hypothetical open reading frames (ORFs). Among them, 38 are baculovirus core genes, 18 are lepidopteran baculovirus conserved genes, and 5 are unique to Baculoviridae. MaphGV has one baculovirus repeat ORF (bro) and three inhibitors of apoptosis proteins (iap), including a newfound iap-6. We found two atypical baculoviral homologous regions (hrs) and four direct repeats (drs) in the MaphGV genome. Based on phylogenetic analysis, MaphGV belongs to Clade b of Betabaculovirus and is closely related to Cydia pomonellagranulovirus (CpGV) and Cryptophlebia leucotretagranulovirus (CrleGV). This novel baculovirus discovery and sequencing are invaluable in understanding the evolution of baculovirus and MaphGV may be a potential biocontrol agent against the bean ravaging pest.

Understanding the Evolutionary Ecology of host--pathogen Interactions Provides Insights into the Outcomes of Insect Pest Biocontrol

The use of viral pathogens to control thepopulation size of pest insects has produced both successful and unsuccessful outcomes. Here, we investigate whether those biocontrol successes and failures can be explained by key ecological and evolutionary processes between hosts and pathogens. Specifically, we examine how heterogeneity inpathogen transmission, ecological and evolutionary tradeoffs, andpathogen diversity affect insect population density and thus successful control. Wefirst review theexisting literature and then use numerical simulations of mathematical models to further explore these processes. Our results show that thecontrol of insect densities using viruses depends strongly on theheterogeneity of virus transmission among insects. Overall, increased heterogeneity of transmission reduces theeffect of viruses on insect densities and increases thelong-term stability of insect populations. Lower equilibrium insect densities occur when transmission is heritable and when there is atradeoff between mean transmission and insect fecundity compared to when theheterogeneity of transmission arises from non-genetic sources. Thus, theheterogeneity of transmission is akey parameter that regulates thelong-term population dynamics of insects and their pathogens. Wealso show that both heterogeneity of transmission and life-history tradeoffs modulate characteristics of population dynamics such as thefrequency and intensity of ``boom--bust" population cycles. Furthermore, we show that because of life-history tradeoffs affecting thetransmission rate, theuse of multiple pathogen strains is more effective than theuse of asingle strain to control insect densities only when thepathogen strains differ considerably intheir transmission characteristics. By quantifying theeffects of ecology and evolution on population densities, we are able to offer recommendations to assess thelong-term effects of classical biocontrol.

Characterization of a novel Helitron family in insect genomes: insights into classification, evolution and horizontal transfer

Background

Helitrons play an important role in shaping eukaryotic genomes due to their ability to transfer horizontally between distantly related species and capture gene fragments during the transposition. However, the mechanisms of horizontal transfer (HT) and the process of gene fragment capturing of Helitrons still remain to be further clarified.

Results

Here, we characterized a novel Helitron family discontinuously distributed in 27 out of 256 insect genomes. The most prominent characteristic of Hel1 family is its high sequence similarity among species of different insect orders. Related elements were also identified in two spiders, representing the first report of spider Helitrons. All these elements were classified into 2 families, 9 subfamilies and 35 exemplars based on our new classification criteria. Autonomous partners of Helitron were reconstructed in the genomes of three insects and one spider. Integration pattern analysis showed that majority of Hel1A elements in Papilio xuthus and Pieris rapae inserted into introns. Consistent with filler DNA model, stepwise sequence acquisition was observed in Sfru_Hel1Aa, Sfru_Hel1Ab and Sfru_Hel1Ac in Spodoptera frugiperda. Remarkably, the evidence that Prap_Hel1Aa in a Lepdidoptera insect, Pieris rapae, was derived from Cves_Hel1Aa in a parasitoid wasp, Cotesia vestalis, suggested the role of nonregular host-parasite interactions in HT of Helitrons.

Conclusions

We proposed a modified classification criteria of Helitrons based on the important role of the 5′-end of Helitrons in transposition, and provided evidence for stepwise sequence acquisition and recurrent HT of a novel Helitron family. Our findings of the nonregular host-parasite interactions may be more conducive to the HT of transposons.

Electronic supplementary material

The online version of this article (10.1186/s13100-019-0165-4) contains supplementary material, which is available to authorized users.

The complete genome of Rachiplusia nu nucleopolyhedrovirus (RanuNPV) and the identification of a baculoviral CPD-photolyase homolog

We described a novel baculovirus isolated from the polyphagous insect pest Rachiplusia nu. The virus presented pyramidal-shaped occlusion bodies (OBs) with singly-embed nucleocapsids and a dose mortality response of 6.9 × 10³ OBs/ml to third-instar larvae of R. nu. The virus genome is 128,587 bp long with a G + C content of 37.9% and 134 predicted ORFs. The virus is an alphabaculovirus closely related to Trichoplusia ni single nucleopolyhedrovirus, Chrysodeixis chalcites nucleopolyhedrovirus, and Chrysodeixis includens single nucleopolyhedrovirus and may constitute a new species. Surprisingly, we found co-evolution among the related viruses and their hosts at species level. Besides, auxiliary genes with homologs in other baculoviruses were found, e.g. a CPD-photolyase. The gene seemed to be result of a single event of horizontal transfer from lepidopterans to alphabaculovirus, followed by a transference from alpha to betabaculovirus. The predicted protein appears to be an active enzyme that ensures likely DNA protection from sunlight.

New Method for Differentiation of Granuloviruses (Betabaculoviruses) Based on Real-Time Polymerase Chain Reaction (Real-Time PCR)

Baculoviridae is a highly diverse family of rod-shaped viruses with double-stranded DNA. To date, almost 100 species have had their complete genomic sequences deposited in the GenBank database, a quarter of which comprises granuloviruses (GVs). Many of the genomes are sequenced using next-generation sequencing, which is currently considered the best method for characterizing new species, but it is time-consuming and expensive. Baculoviruses form a safe alternative to overused chemical pesticides and therefore there is a constant need for identifying new species that can be active components of novel biological insecticides. In this study, we have described a fast and reliable method for the detection of new and differentiation of previously analyzed granulovirus species based on a real-time polymerase chain reaction (PCR) technique with melting point curve analysis. The sequences of highly conserved baculovirus genes, such as granulin and late expression factors 8 and 9 (lef-8 and lef-9), derived from GVs available to date have been analyzed and used for degenerate primer design. The developed method was tested on a representative group of eight betabaculoviruses with comparisons of melting temperatures to allow for quick and preliminary granulovirus detection. The proposed real-time PCR procedure may be a very useful tool as an easily accessible screening method in a majority of laboratories.

The Autographa californica Multiple Nucleopolyhedrovirus ac51 Gene Is Required for Efficient Nuclear Egress of Nucleocapsids and Is Essential for In Vivo Virulence

Alphabaculoviruses are lepidopteran-specific nucleopolyhedroviruses that replicate within the nucleus; however, the anterograde transport of the nucleocapsids of these viruses, which is an obligatory step for progeny virion production, is not well understood. In the present study, a unique Alphabaculovirus gene with unknown function, namely, the Autographa californica multiple nucleopolyhedrovirus (AcMNPV) ac51 gene, was found to be required for efficient nuclear egress of AcMNPV nucleocapsids. Our results indicate that ac51 is a late gene, and Ac51 protein was detectable from 24 to 72 h postinfection using an antibody raised against Ac51. Ac51 is distributed in both the cytoplasm and nuclei of infected cells. Upon ac51 deletion, budded virion (BV) production by 96 h posttransfection was reduced by approximately 1,000-fold compared with that of wild-type AcMNPV. Neither viral DNA synthesis nor viral gene expression was affected. Ac51 was demonstrated to be a nucleocapsid protein of BVs, and ac51 deletion did not interrupt nucleocapsid assembly and occlusion-derived virion (ODV) formation. However, BV production in the supernatants of transfected cells during a viral life cycle was substantially decreased when ac51 was deleted. Further analysis showed that, compared with wild-type AcMNPV, ac51 deletion decreased nucleocapsid egress, while the numbers of nucleocapsids in the nuclei were comparable. Deletion of ac51 also eliminated the virulence of AcMNPV in vivo Taken together, our results support the conclusion that ac51 plays an important role in the nuclear egress of nucleocapsids during BV formation and is essential for the in vivo virulence of AcMNPV.

Biodiversity, Evolution and Ecological Specialization of Baculoviruses: A Treasure Trove for Future Applied Research

The Baculoviridae, a family of insect-specific large DNA viruses, is widely used in both biotechnology and biological control. Its applied value stems from millions of years of evolution influenced by interactions with their hosts and the environment. To understand how ecological interactions have shaped baculovirus diversification, we reconstructed a robust molecular phylogeny using 217 complete genomes and ~580 isolates for which at least one of four lepidopteran core genes was available. We then used a phylogenetic-concept-based approach (mPTP) to delimit 165 baculovirus species, including 38 species derived from new genetic data. Phylogenetic optimization of ecological characters revealed a general pattern of host conservatism punctuated by occasional shifts between closely related hosts and major shifts between lepidopteran superfamilies. Moreover, we found significant phylogenetic conservatism between baculoviruses and the type of plant growth (woody or herbaceous) associated with their insect hosts. In addition, we found that colonization of new ecological niches sometimes led to viral radiation. These macroevolutionary patterns show that besides selection during the infection process, baculovirus diversification was influenced by tritrophic interactions, explained by their persistence on plants and interactions in the midgut during horizontal transmission. This complete eco-evolutionary framework highlights the potential innovations that could still be harnessed from the diversity of baculoviruses.

Analysis of synonymous codon usage bias in helicase gene from Autographa californica multiple nucleopolyhedrovirus

The helicase gene of Autographa californica multiple nucleopolyhedrovirus (AcMNPV) is not only involved in viral DNA replication, but also plays a role in viral host range. To identify the codon usage bias of helicase of AcMNPV, the codon usage bias of helicase was especially studies in AcMNPV and 41 reference strains of baculoviruses by calculating the codon adaptation index (CAI), effective number of codon (ENc), relative synonymous codon usage (RSCU), and other indices. The helicase of baculovirus is less biased (mean ENc = 50.539 > 40; mean CAI = 0.246). AcMNPV helicase has a strong bias toward the synonymous codons with G and C at the third codon position (GC3s = 53.6%). The plot of GC3s against ENc values revealed that GC compositional constraints are the main factor that determines the codon usage bias of major of helicase. Several indicators supported that the codon usage pattern of helicase is mainly subject to mutation pressure. Analysis of variation in codon usage and amino acid composition indicated AcMNPV helicase shows the significant preference for one or more postulated codons for each amino acid. A cluster analysis based on RSCU values suggested that AcMNPV is evolutionarily closer to members of group I alphabaculovirus. Comparison of the codon usage pattern among E. coli, yeast, mouse, human and AcMNPV showed that yeast is a suitable expression system for AcMNPV helicase. AcMNPV helicase shows weak codon usage bias. This study may help in elucidating the functional mechanism of AcMNPV helicase and the evolution of baculovirus helicases.

Isolation of a natural DNA virus of Drosophila melanogaster, and characterisation of host resistance and immune responses

Drosophila melanogaster has played a key role in our understanding of invertebrate immunity. However, both functional and evolutionary studies of host-virus interaction in Drosophila have been limited by a dearth of native virus isolates. In particular, despite a long history of virus research, DNA viruses of D. melanogaster have only recently been described, and none have been available for experimental study. Here we report the isolation and comprehensive characterisation of Kallithea virus, a large double-stranded DNA virus, and the first DNA virus to have been reported from wild populations of D. melanogaster. We find that Kallithea virus infection is costly for adult flies, reaching high titres in both sexes and disproportionately reducing survival in males, and movement and late fecundity in females. Using the Drosophila Genetic Reference Panel, we quantify host genetic variance for virus-induced mortality and viral titre and identify candidate host genes that may underlie this variation, including Cdc42-interacting protein 4. Using full transcriptome sequencing of infected males and females, we examine the transcriptional response of flies to Kallithea virus infection and describe differential regulation of virus-responsive genes. This work establishes Kallithea virus as a new tractable model to study the natural interaction between D. melanogaster and DNA viruses, and we hope it will serve as a basis for future studies of immune responses to DNA viruses in insects.

Deltabaculoviruses encode a functional type I budded virus envelope fusion protein

Envelope fusion proteins (F proteins) are major constituents of budded viruses (BVs) of alpha- and betabaculoviruses (Baculoviridae) and are essential for the systemic infection of insect larvae and insect cell culture. An f homologue gene is absent in gammabaculoviruses. Here we characterized the putative F-homologue (Cuni-F), encoded by (ORF) 104 of Culex nigripalpus nucleopolyhedrovirus (CuniNPV), the only deltabaculovirus member. When expressed alone, this protein seems to locate on the cell surface and is able to induce cell-cell fusion. When expressed by an alphabaculovirus (Autographa california nucleopolyhedrovirus), it was found to be incorporated into BVs. Western blot analyses detected the uncleaved Cuni-F0 and the furin-cleaved F1 forms. Treatment of infected cells with tunicamycin showed that Cuni-F contains N-glycans. Mutagenesis analysis identified the canonical furin cleavage site 126RARR129 as being responsible for the cleavage of Cuni-F in insect cells. The collective evidence suggests that CuniNPV encodes a functional F protein.

Lacanobia oleracea nucleopolyhedrovirus (LaolNPV): A new European species of alphabaculovirus with a narrow host range

During an insect sampling program in alfalfa crops near Montpellier, France in 2011, Lacanobia oleracea larvae were collected that died due to nucleopolyhedrovirus infection (LaolNPV). This virus was subjected to molecular and biological characterization. The virus was a multiple nucleocapsid NPV that showed similar restriction profiles to Mamestra configurata NPV-A (MacoNPV-A) but with significant differences. Polypeptide analysis demonstrated similar proteins in occlusion bodies and occlusion derived virions, to those observed in NPVs from Mamestra spp. Terminal sequencing revealed that the genome organization shared similarity with that of MacoNPV-A. The most homologous virus was MacoNPV-A 90/2 isolate (95.63% identity and 96.47% similarity), followed by MacoNPV-A 90/4 strain (95.37% and 96.26%), MacoNPV-B (89.21% and 93.53%) and M. brassicae MNPV (89.42% and 93.74%). Phylogenetic analysis performed with lef-8, lef-9, polh and a concatenated set of genes showed that LaolNPV and the Mamestra spp. NPVs clustered together with HaMNPV, but with a closer genetic distance to MacoNPV-A strains. The Kimura 2-parameter (K-2-P) distances of the complete genes were greater than 0.05 between LaolNPV and the MbMNPV/MacoNPV-B/HaMNPV complex, which indicates that LaolNPV is a distinct species. K-2-P distances were in the range 0.015–0.050 for comparisons of LaolNPV with MacoNPV-A strains, such that additional biological characteristics should be evaluated to determine species status. While MacoNPV-A was pathogenic to seven lepidopteran species tested, LaolNPV was only pathogenic to Chrysodeixis chalcites. Given these findings, Lacanobia oleracea nucleopolyhedrovirus should be considered as a new species in the Alphabaculovirus genus.

Live phylogeny with polytomies: Finding the most compact parsimonious trees

Construction of phylogenetic trees has traditionally focused on binary trees where all species appear on leaves, a problem for which numerous efficient solutions have been developed. Certain application domains though, such as viral evolution and transmission, paleontology, linguistics, and phylogenetic stemmatics, often require phylogeny inference that involves placing input species on ancestral tree nodes (live phylogeny), and polytomies. These requirements, despite their prevalence, lead to computationally harder algorithmic solutions and have been sparsely examined in the literature to date. In this article we prove some unique properties of most parsimonious live phylogenetic trees with polytomies, and their mapping to traditional binary phylogenetic trees. We show that our problem reduces to finding the most compact parsimonious tree for n species, and describe a novel efficient algorithm to find such trees without resorting to exhaustive enumeration of all possible tree topologies.

Baculovirus FP25K Localization: Role of the Coiled-Coil Domain

Two types of viruses are produced during the baculovirus life cycle: budded virus (BV) and occlusion-derived virus (ODV). A particular baculovirus protein, FP25K, is involved in the switch from BV to ODV production. Previously, FP25K from the model alphabaculovirus Autographa californica multiple nucleopolyhedrovirus (AcMNPV) was shown to traffic ODV envelope proteins. However, FP25K localization and the domains involved are inconclusive. Here we used a quantitative approach to study FP25K subcellular localization during infection using an AcMNPV bacmid virus that produces a functional AcMNPV FP25K-green fluorescent protein (GFP) fusion protein. During cell infection, FP25K-GFP localized primarily to the cytoplasm, particularly amorphous structures, with a small fraction being localized in the nucleus. To investigate the sequences involved in FP25K localization, an alignment of baculovirus FP25K sequences revealed that the N-terminal putative coiled-coil domain is present in all alphabaculoviruses but absent in betabaculoviruses. Structural prediction indicated a strong relatedness of AcMNPV FP25K to long interspersed element 1 (LINE-1) open reading frame 1 protein (ORF1p), which contains an N-terminal coiled-coil domain responsible for cytoplasmic retention. Point mutations and deletions of this domain lead to a change in AcMNPV FP25K localization from cytoplasmic to nuclear. The coiled-coil and C-terminal deletion viruses increased BV production. Furthermore, a betabaculovirus FP25K protein lacking this N-terminal coiled-coil domain localized predominantly to the nucleus and exhibited increased BV production. These data suggest that the acquisition of this N-terminal coiled-coil domain in FP25K is important for the evolution of alphabaculoviruses. Moreover, with the divergence of preocclusion nuclear membrane breakdown in betabaculoviruses and membrane integrity in alphabaculoviruses, this domain represents an alphabaculovirus adaptation for nuclear trafficking of occlusion-associated proteins.

IMPORTANCE

Baculovirus infection produces two forms of viruses: BV and ODV. Manufacturing of ODV involves trafficking of envelope proteins to the inner nuclear membrane, mediated partly through the FP25K protein. Since FP25K is present in alpha-, beta-, and gammabaculoviruses, it is uncertain if this trafficking function is conserved. In this study, we looked at alpha- and betabaculovirus FP25K trafficking by its localization. Alphabaculovirus FP25K localized primarily to the cytoplasm, whereas betabaculovirus FP25K localized to the nucleus. We found that an N-terminal coiled-coil domain present in all alphabaculovirus FP25K proteins, but absent in betabaculovirus FP25K, was critical for alphabaculovirus FP25K cytoplasmic localization. We believe that this represents an evolutionary process that partly led to the gain of function of this N-terminal coiled-coil domain in alphabaculovirus FP25K to aid in nuclear trafficking of occlusion-associated proteins. Due to betabaculovirus breakdown of the nuclear membrane before occlusion, this function is not needed, and the domain was lost or never acquired.

Codon usage in Alphabaculovirus and Betabaculovirus hosted by the same insect species is weak, selection dominated and exhibits no more similar patterns than expected

Mutations shape synonymous codon usage bias in certain organism genomes, while selection shapes it in others. Lepidopteran-specific Alphabaculovirus and Betabaculovirus are two large genera in the family of Baculoviridae. In this study, we analyzed the codon usage patterns in 17 baculoviruses, including 10 alphabaculoviruses and 7 betabaculoviruses, which were isolated from seven insect species, and we characterized the codon usage patterns between Alphabaculovirus and Betabaculovirus. Our results show that all the baculoviruses possessed a general weak trend of codon bias. The differences of ENc (effective number of codons) values, nucleotide contents and the impacts of nucleotide content on ENc value within alpha-/betabaculovirus pairs were independent of whether the host species are the same or different. Furthermore, the majority of amino acid sequences adopted codons unequally in all viruses, but the numbers of common preferred codons between alpha- and betabaculoviruses hosted by the same insect species were not significantly different from the differences observed between alpha- and betabaculoviruses hosted by different insect species. In addition, the amino acids that adopt the same synonymous codon composition between alpha- and betabaculoviruses hosted by the same insect species were statistically as few as those between alpha- and betabaculoviruses hosted by different insect species. Correspondence analysis revealed that no major factors resulted in the codon bias in these baculoviruses, implying multiple minor influential factors exist. Neutrality plot analysis indicated that selection pressure dominated mutations in shaping the codon usage. However, the levels of selection pressure were not significantly different among viruses hosted by the same insect species. We expect that evolution would cause the alpha- and betabaculoviruses hosted by the same insect species to share more patterns, but this effect was not observed.

A Role for the Anti-Viral Host Defense Mechanism in the Phylogenetic Divergence in Baculovirus Evolution

Although phylogenic analysis often suggests co-evolutionary relationships between viruses and host organisms, few examples have been reported at the microevolutionary level. Here, we show a possible example in which a species-specific anti-viral response may drive phylogenic divergence in insect virus evolution. Two baculoviruses, Autographa californica multiple nucleopolyhedrovirus (AcMNPV) and Bombyx mori nucleopolyhedrovirus (BmNPV), have a high degree of DNA sequence similarity, but exhibit non-overlapping host specificity. In our study of their host-range determination, we found that BmNPV replication in B. mori cells was prevented by AcMNPV-P143 (AcP143), but not BmNPV-P143 (BmP143) or a hybrid P143 protein from a host-range expanded phenotype. This suggests that AcMNPV resistance in B. mori cells depends on AcP143 recognition and that BmNPV uses BmP143 to escapes this recognition. Based on these data, we propose an insect-baculovirus co-evolution scenario in which an ancestor of silkworms exploited an AcMNPV-resistant mechanism; AcMNPV counteracted this resistance via P143 mutations, resulting in the birth of BmNPV.

Linking Virus Genomes with Host Taxonomy

Environmental genomics can describe all forms of organisms—cellular and viral—present in a community. The analysis of such eco-systems biology data relies heavily on reference databases, e.g., taxonomy or gene function databases. Reference databases of symbiosis sensu lato, although essential for the analysis of organism interaction networks, are lacking. By mining existing databases and literature, we here provide a comprehensive and manually curated database of taxonomic links between viruses and their cellular hosts.

The genome of Dasychira pudibunda nucleopolyhedrovirus (DapuNPV) reveals novel genetic connection between baculoviruses infecting moths of the Lymantriidae family

Background

DapuNPV (Dasychira pudibunda nucleopolyhedrovirus), presented in this report, belongs to Alphabaculovirus group Ib. Its full, newly sequenced genome shows close relationship to baculovirus OpMNPV isolated from douglas-fir tussock moth Orgyia pseudotsugata. Baculovirus DapuNPV is a natural limiter of pale tussock moth Dasychira pudibunda L. (syn. Calliteara pudibunda L.)(Lepidoptera, Lymantriidae), which can occur in a form of an outbreak on many species of deciduous trees and may cause significant economic losses in the forests.

Methods

Late instars dead larvae of pale tussock moth were mechanically homogenized and polyhedra were purified during series of ultracentrifugation. Viral DNA was extarcted and sequenced using Miseq Illumina platform. 294,902 paired reads were used for de novo assembling. Genome annotation, multiple allingment to others baculoviruses and phylogegentic analises were perform with the use of multiple bioinformatic tools like: Glimmer3, HMMER web server, Geneious 7 and MEGA6.

Results

The genome of DapuNPV is 136,761 bp long with AT pairs content 45.6 %. The predicted number of encoded putative open reading frames (ORFs) is 161 and six of them demonstrate low or no homology to ORFs previously found in baculoviruses. DapuNPV genome shows very high similarity to OpMNPV in a nucleotide sequence (91.1 % of identity) and gene content (150 homologous ORFs), though some major differences (e.g. lack of he65 in OpMNPV) have also been noted.

Conclusions

Similarly to other members of the Baculoviridae family, DapuNPV baculovirus possesses highly conserved core genes. Among them, there is a second copy of occluded derived virus envelope 27 protein (odv-e27), which was previously found only in a member of Alphabaculovirus group II – LyxyMNPV (Lymantria xylina MNPV). Surprisingly enough, DapuNPV and LyxyMNPV genomes share also another feature. Phylogenetic analysis of chitin binding family protein (cbpl) indicates significant similarity of those two baculoviruses from distinct evolutionary groups which infect the same hosts from Lymantriidae. The ubiquitin like family gene (ubil), which has not been described until now, is another characteristic component of DapuNPV genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1963-9) contains supplementary material, which is available to authorized users.

Comparative proteomics analysis of Spodoptera frugiperda cells during Autographa californica multiple nucleopolyhedrovirus infection

Background

Increasing evidence sugggest that in addition of balculovirus controling insect host, host cells also responds to balculovirus infection. However, compared to existing knowledge on virus gene, host cell responses are relatively poorly understood.

Methods

In this study, Spodoptera frugiperda (Sf9) cells were infected with Autographa californica multiple nucleopolyhedrovirus (AcMNPV). The protein composition and protein changes of Spodoptera frugiperda (Sf9) cells of different infection stages were analysed by isobaric tag for relative and absolute quantification (iTRAQ) techniques.

Results

A total of 4004 Sf9 proteins were identified by iTRAQ and 413 proteins were found as more than 1.5-fold changes in abundance. The 413 proteins were categorised according to GO classification for insects and were categorised into: biological process, molecular function and cellular component.

Conclusions

The determination of the protein changes in infected Sf9 cells would help to better understanding of host cell responses and facilitate better design of this virus-host cell interaction in pest insect control and other related fields.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-015-0346-9) contains supplementary material, which is available to authorized users.

Global distribution of Chelonid fibropapilloma-associated herpesvirus among clinically healthy sea turtles

Background

Fibropapillomatosis (FP) is a neoplastic disease characterized by cutaneous tumours that has been documented to infect all sea turtle species. Chelonid fibropapilloma-associated herpesvirus (CFPHV) is believed to be the aetiological agent of FP, based principally on consistent PCR-based detection of herpesvirus DNA sequences from FP tumours. We used a recently described PCR-based assay that targets 3 conserved CFPHV genes, to survey 208 green turtles (Chelonia mydas). This included both FP tumour exhibiting and clinically healthy individuals. An additional 129 globally distributed clinically healthy individual sea turtles; representing four other species were also screened.

Results

CFPHV DNA sequences were obtained from 37/37 (100%) FP exhibiting green turtles, and 45/300 (15%) clinically healthy animals spanning all five species. Although the frequency of infected individuals per turtle population varied considerably, most global populations contained at least one CFPHV positive individual, with the exception of various turtle species from the Arabian Gulf, Northern Indian Ocean and Puerto Rico.

Haplotype analysis of the different gene markers clustered the CFPHV DNA sequences for two of the markers (UL18 and UL22) in turtles from Turks and Caicos separate to all others, regardless of host species or geographic origin.

Conclusion

Presence of CFPHV DNA within globally distributed samples for all five species of sea turtle was confirmed. While 100% of the FP exhibiting green turtles yielded CFPHV sequences, surprisingly, so did 15% of the clinically healthy turtles. We hypothesize that turtle populations with zero (0%) CFPHV frequency may be attributed to possible environmental differences, diet and/or genetic resistance in these individuals. Our results provide first data on the prevalence of CFPHV among seemingly healthy turtles; a factor that may not be directly correlated to the disease incidence, but may suggest of a long-term co-evolutionary latent infection interaction between CFPHV and its turtle-host across species. Finally, computational analysis of amino acid variants within the Turks and Caicos samples suggest potential functional importance in a substitution for marker UL18 that encodes the major capsid protein gene, which potentially could explain differences in pathogenicity. Nevertheless, such a theory remains to be validated by further research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0206-z) contains supplementary material, which is available to authorized users.

Single-stranded DNA fragments of insect-specific nuclear polyhedrosis virus act as selective DNA insecticides for gypsy moth control

This paper focuses on the DNA insecticides as a novel preparation against gypsy moth (Lymantria dispar) based on DNA fragments of the anti-apoptotic gene of its nuclear polyhedrosis virus. It was found that the external application of a solution with two single-stranded DNA fragments from BIR and RING domains of LdMNPV (L.dispar multicapsid nuclear polyhedrosis virus) IAP-3 (inhibitor of apoptosis) gene induces a significantly higher mortality of gypsy moth caterpillars in comparison with the application of the control solutions. This effect does not depend on the infection of caterpillars with LdMNPV. The results also show that DNA insecticides based on LdMNPV IAP-3 gene fragments can be selective in action, and at least are not harmful to tobacco hornworm (Manduca sexta) and black cutworm (Agrotis ipsilon). Part of the gypsy moth genome cloned with the fragments of BIR and RING domains of LdMNPV IAP-3 gene as primers, has an overlap with the corresponding part of the LdMNPV IAP-3 gene and L.dispar IAP-1 mRNA for an inhibitor of apoptosis protein with the high cover by query, allows assuming that we cloned a part of gypsy moth anti-apoptosis gene. This finding gives the grounding that proposed here DNA insecticides might act through the blocking of the mechanisms involved in post transcriptional expression of insect anti-apoptosis genes. The results show the insecticidal potential of the viral genome fragments that can be used to create safe and relatively fast-acting DNA insecticides to control the quantity of gypsy moth populations, important task for forestry and agriculture.

Similarity/Dissimilarity analysis of protein sequences based on a new spectrum-like graphical representation

Sequence comparison is one of the foundations in bioinformatics, which can be used to study evolutionary relations among the sequences. In this study, a 2D spectrum-like graphical representation of protein sequences is presented based on the hydrophobicity scale of amino acids. The frequencies of amplitudes of 4-subsequences are adopted to characterize a spectrum-like graph, and a 17D vector is used as the descriptor of protein sequence. The χ(2) value of compatibility test is performed. New similarity analysis approach is illustrated on the all protein sequences, which are encoded by the mitochondrion genome of 20 different species. Finally, comparison with the ClustalW method shows the utility of our method.

Brown planthopper nudivirus DNA integrated in its host genome

The brown planthopper (BPH), Nilaparvata lugens (Hemiptera:Delphacidae), is one of the most destructive insect pests of rice crops in Asia. Nudivirus-like sequences were identified during the whole-genome sequencing of BPH. PCR examination showed that the virus sequences were present in all of the 22 BPH populations collected from East, Southeast, and South Asia. Thirty-two of the 33 nudivirus core genes were identified, including 20 homologues of baculovirus core genes. In addition, several gene clusters that were arranged collinearly with those of other nudiviruses were found in the partial virus genome. In a phylogenetic tree constructed using the supermatrix method, the original virus was grouped with other nudiviruses and was closely related to polydnavirus. Taken together, these data indicated that the virus sequences belong to a new member of the family Nudiviridae. More specifically, the virus sequences were integrated into the chromosome of its insect host during coevolution. This study is the first report of a large double-stranded circular DNA virus genome in a sap-sucking hemipteran insect.

IMPORTANCE

This is the first report of a large double-stranded DNA virus integrated genome in the planthopper, a plant sap-sucking hemipteran insect. It is an exciting addition to the evolutionary story of bracoviruses (polydnaviruses), nudiviruses, and baculoviruses. The results on the virus sequences integrated in the chromosomes of its insect host also represent a story of successful coevolution of an invertebrate virus and a plant sap-sucking insect.

Comparative analysis of error-prone replication mononucleotide repeats across baculovirus genomes

Genome replication by the baculovirus DNA polymerase often generates errors in mononucleotide repeat (MNR) sequences due to replication slippage. This results in the inactivation of genes that affects different stages of the cell infection cycle. Here we mapped these MNRs in the 59 baculovirus genomes. We found that the MNR frequencies of baculovirus genomes are different and not correlated with the genome sizes. Although the average A/T content of baculoviruses is 58.67%, the A/T MNR frequency is significantly higher than that of the G/C MNRs. Furthermore, the A7/T7 MNRs are the most frequent of those we studied. Finally, MNR frequencies in different classes of baculovirus genes, such as immediate early genes, show differences between baculovirus genomes, suggesting that the distribution and frequency of different MNRs are unique to each baculovirus species or strain. Therefore, the results of this study can help select appropriate baculoviruses for the development of biological insecticides.

Unraveling the entry mechanism of baculoviruses and its evolutionary implications

The entry of baculovirus budded virus into host cells is mediated by two distinct types of envelope fusion proteins (EFPs), GP64 and F protein. Phylogenetic analysis suggested that F proteins were ancestral baculovirus EFPs, whereas GP64 was acquired by progenitor group I alphabaculovirus more recently and may have stimulated the formation of the group I lineage. This study was designed to experimentally recapitulate a possible major step in the evolution of baculoviruses. We demonstrated that the infectivity of an F-null group II alphabaculovirus (Helicoverpa armigera nucleopolyhedrovirus [HearNPV]) can be functionally rescued by coinsertion of GP64 along with the nonfusogenic F(def) (furin site mutated HaF) from HearNPV. Interestingly, HearNPV enters cells by endocytosis and, less efficiently, by direct membrane fusion at low pH. However, this recombinant HearNPV coexpressing F(def) and GP64 mimicked group I virus not only in its EFP composition but also in its abilities to enter host cells via low-pH-triggered direct fusion pathway. Neutralization assays indicated that the nonfusogenic F proteins contribute mainly to binding to susceptible cells, while GP64 contributes to fusion. Coinsertion of GP64 with an F-like protein (Ac23) from group I virus led to efficient rescue of an F-null group II virus. In summary, these recombinant viruses and their entry modes are considered to resemble an evolutionary event of the acquisition of GP64 by an ancestral group I virus and subsequent adaptive inactivation of the original F protein. The study described here provides the first experimental evidence to support the hypothesis of the evolution of baculovirus EFPs.

Pseudoplusia includens single nucleopolyhedrovirus: genetic diversity, phylogeny and hypervariability of the pif-2 gene

The soybean looper (Pseudoplusia includens Walker, 1857) has become a major pest of soybean crops in Brazil. In order to determine the genetic diversity and phylogeny of variants of Pseudoplusia includens single nucleopolyhedrovirus (PsinSNPV-IA to -IG), partial sequences of the genes lef-8, lef-9, pif-2, phr and polh were obtained following degenerate PCR and phylogenetic trees constructed using maximum parsimony and Bayesian methods. The aligned sequences showed polymorphisms among the isolates, where the pif-2 gene was by far the most variable and is predicted to be under positive selection. Furthermore, some of the pif-2 DNA sequence mutations are predicted to result in significant amino acid substitutions, possibly leading to changes in oral infectivity of this baculovirus. Cladistic analysis revealed two closely related monophyletic groups, one containing PsinNPV isolates IB, IC and ID and another containing isolates IA, IE, IF and IG. The phylogeny of PsinSNPV in relation to 56 other baculoviruses was also determined from the concatenated partial LEF-8, LEF-9, PIF-2 and POLH/GRAN deduced amino acid sequences, using maximum-parsimony and Bayesian methods. This analysis clearly places PsinSNPV with the Group II Alphabaculovirus, where PsinSNPV is most closely related to Chrysodeixis chalcites NPV and Trichoplusia ni SNPV.

Comparative genome sequence analysis of Choristoneura occidentalis Freeman and C. rosaceana Harris (Lepidoptera: Tortricidae) alphabaculoviruses

The complete genome sequences of Choristoneura occidentalis and C. rosaceana nucleopolyhedroviruses (ChocNPV and ChroNPV, respectively) (Baculoviridae: Alphabaculovirus) were determined and compared with each other and with those of other baculoviruses, including the genome of the closely related C. fumiferana NPV (CfMNPV). The ChocNPV genome was 128,446 bp in length (1147 bp smaller than that of CfMNPV), had a G+C content of 50.1%, and contained 148 open reading frames (ORFs). In comparison, the ChroNPV genome was 129,052 bp in length, had a G+C content of 48.6% and contained 149 ORFs. ChocNPV and ChroNPV shared 144 ORFs in common, and had a 77% sequence identity with each other and 96.5% and 77.8% sequence identity, respectively, with CfMNPV. Five homologous regions (hrs), with sequence similarities to those of CfMNPV, were identified in ChocNPV, whereas the ChroNPV genome contained three hrs featuring up to 14 repeats. Both genomes encoded three inhibitors of apoptosis (IAP-1, IAP-2, and IAP-3), as reported for CfMNPV, and the ChocNPV IAP-3 gene represented the most divergent functional region of this genome relative to CfMNPV. Two ORFs were unique to ChocNPV, and four were unique to ChroNPV. ChroNPV ORF chronpv38 is a eukaryotic initiation factor 5 (eIF-5) homolog that has also been identified in the C. occidentalis granulovirus (ChocGV) and is believed to be the product of horizontal gene transfer from the host. Based on levels of sequence identity and phylogenetic analysis, both ChocNPV and ChroNPV fall within group I alphabaculoviruses, where ChocNPV appears to be more closely related to CfMNPV than does ChroNPV. Our analyses suggest that it may be appropriate to consider ChocNPV and CfMNPV as variants of the same virus species.

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

Poxviruses are nucleocytoplasmic large DNA viruses encompassing two subfamilies, the Chordopoxvirinae and the Entomopoxvirinae, infecting vertebrates and insects, respectively. While chordopoxvirus genomics have been widely studied, only two entomopoxvirus (EPV) genomes have been entirely sequenced. We report the genome sequences of four EPVs of the Betaentomopoxvirus genus infecting the Lepidoptera: Adoxophyes honmai EPV (AHEV), Choristoneura biennis EPV (CBEV), Choristoneura rosaceana EPV (CREV), and Mythimna separata EPV (MySEV). The genomes are 80% AT rich, are 228 to 307 kbp long, and contain 247 to 334 open reading frames (ORFs). Most genes are homologous to those of Amsacta moorei entomopoxvirus and encode several protein families repeated in tandem in terminal regions. Some genomes also encode proteins of unknown functions with similarity to those of other insect viruses. Comparative genomic analyses highlight a high colinearity among the lepidopteran EPV genomes and little gene order conservation with other poxvirus genomes. As with previously sequenced EPVs, the genomes include a relatively conserved central region flanked by inverted terminal repeats. Protein clustering identified 104 core EPV genes. Among betaentomopoxviruses, 148 core genes were found in relatively high synteny, pointing to low genomic diversity. Whole-genome and spheroidin gene phylogenetic analyses showed that the lepidopteran EPVs group closely in a monophyletic lineage, corroborating their affiliation with the Betaentomopoxvirus genus as well as a clear division of the EPVs according to the orders of insect hosts (Lepidoptera, Coleoptera, and Orthoptera). This suggests an ancient coevolution of EPVs with their insect hosts and the need to revise the current EPV taxonomy to separate orthopteran EPVs from the lepidopteran-specific betaentomopoxviruses so as to form a new genus.

Genomic diversity of Bombyx mori nucleopolyhedrovirus strains

Bombyx mori nucleopolyhedrovirus (BmNPV) is a baculovirus that selectively infects the domestic silkworm. In this study, six BmNPV strains were compared at the whole genome level. We found that the number of bro genes and the composition of the homologous regions (hrs) are the two primary areas of divergence within these genomes. When we compared the ORFs of these BmNPV variants, we noticed a high degree of sequence divergence in the ORFs that are not baculovirus core genes. This result is consistent with the results derived from phylogenetic trees and evolutionary pressure analyses of these ORFs, indicating that ORFs that are not core genes likely play important roles in the evolution of BmNPV strains. The evolutionary relationships of these BmNPV strains might be explained by their geographic origins or those of their hosts. In addition, the total number of hr palindromes seems to affect viral DNA replication in Bm5 cells.

Identification of a previously undescribed divergent virus from the Flaviviridae family in an outbreak of equine serum hepatitis

Theiler's disease is an acute hepatitis in horses that is associated with the administration of equine blood products; its etiologic agent has remained unknown for nearly a century. Here, we used massively parallel sequencing to explore samples from a recent Theiler's disease outbreak. Metatranscriptomic analysis of the short sequence reads identified a 10.5-kb sequence from a previously undescribed virus of the Flaviviridae family, which we designate "Theiler's disease-associated virus" (TDAV). Phylogenetic analysis clusters TDAV with GB viruses of the recently proposed Pegivirus genus, although it shares only 35.3% amino acid identity with its closest relative, GB virus D. An epidemiological survey of additional horses from three separate locations supports an association between TDAV infection and acute serum hepatitis. Experimental inoculation of horses with TDAV-positive plasma provides evidence that several weeks of viremia preceded liver injury and that liver disease may not be directly related to the level of viremia. Like hepatitis C virus, the best characterized Flaviviridae species known to cause hepatitis, we find TDAV is capable of efficient parenteral transmission, engendering acute and chronic infections associated with a diversity of clinical presentations ranging from subclinical infection to clinical hepatitis.

Analysis of the genomic sequence of Philosamia cynthia nucleopolyhedrin virus and comparison with Antheraea pernyi nucleopolyhedrin virus

Background

Two species of wild silkworms, the Chinese oak silkworm (Antheraea pernyi) and the castor silkworm Philosamia cynthia ricini, can acquire a serious disease caused by Nucleopolyhedrin Viruses (NPVs) (known as AnpeNPV and PhcyNPV, respectively). The two viruses have similar polyhedral morphologies and their viral fragments share high sequence similarity. However, the physical maps of the viral genomes and cross-infectivity of the viruses are different. The genome sequences of two AnpeNPV isolates have been published.

Results

We sequenced and analyzed the full-length genome of PhcyNPV to compare the gene contents of the two viruses. The genome of PhcyNPV is 125, 376 bp, with a G + C content of 53.65%, and encodes 138 open reading frames (ORFs) of at least 50 amino acids (aa) (GenBank accession number: JX404026). Between PhcyNPV and AnpeMNPV-L and -Z isolates, 126 ORFs are identical, including 30 baculovirus core genes. Nine ORFs were only found in PhcyNPV. Four genes, cath, v-chi, lef 10 and lef 11, were not found in PhcyNPV. However, most of the six genes required for infectivity via the oral route were found in PhcyNPV and in the two AnpeNPV isolates, with high sequence similarities. The pif-3 gene of PhcyNPV contained 59 aa extra amino acids at the N-terminus compared with AnpeNPV.

Conclusions

Most of the genes in PhcyNPV are similar to the two AnpeNPV isolates, including the direction of expression of the ORFs. Only a few genes were missing from PhcyNPV. These data suggest that PhcyNPV and AnpeNPV might be variants of each other, and that the differences in cross-infection might be caused by gene mutations.

Protein tyrosine phosphatase-induced hyperactivity is a conserved strategy of a subset of baculoviruses to manipulate lepidopteran host behavior

Many parasites manipulate host behavior to increase the probability of transmission. To date, direct evidence for parasitic genes underlying such behavioral manipulations is scarce. Here we show that the baculovirus Autographa californica nuclear polyhedrovirus (AcMNPV) induces hyperactive behavior in Spodoptera exigua larvae at three days after infection. Furthermore, we identify the viral protein tyrosine phosphatase (ptp) gene as a key player in the induction of hyperactivity in larvae, and show that mutating the catalytic site of the encoded phosphatase enzyme prevents this induced behavior. Phylogenetic inference points at a lepidopteran origin of the ptp gene and shows that this gene is well-conserved in a group of related baculoviruses. Our study suggests that ptp-induced behavioral manipulation is an evolutionarily conserved strategy of this group of baculoviruses to enhance virus transmission, and represents an example of the extended phenotype concept. Overall, these data provide a firm base for a deeper understanding of the mechanisms behind baculovirus-induced insect behavior.

Cross-species transmission in the speciation of the currently known murinae-associated hantaviruses

To gain more insight into the phylogeny of Dabieshan virus (DBSV), carried by Niviventer confucianus and other Murinae-associated hantaviruses, genome sequences of novel variants of DBSV were recovered from Niviventer rats trapped in the mountainous areas of Wenzhou, China. Genetic analyses show that all known genetic variants of DBSV, including the ones identified in this study, are distinct from other Murinae-associated hantaviruses. DBSV variants show geographic clustering and high intraspecies diversity. The data suggest that DBSV is a distinct species in the genus Hantavirus. Interestingly, DBSV shows the highest sequence identity to Hantaan virus (HTNV), with a >7% difference in the sequences of the N, GPC, and L proteins, while N. confucianus is more closely related to Rattus norvegicus (the host of Seoul virus [SEOV]) than to Apodemus agrarius (the host of HTNV and Saaremaa virus [SAAV]). Further genetic analyses of all known Murinae-associated hantaviruses (both established and tentative species) show that many of them, including DBSV, may have originated from host switching. The estimation of evolutionary rates and divergence time supports the role of cross-species transmission in the evolution of Murinae-associated hantaviruses. The detection of positive selection suggests that genetic drift may contribute to the speciation of Murinae-associated hantaviruses and that adaptation has a role as well.

Functional and biochemical characterization of the baculovirus caspase inhibitor MaviP35

Many viruses express proteins which prevent the host cell death that their infection would otherwise provoke. Some insect viruses suppress host apoptosis through the expression of caspase inhibitors belonging to the P35 superfamily. Although a number of P35 relatives have been identified, Autographa californica (Ac) P35 and Spodoptera littoralis (Spli) P49 have been the most extensively characterized. AcP35 was found to inhibit caspases via a suicide substrate mechanism: the caspase cleaves AcP35 within its ‘reactive site loop' then becomes trapped, irreversibly bound to the cleaved inhibitor. The Maruca vitrata multiple nucleopolyhedrovirus encodes a P35 family member (MaviP35) that exhibits 81% identity to AcP35. We found that this relative shared with AcP35 the ability to inhibit mammalian and insect cell death. Caspase-mediated cleavage within the MaviP35 reactive site loop occurred at a sequence distinct from that in AcP35, and the inhibitory profiles of the two P35 relatives differed. MaviP35 potently inhibited human caspases 2 and 3, DCP-1, DRICE and CED-3 in vitro, but (in contrast to AcP35) only weakly suppressed the proteolytic activity of the initiator human caspases 8, 9 and 10. Although MaviP35 inhibited the AcP35-resistant caspase DRONC in yeast, and was sensitive to cleavage by DRONC in vitro, MaviP35 failed to inhibit the proteolytic activity of bacterially produced DRONC in vitro.

Complete sequence, analysis and organization of the Orgyia leucostigma nucleopolyhedrovirus genome

The complete genome of the Orgyia leucostigma nucleopolyhedrovirus (OrleNPV) isolated from the whitemarked tussock moth (Orgyia leucostigma, Lymantridae: Lepidoptera) was sequenced, analyzed, and compared to other baculovirus genomes. The size of the OrleNPV genome was 156,179 base pairs (bp) and had a G+C content of 39%. The genome encoded 135 putative open reading frames (ORFs), which occupied 79% of the entire genome sequence. Three inhibitor of apoptosis (ORFs 16, 43 and 63), and five baculovirus repeated ORFs (bro-a through bro-e) were interspersed in the OrleNPV genome. In addition to six direct repeat (drs), a common feature shared among most baculoviruses, OrleNPV genome contained three homologous regions (hrs) that are located in the latter half of the genome. The presence of an F-protein homologue and the results from phylogenetic analyses placed OrleNPV in the genus Alphabaculovirus, group II. Overall, OrleNPV appears to be most closely related to group II alphabaculoviruses Ectropis obliqua (EcobNPV), Apocheima cinerarium (ApciNPV), Euproctis pseudoconspersa (EupsNPV), and Clanis bilineata (ClbiNPV).