The genome sequence of Pseudoplusia includens single nucleopolyhedrovirus and an analysis of p26 gene evolution in the baculoviruses
© Craveiro et al.; licensee BioMed Central. 2015
Received: 1 August 2014
Accepted: 4 February 2015
Published: 25 February 2015
Pseudoplusia includens single nucleopolyhedrovirus (PsinSNPV-IE) is a baculovirus recently identified in our laboratory, with high pathogenicity to the soybean looper, Chrysodeixis includens (Lepidoptera: Noctuidae) (Walker, 1858). In Brazil, the C. includens caterpillar is an emerging pest and has caused significant losses in soybean and cotton crops. The PsinSNPV genome was determined and the phylogeny of the p26 gene within the family Baculoviridae was investigated.
The complete genome of PsinSNPV was sequenced (Roche 454 GS FLX – Titanium platform), annotated and compared with other Alphabaculoviruses, displaying a genome apparently different from other baculoviruses so far sequenced. The circular double-stranded DNA genome is 139,132 bp in length, with a GC content of 39.3 % and contains 141 open reading frames (ORFs). PsinSNPV possesses the 37 conserved baculovirus core genes, 102 genes found in other baculoviruses and 2 unique ORFs. Two baculovirus repeat ORFs (bro) homologs, bro-a (Psin33) and bro-b (Psin69), were identified and compared with Chrysodeixis chalcites nucleopolyhedrovirus (ChchNPV) and Trichoplusia ni single nucleopolyhedrovirus (TnSNPV) bro genes and showed high similarity, suggesting that these genes may be derived from an ancestor common to these viruses. The homologous repeats (hrs) are absent from the PsinSNPV genome, which is also the case in ChchNPV and TnSNPV. Two p26 gene homologs (p26a and p26b) were found in the PsinSNPV genome. P26 is thought to be required for optimal virion occlusion in the occlusion bodies (OBs), but its function is not well characterized. The P26 phylogenetic tree suggests that this gene was obtained from three independent acquisition events within the Baculoviridae family. The presence of a signal peptide only in the PsinSNPV p26a/ORF-20 homolog indicates distinct function between the two P26 proteins.
PsinSNPV has a genomic sequence apparently different from other baculoviruses sequenced so far. The complete genome sequence of PsinSNPV will provide a valuable resource, contributing to studies on its molecular biology and functional genomics, and will promote the development of this virus as an effective bioinsecticide.
Baculoviruses are specific pathogens of the insect orders Lepidoptera, Diptera and Hymenoptera and exhibit rod-shaped nucleocapsids embedded in a crystalline protein matrix (occlusion bodies – OBs) composed of polyhedrin in nucleopolyhedroviruses (NPVs) and granulin in granuloviruses (GVs) [1-3]. The replication cycle of the baculoviruses is characterized by production of two viral phenotypes: occlusion derived viruses (ODVs) and budded viruses (BVs). These particles are genotypically identical, but they are morphologically and functionally distinct, with the BVs involved in systemic infection within host larvae (produced in an early phase of infection) and the ODVs involved in the horizontal transmission of the virus in the host population (produced in the late phase of infection) . The Baculoviridae family consists of four genera: Alphabaculovirus (lepidopteran-specific NPV), Betabaculovirus (lepidopteran-specific GV), Gammabaculovirus (hymenopteran-specific NPV) and Deltabaculovirus (dipteran-specific NPV) [2,5,6]. Alphabaculoviruses can be further divided into Groups I and II based on DNA sequence data and differences in BVs, where the envelope fusogenic protein in Group I is GP64 and in Group II is the fusion (F) protein [7-10].
So far, 64 complete baculovirus genomes are present in GenBank, including many of the Alphabaculoviruses (45) followed by 15 Betabaculoviruses, 3 Gammabaculoviruses and the Culex nigripalpus Deltabaculovirus (CuniNPV) (http://www.ncbi.nlm.nih.gov/genomes/GenomesGroup.cgi?taxid=10442). Baculovirus genomes range in size from 81.7 (Neodiprion lecontei nucleopolyhedrovirus, NeleNPV) to 178.7 kbp (Xestia c-nigum granulovirus, XcGV) with GC content below 50% and containing from 89 (NeleNPV) to 183 (Pseudaletia unipuncta granulovirus, PsunGV) predicted ORFs . The gene diversity in baculoviruses has been estimated to be about 900 genes, among which 37 (core genes) may play essential biological functions in the replication cycle . The common genomic features of the Baculoviridae family include large double-stranded circular DNA, bidirectionally oriented open reading frames (ORFs) which are distributed on both DNA strands, 37 genes common to all baculoviruses (core genes), promoters that regulate the temporal cascade of gene expression and viral genome replication in the host cell nucleus .
The soybean looper, Chrysodeixis includens (syn., Pseudoplusia includens) (Walker, 1858) (Lepidoptera: Noctuidae, Plusiinae) is a lepidopteran pest with restricted distribution in the Western Hemisphere, occurring from the northern United States to southern South America [14,15]. Soybean, cotton, beans, potatoes, tomatoes, tobacco, sunflower, lettuce, cauliflower, cabbage and okra are the most common crops attacked by C. includens [16-21]. However, the polyphagous C. includens was found feeding on 73 plant species from 29 different families in Brazil . Until 2003, Anticarsia gemmatalis was considered one of the most important pests on soybean and the baculovirus Anticarsia gemmatalis MNPV was widely used as a bioinsecticide on approximately two million hectares of soybeans . Recently, C. includens has begun to have an economic impact due to its population growth, causing significant losses in soybean production. Among other factors, this was attributed to a decline in natural enemies, which previously controlled the pest, and to development of resistance due to indiscriminate use of chemical pesticides in soybean fields . Other forms of control are therefore required, and for this, new baculoviruses may be strong candidates for the biocontrol of this emerging pest.
Pseudoplusia includens single nucleopolyhedrovirus (PsinSNPV) is a Group II Alphabaculovirus pathogenic to C. includens . Seven PsinSNPV (IA to IG) isolates collected on cotton and soybean crops from Guatemala and Brazil were reported to cause fatal infections in C. includens larvae . Evidence of significant genetic variations and different degrees of pathogenicity were observed among the isolates analyzed in our previous studies [24,25]. Other PsinNPV isolates, PsinNPV-USA and PsinNPV-GT, have been reported, but little is known about them .
The isolate PsinSNPV-IE was obtained from C. includens larvae collected on Brazilian soybean crops and was found to be one of the most virulent against C. includens among seven isolates analyzed . In this manuscript, we report the complete sequence and organization of the PsinSNPV-IE genome and speculate on the origin of the p26 gene within the Baculoviridae family by potentially distinct acquisition events. The analysis of the PsinSNPV genome will provide important information for a better understanding of its virulence, evolution and molecular biology. These findings may also contribute to the development of a PsinSNPV bioinsecticide for the control of C. includens.
Results and discussion
Nucleotide sequence and gene content of the PsinSNPV genome
Comparison of PsinSNPV with others Alphabaculoviruses
Characteristics of the PsinSNPV genome compared with other Alphabaculoviruses
Genome size (bp)
GC content (%)
Coding sequence (%)
Mean % aa ID with PsinSNPV
Homologs in PsinSNPV
ORFs unique to PsinSNPV
GenBank accession number
Replication, transcription and structural genes
The baculovirus genes are categorized based on their functions during the viral cycle as follows: DNA replication, RNA transcription, ODV and BV structural proteins or oral infectivity proteins . Baculovirus genome replication mechanisms are still not fully understood. Several studies have been developed to try to identify the genes responsible for DNA replication and translation. The essential DNA replication factors late expression factor 1 (lef-1), lef-2, lef-3, DNA polymerase (dnapol), p6.9, 38 k, helicase (hel) and immediate early 1 (ie-1) homologs are all present in the PsinSNPV genome. In addition, the PsinSNPV genome contains genes homologous to proliferating cell nuclear antigen (pcna), major early-transcribed protein 53 (me53), DNA binding protein (dbp), alkaline exonuclease (alkexo) and exon-0/ie-0, which were not identified in all baculoviruses but may influence viral DNA replication [11,27].
The AcMNPV transcription system is activated in two main stages. At first, lef-4, lef-8, lef-9 and p47 are transcribed to encode the 4 subunits of the viral RNA polymerase complex . This complex acts on gene transcription and mRNA processing, including capping and polyadenylation. Then, the transcription enhancers lef-5 and very late factor 1 (vlf-1) are transcribed . All these genes were also found in the PsinSNPV genome. In addition, some supposedly non-essential genes involved in AcMNPV transcription regulation are also present in the PsinSNPV genome: lef-6, lef-11, 39 K, lef-10 and protein kinase 1 (pk1).
The PsinSNPV genome has 28 known baculovirus genes coding for structural proteins. Genes for ODV and BV structural proteins include polyhedrin (polh), orf1629, pk1, occlusion derived virus envelope protein 18 (odv-e18), occlusion derived virus enveloped capsid protein 27 (odv-ec27), p10, viral protein 1054 (vp1054), few polyhedra protein/25 k (fp25k), desmoplakin, 41-kDa glycoprotein (gp41), telokinin-like peptide 20 (tlp20), viral protein 91 (vp91/p95), vp39, p33, odv-e25, p87/vp80, odv-ec43, odv-e66, p13, calyx/polyhedrin enveloped protein (calyx/pep), p24, per os infectivity factor 0 (p74/pif-0), pif-1, pif-2, pif-3, odv-e28/pif-4, odv-e56/pif-5 and fusion (f) protein. The six PIF genes which are components of ODVs and are involved in oral infectivity exhibited high sequence similarity to ChchNPV and TnSNPV PIFs. The genus Alphabaculovirus is divided into Groups I and II based on gene content, and in particular, the BVs fusion protein: GP64 and F protein, respectively. PsinSNPV, a Group II Alphabaculovirus, possesses the expected F protein homolog.
Nucleotide metabolism and DNA repair
Several Group II Alphabaculoviruses and Betabaculoviruses encode genes involved in nucleotide biosynthesis. PsinSNPV possesses the ribonucleotide reductase (RR) large (RR1) and small (RR2) subunits and the dUTPase protein. These RR proteins are enzymes involved in the formation of deoxyribonucleotides from ribonucleotides . The dUTPase protein is responsible for preventing incorporation of mutagenic dUTP into DNA [3,30]. Poly (ADP-ribose) polymerase (PARP) and poly ADP-ribose glycohydrolase (PARG) are enzymes involved in synthesis of ADP riboses that activate and recruit DNA repair enzymes [31-33]. Although it has been reported that all Group II genomes encode PARG homologs , the PsinSNPV genome is notable for its absence.
The CPD photolyase, encoded by the DNA photolyase (phr) gene, acts at cyclobutane pyrimidine dimers to repair ultraviolet (UV) -induced DNA damage. The phr gene was identified in ChchNPV [34,35], TnSNPV , Plusia acuta NPV (PlacNPV)  and Thysanoplusia orichalcea NPV-B9 (ThorNPV-B9) . Studies suggest that the phr gene is conserved in Group II Alphabaculoviruses that infect lepidopteran insects in the Plusiinae subfamily of the Noctuidae family . However, the phr gene was also identified in baculoviruses that infect insects of other subfamilies, such as Spodoptera litura GV (SpliGV) (subfamily Hadeninae) , Clanis bilineata NPV (ClbiNPV) (Sphingidae family) , Apocheima cinerarium NPV (ApciNPV) (Geometridae family)  and Ampelophaga rubiginosa NPV (AmpeNPV) (Sphingidae family) (unpublished, 2008) .
PsinSNPV belongs to a group where the phr gene is conserved and, as expected, its genome encodes a CPD photolyase protein (Psin68). The complete nucleotide sequence of the PsinSNPV phr gene is of 1,512 bp with GC % = 36.2%. The deduced PHR amino acid sequence of PsinSNPV, ChchNPV- PHR1, −PHR2, and TnSNPV were aligned, revealing that the PsinSNPV photolyase possesses high identity to TnSNPV and ChchNPV- PHR1 (Additional file 4). Previous studies showed that ChchNPV - PHR1 is not active when tested in an Escherichia coli photolyase deficient strain . The active copy (PHR-2) is distinct in its possession of two conserved tryptophan residues, which may be involved in an electron transfer mechanism . In the PsinSNPV photolyase protein, the tryptophan residues are replaced by histidine and tyrosine in positions 368 and 370 aa, respectively (Additional file 4). Therefore, both tryptophans are absent in the PsinSNPV photolyase, suggesting that this protein might not be active. The partial PHR amino acid sequences of PsinNPV –GT1 (EU401912); −GT2 [GenBank: EU682272], PsinNPV - USA [GenBank: EU401913] and PsinSNPV-IA to -IG isolates described in the literature [24,39] were aligned, where PsinNPV-GT1 and PsinNPV – USA isolates showed high similarity to the PsinSNPV-IA to -IG isolates. Interestingly, in contrast to other PsinNPV isolates reported so far, PsinNPV - GT2 possesses a tryptophan residue at position 368, which is thought to be essential for enzyme activity (data not shown). Further studies are needed to confirm and investigate the activity of the PsinSNPV - GT2 photolyase.
The auxiliary genes viral ubiquitin (vubi), viral cathepsin (v-cath), chitinase (chiA), 37-kDa glycoprotein (gp37), conotoxin (ctl), superoxide dismutase 29 (sod29), fibroblast growth factors (fgf), phosphotyrosine phosphatase (ptp), ecdysone glucose transferase (egt), actin rearrangement infectivity factor 1 (arif-1), inhibitor of apoptosis 2 (iap-2), iap-3 and p35/p49 were found in the PsinSNPV genome. The auxiliary genes are non-essential for DNA replication, translation or viral particle formation. However, these genes confer selective advantages to viruses as has been observed in homologs of PsinSNPV auxiliary genes described in the literature. The activity of the v-cath and chiA genes is notable in P. includens larvae infected with PsinSNPV, where the encoded proteins cause degradation and liquefaction of the host cadaver [41,42]. The fgf, ptp and egt genes were reported to be involved in host hyperactive behaviors, increasing larval motility and preventing the molt to extend insect life, respectively . The ptp gene was previously reported to be only present in Group I NPVs [3,43], however this gene is present in the PsinSNPV genome.
Homologous regions (hrs) are absent from the PsinSNPV genome
Homologous regions (hrs) are repeated sequences with an imperfect palindromic core that are distributed in the genome as singletons or arranged in tandem. These repeat sequences are present in baculovirus genomes and other closely related invertebrate viruses . These regions act as enhancers of early gene transcription in NPVs and may serve as origins of replication in NPVs and GVs . Homologous regions are a common feature found in genomes of the four genera of the Baculoviridae family. However, no typical baculoviral hrs were found in the PsinSNPV genome, which is also the case in Buzura supressaria NPV (BusuNPV) , ChchNPV, TnSNPV and Agrotis segetum GV (AgseGV) [35,36,44,45].
PsinSNPV bro genes
The PsinSNPV genome sequence contains two baculovirus repeated ORFs (bro genes), named according to their order in the genome: bro-a (ORF-33) and bro-b (ORF-69). The bro genes commonly occur in Alpha-, Beta- and Gammabaculoviruses, varying in number of copies and length among the viruses [46-50]. These genes were first reported from baculoviruses, but bro gene homologs were subsequently identified in other insect dsDNA viruses, such as entomopoxvirus and entomoiridovirus [51-53]. BRO proteins exhibit a highly conserved N-terminal DNA binding domain (BRO-N) in the first 100–150 aa and a variable C-terminal domain (BRO-C) [48,54]. The functions of BRO proteins are not clear, but were proposed to be involved in host DNA replication and/or transcriptional regulation and as a viral replication enhancer in the late phase [46,48,49,54]. Although a deletion of 425 bp (386–811) (~140 aa) is present in the PsinSNPV bro-b gene compared with the ChchNPV bro-b gene, the genes share 72% identity. The PsinSNPV bro-a gene showed higher similarity to the Lymantria xylina MNPV bro-m gene and the Mamestra brassicae MNPV bro-a gene with 58 and 53% identity, respectively. In contrast to the PsinSNPV BRO-B protein with one BRO-N domain, PsinSNPV BRO-A protein contains two BRO-N domains [Pfam: PF02498] at amino acid position 14 to 118 and 139 to 234. In addition, the PsinSNPV BRO-A protein contains a domain of unknown function DUF3627 [Pfam: PF12299] in amino acid position 334 to 423. Although PsinSNPV, ChchNPV and TnSNPV are closely related, their bro genes do not show high similarity.
ORFs unique to PsinSNPV
Two putative ORFs, Psin5 and Psin8, were found to be unique to the PsinSNPV genome. These ORFs do not show significant similarity to other previously described baculovirus ORFs and exhibit signature sequences that describe domains predicted by InterProScan 5 . Psin5 is predicted to encode a 172 amino acid (aa) protein with molecular weight of 18.73 kDa and shows low homology to a signal transducer and activator of transcription protein in the avian species, Pseudopodoces humilis [GenBank: XP_005533966] (%ID = 80%, cover = 40% and e-value = 3E-12). Psin8 is predicted to encode a 150 aa, 11.9 kDa protein, but shows no significant similarity to any genes in GenBank databases (P > 0.01). The TAAG late promoter motif, combined with a TATA early promoter (TATAAGG motif), was identified about 100 bp upstream of both the Psin5 and Psin8 start codons (5,031 and 7,219 nt, respectively). These promoters are thought to be transcribed both by the host RNA polymerase II and viral RNA polymerases , suggesting that these genes could be expressed both early and late in infection. A search was made for protein families, domains and functional sites found in transmembrane domains in the predicted Psin5 and Psin8 proteins. Using TMHMM Server v 2.0, transmembrane helices from amino acid position 79 to 101 and 148 to 170 in Psin5 and 69 to 91 in Psin8 hypothetical proteins were predicted (Additional file 5).
Two p26 homologs in PsinSNPV
Two p26 (Ac136) gene homologs were identified in the PsinSNPV genome. The function of the p26 gene is not well understood, but studies have shown that deletion of the AcMNPV p26 gene produced no differences in phenotype from wild-type AcMNPV in cells and in larvae [3,57,58]. However, a combined deletion of p26, p10 and p74 genes in AcMNPV resulted in polyhedra containing few or no virions . For this reason, p26 is thought to be required for optimal virion occlusion in the OBs.
Amino acid sequences used in phylogenetic analysis of the P26 amino acid sequence
Position (ORF/Signal peptide)
Autographa californica NPV
136 / No
Bombyx mandarina NPV
120 / No
Bombyx mori NPV
121 / No
Rachiplusia ou MNPV
129 / No
Antheraea pernyi NPV
127 / No
Anticarsia gemmatalis MNPV
132 / No
Choristoneura fumiferana DEF MNPV
130 / No
Choristoneura fumiferana MNPV
128 / No
7 / Yes
Choristoneura murinana NPV
22 / No
Choristoneura occidentalis NPV
21 / No
143 / Yes
Choristoneura rosaceana NPV
22 / No
145 / Yes
Epiphyas postvittana NPV
119 / No
Hyphantria cunea NPV
21 / No
Maruca vitrata NPV
104 / No
Orgyia pseudotsugata MNPV
132 / No
Plutella xylostella MNPV
132 / No
Thysanoplusia orichalcea NPV
133 / No
Adoxophyes honmai NPV
31 / Yes
Adoxophyes orana NPV
30 / Yes
Agrotis ipsilon MNPV
149 / Yes
100 / No
Agrotis segetum NPV
142 / Yes
94 / No
Apocheima cinerarium NPV
102 / Yes
40 / No
Buzura suppressaria NPV
8 / Yes
52 / No
Chrysodeixis chalcites NPV
19 / No
63 / No
Clanis bilineata NPV
19 / Yes
61 / No
Ectropis obliqua NPV
18 / No
52 / No
Euproctis pseudoconspersa NPV
51 / Yes
61 / No
Helicoverpa armigera MNPV
151 / Yes
101 / No
Helicoverpa armigera NPV
22 / Yes
Helicoverpa armigera NPVG4
22 / Yes
Helicoverpa armigera NPV NNg1
21 / Yes
Helicoverpa zea SNPV
21 / Yes
Hemileuca sp. NPV
22 / Yes
60 / No
Leucania separata NPV
20 / No
Lymantria dispar MNPV
40 / Yes
Lymantria xylina MNPV
36 / Yes
Mamestra brassicae MNPV
147 / Yes
99 / No
Mamestra configurata NPV A
158 / Yes
109 / No
Mamestra configurata NPV B
157 / Yes
108 / No
Orgyia leucostigma NPV
20 / No
62 / No
Spodoptera exigua MNPV
129 / Yes
87 / No
Spodoptera frugiperda MNPV
131 / Yes
86 / No
Spodoptera litura NPV II
135 / Yes
91 / No
Trichoplusia ni SNPV
19 / Yes
59 / No
The p26 gene is conserved in position, adjacent to the p10 gene in all the Alphabaculoviruses containing a single copy. The PsinSNPV p26a/ORF-20 and p26b/ORF-62 are positioned adjacent to the p10 gene and adjacent to the iap-2 gene, respectively. The copy adjacent to the iap-2 gene is also positionally conserved in all Group II Alphabaculoviruses containing two p26 copies. However, the second p26 copy in Group I Alphabaculoviruses is positioned adjacent to the ptp1 and ptp2 genes.
The acquisition of baculovirus genes may be the results of duplication and horizontal gene transfer by transposable element and/or homologous recombination. The second acquisition event probably occurred by horizontal gene transfer. In this case, the gene duplication hypothesis can be refuted, since the similarity between p26 copies in the same virus is low (less than 30% identity). Furthermore, clades I and II are clearly separated, indicating that the p26 copies did not originate from a common, recent ancestor. In the third acquisition event, there is high similarity between the second p26 copy of the Group I Alphabaculoviruses and first p26 copy of the Group II Alphabaculoviruses, grouping these genes in the same clade (clade IB). Therefore, the second p26 copy of the Group I Alphabaculoviruses may have been acquired from a Group II Alphabaculovirus.
The isoelectric point (pI) and molecular weight (MW) of the deduced P26 protein from Alphabaculoviruses with two p26 copies were calculated, and are presented in the context of their genomic positioning (Table 2). P26 protein from position 1 showed an average pI of 6.7 ± 0.76 and average MW of 30,526 ± 2,387; position 2 P26 proteins have an average pI of 8.2 ± 0.71 and average MW of 27,354 ± 1,306; and position 3 P26 proteins have an average pI of 9.5 ± 0.07 and an average MW of 30,113 ± 24. The mean scores were examined using the Student’s t-test and the pI average showed a significant difference (p < 0.05) between the three positions. The isoelectric point difference suggests that the P26 sequence in position 3 exhibits more basic amino acids than in those in positions 1 and 2.
The presence and location of signal peptide cleavage sites in P26A and P26B amino acid sequences from PsinSNPV were analysed using SignalP v.4.1. The P26A protein showed a signal peptide cleavage site at amino acid residues 1 to 21 and the cleavage site (IMS-T) between amino acids 21 and 22 (Additional file 6). However, the signal peptide cleavage site was absent in the P26B protein. Signal peptides direct the proteins to their proper cellular and extracellular locations. The export of proteins occurs via the secretory pathway, where proteins labeled by an N-terminal signal sequence are translocated across the cytoplasmic membrane, whereafter the N-terminal signal peptide is usually cleaved by an extracellular signal peptidase.
The signal peptide cleavage sites were predicted for other P26 proteins of Alphabaculoviruses with complete sequenced genomes and the results are shown in Table 2. The presence or absence of signal peptides in P26 proteins correlates with the clustering obtained in the phylogenetic analysis, where only P26 proteins from clade IB possess signal peptides. However, four sequences belonging to this clade, LeseNPV_ORF20, EcobNPV_ORF18, OrleNPV_ORF20 and ChchNPV_ORF19, do not possess a signal peptide. The presence or absence of the signal peptide may have led to the differential results found in predicted molecular weight and isoelectric points between the P26 proteins analyzed. The presence of a signal peptide in the first p26 copy of Group II Alphabaculoviruses and in the second p26 copy of Group I Alphabaculoviruses suggests that this domain was acquired from a common ancestor of these viruses. Although the function of P26 is not well understood, the signal peptide may lead to differences in activity of the clade IB proteins compared to the P26 from other clades, which warrants further investigation.
In summary, the complete PsinSNPV-IE genome sequence is apparently different from other baculoviruses sequenced so far. The genome does not contain the typical baculovirus hrs and contains two ORFs with predicted transmembrane domains that are unique to PsinSNPV. The PsinSNPV genome, however, exhibits high sequence similarities and co-linearity to the closely related ChchNPV and TnSNPV. The PsinSNPV genome contains two p26 copies and a phylogenetic analysis of P26 sequences of Alphabaculoviruses showed three potential acquisition events of these genes within the Baculoviridae family. One of the clades comprises P26 protein with a signal peptide, indicating a possible distinct function from other classes of P26 protein. However, further investigations are needed for a better understanding of this protein in baculoviruses. This research reports the first completely sequenced genome of a strain of PsinSNPV, a currently little known baculovirus. It is anticipated that this data will both promote advances in investigations of its molecular biology and gene function and accelerate its development as a biocontrol agent.
Virus and viral DNA extraction
The Pseudoplusia includens SNPV – IE isolate, donated by Dr. Flávio Moscardi, Embrapa Soja (Londrina-PR), was obtained from an infected C. includens larva collected on soybean from a farm in Iguaraçu - PR-Brazil in 2007, and has been deposited in the Invertebrate Virus Collection at Embrapa Genetic Resources and Biotechnology. This isolate is listed in the Brazilian AleloMicro Information System under accession code BRM 005106. Viral OBs from C. includens larval cadavers were purified by differential centrifugation according to procedures described by Maruniak . DNA was extracted from ODVs as described previously [24,25]. The quality of the extracted DNA was determined by 0.5% agarose gel electrophoresis and quantified using a Qubit v. 2.0 Fluorometer (Invitrogen) according to the manufacturer’s instructions.
DNA sequence determination
The genome DNAs of seven PsinSNPV isolates (IA to IG), which have been investigated in our laboratory [24,25], were sequenced using the shotgun approach and were performed using the 454 Roche GS FLX – Titanium instrument at the Federal District (DF, Brazil) High-Performance Genome Center. Raw reads were processed using Newbler v. 2.8 (Roche Applied Science) and Biopieces scripts were used to create the fastq files. FastQC was used for quality assessment and Coral v. 1.4  was used to correct sequencing errors. PrinSeq v. 0.20.3  (Preprocessing and Information of Sequences) was applied to trim low quality reads (Phred ≤ 20) and to remove short sequences (length ≤ 50 bp). An error probability of 0.1% was allowed and 0.27% of the overall reads allowed to contain the ambiguous base ‘N’. The Phred score was measured and the mean sequence quality >30 was estimated, exhibiting an accuracy of 99.9% (Additional file 1). De novo assembly of reads from all isolates were carried out together using the MIRA assembler v. 22.214.171.124  resulting in a single contig of 148,729 bp. This scaffold was used to map the trimmed reads from the isolate PsinSNPV-IE, resulting in a final assembly for this single isolate, with a minimum coverage of 30X, representing its complete genomic sequence. SNPS and indels present in the assembled sequence were observed, since the sequenced isolate was not plaque purified, representing natural genotypic variation within PsinSNPV-IE . The final sequence was the 50% majority rule consensus of the PsinSNPV-IE reads. The polyhedrin gene was identified and the PsinSNPV circular nucleotide sequence was determined. An in silico BamHI, EcoRI, HindIII and PstI endonuclease restriction map was constructed using Geneious R6 v. 6.0.5 (Biomatters, Auckland, New Zealand) and was compared with DNA restriction profiles of PsinSNPV-IE determined previously .
Sequence data bioinformatics analysis
ORF prediction was carried out with ORF Finder (National Center for Biotechnology Information -NCBI) and Geneious (v. 6.0.5.). ATG-initiated ORFs encoding more than 50 amino acids with minimal overlaps were selected. Relevant ORFs were aligned against ChchNPV and TnSNPV genomes, and PsinSNPV-IE gene models confirmed using Artemis software . The BLASTx algorithm  was used to annotate the predicted ORFs. Percentage identities between homologous genes were obtained by alignment of the proteins from whole genomes using the tBLASTn program . Global alignment of the PsinSNPV genome against other baculovirus genomes was performed and the syntenic map constructed using Mauve alignment v. 2.0 implemented in the Geneious v. 6.0.3 package [65,66]. A dot-plot analysis was applied to compare the PsinSNPV genome against ChchNPV, TnSNPV, Mamestra configurata (Maco) NPV-B and Autographa californica (Ac) MNPV using LBDotView v. 1.0 software .
Deduced protein sequences were also analysed using SignalP 4.1 Server  (http://www.cbs.dtu.dk/services/SignalP/) and TMHMM (TransMembrane prediction using Hidden Markov Models) Server v. 2.0 [69,70] for prediction of signal peptide cleavage sites and Transmembrane (TM) helices, respectively.
P26 phylogenetic analysis
PsinSNPV P26a and P26b amino acid sequences were aligned using MUSCLE v. 3.5 software  against the corresponding amino acid sequences from other baculoviruses with sequenced genome (Table 2). A statistical model-fitting approach was conducted using ProtTest  and the LG model  was selected as best-fit model for the P26 alignment. Bayesian phylogenetic inference (BPI) was conducted using MrBayes v. 3.0b4 . Because MrBayes does not support the LG model of evolution, likelihood settings were set to aamodel = mixed rates = invgamma, which allowed the best model of substitution to be selected as a parameter of the analysis . Five Markov chains were run for 600,000 generations (p < 0.01), sampling every 100 generations. The first 25% of the trees obtained in the analysis were discarded as burn-in before computing the consensus tree.
Availability of supporting data
The authors would like to thank Orzenil Bonfim da Silva Jr for initial advise regarding genome assembly and Debora P. Paula for valuable discussions of results and suggestions. This work was supported by the following Brazilian Agencies: EMBRAPA (Empresa Brasileira de Pesquisa Agropecuária), FAPDF/CNPq (Fundação de Apoio à Pesquisa do Distrito Federal/Conselho Nacional de Desenvolvimento Científico e Tecnológico) and CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior).
- Jehle JA, Blissard GW, Bonning BC, Cory JS, Herniou EA, Rohrmann GF, et al. On the classification and nomenclature of baculoviruses: a proposal for revision. Arch Virol. 2006;151(7):1257–66.View ArticlePubMedGoogle Scholar
- Herniou EA, Arif BM, Becnel JJ, Blissard GW, Bonning B, Harrison R, et al. Baculoviridae. In: King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ, editors. Virus Taxonomy. Oxford: Elsevier; 2012. p. 163–74.Google Scholar
- Rohrmann GF. Baculovirus molecular biology. 3rd ed. Bethesda: National Center for Biotechnology Information; 2013.Google Scholar
- Keddie BA, Aponte GW, Volkman LE. The pathway of infection of Autographa californica nuclear polyhedrosis virus in an insect host. Science. 1989;243(4899):1728–30.View ArticlePubMedGoogle Scholar
- Jehle JA, Lange M, Wang H, Hu Z, Wang Y, Hauschild R. Molecular identification and phylogenetic analysis of baculoviruses from Lepidoptera. Virology. 2006;346(1):180–93.View ArticlePubMedGoogle Scholar
- Carstens EB, Ball LA. Ratification vote on taxonomic proposals to the International Committee on Taxonomy of Viruses (2008). Arch Virol. 2009;154(7):1181–8.View ArticlePubMedGoogle Scholar
- Monsma SA, Oomens AG, Blissard GW. The GP64 envelope fusion protein is an essential baculovirus protein required for cell-to-cell transmission of infection. J Virol. 1996;70(7):4607–16.PubMed CentralPubMedGoogle Scholar
- Hefferon KL, Oomens AGP, Monsma SA, Finnerty CM, Blissard GW. Host cell receptor binding by baculovirus GP64 and kinetics of virion entry. Virology. 1999;258(2):455–68.View ArticlePubMedGoogle Scholar
- Pearson MN, Groten C, Rohrmann GF. Identification of the Lymantria dispar nucleopolyhedrovirus envelope fusion protein provides evidence for a phylogenetic division of the Baculoviridae. J Virol. 2000;74(13):6126–31.View ArticlePubMed CentralPubMedGoogle Scholar
- Westenberg M, Uijtdewilligen P, Vlak JM. Baculovirus envelope fusion proteins F and GP64 exploit distinct receptors to gain entry into cultured insect cells. J Gen Virol. 2007;88(12):3302–6.View ArticlePubMedGoogle Scholar
- Ferrelli ML, Berretta MF, Belaich MN, Ghiringhelli PD, Sciocco-Cap A, Romanowski V. The baculoviral genome. In: Garcia M, editor. Viral genomes - molecular structure, diversity, gene expression mechanisms and host-virus interactions. Rijeka, Croatia: InTech; 2012.Google Scholar
- Miele SA, Garavaglia MJ, Belaich MN, Ghiringhelli PD. Baculovirus: molecular insights on their diversity and conservation. Int J Evol Biol. 2011;2011:379–424.View ArticleGoogle Scholar
- Thumbi DK, Béliveau C, Cusson M, Lapointe R, Lucarotti CJ. Comparative genome sequence analysis of Choristoneura occidentalis Freeman and C. rosaceana Harris (Lepidoptera: Tortricidae) alphabaculoviruses. PLoS ONE. 2013;8(7):e68968.View ArticlePubMed CentralPubMedGoogle Scholar
- Kogan M. Dynamics of insect adaptations to soybean: impact of integrated pest management. Environ Entomol. 1981;10(3):363–71.View ArticleGoogle Scholar
- Alford AR, Hammond AM. Temperature modification of female sex pheromone release in Trichoplusia ni (Hübner) and Pseudoplusia includens (Walker) (Lepidoptera: Noctuidae). Environ Entomol. 1982;11(4):889–92.View ArticleGoogle Scholar
- Bottimer LJ. Notes on some Lepidoptera from eastern Texas. J Agric Res. 1926;33:797–819.Google Scholar
- Folsom JW. Notes on little-known insects. J Econ Entomol. 1936;29:282–5.View ArticleGoogle Scholar
- Wolcott GN. A revised annotated check-list of the insects of Puerto Rico. J Agric Univ Puerto Rico. 1936;20:1–627.Google Scholar
- Hensley SD, Newson LD, Chapin J. Observations on the looper complex of the Noctuidae subfamily Plusiinae. J Econ Entomol. 1964;57:1006–7.View ArticleGoogle Scholar
- Herzog DC, Todd JH. Sampling velvetbean carterpillar on soybean. In: Kogan M, Herzog DC, editors. Sampling methods in soybean entomology. New York: Springer-Verlag; 1980. p. 107–40.View ArticleGoogle Scholar
- Bueno RCOF, Parra JRP, Bueno AF, Haddad ML. Desempenho de tricogramatídeos como potenciais agentes de controle de Pseudoplusia includens Walker (Lepidoptera: Noctuidae). Neotrop Entomol. 2009;38:389–94.View ArticlePubMedGoogle Scholar
- Bernardi O, Malvestiti GS, Dourado PM, Oliveira WS, Martinelli S, Berger GU, et al. Assessment of the high-dose concept and level of control provided by MON 87701 X MON 89788 soybean against Anticarsia gemmatalis and Pseudoplusia includens (Lepidoptera: Noctuidae) in Brazil. Pest Manag Sci. 2012;68(7):1083–91.View ArticlePubMedGoogle Scholar
- Moscardi F, de Souza M, de Castro M, Lara Moscardi M, Szewczyk B. Baculovirus pesticides: present state and future perspectives. In: Ahmad I, Ahmad F, Pichtel J, editors. Microbes and microbial technology. New York, USA: Springer; 2011. p. 415–45.View ArticleGoogle Scholar
- Craveiro SR, Melo FL, Ribeiro ZMA, Ribeiro BM, Báo SN, Inglis PW, et al. Pseudoplusia includens single nucleopolyhedrovirus: genetic diversity, phylogeny and hypervariability of the pif-2 gene. J Invertebr Pathol. 2013;114(3):258–67.View ArticlePubMedGoogle Scholar
- Alexandre TM, Ribeiro ZMA, Craveiro SR, Cunha F, Fonseca IC, Moscardi F, et al. Evaluation of seven viral isolates as potential biocontrol agents against Pseudoplusia includens (Lepidoptera: Noctuidae) caterpillars. J Invertebr Pathol. 2010;105(1):98–104.View ArticlePubMedGoogle Scholar
- Xu F, Vlak JM, van Oers MM. Conservation of DNA photolyase genes in group II nucleopolyhedroviruses infecting Plusiinae insects. Virus Res. 2008;136(1–2):58–64.View ArticlePubMedGoogle Scholar
- Mikhailov VS, Okano K, Rohrmann GF. Baculovirus alkaline nuclease possesses a 5′–>3′ exonuclease activity and associates with the DNA-binding protein LEF-3. J Virol. 2003;77(4):2436–44.View ArticlePubMed CentralPubMedGoogle Scholar
- Guarino LA, Xu B, Jin J, Dong W. A virus-encoded RNA polymerase purified from baculovirus-infected cells. J Virol. 1998;72(10):7985–91.PubMed CentralPubMedGoogle Scholar
- van Oers MM, Vlak JM. Baculovirus genomics. Curr Drug Targets. 2007;8(10):1051–68.View ArticlePubMedGoogle Scholar
- Herniou EA, Olszewski JA, Cory JS, O’Reilly DR. The genome sequence and evolution of baculoviruses. Annu Rev Entomol. 2003;48(1):211–34.View ArticlePubMedGoogle Scholar
- Miwa M, Tanaka M, Matsushima T, Sugimura T. Purification and properties of glycohydrolase from calf thymus splitting ribose-ribose linkages of poly(adenosine diphosphate ribose). J Biol Chem. 1974;11:3475–82.Google Scholar
- Chen X, Ijkel WFJ, Tarchini R, Sun X, Sandbrink H, Wang H, et al. The sequence of the Helicoverpa armigera single nucleocapsid nucleopolyhedrovirus genome. J Gen Virol. 2001;82(1):241–57.PubMedGoogle Scholar
- Deng F, Wang R, Fang M, Jiang Y, Xu X, Wang H, et al. Proteomics analysis of Helicoverpa armigera single nucleocapsid nucleopolyhedrovirus identified two new occlusion-derived virus-associated proteins, HA44 and HA100. J Virol. 2007;81(17):9377–85.View ArticlePubMed CentralPubMedGoogle Scholar
- van Oers MM, Herniou EA, Usmany M, Messelink GJ, Vlak JM. Identification and characterization of a DNA photolyase-containing baculovirus from Chrysodeixis chalcites. Virology. 2004;330(2):460–70.View ArticlePubMedGoogle Scholar
- van Oers MM, Abma-Henkens MHC, Herniou EA, Groot JCW, Peters S, Vlak JM. Genome sequence of Chrysodeixis chalcites nucleopolyhedrovirus, a baculovirus with two DNA photolyase genes. J Gen Virol. 2005;86(7):2069–80.View ArticlePubMedGoogle Scholar
- Willis LG, Siepp R, Stewart TM, Erlandson MA, Theilmann DA. Sequence analysis of the complete genome of Trichoplusia ni single nucleopolyhedrovirus and the identification of a baculoviral photolyase gene. Virology. 2005;338(2):209–26.View ArticlePubMedGoogle Scholar
- Wang Y, Choi JY, Roh JY, Woo SD, Jin BR, Je YH. Molecular and phylogenetic characterization of Spodoptera litura granulovirus. J Microbiol. 2008;46(6):704–8.View ArticlePubMedGoogle Scholar
- Zhu SY, Yi JP, Shen WD, Wang LQ, He HG, Wang Y, et al. Genomic sequence, organization and characteristics of a new nucleopolyhedrovirus isolated from Clanis bilineata larva. BMC Genomics. 2009;10(1):91.View ArticlePubMed CentralPubMedGoogle Scholar
- Biernat MA, Ros VID, Vlak JM, van Oers MM. Baculovirus cyclobutane pyrimidine dimer photolyases show a close relationship with lepidopteran host homologues. Insect Mol Biol. 2011;20(4):457–64.View ArticlePubMedGoogle Scholar
- van Oers MM, Lampen MH, Bajek MI, Vlak JM, Eker APM. Active DNA photolyase encoded by a baculovirus from the insect Chrysodeixis chalcites. DNA Repair. 2008;7(8):1309–18.View ArticlePubMedGoogle Scholar
- Ohkawa T, Majima K, Maeda S. A cysteine protease encoded by the baculovirus Bombyx mori nuclear polyhedrosis virus. J Virol. 1994;68(10):6619–25.PubMed CentralPubMedGoogle Scholar
- Slack JM, Kuzio J, Faulkner P. Characterization of v-cath, a cathepsin L-like proteinase expressed by the baculovirus Autographa californica multiple nuclear polyhedrosis virus. J Gen Virol. 1995;76(5):1091–8.View ArticlePubMedGoogle Scholar
- Nguyen Q, Nielsen L, Reid S. Genome scale transcriptomics of baculovirus-insect interactions. Viruses. 2013;5(11):2721–47.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhu Z, Yin F, Liu X, Hou D, Wang J, Zhang L, et al. Genome sequence and analysis of Buzura suppressaria nucleopolyhedrovirus: a group II Alphabaculovirus. PLoS ONE. 2014;9(1):e86450.View ArticlePubMed CentralPubMedGoogle Scholar
- Hilton S, Winstanley D. Genomic sequence and biological characterization of a nucleopolyhedrovirus isolated from the summer fruit tortrix, Adoxophyes orana. J Gen Virol. 2008;89(11):2898–908.View ArticlePubMedGoogle Scholar
- Kang W, Suzuki M, Zemskov E, Okano K, Maeda S. Characterization of baculovirus repeated open reading frames (bro) in Bombyx mori nucleopolyhedrovirus. J Virol. 1999;73(12):10339–45.PubMed CentralPubMedGoogle Scholar
- Kuzio J, Pearson MN, Harwood SH, Funk CJ, Evans JT, Slavicek JM, et al. Sequence and analysis of the genome of a baculovirus pathogenic for Lymantria dispar. Virology. 1999;253(1):17–34.View ArticlePubMedGoogle Scholar
- Iyer LM, Koonin EV, Aravind L. Extensive domain shuffling in transcription regulators of DNA viruses and implications for the origins of fungal APSES transcription factors. Genome Biol. 2002;3:1–11.View ArticleGoogle Scholar
- Bideshi DK, Renault S, Stasiak K, Federici BA, Bigot Y. Phylogenetic analysis and possible function of bro-like genes, a multigene family widespread among large double-stranded DNA viruses of invertebrates and bacteria. J Gen Virol. 2003;84(9):2531–44.View ArticlePubMedGoogle Scholar
- Zhou JB, Li XQ, De-Eknamkul W, Suraporn S, Xu JP. Identification of a new Bombyx mori nucleopolyhedrovirus and analysis of its bro gene family. Virus Genes. 2012;44(3):539–47.View ArticlePubMedGoogle Scholar
- Afonso CL, Tulman ER, Lu Z, Oma E, Kutish GF, Rock DL. The genome of Melanoplus sanguinipes entomopoxvirus. J Virol. 1999;73(1):533–52.PubMed CentralPubMedGoogle Scholar
- Bawden AL, Glassberg KJ, Diggans J, Shaw R, Farmerie W, Moyer RW. Complete genomic sequence of the Amsacta moorei entomopoxvirus: analysis and comparison with other poxviruses. Virology. 2000;274(1):120–39.View ArticlePubMedGoogle Scholar
- Jakob NJ, Müller K, Bahr U, Darai G. Analysis of the first complete DNA sequence of an invertebrate iridovirus: coding strategy of the genome of Chilo iridescent virus. Virology. 2001;286(1):182–96.View ArticlePubMedGoogle Scholar
- Zemskov EA, Kang WK, Maeda S. Evidence for nucleic acid binding ability and nucleosome association of Bombyx mori nucleopolyhedrovirus BRO proteins. J Virol. 2000;74(15):6784–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30(9):1236–40.View ArticlePubMed CentralPubMedGoogle Scholar
- Xing K, Deng R, Wang J, Feng J, Huang M, Wang X. Analysis and prediction of baculovirus promoter sequences. Virus Res. 2005;113(1):64–71.View ArticlePubMedGoogle Scholar
- Simòn O, Williams T, Caballero P, Possee RD. Effects of Acp26 on in vitro and in vivo productivity, pathogenesis and virulence of Autographa californica multiple nucleopolyhedrovirus. Virus Res. 2008;136(1–2):202–5.View ArticlePubMedGoogle Scholar
- Wang L, Salem TZ, Campbell DJ, Turney CM, Kumar CMS, Cheng XW. Characterization of a virion occlusion-defective Autographa californica multiple nucleopolyhedrovirus mutant lacking the p26, p10 and p74 genes. J Gen Virol. 2009;90(7):1641–8.View ArticlePubMedGoogle Scholar
- Maruniak JE. Baculovirus structural proteins and protein synthesis. In: Granados RR, Federici BA, editors. The biology of baculoviruses, vol. 1. Boca Raton: CRC; 1986. p. 129–46.Google Scholar
- Salmela L, Schröder J. Correcting errors in short reads by multiple alignments. Bioinformatics. 2011;27(11):1455–61.View ArticlePubMedGoogle Scholar
- Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27(6):863–4.View ArticlePubMed CentralPubMedGoogle Scholar
- Chevreux B, Pfisterer T, Drescher B, Driesel AJ, Müller WEG, Wetter T, et al. Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs. Genome Res. 2004;14(6):1147–59.View ArticlePubMed CentralPubMedGoogle Scholar
- Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, et al. Artemis: sequence visualization and annotation. Bioinformatics. 2000;16(10):944–5.View ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.View ArticlePubMed CentralPubMedGoogle Scholar
- Darling ACE, Mau B, Blattner FR, Perna NT. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14(7):1394–403.View ArticlePubMed CentralPubMedGoogle Scholar
- Aaron DE, Mau B, Perna NT. Progressive Mauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS ONE. 2010;5(6):e11147.View ArticleGoogle Scholar
- Huang Y, Zhang L. Rapid and sensitive dot-matrix methods for genome analysis. Bioinformatics. 2004;20(4):460–6.View ArticlePubMedGoogle Scholar
- Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8(10):785–6.View ArticlePubMedGoogle Scholar
- Sonnhammer ELL, von Heijne G, Krogh A. A Hidden Markov model for predicting transmembrane helices in protein sequences. In: Proceedings of the 6th International Conference on Intelligent Systems for Molecular Biology. Menlo Park, CA: AAAI Press; 1999.Google Scholar
- Krogh A, Larsson B, von Heijne G, Sonnhammer ELL. Predicting transmembrane protein topology with a Hidden Markov model: application to complete genomes. J Mol Biol. 2001;305(3):567–80.View ArticlePubMedGoogle Scholar
- Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Abascal F, Zardoya R, Posada D. ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005;21(9):2104–5.View ArticlePubMedGoogle Scholar
- Le SQ, Gascuel O. An improved general amino acid replacement matrix. Mol Biol Evol. 2008;25(7):1307–20.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17(8):754–5.View ArticlePubMedGoogle Scholar
- Nylander J. Testing models of evolution—MrModeltest version 1.1 b. Computer program and documentation distributed by author, website: http://www.ebc.uu.se/systzoo/staff/nylander.html; 2002.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.