Am. J. Trop. Med. Hyg., 75(3), 2006, pp. 486-490
Copyright © 2006 by The American Society of Tropical Medicine and Hygiene
SHORT REPORT
ISOLATION OF TWO VACCINIA VIRUS STRAINS FROM A SINGLE BOVINE VACCINIA OUTBREAK IN RURAL AREA FROM BRAZIL: IMPLICATIONS ON THE EMERGENCE OF ZOONOTIC ORTHOPOXVIRUSES
GILIANE S. TRINDADE,
ZÉLIA I. P. LOBATO,
BETÂNIA P. DRUMOND,
JULIANA A. LEITE,
RICARDO C. TRIGUEIRO,
MARIA I. M. C. GUEDES,
FLÁVIO G. DA FONSECA,
JOÃO R. DOS SANTOS,
CLÁUDIO A. BONJARDIM,
PAULO C. P. FERREIRA, AND
ERNA G. KROON*
Laboratório de Vírus, Departamento de Microbiologia, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil; Departamento de Medicina Veterinária Preventiva, Escola de Veterinária, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil; Laboratório de Imunologia Celular e Molecular, Centro de Pesquisas René RachouFIOCRUZ, Belo Horizonte, Brazil
ABSTRACT
Outbreaks of bovine vaccinia disease caused by circulation of Vaccinia virus (VACV) strains have been a common occurrence in Brazil in the recent years, being an important emergent zoonosis. During a single outbreak that took place in 2001, two genetically different VACV strains were isolated and named Guarani P1 virus (GP1V) and Guarani P2 virus (GP2V). Molecular diagnosis was done through restriction fragment length polymorphism (RFLP) of ati gene (A26L) and by sequence analysis of a group of five VACV genes including the C11R, J2R, A56R, B18R, and E3L genes. These findings confirmed the co-circulation of two different Vaccinia virus strains during the same outbreak, raising important questions about the origin, emergence, and circulation of VACV strains in Brazil.
Since 1999, an increasing number of episodes of eruptive skin disease affecting dairy cattle and humans have been reported in Brazil.16 Our research team has been working on these zoonotic outbreaks, and in October of 2001, we studied an outbreak that occurred in Guarani town, Minas Gerais State, in the southeast region of the country.
An epidemiologic study was conducted on the affected area involving 72 visited properties, which comprised an area covered by ~5,700 km of roads. The study revealed that a total of 1,020 lactating cows presented lesions on the teats characterized by the presence of papules that evolved to ulcers. After some time, most of these lesions naturally progressed to healing. In 83% of the farms, human cases were registered, and ~110 persons were infected. Most sick humans were milkers that were contaminated after contact with lesions on cows teats. These milkers presented pleiomorphic lesions on the hands (mainly papules and painful ulcers), fever, lymphadenitis with enlarged lymph nodes, and eventually secondary bacterial infection on the lesions. Also, in some farms, milkers reported transmission from person to person. The time-course of the infection was between 15 and 30 days. Initial diagnosis was done based on the observed clinical and epidemiologic features.3
To attempt laboratory diagnosis and virus isolation, two samples of dried scabs from two affected cows were collected. Each cow belonged to a different neighboring farm, and these farms were localized 10 km apart from each other. To avoid any possibility of laboratory cross-contamination, the different samples were never manipulated simultaneously. Two viral isolates were obtained and named Guarani P1 virus (GP1V) and Guarani P2 virus (GP2V). Virus isolation was performed using conventional methods including inoculation onto chorioallantoic membrane (CAM) of chicken embryonated eggs and plaque purification cloning on Vero cells to assure genetic homogeneity. These viruses were grown and purified in sucrose gradient as described elsewhere.2,3
Biologic diagnosis was based on pock morphology in CAM and a serum neutralization test using anti-VACV polyclonal rabbit serum as a positive control.2,3 In addition, atomic force microscopy was also used to look for typical poxviruses particles on the samples (data not shown).7
Molecular diagnosis of the isolated viruses was performed by polymerase chain reaction (PCR) amplification of the A type inclusion body gene (A26L). The A26L gene analysis, based on partial amplification of this gene followed by digestion using XbaI, has been widely used as a tool for rapid screening and taxonomic differentiation of Orthopoxvirus.8,9 The GP2V A26L gene was amplified, and the amplicon was similar in size to those described for Araçatuba virus (ARAV) and Passatempo virus (PSTV), other VACV strains isolated under similar circumstances in 1999 and 2003, respectively.2,6 When the A26L amplicon of GP2V was digested with XbaI, it presented the same restriction pattern as in ARAV, which was also similar but not identical to those patterns obtained for known VACV strains, such as Western Reserve (WR) and Lister (LST) (Figure 1A
). The main difference observed for the GP2V and ARAV A26L gene amplicon digests was that the larger fragment generated after XbaI digestion migrated faster on the agarose gel than the VACV-WR and VACV-LST analogous fragments (Figure 1A
). Attempts to amplify the A26L gene of GP1V, using the ATI up and ATI low primers, were unsuccessful (data not shown).8 The lack of this amplicon could be caused by a deletion at the A26L gene of the GP1V, a fact previously observed for other Brazilian VACV strains, such as BeAn58058 virus (BAV) and Belo Horizonte Virus (VBH).10,11 To verify the extension of the probable GP1V A26L deletion, we performed a PCR amplification using oligonucleotides P4c1 (located within the p4c1 gene) and RNApol, designed based on Cowpox virus Brighton Red (CPXV BR) as described elsewhere.11 For CPXV, VACV-WR, and GP2V, a product of ~4.300 bp was obtained. However, for GP1V, a fragment of 300 bp (Figure 1B
) was generated and sequenced (GenBank accession no. DQ363383), and it proved that, in GP1V, as well in VBH and BAV, a major portion of the A26L gene is missing (Figure 1C
). Additionally, it is important to point out the usefulness of the A26L PCR-restriction fragment length polymorphism (PCR-RFLP) as a molecular diagnosis approach and for differentiation of Brazilian Orthopoxvirus species.10,11

View larger version (43K):
[in this window]
[in a new window]
|
FIGURE 1. (A), PCR-RFLP analysis of A26L gene amplicons of VACVs isolated in Brazil. Lanes were grouped: odd numbers correspond to non-digested PCR and even numbers to RFLP (M, molecular size marker): 12, VACV-WR; 34, VACV-LST; 56, GP2V; 78, ARAV; 910, CPXV-BR. PCR amplicons digestions were performed using XbaI endonuclease. PCR and PCR-PFLP fragments were separated by electrophoresis in a 1% agarose gel stained with EtBr. (B), PCR amplification employing P4c1 and RNApol primers. Products were fractionated by electrophoresis in a 1% agarose gel stained with EtBr. Lanes (M, molecular size marker; 1-Kb DNA ladder; Promega): 1, negative control; 2, GP1V; 3, GP2V; 4, VACV-WR; 5, CPXV-BR. (C), Schematic comparison of sequenced regions of VAC-WR, GP2V, GP1V, VBH, and CPXV-BR based on work by Meyer and others.8 Boxes, genes; thin arrows, primer positions; bold arrows, orientation of genes; dashed lines, deletion.
|
|
Because GP1V and GP2V presented extensive differences in their A26L ORF, we decided to study other genes to evaluate genomic differences between GP1V and GP2V and other Brazilian VACV strains, such as ARAV, PSTV, Cantagalo virus (CTGV), BAV, and VBH.1,2,6,1113 The thymidine kinase (J2R), vaccinia growth factor (C11R), hemagglutinin (A56R), soluble alpha/beta interferon (IFN) receptor (B18R), and dsRNA-binding protein of VACV (E3L) were chosen for sequencing and analysis. J2R, C11R, and A56R are conserved genes in the Orthopoxvirus genus and have been used before as a taxonomic tool to study ungrouped poxvirus.11,12,14 B18R and E3L genes, both involved in the viral mechanism of pathogenesis, were also amplified from GP1V, GP2V, and other Brazilian VACV strains.15,16 Primers and the reactions conditions to amplify the B18R gene were those described by Marques and others.10 The primers based on the E3L nucleotide sequence of VACV-WR were E3L sense (5'-AGGCGGATCCATGTCTAAAATCTATATC-3') and E3L anti-sense (5'-TCGCAAGCTTTCAGAATCTAATG-ATGAC-3'), containing restriction sites for BamHI and HindIII (underlined), respectively, and they were used to amplify a 589-bp DNA fragment. Annealing was carried out at 45°C, and PCR reactions were performed as described by Marques and others.10
Amplicons were cloned into pGEM-T vector (Promega Corp., Madison, WI) and sequenced in both orientations by the dideoxynucleoside method, using M13 universal primers and ET Dynamic Terminator for MegaBACE (GE HEALTHCARE, UK).17 Nucleotide sequences were assembled using the CAP3 Sequence Assembling Program and deposited in GenBank.18 Nucleotide and inferred amino acids sequences were aligned with other Orthopoxvirus sequences using CLUSTAL W and LALIGN.19,20 J2R sequences of GP1V and GP2V (GenBank accession no. DQ206438 and DQ206439, respectively) were identical to each other, presenting high levels of similarity with VACV sequences. These sequences were also identical to VACV-WR and other Brazilian VACV strains, such as BAV, SPAn232, ARAV, and PSTV (data not shown).2,6,12,13 C11R sequences of GP1V and GP2V (GenBank accession no. DQ206440 and DQ206441, respectively) presented 99.4% of identity compared with each other and high levels of similarity among other VACV strains.
B18R, E3L, and A56R sequences of GP1V and GP2V were also similar to other VACV strains sequences, although there was a higher level of genetic variance among them. GP1V and GP2V B18R sequences presented seven nucleotide substitutions that resulted in the two amino acid modifications in the inferred amino acid sequences. For GP1V, GP2V, ARAV, and BAV, in the position 346 of the B18R protein, a serine was found, whereas in VACV-WR, VACV-IOC, and CTGV, there was a threonine in the same position (data not shown). The GP1V and GP2V E3L genes presented 98.6% of identity to each other, with substitutions in eight nucleotide positions (Table 1
). These mutations led to six amino acid differences observed between the GP1V and GP2V inferred amino acid sequences (data not shown).
However, when the GP1V and GP2V A56R sequences were compared, relatively lower levels of similarity were observed in both nucleotide (95.7%) and amino acid (92.9%) levels. Comparing the GP2V A56R sequence to other Brazilian VACVs, similarity values ranging from 99.1% to 99.3% were obtained (Table 1
). As observed for ARAV, PSTV, CTGV, and VACV-IOC, GP2V also presented the same 18-nucleotide deletion signature showed previously at the A56R gene of those VACV strains. This 18-nucleotide deletion has been found in all Brazilian VACV strains isolated to date, and it has been proposed as a molecular signature of new world VACV strains. However, GP1V and VBH did not present this deletion. Additionally, the lack of this deletion in the GP1V and VBH sequences points out the necessity to review the 18-nucleotide deletion as a molecular signature of new world VACV strains.
To establish a phylogenetic relationship among GP1V, GP2V, and other VACV strains, alignments of the A56R, B18R, and E3L genes were used to construct a concatenated phylogenetic tree by the neighbor-joining method using the Tamura Nei model of nucleotide substitutions implemented in MEGA3.21,22 The tree was midpoint-rooted, and 1,000 bootstrap replicates were performed (Figure 2A
).

View larger version (24K):
[in this window]
[in a new window]
|
FIGURE 2. A, Consensus concatenated phylogenetic tree based on the nucleotide sequences of Orthopoxvirus A56R, B18R, and E3L genes. The tree was constructed by the neighbor-joining method using the Tamura Nei model of nucleotide substitutions implemented in MEGA3. The tree was midpoint-rooted, 1,000 bootstrap replicates were performed, and values > 70% are shown. Nucleotide sequences were obtained from GenBank. The GenBank accession numbers of A56R, B18R, and E3L sequences are GP1V (DQ206436, DQ194380, and DQ194385), GP2V (DQ206437, DQ194381, and DQ194386), ARAV (AY523994, DQ194382, and DQ194389), PSTV (DQ070848, DQ530239, and DQ530240), and VBH (DQ206435, DQ194383, and DQ194390), respectively. Brazilian VACV strains are shown in bold. , Brazilian strains that presented partial A26L gene deletion; , Brazilian strains that presented 18 nt deletion in the A56R gene. B, Brazilian states where and when Brazilian VACV strains were isolated. ES, Espírito Santo; GO, Goiás; MG, Minas Gerais; PA, Pará; RJ, Rio de Janeiro; SP, São Paulo. 1, BeAn 58058 virus in 1963; 2, SPAn232 virus in 1965; 3, Belo Horizonte virus in 1993; 4, ARAV in 1999; 5, CTGV in 1999; 6, GP1V and GP2V in 2001; 7, GO isolates after 2001; 8, SP isolates after 2001; 9, MG isolates after 2001; 10, PSTV in 2003; ES isolates in 2004 (unpublished data).
|
|
As expected, the concatenated phenogram tree grouped GP1V and GP2V within VACV strains. Although isolated from the same outbreak, at the same time, and from neighboring farms, GP1V and GP2V presented enough genetic divergence to place them in different branches of the phylogenetic tree. Whereas GP2V was clustered with Brazilian VACV strains isolated in previous bovine vaccinia outbreaks (ARAV, CTGV, and PSTV), GP1V was grouped together with VACV-WR and VBH. Neither GP1V nor GP2V were clustered with VACV strains LST and IOC that were used in Brazil during the smallpox vaccination campaign.1,2
Our findings show that there are genetically different populations of VACV circulating in the country and even in the same outbreak. It has been proposed that the Brazilian VACV isolates obtained from cowpox-like outbreaks form a genetically homogenous population that could have been originated from the spread of a smallpox vaccine strains, particularly the VACV IOC.1 According to this theory, these strains could have escaped to wild and established circulation in some unknown host, becoming feral during this process.1,2,6 However, the characterization of two different VACV strains isolated from a single outbreak, co-circulating in Brazil, points in a different direction. One hypothesis to explain this virus diversity is that different samples, from different sources, and under different circumstances, could have established circulation in nature. Thus, these viruses would have not one but multiple origins. This does not entirely exclude the possibility that the viruses could have evolved from a common origin, derived from a escapee vaccine strain, and differentiated by means of mutation and/or recombination events during circulation in nature. However, our results suggest that multi-factor events may be the best explanation for the VACV genetic diversity observed in Brazil, including the pre-existence and circulation of a possible autochthonous poxvirus.
Received February 15, 2006.
Accepted for publication May 10, 2006.
Acknowledgments: The authors thank Dr. Marieta Cristina Madureira from Instituto Mineiro de Agropecuária for support during the epidemiologic investigation. We also thank Dr. José Mário Vilela and Dr. Margareth Spangler Andrade from Nanoscopy Laboratory, Fundação Centro Tecnológico de Minas Gerais, for the AFM analysis. The aid of Dr. Fabrício dos Santos, Rodrigo Redondo, and Laboratório de Biodiversidade e Evolução Molecular where all sequences were made is gratefully acknowledged.
Financial support: The CNPq, CAPES, and FAPEMIG provided financial support. E. G. Kroon, C. A. Bonjardim, P. C. P. Ferreira, F. G. Da Fonseca, and J. R. Santos received fellowship from CNPq. M. I. M. C. Guedes received fellowship from FAPEMIG.
* Address correspondence to Erna Geessien Kroon, Laboratório de Vírus, Departamento de Microbiologia, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Av. Antônio Carlos, 6627, caixa postal 486, CEP 31270-901 Belo Horizonte, MG, Brazil. E-mail: kroone{at}icb.ufmg.br 
Authors addresses: Giliane de Souza Trindade, Betânia Paiva Drumond, Juliana Almeida Leite, Ricardo Campos Trigueiro, Maria Isabel Maldonado Coelho Guedes, João Rodrigues dos Santos, Cláudio Antônio Bonjardim, Paulo César Peregrino Ferreira, and Erna Geessien Kroon, Laboratório de Vírus, Departamento de Microbiologia, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Av. Antônio Carlos, 6627, caixa postal 486, CEP 31270-901, Belo Horizonte, MG, Brasil, Fax: 55 31 3443-6482, E-mails: gitrindade{at}yahoo.com.br, betaniadrumond{at}uol.com.br, ju{at}icb.ufmg.br, rictrigueiro{at}yahoo.com.br, isabelguedes{at}icb.ufmg.br, santosjr{at}icb.ufmg.br, claubonj{at}icb.ufmg.br, paulocpf{at}icb.ufmg.br, and kroone{at}icb.ufmg.br. Zélia Inês Portela Lobato, Departamento de Medicina Preventiva, Escola de Veterinária, Universidade Federal de Minas Gerais, Av. Antônio Carlos, 6627, CEP 31270-901, Belo Horizonte, MG, Brasil, Fax: 55 31 3443-6482, E-mail: ziplobat{at}vet.ufmg.br. Flávio Guimarães da Fonseca, Departamento de Microbiologia, Instituto de Ciências Biológicas, Universidade Federal de Minas Gerais, Av. Antônio Carlos, 6627, caixa postal 486, CEP 31270-901, Belo Horizonte, MG, Brasil, Fax: 55 31 3443-6482, E-mail: fdafonseca{at}cpqrr.fiocruz.br.
REFERENCES
- Damaso CRA, Esposito JJ, Condit R, Moussatché N, 2000. An emergent poxvirus from humans and cattle in Rio de Janeiro state: Cantagalo virus may derive from Brazilian smallpox vaccine. Virology 277: 439449.[Web of Science][Medline]
- Leite JA, Drumond BP, Trindade GS, Lobato ZIP, da Fonseca FG, dos Santos JR, Madureira MC, Guedes MIMC, Ferreira JMS, Bonjardim CA, Ferreira PC, Kroon EG, 2005. Passatempo virus: a novel Vaccinia virus isolated during a zoonotic outbreak in Brazil. Emerg Infect Dis 11: 19351938.[Web of Science][Medline]
- Lobato ZIP, Trindade GS, Frois MCM, Ribeiro EBT, Dias GRC, Teixeira BM, Lima FA, Almeida GMF, Kroon EG, 2005. Outbreak of exantemal disease caused by Vaccinia virus in human and cattle in Zona da Mata Region, Minas Gerais. Arq Bras Med Vet Zootech 57: 423429.
- Nagasse-Sugahara TK, Kisielius JJ, Ueda-Ito M, & 8 other authors 2004. Human vaccinia-like virus outbreaks in Sao Paulo and Goias States, Brazil: Virus detection, isolation and identification. Rev Inst Med Trop Sao Paulo 46: 315322.[Medline]
- Schatzmayr HG, Lemos ER, Mazur C, Schubach A, Majerowicz S, Rozental T, Schubach TM, Bustamente MC, Barth OM, 2000. Detection of poxvirus in cattle associated with human cases in the State of Rio de Janeiro: Preliminary report. Mem Inst Oswaldo Cruz 95: 625627.[Web of Science][Medline]
- Trindade GS, Fonseca FG, Marques JT, Nogueira ML, Mendes LC, Borges AS, Pituco EM, Bonjardim CA, Ferreira PC, Kroon EG, 2003. Araçatuba virus: A vaccinia-like virus associated with cattle and human infection. Emerg Infect Dis 9: 155160.[Web of Science][Medline]
- Malkin AJ, McPherson A, Gershon PD, 2003. Structure of intra-cellular mature vaccinia virus vizsualized by in situ atomic force microscopy. J Virol 77: 63326340.[Abstract/Free Full Text]
- Meyer H, Roop SL, Esposito JJ, 1997. Gene for A-type inclusion body protein is useful for a polymerase chain reaction assay to differentiate orthopoxvirus. J Virol Meth 64: 217221.[Web of Science][Medline]
- Neubauer H, Pfeffer M, Meyer H, 1997. Specific detection of mousepox virus by polymerase chain reaction. Lab Anim Sci 31: 201205.
- Marques JT, Trindade GS, da Fonseca FG, dos Santos JR, Bonjardim CA, Ferreira PCP, Kroon EG, 2001. Characterization of ATI, TK and IFN-alpha/betaR genes in the genome of the BeAn 58058 virus, a naturally attenuated wild Orthopoxvirus. Virus Genes 23: 291301.[Web of Science][Medline]
- Trindade GS, da Fonseca FG, Marques JT, Diniz S, Leite JA, De Bodt S, Van der Peer Y, Bonjardim CA, Ferreira PCP, Kroon EG, 2004. Belo Horizonte virus: A vaccinia-like virus lacking the A-type inclusion body gene isolated from infected mice. J Gen Virol 85: 20152021.[Abstract/Free Full Text]
- da Fonseca FG, Trindade GS, Silva RLA, Bonjardim CA, Ferreira PCP, Kroon EG, 2002. Characterization of a vaccinia-like virus isolated in a Brazilian forest. J Gen Virol 83: 223228.[Abstract/Free Full Text]
- Fonseca FG, Lanna MCS, Campos MAS, Kitajima EW, Peres JN, Golgher RR, Ferreira PCP, Kroon EG, 1998. Morphological and molecular characterization of the poxvirus BeAn 58058. Arch Virol 143: 11711186.[Web of Science][Medline]
- Roop SL, Jin QI, Knight JC, Massung RF, Esposito JJ, 1995. PCR strategy for identification and differentiation of smallpox and other Orthopoxviruses. J Clin Microbiol 33: 20692076.[Abstract/Free Full Text]
- Goebel SJ, Johson GP, Perkus ME, Davis SW, Winslow JP, Paoletti E, 1990. The complete DNA sequence of vaccinia virus. Virology 179: 247266.[Web of Science][Medline]
- Johnson GP, Goebel SJ, Paoletti E, 1993. An update on the Vaccinia virus genome. Virology 381: 381401.
- Sanger F, Nicklen S, Coulson AR, 1977. DNA sequencing with chain-terminating inhibitors. Proc Natl Acad Sci USA 74: 54635467.[Abstract/Free Full Text]
- Huang X, Madan A, 1999. CAP3: A DNA sequence assembly program. Genome Res 9: 868877.[Abstract/Free Full Text]
- Pearson W. LALIGN-find multiple matching subsegments in two sequences. http://www.ch.embnet.org/software/LALIGN_form.html. Accessed January 2006.
- Thompson JD, Higgins DG, Gibson TJ, 1994. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 46734680.[Abstract/Free Full Text]
- Tamura K, Nei M, 1993. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 10: 512526.[Abstract]
- Kumar S, Tamura K, Nei M, 2004. MEGA3: Integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief Bioinform 5: 150163.[Abstract/Free Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
J. M. S. Ferreira, J. S. Abrahao, B. P. Drumond, F. M. Oliveira, P. A. Alves, M. A. Pascoal-Xavier, Z. I. P. Lobato, C. A. Bonjardim, P. C. P. Ferreira, and E. G. Kroon
Vaccinia virus: shedding and horizontal transmission in a murine model
J. Gen. Virol.,
December 1, 2008;
89(12):
2986 - 2991.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
G. de Souza Trindade, B. P. Drumond, M. I. M. C. Guedes, J. A. Leite, B. E. F. Mota, M. A. Campos, F. G. da Fonseca, M. L. Nogueira, Z. I. P. Lobato, C. A. Bonjardim, et al.
Zoonotic Vaccinia Virus Infection in Brazil: Clinical Description and Implications for Health Professionals
J. Clin. Microbiol.,
April 1, 2007;
45(4):
1370 - 1372.
[Abstract]
[Full Text]
[PDF]
|
 |
|