BioMed Research International

Volume 2013 (2013), Article ID 270898, 7 pages

http://dx.doi.org/10.1155/2013/270898

Research Article

Bovine Papillomavirus in Brazil: Detection of Coinfection of Unusual Types by a PCR-RFLP Method

R. F. Carvalho,¹ S. T. Sakata,1,2 D. N. S. Giovanni,¹ E. Mori,² P. E. Brandão,² L. J. Richtzenhain,² C. R. Pozzi,³ J. R. P. Arcaro,³ M. S. Miranda,³ J. Mazzuchelli-de-Souza,1,4 T. C. Melo,1,5 G. Comenale,¹ S. L. M. R. Assaf,¹ W. Beçak,1,6 and R. C. Stocco1,4

¹Laboratório de Genética, Instituto Butantan, Secretária de Estado da Saúde, Avenida Vital Brasil 1500, Butantan, 05503-900 São Paulo, SP, Brazil

²Laboratório de Biologia Molecular Aplicada e Sorologia, Departamento de Medicina Veterinária Preventiva e Saúde Animal, Faculdade de Medicina Veterinária, e Zootecnia da Universidade de São Paulo, Avenida Prof. Dr. Orlando Marques de Paiva, 87 Cidade Universitária, Butantan, 05508-270 São Paulo, SP, Brazil

³Centro de Pesquisa em Pecuária do Leite, Instituto de Zootecnia, Agencia Paulista de Tecnologia dos Agronegócios, CAPTA Bovinos Leiteiros, Rodovia Luiz de Queiroz km 129, 13460-000 Nova Odessa, SP, Brazil

⁴Programa de Pós-Graduação Interunidades em Biotecnologia, Instituto de Ciências Biomédicas, Edifício ICB-III, Universidade de São Paulo, Avenida Prof. Lineu Prestes, 2415 Cidade Universitária, Butantan, 05508-900 São Paulo, SP, Brazil

⁵Programa de Pós-Graduação em Biologia Estrutural e Funcional, Universidade Federal de São Paulo, Rua Botucatu, 740 Vila Clementino, 04023-900 São Paulo, SP, Brazil

⁶Departamento de Biologia, Universidade Federal da Integração Latino-Americana (UNILA), Avenida Tancredo Neves 6731, Bloco 4, 85867-970 Foz do Iguaçú, PR, Brazil

Received 26 March 2013; Accepted 16 May 2013

Academic Editor: Lubna Nasir

Copyright © 2013 R. F. Carvalho et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Bovine papillomavirus (BPV) is recognized as a causal agent of benign and malignant tumors in cattle. Thirteen types of BPV are currently characterized and classified into three distinct genera, associated with different pathological outcomes. The described BPV types as well as other putative ones have been demonstrated by molecular biology methods, mainly by the employment of degenerated PCR primers. Specifically, divergences in the nucleotide sequence of the L1 gene are useful for the identification and classification of new papillomavirus types. On the present work, a method based on the PCR-RFLP technique and DNA sequencing was evaluated as a screening tool, allowing for the detection of two relatively rare types of BPV in lesions samples from a six-year-old Holstein dairy cow, chronically affected with cutaneous papillomatosis. These findings point to the dissemination of BPVs with unclear pathogenic potential, since two relatively rare, new described BPV types, which were first characterized in Japan, were also detected in Brazil.

1. Introduction

The Bovine papillomavirus (BPV) is recognized as the causal agent of benign and malignant tumors in cattle, such as cutaneous papillomas, benign fibroplasias, urinary bladder, and esophagus cancer, causing significant economic losses. This oncogenic virus has a double-stranded circular DNA genome of approximately eight kilobases [1].

Currently, the Papillomaviridae family is divided into 16 genera according to their genomic organization [2, 3]. The papillomavirus (PV) genome codified functional, early (E) proteins and structural, late (L) proteins, expressed at different stages of the viral cycle. The L1 is the most conserved gene within a PV genome and has therefore been used for the identification of new PVs: one PV isolate is recognized as a new type if the complete genome has been cloned and the DNA sequence of L1 differs by more than 10% from the closest known PV type. Differences between 2% and 10% define a subtype and less than 2%, a variant [3]. Thirteen types of BPVs are currently well characterized and classified into three distinct genera—Delta, Epsilon, and Xi—each one associated with epithelia lesions of specific histological nature [4].

The BPVs-1 and -2 are classified as Delta papillomaviruses [5]. Characteristically, these types induce the appearance of fibropapillomas, associated with the recruitment of the subepithelial fibroblasts [6]. As far as concerned, both types are also unique in their ability to infect different host species, not only bovines, causing the equine sarcoid [7]. Lately, the genome of a new Delta-BPV type (BPV-13) was fully sequenced [8].

A larger number of BPV types (-3, -4, -6, -9, -10, -11 and -12) belong to the Xipapillomavirus genus. These viruses are considered exclusively epitheliotropic, inducing the formation of “true papillomas,” without the involvement of fibroblasts [9–11]. On the other hand, the BPVs-5 and -8 have the potential to induce both fibropapillomas and true papillomas in the course of their infectious cycle, being classified into a third genus, Epsilonpapillomavirus [12, 13]. The BPV-7 represents an exception and is classified separately (unsigned genus). This virus was first isolated from a cutaneous papilloma lesion and also from healthy teat skin samples [14].

The thirteen described BPV types as well as other putative ones have been demonstrated by molecular biology methods, since papillomaviruses are not prone to be replicated or recovered in cell cultures [15–19]. According to the guidelines outlined by the Papil omavirus Nomenclature Committee (14th International Papillomavirus Conference, Quebec City, QC, Canada), it has been specified that the amplified sequences isolated from novel papillomaviruses isolates could indicate only putative new PV types—instead of PV types—since the PCR amplicons represent only part of the L1 gene [15].

The use of PCR assays with degenerated primers, followed by sequencing has allowed the identification of several PV types in human and other animal hosts [15, 20]. The PCR primer FAP set was designed from two relatively conserved regions found in the L1 gene and has been shown to amplify PVs DNA from both papillomas and healthy tissue of many animal species, including BPVs in bovines [15, 16, 20].

Brazil has a cattle herd of approximately 210 million, being a major exporter of meat, milk and leather. BPVs have been previously detected in Brazil [21], but the extent of the impact of BPVs associated diseases, both in dairy and cattle herds, needs further studies. Available reports in different regions of the country indicate a significant diversity of viral types among the Brazilian herd, implying an evident disease burden [18, 22–25]. Thus, the improvement of knowledge concerning the diagnosis and related clinical aspects of different BPV types among the Brazilian herd should be considered in the development of new sanitary measures, aiming to the prevention of BPV infection and its consequences.

Unfortunately, BPV epidemiological surveys are still limited by the availability of high-throughput diagnostic techniques that could discriminate different BPV sequences at the same time in co infected samples [16, 26]. In this context, the present work represents an effort to identify BPV types employing an alternative screening method based on the PCR-RFLP technique and correlating the histological data of the analyzed lesions with the diagnosed viral type.

2. Material and Methods

In silico generation of RFLPs: the L1 FAP segment digestion profiles of the BPVs-1 to -13 could be generated with NEB cutter 2.0 [27] from all L1 complete nucleotide sequences available in Genbank (http://www.ncbi.nlm.nih.gov/genbank/). The restriction enzymes sites were chosen both by its presence (or absence) as well as the generated digestion fragments sizes, in order to differentiate those thirteen different BPV types. Histopathological analysis: wart biopsies were obtained from the trunk of a six-year-old Holstein dairy cow chronically affected with cutaneous papillomatosis. Samples from three different lesions were submitted to macroscopic, histological (hematoxylin and Eosin staining) and molecular analyses. DNA extraction and PCR: DNA was extracted from warts for viral typing (Illustra tissue and cells genomic Prep Mini Spin GE Healthcare) and an approximately 470, base-pairs L1 gene segment was amplified using the following primer sequences: forward: FAP59 (5′-TAA CWG TIG GIC AYC CWT ATT-3′) reverse: FAP64 (5′-CCWATATCWVHCATITCICCATC-3′). The PCR, were performed with slight modifications of a previously described protocol [16]. In detail the amplification reactions were performed in a Corbett CG1-96 thermocycler (Corbett Life Science, Sydney, Australia), with GoTaq Master Mix (Promega, Madison, USA), under the following conditions: 5 min at 95°C, followed by 35 cycles of 1 min and 30s at 95°C, 2 min at 52°C and 1 min and 30s at 72°C and a final extension step of 5 min at 72°C. Restriction analysis: an aliquot of PCR fragments were submitted to digestion reactions for RFLP analysis with four different restriction enzymes (DdeI, HinfI, HindIII, MsLI), following manufacturer’s instructions (New England Biolabs, Ipswich, USA). Cloned BPV-1 and BPV-2 genomes as well as a known typed clinical sample (Mg-19, BPV-2 typed) were used as positive controls. PCR-RFLP products were analyzed in 2.0% agarose gel electrophoresis stained with ethidium bromide (0.5 μg/mL) in TAE buffer and visualized under UV light. Sequencing: An aliquot of all generated PCR fragments were purified with extraction columns (Illustra GFX PCR DNA and Gel Band Purification Kit GE Healthcare). DNA concentration and purity were determined in a spectrophotometer (Eppendorf BioPhotometer, Hamburg, Germany) and submitted to sequencing reactions: three independent sequencing reactions were done for each PCR fragment in an ABI377 PRISM Genetic Analyzer (Life Applied Biosystems, USA). The quality of DNA sequences was checked and overlapping fragments were assembled using the BioEdit package 7.0.9.0 [28]. Assembled sequences with high quality were aligned using ClustalW 1.83 [29] with default gap penalties. Homology analyses were performed with the NCBI database and BLAST [30]. BioEdit software was used to identify the equivalent amino acid sequences. The sequence alignments were performed using the MEGA 5.0 software [31], using full alignment and 2000 total replications on the bootstrap, in order to ensure a higher level of confidence to our analysis [32]. Phylogenetic Analysis: phylogenetic relationship comparing nucleotide sequences was performed with MEGA. Neighbor-joining trees were drawn using TreeView version 1.6.6 [33]. Nucleotide and amino acid sequences from other BPV types and of a HPV-16 were retrieved from the GenBank (http://www.ncbi.nlm.nih.gov/) for comparison with the obtained sequences here.

Ethics Statement. The protocols used in this study were approved by the Ethical Committee for Animal Experimentation of the Instituto de Zootecnia (Protocol no. 109, on July 06th, 2009) assigned by the President of this Committee. All efforts were made to minimize animal suffering.

3. Results

In silico analyses of the restrictionsites in FAP L1 fragments of 20 deposited BPVs sequences (from BPV-1 to BPV-13) revealed no intratype variation associated with the relative cut positions for the four enzymes employed (Table 1). Macroscopic and histological evaluations identified the collected lesions as cutaneous fibrous papillomas (Figures 1 and 2). From all the three lesions studied, two- (IZ1 and IZ3) RFLP profiling suggests the presence of BPV-11. On the other hand, IZ2 sample gene rated a profile which could not match with any of the thirteen characterized viral types (Figure 3, Table 1). In accordance with these results, DNA sequencing and BLAST analysis of IZ2 indicated a rare putative type (BAPV-3), originally described in Japan [16]. Furthermore, the DNA sequencing of other two samples confirms them as the recently described BPV-11 [10]. These sequences were deposited in GenBank (access numbers: HQ435675 and HQ612180). The phylogenetic reconstruction using this partial genomic sequence allowed its comparison with other sequences of Papillomaviridae family members (Figure 4).

Table 1: L1 FAP segment digestion profiles of the BPVs-1 to -13 for Dde  I, Hinf I, MsL Iand HindIII.

L1 FAP Fragment
DdeI
HinfI
MsLI
HindIII

BPV-1
475 bp
264
329
301
475

159
146
174

52

BPV-2
475 bp
316
329
475
475

159
146

BPV-3
473 bp
319
473
373
473

154

87

13

BPV-4
469 bp
324
469
469
230

145

182

87

BPV-5
469 bp
469
254
469
469

145

70

BPV-6
472 bp
321
384
373
472

151
88
99

BPV-7
484 bp
403
271
484
484

81
119

94

BPV-8
469 bp
469
317
469
322

152

147

BPV-9
469 bp
316
344
469
469

153
125

BPV-10
472 bp
319
381
373
472

90
91
99

63

BPV-11
475 bp
420
407
370
475

55
68
105

BPV-12
469 bp
351
317
469
469

118
152

BPV-13
475 bp
316
329
475
475

109
146

50

Reference Sequences:

BPV-1_NC_001522.1, BPV-1_X02346.1, BPV-2_M20219.1, BPV-2_X01768.1,

BPV-3_AF486184.1, BPV-3_AJ620207.1, BPV-3_NC_004197.1, BPV-4_X05817.1,

BPV-5_AF457465.1, BPV-5_NC_004195.1, BPV-6_AJ620208.1, BPV-7_DQ217793.1,

BPV-7_NC_007612.1, BPV-8_EB_DQ098917.1, BPV-8_NC_009752.1, BPV-9_AB331650.1,

BPV-10_AB331651.1, BPV-11_AB543507, BPV-12_JF834523, BPV-13_JQ798171.

(a)

(b)

(c)

Figure 1: (a) Analyzed animal IZ-1214 with disseminated cutaneous papillomatosis. (b) Gross aspect of the papilloma lesions with a cauliflower-like appearance. (c) Histological preparation (hematoxylin and eosin staining, or HE) of the collected lesions (IZ1, IZ2, and IZ3) indicating characteristic hyperkeratosis, acanthosis, and papillomatosis with dermis proliferation in HE preparation (100x).

(a)

(b)

Figure 2: (a) Histopathology of a wart biopsy: detailed aspect of the IZ2 lesion exhibiting characteristic hyperkeratosis, acanthosis and dermal proliferation, indicated by arrows (100x). (b) Presence of koilocytosis.

(a)

(b)

(c)

(d)

Figure 3: RFLP performed for restriction enzymes DdeI, HinfI, MsLI, and HindIII. As BPV-1 and BPV-2 are the commonly associated with fibropapillomas, their L1 FAP amplicons were generated and digested as positive controls. IZ1, IZ2, IZ3 amplicons are from three different lesions of the same animal (IZ-1214). Mg-19 sample is from another bovine, clinically affected with cutaneous papillomatosis. Molecular marker: 100 bp ladder (New England Biolabs, Ipswich, UK).

Figure 4: Neighbor-joining phylogenic tree constructed with 2000 bootstrap replications, using partial L1 nucleotide sequences, indicating the putative type BAPV-3 in IZ2 sample, originally described in Japan (accession: AY300819). The analysis of the other two samples, IZ1 and IZ3, demonstrated the presence of an identical sequence in both lesions, sharing 99% percent of similarity with the recently described BPV-11 type (accession: AB543507.1). As IZ3 sequence is identical to IZ1, it will not be shown here.

4. Discussion

Originally, the recently characterized BPV-11 was described with the employment of a Xipapillomavirus consensus primer [10]. Here, the same type using the FAP generic primer could be detected, indicating this set as an effective alternative for the identification of BPVs. Other than that, we described the simultaneous presence of two BPV types in three different wart samples, obtained from a chronically affected animal with disseminated papillomatosis.

The typed BPV-11 and BAPV-3 samples have a fibrous aspect, with a similar fibropapilloma histology that is not commonly associated with the exclusive infection of the keratinocytes, commonly attributed to the Xi BPVs. In a previous report [10], Hatama and others discussed the “uncertain nature of BPV-11 tumorigenicity” since BPV-11 was first diagnosed in a fibropapilloma lesion in which the BPV-1 was also detected. Despite the eventual limitations of the diagnostic methods employed, we could detect the same viral type in at least two different cauliflower-like lesions, with a suggestive fibrous core. In addition to Hatama’s report, our data also link the BPV-11 with cutaneous papillomas, indicating it as a pathogenic type. Furthermore, the putative type BAPV-3, described by Ogawa and others [16], was detected just once from a skin papilloma sample, without histological description, being genetically associated with the BPV-3, -4, and -6, or Xi BPVs.

The L1 gene has taxonomical relevance due to its high degree of conservation, which can be accessed with the use of generic primer sets. In accordance with this, the restriction sites located in FAP fragments appears to be maintained, without intratype variations for the restriction enzymes employed, indicating that phylogenetic studies comparing BPVs and other members of the Papillomaviridae family are possible using a relatively short DNA sequence.

Since the early nineties, methods based on PCR-RFLP analysis within the L1 gene of the human papillomavirus (HPV) have been used for viral typing and infection diagnosis from a variety of sources, including cervical samples, fresh and paraffin-embedded tissues [34, 35]. Particularly, the PCR-RFLP method is useful to identify coinfections due to its sensitivity and specificity [36].

5. Conclusion

As an easy, rapid, and cost-effective assay, the PCR-RFLP represents a less-laborious approach than DNA cloning and sequencing, being an alternative as a first-line screening test, both for the diagnosis of an already classified virus type as to indicating the needing of DNA sequencing due to mixed and/or unknown digestion profiles. In a previous survey in Paraná State, Brazil, Claus and others [26] suggest that the occurrence of multiple or mixed BPVs infection may be widespread throughout Brazilian cattle herds and may occur in other Brazilian geographical regions. In agreement, our findings support these results and reinforce the notion that multiple papillomaviruses infections, with a significant pathogenic potential, can be as frequent in cattle as in human hosts [26].

To the best of our knowledge, the restriction map employed here is the first to be produced specifically for the screening and typing of BPVs. Our findings also point to the ubiquity of BPVs dissemination since two relatively rare, new described BPV types, which were first characterized in Japan, were also detected in Brazil.

Acknowledgments

The authors thank the Ministério de Ciência, Tecnologia e Inovação/Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq Proc. 402539/2011-7), Coordenação de Aperfeiçoamento Pessoal de Nível Superior (CAPES), and Fundação do Desenvolvimento Administrativo (PAP-FUNDAP) for the financial support and Nadia Jorge Berriel, Carolina da Paz Sabino, and Jim Hesson (http://www.AcademicEnglishSolutions.com/) for the editorial support.

References

M. S. Campo, “Papillomavirus and disease in humans and animals,” Veterinary and Comparative Oncology, vol. 1, pp. 3–14, 2003.

S.-Y. Chan, H. Delius, A. L. Halpern, and H.-U. Bernard, “Analysis of genomic sequences of 95 papillomavirus types: uniting typing, phylogeny, and taxonomy,” Journal of Virology, vol. 69, no. 5, pp. 3074–3083, 1995.

E.-M. de Villiers, C. Fauquet, T. R. Broker, H.-U. Bernard, and H. Zur Hausen, “Classification of papillomaviruses,” Virology, vol. 324, no. 1, pp. 17–27, 2004.

G. Borzacchiello and F. Roperto, “Bovine papillomaviruses, papillomas and cancer in cattle,” Veterinary Research, vol. 39, no. 5, article 45, 2008.

M. S. Campo, W. F. H. Jarrett, R. Barron, B. W. O'Neil, and K. T. Smith, “Association of bovine papillomavirus type 2 and bracken fern with bladder cancer in cattle,” Cancer Research, vol. 52, no. 24, pp. 6898–6904, 1992.

F. Jelínek and R. Tachezy, “Cutaneous papillomatosis in cattle,” Journal of Comparative Pathology, vol. 132, no. 1, pp. 70–81, 2005.

L. Nasir and M. S. Campo, “Bovine papillomaviruses: their role in the aetiology of cutaneous tumours of bovids and equids,” Veterinary Dermatology, vol. 19, no. 5, pp. 243–254, 2008.

M. Lunardi, A. A. Alfieri, R. A. Otonel, et al., “Genetic characterization of a novel bovine papillomavirus member of the Deltapapillomavirus genus,” Veterinary Microbiology, vol. 162, no. 1, pp. 207–213, 2012.

W. F. H. Jarrett, M. S. Campo, B. W. O'Neil, H. M. Laird, and L. W. Coggins, “A novel bovine papillomavirus (BPV-6) causing true epithelial papillomas of the mammary gland skin: a member of a proposed new BPV subgroup,” Virology, vol. 136, no. 2, pp. 255–264, 1984.

S. Hatama, R. Ishihara, Y. Ueda, T. Kanno, and I. Uchida, “Detection of a novel bovine papillomavirus type 11 (BPV-11) using Xipapillomavirus consensus polymerase chain reaction primers,” Archives of Virology, vol. 156, no. 7, pp. 1281–1285, 2011.

W. Zhu, J. Dong, E. Shimizu et al., “Characterization of novel bovine papillomavirus type 12 (BPV-12) causing epithelial papilloma,” Archives of Virology, vol. 157, no. 1, pp. 85–91, 2012.

N. Bloch, R. H. Sutton, and P. B. Spradbrow, “Bovine cutaneous papillomas associated with bovine papillomavirus type 5,” Archives of Virology, vol. 138, no. 3-4, pp. 373–377, 1994.

Y. Tomita, I. Literak, T. Ogawa, Z. Jin, and H. Shirasawa, “Complete genomes and phylogenetic positions of bovine papillomavirus type 8 and a variant type from a European bison,” Virus Genes, vol. 35, no. 2, pp. 243–249, 2007.

T. Ogawa, Y. Tomita, M. Okada, and H. Shirasawa, “Complete genome and phylogenetic position of bovine papillomavirus type 7,” Journal of General Virology, vol. 88, no. 7, pp. 1934–1938, 2007.

A. Antonsson and B. G. Hansson, “Healthy skin of many animal species harbors papillomaviruses which are closely related to their human counterparts,” Journal of Virology, vol. 76, no. 24, pp. 12537–12542, 2002.

T. Ogawa, Y. Tomita, M. Okada et al., “Broad-spectrum detection of papillomaviruses in bovine teat papillomas and healthy teat skin,” Journal of General Virology, vol. 85, no. 8, pp. 2191–2197, 2004.

Y. Maeda, T. Shibahara, Y. Wada et al., “An outbreak of teat papillomatosis in cattle caused by bovine papilloma virus (BPV) type 6 and unclassified BPVs,” Veterinary Microbiology, vol. 121, no. 3-4, pp. 242–248, 2007.

M. P. Claus, M. Lunardi, A. F. Alfieri, L. M. Ferracin, M. H. P. Fungaro, and A. A. Alfieri, “Identification of unreported putative new bovine papillomavirus types in Brazilian cattle herds,” Veterinary Microbiology, vol. 132, no. 3-4, pp. 396–401, 2008.

S. Hatama, K. Nobumoto, and T. Kanno, “Genomic and phylogenetic analysis of two novel bovine papillomaviruses, BPV-9 and BPV-10,” Journal of General Virology, vol. 89, no. 1, pp. 158–163, 2008.

O. Forslund, A. Antonsson, P. Nordin, B. Stenquist, and B. G. Hansson, “A broad range of human papillomavirus types detected with a general PCR method suitable for analysis of cutaneous tumours and normal skin,” Journal of General Virology, vol. 80, no. 9, pp. 2437–2443, 1999.

R. C. Stocco Dos Santos, C. J. Lindsey, O. P. Ferraz et al., “Bovine papillomavirus transmission and chromosomal aberrations: an experimental model,” Journal of General Virology, vol. 79, no. 9, pp. 2127–2135, 1998.

S. R. Wosiacki, M. P. Claus, A. F. Alfieri, and A. A. Alfieri, “Bovine papillomavirus type 2 detection in the urinary bladder of cattle with chronic enzootic haematuria,” Memorias do Instituto Oswaldo Cruz, vol. 101, no. 6, pp. 635–638, 2006.

A. Yaguiu, C. Carvalho, A. C. Freitas, et al., “Papillomatosis in cattle: in situ detection of Bovine Papillomavirus DNA sequences in reproductive tissues,” Brazilian Journal of Morphological Sciences, vol. 23, no. 3, pp. 129–136, 2006.

N. Diniz, T. C. Melo, J. F. Santos et al., “Simultaneous presence of bovine papillomavirus in blood and in short-term lymphocyte cultures from dairy cattle in Pernambuco, Brazil,” Genetics and Molecular Research, vol. 8, no. 4, pp. 1474–1480, 2009.

M. Lunardi, M. P. Claus, A. A. Alfieri, and M. H. P. Fungaro, “Phylogenetic position of an uncharacterized Brazilian strain of bovine papillomavirus in the genus Xipapillomavirus based on sequencing of the L1 open reading frame,” Genetics and Molecular Biology, vol. 33, no. 4, pp. 745–749, 2010.

M. P. Claus, M. Lunardi, A. A. Alfieri, et al., “Multiple bovine papillomavirus infections associated with cutaneous papillomatosis in Brazilian cattle herds,” Brazilian Archives of Biology and Technology, vol. 52, pp. 93–98, 2009.

T. Vincze, J. Posfai, and R. J. Roberts, “NEBcutter: a program to cleave DNA with restriction enzymes,” Nucleic Acids Research, vol. 31, no. 13, pp. 3688–3691, 2003.

T. A. Hall, “BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT,” Nucleic Acid Symposium Series, vol. 41, pp. 95–98, 1999.

J. D. Thompson, D. G. Higgins, and T. J. Gibson, “CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice,” Nucleic Acids Research, vol. 22, no. 22, pp. 4673–4680, 1994.

S. F. Altschul, T. L. Madden, A. A. Schäffer et al., “Gapped BLAST and PSI-BLAST: a new generation of protein database search programs,” Nucleic Acids Research, vol. 25, no. 17, pp. 3389–3402, 1997.

K. Tamura, D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. Kumar, “MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,” Molecular Biology and Evolution, vol. 28, no. 10, pp. 2731–2739, 2011.

B. Efron, E. Halloran, and S. Holmes, “Bootstrap confidence levels for phylogenetic trees,” Proceedings of the National Academy of Sciences of the United States of America, vol. 93, no. 23, pp. 13429–13434, 1996.

R. D.m.page, “TreeView: an application to display phylogenetic trees on personal computers,” Computer Applications in the Biosciences, vol. 12, no. 4, pp. 357–358, 1996.

O. Lungu, T. C. Wright Jr., and S. Silverstein, “Typing of human papillomaviruses by polymerase chain reaction amplification with L1 consensus primers and RFLP analysis,” Molecular and Cellular Probes, vol. 6, no. 2, pp. 145–152, 1992.

G. Astori, A. Arzese, C. Pipan, E.-M. De Villiers, and G. A. Botta, “Characterization of a putative new HPV genomic sequence from a cervical lesion using L1 consensus primers and restriction fragment length polymorphism,” Virus Research, vol. 50, no. 1, pp. 57–63, 1997.

R. J. Nobre, L. P. de Almeida, and T. C. Martins, “Complete genotyping of mucosal human papillomavirus using a restriction fragment length polymorphism analysis and an original typing algorithm,” Journal of Clinical Virology, vol. 42, no. 1, pp. 13–21, 2008.