Mass identification of transcriptional units expressed from the Bombyx mori nucleopolyhedrovirus genome

Loss of p24 from the Bombyx mori nucleopolyhedrovirus genome results in the formation of cuboidal occlusion bodies

Some insect viruses produce the occlusion body (OB), a large crystalline particle comprising a viral protein that occludes virions to protect them from harsh environments. The shapes and sizes of OBs are diverse depending on baculovirus species, but the detailed molecular mechanism determining them has yet to be totally clarified yet. Here we generated Bombyx mori nucleopolyhedrovirus (BmNPV) mutants of the p24 gene that encodes a viral capsid protein and found that p24-mutated BmNPVs produced cuboidal OBs with a slightly larger size than typical truncated octahedral OBs produced by wild-type BmNPVs. Meanwhile, p24 disruption has no significant impact on progeny virus production and viral pathogenicity. In addition, we experimentally demonstrated that a single amino acid substitution found in the P24 protein of the BmNPV Cubic isolate caused cuboidal OB production. These results suggest that p24 has a crucial role in generating the typical shape of OBs.

Bombyx mori nucleopolyhedrovirus Bm96 suppresses viral virulence in Bombyx mori larvae

Bombyx mori nucleopolyhedrovirus (BmNPV) is a severe pathogen for the domestic silkworm, Bombyx mori. BmNPV harbors over 140 protein-coding genes in its 128.4 kilobase pair-long double-stranded genome. However, many BmNPV genes are still uncharacterized. Here we investigated the role of BmNPV Bm96 in both B. mori cultured cells and larvae. We found that Bm96 is mainly expressed at the late stage of infection and accumulation of Bm96 protein peaks at 24 h post infection (hpi) and declines gradually at 48 hpi in B. mori cultured cells. Compared with the wild-type viruses, Bm96-deletion viruses exhibited higher viral propagation and fast-killing phenotype in B. mori larvae. These results strongly suggest that Bm96 negatively regulates the propagation of BmNPV in B. mori larvae. Furthermore, we observed that larvae infected with Bm96-deletion viruses showed lower locomotory activity at the late stage of infection compared with those infected with the wild-type viruses.

Bombyx mori nucleopolyhedrovirus protein Bm11 is involved in occlusion body production and occlusion-derived virus embedding

Bombyx mori nucleopolyhedrovirus (BmNPV) orf11 (bm11) is a highly conserved gene with unknown function. It is homologous to AcMNPV orf19. In this study, a bm11 knockout virus was constructed and its role was investigated. Expression analysis indicated that bm11 is a late gene and confocal microscopy analysis demonstrated that Bm11 localizes predominantly in the nuclear ring zone at the late phase of infection. The bm11 deletion did not affect budded virus (BV) production or viral genome replication, but markedly reduced the production of occlusion bodies (OBs) and the embedding of occlusion-derived viruses (ODVs). Bio-assays showed that Bm11 was involved in BmNPV infectivity in vivo by direct injection. In conclusion, our results demonstrated that although Bm11 is not essential for BV production or mature ODV formation, it affects OB production and ODV occlusion.

MacoNPV baculovirus midgut-specific gene expression during infection of the bertha armyworm, Mamestra configurata

Baculoviruses have two forms, occlusion derived virus (ODV) which is responsible for primary infection in host midgut tissue and budded virus (BV), which infects all other host tissues during secondary infection. This study examined the primary infection by ODV of midgut cells of bertha armyworm Mamestra configurata fourth instar larvae and measured the expression of viral genes over a time course of infection. Both digital PCR and RNA sequencing methods showed the profile of transcription to be different from those produced by AcMNPV BV infection of in vitro cell cultures. This included having unique collections of genes expressed early, as well as much greater late gene expression of p6.9 and much reduced expression of polh and p10. These differences likely reflect characteristics unique to the critical step of in vivo midgut cell infection, and provide insights into the processes that regulate viral gene expression in different host tissues.

Bombyx mori nucleopolyhedrovirus BM5 protein regulates progeny virus production and viral gene expression

Bombyx mori nucleopolyhedrovirus (BmNPV) orf5 (Bm5) is a core gene of lepidopteran baculoviruses and encodes the protein with the conserved amino acid residues (DUF3627) in its C-terminus. Here, we found that Bm5 disruption resulted in lower titers of budded viruses and fewer numbers of occlusion bodies (OBs) in B. mori cultured cells and larvae, although viral genome replication was not affected. Bm5 disruption also caused aberrant expression of various viral genes at the very late stage of infection. Immunocytochemical analysis revealed that BM5 localized to the nuclear membrane. We also found that DUF3627 is important for OB production, transcriptional regulation of viral genes, and subcellular localization of BM5. Compared with wild-type BmNPV infection, larval death was delayed when B. mori larvae were infected with Bm5 mutants. These results suggest that BM5 is involved in progeny virus production and regulation of viral gene expression at the very late stage of infection.

The killing speed of egt-inactivated Bombyx mori nucleopolyhedrovirus depends on the developmental stage of B. mori larvae

Several lines of evidence have shown that the deletion of the ecdysteroid UDP-glucosyltransferase gene (egt) from the nucleopolyhedrovirus (NPV) genome increases the killing speed of host lepidopteran larvae. However, it has not been investigated in detail whether the effects of egt deletion depend on the larval stages of the host insect. In this study, we performed bioassays using 10 continuous larval stages of the 4th- or 5th-instar Bombyx mori larvae and B. mori NPV egt mutants. The fast-killing phenotype was observed in the egt mutants only when the infection process progressed through larval-larval transition. All day-2 4th-instar larvae infected with the egt mutants entered the molting stage and died much earlier than wild-type-infected larvae. Bodies of egt mutant-infected larvae were filled with excessive fluid immediately after head capsule slippage, owing presumably to the degeneration of Malpighian tubules. Fourth- or 5th-instar larvae infected with the egt mutants at early stages of each instar died similarly to those infected with the wild-type virus. Under infection in the middle stages of the 5th-instar, the survival time of egt mutant-infected larvae was significantly longer than that of the wild-type virus-infected larvae. These results clearly show that the effects of egt deletion on killing speed of NPV are largely dependent on the developmental stage of the host larvae infected by the virus.

Modularity and evolutionary constraints in a baculovirus gene regulatory network

BACKGROUND

The structure of regulatory networks remains an open question in our understanding of complex biological systems. Interactions during complete viral life cycles present unique opportunities to understand how host-parasite network take shape and behave. The Anticarsia gemmatalis multiple nucleopolyhedrovirus (AgMNPV) is a large double-stranded DNA virus, whose genome may encode for 152 open reading frames (ORFs). Here we present the analysis of the ordered cascade of the AgMNPV gene expression.

RESULTS

We observed an earlier onset of the expression than previously reported for other baculoviruses, especially for genes involved in DNA replication. Most ORFs were expressed at higher levels in a more permissive host cell line. Genes with more than one copy in the genome had distinct expression profiles, which could indicate the acquisition of new functionalities. The transcription gene regulatory network (GRN) for 149 ORFs had a modular topology comprising five communities of highly interconnected nodes that separated key genes that are functionally related on different communities, possibly maximizing redundancy and GRN robustness by compartmentalization of important functions. Core conserved functions showed expression synchronicity, distinct GRN features and significantly less genetic diversity, consistent with evolutionary constraints imposed in key elements of biological systems. This reduced genetic diversity also had a positive correlation with the importance of the gene in our estimated GRN, supporting a relationship between phylogenetic data of baculovirus genes and network features inferred from expression data. We also observed that gene arrangement in overlapping transcripts was conserved among related baculoviruses, suggesting a principle of genome organization.

CONCLUSIONS

Albeit with a reduced number of nodes (149), the AgMNPV GRN had a topology and key characteristics similar to those observed in complex cellular organisms, which indicates that modularity may be a general feature of biological gene regulatory networks.

Functional characterization of Bombyx mori nucleopolyhedrovirus CG30 protein

The baculovirus cg30 gene is present in the genomes of most alphabaculoviruses. The gene product, CG30, contains two putative functional domains, a RING finger motif and a leucine zipper motif. A gene-knockout study in Autographa californica nucleopolyhedrovirus (AcMNPV) revealed that a cg30-disrupted virus did not show any striking differences compared with wild-type virus. To determine the roles of cg30 in another alphabaculovirus, we constructed two Bombyx mori NPV (BmNPV) mutants lacking a functional cg30 by lacZ cassette insertion and characterized its infectivity in BmN cells and B. mori larvae. The mutants produced fewer budded viruses (BVs) in BmN cells and B. mori larvae compared with wild-type BmNPV. We also observed a decrease in the release of occlusion bodies (OBs) in the hemolymph of the larvae infected with the cg30 mutants. To investigate the role of a RING finger domain of CG30 during virus growth, we further constructed two mutants; one expressed a mutant CG30 where a RING finger motif is disrupted by a single amino acid substitution, whereas the other possessed a CG30 derivative that completely lacked the RING finger domain. Both produced fewer OBs in the hemolymph of B. mori larvae, demonstrating that a RING finger domain of CG30 is involved in maximum OB production in the hemolymph of B. mori larvae. We also revealed that CG30 is expressed as a nuclear protein with a molecular mass of 30kDa in BmNPV-infected cells.

Identification of Bombyx mori bidensovirus VD1-ORF4 reveals a novel protein associated with viral structural component

Bombyx mori bidensovirus (BmBDV) VD1-ORF4 (open reading frame 4, ORF4) consists of 3,318 nucleotides, which codes for a predicted 1,105-amino acid protein containing a conserved DNA polymerase motif. However, its functions in viral propagation remain unknown. In the current study, the transcription of VD1-ORF4 was examined from 6 to 96 h postinfection (p.i.) by RT-PCR, 5'-RACE revealed the transcription initiation site of BmBDV ORF4 to be -16 nucleotides upstream from the start codon, and 3'-RACE revealed the transcription termination site of VD1-ORF4 to be +7 nucleotides downstream from termination codon. Three different proteins were examined in the extracts of BmBDV-infected silkworms midguts by Western blot using raised antibodies against VD1-ORF4 deduced amino acid, and a specific protein band about 53 kDa was further detected in purified virions using the same antibodies. Taken together, BmBDV VD1-ORF4 codes for three or more proteins during the viral life cycle, one of which is a 53 kDa protein and confirmed to be a component of BmBDV virion.

Characterization of the open reading frame 7a from Bombyx mori nucleopolyhedrovirus

The open reading frame 7a of the Bombyx mori nucleopolyhedrovirus (BmNPV orf7a, Bm7a) encodes a predicted polypeptide consisting of 53 amino acid residues. While this polypeptide's homologues were found in certain baculoviruses, its function has not yet been studied. A protein structural assay revealed that the polypeptide expressed from this gene contained a predicted signal peptide at its N-terminus. To investigate the expression of Bm7a, 3'rapid amplification of cDNA ends was used to detect its transcript. A series of recombinant viruses were also constructed to analyze the cellular localization of the Bm7a products. The results indicated that the Bm7a gene is co-transcribed with bv/odv-e26 and orf9. The function of BM7a was studied by using fused enhanced green fluorescence protein (eGFP) as a visual marker to trace the localization of the protein during infection. Fluorescence was observed in the cellular membrane of the infected cells by confocal microscopy. Western blot analysis using an eGFP-antibody also showed that the BM7a-eGFP fusion protein was present in the budded virions (BVs). These data implied that BM7a is a component of the BV. In this work, we identified the function of the Bm7a gene and suggested that it may play a role in BV packaging.

Phenotypic grouping of 141 BmNPVs lacking viral gene sequences

We constructed a series of gene knockout BmNPVs (KOVs) for each of 141 genes (Gomi et al., 1999; Katsuma et al., 2011) using the BmNPV T3 bacmid system (Ono et al., 2007) and lambda red recombination system (Datsenko and Wanner, 2000). In a subsequent analysis of the properties needed for infection using a marker gene, egfp (enhanced green fluorescent protein gene), inserted into the polyhedrin locus, the knockout viruses (KOVs) were subdivided into four phenotypic types, A to D. Type-A (86 KOVs) showed the ability to expand infections equivalent to the control while type-B (8 KOVs) spread infections more slowly. Type-C (37 KOVs) expressed egfp in transfected-BmN cells but the production of infectious viruses was not observed. Type-D (10 KOVs) showed no ability to express egfp even in the transfection experiments. KOVs lacking genes (pkip (Bm15), gp41 (Bm66), bro-d (Bm131), Bm20, 48, 65, 91, 93, or 101) previously identified as being essential, were placed in the viable type-A and B categories.

Baculovirus-encoded protein BV/ODV-E26 determines tissue tropism and virulence in lepidopteran insects

Lepidopteran nucleopolyhedroviruses (NPVs) show distinct tissue tropism in host insect larvae. However, the molecular mechanism of this tropism is largely unknown. We quantitatively investigated NPV tissue tropism by measuring mRNA levels of viral genes in 16 tissues from Bombyx mori NPV (BmNPV)-infected B. mori larvae and found clear tissue tropism, i.e., BmNPV replicates poorly in the silk glands, midgut, and Malpighian tubule compared with other larval tissues. We next identified the viral genes determining tissue tropism in NPV infection by investigating the phenotypes of larvae infected with 44 BmNPV mutants in which one gene was functionally disrupted by a LacZ cassette insertion. We found that occlusion body (OB) production was markedly enhanced compared with that of the wild type in the middle silk glands (MSGs) of larvae infected with three mutants in which one of three tandemly arrayed genes (Bm7, Bm8, and Bm9) was disrupted. We generated additional mutants in which one or two genes of this gene cluster were partially deleted and showed that Bm8, also known as BV/ODV-E26, was solely required for the suppression of OB production in the MSGs of BmNPV-infected B. mori larvae. Western blotting showed that a LacZ cassette insertion in Bm7 or Bm9 resulted in aberrant expression of Bm8, presumably leading to abnormal OB production in the MSGs. Larval bioassays also revealed that disruption of Bm8 accelerated the death of B. mori larvae. These results suggest that the group I NPV-specific protein BV/ODV-E26 determines tissue tropism and virulence in host lepidopteran insects.

Transcriptomic profiling of Microplitis demolitor bracovirus reveals host, tissue and stage-specific patterns of activity

The polydnaviruses (PDVs) are a family of DNA viruses that are symbiotically associated with parasitoid wasps. The transcription of particular genes or gene-family members have been reported for several PDVs, but no studies have characterized the spatio-temporal patterns of expression for the entire complement of predicted genes in the encapsidated genome of any PDV isolate. The braconid wasp Microplitis demolitor carries the PDV Microplitis demolitor bracovirus (MdBV) and parasitizes larval stage Pseudoplusia (Chrysodeixis) includens. The encapsidated genome consists of 15 genomic segments with 51 predicted ORFs encoding proteins ≥100 aa. A majority of these ORFs form four multimember gene families (ptp, ank, glc and egf) while the remaining ORFs consist of single copy (orph) genes. Here we used RT-PCR and quantitative real-time PCR methods to profile the encapsidated transcriptome of MdBV in P. includens and M. demolitor. Our results indicate that most predicted genes are expressed in P. includens. Spatial patterns of expression in P. includens differed among genes, but temporal patterns of expression were generally similar, with transcript abundance progressively declining between 24 and 120 h. A subset of ptp, ank and orph genes were also expressed in adult female but not male M. demolitor. Only one encapsidated gene (ank-H4) was expressed in all life stages of M. demolitor, albeit at much lower levels than in P. includens. However, another encapsidated gene (orph-B1) was expressed in adult M. demolitor at similar levels to those detected in P. includens.