|
|
||||||||
Centre for Immunology and Cancer Research, University of Queensland, Princess Alexandra Hospital, Brisbane, QLD 4102, Australia
Correspondence
Annika Antonsson
aantonsson{at}cicr.uq.edu.au
| ABSTRACT |
|---|
|
|
|---|
The GenBank/EMBL/DDBJ accession number for the sequence reported in this paper is DQ464069.
| INTRODUCTION |
|---|
|
|
|---|
Previous reports (Antonsson et al., 2000
, 2003a
, b
) have shown that HPV is commonly present in the healthy skin of healthy humans as a presumed commensal agent. Recently, we have also shown that skin papillomaviruses can readily be detected in healthy skin from many different animal species (Antonsson & Hansson, 2002
). Furthermore, these animal papillomavirus types were sufficiently genetically related to their human counterparts to be identifiable by a human skin papillomavirus primer set (FAP59 and FAP64) (Forslund et al., 1999
). Given that papillomaviruses appear to be so widespread, we wished to investigate the distribution of papillomaviruses in Australian fauna. Here we analysed skin swab samples from 54 different Australian animal species that have never before been reported to harbour papillomavirus in healthy skin.
| METHODS |
|---|
|
|
|---|
This project was approved by the Animal Ethics Committee, University of Queensland (CICR 158/05 and CICR 297/06).
Samples.
Swab samples were collected with moist (0.9 % NaCl solution) cotton-tipped swabs (Sarstedt), which were drawn back and forth five times over the forehead skin or fur within an area of 4 by 8 cm, and then suspended in 1 ml 0.9 % NaCl solution. No DNA extraction was performed prior to PCR. Samples were kept at 4 °C for a maximum of 72 h before being analysed.
Three of the wart-like lesions were stored at 20 °C (without 0.9 % NaCl) until analysed. DNA was then extracted using the DNeasy tissue kit (Qiagen). The other three wart-like lesions were put in formalin and were used for histopathology analysis.
PCR and HPV type determination.
A PCR test with the primer pair FAP59/FAP64 was used for detection of skin HPV DNA, as previously described (Forslund et al., 1999
). The protocol was followed as described except for the MgCl2 concentration, which was modified to 3.5 mM. HPV 11 (plasmid) was used as a positive control. The PCR products were cloned into the pCR-script SK(+) cloning vector (Stratagene). Between two and four clones per sample were sequenced with both forward and reverse primers (BigDye Terminator v3.1 cycle sequencing kit; Applied Biosystems), and the sequences obtained were compared with available sequences in the GenBank database using the BLAST server (http://www.ncbi.nlm.nih.gov/blast/).
The ability of the FAP59/FAP64 PCR to detect skin papillomavirus from the animal swab samples was compared with that of an animal papillomavirus PCR described by Rector et al. (2005b)
. Briefly, a 5 µl sample was used in the PCR with the AR-L1F1/AR-L1R3 primer pair, which yielded a PCR product of 600 bp. The COPV plasmid was used as a positive control and H2O without template as a negative control. The bands yielded from the AR-L1 PCR were cut out from the gel and purified on a spin column (Perfectprep gel cleanup kit; Eppendorf), eluted in 30 µl H2O, then cloned and sequenced as described above.
A new type-specific primer pair were designed for the new putative koala PV type detected. The new primers were designed using the FAP59/FAP64 primers and the new putative koala papillomavirus type as templates, which gave the following primer sequences: KoPVFAP59 (forward primer) 5'-TAACTGTGGGGCATCCATATT-3' and KoPVFAP64 (reverse primer) 5'-GATGGCGACATGTCTGATATAGG-3'. The same protocol as described above for FAP59/FAP64 was followed.
Since full-length L1 sequences were not obtained, the new PV isolate detected in this study, which was cloned and submitted to GenBank (accession no. DQ464069
[GenBank]
), was designated a putative PV type. The guidelines from the Papillomavirus Nomenclature Committee 1995 (14th International Papillomavirus Conference, Quebec City, Quebec, Canada) for defining a new putative HPV type (de Villiers, 2001
) were followed. This new putative animal PV type has been designated by an initial letter combination for the species it was found in, i.e. Ko for koala, followed by AA and a unique number (KoAA1), as for previously identified putative animal papillomavirus types (Antonsson & Hansson, 2002
). The papillomavirus isolates identified in this paper by direct sequencing were designated potential PV types.
Phylogenetic analysis.
Phylogenetic analysis was based on multiple alignment with ClustalX (version 1.8) (Jeanmougin et al., 1998
), and the alignments were edited with Genedoc (version 2.4.000) (Nicholas et al., 1997
). PHYLIP (version 3.5) (Felsenstein, 1982
, 1985
) was used for neighbour-joining and maximum-likelihood analyses. These programs were obtained from the website http://evolution.genetics.washington.edu/phylip/software.html.
The taxonomic system with papillomavirus genus groups
to
was used (de Villiers et al., 2004
). This system was applied to demonstrate relatedness between the PV type candidates and previously established PV types.
The region of the L1 gene amplified by the primer pair FAP59/FAP64 was used for the phylogenetic analysis, extending from nucleotide 6044 to 6480 relative to the HPV20 sequence. The following papillomaviruses were included in the phylogenetic analysis (GenBank accession nos, in parentheses): HPV 1 (V01116 [GenBank] ), HPV 4 (X70827 [GenBank] ), HPV 5 (M17463 [GenBank] ), HPV 8 (M12737 [GenBank] ), HPV 9 (X74464 [GenBank] ), HPV 12 (X74466 [GenBank] ), HPV 14 (X74467 [GenBank] ), HPV 15 (X74468 [GenBank] ), HPV 17 (X74469 [GenBank] ), HPV 19 (X74470 [GenBank] ), HPV 20 (U31778 [GenBank] ), HPV 21 (U31779 [GenBank] ), HPV 22 (U31780 [GenBank] ), HPV 23 (U31781 [GenBank] ), HPV 24 (NC_001683 [GenBank] ), HPV 25 (X74471 [GenBank] ), HPV 36 (U31785 [GenBank] ), HPV 37 (U31786 [GenBank] ), HPV 38 (U31787 [GenBank] ), HPV 41 (X56147 [GenBank] ), HPV 48 (U31789 [GenBank] ), HPV 49 (X74480 [GenBank] ), HPV 50 (U31790 [GenBank] ), HPV 60 (U31792 [GenBank] ), HPV 63 (X70828 [GenBank] ), HPV 65 (X70829 [GenBank] ), HPV 75 (Y15173 [GenBank] ), HPV 76 (Y15174 [GenBank] ), HPV 80 (Y15176 [GenBank] ), HPV 92 (AF531420 [GenBank] ), HPV 93 (AY382778 [GenBank] ), HPV 95 (AJ620210 [GenBank] ) and HPV 96 (AY382779 [GenBank] ), colobus monkey PV type 2 (CgPV 2; U72630 [GenBank] ), rhesus monkey PV type 1 (RhPV 1; M60184 [GenBank] ), chimpanzee PV type 1 (CPV 1; AF020905 [GenBank] ), pygmy chimpanzee PV type 1 (PCPV 1; X62844 [GenBank] ), bovine (cattle) PV type 1 (BPV 1; X02346 [GenBank] ), BPV 2 (M20219 [GenBank] ), BPV 3 (AF486184 [GenBank] ), BPV 4 (X05817 [GenBank] ), BPV 5 (AF457465 [GenBank] ), BPV 6 (AJ620208 [GenBank] ), ovine (sheep) PV type 1 (OvPV 1; U83594 [GenBank] ), OvPV 2 (U83585 [GenBank] ), European elk PV (EEPV; M15953 [GenBank] ), deer PV (DPV; M11910 [GenBank] ), canine (dog) oral PV (COPV; L22695 [GenBank] ), canine PV type 2 (CPV 2; AY722648 [GenBank] ), rabbit oral PV (ROPV; AF227240 [GenBank] ), Felis domesticus (cat) PV (FdPV; AF480454 [GenBank] ), Equus caballus (horse) PV (EcPV; AF498323 [GenBank] ), Trichechus manatus latirostris (Florida manatee) PV (TmPV; NC_006563 [GenBank] ), Erethizon dorsatum (North American porcupine) PV (EdPV; NC_006951 [GenBank] ), Procyon lotor (raccoon) PV (PlPV; NC_007150 [GenBank] ), Psittacus erithacus timneh (African grey parrot) PV (PePV; NC_003973 [GenBank] ), Fringilla coelebs (chaffinch) PV (FcPV; AY057109 [GenBank] ), hamster oral PV (HOPV; E15110 [GenBank] ), Reindeer PV (RPV; AF443292 [GenBank] ), Phocoena spinipinnis (Burmeister's porpoise) PV (PsPV; AJ238373 [GenBank] ) and Mastomys natalensis (South African mouse) PV (MnPV; U01834 [GenBank] ). Also included in the analysis was the putative new koala PV type KoAA1 (DQ464069 [GenBank] ). The DNA sequences from previously characterized PV types were obtained from GenBank (http://www.ncbi.nlm.nih.gov/).
The potential PV types described in this paper were not included in the phylogenetic analysis as they were only sequenced once with either FAP59 or FAP64.
| RESULTS |
|---|
|
|
|---|
|
|
PCR with the primer pair FAP59/FAP64 did not identify any PV DNA in the samples from the following animal species: red kangaroo, wallaby, possum, wombat, dingo, emu, quoll, great glider, black flying fox, grey-headed flying fox, salt-water crocodile, bearded dragon, shingleback lizard, carpet python, olive python, coastal carpet python, Darwin carpet python, scrub python, spotted python, Children's python, water python, frog, kookaburra, cockatoo, frogmouth, pigeon, purple-crowned lorikeet, eclectus parrot, superb fruit dove, brown falcon, crested hawk, sealion, ray, blue tang fish, shark, turtle, octopus or starfish.
An alternative animal papillomavirus PCR pair, AR-L1F1/AR-L1R3, as described by Rector et al. (2005b)
, was tested but no papillomavirus DNA was detected in any of the 225 swab samples despite our positive control (COPV plasmid) yielding an amplicon of the correct size.
Wart-like lesions from a koala
During sampling a number of wart-like lesions were noted on the chest of a 2-year-old male koala and three of these were removed for analysis. One of the three lesions tested positive for PV DNA and was identified as KoAA1. After histopathology analysis, the lesion was confirmed to be a sebaceous gland hyperplasia rather than a wart. These scent/sebaceous glands are present on the chest and are used for marking territory.
Phylogenetic analysis
Phylogenetic trees obtained by the neighbour-joining, maximum-parsimony and maximum-likelihood algorithms were compared. The tree in Fig. 1
is divided into the genera
to
and shows 34 fully characterized human skin PV types, 28 fully sequenced animal PV types, and the putative new koala PV type (KoAA1) described in the present study. The phylogenetic analysis places KoAA1 in its own branch in the phylogenetic tree. Furthermore, a phylogenetic analysis based on amino acid alignment was performed and the resulting tree was very much the same as that based on nucleic acid alignment.
|
| DISCUSSION |
|---|
|
|
|---|
Since Australia separated from the other continents about 95 million years, its wildlife has been able to evolve independently. Many of the animals of Australia are distinctly primitive, such as the monotremes that lay eggs and suckle their young, and the marsupials. We found putative and potential papillomavirus types in both a monotreme species (echidna) and marsupials (koala and Eastern grey kangaroo). Only one of the four papillomavirus types could be cloned and sequenced, and hence included in a phylogenetic analysis. The other three sequences were obtained by direct sequencing and were not analysed phylogenetically. Interestingly, the new putative koala PV type, KoAA1, forms a new branch or even a putative new genus, close to the root in the phylogenetic tree. Only one putative papillomavirus type has been described in marsupials before and it was isolated from a skin lesion on the tail of a brush-tail possum (Perrott et al., 2000
). However, this PV type was detected with the MY09/MY11 primer pair, which does not overlap the PCR product of the FAP59/FAP64 primers (used here), meaning we could not include it in our phylogenetic analysis.
The alternative PCR did not detect any papillomavirus in our specimens. However, this primer pair was designed to detect PV DNA in dermal proliferations of racoons with high copy numbers of PV virions (Rector et al., 2005b
) so it is not surprising that PV could not be detected in our swab samples with presumably low number of virus particles.
This study shows that even the healthy skin of primitive mammals like monotremes and marsupials harbour subclinical PV infections, and further shows the ubiquity of skin PVs. Our findings suggest that the phylogenetic tree of the family Papillomaviridae is far from complete and this contention is supported by the recent characterization of a large number of new animal PVs that make up new genera of the growing phylogenetic tree (Ghim et al., 2004
; Rector et al., 2004
, 2005a
, b
; Tachezy et al., 2002a
, b
; Van Doorslaer et al., 2006
).
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
Antonsson, A., Forslund, O., Ekberg, H., Sterner, G. & Hansson, B. G. (2000). The ubiquity and impressive genomic diversity of human skin papillomaviruses suggest a commensalic nature of these viruses. J Virol 74, 1163611641.
Antonsson, A., Erfurt, C., Hazard, K., Holmgren, V., Simon, M., Kataoka, A., Hossain, S., Håkangård, C. & Hansson, B. G. (2003a). Prevalence and type spectrum of human papillomaviruses in healthy skin samples collected in three continents. J Gen Virol 84, 18811886.
Antonsson, A., Karanfilovska, S., Lindqvist, P. G. & Hansson, B. G. (2003b). General acquisition of human papillomavirus infections of skin occurs in early infancy. J Clin Microbiol 41, 25092514.
Bernard, H. U. (1994). Coevolution of papillomaviruses with human populations. Trends Microbiol 2, 140143.[CrossRef][Medline]
Chan, S. Y., Bernard, H. U., Ratterree, M., Birkebak, T. A., Faras, A. J. & Ostrow, R. S. (1997). Genomic diversity and evolution of papillomaviruses in rhesus monkeys. J Virol 71, 49384943.[Abstract]
de Villiers, E.-M. (2001). Taxonomic classification of papillomaviruses. Papillomavirus Rep 12, 5763.
de Villiers, E. M., Fauquet, C., Broker, T. R., Bernard, H. U. & zur Hausen, H. (2004). Classification of papillomaviruses. Virology 324, 1727.[CrossRef][Medline]
Felsenstein, J. (1982). Numerical methods for inferring evolutionary trees. Q Rev Biol 57, 379404.[CrossRef]
Felsenstein, J. (1985). Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39, 783791.[CrossRef]
Forslund, O., Antonsson, A., Nordin, P., Stenquist, B. & Hansson, B. G. (1999). A broad range of human papillomavirus types detected with a general PCR method suitable for analysis of cutaneous tumours and normal skin. J Gen Virol 80, 24372443.
Ghim, S. J., Rector, A., Delius, H., Sundberg, J. P., Jenson, A. B. & Van Ranst, M. (2004). Equine papillomavirus type 1: complete nucleotide sequence and characterization of recombinant virus-like particles composed of the EcPV-1 L1 major capsid protein. Biochem Biophys Res Commun 324, 11081115.[CrossRef][Medline]
Jeanmougin, F., Thompson, J. D., Gouy, M., Higgins, D. G. & Gibson, T. J. (1998). Multiple sequence alignment with Clustal X. Trends Biochem Sci 23, 403405.[CrossRef][Medline]
Nicholas, K. B., Nicholas, H. B., Jr & Deerfield, D. W. (1997). Genedoc: analysis and visualization of genetic variation. EMBNEW News 4, 14.
Perrott, M. R., Meers, J., Greening, G. E., Farmer, S. E., Lugton, I. W. & Wilks, C. R. (2000). A new papillomavirus of possums (Trichosurus vulpecula) associated with typical wart-like papillomas. Arch Virol 145, 12471255.[CrossRef][Medline]
Rector, A., Bossart, G. D., Ghim, S. J., Sundberg, J. P., Jenson, A. B. & Van Ranst, M. (2004). Characterization of a novel close-to-root papillomavirus from a Florida manatee by using multiply primed rolling-circle amplification: Trichechus manatus latirostris papillomavirus type 1. J Virol 78, 1269812702.
Rector, A., Tachezy, R., Van Doorslaer, K., MacNamara, T., Burk, R. D., Sundberg, J. P. & Van Ranst, M. (2005a). Isolation and cloning of a papillomavirus from a North American porcupine by using multiply primed rolling-circle amplification: the Erethizon dorsatum papillomavirus type 1. Virology 331, 449456.[CrossRef][Medline]
Rector, A., Van Doorslaer, K., Bertelsen, M., Barker, I. K., Olberg, R. A., Lemey, P., Sundberg, J. P. & Van Ranst, M. (2005b). Isolation and cloning of the raccoon (Procyon lotor) papillomavirus type 1 by using degenerate papillomavirus-specific primers. J Gen Virol 86, 20292033.
Shadan, F. F. & Villarreal, L. P. (1993). Coevolution of persistently infecting small DNA viruses and their hosts linked to host-interactive regulatory domains. Proc Natl Acad Sci U S A 90, 41174121.
Sundberg, J. P. (1987). Papillomavirus infections in animals. In Papillomaviruses and Human Disease, pp. 40103. Edited by K. Syrjänen, L. Gissmann & L. G. Koss. Berlin: Springer.
Tachezy, R., Duson, G., Rector, A., Jenson, A. B., Sundberg, J. P. & Van Ranst, M. (2002a). Cloning and genomic characterization of Felis domesticus papillomavirus type 1. Virology 301, 313321.[CrossRef][Medline]
Tachezy, R., Rector, A., Havelkova, M., Wollants, E., Fiten, P., Opdenakker, G., Jenson, B., Sundberg, J. & Van Ranst, M. (2002b). Avian papillomaviruses: the parrot Psittacus erithacus papillomavirus (PePV) genome has a unique organization of the early protein region and is phylogenetically related to the chaffinch papillomavirus. BMC Microbiol 2, 19.[CrossRef][Medline]
Van Doorslaer, K., Rector, A., Vos, P. & Van Ranst, M. (2006). Genetic characterization of the Capra hircus papillomavirus: a novel close-to-root artiodactyl papillomavirus. Virus Res 118, 164169.[CrossRef][Medline]
Received 11 May 2006;
accepted 27 July 2006.
This article has been cited by other articles:
![]() |
M. D. Bennett, L. Woolford, A. J. O'Hara, K. S. Warren, and P. K. Nicholls In situ hybridization to detect bandicoot papillomatosis carcinomatosis virus type 1 in biopsies from endangered western barred bandicoots (Perameles bougainville) J. Gen. Virol., February 1, 2008; 89(2): 419 - 423. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Woolford, A. Rector, M. Van Ranst, A. Ducki, M. D. Bennett, P. K. Nicholls, K. S. Warren, R. A. Swan, G. E. Wilcox, and A. J. O'Hara A Novel Virus Detected in Papillomas and Carcinomas of the Endangered Western Barred Bandicoot (Perameles bougainville) Exhibits Genomic Features of both the Papillomaviridae and Polyomaviridae J. Virol., December 15, 2007; 81(24): 13280 - 13290. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |