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

Abstract

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

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

A P

C.W. Arago-Silva, M. S. Andrade, D. M. P. Ardisson-Arajo, J. E. A. Fernandes,F. S. Morgado, S. N. Bo, R. H. P. Moraes, J. L. C. , F. L. Melo & B. M. Ribeiro

Lonomia obliqua (Lepidoptera: Saturniidae) is a species of medical importance due to the severityof reactions caused by accidental contact with the caterpillar bristles. Several natural pathogens L. obliqua, and among them the baculovirus Lonomia obliqua multiple nucleopolyhedrovirus (LoobMNPV). The complete genome of LoobMNPV was sequenced and shown LoobMNPV genome showed that it belongs to Alphabaculovirus group I (lepidopteran-infective NPV). Danaus plexippus, suggesting an independent acquisition through horizontal gene transfer. Homologs of cathepsin and chitinase genes, which are involved in host integument liquefaction and viral spread, were not found in this genome. As L. obliqua presents a gregarious behavior during the larvae stage the impact of this deletion might be neglectable.

Although the vast majority of caterpillars species do not represent a threat to human health, members of 12 lepidopteran families may cause serious human injuries1. Among these, the larvae Lonomia obliqua (Walker, 1855) (Lepidoptera: Saturniidae) is included, causing several accidents in the Southern region of Brazil2. This insect has been detected feeding on plants from the families Anacardiaceae. and Meliaceae and in city orchards2. When in contact with the human skin, a toxin is released by the larvae bristles, resulting in reactions varying from local irritation to severe life threatening conditions such as coagulopathy, acute renal failure and hemorrhagic disorders3. The gregarious nature of the larvae contributes to a higher venom exposure, which aggravates the severity of the symptoms. There are indications that the occurrence of accidents with L. obliqua is increasing2, possibly due to several factors such as deforestation, the introduction of exotic plants and reduction of natural enemies4.

The need to control this insect population has led to the identication of several pathogens and predators, including a nematode (Hexamermis sp.)5 and a baculovirus: Lonomia obliqua multiple nucleopolyhedrovirus (LoobMNPV), which was isolated from L. obliqua larvae6. The infected larvae exhibited all the usual symptoms of a baculovirus infection, although lacking liquefaction and melanization post mortem phenotypes6. A previous study determined the sequence of a conserved baculovirus gene (polyhedrin) gene that showed close similarity with members of group I alphabaculovirus6.

The rst baculovirus to be completely sequenced was the Autographa californica multiple nucleopolyhedro-virus (AcMNPV) in 19947. Since then, the number of complete genomes has been growing rapidly, providing a wealth of data that contributes to the understanding of both biology and evolution of baculovirus. However, there are only three reports of complete genomes from baculoviruses isolated from species of the family Saturnidae: Hemileuca sp. nucleopolyhedrovirus8, Philosamia cynthia nucleopolyhedrovirus9, and two isolates of Antheraea

Departamento de Biologia Celular, Instituto de Cincias Biolgicas, Universidade de Braslia, Braslia, DF, Brazil.

Laboratrio de Biologia Molecular e Virologia, Centro de Cincias Biolgicas e da Sade (CCBS), Universidade Presbiteriana Mackenzie, So Paulo, SP, Brazil. Departamento de Entomologia, Instituto Butantan, Av. Vital [email protected])

SCIENTIFIC REPORTS

1

www.nature.com/scientificreports/

Figure 1. Genome map and heat map of the LoobMNPV genome. (a) Genome map showing all the 134 ORFs. The unique ORFs are represented in black. The outer track contains forward orientation ORFs, and the inner track contains reverse orientation ORFs. Hrs are shown on the line below the genome. (b) Heat map identity of the genomes of the species AcMNPV, ThorNPV, MaviMNPV, DekiNPV, and AnpeNPV (from the outside to the inside) compared to ortholog ORFs from LoobMNPV. The darker the blue, the higher is the correlated ORF identity.

pernyi nucleopolyhedrovirus10. Therefore, the complete genome of LoobMNPV provides additional data that may help the understanding of baculovirus genome evolution and adaptations to their hosts.

Results and Discussion

Genome Features. The assembled genome of LoobMNPV (GenBank accession number: KP763670) is 120,023-bp long with a mean coverage of 20.5 times. We have found 134 Open Reading Frames (ORFs) coding for putative proteins with at least 50 amino acid residues (Fig.1a and Table S1). Out of these, 120 ORFs were found in other baculoviruses and among those, three ORFs (loob078, loob100 (he65) and loob113) showed higher identity with orthologs from betabaculoviruses (Table S1), conrming that horizontal transference between alphabaculovirus and betabaculovirus is indeed a common event, as previously documented11. Moreover, some unique ORFs were located within the hr3 (loob038) and hr4 (loob060 and loob061). All the 37 baculovirus core genes were present, however some genes present in almost all Alpha- and Betabaculovirus12 genomes available were not present: Ac106/107, Ac108 and gp37.

LoobMNPV homologous regions (hrs). Seven homologous regions (hrs) were found in the LoobMNPV genome, in agreement to other baculoviruses13. These DNA palindromic repetitive elements are interspersed throughout most baculoviral genomes and are related to DNA replication14, gene transcription15, and possibly homologous recombination16. The hr1 has 684 bp, hr2 has 541 bp, hr3 has 1,037 bp, hr4 has 1,426 bp, hr5 has 679bp, hr6 has 361bp and hr7 has 836bp, and their size and position were conrmed by PCR (data not shown). As shown in Fig.2, all hrs have various copies of a common imperfect palindromic repeat of 38bp (hr1a, hr2a, hr3a, hr3b, hr4a, hr5a, hr5b, hr6a, hr6b and hr7a).

LoobMNPV Phylogeny. We have performed a phylogenetic analysis (Fig.3) based on the concatenated amino acid sequence alignment of 37 baculovirus core genes of 72 baculovirus species available in GenBank (Table S2). LoobMNPV clustered with group I alphabaculovirus. Its closest relative was found to be DekiNPV and both are basally related to AcMNPV-like viruses. This result disagrees with previous analysis based on the

SCIENTIFIC REPORTS

2

www.nature.com/scientificreports/

Figure 2. LoobMNPV hr palindromes. Alignment showing the position and the sequence of 38 hr-like imperfect palindromes found in the LoobMNPV genome, numbered sequentially, where letters designate palindromes within the same hr. The black shaded areas indicate residues that have strict identity, whereas grey shaded areas indicate conservation within the majority of the group (no strict identity). Dots represent gaps to achive a better alignment.

Figure 3. Phylogeny of baculoviruses. Phylogenetic inference of 37 conserved proteins (core genes) present in 72 baculovirus genomes from dierent host species. The genera Gammabaculovirus, Betabaculovirus, and group II Alphaphabaculovirus are collapsed. CuniNPV was used to root the tree. LoobMNPV belongs to the genus Alphabaculovirus, on group I, and clustered with DekiNPV. Both species are a sister clade of AcMNPV-related species.

polh gene, in which LoobMNPV clustered together with both species Amsacta albistriga nucleopolyhedrovirus (AaNPV) and Thysanoplusia orichalcea nucleopolyhedrovirus (ThorNPV)6. However, the phylogenetic signal of highly conserved genes, such as polh, is usually inaccurate17, while phylogenetic inferences based on a large set of genes is more accurate and robust18.

LoobMNPV Comparative Genomics. Based on the phylogenetic results we have selected some of the closest relatives of LoobMNPV, including AcMNPV, Antheraea pernyi nucleopolyhedrovirus (AnpeNPV), Maruca vitrata nucleopolyhedrovirus (MaviNPV), Dendrolimus kikuchii nucleopolyhedrovirus (DekiNPV) and ThorNPV for comparative genomics. Figure1b and Table S1 reveal all the LoobMNPV ORFs compared to orthologs from these genomes by the level of similarity (in terms of percentage of identity and E-values< 103). Overall, some

SCIENTIFIC REPORTS

3

www.nature.com/scientificreports/

Figure 4. Genome sinteny of LoobMNPV. (a) LoobMNPV sinteny comparison to the baculoviruses AnpeNPV (salmon), MaviMNPV (light pink), ThorNPV (blue), AcMNPV (green), and DekiNPV (orange). Each line represent an ortholog ORF to LoobMNPV based on the protein identity. LoobMNPV and DekiNPV are collinear, while the other genomes show an inversion when compared to LoobMNPV. (b) Zoom from (a) on the location where all baculoviruses compared overlap, showing hrs located closely to LoobMNPV unique ORFs (white gaps).

ORFs are more conserved (high similarity), such as the genes that correspond to polh (loob001), p74 (loob010), p49 (loob014), pif-1 (loob041), vlf-1 (loob067), lef-9 (loob080), lef-8 (loob093) and pif-2 (loob123); while others are more variable, presenting a lower level of similarity, but still present in almost all baculovirus genomes, such as ORF1629 (loob002), ie-2 (loob005), vp80 (loob043), desmoplakin (loob075), f protein (loob122), arif-1 (loob124) and bv/odv-e26 (loob128).

Moreover, we have investigated the genome synteny among LoobMNPV and those selected alphabaculovirus genomes (Fig.4). A circular ideogram layout shown in Fig.4 displays lines connecting related genes, where it is possible to observe the relative position compared to LoobMNPV genome. Interestingly, we have observed that LoobMNPV and DekiNPV show high collinearity, and both present genome inversions and rearrangement in contrast to the genomes of AcMNPV, AnpeNPV, MaviNPV, and ThorNPV. The inverted regions were adjacent or anked by the hr1 and hr2 (Fig. S1). These ndings support the phylogenetic results presented in Fig.3 and suggests that the inversions were originated in a common ancestor of both LoobMNPV and DekiNPV. These inversions are an autapomorphy when compared to the AcMNPV-like viruses. Previous studies have shown that rearrangements in the baculovirus genomes reect evolutionary history, with more closely related viruses presenting higher genome collinearity19.

Furthermore, when analyzing the region where all genomes overlap with LoobMNPV (Fig.4b), we have observed that unique ORFs in LoobMNPV (gaps with no correspondence to the other genomes) are interestingly

SCIENTIFIC REPORTS

4

www.nature.com/scientificreports/

Figure 5. Phylogeny of GTAs and TTF2 genes. Unrooted maximum likelihood phylogeny of the data set containing genes that correspond to TTF2 from Insecta (red), TTF2 from other eukaryote (orange), GTA from group I Alphabaculovirus (green), and loob035 (blue).

found near hrs. According to previous reports, gene rearrangements and acquisitions are of common occurrence on proximities of hrs20, conrming the possibility of gene transfers to viruses by homologous recombination21, which could be also facilitated by factors such as the prevalence of various pathogens infecting the same host22, as well as concomitant infections in eld populations23.

LoobMNPV genome showed 12 genes that do not have any match among baculoviruses. These ORFs are loob004, loob006, loob012, loob035, loob038, loob055, loob059, loob060, loob061, loob071, loob084, and loob097. There are three possible mechanisms for gene acquisition: extensive sequence divergence, which could push homolog genes below the similarity threshold; gene recombination, which produces novel protein products; and horizontal gene transfer (HGT). This third possibility is expected to be detectable by gene similarity from phylogenetically distinct species24. For each of these LoobMNPV unique ORFs, we have searched for baculovirus promoter motifs within 200 bp upstream of the start codon. The late promoter motif TAAG, that appears to be necessary for late transcription by the viral RNA polymerase25 was found in loob004, loob012, loob030, loob038, loob055, and loob071. However, further experiments are necessary to conrm whether these ORFs encode bona de proteins. Moreover, these novel ORFs were searched for known domains and eight of theses ORFs (loob004, loob006, loob012, loob038, loob055, loob061, loob071, and loob097) did not match any predicted domain. The remaining will be discussed below.

Insect immune system-associated domain. The loob060 has an immunoglobulin-like domain, which has been found in some insect proteins, such as the hemolin, an hemolymph component that plays a role in bacterial surface binding, forming a protein complex that initiates the immune response26. Hemolin has also been found in the transcripts of L. obliqua bristles2. Several immunomodulators encoded by viruses have been described27 and may be involved in regulating the immune system and protecting virus-infected cells from the attack of other cells from the immune system28,29. For viruses, the expression of these proteins may indicate benecial susceptibility in multiple pathogen infections, by protecting the host against opportunistic pathogens, reducing competition and beneting viral propagation29. For instance, other saturniid-related alphabaculovirus species HespNPV expresses a functional insect-related serine protease inhibitor (serpin) in its genome that is likely related to host immunity modulation and virulence8,28.

Transcription factor-related domain (). One unique ORF demonstrated high correspondence to the eukaryotic transcription terminator factor type 2 (TTF2) from the buttery Danaus plexippus (GenBank: EHJ68439), with 44% pairwise identity and E-value equals to 3e1011. However, when ltering this result, in order to focus only on the family Baculoviridae, the referred gene presented higher similarity to the Global Transactivator (GTA) gene from the AnpeNPV (YP_611073), with 66% of identity and an E-value of 1e106.

Transcription Factors (TFs) in general are fundamental in a broad array of any cellular processes due to its ability of causing changes in downstream gene expression patterns30. GTA genes are observed in members from the group I alphabaculovirus. According to a previous study31, baculovirus GTAs play an important role in transcriptional activation of viral genes and were probably originated by HGT from the host to the common ancestor of the clade that includes AcMNPV, Bombyx mori nucleopolyhedrovirus (BmNPV), Orgyia pseudotsugata multicapsid nucleopolyhedrovirus (OpMNPV), and Epiphyas postvittana nucleopolyhedrovirus (EppoNPV). Katsuma, et al.32 found that a GTA homolog from BmNPV acts as a viral virulence factor in insect larvae, and may be required for activation of host and/or viral genes, increasing the speed of host killing. Based on domain analysis, we have found that the loob035, TTF2, and GTA genes are members of the SNF2 family. This family of genes encodes proteins with sequence motifs similar to those found in many DNA and RNA helicase protein families, and also proteins from a variety of species with roles in cellular processes such as transcriptional regulation, DNA recombination, chromatin unwinding and various other types of DNA repair33. In this regard, the possible acquisition of loob035 might be involved in the inhibition of the host transcriptional machinery in order to benet viral expression.

To analyze whether loob035 has been independently acquired from the host insect through HGT, or is a divergent baculovirus GTA gene, a phylogenetic analysis has been performeddone, based on an amino acid alignment containing both TTF2 and GTA sequences. As shown in Fig.5, loob035 presented a long branch size, which

SCIENTIFIC REPORTS

5

www.nature.com/scientificreports/

Figure 6. GTA genomic context in several alphabaculoviruses. (a) The GTA gene is present in almost all group I Alphabaculovirus and is shown here aligned to the phylogeny previously showed on Fig.3. In AcMNPV, AgMNPV, AnpeNPV, BmNPV, BomaNPV, CfMNPV, ChocNPV, ChmuNPV, ChroNPV, DekiNPV, EppoNPV, HycuMNPV, LoobMNPV, MaviMNPV, OpMNPV, PhcyNPV, PlxyMNPV, RoMNPV e ThorNPV, the GTA gene is always in between p-47 (followed by lef-12) and odv-e56. LoobMNPV, DekiNPV and MaviMNPV lack the GTA gene. Dotted lines represent absence of ORFs. (b) The position of loob035 in the genome of LoobMNPV is located between homologs of both ac110-like (loob034) and ac111-like (loob036). The position of two loob035 homologs (Table S1) from DekiNPV_Orf138 and ThorNPV_Orf117 are also represented in dierent genome contexts.

indicates great divergence from the other sequences, possibly due to positive selection34. Actually, TFs are over-represented among genes predicted to be positively selected in previous genome-wide selection studies35. To conrm whether loob035 is in fact a new acquisition from insect host or a divergent baculovirus GTA gene, we compared the likelihood of a tree constraining loob035 to the GTA alphabaculovirus group (lnL = 14940.61) with the likelihood of a tree constraining loob035 into TTF2 group (lnL = 1455.06). We argue that the likelihood dierences in combination with high non-parametric bootstrap values and high posterior probabilities constitute considerable evidence that loob035 clusters with the group of TTF2 genes, corroborating to the possibility of a novel HGT.

We further conrmed this independent acquisition hypothesis by analyzing the genomic context of baculovirus GTA genes. As shown in Fig.6a, GTA genes are encountered in all group I alphabaculoviruses:Ac-MNPV, Anticarsia gemmatalis nucleopolyhedrovirus (AgMNPV), AnpeNPV, BmNPV, Bombyx mandarina nucleopolyhedrovirus (BomaNPV), Choristoneura fumiferana multicapsid nucleopolyhedrovirus (CfMNPV), Choristoneura occidentalis nucleopolyhedrovirus (ChocNPV), Choristoneura murinana nucleopolyhedrovirus

SCIENTIFIC REPORTS

6

www.nature.com/scientificreports/

(ChmuNPV), Choristoneura nucleopolyhedrovirus roaceana (ChroNPV), EppoNPV, Hyphantria cunea nucleopolyhedrovirus (HycuMNPV), Orgyia pseudotsugata multicapsid nucleopolyhedrovirus (OpMNPV), Philosamia cynthia nucleopolyhedrovirus (PhcyNPV), Plutella xylostella multiple nucleopolyhedrovirus (PlxyMNPV), Rachiplusia ou multicapsid nucleopolyhedrovirus (RoMNPV) and ThorNPV; except in MaviMNPV, LoobMNPV and DekiNPV, within a conserved position between lef-12 and odv-e66. However, in LoobMNPV, the gene located in this position is loob102, that corresponds to an AcMNPV-like gene (ac044). On the other hand, loob035 is inserted in a completely dierent genome context, located between both the ac110- and the ac111-like genes (Fig.6b), conrming that loob035 has probably a dierent origin not related to the GTA gene. Besides, loob035 homologs found in DekiNPV (Orf 138) and in ThorNPV (Orf 117), according to Table S1, are also inserted in a dierent context (Fig.6b).

Notably, loob035 diverges greatly from all the other compared sequences (Fig. S2), showing that these sequences present similarity only because they all contain the SNF2 conserved domain. Hughes & Friedman31

found that SNF2 baculovirus gene family has homologs in cellular organisms, and clustered closer to homologs in insects (Drosophila), according to the reconstruction of the evolutionary relationship among genes that were potentially acquired through HGT in comparison to baculovirus phylogeny.

The absence of cathepsin and chitinase genes in LoobMNPV. Interestingly, LoobMNPV does not encode two common baculovirus genes that are responsible for the post mortem host melanization and liquefaction beneting virus dissemination: the enzymes cathepsin (v-cath) and chitinase (chiA) genes36. Besides the involvement of these genes in the horizontal spread of the virus in the eld37, it has been reported that the chiA gene and the pro-form of v-cath interact directly and are dependent on each other for the promotion of host liquefaction, and therefore, they are usually acquired or lost together38, since they are adjacent genes on the genome. Among alpha-baculovirus from group I, only LoobMNPV, AgMNPV and PhcyNPV lack v-cath and chiA genes. The recombinant introduction of v-cath and chiA genes from Choristoneura fumiferana DEF multiple nucleopolyhedrovirus (CfDEFNPV) into AgMNPV genome improves production of occlusion bodies and insecticidal activity during A. gemmatalis larvae infection39.

Several hosts from the family Saturniidae40, as well as from Noctuidae, present gregarious behavior41, which may facilitate virus dispersion to new susceptible individuals and reduce selective pressure for the maintenance of v-cath and chiA genes, as observed for LoobMNPV, AgMNPV, and PhcyNPV. However, HespNPV and AnpeNPV also infect gregarious hosts from Saturniidae family, but have v-cath and chiA genes in their genomes. Therefore, it remains to be determined whether the host behavior it is indeed a selective factor for cath and chiA genes maintenance.

Conclusions. In this work, we described the rst complete genome sequence of a baculovirus isolated from a species of medical interest. LoobMNPV is located on a basal position of group I alphabaculovirus and presents inversions in large proportions when compared to the other related genomes. During evolution, uxes in the genomes content, such as genes acquisitions and losses, pressured by positive selection, could possibly implicate in shis on the evolutionary dynamics, by the occurrence of events of adaptation. Therefore, the elucidation of novel genomes will help the studies on baculovirus evolution, proportioning deeper knowledge and understanding of baculovirus as a whole group.

Materials and Methods

Virus origin. The LoobMNPV used in this work was collected in the year 2000, from orchards in the state of Santa Catarina, Brazil6.

Occlusion bodies (OBs) from dead L. obliqua larvae were puried and dissolved in an alkaline solution for DNA purication, all according to previously published protocol to OReilly et al.42. The quantity and quality of the isolated DNA was determined by electrophoresis on 0.8% agarose gel (data not shown). The DNA was submitted to genome sequencing throughout the pyrosequencing technique performed by the 454 Genome Sequencer (GS) FLX Standard (Roche) at Macrogen Inc (Seoul,

Republic of Korea).

Genome assembly and annotation. Genome de novo assembly and annotation was performed using Geneious version 7.1.6 (available at http://www.geneious.com/)43 considering the following established parameters to select the Open Reading Frames (ORFs): start codons corresponding to the methionine residue (ATG), minimal overlapping of adjacent ORFs, and ORFs as large as 150 bp. No gaps were found between connected contigs, only a 700 bp region of low coverage, that was conrmed by PCR, In order to validate assembly, an in silico digestion was carried out using BamHI, EcoRI and PstI restriction enzymes and comparing to a previously published restriction prole6, that used the same virus isolate as the one used in this work (Table S3). ORFs were annotated using BLASTx and PSI-BLASTp searches against the NCBI non-redundant protein database44. When compared to all baculoviruses genomes available at Genbank, unique ORFs were considered when no significant database hits (E-values > 103) were found, and for those, a more sensitive search was performed using HMMer-search against the PFAM-database45. The presence of protein specic domains in HMMer-search unique ORFs was investigated using InterProScan database46. The regions within 200bp upstream of the putative unique ORFs were screened for the presence of TATA-box and CAGT47, as well as GATA motifs48.

Homologous regions (hrs). The homologous regions (hrs) were annotated using DOTPLOT analysis and Tandem Repeat Finder (http://tandem.bu.edu/trf/trf.html)49. An alignment of the repeat unit of each hr was performed using the MAFFT method50. To conrm the size and position of the hrs among the genome, PCRs were performed, using the isolated LoobMNPV DNA and the following primers: hr1 (F: AGA GTT GGA AAT TTC

SCIENTIFIC REPORTS

7

www.nature.com/scientificreports/

GCG CTC and R: GTT TTT ACT CTG TCC GCG CG); hr2 (F: CCC GCT AAT GAA CCG TGT GA and R: AAC CGT TTA AAT CCT TCG TGT); hr3 (F: GCT GGA GTA AAT TGT TCA ATC GC and R: TTT CCA TAA CGG GGT GCC AA); hr4 (F: TAG GGC ACA ATA GCA GCA GC and R: ACG TGC CAA GTC GAA TCT GA); hr5 (F: CGC ATA ACC TTT AGC GTG ACT and R: CTG AAA CGC GAC AAC AGT CC); hr6 (F: AGA GTT GGA AAT TTC GCG CTC and R: TCA TGT CGG CCA ATG AGG AC) and hr7 (F: AAT GCG CAA AAG AAC GGG TC and R: AAC AAC TAA ACT GCG CCC CA).

Phylogenetic analysis and genome comparisons. A MAFFT alignment was performed using amino acid sequences from the predicted 37 core genes from 72 baculovirus genomes (one isolate from each baculovirus species) available in the GenBank up to the date. The phylogeny was inferred by three dierent phylogenetic algorithms including FastTree51, which uses SH-like method for branch support, along with RAxML52 and PhyML53 that use the bootstrap method for branch support. The phylogenies were inferred by the maximum likelihood (ML) method, along with GTR model of DNA evolution, estimated with JmodelProttest 2.4 soware program54.

Moreover, a genomic comparison was performed using blastp identity results of a protein dataset with four alpha-baculovirus genomes AcMNPV, AnpeNPV, MaviNPV, DekiNPV and ThorNPV. This dataset was compared using CGView Comparison Tool55 and the results were plotted using CIRCOS56. Likewise, to investigate the evolutionary relationship between loob035 and the GTA and TTF2 genes, a MAFFT alignment was performed using 14 GTA genes from alphabaculoviruses and 62 TTF2 genes from insects and other eukaryotes. This phylogenetic tree was done by ML method implemented in PhyML53, with JTT substitution model of amino acid evolution, and 100 repetitions of a non-parametric bootstrap. A multiple alignment was performed on the T-coee platform57 with the predicted loob035 along with the GTA and TTF2 proteins, and was submitted to the ESPript 3 platform58. The two most similar crystal structures to loob035 were used to predict the secondary structures of the alignment. They were extracted from the Protein Data Bank (PDB ID 3mwy59 and PDB ID 1z6360), and were aligned along with the GTA genes from AnpeNPV and EppoMPV and TTF2 genes from Danaus plexippus, Bombyx mori, Chelonia mydas and Pterotopus alecto. To show statistical evidence concerning the best hypothesis for loob035 position in the GTA/TTF2 tree, the Bayes factor estimator61 from MrBayes program62 was used.

References

1. Diaz, J. H. The evolving global epidemiology, syndromic classication, management, and prevention of caterpillar envenoming. Am. J. Trop. Med. Hyg. 72, 347357 (2005).2. Carrijo-Carvalho, L. C. & Chudzinski-Tavassi, A. M. The venom of the Lonomia caterpillar: an overview. Toxicon 49, 741757 (2007).

3. Gamborgi, G. P., Metcalf, E. B. & Barros, E. J. G. Acute renal failure provoked by toxin from caterpillars of the species Lonomia obliqua. Toxicon 47, 6874 (2006).

4. Abella, H. B. et al. Manual de diagnstico e tratamento de acidentes porLonomia.20 (Centro de Informao Toxicolgica-SS/RSFEPPS, 1999).

5. Moraes, R. H. P. Identicao dos Inimigos Naturais de Lonomia obliqua Walker, 1855 (Lepidptera: Saturniidae) e possveis fatos determinantes do aumento de sua populao. Masters Thesis, Escola Superior de Agricultura Luiz de Queiroz (2002).

6. Wol, J. L., Moraes, R. H. P., Kitajima, E., Leal, E. D. S. & Zanotto, P. M. A. Identication and characterization of a baculovirus from Lonomia obliqua (Lepidoptera: Saturniidae). J. Invertebr. Pathol. 79, 137145 (2002).

7. Ayres, M. D., Howard, S. C., Kuzio, J., Lopez-Ferber, M. & Possee, R. D. The complete DNA sequence of Autographa californica nuclear polyhedrosis virus. Virology 202, 586605 (1994).

8. Rohrmann, G. F., Erlandson, M. A. & Theilmann, D. A. The genome of a baculovirus isolated from Hemileuca sp. encodes a serpin ortholog. Virus Genes 47, 35764 (2013).

9. Qian, H. et al. Analysis of the genomic sequence of Philosamia cynthia nucleopolyhedrin virus and comparison with Antheraea pernyi nucleopolyhedrin virus. BMC Genomics 14, 115 (2013).

10. Nie, Z. M. et al. Complete sequence and organization of Antheraea pernyi nucleopolyhedrovirus, a dr-rich baculovirus. BMC Genomics 8, 248 (2007).

11. Ardisson-Arajo, D. M. P. et al. Genome sequence of Erinnyis ello granulovirus (ErelGV), a natural cassava hornworm pesticide and the rst sequenced sphingid-infecting betabaculovirus. BMC Genomics 15, 856 (2014).

12. Garavaglia, M. J., Miele, S. A., Iserte, J. A., Belaich, M. N. & Ghiringhelli, P. D. Ac53, ac78, ac101 and ac103 are newly discovered core genes in the family Baculoviridae. J. Virol. 86, 1206912079 (2012).

13. Oliveira, J. V. et al. Genome of the most widely used viral biopesticide: Anticarsia gemmatalis multiple nucleopolyhedrovirus. J. Gen. Virol. 87, 32333250 (2006).

14. Pearson, M. N. & Rohrmann, G. F. Lymantria dispar nuclear polyhedrosis virus homologous regions : characterization of their ability to function as replication origins. J. Virol. 69, 213221 (1995).

15. Pang, Y. et al. Sequence analysis of the Spodoptera litura multicapsid nucleopolyhedrovirus genome. Virology 287, 391404 (2001).16. Hayakawa, T. et al. Sequence analysis of the Xestia c-nigrum granulovirus genome. Virology 262, 277297 (1999).17. Jehle, J. A. et al. On the classication and nomenclature of baculoviruses: a proposal for revision. Arch. Virol. 151, 12571266 (2006).18. Herniou, E. A., Olszewski, J. A., Cory, J. S. & OReilly, D. R. The genome sequence and evolution of baculoviruses. Annu. Rev. Entomol. 48, 211234 (2003).

19. Hu, Z. H. et al. Distinct gene arrangement in the Buzura suppressaria single-nucleocapsid nucleopolyhedrovirus genome physical map was constructed for the restriction. J. Gen. Virol. 79, 28412851 (1998).

20. Thumbi, D. K., Bliveau, C., Cusson, M., Lapointe, R. & Lucarotti, C. J. Comparative genome sequence analysis of Choristoneura occidentalis Freeman and C. rosaceana Harris (Lepidoptera: Tortricidae) alphabaculoviruses. PLoS One 8, e68968 (2013).

21. De Jong, J. G. et al. Analysis of the Choristoneura fumiferana nucleopolyhedrovirus genome. J. Gen. Virol. 86, 929943 (2005).22. Eveleigh, E. S. et al. Fluctuations in density of an outbreak species drive diversity cascades in food webs. Proc. Natl. Acad. Sci.USA 104, 1697616981 (2007).

23. Kemp, E. M., Woodward, D. T. & Cory, J. S. Detection of single and mixed covert baculovirus infections in eastern spruce budworm, Choristoneura fumiferana populations. J. Invertebr. Pathol. 107, 202205 (2011).

24. McLysaght, A., Baldi, P. F. & Gaut, B. S. Extensive gene gain associated with adaptive evolution of poxviruses. Proc. Natl. Acad. Sci. USA 100, 1565515660 (2003).

25. Garrity, D. B., Chang, M. J. & Blissard, G. W. Late Promoter Selection in the Baculovirus gp64 Envelope Fusion Protein Gene. Virology 231, 167181 (1997).

26. Sun, S. C., Lindstrom, I., Boman, H. G., Faye, I. & Schmidt, O. Hemolin : An Insect-Immune Protein Belonging to the Immunoglobulin Superfamily. Science 250, 17291732 (1990).

SCIENTIFIC REPORTS

8

www.nature.com/scientificreports/

27. Spriggs, M. One step ahead of the game: viral immuno-modulatory molecules. Annu. Reg. Immunol. 14, 101130 (1996).28. Beersma, M. F., Bijlmakers, M. J. & Ploegh, H. L. Human cytomegalovirus down-regulates HLA class I expression by reducing the stability of class I H chains. J. Immunol. 151, 44554464 (1993).

29. Ardisson-Araujo, D. M. P., Rohrmann, G. F., Ribeiro, B. M. & Clem, R. J. Functional characterization of hesp018, a baculovirus-encoded serpin gene. J. Gen. Virol. 96, 11501160 (2015).

30. Weatherbee, S. D. et al. Ultrabithorax function in buttery wings and the evolution of insect wing patterns. Curr. Biol. 9, 109115 (1999).

31. Hughes, A. L. & Friedman, R. Genome-wide survey for genes horizontally transferred from cellular organisms to baculoviruses. Mol. Biol. Evol. 20, 979987 (2003).

32. Katsuma, S., Fujii, T., Kawaoka, S. & Shimada, T. Bombyx mori nucleopolyhedrovirus SNF2 global transactivator homologue (Bm33) enhances viral pathogenicity in B. mori larvae. J. Gen. Virol. 89, 30393046 (2008).

33. Eisen, J. A., Sweder, K. S. & Hanawalt, P. C. Evolution of the SNF2 family of proteins: subfamilies with distinct sequences and functions. Nucleic Acids Res. 23, 27152723 (1995).

34. Suzuki, Y. & Gojobori, T. A method for detecting positive selection at single amino acid sites. Mol. Biol. Evol. 16, 13151328 (1999).35. Bustamante, C. D. et al. Natural selection on protein-coding genes in the human genome. Nature 437, 11531157 (2005).36. Hawtin, R. E. et al. Liquefaction of Autographa californica nucleopolyhedrovirus-infected insects is dependent on the integrity of virus-encoded chitinase and cathepsin genes. Virology 238, 243253 (1997).

37. Cory, J. S. & Myers, J. H. The Ecology and Evolution of Insect Baculoviruses. Annu. Rev. Ecol. 34, 239272 (2003).38. Hodgson, J. J., Arif, B. M. & Krell, P. J. Interaction of Autographa californica multiple nucleopolyhedrovirus cathepsin protease progenitor (proV-CATH) with insect baculovirus chitinase as a mechanism for proV-CATH cellular retention. J. Virol. 85, 39183929 (2011).

39. Lima, A. A. et al. A recombinant Anticarsia gemmatalis MNPV harboring chiA and v-cath genes from Choristoneura fumiferana defective NPV induce host liquefaction and increased insecticidal activity. PLoS One 8, e74592 (2013).

40. Lemaire, C. The Saturniidae of America-Hemileucinae. Vol. 3, 1388 (Goecke & Evers, 2002).41. Silva, F. W. et al. Twos a crowd: phenotypic adjustments and prophylaxis in Anticarsia gemmatalis larvae are triggered by the presence of conspecics. PLoS One 8, e61582 (2013).

42. OReilly, D., Miller, L. K. & Luckrow, V. A. Baculovirus Expression Vectors: A laboratory manual. (Freeman and Company, 1992).43. Kearse, M. et al. Geneious Basic: an integrated and extendable desktop soware platform for the organization and analysis of sequence data. Bioinformatics 28, 16471649 (2012).

44. Geer, L. Y. et al. The NCBI BioSystems database. Nucleic Acids Res. 38, D492D496 (2010).45. Finn, R. D. et al. Pfam: the protein families database. Nucleic Acids Res. 42, D222D230 (2014).46. Zdobnov, E. M. & Apweiler, R. InterProScan-an integration platform for the signature-recognition methods in InterPro. Bioinformatics 17, 847848 (2001).

47. Blissard, G. W. & Rohrmann, G. F. Baculovirus diversity and molecular biology. Annu. Rev. Entomol. 35, 127155 (1990).48. Ko, L. J. & Engel, J. D. DNA-Binding Specicities of the GATA Transcription Factor Family. Mol Cell Biol. 13, 40114022 (1993).49. Benson, G. Tandem repeats nder: a program to analyze DNA sequences. Nucleic Acids Res. 27, 573580 (1999).50. 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. 30, 30593066 (2002).

51. Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree: computing large minimum evolution trees with proles instead of a distance matrix. Mol. Biol. Evol. 26, 16411650 (2009).

52. Stamatakis, A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22, 26882690 (2006).

53. Guindon, S. & Gascuel, O. PhyML: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol. 52, 696704 (2003).

54. Posada, D. jModelTest : phylogenetic model averaging. Mol. Biol. Evol. 25, 12531256 (2008).55. Grant, J. R. & Stothard, P. The CGView Server: a comparative genomics tool for circular genomes. Nucleic Acids Res. 36, W181 W184 (2008).

56. Krzywinski, M. et al. Circos : An information aesthetic for comparative genomics. Genome Res. 19, 16391645 (2009).57. Notredame, C., Higgins, D. G. & Heringa, J. T-Coee: A novel method for fast and accurate multiple sequence alignment. J. Mol. Biol. 302, 205217 (2000).

58. Robert, X. & Gouet, P. Deciphering key features in protein structures with the new ENDscript server. Nucl. Acids Res. 42, W320 W324 (2014).

59. Hauk, G., McKnight, J. N., Nodelman, I. M. & Bowman, G. D. The chromodomains of the Chd1 chromatin remodeler regulate DNA access to the ATPase motor. Mol.Cell 39, 711723 (2010).

60. Durr, H., Korner, C., Muller, M., Hickmann, V. & Hopfner, K. P. X-ray structures of the Sulfolobus solfataricus SWI2/SNF20052 ATPase core and its complex with DNA. Cell 121, 363373 (2005).

61. Baele, G., Lemey, P. & Vansteelandt, S. Make the most of your samples : Bayes factor estimators for high-dimensional models of sequence evolution. BMC Bioinformatics 14, 85 (2013).

62. Ronquist, F. et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 61, 539542 (2012).

We thank the Ministry of Science, Technology and Innovation, Conselho Nacional de Desenvolvimento Cientco e Tecnolgico (CNPq): Grant number 407908/2013-7; Fundao de Apoio Pesquisa do Distrito Federal -FAPDF: Grant number 193.000.583/2009 for nancial support and Coordenao de Aperfeioamento de Pessoal de Nvel Superior (CAPES).

Author Contributions

C.W.A.S., D.M.P.A., F.L.M., J.L.C.W., M.S.A. and B.M.R. wrote the main manuscript text. D.M.P.A., F.L.M., S.N.B., J.E.A.F. and B.M.R. designed and conceived the experiments. C.W.A.S., D.M.P.A., F.L.M., J.E.A.F., F.S.M., M.S.A. and R.H.P.M. performed experiments. All authors reviewed the manuscript.

Additional Information

Supplementary information accompanies this paper at http://www.nature.com/srep

Competing nancial interests: The authors declare no competing nancial interests.

SCIENTIFIC REPORTS

9

www.nature.com/scientificreports/

How to cite this article: Arago-Silva, C. W. et al. The complete genome of a baculovirus isolated from an insect of medical interest: Lonomia obliqua (Lepidoptera: Saturniidae). Sci. Rep. 6, 23127; doi: 10.1038/srep23127 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

SCIENTIFIC REPORTS

10