Structural divergence among genomes of closely related baculoviruses and its implications for baculovirus evolution

Genome comparison of long-circulating field CnmeGV isolates from the same region

Cnaphalocrocis medinalis granulovirus (CnmeGV), belonging to Betabaculovirus cnamedinalis, can infect the rice pest, the rice leaf roller. In 1979, a CnmeGV isolate, CnmeGV-EP, was collected from Enping County, China. In 2014, we collected another CnmeGV isolate, CnmeGV-EPDH3, at the same location and obtained the complete virus genome sequence using Illumina and ONT sequencing technologies. By combining these two virus isolates, we updated the genome annotation of CnmeGV and conducted an in-depth analysis of its genome features. CnmeGV genome contains abundant tandem repeat sequences, and the repeating units in the homologous regions (hrs) exhibit overlapping and nested patterns. The genetic variations within EPDH3 population show the high stability of CnmeGV genome, and tandem repeats are the only region of high genetic variation in CnmeGV genome replication. Some defective viral genomes formed by recombination were found within the population. Comparison analysis of the two virus isolates collected from Enping showed that the proteins encoded by the CnmeGV-specific genes were less conserved relative to the baculovirus core genes. At the genomic level, there are a large number of SNPs and InDels between the two virus isolates, especially in and around the bro genes and hrs. Additionally, we discovered that CnmeGV acquired a segment of non-ORF sequence from its host, which does not provide any new proteins but rather serves as redundant genetic material integrated into the viral genome. Furthermore, we observed that the host's transposon piggyBac has inserted into some virus genes. Together, dsDNA viruses could acquire non-coding genetic material from their hosts to expand the size of their genomes. These findings provide new insights into the evolution of dsDNA viruses.

NISES-AnPe-428 cell line derived from the Chinese oak silkworm Antheraea pernyi is permissive for multiple nucleopolyhedrovirus species from insects of four different families

The cell line NISES-AnPe-428 (AnPe), derived from the Chinese oak silkworm Antheraea pernyi, was characterized for its permissiveness and productivity for six different nucleopolyhedrovirus (NPV) species. These NPVs included homologous Antheraea pernyi NPV (AnpeNPV) and heterologous Autographa californica multiple NPV (AcMNPV), Bombyx mori NPV (BmNPV), Hyphantria cunea MNPV (HycuMNPV), Spodoptera exigua MNPV (SeMNPV), and Lymantria dispar MNPV (LdMNPV), representing viruses that had been isolated from insect species belonging to five different families (Saturniidae, Noctuidae, Bombycidae, Arctiidae, and Lymantriidae). We found that AnPe cells supported productive replication of AnpeNPV, AcMNPV, BmNPV, HycuMNPV, and SeMNPV to varying degrees. Upon infection with SeMNPV, a subset of AnPe cell population in the culture underwent apoptosis, while remaining cells produced limited amounts of progeny virions and polyhedra. AnPe cells were refractory to LdMNPV infection and failed to support replication of viral DNA, indicating that viral replication was restricted at or prior to the step of viral DNA replication. These results indicated that AnPe cells have the potential to provide excellent systems for studying the molecular mechanisms underlying cellular permissiveness for NPV replication and host-range determination of NPVs.

Population structure of Cydia pomonella granulovirus isolates revealed by quantitative analysis of genetic variation

Genetic diversity of viruses is driven by genomic mutations and selection through its host, resulting in differences in virulence as well as host responses. For baculoviruses, which are naturally occurring pathogens of insects and which are frequently sprayed on hundred thousands to millions of hectares as biocontrol agents of insect pests, the phenomenon of virus-host co-evolution is of particular scientific interest and economic importance because high virulence of baculovirus products is essential and emergence of host resistance needs to be avoided as much as possible. In the present study, the population structure of twenty isolates of the Cydia pomonella granulovirus (CpGV), including twelve isolates from different geographic origins and eight commercial formulations, were studied on the basis of next-generation sequencing data and by analyzing the distribution of single nucleotide polymorphisms (SNPs). An entirely consensus sequence-free quantitative SNP analysis was applied for the identification of 753 variant SNP sites being specific for single as well as groups of CpGV isolates. Based on the quantitative SNP analysis, homogenous, heterogenous as well as mixed isolates were identified and their proportions of genotypes were deciphered, revealing a high genetic diversity of CpGV isolates from around the world. Based on hierarchical clustering on principal components (HCPC), six distinct isolate/group clusters were identified, representing the proposed main phylogenetic lineages of CpGV but comprising full genome information from virus mixtures. The relative location of different isolates in HCPC reflected the proportion of variable compositions of different genotypes. The established methods provide novel analysis tools to decipher the molecular complexity of genotype mixtures in baculovirus isolates, thus depicting the population structure of baculovirus isolates in a more adequate form than consensus based analyses.

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.

Comparative Analysis of HaSNPV-AC53 and Derived Strains

Complete genome sequences of two Australian isolates of H. armigera single nucleopolyhedrovirus (HaSNPV) and nine strains isolated by plaque selection in tissue culture identified multiple polymorphisms in tissue culture-derived strains compared to the consensus sequence of the parent isolate. Nine open reading frames (ORFs) in all tissue culture-derived strains contained changes in nucleotide sequences that resulted in changes in predicted amino acid sequence compared to the parent isolate. Of these, changes in predicted amino acid sequence of six ORFs were identical in all nine derived strains. Comparison of sequences and maximum likelihood estimation (MLE) of specific ORFs and whole genome sequences were used to compare the isolates and derived strains to published sequence data from other HaSNPV isolates. The Australian isolates and derived strains had greater sequence similarity to New World SNPV isolates from H. zea than to Old World isolates from H. armigera, but with characteristics associated with both. Three distinct geographic clusters within HaSNPV genome sequences were identified: Australia/Americas, Europe/Africa/India, and China. Comparison of sequences and fragmentation of ORFs suggest that geographic movement and passage in vitro result in distinct patterns of baculovirus strain selection and evolution.

The complete sequence of the first Spodoptera frugiperda Betabaculovirus genome: a natural multiple recombinant virus

Spodoptera frugiperda (Lepidoptera: Noctuidae) is a major pest in maize crops in Colombia, and affects several regions in America. A granulovirus isolated from S. frugiperda (SfGV VG008) has potential as an enhancer of insecticidal activity of previously described nucleopolyhedrovirus from the same insect species (SfMNPV). The SfGV VG008 genome was sequenced and analyzed showing circular double stranded DNA of 140,913 bp encoding 146 putative ORFs that include 37 Baculoviridae core genes, 88 shared with betabaculoviruses, two shared only with betabaculoviruses from Noctuide insects, two shared with alphabaculoviruses, three copies of own genes (paralogs) and the other 14 corresponding to unique genes without representation in the other baculovirus species. Particularly, the genome encodes for important virulence factors such as 4 chitinases and 2 enhancins. The sequence analysis revealed the existence of eight homologous regions (hrs) and also suggests processes of gene acquisition by horizontal transfer including the SfGV VG008 ORFs 046/047 (paralogs), 059, 089 and 099. The bioinformatics evidence indicates that the genome donors of mentioned genes could be alpha- and/or betabaculovirus species. The previous reported ability of SfGV VG008 to naturally co-infect the same host with other virus show a possible mechanism to capture genes and thus improve its fitness.

Unraveling the entry mechanism of baculoviruses and its evolutionary implications

The entry of baculovirus budded virus into host cells is mediated by two distinct types of envelope fusion proteins (EFPs), GP64 and F protein. Phylogenetic analysis suggested that F proteins were ancestral baculovirus EFPs, whereas GP64 was acquired by progenitor group I alphabaculovirus more recently and may have stimulated the formation of the group I lineage. This study was designed to experimentally recapitulate a possible major step in the evolution of baculoviruses. We demonstrated that the infectivity of an F-null group II alphabaculovirus (Helicoverpa armigera nucleopolyhedrovirus [HearNPV]) can be functionally rescued by coinsertion of GP64 along with the nonfusogenic F(def) (furin site mutated HaF) from HearNPV. Interestingly, HearNPV enters cells by endocytosis and, less efficiently, by direct membrane fusion at low pH. However, this recombinant HearNPV coexpressing F(def) and GP64 mimicked group I virus not only in its EFP composition but also in its abilities to enter host cells via low-pH-triggered direct fusion pathway. Neutralization assays indicated that the nonfusogenic F proteins contribute mainly to binding to susceptible cells, while GP64 contributes to fusion. Coinsertion of GP64 with an F-like protein (Ac23) from group I virus led to efficient rescue of an F-null group II virus. In summary, these recombinant viruses and their entry modes are considered to resemble an evolutionary event of the acquisition of GP64 by an ancestral group I virus and subsequent adaptive inactivation of the original F protein. The study described here provides the first experimental evidence to support the hypothesis of the evolution of baculovirus EFPs.

Baculovirus replication induces the expression of heat shock proteins in vivo and in vitro

A recent handful of studies have linked baculovirus infection with the induction of heat shock proteins, a highly conserved family of cytoprotective proteins. Here, we demonstrate baculovirus-stimulated upregulation of hsp70 transcription in the natural host, Helicoverpa zea. Larvae lethally infected with Helicoverpa zea single nucleopolyhedrovirus (HzSNPV) accumulated hsp70 transcripts throughout the 72-hour course of infection in the midgut, hemocytes, and fat body. While a maximal 17- or 15-fold induction of hsp70 was noted in the midgut and hemocytes, respectively, by 72 hours postinfection, the level of hsp70 transcription in the fat body of larvae was greater than two orders of magnitude higher than in mock-infected larvae. These results were largely mirrored in cultures of infected cells, and a potentiation effect was observed in cells that were both heat shocked and infected. In contrast, Spodoptera frugiperda multiple nucleopolyhedrovirus and ultraviolet-inactivated HzSNPV did not stimulate hsp70 transcription in these non-permissive larvae and in cell culture, respectively. Taken together, this report documents baculovirus-mediated upregulation of hsp70 in the host and demonstrates the requirement for productive infection for hsp70 induction in vitro and in vivo.

Molecular and immunohistochemical characterization of granulin gene encoded in Pieris rapae granulovirus genome

Pieris rapae granulovirus (PiraGV) is highly pathogenic to the cabbage butterfly (P. rapae), an important pest of cultivated cabbages and mustard crops. It therefore holds significant promise towards exploitation as a potent bio-control agent in the field controlling the pest population. Whole-genome elucidation of the Korean isolate of the granulovirus (PiraGV-K), reported the presence of a granulin gene corresponding to ORF 1 in its genome. Comprehensive studies towards functional characterization of the gene, established that it is composed of 744 nucleotides and encodes a peptide of 247 amino acid residues. It possessed significant homology with AoGV and ClanGV with 87% identity at amino acid level. Multiple alignment data suggests that the C-terminus region of the gene had three different conserved regions. Time-course studies conducted in PiraGV-K infected P. rapae larvae revealed a significant upsurge of the transcript (134-fold) at 4 days post infection followed by a significant decline at the most advanced stages of infection. Anti-PiraGV-K granulin antibody was produced and western blot conducted with the infected larvae further confirmed the induction pattern with a protein of 30 kDa. Immunofluorescent staining showed a granulin-specific signal in fat body and integument of the infected larvae. Granulin-specific signals were noticed 2 days post infection with the eventual systemic spread of infection to the associated tracheal matrix witnessed at 4 days post infection. Immunogold labeling and electron microscopic studies further proved the cytopathological effects as the presence of numerous membrane-bound vesicles with nucleocapsids and abruption of intercellular junctions in fat body and hypertrophied cells in the integument.

A baculovirus photolyase with DNA repair activity and circadian clock regulatory function

Cryptochromes and photolyases belong to the same family of flavoproteins but, despite being structurally conserved, display distinct functions. Photolyases use visible light to repair ultraviolet-induced DNA damage. Cryptochromes, however, function as blue-light receptors, circadian photoreceptors, or repressors of the CLOCK/BMAL1 heterodimer, the transcription activator controlling the molecular circadian clock. Here, we present evidence that the functional divergence between cryptochromes and photolyases is not so univocal. Chrysodeixis chalcites nucleopolyhedrovirus possesses 2 photolyase-like genes: phr1 and phr2. We show that PHR1 and PHR2 are able to bind the CLOCK protein. Only for PHR2, however, the physical interaction with CLOCK represses CLOCK/BMAL1-driven transcription. This result shows that binding of photolyase per se is not sufficient to inhibit the CLOCK/BMAL1 heterodimer. PHR2, furthermore, affects the oscillation of immortalized mouse embryonic fibroblasts, suggesting that PHR2 can regulate the molecular circadian clock. These findings are relevant for further understanding the evolution of cryptochromes and photolyases as well as behavioral changes induced in insects by baculoviruses.

Identification of nucleopolyhedrovirus that infect Nymphalid butterflies Agraulis vanillae and Dione juno

Dione juno and Agraulis vanillae are very common butterflies in natural gardens in South America, and also bred worldwide. In addition, larvae of these butterflies are considered as pests in crops of Passiflora spp. For these reasons, it is important to identify and describe pathogens of these species, both for preservation purposes and for use in pest control. Baculoviridae is a family of insect viruses that predominantly infect species of Lepidoptera and are used as bioinsecticides. Larvae of D. juno and A. vanillae exhibiting symptoms of baculovirus infection were examined for the presence of baculoviruses by PCR and transmission electron microscopy. Degenerate primers were designed and used to amplify partial sequences from the baculovirus p74, cathepsin, and chitinase genes, along with previously designed primers for amplification of lef-8, lef-9, and polh. Sequence data from these six loci, along with ultrastructural observations on occlusion bodies isolated from the larvae, confirmed that the larvae were infected with nucleopolyhedroviruses from genus Alphabaculovirus. The NPVs from the two different larval hosts appear to be variants of the same, previously undescribed baculovirus species. Phylogenetic analysis of the sequence data placed these NPVs in Alphabaculovirus group I/clade 1b.