Analysis of the genetic information of a DNA segment of a new virus from silkworm

Transcriptome analysis of the immune response of silkworm at the early stage of Bombyx mori bidensovirus infection

Bombyx mori bidensovirus (BmBDV) infects silkworm midgut and causes chronic flacherie disease; however, the interaction between BmBDV and silkworm is unclear. Twenty-four hours after BmBDV infection, the midgut was extracted for RNA-seq to analyze the factors associated with BmBDV-invasion and the early antiviral immune response in silkworms. The total reads from each sample were more than 16100000 and the number of expressed genes exceeded 8200. There were 334 upregulated and 272 downregulated differentially expressed genes (DEGs). Gene ontology analysis of DEGs showed that structural constituents of cuticle, antioxidant, and immune system processes were upregulated. Further analysis revealed BmBDV-mediated induction of BmorCPR23 and BmorCPR44, suggesting possible involvement in viral invasion. Antioxidant genes that protect host cells from virus-induced oxidative stress, were significantly upregulated after BmBDV infection. Several genes related to peroxisomes, apoptosis, and autophagy-which may be involved in antiviral immunity-were induced by BmBDV. These results provide insights into the mechanism of BmBDV infection and host defense.

Characterization of the RNA Transcription Profile of Bombyx mori Bidensovirus

Bombyx mori bidensovirus (BmBDV) is a single-stranded DNA (ssDNA) virus from the genus Bidensovirus of the Bidnaviridae family, which, thus far, solely infects insects. It has a unique genome that contains bipartite DNA molecules (VD1 and VD2). In this study, we explored the detailed transcription mapping of the complete BmBDV genome (VD1 and VD2) by rapid amplification of cDNA ends (RACE), reverse transcription quantitative real-time PCR (RT-qPCR), and luciferase assays. For the first time, we report the transcription map of VD2. Our mapping of the transcriptional start sites reveals that the NS genes in VD1 have separate transcripts that are derived from overlapping promoters, P5 and P5.5. Thus, our study provides a strategy for alternative promoter usage in the expression of BmBDV genes.

Increased antiviral capacity of transgenic silkworm via knockdown of multiple genes on Bombyx mori bidensovirus

Bombyx mori bidensovirus (BmBDV) causes fatal flacherie disease leading to severe economic losses in sericultures. The BmDNV-Z genome contains two single-stranded DNA molecules, VD1 and VD2. For generating silkworm lines with antiviral properties, two transgenic RNA interference (RNAi) vectors were constructed. Open reading frames (ORFs) 1-4 of VD1 were knockdown by vector pb-BDV1 while ORF1a, ORF1b, and ORF3 of VD2 were knockdown by vector pb-BDV2. Transgenic silkworm lines BDV1-I and BDV2-I were generated via RNAi microinjection. Mortality rates of BDV1-I and BDV2-I were reduced by 45% and 39%, respectively, and quantitative PCR showed that VD1 and VD2 contents in BDV1-I and BDV2-I were significantly lower than in the non-transgenic line. However, economic traits showed no obvious differences. Thus, knockdown of multiple BmDNV-Z genes provides strong resistance to BDV1-I and BDV2-I lines, and these can be used in sericulture without hampering silk production.

Transcriptome analysis of interactions between silkworm and cytoplasmic polyhedrosis virus

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) specifically infects silkworm midgut (MG) and multiplication occurs mainly in posterior midgut (PM). In this study, MG and fat body (FB) were extracted at 0, 3, 24, and 72 h after BmCPV infection. The total sequence reads of each sample were more than 1510000, and the mapping ratio exceeded 95.3%. Upregulated transcripts increased in MG during the infection process. Gene ontology (GO) categories showed that antioxidants were all upregulated in FB but not in MG. BGI001299, BGI014434, BGI012068, and BGI009201 were MG-specific genes with transmembrane transport function, the expression of which were induced by BmCPV. BGI001299, BGI014434, and BGI012068 expressed in entire MG and may be involved in BmCPV invasion. BGI009201 expressed only in PM and may be necessary for BmCPV proliferation. BmPGRP-S2 and BGI012452 (a putative serine protease) were induced by BmCPV and may be involved in immune defense against BmCPV. The expression level of BmCPV S1, S2, S3, S6, and S7 was high and there was no expression of S9 in MG 72 h, implying that the expression time of structural protein coding genes is earlier. These results provide insights into the mechanism of BmCPV infection and host defense.

The progress and future of enhancing antiviral capacity by transgenic technology in the silkworm Bombyx mori

Bombyx mori is a common lepidopteran model and an important economic insect for silk production. B. mori nucleopolyhedrovirus (BmNPV) is a typical pathogenic baculovirus that causes serious economic losses in sericulture. B. mori and BmNPV are a model of insect host and pathogen interaction including invasion of the host by the pathogen, host response, and enhancement of host resistance. The antiviral capacity of silkworms can be improved by transgenic technology such as overexpression of an endogenous or exogenous antiviral gene, RNA interference of the BmNPV gene, or regulation of the immune pathway to inhibit BmNPV at different stages of infection. Antiviral capacity could be further increased by combining different methods. We discuss the future of an antiviral strategy in silkworm, including possible improvement of anti-BmNPV, the feasibility of constructing transgenic silkworms with resistance to multiple viruses, and the safety of transgenic silkworms. The silkworm model could provide a reference for disease control in other organisms.

Direct interactions between bidensovirus BmDNV-Z proteins and midgut proteins from the virus target Bombyx mori

In this study we aimed to identify the protein-protein interactions between Bombyx mori midgut and the bidensovirus BmDNV-Z via a yeast two-hybrid (Y2H) system. To achieve this we constructed a Gal4 activation domain fusion library that expresses the host genes and Gal4 DNA binding domain fusion bait vectors that express BmDNV-Z genes. Y2H assay revealed 15 potential interactions between host and viral proteins. To verify the interactions, we modified and reconstructed a pair of bimolecular fluorescence complementation (BiFC) vectors and achieved the co-expressions of the candidate host genes and viral genes in insect culture cells. The BiFC assay confirmed the specificity of the interactions including B. mori 35 kDa protease and two BmDNV-Z proteins encoded by VD1-ORF2 and VD2-ORF1; B. mori transgelin and BmDNV-Z protein encoded by VD2-ORF3; and B. mori serine protease precursor and BmDNV-Z encoded by VD2-ORF3 in vitro. Our findings revealed that the specific host midgut proteins are involved in the interactions between B. mori and BmDNV-Z, which will facilitate our understanding of the molecular mechanisms of BmDNV-Z infection.

The non-structural protein NS-2 of Bombyx mori parvo-like virus is localized to the nuclear membrane

Bombyx mori parvo-like virus (BmPLV) has two complementary single-stranded DNA genome (VD1 and VD2) and owns a self-encoding DNA polymerase motif, but its replication mechanism is unclear. In our previous research, a protein encoded by VD1-ORF1 was indentified in the midgut of BmPLV China Zhenjiang isolate-(BmPLV-Z) infected silkworm larvae via two-dimensional gel electrophoresis (2-DE). This protein was named as non-structural protein 2 (NS2), which showed no similarity to that of parvoviruses. To date, little is known about it. In this study, sequence alignment results showed that NS2 shared homology with some chromosomal replication initiator protein dnaA and DNA-binding response regulators. The ns2 was cloned and expressed in E. coli, and then a polyclonal antibody of the NS2 protein was prepared successfully. The data from real-time quantitative PCR displayed that the transcription of VD1-ORF1 from BmPLV-Z-infected midguts started from 28-h post inoculation (h p.i.) in low amounts, but in high amounts at late stages of infection. Immunofluorescence showed that NS2 ultimately concentrated on the nuclear membrane in BmN cells at late stages, indicating that NS2 might be associated with integral membrane protein.

Expression of non-structural protein NS3 gene of Bombyx mori densovirus (China isolate)

The invertebrate parvovirus Bombyx mori Densonucleosis Virus type 3 (China isolate), named BmDNV-3, is a kind of bidensovirus. It is a new type of virus with unique replication mechanisms. To investigate the effects of the NS3 gene during viral DNA replication, a pair of primers was designed for amplifying NS3 gene of Bombyx mori densovirus (China isolate). Gene NS3 amplified was cloned into a prokaryotic expression vector pET-30a and the donor plasmid pFastBacHTe, respectively. The NS3 protein was expressed in Escherichia coli BL21. The pFastBacHTe-NS3 was transformed to E. coli DH10Bac. The recombinant bacmid baculoviruses (rBacmid-EGFP-NS3) isolated from the white colonies were transfected into BmN-4 cells using a transfection reagent. BmN-4 cells were infected with recombinant virus to express fusion proteins. The expression of fusion protein around 30 kDa in E. coli BL21 was identified by SDS-PAGE, Western blotting, and mass spectrometry. The expressed NS3 protein by B. mori nucleopolyhedrovirus bacmid system was confirmed by Western blotting using an anti-NS3 polyclonal antibody. And about 45 kDa protein was found. The expressed fusion protein was smaller than the expected size of EGFP-NS3, 55 kDa. Western blotting analysis indicated that EGFP-NS3 protein was expressed in infected larvae with smaller molecular size.

Differentially expressed genes in resistant and susceptible Bombyx mori strains infected with a densonucleosis virus

We investigated variations in the gene expression of Bombyx mori following infection with a densonucleosis virus (BmDNV-Z). Two B. mori near-isogenic lines, Jingsong and Jingsong.nsd-Z.NIL, which are highly susceptible and completely resistant to BmDNV-Z, respectively, were used in this study. The infection profiles of BmDNV-Z in the midguts of the B. mori Jingsong and Jingsong.nsd-Z.NIL larvae revealed that the virus invaded the midguts of both of these strains. However, its proliferation was notably inhibited in the midgut of the resistant strain. By using the suppression subtractive hybridization method, three cDNA libraries were constructed to compare BmDNV-Z responsive gene expression between the two silkworm lines. In total, 151 differentially expressed genes were obtained. Real-time qPCR analysis confirmed that 11 genes were significantly up-regulated in the midgut of the Jingsong.nsd-Z.NIL strain following BmDNV-Z infection. Our results imply that these up-regulated genes might be involved in B. mori immune responses against BmDNV infection.

Characterization of the genome structure of Bombyx mori densovirus (China isolate)

The genome of Bombyx mori densovirus (China isolate), termed as BmDNV-3, is composed of two kinds of different single-stranded linear DNA molecules (VD1 and VD2). In this study, the viral DNA molecules were purified and cloned into pUC119 vector, and the complete nucleotide sequence was determined. Sequence analysis showed that VD1 genome consisted of 6,543 nts including inverted terminal repeats (ITRs) of 224 nts, and VD2 genome consisted of 6,022 nts including ITRs of 524 nts. Comparison of the complete genome sequence between BmDNV-3 and BmDNV-2 (Yamanashi isolate) showed an identity of 98.4% in VD1 and 97.7% in VD2, with a total number of 228 bp substitutions, 11 bp deletions and 3 bp insertions found in BmDNV-3. A single nucleotide "A" deletion at nt 1589 in BmDNV-3 caused a frame shift mutation and brought about a premature stop codon, thus dividing VD2 of BmDNV-3 into two ORFs (named VD2 ORF1a and VD2 ORF1b) within that region, while there was only one ORF (named VD2 ORF1) in the corresponding region of BmDNV-2 (Yamanashi isolate). Comparative polymorphisms of ORFs and ITR regions of the two viral genomes showed that highly variable regions were mainly located in VD1 ORF3, VD1 ORF4, VD2 ORF2, and ITRs of BmDNV-3. Northern blots analysis revealed that VD1 had 1.1 kb and 1.5 kb transcripts from the left half of its plus strand, and one transcript about 3.3 kb from the right half of its minus strand. Sequencing of 3' and 5' RACE products showed that the 1.1 kb transcript started at nt 290 and ended at nt 1437, the 1.5 kb transcript started at nt 1423 and ended at nt 2931, and the 3.3 kb transcript started at nt 6287 and ended at nt 2922. These results help us to further understand the variation between different DNV genera and its possible causes, providing clues for studying the evolutionary history of densoviruses.

Sequence analysis of the Choristoneura occidentalis granulovirus genome

The genome of the Choristoneura occidentalis granulovirus (ChocGV) isolated from the western spruce budworm, Choristoneura occidentalis, was sequenced completely. It was 104,710 bp long, with a 67.3% A+T content and contained 116 potential open reading frames (ORFs) covering 88.4% of the genome. Of these, 29 ORFs were conserved in all fully sequenced baculovirus genomes, 30 were GV-specific, 53 were present in some nucleopolyhedroviruses (NPVs) and/or GVs, three were common to ChocGV and Choristoneura fumiferana GV (ChfuGV) and one was so far unique. To date, ChocGV is the only GV identified that contains a homologue of the apoptosis inhibitor protein P35/P49, present in some group I NPVs. It is also the first GV without a Xestia c-nigrum GV ORF 26 homologue. Five homologous regions (hrs)/repeat regions, lacking typical NPV hr palindromes were identified. ChocGV hrs were similar to each other but not to other GV hrs. A 1.8 kb repeat region with a high A+T content (81%) and multiple repeats of 21-210 bp was found between choc36 and 37. This area resembled the non-homologous region origin of DNA replication (non-hr ori) identified in Cryptophlebia leucotreta GV (CrleGV) and Cydia pomonella GV (CpGV). Based on the mean amino acid identities of homologous proteins, ChocGV was closest to fully sequenced genomes CpGV (52.3%) and CrleGV (52.1%). The closest amino acid identity was to individual ORFs from the partially sequenced ChfuGV genome (97.2% in 38 ORFs). Phylogenetic analysis placed ChocGV in a clade with CrleGV and CpGV.

Expansion of family Reoviridae to include nine-segmented dsRNA viruses: isolation and characterization of a new virus designated Aedes pseudoscutellaris reovirus assigned to a proposed genus (Dinovernavirus)

Family Reoviridae is known, by definition, to contain dsRNA viruses with 10-12 genome segments. We report here the characterization of the first member of this family with a nine-segmented genome. This virus was isolated from Aedes pseudoscutellaris mosquito cells and designated aedes pseudoscutellaris reovirus (APRV). Virions are single-shelled with turrets but are non-occluded by contrast to cypoviruses. APRV replicates in various mosquito cell lines, but not in mice or mammalian cells. Complete sequence analysis showed that APRV is phylogenetically related to cypoviruses, fijiviruses and oryzaviruses. The maximum amino acid identities with cypoviruses, oryzaviruses or fijiviruses in the polymerase, are compatible with values observed between these genera and lower than values within a given genus. This suggests that APRV should be classified within a new genus that we designated Dinovernavirus (sigla from D: Double-stranded, i: insect, nove: nine from the latin "novem", rna: RNA, virus) in family Reoviridae.

The genome of the Cryptophlebia leucotreta granulovirus

The genome of the Cryptophlebia leucotreta granulovirus (CrleGV) was sequenced and analyzed. The double-stranded circular genome contains 110907 bp and potentially encodes 129 predicted open reading frames (ORFs), 124 of which were similar to other baculovirus ORFs. Five ORFs were CrleGV specific and 26 ORFs were common to other granulovirus genomes. One ORF showed a significant similarity to a nonstructural protein of Bombyx mori densovirus-2. A baculovirus chitinase gene was identified, which is most likely not functional, because its central coding region including the conserved chitinase active site signature is deleted. Three gene copies (Crle20, 23, and 24) containing the Baculo PEP N domain of the polyhedron envelope protein were identified in CrleGV and other GV genomes. One of them (Crle23) appeared also to contain a p10-like sequence encoding of a number of leucine-rich heptad repeats and a proline-rich domain. Another striking feature of the genome is the presence of a hypervariable non-hr ori-like region of about 1800 bp consisting of different kinds of repeats and palindromes. Three other repeat-rich regions were identified within the genome and are considered as homologous regions (hrs). CrleGV is most closely related to the Cydia pomonella granulovirus (CpGV) as revealed by genome order comparisons and phylogenetic analyses. However, the AT content of the CrleGV genome, which is 67.6% and the highest found so far in baculoviruses, differed by 12.8% from the AT content of CpGV. This resulted in a major difference in the codon usage of both viruses and may reflect adaptive selection constraints to their particular hosts.

Evolutionary relationships among parvoviruses: virus-host coevolution among autonomous primate parvoviruses and links between adeno-associated and avian parvoviruses

The current classification of parvoviruses is based on virus host range and helper virus dependence, while little data on evolutionary relationships among viruses are available. We identified and analyzed 472 sequences of parvoviruses, among which there were (virtually) full-length genomes of all 41 viruses currently recognized as individual species within the family Parvoviridae. Our phylogenetic analysis of full-length genomes as well as open reading frames distinguished three evolutionary groups of parvoviruses from vertebrates: (i) the human helper-dependent adeno-associated virus (AAV) serotypes 1 to 6 and the autonomous avian parvoviruses; (ii) the bovine, chipmunk, and autonomous primate parvoviruses, including human viruses B19 and V9; and (iii) the parvoviruses from rodents (except for chipmunks), carnivores, and pigs. Each of these three evolutionary groups could be further subdivided, reflecting both virus-host coevolution and multiple cross-species transmissions in the evolutionary history of parvoviruses. No parvoviruses from invertebrates clustered with vertebrate parvoviruses. Our analysis provided evidence for negative selection among parvoviruses, the independent evolution of their genes, and recombination among parvoviruses from rodents. The topology of the phylogenetic tree of autonomous human and simian parvoviruses matched exactly the topology of the primate family tree, as based on the analysis of primate mitochondrial DNA. Viruses belonging to the AAV group were not evolutionarily linked to other primate parvoviruses but were linked to the parvoviruses of birds. The two lineages of human parvoviruses may have resulted from independent ancient zoonotic infections. Our results provide an argument for reclassification of Parvovirinae based on evolutionary relationships among viruses.

Identification and genetic mapping of RAPD markers linked to the densonucleosis refractoriness gene, nsd-2, in the silkworm, Bombyx mori

In the silkworm, Bombyx mori, nonsusceptibility to B. mori densonucleosis virus type-2 (BmDNV-2) is controlled by a recessive gene, nsd-2 (nonsusceptibility to DNV-2). We investigated the genetic linkage between two random-amplified polymorphic DNA (RAPD) markers and the +nsd-2 gene. Initially, we constructed the JSD-2 strain (nsd-2/+), which is congenic to strain J137 (nsd-2/nsd-2) with respect to the chromosome containing the +nsd-2 gene, starting with a female of strain J137 and a male of strain C137 (+nsd-2/+nsd-2). Genomic DNAs were compared between infected individuals of the JSD-2 strain and J137 by a polymerase chain reaction (PCR) with 700 arbitrary 10-mer primers. Two RAPD markers (OPH19R and OPP01R) linked to the +nsd-2 gene were found. For the crossing-over experiment, a female of J137 was crossed with a male (nsd-2/+) of JSD-2. Segregation analysis showed that the most closely linked RAPD marker (OPP01R) was mapped 4.7 cM distant from +nsd-2.

Analysis of proteins encoded in the bipartite genome of a new type of parvo-like virus isolated from silkworm - structural protein with DNA polymerase motif

Bombyx mori densonucleosis virus type 2 (BmDNV-2) is a small, spherical virus containing two complementary single-stranded linear DNA molecules (VD1, VD2). BmDNV-2 is a new type of virus with a unique, yet unspecified replication mechanism which is different from that of parvoviruses (Bando, H., Choi, H., Ito, Y., Nakagaki, M. , Kawase, S., 1992. Structural analysis on the single-stranded genomic DNAs of the virus newly isolated from silkworm: the DNA molecules share a common terminal sequence, Arch. Virol. 124, 187-193; Bando, H., Hayakawa, T., Asano, S., Sahara, K., Nakagaki, M. , Iizuka, T., 1995. Analysis of the genetic information of a DNA segment of a new virus from silkworm, Arch. Virol., 140, 1147-1155; Hayakawa, T., Asano, S., Sahara, K., Iizuka, T., Bando, H., 1997. Detection of replicative intermediate with closed terminus of Bombyx densonucleosis virus. Arch. Virol. 142, 1-7). Recent analyses on the genomic information of BmDNV-2 identified open reading frames which code for three tentative nonstructural proteins and four (VP1 to 4) of the six known structural proteins (Bando, H., Hayakawa, T., Asano, S., Sahara, K., Nakagaki, M., Iizuka, T., 1995. Analysis of the genetic information of a DNA segment of a new virus from silkworm, Arch. Virol., 140, 1147-1155; Nakagaki et al., in preparation). In this report we demonstrate that the two largest ORFs, VD1-ORF1 and VD2-ORF1, code for the two remaining structural proteins. In addition, computer-assisted analysis revealed that the structural protein encoded in VD1-ORF1 contains sequences conserved among various DNA polymerases, and showed an evolutionary relationship with the DNA polymerases involved in protein-primed replication.