Identification of a conserved non-protein-coding genomic element that plays an essential role in Alphabaculovirus pathogenesis

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.

Nucleocapsid Assembly of Baculoviruses

The baculovirus nucleocapsid is formed through a rod-like capsid encapsulating a genomic DNA molecule of 80~180 kbp. The viral capsid is a large oligomer composed of many copies of various protein subunits. The assembly of viral capsids is a complex oligomerization process. The timing of expression of nucleocapsid-related proteins, transport pathways, and their interactions can affect the assembly process of preformed capsids. In addition, the selection of viral DNA and the injection of the viral genome into empty capsids are the critical steps in nucleocapsid assembly. This paper reviews the replication and recombination of baculovirus DNA, expression and transport of capsid proteins, formation of preformed capsids, DNA encapsulation, and nucleocapsid formation. This review will provide a basis for further study of the nucleocapsid assembly mechanism of baculovirus.

Genomic sequencing of Troides aeacus nucleopolyhedrovirus (TraeNPV) from golden birdwing larvae (Troides aeacus formosanus) to reveal defective Autographa californica NPV genomic features

Background

The golden birdwing butterfly (Troides aeacus formosanus) is a rarely observed species in Taiwan. Recently, a typical symptom of nuclear polyhedrosis was found in reared T. aeacus larvae. From the previous Kimura-2 parameter (K-2-P) analysis based on the nucleotide sequence of three genes in this isolate, polh, lef-8 and lef-9, the underlying virus did not belong to any known nucleopolyhedrovirus (NPV) species. Therefore, this NPV was provisionally named “TraeNPV”. To understand this NPV, the nucleotide sequence of the whole TraeNPV genome was determined using next-generation sequencing (NGS) technology.

Results

The genome of TraeNPV is 125,477 bp in length with 144 putative open reading frames (ORFs) and its GC content is 40.45%. A phylogenetic analysis based on the 37 baculoviral core genes suggested that TraeNPV is a Group I NPV that is closely related to Autographa californica nucleopolyhedrovirus (AcMNPV). A genome-wide analysis showed that TraeNPV has some different features in its genome compared with other NPVs. Two novel ORFs (Ta75 and Ta139), three truncated ORFs (pcna, he65 and bro) and one duplicated ORF (38.7 K) were found in the TraeNPV genome; moreover, there are fewer homologous regions (hrs) than there are in AcMNPV, which shares eight hrs within the TraeNPV genome. TraeNPV shares similar genomic features with AcMNPV, including the gene content, gene arrangement and gene/genome identity, but TraeNPV lacks 15 homologous ORFs from AcMNPV in its genome, such as ctx, host cell-specific factor 1 (hcf-1), PNK/PNL, vp15, and apsup, which are involved in the auxiliary functions of alphabaculoviruses.

Conclusions

Based on these data, TraeNPV would be clarified as a new NPV species with defective AcMNPV genomic features. The precise relationship between TraeNPV and other closely related NPV species were further investigated. This report could provide comprehensive information on TraeNPV for evolutionary insights into butterfly-infected NPV.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5713-2) contains supplementary material, which is available to authorized users.

Comprehensive analysis of single molecule sequencing-derived complete genome and whole transcriptome of Hyposidra talaca nuclear polyhedrosis virus

We sequenced the Hyposidra talaca NPV (HytaNPV) double stranded circular DNA genome using PacBio single molecule sequencing technology. We found that the HytaNPV genome is 139,089 bp long with a GC content of 39.6%. It encodes 141 open reading frames (ORFs) including the 37 baculovirus core genes, 25 genes conserved among lepidopteran baculoviruses, 72 genes known in baculovirus, and 7 genes unique to the HytaNPV genome. It is a group II alphabaculovirus that codes for the F protein and lacks the gp64 gene found in group I alphabaculovirus viruses. Using RNA-seq, we confirmed the expression of the ORFs identified in the HytaNPV genome. Phylogenetic analysis showed HytaNPV to be closest to BusuNPV, SujuNPV and EcobNPV that infect other tea pests, Buzura suppressaria, Sucra jujuba, and Ectropis oblique, respectively. We identified repeat elements and a conserved non-coding baculovirus element in the genome. Analysis of the putative promoter sequences identified motif consistent with the temporal expression of the genes observed in the RNA-seq data.

The Autographa californica Multiple Nucleopolyhedrovirus ac83 Gene Contains a cis-Acting Element That Is Essential for Nucleocapsid Assembly

Baculoviridae is a family of insect-specific viruses that have a circular double-stranded DNA genome packaged within a rod-shaped capsid. The mechanism of baculovirus nucleocapsid assembly remains unclear. Previous studies have shown that deletion of the ac83 gene of Autographa californica multiple nucleopolyhedrovirus (AcMNPV) blocks viral nucleocapsid assembly. Interestingly, the ac83-encoded protein Ac83 is not a component of the nucleocapsid, implying a particular role for ac83 in nucleocapsid assembly that may be independent of its protein product. To examine this possibility, Ac83 synthesis was disrupted by insertion of a chloramphenicol resistance gene into its coding sequence or by deleting its promoter and translation start codon. Both mutants produced progeny viruses normally, indicating that the Ac83 protein is not required for nucleocapsid assembly. Subsequently, complementation assays showed that the production of progeny viruses required the presence of ac83 in the AcMNPV genome instead of its presence in trans Therefore, we reasoned that ac83 is involved in nucleocapsid assembly via an internal cis-acting element, which we named the nucleocapsid assembly-essential element (NAE). The NAE was identified to lie within nucleotides 1651 to 1850 of ac83 and had 8 conserved A/T-rich regions. Sequences homologous to the NAE were found only in alphabaculoviruses and have a conserved positional relationship with another essential cis-acting element that was recently identified. The identification of the NAE may help to connect the data of viral cis-acting elements and related proteins in the baculovirus nucleocapsid assembly, which is important for elucidating DNA-protein interaction events during this process.IMPORTANCE Virus nucleocapsid assembly usually requires specific cis-acting elements in the viral genome for various processes, such as the selection of the viral genome from the cellular nucleic acids, the cleavage of concatemeric viral genome replication intermediates, and the encapsidation of the viral genome into procapsids. In linear DNA viruses, such elements generally locate at the ends of the viral genome; however, most of these elements remain unidentified in circular DNA viruses (including baculovirus) due to their circular genomic conformation. Here, we identified a nucleocapsid assembly-essential element in the AcMNPV (the archetype of baculovirus) genome. This finding provides an important reference for studies of nucleocapsid assembly-related elements in baculoviruses and other circular DNA viruses. Moreover, as most of the previous studies of baculovirus nucleocapsid assembly have been focused on viral proteins, our study provides a novel entry point to investigate this mechanism via cis-acting elements in the viral genome.

The genome of Dasychira pudibunda nucleopolyhedrovirus (DapuNPV) reveals novel genetic connection between baculoviruses infecting moths of the Lymantriidae family

Background

DapuNPV (Dasychira pudibunda nucleopolyhedrovirus), presented in this report, belongs to Alphabaculovirus group Ib. Its full, newly sequenced genome shows close relationship to baculovirus OpMNPV isolated from douglas-fir tussock moth Orgyia pseudotsugata. Baculovirus DapuNPV is a natural limiter of pale tussock moth Dasychira pudibunda L. (syn. Calliteara pudibunda L.)(Lepidoptera, Lymantriidae), which can occur in a form of an outbreak on many species of deciduous trees and may cause significant economic losses in the forests.

Methods

Late instars dead larvae of pale tussock moth were mechanically homogenized and polyhedra were purified during series of ultracentrifugation. Viral DNA was extarcted and sequenced using Miseq Illumina platform. 294,902 paired reads were used for de novo assembling. Genome annotation, multiple allingment to others baculoviruses and phylogegentic analises were perform with the use of multiple bioinformatic tools like: Glimmer3, HMMER web server, Geneious 7 and MEGA6.

Results

The genome of DapuNPV is 136,761 bp long with AT pairs content 45.6 %. The predicted number of encoded putative open reading frames (ORFs) is 161 and six of them demonstrate low or no homology to ORFs previously found in baculoviruses. DapuNPV genome shows very high similarity to OpMNPV in a nucleotide sequence (91.1 % of identity) and gene content (150 homologous ORFs), though some major differences (e.g. lack of he65 in OpMNPV) have also been noted.

Conclusions

Similarly to other members of the Baculoviridae family, DapuNPV baculovirus possesses highly conserved core genes. Among them, there is a second copy of occluded derived virus envelope 27 protein (odv-e27), which was previously found only in a member of Alphabaculovirus group II – LyxyMNPV (Lymantria xylina MNPV). Surprisingly enough, DapuNPV and LyxyMNPV genomes share also another feature. Phylogenetic analysis of chitin binding family protein (cbpl) indicates significant similarity of those two baculoviruses from distinct evolutionary groups which infect the same hosts from Lymantriidae. The ubiquitin like family gene (ubil), which has not been described until now, is another characteristic component of DapuNPV genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1963-9) contains supplementary material, which is available to authorized users.