Complete sequence and organization of Antheraea pernyi nucleopolyhedrovirus, a dr-rich baculovirus
- Zuo-Ming Nie†1,
- Zhi-Fang Zhang†2,
- Dan Wang†1,
- Ping-An He1,
- Cai-Ying Jiang1,
- Li Song1,
- Fang Chen1,
- Jie Xu1,
- Ling Yang2,
- Lin-Lin Yu2,
- Jian Chen1,
- Zheng-Bing Lv1,
- Jing-Jing Lu1,
- Xiang-Fu Wu1 and
- Yao-Zhou Zhang1Email author
© Nie et al; licensee BioMed Central Ltd. 2007
Received: 30 January 2007
Accepted: 24 July 2007
Published: 24 July 2007
The completion and reporting of baculovirus genomes is extremely important as it advances our understanding of gene function and evolution. Due to the large number of viral genomes now sequenced it is very important that authors present significantly detailed analyses to advance the understanding of the viral genomes. However, there is no report of the Antheraea pernyi nucleopolyhedrovirus (AnpeNPV) genome.
The genome of AnpeNPV, which infects Chinese tussah silkworm (Antheraea pernyi), was sequenced and analyzed. The genome was 126,629 bp in size. The G+C content of the genome, 53.4%, was higher than that of most of the sequenced baculoviruses. 147 open reading frames (ORFs) that putatively encode proteins of 50 or more amino acid residues with minimal overlap were determined. Of the 147 ORFs, 143 appeared to be homologous to other baculovirus genes, and 4 were unique to AnpeNPV. Furthermore, there are still 29 and 33 conserved genes present in all baculoviruses and all lepidopteran baculoviruses respectively. In addition, the total number of genes common to all lepidopteran NPVs is sill 74, however the 74 genes are somewhat different from the 74 genes identified before because of some new sequenced NPVs. Only 6 genes were found exclusively in all lepidopteran NPVs and 12 genes were found exclusively in all Group I NPVs. AnpeNPV encodes v-trex(Anpe115, a 3' to 5' repair exonuclease), which was observed only in CfMNPV and CfDEFNPV in Group I NPVs. This gene potentially originated by horizontal gene transfer from an ancestral host. In addition, AnpeNPV encodes two conotoxin-like gene homologues (ctls), ctl1 and ctl2, which were observed only in HycuNPV, OpMNPV and LdMNPV. Unlike other baculoviruses, only 3 typical homologous regions (hr s) were identified containing 2~9 repeats of a 30 bp-long palindromic core. However, 24 perfect or imperfect direct repeats (dr s) with a high degree of AT content were found within the intergenic spacer regions that may function as non-hr, ori-like regions found in GrleGV, CpGV and AdorGV. 9 dr s were also found in intragenic spacer regions of AnpeNPV.
AnpeNPV belongs to Group I NPVs and is most similar to HycuNPV, EppoNPV, OpMNPV and CfMNPV based on gene content, genome arrangement, and amino acid identity. In addition, analysis of genes that flank hr s supported the argument that these regions are involved in the transfer of sequences between the virus and host.
The family Baculoviridae consists of viruses that contain circular DNA genomes ranging in size from approximately 80 to 180 kbp. This family of viruses has only been shown to be pathogenic to arthropods particularly insects [1, 2] and have a very variable G+C content that ranges from 28% to 59 % . The virions have a complex construction and consist of an envelope and a nucleocapsid. Baculoviruses are divided into nucleopolyhedrovirus (NPVs) and granulovirus (GVs) generally. Based on phylogenetic analysis, baculoviruses can be classified into five major groups including GVs, Group I and Group II NPVs, a group of the dipteran viruses, Culex nigripalpus NPV (CuniNPV) and a group of the hymenopteran viruses . NPVs are mainly found in Lepidoptera and other insects such as Hymenoptera, Diptera, Coleoptera, Thysanura and Trichoptera. The virions contain either single or multiple nucleocapsids. All hymenopteran NPVs virions contain single nucleocapsid. GVs have been found only in Lepidoptera [6, 7].
Previous reports demonstrated that there are 29 conserved genes present in all baculoviruses and 33 conserved genes present in all lepidopteran baculoviruses [7–9], and a total of 74 genes are present in all lepidopteran NPVs. To date, the complete genomes of 31 NPVs and 8 GVs are published or available in GenBank. 31 NPV genomes have been published from Autographa californica NPV (AcMNPV)  to Neodiprion abietis NPV (NeabNPV) . 8 GV genomes published include Xestia c-nigrum GV (XecnGV), Plutella xylostella GV (PlxyGV), Cydia pomonella GV (CpGV), Adoxophyes orana GV(AdorGV), Cryptophlebia leucotreta GV (CrleGV) and Choristoneura occidentalis GV(ChocGV).
The Chinese tussah silkworm (also known as Chinese oak silkworm), Antheraea pernyi (Lipidoptera, Saturniidae), is an insect that can spin silk cocoons and eats the leaves of oak trees. As an important economic insect, Antheraea pernyi is commercially cultivated mainly in Middle and Northeastern China. It is also used as a source of food and for cosmetics. Antheraea pernyi NPV (AnpeNPV) can infect Antheraea pernyi resulting in nuclear polyhedrosis that can potentially result in an outbreak of infectious disease. The average incidence rate of this disease reaches approximately 30% and occasionally reaches above 70%. The yield of cocoons can be decreased by 30% when AnpeNPV occurs. This viral disease has brought about huge economic losses to sericulturist in China.
AnpeNPV belongs to the lepidopteran family of NPV. The physical map for AnpeNPV genome has been constructed and the genome was estimated at approximately 130 kbp and 128 kbp in size, respectively [18, 19]. Our group and others have reported the sequences of several AnpeNPV genes including lef-8, lef-9, polyhedrin, lef-7, cathepsin, chitinase, genes ie-2, pe38 (GenBank accession No.DQ372717),truncated egt, lef-1 (GenBank accession no.AY846867), odv-e56 (GenBank accession no.AY846866), lef-3 (GenBank accession no.AY846749) and gp64. A total of approximately 30 kbp of the genome of AnpeNPV has been sequenced and released in GenBank. In this paper, we report the complete sequence and organization of the AnpeNPV genome. The phylogeny of AnpeNPV in comparison with previously published baculovirus genomes is also addressed. Sequence and analysis of the AnpeNPV genome suggest significance for the epidemiological studies in the research community.
Results and discussion
Sequencing, assembly, and analysis of the AnpeNPV genome
A total of 536 recombinant plasmids containing cloned viral DNA fragments (with Sau 3AI partial digestion) were obtained. The size of most of the inserted fragments ranged from 2.0 to 5.0 kbp. Both ends of the inserted fragments were sequenced. The complete sequence of AnpeNPV genome was assembled from the sequenced restriction fragments and PCR amplified products. Approximately over 6 times coverage of the genome was achieved.
The entire AnpeNPV dsDNA genome was sequenced and assembled into a contiguous sequence of 126,629 bp with a G+C content of 53.5% [GenBank: NC_008035]. Of the previously published sequences of AnpeNPV, only gp64 [GenBank:AY854039] showed sequence differences (nt155, nt192, nt989, nt1240, nt1327, nt1380, respectively). The reason may be the variation of different virus isolates. The cytosine of the complement of stop codon of polyhedrin was designated as the first nucleotide of the genome. Compared with other baculovirus genomes, the size of AnpeNPV genome was larger than those of NeleNPV, NeseNPV, CuniNPV, EppoNPV, AdhoNPV, PlxyGV, AdorGV, PhopGV and CrleGV. The G+C content of the OpMNPV and LdMNPV genomes was higher than that of the AnpeNPV genome. Only three hr s were identified in AnpeNPV genome, but twenty-four perfect or imperfect dr s were found in the intergenic spacer regions. In addition, 9 dr s were also found in the intragenic spacer regions of AnpeNPV genome. The number of dr s in AnpeNPV is by far the most in any baculovirus genome. The results of the comparisons between AnpeNPV and other baculovirus genomes are listed in additional files [see Additional file 1].
Identification of AnpeNPV ORFs
Classification of genes in AnpeNPV and other baculoviruses
The classification of genes encoded in different baculovirus genomes may help to determine essential genes for virus survival and to understand features including host range, virulence and morphology. The genes found in all baculoviruses are more likely to be essential, whereas auxiliary genes found in only some baculoviruses may give viruses a selective advantage in nature .
We classified 147 ORFs into six groups: 1) ORFs found in all baculoviruses; 2) ORFs found in all lepidopteran baculoviruses; 3) ORFs found in all lepidopteran NPVs; 4) ORFs found in all Group I NPVs; 5) ORFs found in partial Group I NPVs; 6) ORFs that are unique to AnpeNPV. The detailed results were shown in additional files [see Additional file 4] and Figure 1.
In first group, the 29 genes present in all baculovirus genomes reported previously were also found in AnpeNPV genome [see Additional file 4]. Comparison of these 29 genes between AnpeNPV and the baculovirus type species AcMNPV dedicated that their organization was consistent with Group I NPVs. The 29 genes are more likely to be essential for all baculoviruses survival.
Lauzon et al. (2004) identified 33 ORFs present in all lepidopteran baculoviruses. All of these ORFs were found in AnpeNPV [see Additional file 4]. The 33 conserved genes are more likely to be essential for all lepidopteran baculoviruses survival.
In addition to the genes present in all baculovirus and all lepidopteran baculoviruses described above, 12 ORFs in AnpeNPV were identified in all lepidopteran NPVs [see Additional file 4]. Some of these ORFs were also present in one or more GVs. Of these 12 genes, 6 were found exclusively in all lepidopteran NPVs (Anpe26 (ac34), Anpe54 (ac55), Anpe56 (ChaB), Anpe96 (p87), Anpe100 (ac108), Anpe122 (calyx)). De Jong et al.(2005) identified 8 ORFs exclusive to all lepidopteran NPVs including ac17, ac21(arif-1) and ac57. However, due to the absence of ac17 in ChchNPV, ac21(arif-1) and ac57 in AdhoNPV, these 3 ORFs are now excluded from this class of ORFs. Including these 12 genes, there were still 74 (29+33+12) lepidopteran NPV specific ORFs. In the work of de Jong et al.(2005), 74 baculovirus ORFs were also identified in all lepidopteran NPVs, however these 74 are not the same as those identified here [see Additional file 4]. In this class of genes, Anpe14 had only a partial egt (the 5'end), which matched to 5' end of egt gene of other lepidopteran NPVs by BLASTP, and the amino acid identity ranged from 43% (Bm7, e = 1.0e-16) to 70% (Op14, e = 1.0e-24) for Group I NPVs. egt is a auxiliary gene in lepidopteran NPVs, which is not essential for replication but may facilitate a selective advantage in nature .
Another 41 ORFs in AnpeNPV were found in all Group I NPVs. Of the 41 ORFs, 12 were found exclusively in all Group I NPVs by de Jong et al.(2005) and Hyink et al.(2002) [5, 25], all of which were also found in AnpeNPV including Anpe5(ac5), Anpe9 (ptp-1), Anpe15(odv-e26), Anpe33 (ac30), Anpe42 (gta), Anpe67 (ac72), Anpe68 (ac73), Anpe106 (ac114), Anpe114 (ac124), Anpe119 (gp64), Anpe123(ac132) and Anpe141 (ie-2) [see Additional file 4].
4 ORFs of AnpeNPV (Anpe31, Anpe71, Anpe140 and Anpe142) had no recognizable baculovirus homologs [see Additional file 4]. A homologous region between the Anpe142 and Microtuble associated protein of Drosophila melanogaster was found by BLAST search. The microtuble associated protein is required for the integrity of mitotic spindle. Additionally, a region of Anpe31 showed homology to a carbamate kinase gene of Streptococcus ratti.
In addition, another 27 ORFs [see Additional file 4] were considered in the class of ORFs found in some Group I NPVs. These genes may be auxiliary genes and give viruses a selective advantage in nature. In these ORFs, Anpe89 and Anpe108 were determined as bro, a baculovirus repeated ORF, which are found in almost all baculoviruses. CfdefORF142 and HycuORF148 were identified as genes unique to CfDEFNPV and HycuNPV, respectively, but alignment of Anpe143 to Cfdef142 and Anpe3 to Hycu148 revealed high conservation (34% and 62% amino acid identity, respectively). Anpe110 was a homologue of a gene only found in HycuNPV and OpMNPV (Hycu36 and Op118, 35% and 27% amino acid identity, respectively). Anpe98 and Anpe126 were homologs of genes observed only in AcMNPV and RoNPV. Alignment of the deduced amino acid sequence of Anpe98 (pnk/pnl) to ac86 (63% amino acid identity, e = 0) and Ro83 (63% amino acid identity, e = 0) showed high conservation in their amino acid sequences. At the same time, Anpe98, ac86, and Ro83 shared a 175 residue-long NK (Nucleoside/nucleotide kinase) motif. The percentage identity within these regions (98.9%) suggested that they are functionally important regions. Additionally, multiple alignment of Anpe126 (94 k) to ac134 and Ro127 revealed that Anpe126 matched to two regions of ac134 and Ro127 in the 5' and 3'ends (aa1-273 and aa 419–505 in Anpe126 were 53% and 45% identical at the amino acid level to aa 1–277 and aa 712–800 in ac134 and Ro127, respectively. aa278-aa711 in ac134 and Ro127 was truncated in Anpe126). On the other hand, Anpe 126 also matched to 94 k gene of some Group II NPVs and GVs in the two terminal regions. This implied that the function region of 94 k may be in the two regions. Similarly, Anpe115 was homologous to eukaryotic trex gene, which was observed only in CfMNPV, CfDEFNPV and AgMNPV. The trex gene encodes a 3' to 5' exonuclease [26–28]. Also, multiple alignment of Anpe 115 to Cf114 (65% amino acid identity, e = 4e-84) and Cfdef119 (87% amino acid identity, e = 5e-116) appeared to be highly conservative in their amino acid sequences. Furthermore, Anpe115 was 87% (e = 7e-117) identical at the amino acid level to v-trex gene of AgMNPV. Anpe115, Cf114, Cfdef119 and the v-trex gene of AgMNPV had a 171 residue-long EXOIII motif (exonuclease domain in DNA-polymerase alpha and epsilon chain, ribonuclease T and other exonucleases). The percentage identity (100%) within these regions indicated that this region was also functionally important. This kind of genes may give AnpeNPV a selective advantage in nature.
Baculoviruses contain a lot of genes involved in essential and selective functions, such as viral replication, transcription, inhibition of apoptosis, protein structure, assembly and auxiliary. Besides, there were some conserved genes, of which functions were not clearly identified.
All of the conserved function genes, previously identified in all baculoviruses and all lepidopteran-specific baculoviruses , were also found in AnpeNPV. Several variable function genes identified in some baculoviruses and some conserved baculovirus ORFs of unknown function were also found in AnpeNPV. The detailed functions of these genes are listed in additional files [see Additional file 5].
Whereas, several variable DNA replication genes, including helicase-2(ld50), dna-ligase(ld22), dutpase (op31), RNase reductase-1 (op32) and RNase reductase-2 (op34), were not found in AnpeNPV. Five genes have so far been implicated as inhibitors of apoptosis, including iap-1, iap-2, iap-3, iap-4 and p35. However, only iap-1 and iap-2 were identified in AnpeNPV and considered as genes present in all Group I NPVs. Baculovirus genomes typically encode auxiliary genes that are non-essential for replication but provide a selective advantage to the virus . AnpeNPV contained some auxiliary genes including ctl-1 and ctl-2 (conotoxin-like homologue 1 and 2, Anpe8 and Anpe28, respectively), which are small disulfide-rich ion channel antagonists. ctl1 was found in AcMNPV, OpMNPV, LdMNPV, CfMNPV and HycuNPV, but ctl2 was only found in OpMNPV, HycuNPV and LdMNPV. Numerous structural genes were also found in AnpeNPV, which composed the structure of virus. Three genes in AnpeNPV, pif-1(Anpe111), pif-2(Anpe20) and p74 (Anpe129), were identified as structural genes as well as per os infectivity factor involved in baculovirus infection to midgut epithelial cells of the host. Other structural genes identified in some baculoviruses but not in AnpeNPV included homologues to vef-1(Agse75), vef-2(Agse76), vef-3(Agse128), pk2(ac123). Vefs have been found to dramatically increase the infectivity of baculoviruses in non-natural lepidopteran hosts  and only found in GVs and Group II NPV.
Hrs repeats and the genes around them
Most baculoviruses contain hr regions. An individual NPV hr is typically comprised of direct repeats usually centered around a perfect or imperfect palindrome. hr s have been implicated as origins of DNA replication [31–33] and as enhancers of RNA polymerase II-mediated transcription . GV hr s are different from NPV hr s except for PlxyGV, whose hr s appear more similar to those from NPVs. GV hr s usually lack a palindrome and are homologous to all other hr s within the same genome, but there are no two GVs that have similar hr regions. GV hr s are characterized to be repeat regions that are variable and AT-rich [12, 13, 15, 16].
Direct repeats (dr s) as described by Garcia-Maruniak et al.(2004) are suggested to be a putative non-hr, ori-like regions containing an unusual AT-rich sequences and perfect or imperfect short palindromes in intergenic spacer or within ORFs. These dr s are not repeated elsewhere in the genome. With the exception of NeleNPV, the majority of NPVs have hr s and many of these hr s have homologous palindromes core. The NeleNPV has only nine dr s shared a closer similarity to the organization of repeat regions in GVs than to those in NPVs and showed little similarity to typical NPV hr s, which are characterized by the presence of perfect or imperfect palindrome and AT rich . In comparison NeseNPV not only has the hr s, but also has four dr s. Similar to NeseNPV, 24 perfect or imperfect dr s were identified by Tandem Repeats Finder(TRF) software in addition to the hr s in AnpeNPV [see Additional file 7]. These dr s contained from 2 to 7 repeats. Except for dr 11, they were found to be AT-rich, especially dr 13(94.4%), dr 14(92.9%), dr 15(92.2%), dr 10(91.4%), dr 16(91.4%), dr 12(89.6%), dr 9(88.1%), dr 24(87.2%), dr 17(86.2%), dr 5(86.1%), dr 20(85.0%), dr 8(84.8%), dr 2(82.5%) and dr 6(81.1%). Palindromes ranging in size from 10 to 18 bp were found in several of the dr s. All of the dr s were located inside intergenic spacer regions of the AnpeNPV genome. The existence of twenty-four dr s in the AnpeNPV intergenic spacer regions may reflect a different genome replication strategy for the NPVs, and these dr s may be a new type putative non-hr, ori-like region as is found in GrleGV, CpGV, and AdorGV. In addition, nine tandem repeats were also identified in ORFs of AnpeNPV genome with a high GC content [see Additional file 8]. At the same time, a specific GT-repeat was located from 122228 to122259nt. The sequences contained continuous GT residues with dinucleotide (GT) repeats. The GT-repeat was only found in OpMNPV. Similar sequences of GT repeats have been characterized in other organisms, e.g., Drosophila virilis and Caenorhabditis elegans .
Comparison of AnpeNPV and other baculoviruses
The gene order of AnpeNPV was compared with those of other Group I NPVs (HycuNPV, EppoNPV, CfMNPV, OpMNPV, AcMNPV, RoNPV, BmNPV and CfDEFNPV), two Group II NPVs (MacoNPV A, HearNPV) and a GV (XcGV) by gene parity plots . AnpeNPV shared 133, 125, 130, 130, 127, 125, 123 and 133 ORFs with HycuNPV, EppoNPV, CfMNPV, OpMNPV, AcMNPV, RoNPV, BmNPV and CfDEFNPV, respectively [see Additional file 2]. Among the Group II NPVs, AnpeNPV shared the highest number of ORFs with MacoNPV A (92 ORFs) and the least with AdohNPV (83 ORFs).
Similarity between the 29 core baculovirus genes between AnpeNPV and other Group I NPV was obviously higher than the similarities between AnpeNPV and other groups [see Additional file 9]. CuniNPV was particularly dissimilar, of which similarities were 15.1%–51.1% range and average of 24.3%.
Phylogeny of AnpeNPV
In summary, we found that AnpeNPV is a member of Group I NPV and is most closely related to HycuNPV OpMNPV, EppoNPV, CfMNPV, CfDEFNPV and to a lesser extent to BmNPV, RoNPV and AcMNPV based on the results of multi-alignment, phylogenetic analysis and gene parity plots. Analysis of the genome organization of AnpeNPV and other Group I NPVs demonstrated that there are one gene rearrange region in AnpeNPV with CfDEFNPV, and two gene rearrange regions in AnpeNPV with BmNPV, RoNPV and AcMNPV. Despite its close relation to OpMNPV, CfMNPV, HycuNPV, EppoNPV and CfDEFNPV, the AnpeNPV genome has several distinct features including the presence of 4 genes unique to AnpeNPV(Anpe31, Anpe71, Anpe140 and Anpe142) and 2 genes previously identified as new genes unique to CfDEFNPV and HycuNPV, respectively(Anpe3/Hycu148 and Anpe143/Cfdef142), and that, unlike other baculoviruses there are fewer hr s (only 3 in AnpeNPV), more direct repeats (24 in AnpeNPV intergenic spacer regions, 9 in AnpeNPV intragenic spacer regions) and a GT-repeat previously found only in OpMNPV. At the same time, analysis of the genes around the hr s of AnpeNPV supported the argument that an ancestral baculovirus acquired genes from ancestral insect hosts that were inserted into the regions around the hr s. There are still 29 conserved genes present in all baculoviruses and 33 genes present in all lepidopteran baculoviruses. Besides genes present in all baculoviruses and lepidopteran baculoviruses, 12 ORFs in AnpeNPV were identified in all lepidopteran NPVs. Of the 12 ORFs, 6 were found exclusively in all lepidopteran NPVs (Anpe26/ac34, Anpe54/ac55, Anpe56/ChaB, Anpe96/p87, Anpe100/ac108, Anpe122/calyx), therefore the number of genes common to all lepidopteran NPVs is 74, however these 74 genes are not the same as those identified before. In addition, 12 ORFs are now exclusive to and present in all Group I NPVs(Anpe5/ac5,, Anpe9/ptp-1, Anpe15/odv-e26, Anpe33/ac30, Anpe42/gta, Anpe67/ac72, Anpe68/ac73, Anpe106/ac114, Anpe114/ac124, Anpe119/gp64, Anpe123/ac132, Anpe141/ie-2). Only two bro genes (Anpe89/bro-b and Anpe108/bro-a) and two inhibitors of apoptosis (Anpe36/iap-1 and Anpe66/iap-2) were found in AnpeNPV. So far, two ctls,ctl1 and ctl2, are found only in the genomes of AnpeNPV, HycuNPV, OpMNPV and LdMNPV. These differences identified between AnpeNPV and other baculovirus genomes might provide insight into their evolution, host specificity and pathogenicity. The sequence and analysis of AnpeNPV complete genome can provide a continuing resource to better follow the evolution and phylogeny of baculovirus.
Virus isolation, construction of genomic DNA libraries
The AnpeNPV virus (Liaoning strain) was isolated from the natural infected tussah silkworm from a rearing yard in Liaoning province, China by the Liaoning Sericultural Research Institute. To generate a large number of polyhedra, healthy fifth instar tussah silkworms were inoculated in the laboratory and the haemolymph was collected from the sick tussah silkworms, then washed with ddH2O and centrifuged several rounds to get pure AnpeNPV polyhedra. The genomic DNA of AnpeNPV was purified according to following protocol: about 5 × 108 polyhedra was dissolved in solution (0.1 M Na2CO3, 0.15 M NaCl, pH10.4) on ice for 10 minutes, SDS was then added to the final concentration 0.5%, and the solution was kept on the ice for another 10 minutes. The genomic DNA was extracted by the addition of an equal volume phenol (pH8.0) two times and chloroform one time. The DNA was precipitated with 2-volumes ethanol, washed with 70% ethanol, and dissolved in 0.1 × TE buffer (pH8.0) . The random viral genomic DNA library was constructed as follows: The viral DNA was partially digested by Sau 3A I, a restrict endonuclease with a 4-bp recognition site, and the fragments ranging from 2.0 to 5.0 kbp were recovered from an agarose gel. The cohesive ends were filled in partially by incubating with dATP and dGTP in the presence of Klenow fragment of DNA polymerase I. The plasmid vector pUC19 was fully digested by Sal I, and subsequently filled in by adding dTTP, dCTP and the Klenow fragment. The vector and genomic fragments were mixed and ligated. The ligation products were transformed into E. coli competent cells TG1 and then cultured on LB plates with X-gal/IPTG and Amp (100 μg/ml) . White colonies were selected and cultured. A one-step screening of recombinant plasmids was conducted by assessing the size of insertion . The recombinant plasmids were sequenced with plasmid specific primers and 'primer nesting' from both strands by using BigDye Terminator v3.1 (ABI) on Genetic analyzer 3130XL(ABI).
Assembly and sequence analysis
Restriction fragments from recombinant plasmids were sequenced and assembled into contigs using ContigExpress9.1.0 and SeqMan5.0 from the DNASTAR software package. PCR was used to generate gap-spanning fragments and low quality data regions after preliminary assembly. Open reading frames (ORFs) were identified using ORF finder . The criterion for defining an ORF was a size of at least 50 amino acids and minimal overlap. All BLAST searches were done through the national center for biotechnology information (NCBI) website using BLAST2.2.14  and local BLAST using PowerBlast1.2.0. Multiple alignments and percentage identities were performed using the Clustal W and MegAlign from the DNASTAR software package. Genome homology plots were generated using GeneQuest5.0 from the DNASTAR software package. Dr s were identified by Tandem Repeats Finder (TRF) software with default settings. hr s containing palindromes were identified by searching the genome using hr consensus sequences reported previously . Phylogenetic trees were constructed using Clustal W protein alignment based on the 29 core baculovirus genes.
This work was supported by financial grants from the National Basic Research Program of China (973 Program, No. 2005CB121006), the National High Technology Research and Development Program of China (863 Major Program, 2005AA206120), the National Natural Science Foundation of China (No.30670082, No.30470034) and the Zhejiang Natural Science Foundation (Y305275). We thank Shizuo G. Kamita for English corrections.
- Ahrens CH, Russell RLQ, Funk CJ, Evans JT, Harwood SH, Rohrmann GF: The sequence of the Orgyia pseudotsugata multinucleocapsid nuclear polyhedrosis virus genome. Virology. 1997, 229: 381-399. 10.1006/viro.1997.8448.PubMedView ArticleGoogle Scholar
- Blissard GW, Black B, Crook N, Keddie BA, Possee R, Rohrmann G, Theilmann D, Volkman L: Family Baculoviridae. In Virus Taxonomy. The Seventh Report of the International Committee on Taxonomy of Viruses. 2000, San Diego, Academic PressGoogle Scholar
- ICTVdB Management: 00.006 Baculoviridae. In ICTVdB – The Universal Virus Database, version 3. The Earth Institute and Department of Epidemiology. 2003, New York, Mailman School of Public Health, Columbia UniversityGoogle Scholar
- Afonso CL, Tulman ER, Lu Z, Balinsky CA, Moser BA, Becnel JJ, Rock DL, Kutish GF: Genome sequence of a baculovirus pathogenic for Culex nigripalpus . J Virol. 2001, 75: 11157-11165. 10.1128/JVI.75.22.11157-11165.2001.PubMed CentralPubMedView ArticleGoogle Scholar
- de Jong JG, Lauzon HA, Dominy C, Poloumienko A, Carstens EB, Arif BM, Krell PJ: Analysis of the Choristoneura fumiferana nucleopolyhedrovirus genome. J Gen Virol. 2005, 86: 929-943. 10.1099/vir.0.80490-0.PubMedView ArticleGoogle Scholar
- Federici BA: Baculovirus pathogenesis. In The baculoviruses. 1997, New York, PlenumGoogle Scholar
- Lauzon HA, Lucarotti CJ, Krell PJ, Feng Q, Retnakaran A, Arif BM: Sequence and organization of the Neodiprion lecontei nucleopolyhedrovirus genome. J Virol. 2004, 78: 7023-7035. 10.1128/JVI.78.13.7023-7035.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Herniou EA, Olszewski JA, Cory JS, O'Reilly DR: The genome sequence and evolution of baculoviruses. Annu Rev Entomol. 2003, 48: 211-234. 10.1146/annurev.ento.48.091801.112756.PubMedView ArticleGoogle Scholar
- Garcia-Maruniak A, Maruniak JE, Zanotto PM, Doumbouya AE, Liu JC, Merritt TM, Lanoie JS: Sequence analysis of the genome of the Neodiprion sertifer nucleopolyhedrovirus. J Virol. 2004, 78: 7036-7051. 10.1128/JVI.78.13.7036-7051.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Ayres MD, Howard SC, Kuzio J, Lopez-Ferber M, Possee RD: The complete DNA sequence of Autographa californica nuclear polyhedrosis virus. Virology. 1994, 202: 586-605. 10.1006/viro.1994.1380.PubMedView ArticleGoogle Scholar
- Duffy SP, Young AM, Morin B, Lucarotti CJ, Koop BF, Levin DB: Sequence analysis and organization of the Neodiprion abietis nucleopolyhedrovirus genome. J Virol. 2006, 80: 6952-6963. 10.1128/JVI.00187-06.PubMed CentralPubMedView ArticleGoogle Scholar
- Hayakawa T, Ko R, Okano K, Seong SI, Goto C, Maeda S: Sequence analysis of the Xestia c-nigrum granulovirus genome. Virology. 1999, 262: 277-297. 10.1006/viro.1999.9894.PubMedView ArticleGoogle Scholar
- Hashimoto Y, Hayakawa T, Ueno Y, Fujita T, Sano Y, Matsumoto T: Sequence analysis of the Plutella xylostella granulovirus genome. Virology. 2000, 275: 358-372. 10.1006/viro.2000.0530.PubMedView ArticleGoogle Scholar
- Luque T, Finch R, Crook N, O'Reilly DR, Winstanley D: The complete sequence of the Cydia pomonella granulovirus genome. J Gen Virol. 2001, 82: 2531-2547.PubMedView ArticleGoogle Scholar
- Wormleaton S, Kuzio J, Winstanley D: The complete sequence of the Adoxophyes orana granulovirus genome. Virology. 2003, 311: 350-365. 10.1016/S0042-6822(03)00149-1.PubMedView ArticleGoogle Scholar
- Lange M, Jehle JA: The genome of the Cryptophlebia leucotreta granulovirus. Virology. 2003, 317: 220-236. 10.1016/S0042-6822(03)00515-4.PubMedView ArticleGoogle Scholar
- Escasa SR, Lauzon HA, Mathur AC, Krell PJ, Arif BM: Sequence analysis of the Choristoneura occidentalis granulovirus genome. J Gen Virol. 2006, 87: 1917-1933. 10.1099/vir.0.81792-0.PubMedView ArticleGoogle Scholar
- Huang YJ, Kobayashi J, Yoshimura T: Genome mapping and gene analysis of Antheraea pernyi nucleopolyhedrovirus for improvement of baculovirus expression vector system. J Biosci Bioeng. 2002, 93 (2): 183-191. 10.1263/jbb.93.183.PubMedView ArticleGoogle Scholar
- Wang WB, Zhu SY, Wang LQ, Yu F, Shen WD: Cloning and sequence analysis of the Antheraea pernyi nucleopolyhedrovirus gp64 gene. J Biosci. 2005, 30: 605-610.PubMedView ArticleGoogle Scholar
- Jehle JA, Lange M, Wang H, Hu Z, Wang Y, Hauschild R: Molecular identification and phylogenetic analysis of baculoviruses from Lepidoptera. Virology. 2006, 346: 180-193. 10.1016/j.virol.2005.10.032.PubMedView ArticleGoogle Scholar
- Theilmann DA, Stewart S: Tandemly repeated sequence at the 3' end of the IE-2 gene of the baculovirus Orgyia pseudotsugata multicapsid nuclear polyhedrosis virus is an enhancer element. Virology. 1992, 187: 97-106. 10.1016/0042-6822(92)90298-4.PubMedView ArticleGoogle Scholar
- Lauzon HAM, Lucarotti CJ, Krell PJ, Feng Q, Retnakaran A, Arif BM: Sequence and organization of the Neodiprion lecontei nucleopolyhedrovirus genome. J Virol. 2004, 78: 7023-7035. 10.1128/JVI.78.13.7023-7035.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Van Oers MM, Abma-Henkens MH, Herniou EA, de Groot JCW, Peters S, Vlak JM: Genome sequence of Chrysodeixis chalcites nucleopolyhedrovirus, a baculovirus with two DNA photolyase genes. J Gen Virol. 2005, 86: 2069-2080. 10.1099/vir.0.80964-0.PubMedView ArticleGoogle Scholar
- Nakai M, Goto C, Kang WK, Shikata M, Luque T, Kunimi Y: Genome sequence and organization of a nucleopolyhedrovirus isolated from the smaller tea tortrix, Adoxophyes honmai. Virology. 2003, 316: 171-183. 10.1016/j.virol.2003.08.002.PubMedView ArticleGoogle Scholar
- Hyink O, Dellow RA, Olsen MJ, Caradoc-Davies KMB, Drake K, Herniou EA, Cory JS, O'Reilly DR, Ward VK: Whole genome analysis of the Epiphyas postvittana nucleopolyhedrovirus. J Gen Virol. 2002, 83: 957-971.PubMedView ArticleGoogle Scholar
- Yang DH, de Jong JG, Makhmoudova A, Arif BM, Krell PJ: Choristoneura fumiferana nucleopolyhedrovirus encodes a functional 3'-5' exonuclease. J Gen Virol. 2004, 85: 3569-3573. 10.1099/vir.0.80592-0.PubMedView ArticleGoogle Scholar
- Slack JM, Shapiro M: Anticarsia gemmatalis multicapsid nucleopolyhedrovirus v-trex gene encodes a functional 3' to 5' exonuclease. J Gen Virol. 2004, 85: 2863-2871. 10.1099/vir.0.80109-0.PubMedView ArticleGoogle Scholar
- Slack JM, Ribeiro BM, de Souza ML: The gp64 locus of Anticarsia gemmatalis multicapsid nucleopolyhedrovirus contains a 3' repair exonuclease homologue and lacks v-cath and ChiA genes. J Gen Virol. 2004, 85: 211-219. 10.1099/vir.0.19617-0.PubMedView ArticleGoogle Scholar
- Miller LK: Introduction to the Baculoviruses. In The Baculoviruses. 1997, New York, PlenumView ArticleGoogle Scholar
- Derksen ACG, Granados RR: Alteration of a lepidopteran peritrophic membrane by baculoviruses and enhancement of viral infectivity. Virology. 1988, 167: 242-250. 10.1016/0042-6822(88)90074-8.PubMedView ArticleGoogle Scholar
- Ahrens CH, Rohrmann GF: Replication of Orgyia pseudotsugata baculovirus DNA: lef-2 and ie-1 are essential and ie-2,p34, and Op-iap are stimulatory genes. Virology. 1995, 212: 650-662. 10.1006/viro.1995.1523.PubMedView ArticleGoogle Scholar
- Kool M, Ahrens CH, Vlak JM, Rohrmann GF: Replication of baculovirus DNA. J Gen Virol. 1995, 76: 2103-2118.PubMedView ArticleGoogle Scholar
- Lee H, Krell PJ: Reiterated DNA fragments in defective genomes of Autographa californica nuclear polyhedrosis virus are competent for AcMNPV-dependent DNA replication. Virology. 1994, 202: 418-429. 10.1006/viro.1994.1358.PubMedView ArticleGoogle Scholar
- Huang QH, Deveraux QL, Maeda S, Salvesen GS, Stennicke HR, Hammock BD, Reed JC: Evolutionary conservation of apoptosis mechanisms:Lepidopteran and baculoviral inhibitor of apoptosis proteins are inhibitors of mammalian caspase-9. PNAS. 2000, 97: 1427-1432. 10.1073/pnas.97.4.1427.PubMed CentralPubMedView ArticleGoogle Scholar
- Verhagen MA, Coulson EJ, Vaux DL: Inhibitor of apoptosis proteins and their relatives: IAPs and other BIRPs. Genome Biol. 2001, 2 (7): REVIEWS3009-10.1186/gb-2001-2-7-reviews3009.PubMed CentralPubMedView ArticleGoogle Scholar
- Tautz D, Renz M: Simple DNA sequences of Drosophila virilis isolated by screening with RNA. J Mol Biol. 1984, 172: 229-235. 10.1016/S0022-2836(84)80041-8.PubMedView ArticleGoogle Scholar
- Wilson R, Ainscough R, Anderson K, Baynes C, Berks M, Bonfield J, Burton JL: 2.2 Mb of contiguous nucleotide sequence from chromosome III of C. elegans. Nature. 1994, 368: 32-38. 10.1038/368032a0.PubMedView ArticleGoogle Scholar
- Hu ZH, Arif BM, Jin F, Martens JW, Chen XW, Sun JS, Zuidema D, Goldbach RW, Vlak JM: Distinct gene arrangement in the Buzura suppressaria single-nucleocapsid nucleopolyhedrovirus genome. J Gen Virol. 1998, 79: 2841-2851.PubMedView ArticleGoogle Scholar
- Ikeda M, Shikata M, Shirata N, Chaeychomsri S, Kobayashi M: Gene organization and complete sequence of the Hyphantria cunea nucleopolyhedrovirus genome. J Gen Virol. 2006, 87: 2549-2562. 10.1099/vir.0.81930-0.PubMedView ArticleGoogle Scholar
- Xie Y, Li ZY, Wang S: Review of species invasion into China in Protecting the biological diversity of China. 1996, Beijing, China Environmental Science Press, 91-106.Google Scholar
- Aruga H, Yoshitake N, Watanabe H, Hukuhara T: Studies on nuclear polyhedroses and their inductions in some Lepidoptera. Jpn J Appl Entomol Zool. 1960, 4: 51-56.View ArticleGoogle Scholar
- Xia AF, Zhang ZF, Li JS, Li WG, Mu ZM: Cloning and Analysis of occlusion-derived Virus Envelope-56 Gene of Antheraea pernyi Nucleopolyhedrovirus. Canye Kexue(china). 2005, 31: 300-305.Google Scholar
- Sambrook J, Russell DW: Plasmids and Their Usefulness in Molecular Cloning. In Molecular Cloning:a laboratory manual. 2000, Cold Spring Harbor, Cold Spring Harbor pressGoogle Scholar
- Beuken E, Vink C, Bruggeman CA: One-step procedure for screening recombinant plasmids by size. Biotechniques. 1998, 24: 748-750.PubMedGoogle Scholar
- ORF Finder website. [http://www.ncbi.nlm.nih.gov/gorf]
- NCBI Blast website. [http://www.ncbi.nlm.nih.gov/BLAST]
- Lauzon HA, Jamieson PB, Krell PJ, Arif BM: Gene organization and sequencing of the Choristoneura fumiferana defective nucleopolyhedrovirus genome. J Gen Virol. 2005, 86: 945-961. 10.1099/vir.0.80489-0.PubMedView ArticleGoogle Scholar