Previous Article | Next Article ![]()
Journal of Virology, January 2004, p. 168-177, Vol. 78, No. 1
0022-538X/04/$08.00+0 DOI: 10.1128/JVI.78.1.168-177.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Plum Island Animal Disease Center, Agricultural Research Service, United States Department of Agriculture, Greenport, New York 11944,1 Area of Virology, School of Veterinary Sciences, University of Buenos Aires, 1427 Buenos Aires, Argentina,2 Department of Veterinary Pathobiology, College of Veterinary Medicine, Texas A&M University, College Station, Texas 77843-4467,3 National Animal Disease Center, Agricultural Research Service, United States Department of Agriculture, Ames, Iowa 500104
Received 19 August 2003/ Accepted 22 September 2003
|
|
|---|
|
|
|---|
PPVs cause nonsystemic, eruptive skin disease in domestic and wild mammals. ORFV, the prototype species of PPV, is responsible for contagious ecthyma, an acute disease of sheep and goats. The disease, also known as orf, contagious pustular dermatitis, or scabby mouth, is characterized by proliferative lesions in the skin of the lips, around the nostrils, and in the oral mucosa (27). Lesions progress through a typical pattern of erythema, papula, pustule, and scab and usually resolve in 1 to 2 months (45). Although considered a mild disease, mortality rates up to 93% have been reported in kids (41). High mortality rates occur when lesions in lips and udders prevent infected animals from suckling and grazing, resulting in rapid emaciation (13, 41, 58). Sheep can be repeatedly infected with ORFV, albeit with shorter times to recovery and less pronounced pathological changes than in a primary infection (45). A Th1-type immune response has been implicated in host immunity to ORFV infection (reviewed in reference 32). Attenuated orf vaccines can limit the severity of the infection but they are unable to completely prevent the disease (30).
BPSV infects cattle of all ages but clinical signs are usually seen in calves. The disease has a worldwide distribution and is characterized by papules, often mildly erosive, on the muzzle, oral mucosa, and udder and occasionally in the esophagus and forestomach (40). Like ORFV in sheep and goats, reinfection of cattle with BPSV is commonly observed, suggesting that virus infection does not confer significant immunity. Because of its clinical resemblance to foot-and-mouth disease, BPSV infections are reported to veterinary authorities for differential diagnosis.
Cocirculation of BPSV and ORFV in wild ruminants has been described (35), and PPV isolates from wild ruminants have been experimentally transmitted to sheep, goats, and cattle (59, 60). Both ORFV and BPSV cause occupational infections in humans with lesions characterized by large, painful nodules in the hands and, less frequently, the face (8, 47, 69).
Classification of PPVs has relied on natural host range, pathology, and viral DNA restriction enzyme analysis. The latter revealed considerable genomic heterogeneity, even between isolates of the same virus (26, 35, 63, 64). Hybridization analysis of viral DNA indicates that internal but not terminal genomic regions are conserved between PPVs (26). Data concerning PPV genomics, largely obtained by using ORFV strain NZ2, has revealed (i) colinearity between the ORFV and other poxvirus genera genomes (21, 49, 50), (ii) the presence of a set of genes mostly located at the termini of the viral genome with putative or known roles in virulence or immunomodulation (15, 23, 38, 42, 76), and (iii) the occurrence of genomic rearrangements of the termini upon serial virus passage in vitro (12, 22). Less is known about the gene complement and genomic organization of other PPV. DNA sequencing of the right end of the BPSV strain B177 genome indicated an organization similar to that of the right end of the ORFV genome, except for the lack of a vascular endothelial growth factor (VEGF) gene in BPSV (67). Here we present the complete DNA sequences of two ORFV isolates and one BPSV isolate, thus providing an overview of PPV genomic organization and gene content as well as a comparison between the two viruses.
|
|
|---|
Viral DNA isolation, cloning, sequencing, and sequence analysis. Viral genomic DNA was extracted from infected primary lamb kidney cell cultures as previously described (82). Random DNA fragments were obtained by incomplete enzymatic double digestion with AciI and TaqI endonucleases (New England Biolabs, Beverly, Mass.), and DNA fragments larger than 1.0 kbp were cloned and used in dideoxy sequencing reactions as previously described (2). Reaction products were analyzed on an ABI Prism 3700 automated DNA sequencer (Applied Biosystems, Foster City, Calif.). Sequence data were assembled with the Phrap and CAP3 software programs (18, 19, 33), and gaps were closed as described previously (1). The final DNA consensus sequences for each genome represented on average seven- to ninefold redundancy at each base position and a Consed estimated error rate of 0.01 per 10 kbp (18, 19, 28).
Genome DNA composition, structure, repeats, and restriction enzyme patterns were analyzed as previously described (1) with the Genetics Computer Group GCG version 10 software package (16). Pairwise genomic alignments were done with WABA (http://www.cse.ucsc.edu/~kent/), and multiple genomic alignments were done with Dialign (54) and Clustal (77) alignment programs. Open reading frames (ORFs) longer than 30 codons were evaluated for coding potential and ORFs greater than 60 codons were subjected to homology searches as previously described (1, 2). In addition, Framefinder (http://www.hgmp.mrc.ac.uk/~gslater) was used to evaluate coding potential. Based on these criteria, 131 (BPSV) and 130 (ORFV) putative genes were annotated and orthologous ORFs were similarly numbered. Phylogenetic comparisons were done with the PHYLO_WIN software package (25) and Puzzle (75).
Nucleotide sequence accession number. The genome sequences of ORFV strains IA82 and SA00 and BPSV strain AR02 have been deposited in GenBank under accession no. AY386263, AY386264, and AY386265, respectively.
|
|
|---|
Like other poxviruses, BPSV and ORFV genomes contain a large central coding region bounded by two identical inverted terminal repeat (ITR) regions (12, 26, 48). Assembled ITRs of BV-AR02, OV-SA00, and OV-IA82 contain 1,161, 3,936, and 3,092 bp, respectively. The differences in length between the ITRs of OV-SA00 and OV-IA82 strains are in agreement with previous work, indicating natural intrastrain variations in this genomic region (64). Only one ORF (001), previously described for ORFV strains NZ2 and NZ7 (20, 24) and in BPSV strain B177 (67), initiates and is completely located within the ITRs in the three virus isolates. This ORF of unknown function is unique to BPSV and ORFV, sharing 63% amino acid identity. Putative transcription control elements similar to those described for the ORFV strain NZ2 homologue are found flanking BPSV 001, suggesting early gene expression, as is the case for ORFV NZ2 (20). A second ORFV gene of unknown function (002), not present in BPSV, initiates within the unique region and terminates within the ITR.
Despite the high G+C content and paucity of stop codons, which yield 362 and 345 methionine-initiated ORFs of at least 60 codons in BPSV and ORFV genomes, respectively, coding potential analysis and similarity to known proteins led us to conservatively predict 131 genes in BPSV and 130 genes in ORFV. These genes, which encode proteins of 53 to 1289 amino acids, represent an approximate coding density of 90% for BPSV and 95% for ORFV (Table 1). The central genomic core region (ORFs 009 to 111) contains homologues of conserved poxvirus genes involved in basic replicative mechanisms and structure and morphogenesis of intracellular mature and extracellular enveloped virions (EEV) (55) (Table 1). Homologues of vaccinia virus (VACV) F9L and F10L, which are located at the left end of the conserved core in most ChPVs, are located at the right end of PPV genomes (ORFs 130 and 131). Terminal genomic regions (ORFs 001 to 008 and 112 to 134) represent approximately 20% of the viral genome and contain genes likely affecting pathogenesis. These include genes potentially involved in host range (ankyrin repeat proteins; ORFs 003, 004, 008, 118, 123, 126, 128, and 129), immune evasion (ORF 127), and immune modulation (ORF 117) and genes affecting virulence (ORF 006/132). Notably, PPVs contain a dUTPase gene previously characterized in ORFV (22, 43) but lack homologues of other ChPV genes likely involved in nucleotide metabolism, making this class of genes underrepresented in PPVs.
|
View this table: [in a new window] |
TABLE 1. ORFV and BPSV ORFs
|
BPSV and ORFV contain 15 and 16 ORFs, respectively, which share no significant homology to known proteins and are primarily located at the right end of the genome. Fourteen of these ORFs (ORFs 001, 005, 012, 013, 024, 073, 113, 115, 116, 119 to 121, 124, and 125) are present in both BPSV and ORFV, with amino acid identities ranging from 29 to 64%, two (ORFs 002 and 118) are present only in ORFV, and one (ORF 133) is unique to BPSV. Consistent with a possible host range function, homologues of six of these (ORFs 012, 024, 118, 119, 121, and 125) are transcribed at early times in cells infected with ORFV strain orf-11 (Table 1).
Of the 26 most distantly related ORFs between BPSV and ORFV (amino acid identity < 60%), 10 are unique to PPVs (ORFs 005, 012, 013, 113, 115, 116, 119 to 121, and 124), 3 are characterized ORFV NZ2 strain genes with putative (ORFs 020 and 117) or known (ORF 006/132) roles in pathogenesis, 10 show homology with VACV genes of known structure or function (ORFs 061, 080, 088, 103, 109, 110, 112, 126, 128, and 129), and 3 show homology with VACV virus ORFs of unknown function (ORFs 009, 016, and 122).
Highly variable PPV proteins include homologues of the VACV proteins H5R, A4L, A12L, A26L, A33R, A34R, M1L, and B4R. BPSV and ORFV 061, orthologues of VACV H5R gene, are only 51% identical. H5R encodes a late transcription factor (VLTF-4) which is synthesized before and after DNA synthesis, is phosphorylated by viral kinases, and is hypothesized to have multiple roles in the viral replicative cycle (5, 36). Notably, in closely related capripoxviruses, sheeppox virus and goatpox virus (genomes which share 96% nucleotide identity), VLTF-4 homologues are among the least conserved genes. It is tempting to speculate that PPV ORF 061 may play a role in host range during virus infection.
BPSV and ORFV 080 encode homologues of VACV A4L, a gene with significant variability in other poxvirus genera. A4L encodes an immunodominant late protein associated with the virion core and necessary for viral morphogenesis (84). OV-SA00 and OV-IA82 080 encode products that are 84 and 80 amino acids longer than the BPSV 080 product, respectively, due in part to the lack of four Cys-(Pro-Ala)3 motifs separated by additional Pro/Ala-rich sequences in BPSV 080. Similar Pro/Ala-rich repeats are present in the molluscum contagiosum virus (MOCV) orthologue MC107L but not in A4L. Tandem repeat motifs in A4L-like proteins are thought to be involved in protein-protein interactions and antigenicity (7).
BPSV and ORFV 088, orthologues of VACV A12L virion core protein, share only 56% amino acid identity, an unusual degree of intragenus variability for A12L orthologues (e.g., >90% amino acid identity within orthopoxvirus [OPV], leporipoxvirus, and capripoxvirus). Notably, BPSV and ORFV 088 encode proteins of 223 and 260 amino acid, respectively, while non-PPV ChPV A12L orthologues are 156 to 195 amino acids. The difference in length is partially due to a positively charged 20-residue insertion immediately downstream of the predicted VNA/GG cleavage site (position 170 in ORFV 088) and might suggest specific PPV core structure requirements (83).
PPV 102 and 103 are variable ORFs, with PPV 102 being more conserved between species than 103 (67 to 76% versus 57 to 58% amino acid identity, respectively). PPV 102 and 103 share similarity to each other (32 to 37% amino acid identity) and to homologues of OPV genes involved in formation of virus-filled A-type inclusions (ATIs) (46, 80). Both 102 and 103 are most similar to homologues of OPV P4c, an intracellular mature virion (IMV)-specific protein which helps direct IMV into ATIs. The carboxy termini of PPV 102 and OPV P4c proteins share sequences found in the VACV A27L fusion protein. PPV 103 has weak similarity to OPV ATI proteins, which constitute the crystalline matrix of OPV ATIs. Given the variable nature of these genes in different ChPV genera, PPV homologues may affect genera-specific and species-specific mechanisms of retaining or disseminating IMVs.
PPV 109 and 110 are orthologues of VACV A33R and A34R, respectively, genes encoding envelope type II glycoproteins expressed in intracellular enveloped virions and in EEV (44, 66). Mutations in these genes affect EEV formation (A33R and A34R), cell-to-cell spread of virus (A33R), and infectivity and virulence (A34R) (reviewed in reference 74). PPV 109 and 110, although distantly related to the VACV orthologues, are predicted to contain amino-terminal transmembrane regions and external Cys residues suggesting a similar protein topology and structure. Notably, OV-SA00 109 and 110 are as distantly related to OV-IA82 orthologues (56 and 49% amino acid identity) as they are to BPSV proteins (49 to 50% amino acid identity) (Table 1), with amino acids differences being largely concentrated in the predicted external domain. An explanation for this intraspecies variability is not immediately obvious. Alignment of available ORFV 109 sequences grouped OV-IA82, NZ2, and Orf-11 A33R homologues in a single cluster with 80 to 97% amino acid identity, excluding OV-SA00 A33R, which showed 51 to 55% amino acid identity relative to other ORFV sequences. A similar situation is observed when available ORF 110 sequences are compared. OV-SA00 is the only strain originally isolated from goats, whereas OV-IA82, NZ2, and ORF-11 strains were isolated from sheep. Thus, there appears to be a correlation between the species from which the virus was isolated and clustering of ORFV 109 and 110 amino acid sequences. This raises the possibility that differences in the external domain of both A33R and A34R are associated with host-specific requirements during virus infection by EEV. Differences in disease course have been shown following experimental infection with ORFV isolated from sheep or goats (87).
Variable PPV 126, 128, and 129 correspond to the Ank-1, Ank-2, and Ank-3 genes previously described for the BPSV B177 and ORFV D1701 strains (67). These ORFs are 58, 57, and 50% identical between BPSV and ORFV, respectively, and encode ankyrin repeat-containing proteins (ARPs). Cellular ARPs engage other proteins to form regulatory complexes which are involved in the control of processes such as cell cycle and cell differentiation (71). Many poxviruses encode ARPs, several of which have been linked to host range functions, apoptosis inhibition, and virulence (34, 56). Notably, BPSV contains two additional ARPs in the left terminal genomic region (ORFs 003 and 004) which are not present in ORFV (Table 1).
PPV genes involved in pathogenesis. ORFV encode several proteins with known or putative roles in pathogenesis: 020, an orthologue of vaccinia E3L that functions in interferon (IFN) resistance; 127, a viral homologue of IL-10; 117, a secreted granulocyte-macrophage colony-stimulating factor (GM-CSF) inhibitor (GIF); 112, a putative chemokine-binding protein; and 132, a viral homologue of VEGF. These genes are present in BPSV with predicted amino acid identities ranging from 37 to 77%.
E3L encodes a double-stranded-RNA (dsRNA)-binding protein kinase inhibitor with orthologues in all ChPV genera except Avipoxvirus and Molluscipoxvirus. The E3L gene product provides IFN resistance to VACV-infected cells and broad host range to virus infection in tissue culture and is a virulence determinant (4, 9). The ORFV NZ2 strain E3L homologue is also involved in IFN resistance (31) and is 93% and 53% identical to its OV-SA00 and BPSV counterparts, respectively. ORFV and BPSV 020 are most similar at the carboxy-terminal half of the protein, which is predicted to bind dsRNA through a conserved binding motif (10). The amino-terminal half of the protein is less conserved (45% amino acid identity) and includes two deletions of six and four amino acids in ORFV. In VACV, the amino-terminal half of E3L is dispensable for replication in cell culture and is not required for IFN resistance. However, this region is required for full virulence in a mouse intranasal inoculation model (9). It is thus possible that differences between ORFV and BPSV in the amino-terminal half of 020 are significant for host range and pathogenesis in their respective hosts.
OV-SA00 and OV-IA82 127 are orthologues (95% amino acid identity) of the previously described ORFV NZ2 and NZ7 strain IL-10 genes (23). BV-AR02 127 is divergent from ORFV homologues (77% amino acid identity) with most amino acid differences concentrated in the amino terminal third of the protein (27% amino acid identity in the first 50 amino acids). Nevertheless, a putative signal peptide is present in the amino terminus of all three IL-10 homologues presented here. The carboxy two-thirds of PPV IL-10 are highly conserved with cellular IL-10, with all ORFV interleukin 10 (IL-10) proteins sequenced to date being most similar to ovine IL-10 (23; this work). Notably, and in agreement with previous results, BV-AR02 127 shares in this conserved region six residues identical to bovine but not to ovine and ORFV IL-10, including a His residue at position 75 predicted to interact with the IL-10 receptor chain 1 (67). These features may reflect specific adaptation to the natural host.
PPVs 117 are orthologues of ORFV GIF, a protein which binds and inhibits GM-CSF and IL-2 in vitro and may function as an immunomodulatory factor in vivo (15). OV-SA00, OV-IA82, and NZ2 strain GIF homologues are very similar to each other (94% amino acid identity), containing at least six potentially structurally significant Cys residues and an amino-terminal signal peptide. While BV-AR02 GIF shares these structural features, it is only 38 to 40% identical to ORFV GIF, an unexpected divergence considering the similarity between ovine and bovine GM-CSF and IL-2 (83 and 96% amino acid identity, respectively). Notably, PPV 117 shares 22 to 25% amino acid identity with PPV 112, a gene also predicted to encode a signal peptide and expressed at early times postinfection in ORFV strain Orf-11 (Table 1). Both PPV 117 and 112 share limited sequence similarity and/or Cys patterns with VACV C23L and myxoma virus MT-1 chemokine binding proteins and VACV A41L virulence factor. Taken together, these data suggest that PPV 117 and 112 may be members of a divergent family of poxviral chemokine- and cytokine-binding virulence factors.
ORFV 132 and BPSV 006 encode homologues of mammalian VEGFs, angiogenic factors that bind receptor tyrosine kinases to affect embryonic development and tumor neovascularization (11, 53, 62, 85). ORFV, PCPV, and BPSV encode the only known viral VEGFs (vVEGFs), all of which contain a characteristic cystine knot motif, a potential signal sequence, potential N- and O-linked glycosylation sites, a carboxy-terminal Thr/Pro-rich motif unique to vVEGFs, and an Asp residue (position 85 in BPSV VEGF) associated with specific VEGF receptor binding (38, 79; this paper). All vVEGFs are flanked by similar putative transcriptional control elements (38, 79; this paper) suggesting that, as is the case for ORFV (38), these genes are expressed at early times postinfection. ORFV VEGF is known to play a role in ORFV pathogenesis associated with vascularization and epidermal lesion proliferation (70).
OV-SA00 and OV-IA82 VEGF are 38% identical to each other and most similar to NZ7 and NZ2-like VEGFs (90 and 80% amino acid identity, respectively). Previous sequence analysis of ORFV isolates from diverse geographic origins segregated VEGFs into two divergent groups, a more conserved NZ7-like group and a more variable NZ-2 group (52). The data presented here for U.S. ORFV isolates further supports the notion that VEGF type is independent of the geographic origin (52).
BV-AR02 006 represents a novel variant of PPV VEGF, previously not found in BPSV strain B177 by DNA hybridization (67). BV-AR02 006 is located in a BPSV-specific, left terminal genomic region contrasting the right terminal location of ORFV and PCPV VEGFs. BV-AR02 VEGF is 35 to 50% identical to other vVEGFs and contains a unique charged pentapeptide located at positions 34 to 39. This suggests that, as for PCPV, BPSV VEGF is distinct from ORFV VEGF prototypes. Sequence divergence revealed here may explain the lack of hybridization when BPSV B177 strain DNA was probed with ORFV VEGF probe (67). Alternatively, terminal genomic variability observed for BPSV isolates (26) may have resulted in loss of the gene from the B177 strain. The presence of vVEGF in BPSV suggests its importance in PPV pathogenesis; however, the divergent nature of BPSV VEGF may imply functions or binding specificities distinct from other vVEGFs.
Comparison of PPV with other poxvirus genera.
BPSV and ORFV contain 16 and 17 ORFs, respectively, which have no significant homology to genes from other poxvirus genera, and with the exception of VEGF, to other known proteins (Table 1). Although six of them are transcribed in ORFV strain Orf-11-infected cell cultures, their functions are not known. BPSV and ORFV contain a total of 113 and 111 genes, respectively, with homology to genes from other poxvirus genera. These include homologues of 88 of the 90 genes conserved within all other ChPVs, with 7 of the 11 most similar (
60% amino acid identity; Table 1) involved in transcription, indicating that PPVs utilize basic ChPV replicative mechanisms (81). PPVs are unique within the ChPV subfamily in that they lack homologues of VACV F15L, a gene of unknown function, and VACV D9R, a gene encoding a putative nucleoside triphosphate pyrophosphohydrolase containing a mutT motif similar to that in VACV D10R, a protein affecting viral transcription (55).
PPVs, although distinct, share a number of features with MOCV, the sole member of the molluscipoxvirus genus. PPV and MOCV are the only characterized poxviruses with genomes rich in G+C (64%), while others are rich in A+T. Homology searches revealed that 46 of 104 PPV proteins were most similar to MOCV orthologues, while 26 proteins were more similar to OPV orthologues (Table 1 and data not shown).
PPV 014, 015, and 029 are putative orthologues of MOCV MC026L, 027L, and 043L, respectively, based on amino acid identity and similar genomic location (72, 73). These ORFs of unknown function have no counterparts outside PPV and MOCV and are 61 to 68% identical between ORFV and BPSV. PPV 029 is transcribed at early times postinfection in Orf-11-infected cells (Table 1). PPV 014 and MC026L contain a RING-H2 motif which is present in proteins from diverse organisms. RING-H2 proteins are subunits of heteromeric ubiquitin ligases (E3s) which affect multiple cellular processes including cell cycle regulation and immune response (37).
PPVs and MOCV both lack genes present or conserved in other poxviruses. These comprise homologues of most poxviral genes likely involved in nucleotide metabolism, including homologues of OPV ribonucleotide reductase, thymidine kinase, guanylate kinase, thymidylate kinase, and a putative ribonucleotide reductase cofactor, VACV glutaredoxin O2L (3). PPVs, MOCV, and Melanoplus sanguinipes entomopoxvirus are the only poxviruses known to lack a thymidine kinase gene. In contrast, PPV 007, a gene not essential for growth of ORFV in vitro, is a dUTPase gene missing in MOCV (22, 43, 72, 73). Notably, PPVs and MOCV are the only ChPVs lacking homologues of VACV B1R, a Ser/Thr protein kinase similar to cellular VRK1 homologues and giving a temperature-sensitive DNA-negative phenotype (61). PPVs and MOCV lack homologues of VACV A50R DNA ligase, a gene also absent in yatapoxviruses (81). Also absent in PPV and MOCV are serine protease inhibitor and kelch-like gene families present in other ChPVs. These gene families are known to affect host responses, including inflammation, apoptosis, complement activation, and coagulation (39), and are associated with virulence (78). The lack of ChPV-like genes in both PPVs and MOCV may reflect adaptation for specific tissue tropism, which is notable considering that PPVs and MOCV appear to replicate in cycling cells (17, 45).
Features similar in both PPV and MOCVnucleotide composition, amino acid similarity, and gene contentsuggest that they are distinct from other known mammalian poxvirus genera. Phylogenetic analysis of protein sequences also supports the idea that, although divergent, PPVs and MOCV are distinct from other known mammalian ChPVs (Fig. 1).
![]() View larger version (14K): [in a new window] |
FIG. 1. Phylogenetic analysis of PPV proteins. PPV025 DNA polymerase and homologous sequences were aligned using ClustalW. Unrooted trees were generated using the neighbor-joining algorithm with Poisson correction for multiple substitutions. Bootstrap values greater than 95% after 1,000 replicates are indicated at appropriate nodes. Homologous protein sequences from the following viruses and accession numbers were compared: ORFV, AY386264; MOCV, MCU60315; VACV, M35027; yaba-like disease virus (YLDV), YDI293568; lumpy skin disease virus (LSDV), AF325528; myxoma virus (MYXV), AF170726; and swinepox virus (SWPV), AF410153. Similar results were obtained for 16 additional conserved PPV proteins, with 15 maintaining 100% bootstrap support for separation of PPV and MOCV from other mammalian ChPVs. Similar results were also obtained for these 17 proteins using the maximum likelihood algorithm with Dayhoff correction for multiple substitutions and for whole genomic nucleotide alignment utilizing amino acid translation as implemented by using Dialign (54).
|
|
|
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»