Laboratory of Virology, Wageningen University and Research Centre, Binnenhaven 11, 6709 PD Wageningen, The Netherlands1
Author for correspondence: Just Vlak. Fax +31 317 484820. e-mail just.vlak{at}medew.viro.wau.nl
![]() |
Abstract |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Baculoviruses almost exclusively infect insects, belonging mainly to the orders Lepidoptera, Hymenoptera and Diptera (Adams & McClintock, 1991 ). SeMNPV infects only a single insect species, the beet army worm Spodoptera exigua (Lepidoptera: Noctuidae), a worldwide insect pest of agricultural importance. SeMNPV differs from many other baculoviruses in that it is monospecific and highly virulent for S. exigua larvae (Smits, 1987
). However, the molecular mechanism associated with these properties is unknown. Therefore, it is important at this point to study the genetic information available for the virus and the expression of its genes.
The best-characterized baculoviruses are Autographa californica (Ac) MNPV (Ayres et al., 1994 ), Bombyx mori (Bm) NPV (Gomi et al., 1999
), Orgyia pseudotsugata (Op) MNPV (Ahrens et al., 1997
) and Lymantria dispar (Ld) MNPV (Kuzio et al., 1999
). The genome of AcMNPV is composed of 133894 bp, potentially encoding 154 proteins (Ayres et al., 1994
). Fourteen ORFs are unique to AcMNPV, whereas most of the other ORFs have baculovirus homologues. Eight homologous regions (hr), implicated as transcriptional activators and as putative origins of replication, are present in the AcMNPV genome. BmNPV is 128413 bp in size and contains 136 putative genes (Gomi et al., 1999
). Only four ORFs are unique to BmNPV, all other ORFs possessing a baculovirus homologue. Five copies of an AcMNPV ORF2 homologue, named bro, and seven hr sequences are present, dispersed along the BmNPV genome. The OpMNPV genome contains 131990 bp and potentially encodes 152 proteins (Ahrens et al., 1997
). Twenty-six genes are unique to OpMNPV. OpMNPV contains one complete bro gene and two truncated ORFs that show homology to the bro genes. Only five hr sequences are present, dispersed throughout the OpMNPV genome. AcMNPV, BmNPV and OpMNPV belong to the group I NPVs (Zanotto et al., 1993
). The genome of LdMNPV (group II) is composed of 161046 bp and contains 163 ORFs (Kuzio et al., 1999
). Forty-seven genes are unique to LdMNPV. The large size of the LdMNPV genome is largely due to the presence of 13 hr sequences and 16 bro gene homologues.
A detailed physical map was recently constructed for an American isolate of SeMNPV (US1) in order to map a mutant SeMNPV. This mutant was obtained within the first passage in insect cell culture. It has a single deletion of approximately 25 kb and is unable to infect S. exigua larvae (Heldens et al., 1996 ). The sequences of a number of SeMNPV genes, including polyhedrin (van Strien et al., 1992
), p10 (Zuidema et al., 1993
), ubiquitin (van Strien et al., 1996
), ribonucleotide reductase large subunit (van Strien et al., 1997
) and p143 (Heldens et al., 1997b
), have been elucidated and characterized. The locations of these genes on the SeMNPV genome differ considerably from those of other baculoviruses, such as AcMNPV, BmNPV, OpMNPV and LdMNPV, suggesting that the genetic organization is markedly different. It was noted, however, that among baculoviruses the genomic region located in the centre of the linearized genome was highly conserved (Heldens et al., 1997b
). Six hr sequences have been identified on the SeMNPV genome, which are similar in structure to those of other baculoviruses (Broer et al., 1998
). Here, we present the complete sequence and gene organization of the SeMNPV genome and compare them to other baculoviruses by genomic and phylogeny analysis.
![]() |
Methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
The SeMNPV XbaI plasmid and Sau3AI cosmid libraries were described previously (Heldens et al., 1996 ). The XbaI-A and XbaI-B fragments were too large to be cloned into pUC19. The XbaI-A fragment was subcloned from SeMNPV cosmids 24 and 17 into plasmids SeBglII-H, SePstI-M (cosmid 24) and SeBSpeI-5.4, SeSpeI-H and SeBP-5.6 (cosmid 17). In addition, the XbaI-B fragment was subcloned from SeMNPV cosmid 22 into plasmids SeSpeIH-3.2, SeBSpeI-6.3, SeSpeIH-2.8, SeHBglII-6.2 and SeEcoRI-2.2. Some of these clones were described previously (Broer et al., 1998
; van Strien et al., 1996
).
Four regions of the SeMNPV genome were difficult to clone, as restriction fragments or sequencing attempts resulted in premature termination. The four regions were located at positions 1653019106, 3283136808, 4483946417 and 5475854937 on the SeMNPV genome. These regions were amplified by PCR and cloned into pGEM-T vectors (Promega). Template DNA for sequencing was purified from plasmids by using Jetstar columns according to manufacturers protocol (ITK Diagnostics).
Sequencing was done by using plasmid, cosmid and PCR products from both strands of the viral genomic DNA as templates. Sequence reactions were performed at the Sequencing Core Facilities of Wageningen Agricultural University and Queens University (Kingston, Ontario, Canada) by primer walking.
DNA sequence analysis.
Genomic DNA composition, structure, repeats and restriction enzyme pattern were analysed with the Wisconsin Genetics Computer Group programs (Devereux et al., 1984 ) and DNASTAR. ORFs consisting of more than 50 amino acids were considered to encode proteins. Relevant ORFs (119 of 139) were checked for maximum alignment with known baculovirus gene homologues from GenBank; ORFs with significant overlap of hr sequences were excluded. The overlap between any two ORFs with known baculovirus homologues was set to a maximum of 25 amino acids; otherwise the largest ORF was selected. DNA and protein comparisons with entries in the genetic databases were performed with FASTA and BLAST programs (Pearson, 1990
; Altschul et al., 1990
). Multiple sequence alignments were performed with the GCG PileUp and Gap computer programs with gap creation and extension penalties set to 9 and 2, respectively (Devereux et al., 1984
). Percentage identity indicates the percentage of identical residues between two complete sequences. Motif searches were done against the Prosite release 14 database (Fabian et al., 1997
; Bairoch et al., 1997
). Prediction of transmembrane domains was accomplished with SignalP and PHD software (Nielsen et al., 1997
; Rost & Sander, 1993
). GeneParityPlot analysis was performed on the SeMNPV genome versus the genomes of AcMNPV, BmNPV, OpMNPV and LdMNPV as described previously (Hu et al., 1998
).
![]() |
Results and Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
One hundred and thirty-nine ORFs, defined as methionine-initiated ORFs encoding more than 50 amino acids and with minimal overlap with other ORFs, were present on the SeMNPV genome (Fig. 1). The SeMNPV ORFs were, in general, tightly packed with minimal intergenic distances and their orientation was distributed almost evenly along the genome (55% clockwise, 45% anticlockwise; Fig. 1
). The locations, orientations and sizes of the predicted ORFs are shown in detail in Table 1
. The distribution of the ATG, TAG and TGA codons in the SeMNPV sequence was not random, while the TAA frequency (1·58%) was not significantly different from the expected random distribution (1·56%). The ATG codon (1·77%) and TGA stop codon (1·78%) occurred more frequently in the SeMNPV sequence, while there was paucity of TAG stop codons (0·83%), as is the case for AcMNPV (Ranjan & Hasnain, 1995
). Predicted ORFs represented 90% coding density, with a mean ORF length of 875 nucleotides. Twenty ORFs had small (<25 aa) overlaps with adjacent ORFs. One hundred and nineteen (86%) of the 139 SeMNPV ORFs had an assigned function or had homologues among other baculovirus genes (Table 1
). Twenty ORFs are so far unique to SeMNPV. These ORFs accounted for 7% (9·3 kb) of the genome. Six hr sequences similar in structure to those of other baculovirus hr sequences have been identified previously (Broer et al., 1998
) and no further hr sequences were detected in the complete sequence. Furthermore, one non-hr, putative origin of replication, is present on the XbaI-F fragment (Heldens et al., 1997a
). The positions of the hr sequences are presented in Fig. 1
.
|
|
|
The overall characteristics of the different baculovirus genomes are shown in Table 2. The G+C content of the SeMNPV genome was 43·8 mol%, which is similar to that of AcMNPV (Ayres et al., 1994
) and BmNPV (Gomi et al., 1999
) but much lower than that of OpMNPV (Ahrens et al., 1997
) and LdMNPV (Kuzio et al., 1999
) (Table 2
). The smaller number of SeMNPV ORFs compared with AcMNPV and OpMNPV, the genomes of which are similar in size, is caused by the absence of a number of small putative ORFs. In AcMNPV and OpMNPV, these small putative ORFs are located between larger ORFs. The larger ORFs have homologues in many other baculoviruses. SeMNPV does not possess these smaller ORFs, but does possess the larger ORFs. This suggests that the numbers of ORFs in AcMNPV and OpMNPV may be overestimates and that these smaller putative ORFs may not be functional. The frequency of different temporal consensus promoter elements showed considerable variation between baculoviruses, except for the late promoter motif (Table 2
).
|
|
Fifty-three AcMNPV genes had no homologues in SeMNPV (Table 4). Most of these genes are also absent in at least one of the other three baculoviruses compared. However, SeMNPV also lacks AcMNPV ORFs 4, 11, 38, 111 and 115, which are present in the four other baculoviruses. To date, no functions have been assigned to these ORFs. In addition to the AcMNPV ORFs that have no SeMNPV homologues, a number of BmNPV (1), OpMNPV (17) and LdMNPV (33) ORFs without AcMNPV homologues are also absent in SeMNPV (Table 4
). The total number of AcMNPV, BmNPV, OpMNPV and LdMNPV ORFs without homologues in SeMNPV is 53, 7, 40 and 78, respectively. Although most of the baculovirus ORFs not found in SeMNPV have not yet been characterized, information is available for some (Ayres et al., 1994
; Ahrens et al., 1997
; Kuzio et al., 1999
). The presence or absence of baculovirus ORF homologues in SeMNPV and their implications for SeMNPV characteristics are discussed below.
|
Surprisingly, SeMNPV lacks a homologue of the budded virus (BV) surface glycoprotein gp64 (Ac128). A similar situation exists for LdMNPV (Kuzio et al., 1999 ). GP64 is a major envelope glycoprotein that is acquired by virions during budding through the plasma membrane. GP64 is required for efficient virion budding in AcMNPV; deletion of the cytoplasmic tail domain resulted in a reduction in progeny BV and in a virus that was incapable of efficient propagation in cell culture (Oomens & Blissard, 1999
). It has been suggested that Ld130, which shares 22% identity with Ac23, could substitute for the lack of GP64 (Kuzio et al., 1999
). This supposition was based on the presence of N-terminal signal and transmembrane domains, which are indicative of transmembrane receptor-like proteins. SeMNPV has a homologue of Ld130, Se8, which has twice the identity to LdMNPV (41%) than to Ac23, Bm14 and Op21 (~22%).
DNA replication genes
The genes essential for DNA replication were only moderately conserved: DNA pol, helicase, lef-2 and lef-1 were 44% identical, whereas lef-3 showed approximately 30% identity. Se126, a homologue of Ac25 and Bm16, which encodes a putative DNA-binding protein (DBP) (Okano et al., 1999 ; Mikhailov et al., 1998
), showed approximately 30% identity to its homologues. The relatively low identity of these SeMNPV proteins to their homologues may explain the specificity of the virus DNA replication process (Heldens et al., 1997b
).
None of the ie-2, pe38, lef-7 or p35 genes, found to be stimulatory in AcMNPV and BmNPV DNA replication assays, had a homologue in the SeMNPV genome or the LdMNPV genome (Kuzio et al., 1999 ). The ie-2, lef-7 and p35 genes are also non-essential for BmNPV virus replication, since functional deletion by insertion experiments resulted in viable virus mutants (Gomi et al., 1997
). A reduction of viral DNA synthesis, however, was demonstrated in only two of the three cell lines infected with Ac
lef-7 (Chen & Thiem, 1997
).
Genes regulating gene expression
The genes required for transactivation of early gene transcription, such as ie-1, ie-0 and me53, were poorly conserved in their amino acid sequence (~35%) among baculoviruses, whereas the late transcription activators including the RNA polymerase, lef-4 (identity to AcMNPV homologue 51%), lef-8 (64%), lef-9 (63%), p47 (56%) and vlf-1 (65%) were, in contrast, very well conserved. This is compatible with the supposition that specificity is already displayed early in infection because baculoviruses have to adapt to the host transcription system. This could also be explained as the result of a higher constraint in the late and virus-encoded transcription system.
SeMNPV had a CG30 (Ac88) homologue, Se76, that is absent from the LdMNPV genome. This ORF contains features characteristic of a transcription-regulatory protein: (i) two nucleic acid-binding sites, (ii) a zinc finger and (iii) a leucine zipper. It is considered to be a prime candidate in the regulation of genes at late times in infection. The SeMNPV homologue was extended compared with the AcMNPV, BmNPV and OpMNPV genes, as the sequence ACC(G/A)TCGACATCGGC(C/T)GG was repeated seven times. The zinc finger and the leucine zipper were present in Se76, although one of the four leucines was changed to a methionine, as was also found for the OpMNPV homologue.
Inhibitors of apoptosis
Baculoviruses have genes involved in the inhibition of apoptosis: p35-like genes and iap-like genes. SeMNPV lacked homologues of iap-1, iap-4 and p35. The AcMNPV annihilator mutant (Acp35) causes cell line-specific apoptosis after infection. This is in contrast to AcMNPV wild-type and iap-1 or iap-2 deletion recombinants. This suggests that iap-1 and iap-2 are not required for prevention of apoptosis in these cell lines (Griffiths et al., 1999
). OpMNPV and Cydia pomonella GV iap-3 have proven to be inhibitors of apoptosis in different cell lines upon infection with Acp35
recombinant virus (Vucic et al., 1998
; Seshagiri & Miller, 1997
; Ahrens & Rohrmann, 1995
; Lu & Miller, 1995
; Clem et al., 1994
; Clem & Miller, 1994
; Birnbaum et al., 1994
; Crook et al., 1993
). SeMNPV possessed iap-2 and iap-3 homologues. The iap-3 gene product shared 48% similarity with its OpMNPV homologue (Table 1
) and may be involved in the prevention of apoptosis in S. exigua larvae and different S. exigua cell lines.
Nucleotide metabolism
SeMNPV possesses a number of previously described baculovirus genes involved in nucleotide metabolism. Genes encoding the large and small subunits of ribonucleotide reductase (rr1 and rr2) and a dUTPase were present in SeMNPV, as well as in OpMNPV and LdMNPV. By means of these proteins, SeMNPV may promote deoxyribonucleotide synthesis in non-dividing cells and conversion of dUTP to dUMP, which serves as a precursor for dTTP (Elledge et al., 1992 ). The rr1 gene of SeMNPV, Se139, has been described previously (van Strien et al., 1997
). The SeMNPV rr2 gene, Se45, which was found distal from rr1, was more closely related to Ld120 (rr2b) than to the Op34 (rr2) or Ld147 (rr2a) homologues in terms of protein identity and location on the genome (Fig. 2C
, D
). This is in contrast to the rr1 gene, which is equally related to its OpMNPV and LdMNPV homologues. In contrast to OpMNPV (ORF32 and 34) and LdMNPV (ORF 147 and 148), the SeMNPV rr1 and rr2 as well as the LdMNPV rr2b genes appear to have been acquired from a source more closely related to eukaryote than prokaryote homologues (van Strien et al., 1997
; Kuzio et al., 1999
).
The dutpase gene of SeMNPV, Se55, differed more from the dutpase of OpMNPV (Op31) than from the LdMNPV (Ld116) homologue. The locations of the dutpase gene on the genomes of LdMNPV and SeMNPV are quite similar, whereas the OpMNPV homologue is located in a different region on the genome (Fig. 2C, D
). This suggests that the SeMNPV and LdMNPV dutpase genes were acquired from the same source, whereas the OpMNPV dutpase may have been acquired independently from a different source. The different location of OpMNPV dutpase could also have resulted from gene duplication and rearrangement.
Genes with auxiliary functions
The auxiliary genes (OReilly, 1997 ) superoxide dismutase (sod) (identity to AcMNPV homologue 67%), chitinase (65%), cathepsin (55%) and ecdysteroid UDP-glucosyltransferase (50%) were quite well conserved, whereas the fibroblast growth factor (fgf) gene (identity to AcMNPV homologue 24%) was quite different from the other baculovirus fgf genes.
Ac1 encodes a protein tyrosine/serine phosphatase with dual specificity (dsPTP) (Tilakaratne et al., 1991 ; Kim & Weaver, 1993
). This protein removes phosphates specifically from both tyrosine and serine/threonine residues and regulates the phosphorylation status of a variety of proteins, including growth factors, which in turn regulate developmental processes in living cells (Wishart et al., 1995
). The absence of a ptp-1 homologue in SeMNPV may not necessarily result in loss of PTPase function, however, because a ptp-2 (Op9) homologue was present in the SeMNPV genome. The SeMNPV PTP-2 homologue, Se26, contained the conserved domain [HCXXGXXR(S/T)] encoding the dsPTP catalytic loop. It is therefore likely that the PTP-2 homologue encodes an active tyrosine/serine phosphatase.
A protein kinase enzyme activity also appeared to have been retained in the SeMNPV genome, since a pk-1 homologue of AcMNPV was present (Se3). However, a pk-2 homologue was absent. pk-2 was shown to be non-essential for AcMNPV, since a pk-2 deletion mutant had no detectable effect on AcMNPV replication in cell cultures (Chen & Thiem, 1997 ). Although AcMNPV pk-2 is non-essential, its presence favours virus gene expression by inhibiting a host stress response in infected cells (Dever et al., 1998
).
A homologue of the actin rearrangement-inducing factor-1 (arif-1) was present in the SeMNPV genome. This arif-1 gene induces rearrangements of the actin cytoskeleton after infection, but the functional significance of these conformational changes remains to be elucidated (Roncarati & Knebel-Mörsdorf, 1997 ).
No homologues of the LdMNPV viral enhancing factors (vef)-1 and -2 have been identified in SeMNPV (Bischoff & Slavicek et al., 1997 ; Hashimoto et al., 1991
). In GVs, the vef gene products increase virus potency by disrupting the peritrophic membrane, thereby allowing virions access to the surface of midgut epithelial cells (Wang & Granados, 1998
; Derksen & Granados, 1988
). The vef genes encode metalloproteases that specifically degrade the mucin protein component of the peritrophic membrane (Wang & Granados, 1997
; Lepore et al., 1996
). vef homologues are absent in group I baculoviruses and may, therefore, be unique to LdMNPV among the NPVs.
SeMNPV homologue ORFs of unknown function: two p26 homologues
ORFs without assigned functions, but well conserved among the four baculoviruses, include Ac106/107 (identity to AcMNPV homologue 58 %), Ac38 (58%), Ac22 (57%), Ac92 (55%) and Ac103 (54%). The high percentage identities between the baculovirus homologues suggest that these ORFs have essential functions in virus multiplication and pathology, for which a certain degree of conservation is required.
Notable is the pairwise conservation of Ac76 between AcMNPV and OpMNPV (81%) and between LdMNPV and SeMNPV (71%), which may suggest that Ac76 homologues have been acquired twice during evolution from two different sources. Other pairwise alignments yielded identities no higher than 45%. To this end, it can be speculated that SeMNPV and LdMNPV have a more recent baculovirus ancestor in common than SeMNPV and AcMNPV or OpMNPV.
Some ORFs that were previously unique to LdMNPV have homologues in the SeMNPV genome (SeMNPV ORFs 15, 28, 30, 33, 49, 51, 52 and 107) (Table 1). The SeMNPV genome also contained a homologue of the previously described LdMNPV ORF4 (Bjornson & Rohrmann, 1992
), although this ORF was not included in the LdMNPV genome analysis (Kuzio et al., 1999
) due to overlap with Ld137. A similar situation occurs for the LdMNPV homologue of Se37 that overlaps with Ld155.
gp37 (Se25), named spindle-like protein or fusolin because of its obvious homology to the entomopoxvirus spindle-shaped proteins, is a conserved NPV gene (identity to AcMNPV homologue 56%; Liu & Carstens, 1996 ). The gp37/fusolin gene family may be essential for virus replication, based on the failure to construct an insertion mutant for this gene in AcMNPV (Wu & Miller, 1989
). Furthermore, studies have suggested that gp37/fusolin is involved in enhancement of virus infection in vivo (Phanis et al., 1999
).
Unlike any baculovirus genome so far analysed, SeMNPV possessed two copies of p26 (Se87 and 129). Se87 was located in the proximity of the non-hr (Fig. 1). This region is organized differently compared with AcMNPV, in contrast to the Se129 region (Fig. 2
; see position of p26 in all GeneParityPlots). It is possible that Se87 was acquired independently from a different source than Se129. This view is further supported by the 8% less identity of Se87 to its AcMNPV homologue than Se129. It is equally possible that Se129 has diverged from Se87 and has been rearranged following duplication. Transcripts have been identified for the AcMNPV homologue, which are synthesized by the host polymerase II both early and late in infection (Huh & Weaver, 1990
). The P26 protein was localized primarily to the cytoplasm and is present in the membrane fraction of BV (Goenka & Weaver, 1996
). A function is not yet been assigned to the P26 protein, but its conservation in all MNPV genomes analysed so far suggests a function basic to baculoviruses.
Baculovirus repeated ORFs (bro genes)
bro genes, present in a number of other baculoviruses and to date of no known function, were not identified in SeMNPV. Five copies of a homologue of Ac2 were identified in BmNPV (Gomi et al., 1999 ) and there were 16 copies in LdMNPV (Kuzio et al., 1999
). In OpMNPV, a truncated version and two smaller bro-related ORFs are present (Ahrens et al., 1997
). Similarity searches revealed that Se13 showed weak homology (~25%) to some bro genes, particularly to BmNPV bro-d (Gomi et al., 1999
) and LdMNPV bro-j (Kuzio et al., 1999
) (Table 5
). However, Se13 had higher homology (~33%) to Ac13 and its homologues Bm5, Op12 and Ld122 (Table 1
). Furthermore, Se13 is located adjacent to a homologue of Ac14 (Se14) and these two genes are clustered in all baculoviruses compared. Therefore, we consider that Se13 is not a bro gene sensu stricto.
|
The SeMNPV genome contained an ORF (Se4) with homology to the hoar ORF of HzSNPV and Helicoverpa armigera NPV (Le et al., 1997 ). In SeMNPV and HzSNPV, the pk-1 gene is downstream of the hoar gene (Table 1
; Le et al., 1997
). However, the upstream flanking ORFs of the HzSNPV hoar gene had no homologues in SeMNPV (HzSNPV ORF480 and ORF321) or were present in different locations in SeMNPV (HzSNPV ORF16 corresponding to Se138132) (Le et al., 1997
; Table 1
). The SeMNPV hoar ORF, like the Heliothis sp. homologues, contained a complex, A+T-rich, triplet repeat region (RAT-repeats) distributed over 330 bp and a C3HC4 (RING finger) zinc-binding motif.
Unique SeMNPV ORFs
Twenty ORFs in the SeMNPV genome were completely unique to this virus and did not exhibit significant homology to any sequence in the GenBank. Hence no putative functions could be assigned to these ORFs. The functions of these ORFs are being investigated. This number is roughly proportional to the size of the genome (Table 2).
Organization of the SeMNPV genome
The genomic organization, i.e. the order of genes, is similar in AcMNPV, BmNPV and OpMNPV, except for a small number of rearrangements (Ahrens et al., 1997 ; Hu et al., 1998
; Gomi et al., 1999
). To investigate whether the organization in SeMNPV was collinear with these viruses and to the recently sequenced LdMNPV (Kuzio et al., 1999
), a comparison was made between the SeMNPV genome organization and those of AcMNPV, BmNPV, OpMNPV and LdMNPV by using GeneParityPlot analysis (Fig. 2
; Hu et al., 1998
). The gene organization was most conserved in the central region (ORFs 3070) of the linearized baculovirus genomes, confirming the assumption of Heldens et al. (1997b)
. The left part of the SeMNPV genome (ORFs 130) displayed a considerable number of gene inversions and translocations in the GeneParityPlot analyses. The right part (ORFs 70100) showed a high degree of gene scrambling (Fig. 2
). From these analyses, it is concluded that the organization of SeMNPV is highly characteristic and distinct from those of AcMNPV, BmNPV, OpMNPV and LdMNPV.
By convention, the orientation of a circular baculovirus genome is determined by the relative position of two genes, polyhedrin at map unit 0 and p10 approximately at map unit 90 (Vlak & Smith, 1982 ). In the initial GeneParityPlot analysis, the orientation of the SeMNPV genome appeared to be inverted for more than 50% of the ORFs compared with AcMNPV, BmNPV, OpMNPV and LdMNPV. This led to perpendicularity in the graph where collinearity was known to occur, i.e. in the conserved, central part of the genome. To facilitate convenient comparison and interpretation of the different genomes, the SeMNPV gene order was reversed before it was subjected to GeneParityPlot analysis. The previously satisfactory choices of polyhedrin and p10 for the zero point and directional orientation, respectively, were not convenient for GeneParityPlot analysis in this case because both genes are located in regions that show extensive rearrangements (Fig. 2
).
Comparison of the relative gene order between SeMNPV and AcMNPV, BmNPV, OpMNPV and LdMNPV revealed the presence of certain gene clusters that are conserved in all baculovirus genomes compared. These clusters were numbered according to their sequential appearance in the GeneParityPlots. Fourteen clusters conserved in all five baculoviruses were identified (Fig. 2, Table 1
). Cluster 3 was interrupted in SeMNPV by the insertion of Se17 and Se18, which are unique to SeMNPV. Cluster 12 is discontinuous in LdMNPV because five copies of the bro gene and two other ORFs are inserted in this cluster. Four additional clusters were identified in the GeneParityPlot of SeMNPV versus LdMNPV (Fig. 2D
; Table 1
; clusters 1518). Furthermore, clusters 2 and 5 were extended to include genes Se15 and Se38+Se41, respectively. Clusters 7 and 8 and clusters 9 and 10 were present as two contiguous clusters in LdMNPV and SeMNPV. This is in contrast to the other three baculoviruses, where the positions of genes of these clusters in the gene parity plot were perpendicular to each other due to inversion of one of the clusters (Ayres et al., 1994
; Gomi et al., 1999
; Ahrens et al., 1997
, Kuzio et al., 1999
). The additional and the enlarged clusters of SeMNPV and LdMNPV suggest that the genomic organization of SeMNPV is more closely related to that of LdMNPV than to that of AcMNPV, BmNPV and OpMNPV. This agrees with the phylogenetic analysis of single genes such as egt, lef2 and rr1, which shows that SeMNPV is more closely related to LdMNPV than to AcMNPV, BmNPV or OpMNPV (Chen et al., 1997
, 1999
; Hu et al., 1997
; van Strien et al., 1997
). Thus, juxtaposition of ORFs can be used as a phylogenetic marker to study the ancestral relationship of baculoviruses, independent of the evolution of individual genes.
Between- and within-baculovirus genome rearrangement
Comparison of SeMNPV with AcMNPV, BmNPV, OpMNPV and LdMNPV showed that baculovirus genomes may vary due to deletions, (gene) insertions, inversions and duplications (Ayres et al., 1994 ; Gomi et al., 1999
; Ahrens et al., 1997
; Kuzio et al., 1999
). The mechanisms underlying these rearrangements are still unclear. Transposable elements that may play a role in rearrangements of baculovirus genomes have been identified in several baculoviruses (Friesen, 1993
; Jehle, 1996
; Jehle et al., 1997
). Furthermore, there is evidence to suggest that hr sequences are related to the generation of variant baculovirus genotypes (Muñoz et al., 1999
).
Genome rearrangements also occur within one baculovirus species, as is the case for SeMNPV. A mutant SeMNPV, containing a single deletion of approximately 25 kb, was obtained within the first passage in cell culture (Heldens et al., 1996 ). This deletion is located approximately between 17·5 and 42·0 kb (±1·0 kb) and encompasses Se15 to 41 (Table 1
). So far, none of these ORFs has been shown to be essential for virus replication. However, deletion mutant SeMNPV polyhedra produced in vitro do not cause any pathological effect in vivo nor does the injection of BV into the haemolymph. In contrast, mutant SeMNPV BV was highly infectious for Se-UCR1 cells and resulted in polyhedron production (Heldens et al., 1996
). Therefore, at least one gene located in the deleted sequences contains information that is important for virulence in vivo.
In conclusion, sequencing revealed that the genome of SeMNPV is distinct from those of other baculoviruses both in gene content and arrangement. Two, probably independently acquired, p26 and odv-e66 genes are present. Notably, SeMNPV lacks homologues of the gp64, ie-2 and multiple bro genes. Furthermore, SeMNPV and LdMNPV may have a recent common ancestor, whereas they are more distantly related to AcMNPV, BmNPV and OpMNPV on the basis of gene homology and genomic organization. The gene order in the central part of baculovirus genomes is highly conserved, whereas the gene order in the other segments has been subjected to multiple rearrangements. The GeneParityPlot analyses demonstrate that this method can be used as an independent means of phylogenetic study and can provide an initial view of the conservation of gene clusters and how viruses may have obtained additional genes. The genomic sequences absent in the deletion mutant of SeMNPV contain information that is important for virulence in vivo. Further studies will concentrate on the functional analysis of the ORFs that are unique to SeMNPV. These studies will provide insight in the roles these ORFs may play in the high virulence and narrow host-range of SeMNPV.
![]() |
Acknowledgments |
---|
![]() |
Footnotes |
---|
b Present address: Fort Dodge Animal Health Holland, Biological Research and Development, C. J. van Houtenlaan 36, 1381 CP Weesp, The Netherlands.
c Present address: Department of Virology, Leiden University, 2300 RC Leiden, The Netherlands.
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() |
---|
Ahrens, C. H. & Rohrmann, G. F. (1995). Replication of Orgyia pseudotsugata baculovirus DNA: lef-2 and ie-1 are essential and ie-2, p34, and Op-iap are stimulatory genes. Virology 212, 650-662.[Medline]
Ahrens, C. H., Russell, R. L., Funk, C. J., Evans, J. T., Harwood, S. H. & Rohrmann, G. F. (1997). The sequence of the Orgyia pseudotsugata multinucleocapsid nuclear polyhedrosis virus genome. Virology 229, 381-399.[Medline]
Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. (1990). Basic local alignment search tool. Journal of Molecular Biology 215, 403-410.[Medline]
Ayres, M. D., Howard, S. C., Kuzio, J., Lopez-Ferber, M. & Possee, R. D. (1994). The complete DNA sequence of Autographa californica nuclear polyhedrosis virus. Virology 202, 586-605.[Medline]
Bairoch, A., Bucher, P. & Hofmann, K. (1997). The PROSITE database, its status in 1997. Nucleic Acids Research 25, 217-221.
Birnbaum, M. J., Clem, R. J. & Miller, L. K. (1994). An apoptosis-inhibiting gene from a nuclear polyhedrosis virus encoding a polypeptide with Cys/His sequence motifs. Journal of Virology 68, 2521-2528.[Abstract]
Bischoff, D. S. & Slavicek, J. M. (1997). Molecular analysis of an enhancin gene in the Lymantria dispar nuclear polyhedrosis virus. Journal of Virology 71, 8133-8140.[Abstract]
Bjornson, R. M. & Rohrmann, G. F. (1992). Nucleotide sequence of the polyhedron envelope protein gene region of the Lymantria dispar nuclear polyhedrosis virus. Journal of General Virology 73, 1499-1504.[Abstract]
Broer, R., Heldens, J. G. M., van Strien, E. A., Zuidema, D. & Vlak, J. M. (1998). Specificity of multiple homologous genomic regions in Spodoptera exigua nucleopolyhedrovirus DNA replication. Journal of General Virology 79, 1563-1572.[Abstract]
Chen, C. J. & Thiem, S. M. (1997). Differential infectivity of two Autographa californica nucleopolyhedrovirus mutants on three permissive cell lines is the result of lef-7 deletion. Virology 227, 88-95.[Medline]
Chen, X., Hu, Z., Jehle, J. A., Zhang, Y. & Vlak, J. M. (1997). Analysis of the ecdysteroid UDP-glucosyltransferase gene of Heliothis armigera single-nucleocapsid baculovirus. Virus Genes 15, 219-225.[Medline]
Chen, X., IJkel, W. F. J., Dominy, C., Zanotto, P. M. de A., Hashimoto, Y., Faktor, O., Hayakawa, T., Wang, C. H., Premkumar, A., Mathavan, S., Krell, P. J., Hu, Z. & Vlak, J. M. (1999). Identification, sequence and phylogeny of the lef-2 gene of Helicoverpa armigera single-nucleocapsid baculovirus. Virus Research (in press).
Clem, R. J. & Miller, L. K. (1994). Control of programmed cell death by the baculovirus genes p35 and iap. Molecular and Cellular Biology 14, 5212-5222.[Abstract]
Clem, R. J., Robson, M. & Miller, L. K. (1994). Influence of infection route on the infectivity of baculovirus mutants lacking the apoptosis-inhibiting gene p35 and the adjacent gene p94. Journal of Virology 68, 6759-6762.[Abstract]
Crook, N. E., Clem, R. J. & Miller, L. K. (1993). An apoptosis-inhibiting baculovirus gene with a zinc finger-like motif. Journal of Virology 67, 2168-2174.[Abstract]
Derksen, A. C. & Granados, R. R. (1988). Alteration of a lepidopteran peritrophic membrane by baculoviruses and enhancement of viral infectivity. Virology 167, 242-250.[Medline]
Dever, T. E., Sripriya, R., McLachlin, J. R., Lu, J., Fabian, J. R., Kimball, S. R. & Miller, L. K. (1998). Disruption of cellular translational control by a viral truncated eukaryotic translation initiation factor 2 kinase homolog. Proceedings of the National Academy of Sciences, USA 95, 4164-4169.
Devereux, J., Haeberli, P. & Smithies, O. (1984). A comprehensive set of sequence analysis programs for the VAX. Nucleic Acids Research 12, 387-395.[Abstract]
Elledge, S. J., Zhou, Z. & Allen, J. B. (1992). Ribonucleotide reductase: regulation, regulation, regulation. Trends in Biochemical Sciences 17, 119-123.[Medline]
Fabian, P., Murvai, J., Hatsagi, Z., Vlahovicek, K., Hegyi, H. & Pongor, S. (1997). The SBASE protein domain library, release 5.0: a collection of annotated protein sequence segments. Nucleic Acids Research 25, 240-243.
Friesen, P. D. (1993). Invertebrate transposable elements in the baculovirus genome: characterization and significance. In Parasites and Pathogens of Insects: Parasites, pp. 147-178. Edited by N. E. Beckage, S. N. Thompson & B. A. Federici. London: Academic Press.
Friesen, P. D. & Miller, L. K. (1987). Divergent transcription of early 35- and 94-kilodalton protein genes encoded by the HindIII K genome fragment of the baculovirus Autographa californica nuclear polyhedrosis virus. Journal of Virology 61, 2264-2272.[Medline]
Gelernter, W. D. & Federici, B. A. (1986). Isolation, identification and determination of virulence of a nuclear polyhedrosis virus from the beet armyworm, Spodoptera exigua (Lepidoptera: Noctuidae). Environmental Entomology 15, 240-245.
Goenka, S. & Weaver, R. F. (1996). Cellular localization of p26, an early AcMNPV gene product. In Abstracts of the 15th Annual Meeting of the American Society of Virology, p. 154.
Gomi, S., Zhou, C. E., Yih, W., Majima, K. & Maeda, S. (1997). Deletion analysis of four of eighteen late gene expression factor gene homologues of the baculovirus, BmNPV. Virology 230, 35-47.[Medline]
Gomi, S., Majima, K. & Maeda, S. (1999). Sequence analysis of the genome of Bombyx mori nucleopolyhedrovirus. Journal of General Virology 80, 1323-1337.[Abstract]
Griffiths, C. M., Barnett, A. L., Ayres, M. D., Windass, J., King, L. A. & Possee, R. D. (1999). In vitro host range of Autographa californica nucleopolyhedrovirus recombinants lacking functional p35, iap1 or iap2. Journal of General Virology 80, 1055-1066.[Abstract]
Hashimoto, Y., Corsaro, B. G. & Granados, R. R. (1991). Location and nucleotide sequence of the gene encoding the viral enhancing factor of the Trichoplusia ni granulosis virus. Journal of General Virology 72, 2645-2651.[Abstract]
Heldens, J. G. M., van Strien, E. A., Feldmann, A. M., Kulcsár, P., Munoz, D., Leisy, D. J., Zuidema, D., Goldbach, R. W. & Vlak, J. M. (1996). Spodoptera exigua multicapsid nucleopolyhedrovirus deletion mutants generated in cell culture lack virulence in vivo. Journal of General Virology 77, 3127-3134.[Abstract]
Heldens, J. G. M., Broer, R., Zuidema, D., Goldbach, R. W. & Vlak, J. M. (1997a). Identification and functional analysis of a non-hr origin of DNA replication in the genome of Spodoptera exigua multicapsid nucleopolyhedrovirus. Journal of General Virology 78, 1497-1506.[Abstract]
Heldens, J. G. M., Liu, Y., Zuidema, D., Goldbach, R. W. & Vlak, J. M. (1997b). Characterization of a putative Spodoptera exigua multicapsid nucleopolyhedrovirus helicase gene. Journal of General Virology 78, 3101-3114.[Abstract]
Hong, T., Braunagel, S. C. & Summers, M. D. (1994). Transcription, translation, and cellular localization of PDV-E66: a structural protein of the PDV envelope of Autographa californica nuclear polyhedrosis virus. Virology 204, 210-222.[Medline]
Hu, Z. H., Broer, R., Westerlaken, J., Martens, J. W. M., Jin, F., Jehle, J. A., Wang, L. M. & Vlak, J. M. (1997). Characterization of the ecdysteroid UDP-glucosyltransferase gene of a single nucleocapsid nucleopolyhedrovirus of Buzura suppressaria. Virus Research 47, 91-97.[Medline]
Hu, Z. H., Arif, B. M., Jin, F., Martens, J. W. M., Chen, X. W., Sun, J. S., Zuidema, D., Goldbach, R. W. & Vlak, J. M. (1998). Distinct gene arrangement in the Buzura suppressaria single-nucleocapsid nucleopolyhedrovirus genome. Journal of General Virology 79, 2841-2851.[Abstract]
Huh, N. E. & Weaver, R. F. (1990). Identifying the RNA polymerases that synthesize specific transcripts of the Autographa californica nuclear polyhedrosis virus. Journal of General Virology 71, 195-201.[Abstract]
Jehle, J. A. (1996). Transmission of insect transposons into baculovirus genomes: an unusual hostpathogen interaction. In Transgenic Organisms Biological and Social Implications, pp. 81-97. Edited by J. Tomiuk, K. Wohrmann & A. Sentker. Basel: Birlhauser Verlag.
Jehle, J. A., van der Linden, I. F., Backhaus, H. & Vlak, J. M. (1997). Identification and sequence analysis of the integration site of transposon TCp3.2 in the genome of Cydia pomonella granulovirus. Virus Research 50, 151-157.[Medline]
Kim, D. & Weaver, R. F. (1993). Transcription mapping and functional analysis of the protein tyrosine/serine phosphatase (PTPase) gene of the Autographa californica nuclear polyhedrosis virus. Virology 195, 587-595.[Medline]
Kuzio, J., Pearson, M. N., Harwood, S. H., Funk, C. J., Evans, J. T., Slavicek, J. M. & Rohrmann, G. F. (1999). Sequence and analysis of the genome of a baculovirus pathogenic for Lymantria dispar. Virology 253, 17-34.[Medline]
Le, T. H., Wu, T., Robertson, A., Bulach, D., Cowan, P., Goodge, K. & Tribe, D. (1997). Genetically variable triplet repeats in a RING-finger ORF of Helicoverpa species baculoviruses. Virus Research 49, 67-77.[Medline]
Lepore, L. S., Roelvink, P. R. & Granados, R. R. (1996). Enhancin, the granulosis virus protein that facilitates nucleopolyhedrovirus (NPV) infections, is a metalloprotease. Journal of Invertebrate Pathology 68, 131-140.[Medline]
Liu, J. J. & Carstens, E. B. (1996). Identification, molecular cloning, and transcription analysis of the Choristoneura fumiferana nuclear polyhedrosis virus spindle-like protein gene. Virology 223, 396-400.[Medline]
Lu, A. & Miller, L. K. (1995). The roles of eighteen baculovirus late expression factor genes in transcription and DNA replication. Journal of Virology 69, 975-982.[Abstract]
Mikhailov, V. S., Mikhailova, A. L., Iwanaga, M., Gomi, S. & Maeda, S. (1998). Bombyx mori nucleopolyhedrovirus encodes a DNA-binding protein capable of destabilizing duplex DNA. Journal of Virology 72, 3107-3116.
Muñoz, D., Castillejo, J. I. & Caballero, P. (1998). Naturally occurring deletion mutants are parasitic genotypes in a wild-type nucleopolyhedrovirus population of Spodoptera exigua. Applied and Environmental Microbiology 64, 4372-4377.
Muñoz, D., Murillo, R., Krell, P. J., Vlak, J. M. & Caballero, P. (1999). Four genotypic variants of a Spodoptera exigua Nucleopolyhedrovirus (Se-SP2) are distinguishable by a hypervariable genomic region. Virus Research 59, 61-74.[Medline]
Nielsen, H., Engelbrecht, J., Brunak, S. & von Heijne, G. (1997). Identification of prokaryotic and eukaryotic signal peptides and prediction of their cleavage sites. Protein Engineering 10, 1-6.[Abstract]
Okano, K., Mikhailov, V. S. & Maeda, S. (1999). Colocalization of baculovirus IE-1 and two DNA-binding proteins, DBP and LEF-3, to viral replication factories. Journal of Virology 73, 110-119.
Oomens, A. G. P. & Blissard, G. W. (1999). Requirement for GP64 to drive efficient budding of Autographa californica multicapsid nucleopolyhedrovirus. Virology 254, 297-314.[Medline]
OReilly, D. R. (1997). Auxiliary genes of baculoviruses. In The Baculoviruses, pp. 267-295. Edited by L. K. Miller. New York: Plenum Press.
Pearson, W. R. (1990). Rapid and sensitive sequence comparison with FASTP and FASTA. Methods in Enzymology 183, 63-98.[Medline]
Phanis, C. G., Miller, D. P., Cassar, S. C., Tristem, M., Thiem, S. M. & OReilly, D. R. (1999). Identification and expression of two baculovirus gp37 genes. Journal of General Virology 80, 1823-1831.[Abstract]
Ranjan, A. & Hasnain, S. E. (1995). Codon usage in the prototype baculovirus Autographa californica nuclear polyhedrosis virus. Indian Journal of Biochemistry and Biophysics 32, 424-428.[Medline]
Roncarati, R. & Knebel-Mörsdorf, D. (1997). Identification of the early actin-rearrangement-inducing factor gene, arif-1, from Autographa californica multicapsid nuclear polyhedrosis virus. Journal of Virology 71, 7933-7941.[Abstract]
Rost, B. & Sander, C. (1993). Prediction of protein secondary structure at better than 70% accuracy. Journal of Molecular Biology 232, 584-599.[Medline]
Seshagiri, S. & Miller, L. K. (1997). Baculovirus inhibitors of apoptosis (IAPs) block activation of Sf-caspase-1. Proceedings of the National Academy of Sciences, USA 94, 13606-13611.
Smits, P. H. (1987). Nuclear polyhedrosis virus as biological control agent of Spodoptera exigua. PhD thesis. Wageningen Agricultural University, Wageningen, The Netherlands.
Smits, P. H., van de Vrie, M. & Vlak, J. M. (1988). Nuclear polyhedrosis virus control of Spodoptera exigua on glasshouse crops. Entomologia Experimentalis et Applicata 43, 73-80.
Tilakaratne, N., Hardin, S. E. & Weaver, R. F. (1991). Nucleotide sequence and transcript mapping of the HindIII F region of the Autographa californica nuclear polyhedrosis virus genome. Journal of General Virology 72, 285-291.[Abstract]
van Strien, E. A., Zuidema, D., Goldbach, R. W. & Vlak, J. M. (1992). Nucleotide sequence and transcriptional analysis of the polyhedrin gene of Spodoptera exigua nuclear polyhedrosis virus. Journal of General Virology 73, 2813-2821.[Abstract]
van Strien, E. A., Jansen, B. J. H., Mans, R. M. W., Zuidema, D. & Vlak, J. M. (1996). Sequence and transcriptional analysis of the ubiquitin gene cluster in the genome of Spodoptera exigua nucleopolyhedrovirus. Journal of General Virology 77, 2311-2319.[Abstract]
van Strien, E. A., Faktor, O., Hu, Z. H., Zuidema, D., Goldbach, R. W. & Vlak, J. M. (1997). Baculoviruses contain a gene for the large subunit of ribonucleotide reductase. Journal of General Virology 78, 2365-2377.[Abstract]
Vlak, J. M. & Smith, G. E. (1982). Orientation of the genome of Autographa californica nuclear polyhedrosis virus: a proposal. Journal of Virology 41, 1118-1121.
Volkman, L. E., Blissard, G. W., Friesen, P., Keddie, B. A., Possee, R. & Theilmann, D. A. (1995). Family Baculoviridae. In Virus Taxonomy. Sixth Report of the International Committee on Taxonomy of Viruses, pp. 104-113. Edited by F. A. Murphy, C. M. Fauquet, D. H. L. Bishop, S. A. Ghabrial, A. W. Jarvis, G. P. Martelli, M. A. Mayo & M. D. Summers. Vienna & New York: Springer-Verlag.
Vucic, D., Kaiser, W. J. & Miller, L. K. (1998). A mutational analysis of the baculovirus inhibitor of apoptosis Op-IAP. Journal of Biological Chemistry 273, 33915-33921.
Wang, P. & Granados, R. R. (1997). An intestinal mucin is the target substrate for a baculovirus enhancin. Proceedings of the National Academy of Sciences, USA 94, 6977-6982.
Wang, P. & Granados, R. R. (1998). Observations on the presence of the peritrophic membrane in larval Trichoplusia ni and its role in limiting baculovirus infection. Journal of Invertebrate Pathology 72, 57-62.[Medline]
Wang, J. W., Qi, Y. P., Huang, Y. X. & Li, S. D. (1995). Nucleotide sequence of a 1446 base pair SalI fragment and structure of a novel early gene of Leucania separata nuclear polyhedrosis virus. Archives of Virology 140, 2283-2291.[Medline]
Wishart, M. J., Denu, J. M., Williams, J. A. & Dixon, J. E. (1995). A single mutation converts a novel phosphotyrosine binding domain into a dual-specificity phosphatase. Journal of Biological Chemistry 270, 26782-26785.
Wu, J. & Miller, L. K. (1989). Sequence, transcription and translation of a late gene of the Autographa californica nuclear polyhedrosis virus encoding a 34·8K polypeptide. Journal of General Virology 70, 2449-2459.[Abstract]
Zanotto, P. M., Kessing, B. D. & Maruniak, J. E. (1993). Phylogenetic interrelationships among baculoviruses: evolutionary rates and host associations. Journal of Invertebrate Pathology 62, 147-164.[Medline]
Zuidema, D., van Oers, M. M., van Strien, E. A., Caballero, P. C., Klok, E.-J., Goldbach, R. W. & Vlak, J. M. (1993). Nucleotide sequence and transcriptional analysis of the p10 gene of Spodoptera exigua nuclear polyhedrosis virus. Journal of General Virology 74, 1017-1024.[Abstract]
Received 30 June 1999;
accepted 9 September 1999.