J Gen Virol Tips for Better Browsing
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Gen Virol 87 (2006), 2191-2201; DOI 10.1099/vir.0.81738-0

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Jarmin, S.
Right arrow Articles by Skinner, M. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Jarmin, S.
Right arrow Articles by Skinner, M. A.
Agricola
Right arrow Articles by Jarmin, S.
Right arrow Articles by Skinner, M. A.
© 2006 Society for General Microbiology

Avipoxvirus phylogenetics: identification of a PCR length polymorphism that discriminates between the two major clades

Susan Jarmin1,{dagger}, Ruth Manvell2, Richard E. Gough2, Stephen M. Laidlaw1,{dagger} and Michael A. Skinner1,{dagger}

1 Institute for Animal Health, Compton, Newbury, Berkshire RG20 7NN, UK
2 Veterinary Laboratory Agency, New Haw, Addlestone, Surrey KT15 3NB, UK

Correspondence
Michael A. Skinner
m.skinner{at}imperial.ac.uk


   ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Avipoxvirus infections have been observed in an extensive range of wild, captive and domesticated avian hosts, yet little is known about the genome diversity and host-range specificity of the causative agent(s). Genome-sequence data are largely restricted to Fowlpox virus (FWPV) and Canarypox virus (CNPV), which have been sequenced completely, showing considerable divergence between them. It is therefore proving difficult, by empirical approaches, to identify pan-genus, avipoxvirus-specific oligonucleotide probes for PCR and sequencing to support phylogenetic studies. A previous preliminary study used the fpv167 locus, which encodes orthologues of vaccinia virus core protein P4b (A3). PCR per se did not discriminate between viruses, but restriction-enzyme or sequence analysis indicated that the avipoxviruses clustered either with FWPV or with CNPV. Here, further study of the P4b locus demonstrated a third cluster, from psittacine birds. A newly identified locus, flanking fpv140 (orthologue of vaccinia virus H3L), confirms the taxonomic structure. This locus is particularly useful in that viruses from the fowlpox-like and canarypox-like clusters can be discriminated by PCR on the basis of fragment size, whilst sequence comparison allows discrimination for the first time between Pigeonpox virus and Turkeypox virus. Except within the psittacines, virus and avian host taxonomies do not show tight correlation, with viruses from the same species located in very different clades. Nor are all the existing recognized avipoxvirus species, defined primarily by avian host species (such as CNPV and Sparrowpox virus), resolved within the present structure.

The GenBank/EMBL/DDBJ accession numbers for the sequences reported in this paper are shown in Table 1Up.

{dagger}Present address: Department of Virology, Imperial College London, Faculty of Medicine, St Mary's Campus, Norfolk Place, London W2 1PG, UK. Back


   INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Avipoxvirus infections have been observed in more than 230 of the known 9000 species of bird, spanning 23 orders (Bolte et al., 1999Down), yet little is known about the genome diversity and host-range specificity of the causative agent(s). Avian poxviruses are all classified in the genus Avipoxvirus of the family Poxviridae within the subfamily Chordopoxvirinae. In common with other poxviruses, they contain a double-stranded DNA genome, ranging from 230 to >300 kbp, and replicate within the cytoplasm of the cells that they infect. Fowlpox virus (FWPV) is the prototypic member of the avipoxviruses. It is well known that the diseases present in two main forms, which may be the result of different routes of infection. The cutaneous form occurs following mechanical trauma, including bites of arthropods that serve as vectors for mechanical transmission of the virus. Lesions normally occur on unfeathered areas around the eyes and nares, on the comb, wattle or legs. In poultry, this form is associated with low mortality and, frequently, subclinical infection. The diphtheritic form, following inhalation or ingestion of virus, involves mucous membranes of the mouth, pharynx, larynx and sometimes trachea. It is associated with a higher mortality rate, usually due to occlusion of the oropharynx. Avipoxvirus infections of canaries (Johnson & Castro, 1986Down) and finches can cause significantly higher rates of mortality than those of chickens or turkeys. A third, pneumonia-like form of the disease has been described in canaries (Johnson & Castro, 1986Down).

Conventional diagnosis of avipoxviruses is carried out by histopathological examination, electron microscopy, virus isolation on chorioallantoic membranes of embryonated chicken eggs, serological methods and restriction-fragment length polymorphism. These detection methods are being supplemented by the use of PCR using primers for a 578 bp product of fpv167, the orthologue of vaccinia virus A3L that encodes P4b (Huw Lee & Hwa Lee, 1997Down). [To distinguish genes or open reading frames (ORFs) from their encoded protein products, italicized lower case is used for genes and ORFs (e.g. fpv030) and plain text (with capitalized first letter) for the products that they encode (e.g. Fpv030)].

Vaccines are available that utilize FWPV, Canarypox virus (CNPV), Pigeonpox virus (PGPV) or Quailpox virus (QUPV), but, due to new trends in avian farming and conservation, there is a demand for vaccines for use in other bird species, such as rare, and often expensive, falcons, eagles and ostriches. Vaccine application is by trial and error, with some viruses inducing cross-protection. For instance, vaccination with PGPV-based vaccines, but not those based on QUPV or Turkeypox virus (TKPV) (Winterfield & Reed, 1985Down), protects chickens from fowlpox. Such antigenic differences have served as the basis for avipoxvirus classification (Tripathy, 1984Down).

Little is known about the host ranges of avipoxviruses, although they are generally assumed to be limited. For instance, a virus isolated recently from an Andean condor (Vultur gryphus) caused an aggressive, diphtheritic form of the disease in the condor, but produced only small, localized, cutaneous lesions in inoculated chickens, inducing no cross-protection against challenge with FWPV (Kim et al., 2003Down).

Partial sequencing of the CNPV genome around the thymidine kinase locus (Amano et al., 1999Down) first showed that it was surprisingly highly diverged from that of FWPV. This was confirmed when the complete genome sequence of CNPV was derived (Tulman et al., 2004Down), allowing it to be compared in its entirety with that of FWPV (Afonso et al., 2000Down). The extent of divergence is illustrated by the fact that CNPV and FWPV orthologues of P4b, a conserved 75.2 kDa virion core protein found in all poxviruses (Binns et al., 1989Down; Huw Lee & Hwa Lee, 1997Down), share only 64.2 % amino acid identity.

Despite sequencing of the complete genomes of FWPV and CNPV, little is known about the phylogenetic relationship between these and other avipoxviruses and how this might relate to host range, pathogenesis and epidemiology. Concern has been raised that poxvirus infections are responsible for the decline in populations of native birds in Hawaii and the Canary Islands (Tripathy et al., 2000Down) and that they might threaten marginal species, such as the flightless, nocturnal New Zealand owl parrot or kakapo (Stone & Forbes, 2002–2003Down). Recent preliminary phylogenetic studies of avipoxvirus isolates using the fpv167 locus (Luschow et al., 2004Down; Weli et al., 2004Down) indicated that most of the isolates clustered around either CNPV or FWPV.

Corroboration of the phylogenetics of the avipoxviruses will require data from additional loci. For diagnostic purposes, it would also be useful to have fragment-length polymorphisms that could be detected directly by PCR rather than having to resort to the restriction-enzyme digestion that was used at the fpv167 locus by Luschow et al. (2004)Down.

In an attempt to address these issues, a collection of 22 avipoxviruses, comprising isolates, vaccines and laboratory strains, was examined, initially at the fpv167 locus, then subjected to PCR at the loci of genes encoding the three immunodominant FWPV antigens Fpv140, Fpv168 and Fpv191 (Boulanger et al., 2002Down), using a variety of primers with sequences specific for targets conserved between FWPV and CNPV. Sequence data obtained successfully from the fpv140 locus, as well as from the fpv167 locus, were subjected to phylogenetic analysis. The implications of the results for avipoxvirus epidemiology, virus host range and nomenclature are considered.


   METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Viruses and cells.
Twenty-two avipoxviruses were used in this study (Table 1Down). Five were commercial fowlpox vaccines, Diftosec CT (Rhone Merieux), Nobilis Variole W (Intervet), Chick-N-Pox, Websters FPV M (Salsbury Laboratories Inc.) and Poxine (Duphar). One was a canarypox vaccine, Canarypox V (Duphar). One was the attenuated FWPV FP9, the genome of which has been sequenced completely (Laidlaw & Skinner, 2004Down). The remaining 16 viruses were isolates obtained from outbreaks and samples collected at the Central Veterinary Laboratory (CVL) of the Veterinary Laboratory Agency (VLA), Weybridge, UK. These poxviruses from chicken (Gallus gallus), canary (Serinus canaria), pigeon (Columba livia), turkey (Meleagris gallopavo), sparrow (Passer domesticus), starling (Sturnus vulgaris), houbara bustard (Chlamydotis undulata), macaw (Ara spp.), parrot (Amazona spp.), falcon (Falco spp.) and albatross (Diomedea melanophis) were isolated from birds with clinical signs of infection. Most of the isolates had been passaged on chick embryo fibroblasts (CEFs) and/or specific-pathogen-free (SPF) embryonated chicken eggs between five and 30 times.


View this table:
[in this window]
[in a new window]
 
Table 1. Details of avipoxviruses used in this study

CVL, Central Veterinary Laboratories, Weybridge, UK.

 
Samples were propagated for PCR on CEFs under E199 medium (1x E199 with Earle's salts, 10 % tryptose phosphate broth, 3 % sodium bicarbonate, 1 % L-glutamine, 0.1 % penicillin/streptomycin, 2.5 % nystatin, 2 % newborn calf serum). The samples were incubated for 4 days under 5 % CO2 at 37 °C or until cytopathic effects were observed.

Virus purification and DNA extraction.
Virus DNA was isolated by using standard methods, as described previously (Boulanger et al., 1998Down).

PCR amplification.
Primers described previously (Huw Lee & Hwa Lee, 1997Down) were used to amplify the FWPV P4b gene (fpv167). The sequences of the primers were M2925, 5'-CAGCAGGTGCTAAACAACAA-3' (fpv167 nt 459–478), and M2926, 5'-CGGTAGCTTAACGCCGAATA-3' (complementary to nt 1016–1035).

The 3' end of fpv139 through to the 5' end of fpv141 was amplified by PCR using primers designed on sequence conserved between CNPV and FWPV. M2904, 5'-GAAGTAGAGTTACGGTTC-3', is identical to nt 171290–171311 of FP9 (Laidlaw & Skinner, 2004Down) in fpv139, and M2912, 5'-GGTGATCCATTTCCATTTC-3'), is complementary to nt 173254–173272 in fpv141. PCR amplification was performed essentially as described by Huw Lee & Hwa Lee (1997)Down, but using Expand High Fidelity PCR Enzyme mix (2.6 U; Roche).

Electrophoresis, purification and sequencing of PCR product.
PCR samples (5 µl) were loaded onto 1 % agarose gels [90 mA for 30 min in TBE buffer (0.1 M Tris, 0.09 M boric acid, 0.01 M EDTA)]. PCR products were purified by using a QIAquick PCR Purification kit (Qiagen) and resuspended in 50 µl H2O. These products were then sequenced on both strands, using the amplification primers and internal primers in an automated sequencer (Beckman CEQ 8000).

For FWPV-like viruses, primers M3238 (5'-GTTAGAAGAGGAGATGGAGATGG-3'), M3239 (5'-CATCATAGTTATTTAACTAAATC-3'), M3242 (5'-GTCTACAGACGTCAGAATAATAACC-3') and M3243 (5'-GGTTATTATTCTGACGTCTGTAGAC-3') were used.

For CNPV-like viruses, primers M3240 (5'-ATAAGAAGAGGAGAAGGAGACGG-3'), M3241 (5'-GTATACGCTATTAAATGAACATTC-3'), M3244 (5'-GCGATAGAATTTTTTGTTTC-3'), M3245 (5'-GGATGACAGATTCTATGACG-3'), M3251 (5'-GGAGAGCTTAGTTTCTCTATC-3'), M3252 (5'-GCATTGATTAAAATCTATATCACC-3'), M3253 (5'-GCACTAATATCATACTTTCCAAAC-3'), M3254 (5'-GGTATAAATAGATATGTTATGAGTAACG-3'), M3255 (5'-CGCGTAAAAAAAGATACGGG-3'), M3256 (5'-GGTTTACAGTATCGGGTATATGG-3'), M3257 (5'-CGATACTAGTCATCTATTATTAG-3') and M3258 (5'-GCCGAAATCAACACATTTAAATC-3') were used to extend the reads. Sequencing PCR was carried out in a total volume of 10 µl, containing 2 µl purified PCR product, 2 µl Quick Start Enzyme mix (Beckman), 0.5 µl buffer and 5 pmol primer. Samples were subjected to 30 cycles (20 s at 96 °C, 20 s at 40 °C, 4 min at 60 °C) in a TC-412 (Techne).

The sequencing samples were run by using the LFR program (Beckman). Sequence analysis was carried out with STADEN (Staden, 1982Down) and alignments were made with CLUSTAL_X (Thompson et al., 1997Down), manipulated with Se-Al (Rambaut, 1996Down; http://evolve.zoo.ox.ac.uk/) and displayed with MacGDE (http://www.msu.edu/~lintone/macgde/). Phylogenetic analysis was performed by using PAUP, with approaches and methods as described by Hall (2001)Down.


   RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
P4b gene (fpv167) locus
Primers for fpv167 were used to amplify sequences from DNA extracted from isolates propagated on CEFs. These primers gave products of about 600 bp (predicted size, 578 bp) on all viruses tested (data not shown). The products were purified and sequenced; a common sequence of 450 nt was obtained from all of the isolates, allowing comparison with equivalent, previously published sequences in the databases (Table 2Down). The derived and aligned protein sequences (Fig. 1aDown) were used to generate phylogenetic trees for P4b, including viruses from other chordopoxvirus genera, by using neighbour-joining analysis (with bootstrapping) and parsimony. Trees were also derived for the P4b-encoding sequences of the avipoxviruses, with Molluscum contagiosum virus (MOCV) as an outgroup, using neighbour-joining analysis (with bootstrapping), maximum likelihood and parsimony.


View this table:
[in this window]
[in a new window]
 
Table 2. Details of poxvirus sequences obtained from GenBank

 

Figure 1
View larger version (10K):
[in this window]
[in a new window]
 
Fig. 1. (a) Fpv167 (P4b) amino acid sequences aligned by using CLUSTAL_X and displayed by using MacGDE and Se-Al. See Table 2Up for GenBank accession numbers. FWPV-like: FWPVD, FWPVN, FWPVM, FWPV174, FWPVHPB, TKPV13401, SRPVDD1258, FWPVFP9, FWPVV250, PGPVTP2, OSPV, PGPVP, TKPV66, TKPV98, ABPV353 and FLPV138196; CNPVATCC-like: CNPVATCC, SCPV182, GTPVA311, GTV256, SPRV32002, HOPV1252 and SRPV9037. In the consensus line, strictly conserved residues are in upper case. (b) Fpv139 amino acid alignment. FWPV-like: FWPVFP9, FWPVD, FWPVN, FWPVM, FWPV174 and FWPVC; CNPVATCC-like: CNPVATCC, CNPV1445, HOPV1165, HOPV1252 and SRPV9037.

 
The grouping of the mammalian poxviruses in the P4b protein trees differs slightly from that in trees based on consensus sequences of many conserved proteins (McLysaght et al., 2003Down; Gubser et al., 2004Down). In the consensus trees, yatapoxviruses form the most proximal branch of the group containing leporipox-, suipox- and capripoxviruses after that group's divergence from the orthopoxviruses. In the P4b protein trees, the yatapoxviruses are in the same group, but are in a distal node with Swinepox virus. It must be acknowledged that the single-protein trees are more susceptible to functional pressures upon the particular protein than are consensus trees. However, our main interest in the protein trees was to compare the relative divergence of avian and mammalian viruses (acknowledging that the rates of divergence may vary among viruses in these different classes of vertebrates, thereby preventing any conclusions about the timescales involved). Such comparison, by neighbour-joining and parsimony analysis, indicates that the extent of protein divergence between avipoxvirus clades is equivalent to that between the various genera of mammalian poxviruses (Fig. 2Down).


Figure 2
View larger version (13K):
[in this window]
[in a new window]
 
Fig. 2. Phylogenetic analysis of avipoxvirus P4b orthologues. Abbreviations are as in Tables 1 and 2UpUp. (a) Unrooted, neighbour-joining phylogram of partial protein sequences for all genera in the subfamily Chordopoxvirinae. Bootstrap values (1000 replicates) of >60 % are indicated. (b) Neighbour-joining phylogram of DNA sequences from genes encoding avipoxvirus and MOCV P4b orthologues, rooted on MOCV. Bootstrap values (1000 replicates) of >80 % are shown. Sequences derived as part of this study are indicated by a dot. Avipoxvirus clades A–C and subclades are labelled. Isolates from the same genus of host located in different clades or subclades are as indicated: *, pigeons; +, falcons. SRPVDD1258 and TKPV13401 are not thus marked as it could not be excluded that they were really FWPV (Luschow et al., 2004Down).

 
Protein and DNA trees show the avipoxviruses forming two major clades (A and B), with strong support by neighbour-joining bootstrap, parsimony and (for DNA) maximum likelihood (Fig. 2Up). In protein and DNA trees, clade B diverges into two subclades, B1 and B2, also with strong support. B1 consists of the isolates from canaries, great tits, houbara, stone curlew and four of the five sparrow isolates, whilst B2 comprises the starling isolate as well as two of the four viruses from pigeons. In protein and DNA trees, clade A comprises FWPV, TKPV, Falconpox virus (with one of the two, FLPV36202, as a proximal outlier), Ospreypox virus and Albatrosspox virus, as well as the remaining isolates from pigeons and sparrow. In the DNA trees, the clade A viruses are further resolved (with strong support by neighbour-joining bootstrap, parsimony and maximum likelihood). FLPV36202 still forms a proximal outlier (subclade A4) to clade A. Subclade A1, with strong support, comprises essentially the FWPVs plus two single isolates from a turkey and a sparrow. Subclade A2, with strong support, consists of the two (of four) isolates from pigeons, two viruses from turkeys (isolated 30 years apart in the UK and Italy) and an isolate from an osprey. Subclade A3, also with strong support, consists of an albatross virus and the second falcon virus.

In the DNA trees, a third major clade, C, which consists of the psittacine viruses, branches proximal to the MOCV outgroup, with strong support by neighbour-joining bootstrap, parsimony and maximum likelihood. In the protein trees, however, placement of the psittacine group is more ambiguous. By neighbour joining, the psittacine group clusters (at a bootstrap value of 71 %) with FWPV and the FWPV-like viruses of clade A. By parsimony, it clusters with CNPV and the CNPV-like viruses of clade B in 75 % of trees but, as in the DNA trees, it forms a distinct clade, proximal to MOCV, in 25 % of trees. Detailed examination of the various coding sequences shows that a number of nucleic acid changes unique to the psittacine viruses are masked in the protein sequence, due to codon redundancy. We therefore conclude that the DNA trees represent a more accurate description of the status of the psittacine viruses as a distinct clade, in accord with in vivo cross-protection studies (Boosinger et al., 1982Down; Winterfield & Reed, 1985Down).

In the DNA trees, clade A represents the clusters of FWPV and FWPV-like viruses seen previously [clusters I–III of Luschow et al. (2004)Down]. Clade B represents CNPV and the CNPV-like viruses seen previously [cluster IV of Luschow et al. (2004)Down]. Branch V observed by Luschow et al. (2004)Down consisted of only one virus, an isolate from agapornis. We can now see that it is joined by viruses from two other psittacines, parrot and macaw, to form clade C. This clade is the least divergent of the three, although this may be due to it consisting of only three isolates. Branch IV observed by Luschow et al. (2004)Down is represented by subclade B1.

Clade A.
Sequence analysis at the P4b locus revealed that the FWPVs (FWPVFP9, FWPV174, FWPVD, FWPVHPB and FWPVN), in subclade A1, were identical to each other as well as to SRPVDD1258 and TKPV13401, both German clinical isolates (Fig. 1aUp). The subclade A1 viruses showed 95.3 % nucleic acid identity (95 % amino acid identity) to FLPV36202 in subclade A4. Two other TKPV isolates (66 and 98) could be distinguished from FWPV (91 % nucleic acid identity), but not from PGPVP and PGPVTP2, as they clustered on the same branch of subclade A2, with 100 % nucleic acid identity, as they did with OSPV 724/01-20. ABPV and FLPV138196 (subclade A3) show 99 % nucleic acid identity to the subclade A2 viruses.

Clade B.
Viruses in clade B, all with 98 % or more nucleic acid identity, show about 80 % nucleic and amino acid identity to viruses in clades A and C.

Clade C.
MCPV and PRPV have identical sequences and show 99 % nucleic acid identity to AGPV; thus, these viruses form clade C, with 86 % amino acid identity to viruses in clade A1 and 80 % to those in clade B1. Although the low divergence within this clade may be attributable to the small number of isolates, the isolates were apparently well separated in terms of host species, location and time: MCPV and PRPV were isolated from the different species in UK quarantine facilities some 3 years apart and AGPV was isolated in Germany (Luschow et al., 2004Down).

H3L gene (fpv140) locus
PCR products for the H3L (fpv140) gene locus, amplified from primers in fpv139 and fpv141, were obtained for viruses from clades A (FWPV-like) and B (CNPV-like), but not C (psittacine viruses). A clade-specific fragment-length polymorphism was observed for the PCR product, with each virus producing a product of either 1800 bp (for clade A viruses) or 2400 bp (for clade B viruses). Comparison of the known genome sequences for FWPV and CNPV indicated that the polymorphism was probably due to the absence (from clade A viruses) or presence (in clade B viruses) of a gene, CNPV185 or its equivalents, inserted between the orthologues of fpv139 (CNPV184) and fpv140 (CNPV186; Fig. 3Down).


Figure 3
View larger version (31K):
[in this window]
[in a new window]
 
Fig. 3. (a) Schematic (to scale) showing the arrangement of FWPV genes fpv138–fpv141 aligned with their CNPV orthologues. CNPV185 has no FWPV orthologue.Primers M2904 and M2912 are shown below the CNPV sequence to which they bind. (b) PCR amplification of the fpv140 gene from avipoxvirus isolates. M, {lambda}X174 DNA digested with HaeIII. See Table 1Up for virus abbreviations.

 
Because of the presence of the insertion, sequences were determined from each end of the products by using internal primers in combination with the flanking primers M2904 and M2912. Sequence analysis of the PCR products confirmed the anticipated correlation between sequence-based phylogeny and the size of PCR product. All clade A viruses, like FWPV, encode a product for fpv139 consisting of 190 aa, whereas clade B viruses encode products, like CNPV184, of 189 aa. Within clade A, genes fpv141 and fpv140 are translated from the same reading frame, the two ORFs being separated only by the stop codon of fpv141. In clade B, in contrast, the orthologous proteins are translated from different reading frames, the start codon of CNPV186 being separated from the stop codon of CNPV187 by 91 bp.

The phylogenetic trees for fpv139 (Fig. 4Down) and fpv141/140 (data not shown) are essentially the same as each other and as those seen for clades A and B with P4b. However, this locus also provides additional resolution within subclade A2, allowing us for the first time to discriminate between the turkey and the pigeon isolates from within that subclade. Bootstrap analysis indicates that the pigeon isolate available to us from that subclade, PGPVP, is now clearly distinct from the turkey isolates remaining in subclade A2, clustering either with (but proximal to) subclade A3 (in 55–60 % of trees) or proximal to both subclades A2 and A3 (in 40–45 % of trees).


Figure 4
View larger version (10K):
[in this window]
[in a new window]
 
Fig. 4. Phylogenetic analysis of avipoxvirus and MOCV orthologues of fpv139. Abbreviations are as in Tables 1 and 2UpUp. Neighbour-joining phylogram of DNA sequences (480 bp) from avipoxvirus and MOCV genes, rooted on MOCV. All avipoxvirus sequences, except FP9 and CNPVATCC, were derived as part of this study. Avipoxvirus clades A1, A2, A3, B1 and B2 are labelled. Bootstrap values (1000 replicates) of >75 % are shown, as is the value (in parentheses) for PGPVP. It was found as shown in a node with, but proximal to, subclade A3 in 55 % of trees or, alternatively, proximal to both subclades A2 and A3 in 45 % of trees.

 
Clade A.
At this locus, viruses in subclade A1 all show 100 % amino acid identity to each other and 97–98 % amino acid identity to viruses in subclade A2. Within subclade A2, the TKPV isolates (66 and 98) are identical to each other, with 97.9 and 96.5 % amino acid identity to ABPV and PGPVP isolates, respectively. The latter show 97.4 % identity to each other.

Clade B.
At this locus, viruses in subclade B1 show 86–87 % amino acid identity to those in clade A. Within subclade B1, divergence between viruses isolated from different species (canary, houbara and sparrow) is barely more than that observed within CNPV.


   DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
The phylogenetics and epidemiology of few large groups of viruses are understood as poorly as are those of the avian poxviruses. Yet this group of viruses has, by virtue of the lifestyle of its hosts, an ability to cross natural and political boundaries, reaching new hosts as their migration and habitation patterns change. This feature of avian viruses is well illustrated by recent developments with West Nile virus and avian influenza virus strain H5N1.

Without knowing more about the avipoxviruses, it is difficult to assess the threat to reintroduction programmes, such as those for the great bustard on Salisbury Plain, Wiltshire, UK. The risks involved with reintroduction programmes in the Middle East for the closely related houbara bustard have been recognized (Samour et al., 1996Down; Bailey et al., 2002Down; Smits et al., 2005Down).

Most of our knowledge concerning the avipoxviruses is based on clinical isolates from diseased birds. The infectious history of such isolates is rarely clear. It is important, therefore, that large-scale surveys, such as that performed in the Canary Islands (Huw Lee & Hwa Lee, 1997Down), are also conducted to understand the nature of viruses carried, sometimes as subclinical infections, by natural populations. Comparison at the molecular level of viruses from these two situations could then help to identify novel threats faced by birds, whether captive, farmed or free-living, in response to natural or man-made changes to their habitats and avian contacts. Unfortunately, the divergence among avipoxviruses is, as reinforced by this study, so great that it is difficult to identify PCR primers, whether pan-genus or species-specific. As a consequence, the only PCR locus used until now has been the P4b locus.

It is in this context that this study sought to obtain a better phylogenetic framework for the avipoxviruses, by surveying more viruses with the previously used P4b locus and to identify new loci, such as H3L, that would permit a more robust classification.

P4b gene (fpv167) locus
PCR analysis of the highly conserved P4b gene was initially used as a diagnostic marker for FWPV infections, generating a 578 bp product (Tadese & Reed, 2003Down). Amplification of an overlapping 1409 bp product (by changing one primer) allowed subsequent discrimination of FWPV from PGPV and TKPV by EcoRV restriction-enzyme digestion profile (Luschow et al., 2004Down). MseI or EcoRV digestion of the original 578 bp P4b PCR product identified six different profiles, broadly discriminating between FWPV-like viruses, CNPV-like viruses, SRPV, PGPV, FLPV and AGPV. After sequence and phylogenetic analysis, the viruses were clustered into three groups equivalent to our subclades A1, A2 and B1, with two more single-virus branches, equivalent to subclade A4 and clade C.

A problem with the host species-based approach to the taxonomy of avipoxviruses is that sequences of isolates taken from particular species can be found in different subclades (e.g. FLPV) or even different clades (e.g. PGPV). Luschow et al. (2004)Down had already reported that one of two SRPVs that they had sequenced (SRPV32002) grouped with CNPV-like viruses (clade B), but that the other (SRPVDD1258) grouped with FWPV-like viruses (clade A). Two others, examined only by restriction-enzyme analysis, were CNPV-like. Both of our SRPV isolates (SRPV23 and SRPV9037), as well as one from another study (SRPV32002; Weli et al., 2004Down), also cluster with CNPV-like viruses. SRPVDD1258 is actually identical to all of the FWPV strains at the P4b locus. As the majority of SRPV isolates are therefore CNPV-like, it is not clear whether SRPVDD1258 represents (i) a rarer form of the disease in sparrows, (ii) a sporadic species jump, perhaps into an already sick bird, or (iii) an emerging disease.

The phylogenetic analysis conducted by Luschow et al. (2004)Down included only one PGPV (PGPVTP2), placing it in a cluster with OSPV, although restriction-enzyme analysis would place two more PGPV isolates in the same cluster (subclade A2). Weli et al. (2004)Down, however, obtained the sequence of a PGPV isolate (B7) that clusters with CNPV-like viruses in clade B. Of the two PGPV clinical isolates analysed in this study, one (Peekham) falls with PGPVTP2 in subclade A2 and the other (950) in subclade B2 with PGPVB7. It is therefore apparent that viruses from clades A and B can infect and cause disease in pigeons. It is not clear whether there are any significant differences in the disease caused by the two different types of viruses.

In this study, both a 1962 UK TKPV isolate and a 1998 Italian TKPV isolate cluster not with FWPV in subclade A1, but with PGPVTP2 and PGPVP, as well as with OSPV, in subclade A2. However, by examining more recent isolates, Luschow et al. (2004)Down showed that a 2001 TKPV isolate from Germany clusters with FWPV (in subclade A1). By restriction-enzyme analysis, two more German TKPV isolates (from 1999 and 2000) cluster similarly. This therefore reinforces the changing picture of poxvirus infections in turkeys, with FWPV as a recently emerging pathogen of turkeys.

A feature of the host-based approach to the taxonomy of avipoxviruses is that the host genus or species in which a virus causes severe disease and, therefore, from which it is most likely to be isolated need not be the normal reservoir host. Such situations may be common amongst passeriform birds infected by clade B CNPV-like viruses. Indeed, canaries, in which poxviruses cause devastating disease, may not be the natural reservoir for CNPV. In support of this perspective, a survey of wild birds in the Canary Islands found that 50 % of short-toed larks (Calandrella rufescens) and 28 % of Bertholet's pipits (Anthus berthelotti) had mild poxvirus lesions (Smits et al., 2005Down).

Virus isolates from sparrow and stone curlew cluster with CNPV (subclade B1), which correlates with all of the hosts for this subclade belonging to the order of passeriform birds. HOPV is part of the same viral phylogenetic cluster and, taxonomically, houbara bustard (order Gruiformes) is grouped with the order Passeriformes. PGPV clusters with viruses in the FWPV-like cluster (subclade A1), even though pigeon (order Columbiformes) is classified separately from turkey and chicken (order Galliformes) as well as ostrich (order Struthioniformes). However, another group of PGPVs clusters loosely with the order Passeriformes, in subclade B2. Therefore, the evolutionary taxonomy of the host does not appear to be the primary factor in driving avipoxvirus evolution, in so far as samples and sequences used within this study can show.

In this context, it is interesting to note a recent report (Adams et al., 2005Down) describing the closely related P4b sequences from a wide range of different birds presenting to a wildlife centre in Virginia, USA, in 2003–2004. The viruses form a tight cluster resembling our subclade B1, but distinct from it and B2, essentially forming a new subclade B3. The authors admit the possibility of cross-infection in the clinic, but, even excluding those birds that developed lesions during hospitalization, the hosts included birds as diverse as cardinal, blue-grey gnatcatcher, mockingbird, house finch and great blue heron. It will be interesting to see whether the observation is confirmed by other centres, showing whether the viruses are more widespread within the USA, whether they can be found in earlier epornitics and whether they persist over time.

H3L gene (fpv140) locus
The fpv140 locus PCR product allows easy distinction between clades A and B based solely on the size of the primary product, without further restriction-enzyme digestion. Sequence analysis of these products also allows better discrimination between viruses within subclade A2, in particular between the isolates from turkey and pigeon that lie within this subclade.

These primers can therefore be used as diagnostic markers to distinguish between FWPV-like and CNPV-like isolates by PCR and to differentiate between isolates by sequence analysis. They provide a useful addition to the currently used primers for the gene encoding P4b orthologues.

Unfortunately, the marker appeared to fail with clade C viruses and with the one subclade A3 virus available to us. There is, therefore, still a need for other pan-genus markers, such as the P4b locus. The difficulties of finding such a marker for genomes of this size and diversity, in the absence of more extensive genome-sequence data, should not be underestimated.

Implications for phylogeny and systematics of the poxviruses
Not all currently recognized avipoxvirus species are represented in the current tree: Juncopox virus, Mynahpox virus and QUPV (as well as tentative species Crowpox virus, Peacockpox virus and Penguinpox virus) remain to be analysed. However, a recent study (Hsieh et al., 2005Down) demonstrated that an isolate from Chinese jungle mynah (Acridotheres cristatellus), the sequence of which is not available in the databases, grouped with that from a wood pigeon (GenBank accession no. AY453177 [GenBank] , described as PGPVB7 in this manuscript). This is interesting, and consistent with the phylogenetic structure reported here, because our starling isolate (SLPV) also groups with PGPVB7 and because transmission of enzootic poxvirus from starlings to mynahs has been reported (Landolt & Kocan, 1976Down).

It may be premature to suggest changes to the nomenclature of avipoxviruses specifically and of poxviruses generally, particularly as most data are available for only the P4b gene/fpv168 locus (although more limited data for the H3L gene/fpv140 locus are consistent). However, it is appropriate to consider the existing taxonomic structure (Fauquet et al., 2005Down) and to question whether it can accommodate the emerging data without reorganization.

As an interim approach to nomenclature, we recommend reducing the number of recognized species to correspond with the phylogenetic groups. These could be FWPV (subclade A1), TKPV (subclade A2), CNPV (subclade B1), SLPV (subclade B2) and PSPV (clade C). PGPV would not be used, as it remains ambiguous until it can be established whether the existing taxonomic unit represents viruses in subclades A2 or B2. SRPV would be subsumed by CNPV and Mynahpox virus by SLPV.

It is relevant to any longer-term changes in nomenclature that the avipoxviruses are currently represented by only one genus (Avipoxvirus), whereas the mammalian poxviruses are represented by seven genera. It is clear that, for the P4b orthologues, the divergence between the avipoxvirus clades A, B and C is comparable to that seen between the existing genera of mammalian poxviruses.

The family Poxviridae is divided into subfamilies of insect (Entomopoxvirinae) and vertebrate (Chordopoxvirinae) viruses. This division offers little scope for taxonomic nomenclature that would both distinguish mammalian poxviruses from avipoxviruses and allow systematic groupings of the avipoxviruses. Two alternative approaches would offer such scope. The first would envisage replacing the two existing subfamilies with three: for the insecta, the diapsida (or archosauria) and the synapsida (or mammalia). The other alternative would envisage elevating the family Poxviridae to order status and the existing subfamily divisions to family level. This would allow groupings for the insecta, the diapsida (or archosauria) and the synapsida (or mammalia) at the subfamily level.

A third approach would be to maintain the existing nomenclature for family and subfamilies, but to replace the existing genus Avipoxvirus with three (or more) genera based on the major clades, with appropriate nomenclature. The loss of the genus Avipoxvirus might be mitigated by the use of genus names that incorporate ‘avipoxvirus’, such as: Falconiform avipoxvirus (with or without the space), Psittacine avipoxvirus, Passeriform avipoxvirus, Galliform avipoxvirus or Collumbiform avipoxvirus. Such genus names would, of course, also be compatible with the earlier two approaches.

We hope that the data presented here will encourage the acquisition of sequence data from a wider range of enzootic and epizootic isolates of this fascinating and neglected group of viruses. This will allow a rational and effective approach to more extensive genome sequencing of carefully selected members of the group, bringing more insight into the complex relationship between avian viruses and their hosts. We urge those responsible for collecting isolates to be as rigorous as possible in recording details of the host (genus and species, if and where possible supported by archiving of a small sample of host tissue for future typing), relevant epidemiology (date, location, whether captive or wild, any known disease contacts), as well as virus isolation and passage history.


   ACKNOWLEDGEMENTS
 
This work was funded by the Biotechology and Biological Sciences Research Council.


   REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES
 
Adams, C. J., Feldman, S. H. & Sleeman, J. M. (2005). Phylogenetic analysis of avian poxviruses among free-ranging birds of Virginia. Avian Dis 49, 601–605.[Medline]

Afonso, C. L., Tulman, E. R., Lu, Z., Zsak, L., Kutish, G. F. & Rock, D. L. (2000). The genome of fowlpox virus. J Virol 74, 3815–3831.[Abstract/Free Full Text]

Amano, H., Morikawa, S., Shimizu, H., Shoji, I., Kurosawa, D., Matsuura, Y., Miyamura, T. & Ueda, Y. (1999). Identification of the canarypox virus thymidine kinase gene and insertion of foreign genes. Virology 256, 280–290.[CrossRef][Medline]

Bailey, T. A., Silvanose, C., Manvell, R., Gough, R. E., Kinne, J., Combreau, O. & Launay, F. (2002). Medical dilemmas associated with rehabilitating confiscated houbara bustards (Chlamydotis undulata macqueenii) after avian pox and paramyxovirus type 1 infection. J Wildl Dis 38, 518–532.[Abstract]

Binns, M. M., Boursnell, M. E. G., Tomley, F. M. & Campbell, J. (1989). Analysis of the fowlpoxvirus gene encoding the 4b core polypeptide and demonstration that it possesses efficient promoter sequences. Virology 170, 288–291.[CrossRef][Medline]

Bolte, A. L., Meurer, J. & Kaleta, E. F. (1999). Avian host spectrum of avipoxviruses. Avian Pathol 28, 415–432.[CrossRef]

Boosinger, T. R., Winterfield, R. W., Feldman, D. S. & Dhillon, A. S. (1982). Psittacine pox virus: virus isolation and identification, transmission, and cross-challenge studies in parrots and chickens. Avian Dis 26, 437–444.[CrossRef][Medline]

Boulanger, D., Green, P., Smith, T., Czerny, C.-P. & Skinner, M. A. (1998). The 131-amino-acid repeat region of the essential 39-kilodalton core protein of fowlpox virus FP9, equivalent to vaccinia virus A4L protein, is nonessential and highly immunogenic. J Virol 72, 170–179.[Abstract/Free Full Text]

Boulanger, D., Green, P., Jones, B., Henriquet, G., Hunt, L. G., Laidlaw, S. M., Monaghan, P. & Skinner, M. A. (2002). Identification and characterization of three immunodominant structural proteins of fowlpox virus. J Virol 76, 9844–9855.[Abstract/Free Full Text]

Fauquet, C. M., Mayo, M. A., Maniloff, J., Desselberger, U. & Ball, L. A. (editors) (2005). Virus Taxonomy: Eighth Report of the International Committee on Taxonomy of Viruses. San Diego: Elsevier/Academic Press.

Gubser, C., Hué, S., Kellam, P. & Smith, G. L. (2004). Poxvirus genomes: a phylogenetic analysis. J Gen Virol 85, 105–117.[Abstract/Free Full Text]

Hall, B. G. (2001). Phylogenetic Trees Made Easy: A How-To Manual for Molecular Biologists. Sunderland, MA: Sinauer Associates.

Hsieh, Y. C., Chen, S. H., Wang, C. W., Lee, Y. F., Chung, W. C., Tsai, M. C., Chang, T. C., Lien, Y. Y. & Tsai, S. S. (2005). Unusual pox lesions found in Chinese jungle mynahs (Acridotheres cristatellus). Avian Pathol 34, 415–417.[CrossRef][Medline]

Huw Lee, L. & Hwa Lee, K. (1997). Application of the polymerase chain reaction for the diagnosis of fowl poxvirus infection. J Virol Methods 63, 113–119.[CrossRef][Medline]

Johnson, B. J. & Castro, A. E. (1986). Canary pox causing high mortality in an aviary. J Am Vet Med Assoc 189, 1345–1347.[Medline]

Kim, T. J., Schnitzlein, W. M., McAloose, D., Pessier, A. P. & Tripathy, D. N. (2003). Characterization of an avianpox virus isolated from an Andean condor (Vultur gryphus). Vet Microbiol 96, 237–246.[CrossRef][Medline]

Laidlaw, S. M. & Skinner, M. A. (2004). Comparison of the genome sequence of FP9, an attenuated, tissue culture-adapted European strain of Fowlpox virus, with those of virulent American and European viruses. J Gen Virol 85, 305–322.[Abstract/Free Full Text]

Landolt, M. & Kocan, R. M. (1976). Transmission of avian pox from starlings to Rothchild's mynahs. J Wildl Dis 12, 353–356.[Abstract/Free Full Text]

Luschow, D., Hoffmann, T. & Hafez, H. M. (2004). Differentiation of avian poxvirus strains on the basis of nucleotide sequences of 4b gene fragment. Avian Dis 48, 453–462.[Medline]

McLysaght, A., Baldi, P. F. & Gaut, B. S. (2003). Extensive gene gain associated with adaptive evolution of poxviruses. Proc Natl Acad Sci U S A 100, 15655–15660.[Abstract/Free Full Text]

Rambaut, A. (1996). Se-Al: Sequence Alignment Editor. Zoology Department, University of Oxford, UK.

Samour, J. H., Kaaden, O. R., Wernery, U. & Bailey, T. A. (1996). An epornitic of avian pox in houbara bustards (Chlamydotis undulata macqueenii). Zentralbl Veterinarmed B 43, 287–292.[Medline]

Smits, J. E., Tella, J. L., Carrete, M., Serrano, D. & López, G. (2005). An epizootic of avian pox in endemic short-toed larks (Calandrella rufescens) and Berthelot's pipits (Anthus berthelotti) in the Canary Islands, Spain. Vet Pathol 42, 59–65.[Abstract/Free Full Text]

Staden, R. (1982). Automation of the computer handling of gel reading data produced by the shotgun method of DNA sequencing. Nucleic Acids Res 10, 4731–4751.[Abstract/Free Full Text]

Stone, M. & Forbes, V. (2002–2003). MAF keen to hear about parrot disease, Psittacine poxvirus. In Protect, pp. 1–2. Wellington: MAF New Zealand. http://www.biosecurity.org.nz/files/Parrot_Disease.pdf

Tadese, T. & Reed, W. M. (2003). Detection of specific reticuloendotheliosis virus sequence and protein from REV-integrated fowlpox virus strains. J Virol Methods 110, 99–104.[CrossRef][Medline]

Thompson, J. D., Gibson, T. J., Plewniak, F., Jeanmougin, F. & Higgins, D. G. (1997). The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25, 4876–4882.[Abstract/Free Full Text]

Tripathy, D. N. (1984). Pox. In Diseases in Poultry, 9th edn, pp. 583–596. Edited by B. W. Calnek. Ames, IA: Iowa State Press.

Tripathy, D. N., Schnitzlein, W. M., Morris, P. J., Janssen, D. L., Zuba, J. K., Massey, G. & Atkinson, C. T. (2000). Characterization of poxviruses from forest birds in Hawaii. J Wildl Dis 36, 225–230.[Abstract]

Tulman, E. R., Afonso, C. L., Lu, Z., Zsak, L., Kutish, G. F. & Rock, D. L. (2004). The genome of canarypox virus. J Virol 78, 353–366.[Abstract/Free Full Text]

Weli, S. C., Traavik, T., Tryland, M., Coucheron, D. H. & Nilssen, Ø. (2004). Analysis and comparison of the 4b core protein gene of avipoxviruses from wild birds: evidence for interspecies spatial phylogenetic variation. Arch Virol 149, 2035–2046.[Medline]

Winterfield, R. W. & Reed, W. (1985). Avian pox: infection and immunity with quail, psittacine, fowl, and pigeon pox viruses. Poult Sci 64, 65–70.[Medline]

Received 6 December 2005; accepted 18 April 2006.


This article has been cited by other articles:


Home page
J. Gen. Virol.Home page
S. A. Jarmin, R. Manvell, R. E. Gough, S. M. Laidlaw, and M. A. Skinner
Retention of 1.2 kbp of 'novel' genomic sequence in two European field isolates and some vaccine strains of Fowlpox virus extends open reading frame fpv241
J. Gen. Virol., December 1, 2006; 87(12): 3545 - 3549.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Jarmin, S.
Right arrow Articles by Skinner, M. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Jarmin, S.
Right arrow Articles by Skinner, M. A.
Agricola
Right arrow Articles by Jarmin, S.
Right arrow Articles by Skinner, M. A.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS