Genome sequence of Perigonia lusca single nucleopolyhedrovirus: insights into the evolution of a nucleotide metabolism enzyme in the family Baculoviridae

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

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

Genome analysis of Psilogramma increta granulovirus and its intrapopulation diversity

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

Viral dUTPases: Modulators of Innate Immunity

Most free-living organisms encode for a deoxyuridine triphosphate nucleotidohydrolase (dUTPase; EC 3.6.1.23). dUTPases represent a family of metalloenzymes that catalyze the hydrolysis of dUTP to dUMP and pyrophosphate, preventing dUTP from being incorporated into DNA by DNA polymerases, maintaining a low dUTP/dTTP pool ratio and providing a necessary precursor for dTTP biosynthesis. Thus, dUTPases are involved in maintaining genomic integrity by preventing the uracilation of DNA. Many DNA-containing viruses, which infect mammals also encode for a dUTPase. This review will summarize studies demonstrating that, in addition to their classical enzymatic activity, some dUTPases possess novel functions that modulate the host innate immune response.

Entomopathogenic Viruses in the Neotropics: Current Status and Recently Discovered Species

The market for biological control of insect pests in the world and in Brazil has grown in recent years due to the unwanted ecological and human health impacts of chemical insecticides. Therefore, research on biological control agents for pest management has also increased. For instance, insect viruses have been used to protect crops and forests around the world for decades. Among insect viruses, the baculoviruses are the most studied and used viral biocontrol agent. More than 700 species of insects have been found to be naturally infected by baculoviruses, with 90% isolated from lepidopteran insects. In this review, some basic aspects of baculovirus infection in vivo and in vitro infection, gene content, viral replication will be discussed. Furthermore, we provide examples of the use of insect viruses for biological pest control and recently characterized baculoviruses in Brazil.

Genome Analysis of a Novel Clade II.b Alphabaculovirus Obtained from Artaxa digramma

Artaxa digramma is a lepidopteran pest distributed throughout southern China, Myanmar, Indonesia, and India. Artaxa digramma nucleopolyhedrovirus (ArdiNPV) is a specific viral pathogen of A. digramma and deemed as a promising biocontrol agent against the pest. In this study, the complete genome sequence of ArdiNPV was determined by deep sequencing. The genome of ArdiNPV contains a double-stranded DNA (dsDNA) of 161,734 bp in length and 39.1% G+C content. Further, 149 hypothetical open reading frames (ORFs) were predicted to encode proteins >50 amino acids in length, covering 83% of the whole genome. Among these ORFs, 38 were baculovirus core genes, 22 were lepidopteran baculovirus conserved genes, and seven were unique to ArdiNPV, respectively. No typical baculoviral homologous regions (hrs) were identified in the genome. ArdiNPV had five multi-copy genes including baculovirus repeated ORFs (bros), calcium/sodium antiporter B (chaB), DNA binding protein (dbp), inhibitor of apoptosis protein (iap), and p26. Interestingly, phylogenetic analyses showed that ArdiNPV belonged to Clade II.b of Group II Alphabaculoviruses, which all contain a second copy of dbp. The genome of ArdiNPV was the closest to Euproctis pseudoconspersa nucleopolyhedrovirus, with 57.4% whole-genome similarity. Therefore, these results suggest that ArdiNPV is a novel baculovirus belonging to a newly identified cluster of Clade II.b Alphabaculoviruses.

A Nymphalid-Infecting Group I Alphabaculovirus Isolated from the Major Passion Fruit Caterpillar Pest Dione juno juno (Lepidoptera: Nymphalidae)

Baculoviruses are capable of infecting a wide diversity of insect pests. In the 1990s, the Dione juno nucleopolyhedrovirus (DijuNPV) was isolated from larvae of the major passionfruit defoliator pest Dione juno juno (Nymphalidae) and described at ultrastructural and pathological levels. In this study, the complete genome sequence of DijuNPV was determined and analyzed. The circular genome presents 122,075 bp with a G + C content of 50.9%. DijuNPV is the first alphabaculovirus completely sequenced that was isolated from a nymphalid host and may represent a divergent species. It appeared closely related to Orgyia pseudotsugata multiple nucleopolyhedrovirus (OpMNPV) and other Choristoneura-isolated group I alphabaculoviruses. We annotated 153 open reading frames (ORFs), including a set of 38 core genes, 26 ORFs identified as present in lepidopteran baculoviruses, 17 ORFs unique in baculovirus, and several auxiliary genes (e.g., bro, cathepsin, chitinase, iap-1, iap-2, and thymidylate kinase). The thymidylate kinase (tmk) gene was present fused to a dUTPase (dut) gene in other baculovirus genomes. DijuNPV likely lost the dut portion together with the iap-3 homolog. Overall, the genome sequencing of novel alphabaculoviruses enables a wide understanding of baculovirus evolution.

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

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

Biological and molecular characterization of two Anticarsia gemmatalis multiple nucleopolyhedrovirus clones exhibiting contrasting virulence

Baculovirus natural populations are known to be genetically heterogeneous and such genotypic diversity could have implications in the performance of biocontrol agents. The Anticarsia gemmatalis nucleopolyhedrovirus (AgMNPV) has been widely used to control the velvetbean caterpillar, Anticarsia gemmatalis, in Brazil. In the present work, morphological and molecular analyses as well as the biological activity of AgMNPV genotypes derived from a Brazilian field isolate (AgMNPV-79) were carried out. The existence of genotypic variants in the population was confirmed by DNA restriction analysis. Although difference in virulence was observed among the variants, the most (Ag79-01) and the least (AgL-16) virulent clones do not show any morphological and cytopathological changes when compared to the most studied isolate (AgMNPV-2D). The complete genome analysis of the two viral clones showed the presence of single open reading frames (ORFs) of the pe-38 and he65 genes, which contrasts with the two split ORFs present in the genome of the AgMNPV-2D isolate. The viral clone AgL-16 has many variations in the ie-2 and pe-38 genes, which are transcription regulatory genes responsible for the regulation of viral early gene expression during insect cell infection. Furthermore, other genes showed alterations like the odv-e56, which have an essential role in the maturation and envelopment of the ODVs, and bro-a and bro-b genes which were fused to form a single ORF. For the Ag79-01, although the total number of single nucleotide variants (SNVs) was more prominent in the pe-38 gene, its genome showed very few modifications in comparison to the AgMNPV-2D genome.

The complete genome sequence of a second alphabaculovirus from the true armyworm, Mythimna unipuncta: implications for baculovirus phylogeny and host specificity

The Mythimna unipuncta nucleopolyhedrovirus isolate KY310 (MyunNPV-KY310) is an alphabaculovirus isolated from a true armyworm (Mythimna unipuncta) population in Kentucky, USA. Occlusion bodies of this virus were examined by electron microscopy and the genome sequence was determined by 454 pyrosequencing. MyunNPV-KY310 occlusion bodies consisted of irregular polyhedra measuring 0.8-1.8 µm in diameter and containing multiple virions, with one to six nucleocapsids per virion. The genome sequence was determined to be 156,647 bp with a nucleotide distribution of 43.9% G+C. 152 ORFs and six homologous repeat (hr) regions were annotated for the sequence, including the 38 core genes of family Baculoviridae and an additional group of 26 conserved alphabaculovirus genes. BLAST queries and phylogenetic inference confirmed that MyunNPV-KY310 is most closely related to the alphabaculovirus Leucania separata nucleopolyhedrovirus isolate AH1, which infects Mythimna separata. In contrast, MyunNPV-KY310 did not exhibit a close relationship with Mythimna unipuncta nucleopolyhedrovirus isolate #7, an alphabaculovirus from the same host species. MyunNPV-KY310 lacks the gp64 envelope glycoprotein, which is a characteristic of group II alphabaculoviruses. However, this virus and five other alphabaculoviruses lacking gp64 are placed outside the group I and group II clades in core gene phylogenies, further demonstrating that viruses of genus Alphabaculovirus do not occur in two monophyletic clades. Potential instances of MyunNPV-KY310 ORFs arising by horizontal transfer were detected. Although there are now genome sequences of four different baculoviruses from M. unipuncta, comparison of their genome sequences provides little insight into the genetic basis for their host specificity.

Genome-wide diversity in temporal and regional populations of the betabaculovirus Erinnyis ello granulovirus (ErelGV)

Background

Erinnyis ello granulovirus (ErelGV) is a betabaculovirus infecting caterpillars of the sphingid moth E. ello ello (cassava hornworm), an important pest of cassava crops (Manihot esculenta). In this study, the genome of seven field isolates of the virus ErelGV were deep sequenced and their inter- and intrapopulational sequence diversity were analyzed.

Results

No events of gene gain/loss or translocations were observed, and indels were mainly found within highly repetitive regions (direct repeats, drs). A naturally occurring isolate from Northern Brazil (Acre State, an Amazonian region) has shown to be the most diverse population, with a unique pattern of polymorphisms. Overall, non-synonymous substitutions were found all over the seven genomes, with no specific gathering of mutations on hotspot regions. Independently of their sizes, some ORFs have shown higher levels of non-synonymous changes than others. Non-core genes of known functions and structural genes were among the most diverse ones; and as expected, core genes were the least variable genes. We observed remarkable differences on diversity of paralogous genes, as in multiple copies of p10, fgf, and pep. Another important contrast on sequence diversity was found on genes encoding complex subunits and/or involved in the same biological processes, as late expression factors (lefs) and per os infectivity factors (pifs). Interestingly, several polymorphisms in coding regions lie on sequences encoding specific protein domains.

Conclusions

By comparing and integrating information about inter- and intrapopulational diversity of viral isolates, we provide a detailed description on how evolution operates on field isolates of a betabaculovirus. Our results revealed that 35–41% of the SNPs of ErelGV lead to amino acid changes (non-synonymous substitutions). Some genes, especially non-core genes of unknown functions, tend to accumulate more mutations, while core genes evolve slowly and are more conserved. Additional studies would be necessary to understand the actual effects of such gene variations on viral infection and fitness.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5070-6) contains supplementary material, which is available to authorized users.

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

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