Skip to main content
  • Research article
  • Open access
  • Published:

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



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.


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.


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.


Cassava (Manihot esculenta) is the basic source for dietary energy of 500 million people in tropical and subtropical areas of Africa, Asia, and Latin America [1]. In Brazil, the hornworm Erinnyis ello ello (Lepidoptera: Sphingidae) is one of the most important pests [2] occurring throughout the year and greatly impacting cassava production [3, 4]. This pest has been observed in 35 plant species, especially in the Euphorbiaceae family [5, 6]. In large infestations, the cassava pest may reduce by 50% the roots yield. In the 1980s, a naturally occurring baculovirus was isolated from this pest and applied as a bio-pesticide in Brazil [5]. The biological control program has proven to be safe and economical [5, 6]. However, genomic and structural information about this virus is lacking.

The Baculoviridae is a family of insect viruses with circular double-stranded genomic DNA [79] that have been successfully applied for the control of agricultural and forest pests [10]. So far, Alpha and Betabaculovirus are the most studied baculovirus genera; both infect Lepidoptera [7]. The infection is initiated when larvae feed on foliage contaminated with orally infectious occlusion bodies (OBs) [11] that release occlusion derived-virions (ODVs) in the midgut [12]. Early after primary midgut epithelial cell infection, budded virions (BV) are produced and cause systemic infection. Infection symptoms include cuticle discoloration, movement loss, and incapability for feeding [13, 14].

Few full-length betabaculovius genome sequences are available compared to those from Alphabaculovirus and none of them was isolated from sphingid host. In this context, identification and sequencing of virus species from different lepidopteran families will provide a wider empirical database to help understand baculovirus evolution [15, 16]. Here, we presented the morphological characterization, the complete genome sequence, and the phylogenetic analyses of the natural cassava hornworm pesticide, the first completely sequenced betabaculovirus isolated from a sphingid host.

Results and discussion

Virus characterization and genome features

A naturally occurring baculovirus was isolated from dead cassava hornworm (E. ello ello) caterpillars in crops from the South of Brazil in 1986. As shown in Figure 1A, the larvae is usually found hanged in cassava apical leaves, which is a characteristic symptom of the baculovirus infections [17]. Neither cuticle melanization nor post-mortem melting phenotypes were observed among the caterpillar cadavers, an attribute which probably facilitated virus collection and use for pesticide production as previously observed in another baculovirus (Anticarsia gemmatalis multiple nucleopolyhedrovirus - AgMNPV) [10]. Ultrastructural analyses revealed a granular OB with irregular form and size (Figure 1B) containing single rod-shaped nucleocapsid (Figure 1C). Both of these structural features, i.e. granular form and nucleocapsid shape, are typical of viruses from the genus Betabaculovirus [8, 18] and thus, we named it Erinnyis ello granulovirus (ErelGV) isolate Br-S86 (Brazil/South/1986). Two other cassava hornworm-isolated granuloviruses were previously reported, one isolated in Colombia [19] and another from an undisclosed geographical source [20]. Restriction endonuclease profile analyses (Figure 1D) suggest that the Brazilian and the Colombian viruses (previously published in [19]) are either variants of the same species or are distinct species infecting the same host. However, the absence of sequence data from the latter prevents establishment of any phylogenetic relationship.

Figure 1
figure 1

Erinnyis ello granulovirus (ErelGV) infection and virus characterization. (A) Cassava hornworm cadaver found hanging in the field due to terminal baculovirus infection (Source: José Osmar Lorenzi). (B) Scanning and (C) transmission electron micrographs reveal granular occlusion bodies containing singly embedded rod-shaped nucleocapsid (nc) (scale bars = 0.5 μm). (C) Restriction enzyme profile of Brazilian isolate genomic DNA. Agarose gel electrophoresis-resolved DNA fragments digested with HindIII (lane 1), EcoRI (lane 2), BamHI (lane 3).

We sequenced the genome of ErelGV, the first completely sequenced sphingid host-isolated betabaculovirus (Genbank accession number KJ406702). The genome is 102,759 bp long with a G + C content of 38.7% (Table 1). We found 130 putative genes coding for polypeptides of at least 50 amino acid residues. Additional file 1: Table S1 summarizes the ErelGV genes and compares each predicted protein sequence with its orthologs in other baculoviruses. Eight of these were shown to be unique (ErelOrf-11, ErelOrf-15, ErelOrf-27, ErelOrf-53, ErelOrf-59, ErelOrf-70, ErelOrf-90, ErelOrf-102), and all of them are peptides with no significant similarity to any other sequence in GenBank. All 37 Baculoviridae core genes were found and no typical homologous regions (hrs) were detected. However, we identified five putative homologous regions (hrs)/repeat regions lacking typical alphabaculovirus hr palindromes. This feature is also found in both Choristoneura occidentalis granulovirus (ChocGV) and Pieris rapae granulovirus (PiraGV) genomes. As observed in ChocGV [21], ErelGV lacks both gp37 and exon0, which was previously predicted for being shared among all Alpha and Betabaculovirus [22].

Table 1 All species from the genus Betabaculovirus completely sequenced to date

Phylogenetic analysis

In order to better understand the evolutionary history of ErelGV and the genus Betabaculovirus, we carried out a maximum likelihood phylogenetic analysis using the 37 baculovirus core gene alignment from all baculovirus genome available. ErelGV clustered with ChocGV and both viruses share the same ancestor with PiraGV isolates (Figure 2). Since the Chinese and Korean PiraGV isolates are very similar to each other (99.5%), we have included only the Chinese isolate in our analyses. Using Mauve alignment [34], we found that ChocGV and PiraGV genomes have respectively 38.5% and 34.5% of global pairwise identity when compared to ErelGV genome. Additionally, our phylogenetic analyses did not find support for Betabaculovirus division in two clades (A and B), as described previously using neighbor joining clustering method [25, 28]. Phylogenetic relationships in Baculoviridae, in particular in the genus Betabaculovirus, are difficult to discern due to the limited number of sequenced genomes available (Table 1). Therefore, we further evaluated ErelGV phylogenetic relationships using granulin, lef-8, and lef-9 partial gene dataset as previously carried out [20, 35] (28 partial sequences), but including new sequences publicly available (seven sequences from completely sequenced baculovirus) totalizing 35 granulovirus sequences. This analysis revealed that ErelGV isolate Br-S86 is closely related to another ErelGV (also called EeGV) from the Steinhaus collection [20] and that both are closer to Andraca bipunctata granulovirus (AnbiGV) (data not shown).

Figure 2
figure 2

Maximum likelihood tree for Betabaculovirus . The phylogenetic inference was based on the concatenated amino acid sequences of the 37 core genes identified in all complete baculovirus genome sequences. We collapsed all the Gammabaculovirus and Alphaphabaculovirus. The CuniNPV was used as root. ErelGV (boldface) clustered with ChocGV and both were closely related to PiraGV isolates.

Betabaculovirusgene comparison

We performed BLAST comparisons between ErelGV and all other full betabaculovirus genomes available in Genbank using the CGView Comparison Tool [36] and CIRCOS [37]. As shown in Figure 3, most of the ErelGV-encoded ORFs are conserved among all betabaculovirus, but protein similarity varies widely across the species. Some structural proteins, such as granulin and the per os infectivity factors (PIFs), were the most conserved genes. Conversely, F protein, the major Betabaculovirus envelope fusion protein (EFP, encoded by ErelOrf-28) and matrix metalloproteinase (MMP, a stromelysin-1-like protein, encoded by ErelOrf-39) were particularly variable despite of both being present in every betabaculovirus sequenced to date. The EFP is essential for cell-to-cell movement and systemic virus spread [7]. GP64 is the EFP found in Group I Alphabaculovirus and all orthologs are closely related to each other (81% of protein sequence identity), whereas the F protein, found in both Alpha and Betabaculovirus [38], is very diverse (20 to 40% sequence identity). Interestingly, deletion of the gp64 or f protein genes is lethal for BV propagation in Autographa californica multiple nucleopolyhedrovirus (AcMNPV) [39] and Helicoverpa armigera nuclepolyhedrovirus (HaNPV) [40, 41], respectively. The deficiency can be rescued by efp homologs from many different viruses in the case of AcMNPV [42], but the opposite is not true; AcMNPV gp64 is not able to completely rescue an efp-deleted HaNPV. However, it is not clear why the F protein from Plutella xylostella granulovirus (PlxyGV) is not able to rescue the infectivity of gp64-null AcMNPV [42] but that from AgseGV can. PlxyGV causes systemic infection to the diamondback moth P. xylostella (Plutellidae) larvae [43] and AgseGV infects the cutworm A. segetum (Noctuidae) [44]. Thus, the betabaculovirus EFP variability might reflect the cell machinery adjustment at the insect family level considering that AcMNPV infects caterpillar from the same insect family of A. segetum. A second highly variable gene, MMP, is a proteinase able to produce a distinct pattern of melanization in Bombyx mori larvae infected with the Xestia c-nigrum granulovirus (XecnGV) metalloproteinase-expressing Bombyx mori nucleopolyhedrovirus [45]. The enzyme is thought to enhance, replace, or act synergistically with proteins from virus or host playing an important role in the virus spread [46]. This variability is not unexpected since granulovirus genomes vary in content with respect to the presence or absence of the proteases cathepsin and enhancin genes and also the chitinase gene [12, 4547].

Figure 3
figure 3

Gene comparison of ErelGV genome and all completely sequenced betabaculoviruses available in Genbank. CDS identities were acquired by BLAST analysis and ranked from 0 to 100%. From the outermost ring: ChocGV, PiraGV-E3, PiraGV-China, ClanGV, CaLGV, CrleGV, CypoGV, AdorGV, PhopGV, EpapGV, AgseGV, PlxyGV, PsunGV, XecnGV, HearGV, and SpliGV-K1. For this representation, gene synteny is not taken into account. CDS that were absent in the ErelGV genome but present in the query sequences were not displayed. To prevent the missing of known homologues, like p6.9 and odv-e18 (asterisk), all the low identity hits (bellow 20%) were plotted as well. Unique genes are shown in red, core genes are in boldface, and Densovirus-related genes are shown in green.

Lack of cathepsin and chitinasegenes

ErelGV lacks cathepsin and chitinase genes, despite of their importance for promoting baculovirus horizontal transmission [48]. This feature can explain the integrity of caterpillar flesh and light color after death (Figure 1A). Other betabaculovirus genomes also lack both enzymes: complete deletion in ChocGV [21], Adoxophyes orana granulovirus (AdorGV) [23], Phthorimaea opercullela granulovirus (PhopGV) (unpublished), PlxyGV [31] and Spodoptera litura granulovirus (SpliGV) [32]; Cryptophlebia leucotreta granulovirus (CrleGV) [26chitinase has an interruption; and in Helicoverpa armigera granulovirus (HearGV) [29] only cathepsin is absent. Interestingly, most of these deletions seem to have occurred independently of each other within Betabaculovirus (data not shown), aside from ChocGV and ErelGV in which is strongly supported an ancestral lacking. Thus, it is reasonable to expect that AnbiGV, the closest relative to ErelGV, might also lack both cathepsin and chitinase. Taken together, these results reinforce the notion that both genes are most likely non-essential for the persistence of baculoviruses in the environment. Conversely, previous work from our research team has shown that introduction of cathepsin and chitinase from Choristoneura fumiferana defective nucleopolyhedrovirus into AgMNPV (which naturally lacks both genes) increases pathogenicity and occlusion body production relative to the wild type virus [49].


ErelOrf-5 codes for a nucleotide metabolism-related gene homologous to Orgyia pseudotsugata multiple nucleopolyhedrovirus (OpMNPV) Orf-31. The gene seems to be composed of a fusion between two distinct ORFs; the C-terminal portion is related to a baculovirus thymidylate kinase-like gene and the N-terminal portion is related to several dUTPase-like genes. The thymidylate kinase enzyme catalyzes a critical step in the biosynthesis of deoxythymidine triphosphate [50]. dUTPase catalyses dUTP dephosphorylation to generate dUMP [51]. High levels of dUTP can be deleterious for virus genomic DNA replication since dTTP can be substituted for dUTP during DNA synthesis [52]. A high dUTP/dTTP ratio promotes uracil incorporation into DNA. Uracils in DNA are then targeted by uracil DNA glycosylase and excised, leading to futile repair cycles and DNA breakage and or translesional DNA synthesis [53, 54]. Nucleotide metabolism-related enzyme acquisition is common in baculoviruses [28] and could avoid this deleterious response by decreasing the dUTP/dTTP ratio, however how these genes alter the virus fitness is not clear [55].

The he65-like and p43-like genes

The ErelGV genome contains homologues of the he65 and p43 genes. Homologues of he65 are harbored by several alphabaculoviruses, four betabaculoviruses (Agrotis segetum granulovirus (AgseGV), HearGV, Pseudaletia unipuncta granulovirus (PsunGV), and Xestia c-nigrum granulovirus (XecnGV)), and two betaentomopoxviruses (Amsacta moorei entomopoxvirus - AMV and Mythimna separate entomopoxvirus - MSV). This gene is a member of a distinct RNA ligase family related to the T4 RNA ligase gp63-like gene and is present in all the domains of life (Bacteria, Archaea, and Eukarya) [7, 56]. The alignment of baculovirus and entomopoxvirus he65-like genes revealed large, independent, and recurrent deletions in the C-terminal region (data not shown), which contain five nucleotidyl transferase motifs [56]. The amino-terminal region was highly conserved although no previously characterized motifs were present. We performed a phylogenetic reconstruction based on this conserved domain. The he65 reconstruction revealed distinct horizontal gene transfer (HGT) events from Alphabaculovirus to Betabaculovirus and Betaentomopoxvirus (Figure 4A). Betabaculovirus likely endured two independent acquisitions from Group II Alphabaculovirus in distinct genomic regions: (i) a synapomorphic introduction for HearGV, PsunGV, and XecnGV (Figure 4A, orange rectangle); and (ii) an additional gain for AgseGV (Figure 4A, dark orange rectangle). Importantly, support for AgseGV branch is low. However, the genomic context of the gene is conserved among HearGV, PsunGV, and XecnGV but not in AgseGV (data not shown), reinforcing our hypothesis that two independent introductions occurred. Likewise, Betaentomopoxvirus homologues were probably acquired from Group II Alphabaculovirus (Figure 4, yellow rectangle). Remarkably, ErelGV is the first Betabaculovirus with a he65-like gene (ErelOrf-36 - Figure 4A, green rectangle) acquired from Group I Alphabaculovirus. It is not clear whether C-terminal deleted he65 remains functional in baculovirus. However the maintenance of the amino-terminal region indicates that this gene region is under positive selection pressure.

Figure 4
figure 4

Phylogeny of he65 and p43 reveals horizontal gene transfer in ErelGV from Alphabaculovirus . (A) The maximum likelihood (ML) tree was inferred using the conserved amino-terminal region alignment of he65-like gene for 36 baculoviruses and two entomopoxviruses. Circles indicate the presence (blue) or absence (red) of the carboxy-terminal region. The postulated horizontal gene transfer (HGT) events are highlighted for Betabaculovirus (light and dark orange), Betaentomopoxvirus (yellow), and ErelGV (green). (B) ML-Phylogenetic reconstruction for p43-like gene found in ErelGV genome. The trees are midpoint rooted for purposes of clarity.

Furthermore, we found in ErelGV genome a p43-like gene (ErelOrf-105) whose homologues were found only in baculovirus species from the genus Alphabaculovirus (Figure 4B) with conserved amino acid sequence and position in the genome [7]. Deletion of p43 in AcMNPV does not affect virus replication in cell culture and the reason for gene acquisition and preservation is not clear [57]. Two hypotheses can be raised for p43 introduction in ErelGV: (i) ErelGV acquired the p43-like gene from Group I Alphabaculovirus, specifically from AcMNPV-related viruses; or (ii) ErelGV acquired from Group II Alphabaculovirus, specifically from a baculovirus (e.g. Clanis bilineata nucleopolyhedrovirus - ClbiNPV [58]) during co-infection of a sphingid host.

Acquisitions of Densovirus-related genes in Betabaculovirus

ErelOrf-57 and ErelOrf-100 are homologues to a non-structural Densovirus gene. Densovirus-related genes were previously described in two betabaculoviruses, ChocGV (ChocOrf-25) [21] and CrleGV (CrleOrf-9) [26], and one gammabaculovirus (baculovirus infective to hymenoptera), Neodiprion lecontei nucleopolyhedrovirus (NeleNPV - NeleOrf-81) [59]. The latter did not match the other two homologues (data not shown), suggesting these genes resulted from at least two HGT events between densoviruses and baculoviruses. Despite the limited number of Densovirus genomes available, we performed a phylogenetic analysis to help understand the origins of Betabaculovirus homologues. We found that the genes were dispersed over the phylogenetic tree, suggestive of multiple HGT events. As shown in Figure 5, Betabaculovirus homologues did not form a monophyletic cluster. To further substantiate our findings, we compared the likelihood of the observed tree to that estimated assuming a Betabaculovirus monophyletic clade (single-HGT event). Indeed, the likelihood ratio test rejected the monophyletic hypothesis favoring the multiple-HGT scenario, which was also supported by the distinct genomic context observed for the homologous betabaculovirus genes (data not shown). Moreover, both ErelOrf-57 and ErelOrf-100 form a well-supported clade, indicating that they probably represent a gene duplication event during ErelGV evolution.

Figure 5
figure 5

Densovirus -related genes in betabaculovirus and phylogenetic relationship. ML tree was inferred using the alignment of ErelOrf-57 and ErelOrf-100 from ErelGV with non-structural protein (NS) from Bombyx mori densovirus 2 and 3 (BmbDENSV-2 / Genbank YP_007714627.1 and BmDENSV-3/Genbank YP_007714627), NS3 from Diatraea saccharalis densovirus (DisaDENSV/Genbank NP_046812.1), Mythimna loreyi densovirus (MyloDENSV/Genbank NP_958098.1), Helicoverpa armigera densovirus (HaDENSV/Genbank AFK91982.1), Galleria mellonella densovirus (GameDENSV/Genbank NP_899649.1), Junonia coenia densovirus (JucoDENSV/Genbank AGO32182.1), and Pseudoplusia includens densovirus (PsinDENSV/Genbank YP_007003822.1), Orf-25 from ChocGV, and Orf-9 from CrleGV. The tree is midpoint rooted for purposes of clarity only. We hypothesized gene duplication for both ErelGV genes (boldface).


ErelGV is a new betabaculovirus species closely related to ChocGV and PiraGV isolates. Its genome encodes 130 ORFs, eight of which are unique. We found evidence suggesting horizontal gene transfers from Alphabaculovirus and Densovirus to Betabaculovirus. The he65-like gene was independently acquired tree times from Alphabaculovirus. We found a dUTPase-like gene homologous to OpMNPV Orf-31 and two Densovirus-related genes. The contribution of these genes to baculovirus fitness is not clear and is being experimentally tested in our lab. Importantly, recombination events between insect viruses from unrelated families and genera might drive baculovirus genomic evolution.

Nucleotide sequence accession number

The ErelGV genome sequence was submitted to GenBank under accession number KJ406702.

Availability of supporting data

The complete ErelGV genome sequence has been submitted to GenBank (accession number KJ406702). All supporting data is included as additional files.


Virus purification

Insect cadavers of the hornworm E. ello ello with baculovirus infection symptoms were collected in cassava crops in the South of Brazil (Itajaí, Santa Catarina) in 1986. They were kindly provided by Dr. Renato Arcanjo Pegoraro (EPAGRI). The cadavers were kept in the freezer and later used for OBs purification. Insect cadavers were homogenized with ddH2O (w/v), filtered through three layers of gauze, and centrifuged at 7,000 x g for 10 min. The pellet was resuspended in 0.5% (w/v) SDS and again centrifuged at 7,000 x g for 10 min. The dilution and centrifugation steps were repeated four times, and the final pellet was washed in 0.5 M NaCl. The pellet was resuspended in ddH2O, loaded onto a continuous 40-65% sucrose gradient, and centrifuged at 104.000 x g for 40 min at 4°C. The OB band was collected, diluted 4-fold in ddH2O, and centrifuged at 7,000 x g for 15 min at 4°C.

Electron microscopy

For scanning electron microscopy (SEM), 100 μl of the OB-containing solution (109 OBs/ml) were incubated with 300 μl of acetone at 25°C for 1 hour. The solution was loaded in a metallic stub, dried overnight at 37°C, coated with gold in a Sputter Coater (Balzers) for 3 min, and observed in a scanning electron microscope Jeol JSM 840 A at 10 kV. For transmission electron microscopy (TEM) pellets of purified granules were fixed in Karnovsky fixative (2.5% glutaraldehyde, 2% paraformaldehyde, in 0.1 M, cacodylate buffer, pH 7.2) for 2 h, post-fixed in 1% osmium tetroxide in the same buffer for 1 h and then stained en bloc with 0.5% aqueous uranyl acetate, dehydrated in acetone, and embedded in Spurr’s low viscosity embedding medium. The ultrathin sections were contrasted with 2% uranyl acetate and observed in a ZEISS TEM 109 at 80 kV.

Genomic DNA restriction analyses

Purified granules (109 OBs/ml) were dissolved in an alkaline solution and used to extract DNA according to O’Reilly et al. [60]. The quantity and quality of the isolated DNA were determined by electrophoresis on 0.8% agarose (data not shown). The viral DNA (1–2 μg) was individually cleaved with the restriction enzymes HindIII, EcoRI, and BamHI (Promega) according to manufacturer’s instructions. The DNA fragments generated were analyzed by 0.8% agarose gel electrophoresis [61], visualized, and photographed in AlphaImager® Mini (Alpha Innothech).

Genome sequencing, assembly, and annotation

ErelGV genomic DNA was sequenced with the 454 Genome Sequencer (GS) FLX™ Standard (Roche) at the Centro de Genômica de Alto Desempenho do Distrito Federal (Brasília, Brazil). The genome was assembled de novo using Geneious 6.0 [62] and confirmed using restriction enzyme digestion profile. The annotation was performed using Geneious 6.0 to identify the open reading frames (ORFs) that started with a methionine codon (ATG) encoding at least 50 amino acids and blastp [63] to identify homologues.

Phylogeny, genome, and gene comparisons

For Baculoviridae phylogenetic analysis, a MAFFT alignment [64] was carried out with concatenated amino acid sequences predicted for 37 baculovirus core genes. A maximum likelihood tree was inferred using PhyML with 100 repetitions of a non parametric bootstrap [65], implemented in Geneious, with LG + I + G + F model selected by Prottest 2.4 [66]. Moreover, a genomic comparison was performed using the protein dataset of all the complete Betabaculovirus genomes available in Genbank. The dataset was compared using CGView Comparison Tool [36] and the results were plotted using CIRCOS [37]. We also compared ChocGV and PiraGV genomes with ErelGV genome using Mauve alignment [34]. The horizontal gene transfer (HGTs) events were investigated comparing the maximum likelihood phylogenetic tree inferred using the RAxML method [67] and a MAFT alignment of homologues for he65-like and p43-like, and Densovirus-related genes with 100 repetitions of a non parametric bootstrap for branch support.

Author’s contributions

Conceived and designed the experiments: DMPAA, FLM, BMR, MLS; Performed the experiments: DMPAA, FLM, MSA, WS; Analyzed the data: DMPAA, FLM, BMR; Contributed reagents/materials/analysis tools: BMR, DMPAA, FLM, MSA, SNB, MLS; Wrote the paper: DMPAA, FLM, BMR, MLS. All authors read and approved the final manuscript.


  1. El-Sharkawy MA: Cassava biology and physiology. Plant Mol Biol. 2004, 56 (4): 21-

    Article  Google Scholar 

  2. Pietrowski V, Ringenberger R, Rheinheimer AR, Bellon PP, Gazola D, Miranda AM: Insetos-praga da cultura da mandioca na região Centro-Sul do Brasil. 2010, Marechal Cândido Rondon - Parana: EMBRAPA, 1:

    Google Scholar 

  3. Bellotti AC, Arias B, Guzman OL: Biological control of the cassava hornworm erinnyis ello (Lepidoptera: Sphingidae). The Florida Entomologist. 1992, 75 (4): 10-

    Article  Google Scholar 

  4. Fazolin M, Estrela JLV, Filho MDC, Santiago ACC, Frota FS: Manejo Integrado do Mandarová-da-Mandioca Erinnyis ello (L.) (Lepidoptera: Sphingidae): Conceitos e Experiências na Região do Vale do Rio Juruá, Acre. Edited by: Embrapa. 2007, Rio Branco, Brasil: Embrapa

    Google Scholar 

  5. Schmitt AT: Eficiência da aplicação de Baculovirus erinnyis no controle do mandarová da mandioca. 1985, Florianópolis, Santa Catarina - Brasil: EMPASC Comunicado Técnico edn, 88:

    Google Scholar 

  6. Schmitt AT: Principais insetos pragas da mandioca e seu controle. Cultura de tuberosas amiláceas Latino Americanas. Edited by: Cereda MP. 2002, São Paulo: Fundação Cargill, 2: 7-

    Google Scholar 

  7. Rohrmann GF: Baculovirus Molecular Biology, Third Edition [Internet] edn. 2011, Bethesda (MD): National Center for Biotechnology Information (US), Available from:

    Google Scholar 

  8. Jehle JA, Blissard GW, Bonning BC, Cory JS, Herniou EA, Rohrmann GF, Theilmann DA, Thiem SM, Vlak JM: On the classification and nomenclature of baculoviruses: a proposal for revision. Arch Virol. 2006, 151 (7): 1257-1266. 10.1007/s00705-006-0763-6.

    Article  CAS  PubMed  Google Scholar 

  9. Herniou EA, Arif BM, Becnel JJ, Blissard GW, Bonning B, Harrison R, Jehle JA, Theilmann DA, Vlak JM: Baculoviridae. Virus taxonomy Ninth Report of the International Committee on Taxonomy of Viruses. Edited by: King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ. 2012, San Diego: Elsevier-Academic Press, 163-173.

    Google Scholar 

  10. Moscardi F: Assessment of the application of baculoviruses for control of Lepidoptera. Annu Rev Entomol. 1999, 44: 257-289. 10.1146/annurev.ento.44.1.257.

    Article  CAS  PubMed  Google Scholar 

  11. Ji X, Sutton G, Evans G, Axford D, Owen R, Stuart DI: How baculovirus polyhedra fit square pegs into round holes to robustly package viruses. EMBO J. 2010, 29 (2): 505-514. 10.1038/emboj.2009.352.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  12. Slack J, Arif BM: The baculoviruses occlusion-derived virus: virion structure and function. Adv Virus Res. 2007, 69: 99-165.

    Article  CAS  PubMed  Google Scholar 

  13. Wang R, Deng F, Hou D, Zhao Y, Guo L, Wang H, Hu Z: Proteomics of the Autographa californica nucleopolyhedrovirus budded virions. J Virol. 2010, 84 (14): 7233-7242. 10.1128/JVI.00040-10.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  14. Washburn JO, Chan EY, Volkman LE, Aumiller JJ, Jarvis DL: Early synthesis of budded virus envelope fusion protein GP64 enhances Autographa californica multicapsid nucleopolyhedrovirus virulence in orally infected Heliothis virescens. J Virol. 2003, 77 (1): 280-290. 10.1128/JVI.77.1.280-290.2003.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  15. Herniou EA, Olszewski JA, Cory JS, O'Reilly DR: The genome sequence and evolution of baculoviruses. Annu Rev Entomol. 2003, 48: 211-234. 10.1146/annurev.ento.48.091801.112756.

    Article  CAS  PubMed  Google Scholar 

  16. Herniou EA, Luque T, Chen X, Vlak JM, Winstanley D, Cory JS, O'Reilly DR: Use of whole genome sequence data to infer baculovirus phylogeny. J Virol. 2001, 75 (17): 8117-8126. 10.1128/JVI.75.17.8117-8126.2001.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  17. Hoover K, Grove M, Gardner M, Hughes DP, McNeil J, Slavicek J: A gene for an extended phenotype. Science. 2011, 333 (6048): 1401-10.1126/science.1209199.

    Article  CAS  PubMed  Google Scholar 

  18. Ackermann H, Smirnoff W: A morphological investigation of 23 baculoviruses. J Invertebr Pathol. 1983, 41: 12-

    Article  Google Scholar 

  19. Finnerty CM, Li G, Granados RR: Characterization of a granulovirus from the cassava hornworm (Erinnyis ello: Sphingidae). J Invertebr Pathol. 2000, 75 (4): 273-278. 10.1006/jipa.2000.4929.

    Article  CAS  PubMed  Google Scholar 

  20. 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-193. 10.1016/j.virol.2005.10.032.

    Article  CAS  PubMed  Google Scholar 

  21. Escasa SR, Lauzon HA, Mathur AC, Krell PJ, Arif BM: Sequence analysis of the Choristoneura occidentalis granulovirus genome. J Gen Virol. 2006, 87 (Pt 7): 1917-1933.

    Article  CAS  PubMed  Google Scholar 

  22. Garavaglia MJ, Miele SA, Iserte JA, Belaich MN, Ghiringhelli PD: The ac53, ac78, ac101, and ac103 genes are newly discovered core genes in the family Baculoviridae. J Virol. 2012, 86 (22): 12069-12079. 10.1128/JVI.01873-12.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  23. Wormleaton S, Kuzio J, Winstanley D: The complete sequence of the Adoxophyes orana granulovirus genome. Virology. 2003, 311 (2): 350-365. 10.1016/S0042-6822(03)00149-1.

    Article  CAS  PubMed  Google Scholar 

  24. Zhang X, Liang Z, Yin X, Wang J, Shao X: Complete genome sequence of Agrotis segetum granulovirus Shanghai strain. Arch Virol. 2014,,

    Google Scholar 

  25. Liang Z, Zhang X, Yin X, Cao S, Xu F: Genomic sequencing and analysis of Clostera anachoreta granulovirus. Arch Virol. 2011, 156 (7): 1185-1198. 10.1007/s00705-011-0977-0.

    Article  CAS  PubMed  Google Scholar 

  26. Lange M, Jehle JA: The genome of the Cryptophlebia leucotreta granulovirus. Virology. 2003, 317 (2): 220-236. 10.1016/S0042-6822(03)00515-4.

    Article  CAS  PubMed  Google Scholar 

  27. Luque T, Finch R, Crook N, O'Reilly DR, Winstanley D: The complete sequence of the Cydia pomonella granulovirus genome. J Gen Virol. 2001, 82 (Pt 10): 2531-2547.

    Article  CAS  PubMed  Google Scholar 

  28. Ferrelli ML, Salvador R, Biedma ME, Berretta MF, Haase S, Sciocco-Cap A, Ghiringhelli PD, Romanowski V: Genome of Epinotia aporema granulovirus (EpapGV), a polyorganotropic fast killing betabaculovirus with a novel thymidylate kinase gene. BMC Genomics. 2012, 13: 14-10.1186/1471-2164-13-14.

    Article  Google Scholar 

  29. Harrison RL, Popham HJ: Genomic sequence analysis of a granulovirus isolated from the Old World bollworm, Helicoverpa armigera. Virus Genes. 2008, 36 (3): 565-581. 10.1007/s11262-008-0218-0.

    Article  CAS  PubMed  Google Scholar 

  30. Zhang BQ, Cheng RL, Wang XF, Zhang CX: The Genome of Pieris rapae Granulovirus. J Virol. 2012, 86 (17): 9544-10.1128/JVI.01431-12.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  31. Hashimoto Y, Hayakawa T, Ueno Y, Fujita T, Sano Y, Matsumoto T: Sequence analysis of the Plutella xylostella granulovirus genome. Virology. 2000, 275 (2): 358-372. 10.1006/viro.2000.0530.

    Article  CAS  PubMed  Google Scholar 

  32. Wang Y, Choi JY, Roh JY, Liu Q, Tao XY, Park JB, Kim JS, Je YH: Genomic sequence analysis of granulovirus isolated from the tobacco cutworm, Spodoptera litura. PLoS One. 2011, 6 (11): e28163-10.1371/journal.pone.0028163.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  33. Hayakawa T, Ko R, Okano K, Seong SI, Goto C, Maeda S: Sequence analysis of the Xestia c-nigrum granulovirus genome. Virology. 1999, 262 (2): 277-297. 10.1006/viro.1999.9894.

    Article  CAS  PubMed  Google Scholar 

  34. Darling AC, Mau B, Blattner FR, Perna NT: Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004, 14 (7): 1394-1403. 10.1101/gr.2289704.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. Lange M, Wang H, Zhihong H, Jehle JA: Towards a molecular identification and classification system of lepidopteran-specific baculoviruses. Virology. 2004, 325 (1): 36-47. 10.1016/j.virol.2004.04.023.

    Article  CAS  PubMed  Google Scholar 

  36. Grant JR, Arantes AS, Stothard P: Comparing thousands of circular genomes using the CGView Comparison Tool. BMC Genomics. 2012, 13: 8-10.1186/1471-2164-13-8.

    Article  Google Scholar 

  37. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics. Genome Res. 2009, 19 (9): 1639-1645. 10.1101/gr.092759.109.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Pearson MN, Rohrmann GF: Transfer, incorporation, and substitution of envelope fusion proteins among members of the Baculoviridae, Orthomyxoviridae, and Metaviridae (insect retrovirus) families. J Virol. 2002, 76 (11): 5301-5304. 10.1128/JVI.76.11.5301-5304.2002.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  39. Oomens AG, Blissard GW: Requirement for GP64 to drive efficient budding of Autographa californica multicapsid nucleopolyhedrovirus. Virology. 1999, 254 (2): 297-314. 10.1006/viro.1998.9523.

    Article  CAS  PubMed  Google Scholar 

  40. Wang M, Yin F, Shen S, Tan Y, Deng F, Vlak JM, Hu Z, Wang H: Partial functional rescue of Helicoverpa armigera single nucleocapsid nucleopolyhedrovirus infectivity by replacement of F protein with GP64 from Autographa californica multicapsid nucleopolyhedrovirus. J Virol. 2010, 84 (21): 11505-11514. 10.1128/JVI.00862-10.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  41. Wang M, Tan Y, Yin F, Deng F, Vlak JM, Hu Z, Wang H: The F protein of Helicoverpa armigera single nucleopolyhedrovirus can be substituted functionally with its homologue from Spodoptera exigua multiple nucleopolyhedrovirus. J Gen Virol. 2008, 89 (Pt 3): 791-798.

    Article  CAS  PubMed  Google Scholar 

  42. Lung O, Westenberg M, Vlak JM, Zuidema D, Blissard GW: Pseudotyping Autographa californica multicapsid nucleopolyhedrovirus (AcMNPV): F proteins from group II NPVs are functionally analogous to AcMNPV GP64. J Virol. 2002, 76 (11): 5729-5736. 10.1128/JVI.76.11.5729-5736.2002.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  43. Harrison RL, Lynn DE: Genomic sequence analysis of a nucleopolyhedrovirus isolated from the diamondback moth, Plutella xylostella. Virus Genes. 2007, 35 (3): 857-873. 10.1007/s11262-007-0136-6.

    Article  CAS  PubMed  Google Scholar 

  44. Wennmann JT, Jehle JA: Detection and quantitation of Agrotis baculoviruses in mixed infections. J Virol Methods. 2014, 197: 39-46.

    Article  CAS  PubMed  Google Scholar 

  45. Ko R, Okano K, Maeda S: Structural and functional analysis of the Xestia c-nigrum granulovirus matrix metalloproteinase. J Virol. 2000, 74 (23): 11240-11246. 10.1128/JVI.74.23.11240-11246.2000.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  46. Means JC, Passarelli AL: Viral fibroblast growth factor, matrix metalloproteases, and caspases are associated with enhancing systemic infection by baculoviruses. Proc Natl Acad Sci U S A. 2010, 107 (21): 9825-9830. 10.1073/pnas.0913582107.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  47. Lepore LS, Roelvink PR, Granados RR: Enhancin, the granulosis virus protein that facilitates nucleopolyhedrovirus (NPV) infections, is a metalloprotease. J Invertebr Pathol. 1996, 68 (2): 131-140. 10.1006/jipa.1996.0070.

    Article  CAS  PubMed  Google Scholar 

  48. D'Amico V, Slavicek J, Podgwaite JD, Webb R, Fuester R, Peiffer RA: Deletion of v-chiA from a baculovirus reduces horizontal transmission in the field. Appl Environ Microbiol. 2013, 79 (13): 4056-4064. 10.1128/AEM.00152-13.

    Article  PubMed Central  PubMed  Google Scholar 

  49. Lima AA, Aragao CW, de Castro ME, Oliveira JV, Sosa Gomez DR, Ribeiro BM: A Recombinant MNPV Harboring and Genes from Defective NPV Induce Host Liquefaction and Increased Insecticidal Activity. PLoS One. 2013, 8 (9): e74592-10.1371/journal.pone.0074592.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  50. Cui Q, Shin WS, Luo Y, Tian J, Cui H, Yin D: Thymidylate kinase: an old topic brings new perspectives. Curr Med Chem. 2013, 20 (10): 1286-1305. 10.2174/0929867311320100006.

    Article  CAS  PubMed  Google Scholar 

  51. Penades JR, Donderis J, Garcia-Caballer M, Tormo-Mas MA, Marina A: dUTPases, the unexplored family of signalling molecules. Curr Opin Microbiol. 2013, 16 (2): 163-170. 10.1016/j.mib.2013.02.005.

    Article  CAS  PubMed  Google Scholar 

  52. Priet S, Sire J, Querat G: Uracils as a cellular weapon against viruses and mechanisms of viral escape. Curr HIV Res. 2006, 4 (1): 31-42. 10.2174/157016206775197673.

    Article  CAS  PubMed  Google Scholar 

  53. Castillo-Acosta VM, Aguilar-Pereyra F, Bart JM, Navarro M, Ruiz-Perez LM, Vidal AE, Gonzalez-Pacanowska D: Increased uracil insertion in DNA is cytotoxic and increases the frequency of mutation, double strand break formation and VSG switching in Trypanosoma brucei. DNA Repair (Amst). 2012, 11 (12): 986-995. 10.1016/j.dnarep.2012.09.007.

    Article  CAS  Google Scholar 

  54. Guillet M, Van Der Kemp PA, Boiteux S: dUTPase activity is critical to maintain genetic stability in Saccharomyces cerevisiae. Nucleic Acids Res. 2006, 34 (7): 2056-2066. 10.1093/nar/gkl139.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  55. Herniou EA, Jehle JA: Baculovirus phylogeny and evolution. Curr Drug Targets. 2007, 8 (10): 1043-1050. 10.2174/138945007782151306.

    Article  CAS  PubMed  Google Scholar 

  56. Ho CK, Shuman S: Bacteriophage T4 RNA ligase 2 (gp24.1) exemplifies a family of RNA ligases found in all phylogenetic domains. Proc Natl Acad Sci U S A. 2002, 99 (20): 12709-12714. 10.1073/pnas.192184699.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  57. Yu M, Carstens EB: Characterization of an Autographa californica multiple nucleopolyhedrovirus mutant lacking the ac39(p43) gene. Virus Res. 2011, 155 (1): 300-306. 10.1016/j.virusres.2010.10.025.

    Article  CAS  PubMed  Google Scholar 

  58. Zhu SY, Yi JP, Shen WD, Wang LQ, He HG, Wang Y, Li B, Wang WB: Genomic sequence, organization and characteristics of a new nucleopolyhedrovirus isolated from Clanis bilineata larva. BMC Genomics. 2009, 10: 91-10.1186/1471-2164-10-91.

    Article  PubMed Central  PubMed  Google Scholar 

  59. Lauzon HA, Lucarotti CJ, Krell PJ, Feng Q, Retnakaran A, Arif BM: Sequence and organization of the Neodiprion lecontei nucleopolyhedrovirus genome. J Virol. 2004, 78 (13): 7023-7035. 10.1128/JVI.78.13.7023-7035.2004.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  60. O'Reilly D, Miller LK, Luckow VA: Baculovirus Expression Vectors: a laboratory manual. 1992, New York: Freeman and Company

    Google Scholar 

  61. Sambrook J, Russel DW: Molecular Cloning: a laboratory manual, 3rd ed. edn. 2001, New York: Cold Spring Harbor

    Google Scholar 

  62. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, Thierer T, Ashton B, Meintjes P, Drummond A: Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012, 28 (12): 1647-1649. 10.1093/bioinformatics/bts199.

    Article  PubMed Central  PubMed  Google Scholar 

  63. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  64. Katoh K, Misawa K, Kuma K, Miyata T: MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30 (14): 3059-3066. 10.1093/nar/gkf436.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  65. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59 (3): 307-321. 10.1093/sysbio/syq010.

    Article  CAS  PubMed  Google Scholar 

  66. Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21 (9): 2104-2105. 10.1093/bioinformatics/bti263.

    Article  CAS  PubMed  Google Scholar 

  67. Stamatakis A, Hoover P, Rougemont J: A rapid bootstrap algorithm for the RAxML Web servers. Syst Biol. 2008, 57 (5): 758-771. 10.1080/10635150802429642.

    Article  PubMed  Google Scholar 

Download references


We thank Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), Fundação de Apoio à Pesquisa do Distrito Federal (FAPDF), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), and EMBRAPA Recursos Genéticos e Biotecnologia for the financial support; José Osmar Lorenzi for the picture of E. ello ello infected by ErelGV; Ingrid Gracielle Martins da Silva for kindly helping with the sample for microscopy; and Jeffrey J. Hodgson for kindly reviewing the English writing style.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Bergmann Morais Ribeiro.

Additional information

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material


Additional file 1: Table S1: Characteristics of the Erinnyis ello granulovirus (ErelGV) genome: analysis and homology search. Predicted ORFs are compared with homologous genes in three related genomes. (DOCX 67 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ardisson-Araújo, D.M.P., de Melo, F.L., Andrade, M.d.S. et al. Genome sequence of Erinnyis ello granulovirus (ErelGV), a natural cassava hornworm pesticide and the first sequenced sphingid-infecting betabaculovirus. BMC Genomics 15, 856 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: