Biodiversity, Evolution and Ecological Specialization of Baculoviruses: A Treasure Trove for Future Applied Research

Assembly and Genome Annotation of Different Strains of Apple Fruit Moth Virus (Cydia pomonella granulovirus)

Cydia pomonella granulovirus is a natural pathogen for Cydia pomonella that is used as a biocontrol agent of insect populations. The study of granulovirus virulence is of particular interest since the development of resistance in natural populations of C. pomonella has been observed during the long-term use of the Mexican isolate CpGV. In our study, we present the genomes of 18 CpGV strains endemic to southern Russia and from Kazakhstan, as well as a strain included in the commercial preparation "Madex Twin", which were sequenced and analyzed. We performed comparative genomic analysis using several tools. From comparisons at the level of genes and protein products that are involved in the infection process of virosis, synonymous and missense substitution variants have been identified. The average nucleotide identity has demonstrated a high similarity with other granulovirus genomes of different geographic origins. Whole-genome alignment of the 18 genomes relative to the reference revealed regions of low similarity. Analysis of gene repertoire variation has shown that BZR GV 4, BZR GV 6, and BZR GV L-7 strains have been the closest in gene content to the commercial "Madex Twin" strain. We have confirmed two deletions using read depth coverage data in regions lacking genes shown by homology analysis for granuloviruses BZR GV L-4 and BZR GV L-6; however, they are not related to the known genes causing viral pathogenicity. Thus, we have isolated novel CpGV strains and analyzed their potential as strains producing highly effective bioinsecticides against C. pomonella.

Evolutionary analysis and biological characterization of a novel alphabaculovirus isolated from Mythimna separata

Baculoviruses are insect-specific pathogens. Novel baculovirus isolates provide new options for the biological control of pests. Therefore, research into the biological characteristics of newly isolated baculoviruses, including accurate classification and nomenclature, is important. In this study, a baculovirus was isolated from Mythimna separata and its complete genome sequence was determined by next-generation sequencing. The double-stranded DNA genome was 153 882 bp in length, encoding 163 open reading frames. The virus was identified as a variant of Mamestra brassicae multiple nucleopolyhedrovirus (MbMNPV) and designated Mamestra brassicae multiple nucleopolyhedrovirus CHN1 (MbMNPV-CHN1) according to ultrastructural analysis, genome comparison and phylogenetic analysis. Phylogenetic inference placed MbMNPV-CHN1 in a clade containing isolates of MacoNPV-A, MacoNPV-B and MbMNPV, which we have designated the Mb-McNPV group. The genomes of isolates in the Mb-McNPV group exhibited a high degree of collinearity with relatively minor differences in the content of annotated open reading frames. The development of codon usage bias in the Mb-McNPV group was affected mainly by natural selection. MbMNPV-CHN1 shows high infectivity against seven species of Lepidoptera. The yield of MbMNPV-CHN1 in the fourth- and fifth-instar M. separata larvae was 6.25×109-1.23×1010 OBs/cadaver. Our data provide insights into the classification, host range and virulence differences among baculoviruses of the Mb-McNPV group, as well as a promising potential new baculoviral insecticide.

The Parapoynx stagnalis Nucleopolyhedrovirus (PastNPV), a Divergent Member of the Alphabaculovirus Group I Clade, Encodes a Homolog of Ran GTPase

We report the analysis of the genome of a novel Alphabaculovirus, Parapoynx stagnalis nucleopolyhedrovirus isolate 473 (PastNPV-473), from cadavers of the rice case bearer, Parapoynx stagnalis Zeller (Lepidoptera: Crambidae), collected in rice fields in Kerala, India. High-throughput sequencing of DNA from PastNPV occlusion bodies and assembly of the data yielded a circular genome-length contig of 114,833 bp with 126 annotated opening reading frames (ORFs) and six homologous regions (hrs). Phylogenetic inference based on baculovirus core gene amino acid sequence alignments indicated that PastNPV is a member of the group I clade of viruses in genus Alphabaculovirus, but different phylogenetic methods yielded different results with respect to the placement of PastNPV and four similarly divergent alphabaculoviruses in the group I clade. Branch lengths and Kimura-2-parameter pairwise nucleotide distances indicated that PastNPV-473 cannot be classified in any of the currently listed species in genus Alphabaculovirus. A unique feature of the PastNPV genome was the presence of an ORF encoding a homolog of Ran GTPase, a regulator of nucleocytoplasmic trafficking. PastNPV appears to have acquired a homolog of Ran relatively recently from a lepidopteran host via horizontal gene transfer.

Biological, morphological, and molecular characterization of the baculovirus PlxyMNPV_LBIV-11, and its virulence towards Plutella xylostella, Trichoplusia ni, and Spodoptera frugiperda larvae

PlxyMNPV_LBIV-11 is an alphabaculovirus strain, isolated from Plutella xylostella larvae. This work characterized this strain at a biological, morphological, and molecular level to evaluate its similarity with other baculoviruses. Its ultrastructure showed a multiple arrangement of nucleocapsids within enveloped virions, all occluded within large cubical polyhedra. PlxyMNPV_LBIV-11 showed infectivity on the Hi5 and Sf9 cell lines, despite these being from heterologous origin. This in vitro infectivity was observed using either BVs or by transfection with genomic DNA. Restriction fragment patterns of PlxyMNPV_LBIV-11, using the enzymes EcoRI, BamHI and HindIII, showed a high relationship with those patterns shown by AcMNPV, except for one or two differential bands with each enzyme. Sequences of core genes lef-8 and lef-9 and the conserved polh gene showed identities ranging from 98 to 100% when compared with those of AcMNPV. Somewhat lower was the sequence identity of the gp64 gene (94%) as compared with those of AcMNPV and PlxyMNPV_CL3, which might be related to the difference in virulence. Besides, the presence of this gene in PlxyMNPV_LBIV-11 indicates that it belongs to group 1 of alphabaculoviruses. A phylogram was estimated with the core and conserved gene sequences, corroborating its high relationship with AcMNPV and PlxyMNPV_CL3. Bioassays were performed with P. xylostella larvae reared on a meridic diet, whose LC₅₀ values indicated lower virulence than AcMNPV when tested against P. xylostella, Spodoptera frugiperda, and Trichoplusia ni larvae. Its virulence against S. frugiperda was only seven times lower than AcMNPV. Its potential as a biological control agent is discussed.

The complete genome sequence of an alphabaculovirus from the brown tussock moth, Olene mendosa Hübner, expands our knowledge of lymantriine baculovirus diversity and evolution

The complete genome sequence was determined for an apparent alphabaculovirus isolated from larval cadavers of the brown tussock moth, Olene mendosa Hübner, collected during an epizootic in Coimbatore, India. The genome was determined to be a circular 142,291 bp molecule, and 147 ORFs and nine homologous regions were annotated for the sequence. Analysis of the sequence confirmed that this virus, Olene mendosa nucleopolyhedrovirus (OlmeNPV), was a member of genus Alphabaculovirus in family Baculoviridae. Phylogenies inferred from nucleotide and amino acid alignments indicated that OlmeNPV was part of a group of viruses that infect moths of genus Lymantria, suggesting that OlmeNPV may have shifted hosts from a Lymantria species to an ancestral Olene species at some point during its evolutionary history. OlmeNPV was most closely related to Lymantria xylina multiple nucleopolyhedrovirus isolate 5 (LyxyMNPV-5). The genomes of OlmeNPV and LyxyMNPV-5 were distinguished not only by differences in ORF content, but by a 27 kbp region of the genome that is inverted in LyxyMNPV-5 relative to OlmeNPV. Pairwise nucleotide distances between OlmeNPV and other Lymantria spp. alphabaculoviruses indicate that OlmeNPV represents a new baculovirus species.

Genomic and biological characteristics of an alphabaculovirus isolated from Trabala vishnou gigantina

The oak lappet moth, Trabala vishnou gigantina is a forest insect pest that damages broad-leaf trees severely. Trabala vishnou gigantina nucleopolyhedrovirus (TrviNPV) has been isolated from a naturally infected T. vishnou gigantina larva and investigated for its biology and the potential to be a biological control agent against its insect host. TrviNPV was characterized by electron microscope of occlusion bodies (OBs), genomic sequencing and field control efficacy. TrviNPV OBs exhibited an irregular polyhedral shape varying in size from 0.99 to 3.99 μm with multiple nucleocapsids per virion. The genome of this virus was 165 657 bp in length with 40.33% GC content and encoded 146 putative ORFs including the 38 baculovirus core genes. TrviNPV is a group II alphabaculovirus that encodes F protein and lacks the gp64 gene specific to group I alphabaculoviruses. Phylogeny and Kimura-2 parameter analysis revealed TrviNPV to be a novel species and closest to ArdiNPV, EupsNPV and OrleNPV. Bioassays and field trials in a shrubland revealed that TrviNPV was virulent and effective to control T. vishnou gigantina in arid semi-desert region. This work firstly reported the whole genome of TrviNPV as well as its biological characters for a possibility to develop this virus as bio-pesticide in the future.

Genome Analysis of an Alphabaculovirus Isolated from the Larch Looper, Erannis ankeraria

The larch looper, Erannis ankeraria Staudinger (Lepidoptera: Geometridae), is one of the major insect pests of larch forests, widely distributed from southeastern Europe to East Asia. A naturally occurring baculovirus, Erannis ankeraria nucleopolyhedrovirus (EranNPV), was isolated from E. ankeraria larvae. This virus was characterized by electron microscopy and by sequencing the whole viral genome. The occlusion bodies (OBs) of EranNPV exhibited irregular polyhedral shapes containing multiple enveloped rod-shaped virions with a single nucleocapsid per virion. The EranNPV genome was 125,247 bp in length with a nucleotide distribution of 34.9% G+C. A total of 131 hypothetical open reading frames (ORFs) were identified, including the 38 baculovirus core genes and five multi-copy genes. Five homologous regions (hrs) were found in the EranNPV genome. Phylogeny and pairwise kimura 2-parameter analysis indicated that EranNPV was a novel group II alphabaculovirus and was most closely related to Apocheima cinerarium NPV (ApciNPV). Field trials showed that EranNPV was effective in controlling E. ankeraria in larch forests. The above results will be relevant to the functional research on EranNPV and promote the use of this virus as a biocontrol agent.

Natural Coinfection between Novel Species of Baculoviruses in Spodoptera ornithogalli Larvae

Spodoptera ornithogalli (Guenée) (Lepidoptera: Noctuidae) is an important pest in different crops of economic relevance in America. For its control, strategies that include chemicals are usually used; so, the description of entomopathogens would be very useful for the formulation of biopesticides. In this regard, two different baculoviruses affecting S. ornithogalli were isolated in Colombia, with one of them being an NPV and the other a GV. Ultrastructural, molecular, and biological characterization showed that both isolates possess the 38 core genes and are novel species in Baculoviridae, named as Spodoptera ornithogalli nucleopolyhedrovirus (SporNPV) and Spodoptera ornithogalli granulovirus (SporGV). The bioassays carried out in larvae of S. ornithogalli and S. frugiperda showed infectivity in both hosts but being higher in the first. In addition, it was observed that SporGV potentiates the insecticidal action of SporNPV (maximum value in ratio 2.5:97.5). Both viruses are individually infective but coexist in nature, producing mixed infections with a synergistic effect that improves the performance of the NPV and enables the transmission of the GV, which presents a slowly killing phenotype.

Monitoring insect transposable elements in large double-stranded DNA viruses reveals host-to-virus and virus-to-virus transposition

The mechanisms by which transposable elements (TEs) can be horizontally transferred between animals are unknown, but viruses are possible candidate vectors. Here, we surveyed the presence of host-derived TEs in viral genomes in 35 deep sequencing data sets produced from 11 host–virus systems, encompassing nine arthropod host species (five lepidopterans, two dipterans, and two crustaceans) and six different double-stranded (ds) DNA viruses (four baculoviruses and two iridoviruses). We found evidence of viral-borne TEs in 14 data sets, with frequencies of viral genomes carrying a TE ranging from 0.01% to 26.33% for baculoviruses and from 0.45% to 7.36% for iridoviruses. The analysis of viral populations separated by a single replication cycle revealed that viral-borne TEs originating from an initial host species can be retrieved after viral replication in another host species, sometimes at higher frequencies. Furthermore, we detected a strong increase in the number of integrations in a viral population for a TE absent from the hosts’ genomes, indicating that this TE has undergone intense transposition within the viral population. Finally, we provide evidence that many TEs found integrated in viral genomes (15/41) have been horizontally transferred in insects. Altogether, our results indicate that multiple large dsDNA viruses have the capacity to shuttle TEs in insects and they underline the potential of viruses to act as vectors of horizontal transfer of TEs. Furthermore, the finding that TEs can transpose between viral genomes of a viral species sets viruses as possible new niches in which TEs can persist and evolve.

Mixtures of Insect-Pathogenic Viruses in a Single Virion: towards the Development of Custom-Designed Insecticides

Alphabaculoviruses (Baculoviridae) are pathogenic DNA viruses of Lepidoptera that have applications as the basis for biological insecticides and expression vectors in biotechnological processes. These viruses have a characteristic physical structure that facilitates the transmission of groups of genomes. We demonstrate that coinfection of a susceptible insect by two different alphabaculovirus species results in the production of mixed-virus occlusion bodies containing the parental viruses. This occurred between closely related and phylogenetically more distant alphabaculoviruses. Approximately half the virions present in proteinaceous viral occlusion bodies produced following coinfection of insects with a mixture of two alphabaculoviruses contained both viruses, indicating that the viruses coinfected and replicated in a single cell and were coenveloped within the same virion. This observation was confirmed by endpoint dilution assay. Moreover, both viruses persisted in the mixed-virus population by coinfection of insects during several rounds of insect-to-insect transmission. Coinfection by viruses that differed in genome size had unexpected results on the length of viral nucleocapsids, which differed from those of both parental viruses. These results have unique implications for the development of alphabaculoviruses as biological control agents of insect pests.IMPORTANCE Alphabaculoviruses are used as biological insecticides and expression vectors in biotechnology and medical applications. We demonstrate that in caterpillars infected with particular mixtures of viruses, the genomes of different baculovirus species can be enveloped together within individual virions and occluded within proteinaceous occlusion bodies. This results in the transmission of mixed-virus populations to the caterpillar stages of moth species. Once established, mixed-virus populations persist by coinfection of insect cells during several rounds of insect-to-insect transmission. Mixed-virus production technology opens the way to the development of custom-designed insecticides for control of different combinations of caterpillar pest species.

Comparative genomic analysis of three geographical isolates from China reveals high genetic stability of Plutella xylostella granulovirus

In this study, the genomes of three Plutella xylostella granulovirus (PlxyGV) isolates, PlxyGV-W and PlxyGV-Wn from near Wuhan and PlxyGV-B from near Beijing, China were completely sequenced and comparatively analyzed to investigate genetic stability and diversity of PlxyGV. PlxyGV-W, PlxyGV-B and PlxyGV-Wn consist of 100,941bp, 100,972bp and 100,999bp in length with G + C compositions of 40.71–40.73%, respectively, and share nucleotide sequence identities of 99.5–99.8%. The three individual isolates contain 118 putative protein-encoding ORFs in common. PlxyGV-W, PlxyGV-B and PlxyGV-Wn have ten, nineteen and six nonsynonymous intra isolate nucleotide polymorphisms (NPs) in six, fourteen and five ORFs, respectively, including homologs of five DNA replication/late expression factors and two per os infectivity factors. There are seventeen nonsynonymous inter isolate NPs in seven ORFs between PlxyGV-W and PlxyGV-B, seventy three nonsynonymous NPs in forty seven ORFs between PlxyGV-W and PlxyGV-Wn, seventy seven nonsynonymous NPs in forty six ORFs between PlxyGV-B and PlxyGV-Wn. Alignment of the genome sequences of nine PlxyGV isolates sequenced up to date shows that the sequence homogeneity between the genomes are over 99.4%, with the exception of the genome of PlxyGV-SA from South Africa, which shares a sequence identity of 98.6–98.7% with the other ones. No events of gene gain/loss or translocations were observed. These results suggest that PlxyGV genome is fairly stable in nature. In addition, the transcription start sites and polyadenylation sites of thirteen PlxyGV-specific ORFs, conserved in all PlxyGV isolates, were identified by RACE analysis using mRNAs purified from larvae infected by PlxyGV-Wn, proving the PlxyGV-specific ORFs are all genuine genes.

Genomic analyses of Biston suppressaria nucleopolyhedrovirus: a viral isolate obtained from the tea looper caterpillar, Biston suppressaria (Guenée, 1857)

We described the complete genome sequence of a novel baculovirus isolate of species Buzura suppressaria nucleopolyhedrovirus, called by isolate CNPSo-25. The occlusion bodies were found to be polyhedral in shape and to contain virions with singly embedded nucleocapsids. The size of the genome is 121,377 bp with a G+C content of 36.7%. We annotated 131 ORFs that cover 90.42% of the genome. Moreover, phylogenetic inference indicated that CNPSo-25 is a member of genus Alphabaculovirus that clustered together with two other Chinese isolates of the same species. We called the virus by Biston suppressaria nucleopolyhedrovirus isolate CNPSo-25 (BisuNPV-CNPSo-25), as Buzura was placed inside the lepidopteran genus Biston. As expected, we detected intra-population variability in the virus sample when the novel isolate was compared to the Chinese isolates: 292 single nucleotide variants were found in the genome, with 181 affecting the protein product. The closest representatives of other species to BisuNPV-CNPSo-25 was found to be Sucra jujuba nucleopolyhedrovirus and Hyposidra talaca nucleopolyhedrovirus, two other virus isolates of geometrid caterpillars. The study of baculovirus genomes is of importance for the development of tools for insect pest biological control and biotechnology.

Mosaic genome evolution and phylogenetics of Chrysodeixis includens nucleopolyhedrovirus (ChinNPV) and virulence of seven new isolates from the Brazilian states of Minas Gerais and Mato Grosso

In a comparative analysis of genome sequences from isolates of the baculovirus Chrysodeixis includens nucleopolyhedrovirus (ChinNPV) from Brazil and Guatemala, we identified a subset of isolates possessing chimeric genomes. We identified six distinct phylogenetically incongruous regions (PIRs) dispersed in the genomes, of between 279 and 3345 bp in length. The individual PIRs possessed high sequence similarity among the affected ChinNPV isolates but varied in coverage in some instances. The donor for four of the PIRs implicated in horizontal gene transfer (HGT) was identified as Trichoplusia ni single nucleopolyhedrovirus (TnSNPV), an alphabaculovirus closely related to ChinNPV, or another unknown but closely related virus. BLAST searches of the other two PIRs returned only ChinNPV sequences, but HGT from an unknown donor baculovirus cannot be excluded. Although Chrysodeixis includens and Trichoplusia ni are frequently co-collected from soybean fields in Brazil, pathogenicity data suggest that natural coinfection of C. includens larvae with ChinNPV and TnSNPV is probably uncommon. Additionally, since the chimeric ChinNPV genomes with tracts of TnSNPV sequence were restricted to a single monophyletic lineage of closely related isolates, a model of progressive restoration of the native DNA sequence by recombination with ChinNPV possessing a fully or partially non-chimeric genome is reasonable. However, multiple independent HGT from TnSNPV to ChinNPV during the evolution of these isolates cannot be excluded. Mortality data suggest that the ChinNPV isolates with chimeric genomes are not significantly different in pathogenicity towards C. includens when compared to most other ChinNPV isolates. Exclusion of the PIRs prior to phylogenetic analysis had a large impact on the topology of part of the maximum-likelihood tree, revealing a homogenous clade of three isolates (IB, IC and ID) from Paraná state in Brazil collected in 2006, together with an isolate from Guatemala collected in 1972 (IA), comprising the lineage uniquely affected by HGT from TnSNPV. The other 10 Brazilian ChinNPV isolates from Paraná, Mato Grosso, and Minas Gerais states showed higher variability, where only three isolates from Paraná state formed a monophyletic group correlating with geographical origin.

Advances in the Bioinformatics Knowledge of mRNA Polyadenylation in Baculovirus Genes

Baculoviruses are a group of insect viruses with large circular dsDNA genomes exploited in numerous biotechnological applications, such as the biological control of agricultural pests, the expression of recombinant proteins or the gene delivery of therapeutic sequences in mammals, among others. Their genomes encode between 80 and 200 proteins, of which 38 are shared by all reported species. Thanks to multi-omic studies, there is remarkable information about the baculoviral proteome and the temporality in the virus gene expression. This allows some functional elements of the genome to be very well described, such as promoters and open reading frames. However, less information is available about the transcription termination signals and, consequently, there are still imprecisions about what are the limits of the transcriptional units present in the baculovirus genomes and how is the processing of the 3′ end of viral mRNA. Regarding to this, in this review we provide an update about the characteristics of DNA signals involved in this process and we contribute to their correct prediction through an exhaustive analysis that involves bibliography information, data mining, RNA structure and a comprehensive study of the core gene 3′ ends from 180 baculovirus genomes.

Whole-genome-based phylogeny of African swine fever virus

Aim

A genome-scale phylogenetic analysis was used to infer the evolutionary dynamics of Asfarviridae - African swine fever virus (ASFV) - and better define its genetic diversity.

Materials and Methods

All complete ASFV genomes from NCBI's resource as of March 2020 were used. The phylogenetic analysis used maximum likelihood and neighbor-joining methods. The evolutionary models detection was done with the help of the package of programs MEGA-X. Algorithms were used to build phylogenetic trees for type B DNA polymerases of ASFV (n=34) and HcDNAV (n=2), as an external group.

Results

An expedient categorization of the Asfarviridae family uses five clades. Genotype 1 (except for LIV 5/40 virus isolate) as well genotype 7 are assigned to the alpha clade; genotype 2 to the beta clade; genotypes 8, 9, and 10 to the gamma clade; genotype 5 to the delta clade; and genotypes 3, 4, and 20, as well as genotype 22 and the LIV 5/40 isolate to the epsilon clade. Branch lengths on the phylogenetic tree are proportional to genetic distance along the branch. Branches at the phylogenetic tree of Asfarviridae are much shorter than branches for Baculoviridae. Shorter branches in ASFVs population suggest that Asfarviridae evolved relatively recently and remain more closely related.

Conclusion

We suggest applying more robust standards using whole genomes to ensure the correct classification of ASFV and maintain phylogeny as a useful tool.

Identification and Characterization of the Nucleolar Localization Signal of Autographa californica Multiple Nucleopolyhedrovirus LEF5

Autographa californica multiple nucleopolyhedrovirus (AcMNPV) late expression factor 5 (LEF5) is highly conserved in all sequenced baculovirus genomes and plays an important role in production of infectious viral progeny. In this study, nucleolar localization of AcMNPV LEF5 was characterized. Through transcriptome analysis, we identified two putative nucleolar proteins, Spodoptera frugiperda nucleostemin (SfNS) and fibrillarin (SfFBL), from Sf9 cells. Immunofluorescence analysis demonstrated that SfNS and SfFBL were localized to the nucleolus. AcMNPV infection resulted in reorganization of the nucleoli of infected cells. Colocalization of LEF5 and SfNS showed that AcMNPV LEF5 was localized to the nucleolus in Sf9 cells. Bioinformatic analysis revealed that basic amino acids of LEF5 are enriched at residues 184 to 213 and may contain a nucleolar localization signal (NoLS). Green fluorescent protein (GFP) fused to NoLS of AcMNPV LEF5 localized to the nucleoli of transfected cells. Multiple-point mutation analysis demonstrated that amino acid residues 197 to 204 are important for nucleolar localization of LEF5. To identify whether the NoLS in AcMNPV LEF5 is important for production of viral progeny, a lef5-null AcMNPV bacmid was constructed; several NoLS-mutated LEF5 proteins were reinserted into the lef5-null AcMNPV bacmid with a GFP reporter. The constructs containing point mutations at residues 185 to 189 or 197 to 204 in AcMNPV LEF5 resulted in reduction in production of infectious viral progeny and occlusion body yield in bacmid-transfected cells. Together, these data suggested that AcMNPV LEF5 contains an NoLS, which is important for nucleolar localization of LEF5, progeny production, and occlusion body production.IMPORTANCE Many viruses, including human and plant viruses, target nucleolar functions as part of their infection strategy. However, nucleolar localization for baculovirus proteins has not yet been characterized. In this study, two nucleolar proteins, SfNS and SfFBL, were identified in Sf9 cells. Our results showed that Autographa californica multiple nucleopolyhedrovirus (AcMNPV) infection resulted in redistribution of the nucleoli of infected cells. We demonstrated that AcMNPV late expression factor 5 (LEF5) could localize to the nucleolus and contains a nucleolar localization signal (NoLS), which is important for nucleolar localization of AcMNPV LEF5 and for production of viral progeny and yield of occlusion bodies.

Automated, phylogeny-based genotype delimitation of the Hepatitis Viruses HBV and HCV

Background

The classification of hepatitis viruses still predominantly relies on ad hoc criteria, i.e., phenotypic traits and arbitrary genetic distance thresholds. Given the subjectivity of such practices coupled with the constant sequencing of samples and discovery of new strains, this manual approach to virus classification becomes cumbersome and impossible to generalize.

Methods

Using two well-studied hepatitis virus datasets, HBV and HCV, we assess if computational methods for molecular species delimitation that are typically applied to barcoding biodiversity studies can also be successfully deployed for hepatitis virus classification. For comparison, we also used ABGD, a tool that in contrast to other distance methods attempts to automatically identify the barcoding gap using pairwise genetic distances for a set of aligned input sequences.

Results—Discussion

We found that the mPTP species delimitation tool identified even without adapting its default parameters taxonomic clusters that either correspond to the currently acknowledged genotypes or to known subdivision of genotypes (subtypes or subgenotypes). In the cases where the delimited cluster corresponded to subtype or subgenotype, there were previous concerns that their status may be underestimated. The clusters obtained from the ABGD analysis differed depending on the parameters used. However, under certain values the results were very similar to the taxonomy and mPTP which indicates the usefulness of distance based methods in virus taxonomy under appropriate parameter settings. The overlap of predicted clusters with taxonomically acknowledged genotypes implies that virus classification can be successfully automated.

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.

Genetic Variation and Biological Activity of Two Closely Related Alphabaculoviruses during Serial Passage in Permissive and Semi-Permissive Heterologous Hosts

Phylogenetic analyses suggest that Mamestra brassicae multiple nucleopolyhedrovirus (MbMNPV) and Helicoverpa armigera multiple nucleopolyhedrovirus (HearMNPV) may be strains of the same virus species. Most of the studies comparing their biological activities have been performed in their homologous hosts. A comparison of host range and stability in alternative hosts was performed. The host range of these viruses was compared using high concentrations of inoculum to inoculate second instars of six species of Lepidoptera. One semi-permissive host (Spodoptera littoralis) and one permissive host (S. exigua) were then selected and used to perform six serial passages involving a concentration corresponding to the ~25% lethal concentration for both viruses. Restriction endonuclease analysis showed fragment length polymorphisms in every host-virus system studied. In S. littoralis, serial passage of MbMNPV resulted in decreased pathogenicity and an increase in speed-of-kill, whereas no significant changes were detected for HearMNPV with respect to the initial inoculum. In contrast, both viruses showed a similar trend in S. exigua. These results highlight the low genetic diversity and a high phenotypic stability of HearMNPV with respect to the original inoculum after six successive passages in both insect hosts. This study concludes that host-baculovirus interactions during serial passage are complex and the process of adaptation to a novel semi-permissive host is far from predictable.

Identification of Multiple Replication Stages and Origins in the Nucleopolyhedrovirus of Anticarsia gemmatalis

To understand the mechanism of replication used by baculoviruses, it is essential to describe all the factors involved, including virus and host proteins and the sequences where DNA synthesis starts. A lot of work on this topic has been done, but there is still confusion in defining what sequence/s act in such functions, and the mechanism of replication is not very well understood. In this work, we performed an AgMNPV replication kinetics into the susceptible UFL-Ag-286 cells to estimate viral genome synthesis rates. We found that the viral DNA exponentially increases in two different phases that are temporally separated by an interval of 5 h, probably suggesting the occurrence of two different mechanisms of replication. Then, we prepared a plasmid library containing virus fragments (0.5-2 kbp), which were transfected and infected with AgMNPV in UFL-Ag-286 cells. We identified 12 virus fragments which acted as origins of replication (ORI). Those fragments are in close proximity to core genes. This association to the core genome would ensure vertical transmission of ORIs. We also predict the presence of common structures on those fragments that probably recruit the replication machinery, a structure also present in previously reported ORIs in baculoviruses.

A Novel Alphabaculovirus from the Soybean Looper, Chrysodeixis includens, that Produces Tetrahedral Occlusion Bodies and Encodes Two Copies of he65

Isolates of the alphabaculovirus species, Chrysodeixis includens nucleopolyhedrovirus, have been identified that produce polyhedral occlusion bodies and infect larvae of the soybean looper, Chrysodeixis includens. In this study, we report the discovery and characterization of a novel C. includens-infecting alphabaculovirus, Chrysodeixis includens nucleopolyhedrovirus #1 (ChinNPV#1), that produces tetrahedral occlusion bodies. In bioassays against C. includens larvae, ChinNPV #1 exhibited a degree of pathogenicity that was similar to that of other ChinNPV isolates, but killed larvae more slowly. The host range of ChinNPV#1 was found to be very narrow, with no indication of infection occurring in larvae of Trichoplusia ni and six other noctuid species. The ChinNPV#1 genome sequence was determined to be 130,540 bp, with 126 open reading frames (ORFs) annotated but containing no homologous repeat (hr) regions. Phylogenetic analysis placed ChinNPV#1 in a clade with other Group II alphabaculoviruses from hosts of lepidopteran subfamily Plusiinae, including Chrysodeixis chalcites nucleopolyhedrovirus and Trichoplusia ni single nucleopolyhedrovirus. A unique feature of the ChinNPV#1 genome was the presence of two full-length copies of the he65 ORF. The results indicate that ChinNPV#1 is related to, but distinct from, other ChinNPV isolates.

Genetic Variability of Chrysodeixis Includens Nucleopolyhedrovirus (ChinNPV) and the Insecticidal Characteristics of Selected Genotypic Variants

Genetic variation in baculoviruses is recognized as a key factor, not only due to the influence of such variation on pathogen transmission and virulence traits, but also because genetic variants can form the basis for novel biological insecticides. In this study, we examined the genetic variability of Chrysodeixis includens nucleopolyhedrovirus (ChinNPV) present in field isolates obtained from virus-killed larvae. Different ChinNPV strains were identified by restriction endonuclease analysis, from which genetic variants were isolated by plaque assay. Biological characterization studies were based on pathogenicity, median time to death (MTD), and viral occlusion body (OB) production (OBs/larva). Nine different isolates were obtained from eleven virus-killed larvae collected from fields of soybean in Mexico. An equimolar mixture of these isolates, named ChinNPV-Mex1, showed good insecticidal properties and yielded 23 genetic variants by plaque assay, one of which (ChinNPV-R) caused the highest mortality in second instars of C. includens. Five of these variants were selected: ChinNPV-F, ChinNPV-J, ChinNPV-K, ChinNPV-R, and ChinNPV-V. No differences in median time to death were found between them, while ChinNPV-F, ChinNPV-K, ChinNPV-R and ChinNPV-V were more productive than ChinNPV-J and the original mixture of field isolates ChinNPV-Mex1. These results demonstrate the high variability present in natural populations of this virus and support the use of these new genetic variants as promising active substances for baculovirus-based bioinsecticides.

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.

The complete genome sequence of an alphabaculovirus from Spodoptera exempta, an agricultural pest of major economic significance in Africa

Spodoptera exempta nucleopolyhedrovirus (SpexNPV) is a viral pathogen of the African armyworm, Spodoptera exempta (Lepidoptera: Noctuidae), a significant agricultural pest of cereal crops in Africa. SpexNPV has been evaluated as a potential insecticide for control of this pest and has served as the subject of research on baculovirus pathology and transmission. Occlusion bodies (OBs) of SpexNPV isolate 244.1 were examined, and the nucleotide sequence of the genome was determined and characterized. SpexNPV-244.1 OBs consisted of irregular polyhedra with a size and appearance typical for alphabaculoviruses. Virions within the polyhedra contained 1–8 nucleocapsids per unit envelope. The SpexNPV-244.1 genome was comprised of a 129,528 bp circular sequence, in which 139 ORFs were annotated. Five homologous regions (hrs) consisting of a variable number of 28-bp imperfect palindromes were identified in the genome. The genome sequence contained the 38 core genes of family Baculoviridae, as well as three ORFs unique to the SpexNPV sequence and one ORF that was apparently acquired by horizontal gene transfer with a betabaculovirus ancestor. Phylogenetic inference with core gene amino acid sequence alignments placed SpexNPV-244.1 in a lineage containing alphabaculoviruses of Spodoptera frugiperda and Spodopotera exigua which in turn is part of a larger group of alphabaculoviruses from the subfamily Noctuinae in the lepidopteran family Noctuidae. Kimura-2-parameter pairwise nucleotide distances indicated that SpexNPV-244.1 represented a different and previously unlisted species in the genus Alphabaculovirus. Gene parity plots indicated that the gene order of SpexNPV-244.l was extensively collinear with that of Spodoptera exigua NPV (SeMNPV). These plots also revealed a group of 17 core genes whose order was conserved in other alpha- and betabaculoviruses.

Viral and metazoan poxins are cGAMP-specific nucleases that restrict cGAS-STING signalling

Cytosolic DNA triggers innate immune responses through the activation of cyclic GMP-AMP synthase (cGAS) and production of the cyclic dinucleotide second messenger 2',3'-cyclic GMP-AMP (cGAMP)1-4. 2',3'-cGAMP is a potent inducer of immune signalling; however, no intracellular nucleases are known to cleave 2',3'-cGAMP and prevent the activation of the receptor stimulator of interferon genes (STING)5-7. Here we develop a biochemical screen to analyse 24 mammalian viruses, and identify poxvirus immune nucleases (poxins) as a family of 2',3'-cGAMP-degrading enzymes. Poxins cleave 2',3'-cGAMP to restrict STING-dependent signalling and deletion of the poxin gene (B2R) attenuates vaccinia virus replication in vivo. Crystal structures of vaccinia virus poxin in pre- and post-reactive states define the mechanism of selective 2',3'-cGAMP degradation through metal-independent cleavage of the 3'-5' bond, converting 2',3'-cGAMP into linear Gp[2'-5']Ap[3']. Poxins are conserved in mammalian poxviruses. In addition, we identify functional poxin homologues in the genomes of moths and butterflies and the baculoviruses that infect these insects. Baculovirus and insect host poxin homologues retain selective 2',3'-cGAMP degradation activity, suggesting an ancient role for poxins in cGAS-STING regulation. Our results define poxins as a family of 2',3'-cGAMP-specific nucleases and demonstrate a mechanism for how viruses evade innate immunity.

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.

Autographa Californica Multiple Nucleopolyhedrovirus Enters Host Cells via Clathrin-Mediated Endocytosis and Direct Fusion with the Plasma Membrane

The cell entry mechanism of Autographa californica multiple nucleopolyhedrovirus (AcMNPV) is not fully understood. Previous studies showed that AcMNPV entered host cells primarily through clathrin-mediated endocytosis, and could efficiently infect cells via fusion with the plasma membrane after a low-pH trigger. However, whether AcMNPV enters cells via these two pathways simultaneously, and the exact manner in which AcMNPV particles are internalized into cells remains unclear. In this study, using single-virus tracking, we observed that AcMNPV particles were first captured by pre-existing clathrin-coated pits (CCP), and were then delivered to early endosomes. Population-based analysis of single-virus tracking and quantitative electron microscopy demonstrated that the majority of particles were captured by CCPs and internalized via invagination. In contrast, a minority of virus particles were not delivered to CCPs, and were internalized through direct fusion with the plasma membrane without invagination. Quantitative electron microscopy also showed that, while inhibition of CCP assembly significantly impaired viral internalization, inhibition of endosomal acidification blocked virus particles out of vesicles. Collectively, these findings demonstrated that approximately 90% of AcMNPV particles entered cells through clathrin-mediated endocytosis and 10% entered via direct fusion with the plasma membrane. This study will lead toward a better understanding of AcMNPV infection.