Removal of transposon target sites from the Autographa californica multiple nucleopolyhedrovirus fp25k gene delays, but does not prevent, accumulation of the few polyhedra phenotype

Genome Sequencing and Organization of Three Geographically Different Isolates of Nucleopolyhedrovirus from the Gypsy Moth Reveal Significant Genomic Differences

Background

The gypsy moth (Lymantria dispar L., Lepidoptera: Erebidae) is a worldwide pest of trees and forests. Lymantria dispar nucleopolyhedrovirus (LdMNPV) belongs to the Baculoviridae family and is an insect virus specific to gypsy moth larvae. In this study, we describe the complete genome sequences of three geographically diverse isolates, H2 (China), J2 (Japan), and T3 (Turkey), of Lymantria dispar multiple nucleopolyhedrovirus (LdMNPV).

Methods

The genomes of isolates H2, J2, and T3 were subjected to shotgun pyrosequencing using Roche 454 FLX and assembled using Roche GS De Novo Assembler. Comparative analysis of all isolates was performed using bioinformatics methods.

Results

The genomes of LdMNPV-H2, J2, and T3 were 164,746, 162,249, and 162,614 bp in size, had GC content of 57.25%, 57.30%, and 57.46%, and contained 162, 165, and 164 putative open reading frames (ORFs ≥ 150 nt), respectively. Comparison between the reference genome LdMNPV-5/6 (AF081810) and the genomes of LdMNPV-H2, J2, and T3 revealed differences in gene content. Compared with LdMNPV-5/6, ORF5, 6, 8, 10, 31, and 67 were absent in LdMNPV-H2, ORF5, 13, and 66 were absent in LdMNPV-J2, and ORF10, 13, 31, and 67 were absent in LdMNPV-T3. In addition, the gene encoding the mucin-like protein (ORF4) was split into two parts in isolates H2 and T3 and designated ORF4a and ORF4b. Phylogenetic analysis grouped isolates H2 and J2 in a different cluster than isolate T3, which is more closely related to the Turkish and Polish isolates. In addition, H2 was found to be closely related to a South Korean LdMNPV isolate.

Conclusion

This study provided a more detailed overview of the relationships between different geographic LdMNPV isolates. The results showed remarkable differences between groups at the genome level.

A detailed model and Monte Carlo simulation for predicting DIP genome length distribution in baculovirus infection of insect cells

Baculoviruses have enormous potential for use as biopesticides to control insect pest populations without the adverse environmental effects posed by the widespread use of chemical pesticides. However, continuous baculovirus production is susceptible to DNA mutation and the subsequent production of defective interfering particles (DIPs). The amount of DIPs produced and their genome length distribution are of great interest not only for baculoviruses but for many other DNA and RNA viruses. In this study, we elucidate this aspect of virus replication using baculovirus as an example system and both experimental and modeling studies. The existing mathematical models for the virus replication process consider DIPs as a lumped quantity and do not consider the genome length distribution of the DIPs. In this study, a detailed population balance model for the cell-virus culture is presented, which predicts the genome length distribution of the DIP population along with their relative proportion. The model is simulated using the kinetic Monte Carlo algorithm, and the results agree well with the experimental results. Using this model, a practical strategy to maintain the DIP fraction to near to its maximum and minimum limits has been demonstrated.

Statistical modeling of cell-to-cell variability in viral infection during passaging in suspension cell culture: Application in Monte-Carlo simulation

Packaging during the passaging of viruses in cell cultures yields various phenotypes and is regulated by viral protein expression in infected cells. Although such a packaging mechanism has a profound effect in controlling the virus yield, little is known about the underlying statistical models followed by virus packaging and protein expression among cells infected with the virus. A predictive framework combining identification of the probability density function (PDF) based on log-likelihood and using the PDF for Monte-Carlo simulations is developed. The Birnbaum-Saunders distribution was found to be consistent with all three-virus packaging levels, including nucleocapsids/occlusion-derived virus (ODV), ODVs/polyhedra, and polyhedra/cell for both wild-type and genetically modified AcMNPV. Next, it was demonstrated that PDF fitting could be used to compare two viruses having distinctly different genetic configurations. Finally, the identified PDF can be incorporated in RNA synthesis parameters for baculovirus infection to predict the cell-to-cell variability in protein expression using Monte-Carlo simulations. The proposed tool can be used for the estimation of uncertainty in the kinetic parameter and prediction of cell-to-cell variability for other biological systems.

Cell-to-Cell Variability in Protein Expression during Viral Infection: Monte-Carlo Simulation and Validation based on Confocal Imaging

One of the major challenges is to identify the statistical model underlying the heterogeneity in viral protein expression in single cells. In this endeavor, we propose a computational tool to address the cell-to-cell variability in protein expression by random variate generation following probability distributions. Here, we show that statistical modeling using the probability density function of various distribution offers considerable potential for providing stochastic inputs to Monte Carlo simulation. Specifically, we present the ranking between three distribution families including gamma, normal and Weibull distribution using a comparison of cumulative frequency obtained from experiment and simulation. The major contribution of the proposed simulation method is to identify the underlying statistical model in kinetic parameters that capture the variability in protein expression in single cells obtained through imaging using confocal microscopy.

Baculovirus FP25K Localization: Role of the Coiled-Coil Domain

Two types of viruses are produced during the baculovirus life cycle: budded virus (BV) and occlusion-derived virus (ODV). A particular baculovirus protein, FP25K, is involved in the switch from BV to ODV production. Previously, FP25K from the model alphabaculovirus Autographa californica multiple nucleopolyhedrovirus (AcMNPV) was shown to traffic ODV envelope proteins. However, FP25K localization and the domains involved are inconclusive. Here we used a quantitative approach to study FP25K subcellular localization during infection using an AcMNPV bacmid virus that produces a functional AcMNPV FP25K-green fluorescent protein (GFP) fusion protein. During cell infection, FP25K-GFP localized primarily to the cytoplasm, particularly amorphous structures, with a small fraction being localized in the nucleus. To investigate the sequences involved in FP25K localization, an alignment of baculovirus FP25K sequences revealed that the N-terminal putative coiled-coil domain is present in all alphabaculoviruses but absent in betabaculoviruses. Structural prediction indicated a strong relatedness of AcMNPV FP25K to long interspersed element 1 (LINE-1) open reading frame 1 protein (ORF1p), which contains an N-terminal coiled-coil domain responsible for cytoplasmic retention. Point mutations and deletions of this domain lead to a change in AcMNPV FP25K localization from cytoplasmic to nuclear. The coiled-coil and C-terminal deletion viruses increased BV production. Furthermore, a betabaculovirus FP25K protein lacking this N-terminal coiled-coil domain localized predominantly to the nucleus and exhibited increased BV production. These data suggest that the acquisition of this N-terminal coiled-coil domain in FP25K is important for the evolution of alphabaculoviruses. Moreover, with the divergence of preocclusion nuclear membrane breakdown in betabaculoviruses and membrane integrity in alphabaculoviruses, this domain represents an alphabaculovirus adaptation for nuclear trafficking of occlusion-associated proteins.

IMPORTANCE

Baculovirus infection produces two forms of viruses: BV and ODV. Manufacturing of ODV involves trafficking of envelope proteins to the inner nuclear membrane, mediated partly through the FP25K protein. Since FP25K is present in alpha-, beta-, and gammabaculoviruses, it is uncertain if this trafficking function is conserved. In this study, we looked at alpha- and betabaculovirus FP25K trafficking by its localization. Alphabaculovirus FP25K localized primarily to the cytoplasm, whereas betabaculovirus FP25K localized to the nucleus. We found that an N-terminal coiled-coil domain present in all alphabaculovirus FP25K proteins, but absent in betabaculovirus FP25K, was critical for alphabaculovirus FP25K cytoplasmic localization. We believe that this represents an evolutionary process that partly led to the gain of function of this N-terminal coiled-coil domain in alphabaculovirus FP25K to aid in nuclear trafficking of occlusion-associated proteins. Due to betabaculovirus breakdown of the nuclear membrane before occlusion, this function is not needed, and the domain was lost or never acquired.

Reduction of polyhedrin mRNA and protein expression levels in Sf9 and Hi5 cell lines, but not in Sf21 cells, infected with Autographa californica multiple nucleopolyhedrovirus fp25k mutants

During cell infection, the fp25k gene of baculoviruses frequently mutates, producing the few polyhedra (FP) per cell phenotype with reduced polyhedrin (polh) expression levels compared with wild-type baculoviruses. Here we report that the fp25k gene of the model baculovirus, Autographa californica multiple nucleopolyhedrovirus (AcMNPV), contains two hypermutable seven-adenine (A7) mononucleotide repeats (MNRs) that were mutated to A8 MNRs and a TTAA site that had host DNA insertions, producing fp25k mutants during Sf21 cell infection. The FP phenotype in Sf9 and Hi5 cells was more pronounced than in Sf21 cells. AcMNPV fp25k mutants produced similar levels of polyhedra or enhanced GFP, which were both under the control of the AcMNPV polh promoter for expression, in Sf21 cells but lower levels in Sf9 and Hi5 cells compared with AcMNPV with an intact fp25k gene. This correlated with the polh mRNA levels detected in each cell line. The majority of Sf21 cells infected with fp25 mutants showed high polh promoter-mediated GFP expression levels. Two cell lines subcloned from Sf21 cells that were infected with fp25k mutants showed different GFP expression levels. Furthermore, a small proportion of Hi5 cells infected with fp25k mutants showed higher production of polyhedra and GFP expression than the rest, and the latter was not correlated with increased m.o.i. Therefore, these data suggest that AcMNPV polh promoter-mediated gene expression activities differ in the three cell lines and are influenced by different cells within the cell line.

AcMNPV core gene ac109 is required for budded virion transport to the nucleus and for occlusion of viral progeny

The Autographa californica multiple nucleopolyhedrovirus (AcMNPV) ac109 core gene has been previously characterized as an essential late gene. Our results showed that budded virions could be detected in supernatants of infected Sf-9 cells, even when ac109 knockout viruses displayed a single-cell infection phenotype. Moreover, confocal microscopy analysis revealed that budded virions can enter the cytoplasm but are unable to enter the cell nucleus. This defect could be repaired by complementing ac109 in trans. In addition, polyhedra of normal size could be detected in Sf-9 nuclei infected with ac109 knockout viruses. However, electron microscopy demonstrated that these occlusion bodies were empty. Altogether, these results indicate that ac109 is required for infectivity of both phenotypes of virus.

Production of baculovirus defective interfering particles during serial passage is delayed by removing transposon target sites in fp25k

Accumulation of baculovirus defective interfering particle (DIP) and few polyhedra (FP) mutants is a major limitation to continuous large-scale baculovirus production in insect-cell culture. Although overcoming these mutations would result in a cheaper platform for producing baculovirus biopesticides, little is known regarding the mechanism of FP and DIP formation. This issue was addressed by comparing DIP production of wild-type (WT) Autographa californica multiple nucleopolyhedrovirus (AcMNPV) with that of a recombinant AcMNPV (denoted Ac-FPm) containing a modified fp25k gene with altered transposon insertion sites that prevented transposon-mediated production of the FP phenotype. In addition to a reduction in the incidence of the FP phenotype, DIP formation was delayed on passaging of Ac-FPm compared with WT AcMNPV. Specifically, the yield of DIP DNA in Ac-FPm was significantly lower than in WT AcMNPV up to passage 16, thereby demonstrating that modifying the transposon insertion sites increases the genomic stability of AcMNPV. A critical component of this investigation was the optimization of a systematic method based on the use of pulsed-field gel electrophoresis (PFGE) to characterize extracellular virus DNA. Specifically, PFGE was used to detect defective genomes, determine defective genome sizes and quantify the amount of defective genome within a heterogeneous genome population of passaged virus.