Sequence analysis of the Xestia c-nigrum granulovirus genome

Genome comparison of long-circulating field CnmeGV isolates from the same region

Cnaphalocrocis medinalis granulovirus (CnmeGV), belonging to Betabaculovirus cnamedinalis, can infect the rice pest, the rice leaf roller. In 1979, a CnmeGV isolate, CnmeGV-EP, was collected from Enping County, China. In 2014, we collected another CnmeGV isolate, CnmeGV-EPDH3, at the same location and obtained the complete virus genome sequence using Illumina and ONT sequencing technologies. By combining these two virus isolates, we updated the genome annotation of CnmeGV and conducted an in-depth analysis of its genome features. CnmeGV genome contains abundant tandem repeat sequences, and the repeating units in the homologous regions (hrs) exhibit overlapping and nested patterns. The genetic variations within EPDH3 population show the high stability of CnmeGV genome, and tandem repeats are the only region of high genetic variation in CnmeGV genome replication. Some defective viral genomes formed by recombination were found within the population. Comparison analysis of the two virus isolates collected from Enping showed that the proteins encoded by the CnmeGV-specific genes were less conserved relative to the baculovirus core genes. At the genomic level, there are a large number of SNPs and InDels between the two virus isolates, especially in and around the bro genes and hrs. Additionally, we discovered that CnmeGV acquired a segment of non-ORF sequence from its host, which does not provide any new proteins but rather serves as redundant genetic material integrated into the viral genome. Furthermore, we observed that the host's transposon piggyBac has inserted into some virus genes. Together, dsDNA viruses could acquire non-coding genetic material from their hosts to expand the size of their genomes. These findings provide new insights into the evolution of dsDNA viruses.

Helical reconstruction of VP39 reveals principles for baculovirus nucleocapsid assembly

Baculoviruses are insect-infecting pathogens with wide applications as biological pesticides, in vitro protein production vehicles and gene therapy tools. Its cylindrical nucleocapsid, which encapsulates and protects the circular double-stranded viral DNA encoding proteins for viral replication and entry, is formed by the highly conserved major capsid protein VP39. The mechanism for VP39 assembly remains unknown. We use electron cryomicroscopy to determine a 3.2 Å helical reconstruction of an infectious nucleocapsid of Autographa californica multiple nucleopolyhedrovirus, revealing how dimers of VP39 assemble into a 14-stranded helical tube. We show that VP39 comprises a distinct protein fold conserved across baculoviruses, which includes a Zinc finger domain and a stabilizing intra-dimer sling. Analysis of sample polymorphism shows that VP39 assembles in several closely-related helical geometries. This VP39 reconstruction reveals general principles for baculoviral nucleocapsid assembly.

Helical reconstruction of VP39 reveals principles for baculovirus nucleocapsid assembly

Baculoviruses are insect-infecting pathogens with wide applications as biological pesticides, in vitro protein production vehicles and gene therapy tools. Its cylindrical nucleocapsid, which encapsulates and protects the circular double-stranded viral DNA encoding proteins for viral replication and entry, is formed by the highly conserved major capsid protein VP39. The mechanism for VP39 assembly remains unknown. We determined a 3.2 Å electron cryomicroscopy helical reconstruction of an infectious nucleocapsid of Autographa californica multiple nucleopolyhedrovirus, revealing how dimers of VP39 assemble into a 14-stranded helical tube. We show that VP39 comprises a unique protein fold conserved across baculoviruses, which includes a Zinc finger domain and a stabilizing intra-dimer sling. Analysis of sample polymorphism revealed that VP39 assembles in several closely-related helical geometries. This VP39 reconstruction reveals general principles for baculoviral nucleocapsid assembly.

The History of Baculovirology in Africa

Baculovirology has been studied on the African continent for the development of insect virus-based biopesticides and, to a much lesser extent, vaccine production and delivery, since the 1960s. In this review, we focus only on baculoviruses as biopesticides for agricultural pests in Africa. At least 11 species of baculovirus have been discovered or studied on the African continent, some with several distinct isolates, with the objective in most cases being the development of a biopesticide. These include the nucleopolyhedroviruses of Helicoverpa armigera, Cryptophlebia peltastica, Spodoptera exempta, Spodoptera frugiperda, Spodoptera littoralis, and Maruca vitrata, as well as the granuloviruses of Cydia pomonella, Plutella xylostella, Thaumatotibia (Cryptophlebia) leucotreta, Choristoneura occidentalis, and Phthorimaea operculella. Eleven different baculovirus-based biopesticides are recorded as being registered and commercially available on the African continent. Baculoviruses are recorded to have been isolated, researched, utilised in field trials, and/or commercially deployed as biopesticides in at least 13 different African countries. Baculovirus research is ongoing in Africa, and researchers are confident that further novel species and isolates will be discovered, to the benefit of environmentally responsible agricultural pest management, not only in Africa but also elsewhere.

Genome analysis of Psilogramma increta granulovirus and its intrapopulation diversity

The complete genome of Psilogramma increta granulovirus (PsinGV), isolated from P. increta (Lepidoptera: Sphingidae), was ultra-deep sequenced with a Novaseq PE150 platform and de novo assembled and annotated. The PsinGV genome is a circular double-stranded DNA, 103,721 bp in length, with a G+C content of 33.0%, the third lowest G+C content in present sequenced baculoviruses. It encodes 123 putative open reading frames, including 38 baculovirus core genes, 42 lepidopteran baculovirus conserved genes, 38 betabaculovirus conserved genes, and 5 genes unique to PsinGV. Meanwhile, 3 homologous repeated regions with the core sequence TTGCAA and 3 direct repeated sequences were identified within the PsinGV genome. Kimura two-parameters distance analysis confirmed that Psilogramma increta granulovirus is a representative of a prospective new species of the genus Betabaculovirus. Phylogenetic analysis based on the baculovirus core genes showed that PsinGV is closely related to Choristoneura fumiferana granulovirus, Clostera anastomosis granulovirus-B, and Erinnyis ello granulovirus. These four species therefore share a common ancestor residing in the Betabaculovirus genus. The genome of the PsinGV isolate contained two p10 copies: p10 and p10-2. Phylogenetic reconstruction of P10 suggests a transfer event of the p10-2 gene from an alphabaculovirus to the aforementioned common ancestor. Analysis of genomic diversity showed that 203 intrahost variants, including 182 single nucleotide variants and 21 short insertions/deletions, are present within the PsinGV isolate. Meanwhile, allele frequency indicated that the isolate contains three major genotypes, implying the archived isolate consists of several P. increta carcasses killed by PsinGV with different genotypes.

Genome sequence and organization of the Mythimna (formerly Pseudaletia) unipuncta granulovirus Hawaiian strain

Purified occlusion bodies (OBs) of Mythimna (formerly Pseudaletia) unipuncta (the true armyworm) granulovirus Hawaiian strain (MyunGV-A) were observed, showing typical GV morphological characteristics under scanning and transmission electron microscopy (EM). The genome of MyunGV-A was completely sequenced and analysed. The genome is 176,677 bp in size, with a G+C content of 39.79%. It contains 183 open reading frames (ORFs) encoding 50 or more amino acids with minimal overlap. Comparison of MyunGV-A with TnGV, XcGV, and HearGV genomes revealed extensive sequence similarity and collinearity, and the four genomes contain the same nine homologous regions (hrs) with conserved structures and locations. Three unique genes, 12 baculovirus repeated ORF (bro), 2 helicase, and 3 enhancin genes, were identified. In particular, two repeated genes (ORF39 and 49) are present in the genome, in reverse and complementarily orientations. Twenty-four OB proteins were identified from the putative protein database of MyunGV-A. In addition, MyunGV-A belongs to the Betabaculovirus group and is most closely related to TnGV (99% amino acid identity) according to a phylogenetic tree based on the combined amino acid sequences of 38 core gene contents.

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.

Bombyx mori nucleopolyhedrovirus Bm46 is essential for efficient production of infectious BV and nucleocapsid morphogenesis

Bombyx mori nucleopolyhedrovirus (BmNPV) orf46 (Bm46), the orthologues of Autographa californica multiple nucleopolyhedrovirus (AcMNPV) ac57, is a highly conserved gene in group Ⅰ and group Ⅱ nucleopolyhedroviruses (NPVs). However, its function in viral life cycle is unclear. Our results indicated that Bm46 transcript was detected from infected cells at 12 h post infection, while Bm46 protein was detectable from 24 to 72 h post infection. Upon the deletion of Bm46, fewer infectious BVs were produced by titer assays, but neither viral DNA synthesis nor occlusion bodies (OBs) production was affected. Electron microscopy revealed that Bm46 knockout interrupted nucleocapsid assembly and occlusion-derived virus (ODV) embedding, resulting in aberrant capsid-like tubular structures accumulated in the RZ (ring zone). Interestingly, this abnormally elongated capsid structures were consistent with the immunofluorescence microscopy results showing that VP39 assembled into long filaments and cables in the RZ. Moreover, DNA copies decreased by 30 % in occlusion bodies (OBs) produced by Bm46-knockout virus. qRT-PCR and Western blot analysis showed that the expression of VP39 was affected by Bm46 disruption. Taken together, our findings clearly pointed out that Bm46 played an important role in BV production and the proper formation of nucleocapsid morphogenesis.

Genome sequence analysis of a Helicoverpa armigera single nucleopolyhedrovirus (HearNPV-TR) isolated from Heliothis peltigera in Turkey

The entire genome of Helicoverpa armigera single nucleopolyhedrovirus (HearNPV-TR) was sequenced, and compared to genomes of other existing isolates. HearNPV-TR genome is 130.691 base pairs with a 38.9% G+C content and has 137 open reading frames (ORFs) of ≥ 150 nucleotides. Five homologous repeated sequences (hrs) and two baculovirus repeated ORFs (bro-a and bro-b) were identified. Phylogenetic analysis showed that HearNPV-TR is closer to HaSNPV-C1, HaSNPV-G4, HaSNPV-AU and HasNPV. However, there are significant differences in hr3, hr5 regions and in bro-a gene. Pairwise Kimura-2 parameter analysis of 38 core genes sequences of HearNPV-TR and other Helicoverpa NPVs showed that the genetic distances for these sequences were below 0.015 substitutions/site. Genomic differences as revealed by restriction profiles indicated that hr3, hr5 regions and bro-a gene may play a role in the virulence of HearNPV-TR.

The peritrophic membrane as a target of proteins that play important roles in plant defense and microbial attack

The peritrophic membrane (or peritrophic matrix: PM) is a thin membranous structure that lies along the midgut epithelium in the midgut lumen and consists of chitin and proteins. PM exists between ingested food material and midgut epithelium cells and it is on the frontline of insect-plant and insect-microbe interactions. Therefore, proteins that play major roles in plant defense against herbivorous insects and in microbial attack on insects should penetrate, destroy or modify the PM to accomplish their roles. Recently, it has become clear that some proteins crucial to plant defense or microbial attack have the PM as their primary target. In addition, several plant defense proteins have been reported to affect the PM, although it is still unclear whether the PM is their primary target. This review introduces several of these proteins: fusolin and enhancin, two proteins produced by insect viruses that greatly enhance infection of the viruses by disrupting the PM; the MLX56 family proteins found in mulberry latex as defense proteins against insect herbivores, which modify the PM to a thick structure that inhibits digestive processes; Mir1-CP, a defense cysteine protease from maize that inhibits the growth of insects at very low concentrations and degrades the PM structures; and chitinases and lectins. The importance, necessary characteristics, and modes of action of PM-targeting proteins are then discussed from a strategic point of view, by spotlighting the importance of selective permeability of the PM. Finally, the review discusses the possibility of applying PM-targeting proteins for the control of pest insects.

Genome sequence analysis and organization of the Hyphantria cunea granulovirus (HycuGV-Hc1) from Turkey

The fall webworm (Hyphantria cunea) impacts a wide variety of crops and cultivated broadleaf plant species. The pest is native to North America, was introduced to Europe and has since spread further as far as central Asia. Despite several attempts to control its distribution, the pest continues to spread causing damage all over the world. A naturally occurring baculovirus, Hyphantria cunea granulovirus (HycuGV-Hc1), isolated from the larvae of H. cunea in Turkey appears to have a potential as microbial control agent against this pest. In this report we describe the complete genome sequence and organization of the granulovirus isolate (HycuGV-Hc1) that infects the larval stages and compare it to other baculovirus genomes. The HycuGV-Hc1 genome is a circular double-stranded DNA of 114,825 bp in size with a nucleotide distribution of 39.3% G + C. Bioinformatics analysis predicted 132 putative open reading frames of (ORFs) ≥ 150 nucleotides. There are 24 ORFs with unknown function. Seven homologous repeated regions (hrs) and two bro genes (bro-1 and bro-2) were identified in the genome. Comparison to other baculovirus genomes, HycuGV-Hc1 revealed some differences in gene content and organization. Gene parity plots and phylogenetics confirmed that HycuGV-Hc1 is a Betabaculovirus and is closely related to Plutella xylostella granulovirus. This study expands our knowledge on the genetic variation of HycuGV isolates and provides further novel knowledge on the nature of granuloviruses.

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.

Genome Characteristics of the Cyclophragma Undans Nucleopolyhedrovirus: A Distinct Species in Group I of Alphabaculovirus

The Cyclophragma undans nucleopolyhedrovirus (CyunNPV), a potential pest control agent, was isolated from Cyclophragma undans (Lepidoptera: Lasiocampidae), an important forest pest. In the present study, we performed detailed genome analysis of CyunNPV and compared its genome to those of other Group I alphabaculoviruses. Sequencing of the CyunNPV genome using the Roche 454 sequencing system generated 142,900 bp with a G + C content of 45%. Genome analysis predicted a total of 147 hypothetical open reading frames comprising 38 baculoviral core genes, 24 lepidopteran baculovirus conserved genes, nine Group I Alphabaculovirus conserved genes, 71 common genes, and five genes that are unique to CyunNPV. In addition, the genome contains 13 homologous repeated sequences (hrs). Phylogenetic analysis groups CyunNPV under a distinct branch within clade "a" of Group I in the genus Alphabaculovirus. Unlike other members of Group I, CyunNPV harbors only nine of the 11 genes previously determined to be specific to Group I viruses. Furthermore, the CyunNPV lacks the tyrosine phosphatase gene and the ac30 gene. The CyunNPV F-like protein contains two insertions of continuous polar amino acids, one at the conventional fusion peptide and a second insertion at the pre-transmembrane domain. The insertions are likely to affect the fusion function and suggest an evolutionary process that led to inactivation of the F-like protein. The above findings imply that CyunNPV is a distinct species under Group I Alphabaculovirus.

Complete genome sequence analysis of the Malacosoma neustria nucleopolyhedrovirus from Turkey

The lackey moth, Malacosoma neustria (Linnaeus, 1758), a worldwide pest, causes extensive economic losses particularly on hazelnut, plum, oak, poplar, and willow trees. A baculovirus, Malacosoma neustria nucleopolyhedrovirus (ManeNPV-T2), has been isolated from the larvae collected in Turkey and appears to have a potential as a microbial control agent. In this study, we describe the complete genome sequence of ManeNPV-T2 and compare it to other sequenced baculovirus genomes. The ManeNPV-T2 genome is a circular double-stranded DNA molecule of 130,202 bp, has 38.2% G + C, and is predicted to contain 131 putative open reading frames (ORFs) each with a coding capacity of more then 50 amino acids. There are 27 ORFs with unknown function of which 6 are unique to ManeNPV-T2. Eleven homologous regions (hrs) and two bro genes (bro-a and bro-b) were identified in the genome. There are two homologues of chaB and nicotinamide riboside kinase-1 genes, separated from themselves with a few nucleotides. Additionally, ac145, thought to be per os infectivity factor (pif) gene, is also found as two homologues. All 38 core genes are found in the ManeNPV-T2 genome. The phylogenetic tree of ManeNPV-T2 in relation to 50 other baculoviruses whose genomes have been completely sequenced showed ManeNPV-T2 to be closely related to the group II NPVs. This study expands our knowledge on baculoviruses, describes the characterization ManeNPV, and ultimately contributes to the registration of this virus as a microbial pesticide.

The 38K-Mediated Specific Dephosphorylation of the Viral Core Protein P6.9 Plays an Important Role in the Nucleocapsid Assembly of Autographa californica Multiple Nucleopolyhedrovirus

Encapsidation of the viral genomes, leading to the assembly of the nucleocapsids to form infectious progeny virions, is a key step in many virus life cycles. Baculovirus nucleocapsid assembly is a complex process that involves many proteins. Our previous studies showed that the deletion of the core gene 38K (ac98) interrupted the nucleocapsid assembly by producing capsid sheaths devoid of viral genomes by an unknown mechanism. All homologs of 38K contain conserved motifs of the haloacid dehalogenase superfamily, which are involved in phosphoryl transfer. The requirements of these motifs for nucleocapsid assembly, confirmed in the present study, suggest that 38K may be a functioning haloacid dehalogenase. P6.9 is also encoded by a core gene (ac100) and is required for viral genome encapsidation. It has been reported that multiple phosphorylated species of P6.9 are present in virus-infected cells, while only an unphosphorylated species is detected in the budded virus. Therefore, whether 38K mediates the dephosphorylation of P6.9 was investigated. An additional phosphorylated species of P6.9 in 38K-deleted or -mutated virus-transfected cells was detected, and the dephosphorylated sites mediated by 38K were determined by mass spectrometry. To assess the effects of dephosphorylation of P6.9 mediated by 38K on virus replication, these sites were mutated to glutamic acids (phosphorylation-mimic mutant) or to alanines (phosphorylation-deficient mutant). Studies showed that the nucleocapsid assembly was interrupted in phosphorylation-mimic mutant virus-transfected cells. Taken together, our findings demonstrate that 38K mediates the dephosphorylation of specific sites at the C terminus of P6.9, which is essential for viral genome encapsidation.IMPORTANCE Genome packaging is a fundamental process in the virus life cycle, and viruses have different strategies to perform this step. For several double-stranded DNA (dsDNA) viruses, the procapsid is formed before genome encapsidation, which may require basic proteins that help to neutralize the nucleic acid charge repulsion to facilitate the compaction of the genome within the confined capsid space. Baculovirus encodes a small basic protein, P6.9, which is required for a variety of processes in the virus infection cycle. The phosphorylation of P6.9 is thought to result in nucleocapsid uncoating, while the dephosphorylation of P6.9 is involved in viral DNA encapsidation during nucleocapsid assembly. Here, we demonstrate that a haloacid dehalogenase homolog encoded by baculovirus core gene 38K is involved in nucleocapsid assembly by mediating the dephosphorylation of 5 specific sites at the C terminus of P6.9. This finding contributes to the understanding of the mechanisms of virus nucleocapsid assembly.

Genome analysis of a novel Group I alphabaculovirus obtained from Oxyplax ochracea

Oxyplax ochracea (Moore) is a pest that causes severe damage to a wide range of crops, forests and fruit trees. The complete genome sequence of Oxyplax ochracea nucleopolyhedrovirus (OxocNPV) was determined using a Roche 454 pyrosequencing system. OxocNPV has a double-stranded DNA (dsDNA) genome of 113,971 bp with a G+C content of 31.1%. One hundred and twenty-four putative open reading frames (ORFs) encoding proteins of >50 amino acids in length and with minimal overlapping were predicted, which covered 92% of the whole genome. Six baculoviral typical homologous regions (hrs) were identified. Phylogenetic analysis and gene parity plot analysis showed that OxocNPV belongs to clade “a” of Group I alphabaculoviruses, and it seems to be close to the most recent common ancestor of Group I alphabaculoviruses. Three unique ORFs (with no homologs in the National Center for Biotechnology Information database) were identified. Interestingly, OxocNPV lacks three auxiliary genes (lef7, ie-2 and pcna) related to viral DNA replication and RNA transcription. In addition, OxocNPV has significantly different sequences for several genes (including ie1 and odv-e66) in comparison with those of other baculoviruses. However, three dimensional structure prediction showed that OxocNPV ODV-E66 contain the conserved catalytic residues, implying that it might possess polysaccharide lyase activity as AcMNPV ODV-E66. All these unique features suggest that OxocNPV represents a novel species of the Group I alphabaculovirus lineage.

Genomic analysis of a Trichoplusia ni Betabaculovirus (TnGV) with three different viral enhancing factors and two unique genes

The complete genome of a Trichoplusia ni granulovirus (TnGV) is described and analyzed. The genome contains 175,360 bp (KU752557), becoming the third largest genome within the genus Betabaculovirus, smaller only than the Xestia c-nigrum GV (XecnGV) (178,733 pb) and the Pseudaletia unipuncta GV (PsunGV) (176,677 pb) genomes. The TnGV genome has a 39.81% C+G content and a total of 180 ORFs were identified, 96 of them in the granulin gene direction and 84 in the opposite direction. A total of 94.38% of the ORFs showed high identity with those of ClanGV, HaGV, and SlGV. Eight homologous regions (hrs) were identified as well as one apoptosis inhibitor (IAP-3). Interestingly, three viral enhancing factors (VEFs) were located in TnGV genome: VEF-1 (orf153), VEF-3 (orf155), and VEF-4 (orf164), additional to another metalloprotease (orf37). Two ORFs were unique to TnGV (orf100 and orf101) and another one was shared by only TnGV and AgseGV (orf2). Eleven of the deduced proteins showed high identity with proteins from nucleopolyhedroviruses, three with proteins from ascoviruses, and one with an entomopoxvirus protein. The largest deduced protein contains 1,213 amino acids (orf43) and the smallest deduced protein contains only 50 amino acids (orf143). Sequence identity and phylogenetic analyses showed that the closest related genomes to TnGV are, to date, those of PsunGV and XecnGV. This genome analysis may contribute to functional research on TnGV, and may form the bases for the utilization of this betabaculovirus as a pest control agent.

Whole genome sequence and comparative genomic sequence analysis of Helicoverpa armigera nucleopolyhedrovirus (HearNPV-L1) isolated from India

The whole genome of Helicoverpa armigera nucleopolyhedrovirus (HearNPV) from India, HearNPV-L1, was sequenced and analyzed, with a view to look for genes and/or nucleotide sequences that might be involved in the differences and virulence among other HearNPVs sequenced from other countries like SP1A (Spain), NNg1 (Kenya) and G4 (China). The entire nucleotide sequence of the HearNPV-L1 genome was 136,740 bp in length having GC content of 39.19% and contained 113 ORFs that could encode polypeptides with more than 50 amino acids (GenBank accession number KT013224). Two ORFs, viz., ORF 18 (300 bp) and ORF 19 (401 bp) identified were unique in HearNPV-L1 genome. Most of the HearNPV-L1 ORFs showed high similarity to NNg1, SP1A and G4 genomes. HearNPV-L1 genome contains 5 h (hr1-hr5), these regions were found 84-100% similar to hr region of NNg1, SP1A and G4 genomes. A total of four bro genes were observed in HearNPV-L1 genome, of which bro-a gene was 12 and 351 bp bigger than SP1A and G4 bro-a, respectively, while bro-b was 15 bp bigger SP1A and NNg1 bro-b, whereas 593 bp shorter than G4 bro-b, while bro-c was 12 bp shorter than NNg1, however bro-c was absent in G4 genome. HearNPV-L1 bro-d was 100% homologous to bro-d of SP1A, NNg1 and G4 genomes, respectively. The comparative analysis of HearNPV-L1 genome indicated that there are several other putative genes and nucleotide sequences that may be responsible for insecticidal activity in HearNPV-L1 isolate, however, further functional analysis of the hypothetical (putative) genes may help identifying the genes that are crucial for the virulence and insecticidal activity.

The Complete Genome Sequence of a Second Distinct Betabaculovirus from the True Armyworm, Mythimna unipuncta

The betabaculovirus originally called Pseudaletia (Mythimna) sp. granulovirus #8 (MyspGV#8) was examined by electron microscopy, host barcoding PCR, and determination of the nucleotide sequence of its genome. Scanning and transmission electron microscopy revealed that the occlusion bodies of MyspGV#8 possessed the characteristic size range and morphology of betabaculovirus granules. Barcoding PCR using cytochrome oxidase I primers with DNA from the MyspGV#8 collection sample confirmed that it had been isolated from the true armyworm, Mythimna unipuncta (Lepidoptera: Noctuidae) and therefore was renamed MyunGV#8. The MyunGV#8 genome was found to be 144,673 bp in size with a nucleotide distribution of 49.9% G+C, which was significantly smaller and more GC-rich than the genome of Pseudaletia unipuncta granulovirus H (PsunGV-H), another M. unipuncta betabaculovirus. A phylogeny based on concatenated baculovirus core gene amino acid sequence alignments placed MyunGV#8 in clade a of genus Betabaculovirus. Kimura-2-parameter nucleotide distances suggested that MyunGV#8 represents a virus species different and distinct from other species of Betabaculovirus. Among the 153 ORFs annotated in the MyunGV#8 genome, four ORFs appeared to have been obtained from or donated to the alphabaculovirus lineage represented by Leucania separata nucleopolyhedrovirus AH1 (LeseNPV-AH1) during co-infection of Mythimna sp. larvae. A set of 33 ORFs was identified that appears only in other clade a betabaculovirus isolates. This clade a-specific set includes an ORF that encodes a polypeptide sequence containing a CIDE_N domain, which is found in caspase-activated DNAse/DNA fragmentation factor (CAD/DFF) proteins. CAD/DFF proteins are involved in digesting DNA during apoptosis.

The Complete Genome Sequence of Plodia Interpunctella Granulovirus: Evidence for Horizontal Gene Transfer and Discovery of an Unusual Inhibitor-of-Apoptosis Gene

The Indianmeal moth, Plodia interpunctella (Lepidoptera: Pyralidae), is a common pest of stored goods with a worldwide distribution. The complete genome sequence for a larval pathogen of this moth, the baculovirus Plodia interpunctella granulovirus (PiGV), was determined by next-generation sequencing. The PiGV genome was found to be 112, 536 bp in length with a 44.2% G+C nucleotide distribution. A total of 123 open reading frames (ORFs) and seven homologous regions (hrs) were identified and annotated. Phylogenetic inference using concatenated alignments of 36 baculovirus core genes placed PiGV in the “b” clade of viruses from genus Betabaculovirus with a branch length suggesting that PiGV represents a distinct betabaculovirus species. In addition to the baculovirus core genes and orthologues of other genes found in other betabaculovirus genomes, the PiGV genome sequence contained orthologues of the bidensovirus NS3 gene, as well as ORFs that occur in alphabaculoviruses but not betabaculoviruses. While PiGV contained an orthologue of inhibitor of apoptosis-5 (iap-5), an orthologue of inhibitor of apoptosis-3 (iap-3) was not present. Instead, the PiGV sequence contained an ORF (PiGV ORF81) encoding an IAP homologue with sequence similarity to insect cellular IAPs, but not to viral IAPs. Phylogenetic analysis of baculovirus and insect IAP amino acid sequences suggested that the baculovirus IAP-3 genes and the PiGV ORF81 IAP homologue represent different lineages arising from more than one acquisition event. The presence of genes from other sources in the PiGV genome highlights the extent to which baculovirus gene content is shaped by horizontal gene transfer.

An alphabaculovirus isolated from dead Lymantria dispar larvae shows high genetic similarity to baculovirus previously isolated from Lymantria monacha - An example of adaptation to a new host

A new isolate of baculovirus, Lymantria dispar multiple nucleopolyhedrovirus-BNP (LdMNPV-BNP), was found in dead gypsy moth (L. dispar) caterpillars collected in the Biebrzanski National Park in Poland. Here, we examined its biological activity, structure, genetic content and phylogeny. Multiple nucleocapsids of LdMNPV-BNP are enveloped together in 2-26 virions embedded in occluded bodies (OBs) very similar to the OBs previously described in viruses infecting Lymantriinae. This isolate kills pest larvae in a relatively short time (LT50 of approximately 9days for a dose of 2×10(7)OBs/ml), highlighting the possibility for its use as a biopesticide. Next-generation sequencing of LdMNPV-BNP revealed gene content (e.g. DNA photolyase) that is not present in any LdMNPV isolate sequenced to date. The genome is 157,270 base pairs long and has a notably lower G+C content in comparison to other LdMNPVs (50.3% G+C content compared to an average of 57.4% among other LdMNPVs). According to our phylogenetic analysis based on 37 core genes, LdMNPV-BNP is a member of group II alphabaculoviruses, which are closely related to LdMNPV and LyxyMNPV (Lymantria xylina multiple nucleopolyhedrovirus). Molecular evolution inference based on the partial sequence of lef-8, lef-9 and polh genes shows that LdMNPV-BNP and isolates of Lymantria monacha nucleopolyhedrovirus (LymoNPV) may share a very recent common ancestor or be isolates of the same virus species. LdMNPV-BNP, like other baculoviruses, could be beneficial as an active component of biopesticides that can be used during forest integrated pest management.

The complete genome of a baculovirus isolated from an insect of medical interest: Lonomia obliqua (Lepidoptera: Saturniidae)

Lonomia obliqua (Lepidoptera: Saturniidae) is a species of medical importance due to the severity of reactions caused by accidental contact with the caterpillar bristles. Several natural pathogens have been identified in L. obliqua, and among them the baculovirus Lonomia obliqua multiple nucleopolyhedrovirus (LoobMNPV). The complete genome of LoobMNPV was sequenced and shown to have 120,022 bp long with 134 putative open reading frames (ORFs). Phylogenetic analysis of the LoobMNPV genome showed that it belongs to Alphabaculovirus group I (lepidopteran-infective NPV). A total of 12 unique ORFs were identified with no homologs in other sequenced baculovirus genomes. One of these, the predicted protein encoded by loob035, showed significant identity to an eukaryotic transcription terminator factor (TTF2) from the Lepidoptera Danaus plexippus, suggesting an independent acquisition through horizontal gene transfer. Homologs of cathepsin and chitinase genes, which are involved in host integument liquefaction and viral spread, were not found in this genome. As L. obliqua presents a gregarious behavior during the larvae stage the impact of this deletion might be neglectable.

Genome Sequencing and Analysis of Catopsilia pomona nucleopolyhedrovirus: A Distinct Species in Group I Alphabaculovirus

The genome sequence of Catopsilia pomona nucleopolyhedrovirus (CapoNPV) was determined by the Roche 454 sequencing system. The genome consisted of 128,058 bp and had an overall G+C content of 40%. There were 130 hypothetical open reading frames (ORFs) potentially encoding proteins of more than 50 amino acids and covering 92% of the genome. Among all the hypothetical ORFs, 37 baculovirus core genes, 23 lepidopteran baculovirus conserved genes and 10 genes conserved in Group I alphabaculoviruses were identified. In addition, the genome included regions of 8 typical baculoviral homologous repeat sequences (hrs). Phylogenic analysis showed that CapoNPV was in a distinct branch of clade “a” in Group I alphabaculoviruses. Gene parity plot analysis and overall similarity of ORFs indicated that CapoNPV is more closely related to the Group I alphabaculoviruses than to other baculoviruses. Interesting, CapoNPV lacks the genes encoding the fibroblast growth factor (fgf) and ac30, which are conserved in most lepidopteran and Group I baculoviruses, respectively. Sequence analysis of the F-like protein of CapoNPV showed that some amino acids were inserted into the fusion peptide region and the pre-transmembrane region of the protein. All these unique features imply that CapoNPV represents a member of a new baculovirus species.

The Autographa californica Multiple Nucleopolyhedrovirus ac54 Gene Is Crucial for Localization of the Major Capsid Protein VP39 at the Site of Nucleocapsid Assembly

Baculovirus DNAs are synthesized and inserted into preformed capsids to form nucleocapsids at a site in the infected cell nucleus, termed the virogenic stroma. Nucleocapsid assembly ofAutographa californicamultiple nucleopolyhedrovirus (AcMNPV) requires the major capsid protein VP39 and nine minor capsid proteins, including VP1054. However, how VP1054 participates in nucleocapsid assembly remains elusive. In this study, the VP1054-encoding gene (ac54) was deleted to generate theac54-knockout AcMNPV (vAc54KO). In vAc54KO-transfected cells, nucleocapsid assembly was disrupted, leading to the formation of abnormally elongated capsid structures. Interestingly, unlike cells transfected with AcMNPV mutants lacking other minor capsid proteins, in which capsid structures were distributed within the virogenic stroma,ac54ablation resulted in a distinctive location of capsid structures and VP39 at the periphery of the nucleus. The altered distribution pattern of capsid structures was also observed in cells transfected with AcMNPV lacking BV/ODV-C42 or in cytochalasind-treated AcMNPV-infected cells. BV/ODV-C42, along with PP78/83, has been shown to promote nuclear filamentous actin (F-actin) formation, which is another requisite for nucleocapsid assembly. Immunofluorescence using phalloidin indicated that the formation and distribution of nuclear F-actin were not affected byac54deletion. However, immunoelectron microscopy revealed that BV/ODV-C42, PP78/83, and 38K failed to integrate into capsid structures in the absence of VP1054, and immunoprecipitation further demonstrated that in transient expression assays, VP1054 interacted with BV/ODV-C42 and VP80 but not VP39. Our findings suggest that VP1054 plays an important role in the transport of capsid proteins to the nucleocapsid assembly site prior to the process of nucleocapsid assembly.

IMPORTANCE

Baculoviruses are large DNA viruses whose replication occurs within the host nucleus. The localization of capsids into the capsid assembly site requires virus-induced nuclear F-actin; the inhibition of nuclear F-actin formation results in the retention of capsid structures at the periphery of the nucleus. In this paper, we note that the minor capsid protein VP1054 is essential for the localization of capsid structures, the major capsid protein VP39, and the minor capsid protein 38K into the capsid assembly site. Moreover, VP1054 is crucial for correct targeting of the nuclear F-actin factors BV/ODV-C42 and PP78/83 for capsid maturation. However, the formation and distribution of nuclear F-actin are not affected by the lack of VP1054. We further reveal that VP1054 interacts with BV/ODV-C42 and a capsid transport-related protein, VP80. Taken together, our findings suggest that VP1054 plays a unique role in the pathway(s) for transport of capsid proteins.

Molecular and in Silico Characterization of Achaea janata Granulovirus Granulin Gene

Achaea janata granulovirus (AcjaGV), an insect virus belonging to Baculoviridae, infects semilooper, a widely distributed defoliating pest on castor beans (Ricinus communis L.) and several other plant hosts in India. The propagation and purification of the Hyderabad isolate AcjaGV were performed, granulin gene from this isolate was amplified, cloned and sequenced, and its homology with other known granulin genes was assessed. The 753-bp granulin ORF of AcjaGV encoded for a granulin protein of 250 amino acids with a molecular mass of 29.5 ± 0.7 kDa. This amino acid sequence exhibited significant homology with Spodoptera litura granulovirus (SpliGV) and other GVs infecting insects in the same Noctuidae family of Lepidoptera. Peptide analysis of granulin protein indicated close homology with that of SpliGV. Virtual RFLP patterns from in silico digestions of granulin gene of 18 granuloviruses mapped by 12 restriction enzymes were used for simulated digestions. Implications of the phylogenetic relationships of granulin nucleotide and deduced amino acid sequence are discussed. We have established the sequence identity of granulin gene of AcjaGV and characterized its protein product and the phylogenetic relationship with other known GVs. Our results indicate the presence of unique restriction sites for three restriction enzymes, and this can be used as a tool for identification of AcjaGV from various sources. This is the first report from the Indian subcontinent to describe the complete granulin gene of a GV isolated from A. janata.

Genome of Cnaphalocrocis medinalis Granulovirus, the First Crambidae-Infecting Betabaculovirus Isolated from Rice Leaffolder to Sequenced

Cnaphalocrocis medinalis is a major pest of rice in South and South-East Asia. Insecticides are the major means farmers use for management. A naturally occurring baculovirus, C. medinalis granulovirus (CnmeGV), has been isolated from the larvae and this has the potential for use as microbial agent. Here, we described the complete genome sequence of CnmeGV and compared it to other baculovirus genomes. The genome of CnmeGV is 112,060 base pairs in length, has a G+C content of 35.2%. It contains 133 putative open reading frames (ORFs) of at least 150 nucleotides. A hundred and one (101) of these ORFs are homologous to other baculovirus genes including 37 baculovirus core genes. Thirty-two (32) ORFs are unique to CnmeGV with no homologues detected in the GeneBank and 53 tandem repeats (TRs) with sequence length from 25 to 551 nt intersperse throughout the genome of CnmeGV. Six (6) homologous regions (hrs) were identified interspersed throughout the genome. Hr2 contains 11 imperfect palindromes and a high content of AT sequence (about 73%). The unique ORF28 contains a coiled-coil region and a zinc finger-like domain of 4–50 residues specialized by two C2C2 zinc finger motifs that putatively bound two atoms of zinc. ORF21 encoding a chit-1 protein suggesting a horizontal gene transfer from alphabaculovirus. The putative protein presents two carbohydrate-binding module family 14 (CBM_14) domains rather than other homologues detected from betabaculovirus that only contains one chit-binding region. Gene synteny maps showed the colinearity of sequenced betabaculovirus. Phylogenetic analysis indicated that CnmeGV grouped in the betabaculovirus, with a close relation to AdorGV. The cladogram obtained in this work grouped the 17 complete GV genomes in one monophyletic clade. CnmeGV represents a new crambidae host-isolated virus species from the genus Betabaculovirus and is most closely relative of AdorGV. The analyses and information derived from this study will provide a better understanding of the pathological symptoms caused by this virus and its potential use as a microbial pesticide.

Entomopoxvirus infection induces changes in both juvenile hormone and ecdysteroid levels in larval Mythimna separata

Insect viruses are among the most important pathogens of lepidopteran insects. Virus-infected larvae often show developmental defects including a prolonged larval period and a failure to pupate, but the mechanisms by which insect viruses regulate host development need further investigation. In this study, the regulation of host endocrinology by a lepidopteran entomopoxvirus (EPV), Mythimna separata EPV (MySEV), was examined. When fourth instar M. separata were inoculated with MySEV occlusion bodies, pupation was prevented and the insects died during the final (sixth) larval instar. Liquid chromatography-MS analysis revealed that juvenile hormone (JH) titres in the haemolymph of MySEV-infected sixth instars were higher than those in mock-infected larvae. JH esterase (JHE) activity was also examined by kinetic assay using a colorimetric substrate. The level of JHE activity in the haemolymph of MySEV-infected larvae was generally lower than that found in mock-infected larvae. In contrast, ecdysteroid titre in the haemolymph of final-instar MySEV-infected larvae was lower than that found in mock-infected larvae when measured by radioimmunoassay. A statistically significant difference in the release of ecdysteroids from prothoracic glands (PGs) that were dissected from MySEV- or mock-infected sixth instar Day 3 larvae was not found following prothoracicotropic hormone (PTTH) exposure. Our results indicate that the release of ecdysteroids was reduced not by infection of the PGs by MySEV, but by reduced PTTH production from the brain. Taken together our study suggests that EPVs retard host development by both reducing ecdysone titre and maintaining status quo levels of JH by preventing its metabolism.

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.

The Complete Genome of a New Betabaculovirus from Clostera anastomosis

Clostera anastomosis (Lepidoptera: Notodontidae) is a defoliating forest insect pest. Clostera anastomosis granulovirus-B (ClasGV-B) belonging to the genus Betabaculovirus of family Baculoviridae has been used for biological control of the pest. Here we reported the full genome sequence of ClasGV-B and compared it to other previously sequenced baculoviruses. The circular double-stranded DNA genome is 107,439 bp in length, with a G+C content of 37.8% and contains 123 open reading frames (ORFs) representing 93% of the genome. ClasGV-B contains 37 baculovirus core genes, 25 lepidopteran baculovirus specific genes, 19 betabaculovirus specific genes, 39 other genes with homologues to baculoviruses and 3 ORFs unique to ClasGV-B. Hrs appear to be absent from the ClasGV-B genome, however, two non-hr repeats were found. Phylogenetic tree based on 37 core genes from 73 baculovirus genomes placed ClasGV-B in the clade b of betabaculoviruses and was most closely related to Erinnyis ello GV (ErelGV). The gene arrangement of ClasGV-B also shared the strongest collinearity with ErelGV but differed from Clostera anachoreta GV (ClanGV), Clostera anastomosis GV-A (ClasGV-A, previously also called CaLGV) and Epinotia aporema GV (EpapGV) with a 20 kb inversion. ClasGV-B genome contains three copies of polyhedron envelope protein gene (pep) and phylogenetic tree divides the PEPs of betabaculoviruses into three major clades: PEP-1, PEP-2 and PEP/P10. ClasGV-B also contains three homologues of P10 which all harbor an N-terminal coiled-coil domain and a C-terminal basic sequence. ClasGV-B encodes three fibroblast growth factor (FGF) homologues which are conserved in all sequenced betabaculoviruses. Phylogenetic analysis placed these three FGFs into different groups and suggested that the FGFs were evolved at the early stage of the betabaculovirus expansion. ClasGV-B is different from previously reported ClasGV-A and ClanGV isolated from Notodontidae in sequence and gene arrangement, indicating the virus is a new notodontid betabaculovirus.

Two year field study to evaluate the efficacy of Mamestra brassicae nucleopolyhedrovirus combined with proteins derived from Xestia c-nigrum granulovirus

Japan has only three registered baculovirus biopesticides despite its long history of studies on insect viruses. High production cost is one of the main hindrances for practical use of baculoviruses. Enhancement of insecticidal effect is one possible way to overcome this problem, so there have been many attempts to develop additives for baculoviruses. We found that alkaline soluble proteins of capsules (GVPs) of Xestia c-nigrum granulovirus can increase infectivity of some viruses including Mamestra brassicae nucleopolyhedrovirus (MabrNPV), and previously reported that MabrNPV mixed with GVPs was highly infectious to three important noctuid pests of vegetables in the following order, Helicoverpa armigera, M. brassicae, and Autographa nigrisigna. In this study, small-plot experiments were performed to assess concentrations of MabrNPV and GVPs at three cabbage fields and a broccoli field for the control of M. brassicae. In the first experiment, addition of GVPs (10 µg/mL) to MabrNPV at 10⁶ OBs/mL resulted in a significant increase in NPV infection (from 53% to 66%). In the second experiment, the enhancing effect of GVP on NPV infection was confirmed at 10-times lower concentrations of MabrNPV. In the third and fourth experiments, a 50% reduction in GVPs (from 10 µg/mL to 5 µg/mL) did not result in a lowering of infectivity of the formulations containing MabrNPV at 10⁵ OBs/mL. These results indicate that GVPs are promising additives for virus insecticides.

Reaching the melting point: Degradative enzymes and protease inhibitors involved in baculovirus infection and dissemination

Baculovirus infection of a host insect involves several steps, beginning with initiation of virus infection in the midgut, followed by dissemination of infection from the midgut to other tissues in the insect, and finally culminating in "melting" or liquefaction of the host, which allows for horizontal spread of infection to other insects. While all of the viral gene products are involved in ultimately reaching this dramatic infection endpoint, this review focuses on two particular types of baculovirus-encoded proteins: degradative enzymes and protease inhibitors. Neither of these types of proteins is commonly found in other virus families, but they both play important roles in baculovirus infection. The types of degradative enzymes and protease inhibitors encoded by baculoviruses are discussed, as are the roles of these proteins in the infection process.

The complete sequence of the first Spodoptera frugiperda Betabaculovirus genome: a natural multiple recombinant virus

Spodoptera frugiperda (Lepidoptera: Noctuidae) is a major pest in maize crops in Colombia, and affects several regions in America. A granulovirus isolated from S. frugiperda (SfGV VG008) has potential as an enhancer of insecticidal activity of previously described nucleopolyhedrovirus from the same insect species (SfMNPV). The SfGV VG008 genome was sequenced and analyzed showing circular double stranded DNA of 140,913 bp encoding 146 putative ORFs that include 37 Baculoviridae core genes, 88 shared with betabaculoviruses, two shared only with betabaculoviruses from Noctuide insects, two shared with alphabaculoviruses, three copies of own genes (paralogs) and the other 14 corresponding to unique genes without representation in the other baculovirus species. Particularly, the genome encodes for important virulence factors such as 4 chitinases and 2 enhancins. The sequence analysis revealed the existence of eight homologous regions (hrs) and also suggests processes of gene acquisition by horizontal transfer including the SfGV VG008 ORFs 046/047 (paralogs), 059, 089 and 099. The bioinformatics evidence indicates that the genome donors of mentioned genes could be alpha- and/or betabaculovirus species. The previous reported ability of SfGV VG008 to naturally co-infect the same host with other virus show a possible mechanism to capture genes and thus improve its fitness.

The genome sequence of Agrotis segetum nucleopolyhedrovirus B (AgseNPV-B) reveals a new baculovirus species within the Agrotis baculovirus complex

The genome of Agrotis segetum nucleopolyhedrovirus B (AgseNPV-B) was completely sequenced and compared with whole genome sequences of the Agrotis segetum nucleopolyhedrovirus A (AgseNPV-A) and Agrotis ipsilon nucleopolyhedrovirus (AgipNPV). The AgseNPV-B genome is 148,981 bp in length and encodes 150 putative open reading frames. AgseNPV-B contains two copies of the gene viral enhancing factor (vef), making the Agrotis nucleopolyhedroviruses and A. segetum granulovirus (AgseGV) very rich in vef in comparison to other baculoviruses. Genome alignments of AgseNPV-B, AgseNPV-A and AgipNPV showed a very high genome co-linearity interspersed with variable regions, which are considered as putative sites of genomic recombination. Phylogenetic analyses revealed that all three viruses are distinct. However, AgseNPV-B is more closely related to AgipNPV suggesting that both viruses are at an early stage of phylogenetic divergence. It is proposed that AgseNPV-B belongs to a third Alphabaculovirus species of the Agrotis baculovirus complex. The Agrotis exclamationis nucleopolyhedrovirus (AgexNPV) shared high nucleotide sequence identities with AgseNPV-B, suggesting it is actually an AgseNPV-B isolate.

Genome sequence of Erinnyis ello granulovirus (ErelGV), a natural cassava hornworm pesticide and the first sequenced sphingid-infecting betabaculovirus

Background

Cassava (Manihot esculenta) is the basic source for dietary energy of 500 million people in the world. In Brazil, Erinnyis ello ello (Lepidoptera: Sphingidae) is a major pest of cassava crops and a bottleneck for its production. In the 1980s, a naturally occurring baculovirus was isolated from E. ello larva and successfully applied as a bio-pesticide in the field. Here, we described the structure, the complete genome sequence, and the phylogenetic relationships of the first sphingid-infecting betabaculovirus.

Results

The baculovirus isolated from the cassava hornworm cadavers is a betabaculovirus designated Erinnyis ello granulovirus (ErelGV). The 102,759 bp long genome has a G + C content of 38.7%. We found 130 putative ORFs coding for polypeptides of at least 50 amino acid residues. Only eight genes were found to be unique. ErelGV is closely related to ChocGV and PiraGV isolates. We did not find typical homologous regions and cathepsin and chitinase homologous genes are lacked. The presence of he65 and p43 homologous genes suggests horizontal gene transfer from Alphabaculovirus. Moreover, we found a nucleotide metabolism-related gene and two genes that could be acquired probably from Densovirus.

Conclusions

The ErelGV represents a new virus species from the genus Betabaculovirus and is the closest relative of ChocGV. It contains a dUTPase-like, a he65-like, p43-like genes, which are also found in several other alpha- and betabaculovirus genomes, and two Densovirus-related genes. Importantly, recombination events between insect viruses from unrelated families and genera might drive baculovirus genomic evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-856) contains supplementary material, which is available to authorized users.

LITAF mutations associated with Charcot-Marie-Tooth disease 1C show mislocalization from the late endosome/lysosome to the mitochondria

Charcot-Marie-Tooth (CMT) disease is one of the most common heritable neuromuscular disorders, affecting 1 in every 2500 people. Mutations in LITAF have been shown to be causative for CMT type 1C disease. In this paper we explore the subcellular localization of wild type LITAF and mutant forms of LITAF known to cause CMT1C (T49M, A111G, G112S, T115N, W116G, L122V and P135T). The results show that LITAF mutants A111G, G112S, W116G, and T115N mislocalize from the late endosome/lysosome to the mitochondria while the mutants T49M, L122V, and P135T show partial mislocalization with a portion of the total protein present in the late endosome/lysosome and the remainder of the protein localized to the mitochondria. This suggests that different mutants of LITAF will produce differing severity of disease. We also explored the effect of the presence of mutant LITAF on wild-type LITAF localization. We showed that in cells heterozygous for LITAF, CMT1C mutants T49M and G112S are dominant since wild-type LITAF localized to the mitochondria when co-transfected with a LITAF mutant. Finally, we demonstrated how LITAF transits to the endosome and mitochondria compartments of the cell. Using Brefeldin A to block ER to Golgi transport we demonstrated that wild type LITAF traffics through the secretory pathway to the late endosome/lysosome while the LITAF mutants transit to the mitochondria independent of the secretory pathway. In addition, we demonstrated that the C-terminus of LITAF is necessary and sufficient for targeting of wild-type LITAF to the late endosome/lysosome and the mutants to the mitochondria. Together these data provide insight into how mutations in LITAF cause CMT1C disease.

Genome sequence and analysis of Buzura suppressaria nucleopolyhedrovirus: a group II Alphabaculovirus

The genome of Buzura suppressaria nucleopolyhedrovirus (BusuNPV) was sequenced by 454 pyrosequencing technology. The size of the genome is 120,420 bp with 36.8% G+C content. It contains 127 hypothetical open reading frames (ORFs) covering 90.7% of the genome and includes the 37 conserved baculovirus core genes, 84 genes found in other baculoviruses, and 6 unique ORFs. No typical baculoviral homologous repeats (hrs) were present but the genome contained a region of repeated sequences. Gene Parity Plots revealed a 28.8 kb region conserved among the alpha- and beta-baculoviruses. Overall comparisons of BusuNPV to other baculoviruses point to a distinct species in group II Alphabaculovirus.

Analysis of the genome of a Korean isolate of the Pieris rapae granulovirus enabled by its separation from total host genomic DNA by pulse-field electrophoresis

BACKGROUND

Most traditional genome sequencing projects involving viruses include the culture and purification of the virus particles. However, purification of virions may yield insufficient material for traditional sequencing. The electrophoretic method described here provides a strategy whereby the genomic DNA of the Korean isolate of Pieris rapae granulovirus (PiraGV-K) could be recovered in sufficient amounts for sequencing by purifying it directly from total host DNA by pulse-field gel electrophoresis (PFGE).

METHODOLOGY/PRINCIPAL FINDINGS

The total genomic DNA of infected P. rapae was embedded in agarose plugs, treated with restriction nuclease and methylase, and then PFGE was used to separate PiraGV-K DNA from the DNA of P. rapae, followed by mapping of fosmid clones of the purified viral DNA. The double-stranded circular genome of PiraGV-K was found to encode 120 open reading frames (ORFs), which covered 92% of the sequence. BLAST and ORF arrangement showed the presence of 78 homologs to other genes in the database. The mean overall amino acid identity of PiraGV-K ORFs was highest with the Chinese isolate of PiraGV (~99%), followed up with Choristoneura occidentalis ORFs at 58%. PiraGV-K ORFs were grouped, according to function, into 10 genes involved in transcription, 11 involved in replication, 25 structural protein genes, and 15 auxiliary genes. Genes for Chitinase (ORF 10) and cathepsin (ORF 11), involved in the liquefaction of the host, were found in the genome.

CONCLUSIONS/SIGNIFICANCE

The recovery of PiraGV-K DNA genome by pulse-field electrophoretic separation from host genomic DNA had several advantages, compared with its isolation from particles harvested as virions or inclusions from the P. rapae host. We have sequenced and analyzed the 108,658 bp PiraGV-K genome purified by the electrophoretic method. The method appears to be generally applicable to the analysis of genomes of large viruses.

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.

Molecular and immunohistochemical characterization of the chitinase gene from Pieris rapae granulovirus

The chitinase gene of baculoviruses is expressed in the late phase of virus replication in insects and possesses high exo- and endochitinase activity, which can hydrolyze chitin in the body of the insect, thus promoting terminal host liquefaction. Alphabaculovirus viral chitinases (vChitA) have been well analyzed, but information regarding viral chitinases from betabaculoviruses is limited. Whole-genome sequencing of a Korean isolate of Pieris rapae GV (PiraGV-K) predicted a putative chitinase gene corresponding to ORF10. The PiraGV-K chitinase gene had a coding sequence of 1,761 bp, encoding a protein of 586 amino acid (aa) residues, including an 18-aa putative signal peptide. Time course induction pattern observed by SDS-PAGE and subsequent Western blot with anti-PiraGV-K chitinase antibody revealed the cleavage of the signal peptide from the intact chitinase. Edman sequencing analysis was further conducted to confirm the exact nature of the mature chitinase, and the N-terminal amino acid sequence (KPGAP) exactly matched the sequence following the signal peptide sequence. The transcriptomics of PiraGV-K chitinase in infected P. rapae larvae, examined by real-time PCR, revealed a significant 75-fold increase after four days of feeding with PiraGV-K-treated leaves, with a subsequent decline at the later stages of infection. Confocal microscopic analysis showed that PiraGV-K chitinase possibly exists as a secreted protein, with strong chitinase-specific signals in fat body cells and integument at four days postinfection. Furthermore, immunogold labeling and electron microscopy studies localized the PiraGV-K chitinase in the cytoplasm and sparsely within vacuolar structures in the fat body apart from the extensive aggregation in the cuticular lining of the integument.

Origin of Ecdysosteroid UDP-glycosyltransferases of Baculoviruses through Horizontal Gene Transfer from Lepidoptera

Baculoviruses infecting Lepidoptera (butterflies and moths) encodes an enzyme known as ecdysosteroid UDP-glycosyltransferase (EGT), which inactivates insect host ecdysosteroid hormones, thereby preventing molt and pupation and permitting a build-up of the viral population within the host. Baculovirus EGT shows evidence of homology to insect UDP-glycosyltransferases, and a phylogenetic analysis supported the closest relative of baculovirus EGT are the UGT33 and UGT34 families of lepidopteran UDP-glycosyltransferases. The phylogenetic analysis thus supported that baculovirus EGT arose by horizontal gene transfer of a UDP-glycosyltransferase from a lepidopteran host, an event that occurred 70 million years ago at the earliest but possibly much more recently. Three amino acid replacements unique to baculovirus EGTs and conserved in all available baculovirus sequences were identified in the N-terminal region of the molecule. Because of their conservation, these amino acids are candidates for playing an important functional role in baculovirus EGT function.

Genome of Epinotia aporema granulovirus (EpapGV), a polyorganotropic fast killing betabaculovirus with a novel thymidylate kinase gene

BACKGROUND

Epinotia aporema (Lepidoptera: Tortricidae) is an important pest of legume crops in South America. Epinotia aporema granulovirus (EpapGV) is a baculovirus that causes a polyorganotropic infection in the host larva. Its high pathogenicity and host specificity make EpapGV an excellent candidate to be used as a biological control agent.

RESULTS

The genome of Epinotia aporema granulovirus (EpapGV) was sequenced and analyzed. Its circular double-stranded DNA genome is 119,082 bp in length and codes for 133 putative genes. It contains the 31 baculovirus core genes and a set of 19 genes that are GV exclusive. Seventeen ORFs were unique to EpapGV in comparison with other baculoviruses. Of these, 16 found no homologues in GenBank, and one encoded a thymidylate kinase. Analysis of nucleotide sequence repeats revealed the presence of 16 homologous regions (hrs) interspersed throughout the genome. Each hr was characterized by the presence of 1 to 3 clustered imperfect palindromes which are similar to previously described palindromes of tortricid-specific GVs. Also, one of the hrs (hr4) has flanking sequences suggestive of a putative non-hr ori. Interestingly, two more complex hrs were found in opposite loci, dividing the circular dsDNA genome in two halves. Gene synteny maps showed the great colinearity of sequenced GVs, being EpapGV the most dissimilar as it has a 20 kb-long gene block inversion. Phylogenetic study performed with 31 core genes of 58 baculoviral genomes suggests that EpapGV is the baculovirus isolate closest to the putative common ancestor of tortricid specific betabaculoviruses.

CONCLUSIONS

This study, along with previous characterization of EpapGV infection, is useful for the better understanding of the pathology caused by this virus and its potential utilization as a bioinsecticide.

Genomic sequencing and analyses of HearMNPV--a new Multinucleocapsid nucleopolyhedrovirus isolated from Helicoverpa armigera

Background

HearMNPV, a nucleopolyhedrovirus (NPV), which infects the cotton bollworm, Helicoverpa armigera, comprises multiple rod-shaped nucleocapsids in virion(as detected by electron microscopy). HearMNPV shows a different host range compared with H. armigera single-nucleocapsid NPV (HearSNPV). To better understand HearMNPV, the HearMNPV genome was sequenced and analyzed.

Methods

The morphology of HearMNPV was observed by electron microscope. The qPCR was used to determine the replication kinetics of HearMNPV infectious for H. armigera in vivo. A random genomic library of HearMNPV was constructed according to the “partial filling-in” method, the sequence and organization of the HearMNPV genome was analyzed and compared with sequence data from other baculoviruses.

Results

Real time qPCR showed that HearMNPV DNA replication included a decreasing phase, latent phase, exponential phase, and a stationary phase during infection of H. armigera. The HearMNPV genome consists of 154,196 base pairs, with a G + C content of 40.07%. 162 putative ORFs were detected in the HearMNPV genome, which represented 90.16% of the genome. The remaining 9.84% constitute four homologous regions and other non-coding regions. The gene content and gene arrangement in HearMNPV were most similar to those of Mamestra configurata NPV-B (MacoNPV-B), but was different to HearSNPV. Comparison of the genome of HearMNPV and MacoNPV-B suggested that HearMNPV has a deletion of a 5.4-kb fragment containing five ORFs. In addition, HearMNPV orf66, bro genes, and hrs are different to the corresponding parts of the MacoNPV-B genome.

Conclusions

HearMNPV can replicate in vivo in H. armigera and in vitro, and is a new NPV isolate distinguished from HearSNPV. HearMNPV is most closely related to MacoNPV-B, but has a distinct genomic structure, content, and organization.

Genomic sequence analysis of granulovirus isolated from the tobacco cutworm, Spodoptera litura

BACKGROUND

Spodoptera litura is a noctuid moth that is considered an agricultural pest. The larvae feed on a wide range of plants and have been recorded on plants from 40 plant families (mostly dicotyledons). It is a major pest of many crops. To better understand Spodoptera litura granulovirus (SpliGV), the nucleotide sequence of the SpliGV DNA genome was determined and analyzed.

METHODOLOGY/PRINCIPAL FINDINGS

The genome of the SpliGV was completely sequenced. The nucleotide sequence of the SpliGV genome was 124,121 bp long with 61.2% A+T content and contained 133 putative open reading frames (ORFs) of 150 or more nucleotides. The 133 putative ORFs covered 86.3% of the genome. Among these, 31 ORFs were conserved in most completely sequenced baculovirus genomes, 38 were granulovirus (GV)-specific, and 64 were present in some nucleopolyhedroviruses (NPVs) and/or GVs. We proved that 9 of the ORFs were SpliGV specific.

CONCLUSIONS/SIGNIFICANCE

The genome of SpliGV is 124,121 bp in size. One hundred thirty-three ORFs that putatively encode proteins of 50 or more amino acid residues with minimal overlap were determined. No chitinase or cathepsin genes, which are involved in the liquefaction of the infected host, were found in the SpliGV genome, explaining why SpliGV-infected insects do not degrade in a typical manner. The DNA photolyase gene was first found in the genus Granulovirus. When phylogenic relationships were analyzed, the SpliGV was most closely related to Trichoplusia ni granulovirus (TnGV) and Xestia c-nigrum granulovirus (XecnGV), which belong to the Type I-granuloviruses (Type I-GV).

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).

Biology and genomics of viruses within the genus Gammabaculovirus

Hymenoptera is a very large and ancient insect order encompassing bees, wasps, ants and sawflies. Fossil records indicate that they existed over 200 million years ago and about 100 million years before the appearance of Lepidoptera. Sawflies have been major pests in many parts of the world and some have caused serious forest defoliation in North America. All baculoviruses isolated from sawflies are of the single nucleocapsids phenotype and appear to replicate in midgut cells only. This group of viruses has been shown to be excellent pest control agents and three have been registered in Canada and Britain for this purpose. Sawfly baculoviruses contain the smallest genome of all baculoviruses sequenced so far. Gene orders among sequenced sawfly baculoviruses are co-linear but this is not shared with the genomes of lepidopteran baculoviruses. One distinguishing feature among all sequenced sawfly viruses is the lack of a gene encoding a membrane fusion protein, which brought into question the role of the budded virus phenotype in Gammabaculovirus biology.

Comparative genomic sequence analysis of novel Helicoverpa armigera nucleopolyhedrovirus (NPV) isolated from Kenya and three other previously sequenced Helicoverpa spp. NPVs

A newly cloned Helicoverpa armigera nucleopolyhedrovirus (HearNPV) from Kenya, HearNPV-NNg1, has a higher insecticidal activity than HearNPV-G4, which also exhibits lower insecticidal activity than HearNPV-C1. In the search for genes and/or nucleotide sequences that might be involved in the observed virulence differences among Helicoverpa spp. NPVs, the entire genome of NNg1 was sequenced and compared with previously sequenced genomes of G4, C1 and Helicoverpa zea single-nucleocapsid NPV (Hz). The NNg1 genome was 132,425 bp in length, with a total of 143 putative open reading frames (ORFs), and shared high levels of overall amino acid and nucleotide sequence identities with G4, C1 and Hz. Three NNg1 ORFs, ORF5, ORF100 and ORF124, which were shared with C1, were absent in G4 and Hz, while NNg1 and C1 were missing a homologue of G4/Hz ORF5. Another three ORFs, ORF60 (bro-b), ORF119 and ORF120, and one direct repeat sequence (dr) were unique to NNg1. Relative to the overall nucleotide sequence identity, lower sequence identities were observed between NNg1 hrs and the homologous hrs in the other three Helicoverpa spp. NPVs, despite containing the same number of hrs located at essentially the same positions on the genomes. Differences were also observed between NNg1 and each of the other three Helicoverpa spp. NPVs in the diversity of bro genes encoded on the genomes. These results indicate several putative genes and nucleotide sequences that may be responsible for the virulence differences observed among Helicoverpa spp., yet the specific genes and/or nucleotide sequences responsible have not been identified.

Baculovirus: molecular insights on their diversity and conservation

The Baculoviridae is a large group of insect viruses containing circular double-stranded DNA genomes of 80 to 180 kbp. In this study, genome sequences from 57 baculoviruses were analyzed to reevaluate the number and identity of core genes and to understand the distribution of the remaining coding sequences. Thirty one core genes with orthologs in all genomes were identified along with other 895 genes differing in their degrees of representation among reported genomes. Many of these latter genes are common to well-defined lineages, whereas others are unique to one or a few of the viruses. Phylogenetic analyses based on core gene sequences and the gene composition of the genomes supported the current division of the Baculoviridae into 4 genera: Alphabaculovirus, Betabaculovirus, Gammabaculovirus, and Deltabaculovirus.

The genomic diversity and phylogenetic relationship in the family iridoviridae

The Iridoviridae family are large viruses (∼120–200 nm) that contain a linear double-stranded DNA genome. The genomic size of Iridoviridae family members range from 105,903 bases encoding 97 open reading frames (ORFs) for frog virus 3 to 212,482 bases encoding 211 ORFs for Chilo iridescent virus. The family Iridoviridae is currently subdivided into five genera: Chloriridovirus, Iridovirus, Lymphocystivirus, Megalocytivirus, and Ranavirus. Iridoviruses have been found to infect invertebrates and poikilothermic vertebrates, including amphibians, reptiles, and fish. With such a diverse array of hosts, there is great diversity in gene content between different genera. To understand the origin of iridoviruses, we explored the phylogenetic relationship between individual iridoviruses and defined the core-set of genes shared by all members of the family. In order to further explore the evolutionary relationship between the Iridoviridae family repetitive sequences were identified and compared. Each genome was found to contain a set of unique repetitive sequences that could be used in future virus identification. Repeats common to more than one virus were also identified and changes in copy number between these repeats may provide a simple method to differentiate between very closely related virus strains. The results of this paper will be useful in identifying new iridoviruses and determining their relationship to other members of the family.