Next Article in Journal
Human Endogenous Retrovirus Group E and Its Involvement in Diseases
Next Article in Special Issue
Molecular Genetic Analysis of Orf Virus: A Poxvirus That Has Adapted to Skin
Previous Article in Journal
Simple Mathematical Models Do Not Accurately Predict Early SIV Dynamics
Previous Article in Special Issue
The Origin of the Variola Virus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

From Lesions to Viral Clones: Biological and Molecular Diversity amongst Autochthonous Brazilian Vaccinia Virus

1
Laboratório de Vírus, Departamento de Microbiologia, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais 31270-901, Brazil
2
AQUACEN) Laboratório Nacional Oficial de Referência de Doenças de Animais Aquáticos, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais 31270-901, Brazil
3
Laboratório de Retroviroses, Departamento de Medicina Veterinária Preventiva, Escola de Veterinária da Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais 31270-901, Brazil
4
Centro Agropecuário, Departamento de Ciência Animal, Universidade Federal do Pará, Pará 66075-110, Brazil
*
Author to whom correspondence should be addressed.
Submission received: 17 December 2014 / Revised: 5 March 2015 / Accepted: 9 March 2015 / Published: 16 March 2015
(This article belongs to the Special Issue Poxvirus Evolution)

Abstract

:
Vaccinia virus (VACV) has had an important role for humanity because of its use during the smallpox eradication campaign. VACV is the etiologic agent of the bovine vaccinia (BV), an emerging zoonosis that has been associated with economic, social, veterinary and public health problems, mainly in Brazil and India. Despite the current and historical VACV importance, there is little information about its circulation, prevalence, origins and maintenance in the environment, natural reservoirs and diversity. Brazilian VACV (VACV-BR) are grouped into at least two groups based on genetic and biological diversity: group 1 (G1) and group 2 (G2). In this study, we went to the field and investigated VACV clonal diversity directly from exanthemous lesions, during BV outbreaks. Our results demonstrate that the G1 VACV-BR were more frequently isolated. Furthermore, we were able to co-detect the two variants (G1 and G2) in the same sample. Molecular and biological analysis corroborated previous reports and confirmed the co-circulation of two VACV-BR lineages. The detected G2 clones presented exclusive genetic and biological markers, distinct to reference isolates, including VACV-Western Reserve. Two clones presented a mosaic profile, with both G1 and G2 features based on the molecular analysis of A56R, A26L and C23L genes. Indeed, some SNPs and INDELs in A56R nucleotide sequences were observed among clones of the same virus population, maybe as a result of an increased mutation rate in a mixed population. These results provide information about the diversity profile in VACV populations, highlighting its importance to VACV evolution and maintenance in the environment.

1. Introduction

Vaccinia virus (VACV) has had an important role in human history due to its use as a vaccine during the smallpox vaccination campaign, resulting in the eradication of this deadly disease in 1980 [1]. More recently, VACV has been widely used as a vector for recombinant vaccines [2,3,4]. Following Smallpox eradication, other zoonotic orthopoxviruses (OPV) have emerged worldwide, such as Cowpox virus (CPXV), in Europe, Monkeypox virus (MPXV), endemic in many African countries and recently introduced in the USA, and finally VACV, endemic in Brazil and India [5,6,7,8,9,10,11,12,13]. VACV is the causative agent of bovine vaccinia (BV), an exanthemous disease responsible for outbreaks that affect both bovines and humans, causing public health impacts and economic losses in South America and India [10,11,12,14,15]. The clinical prognosis of BV includes the appearance of papules that progress to vesicles and growing scabs mainly on the teats and udder of infected bovines. In humans, the lesions occur primarily on milkers’ hands and arms. Other symptoms, such as fever, myalgia, headache, arthralgia and lymphadenopathy are also frequently reported [11,16]. During outbreaks, direct contact with infected animals is the main transmission route to humans, but it may also occur by contact with contaminated fomites. Furthermore, reports have suggested a possible human-to-human transmission of VACV [16,17,18].
In Brazil, VACV was first isolated from rodents in the 1960s and since 1999 BV outbreaks have been consistently reported [10,12,15,19,20,21,22,23,24]. Some reports have shown a great genetic and biologic heterogeneity among Brazilian VACV (VACV-BR) isolates. This variability allowed clustering of VACV-BR into at least two distinct groups (Group 1—G1 and Group 2—G2), demonstrating a dichotomy supported by evidence of molecular and biological diversity such as virulence in BALB/c mouse model and plaque phenotype in BSC-40 cells. Group 2 strains display higher plaque sizes and are virulent to mice, unlike Group 1 [24,25,26,27,28]. Moreover, polymorphisms observed in specific VACV genes, such as hemagglutinin gene (A56R), A-type inclusion body gene (A26L), and chemokine binding protein gene (C23L), have been used in phylogenetic studies and further confirmed the dichotomy between G1 and G2 VACV-BR [22,27,29,30].
The VACV-BR G2 comprises of viral isolates Guarani P1 virus (GP1V), Pelotas 1 virus (P1V), Belo Horizonte virus (VBH), the newest isolate Serro human virus 2 (SH2V) and others, while G1, the group most frequently reported during BV outbreaks, includes isolates of Araçatuba virus (ARAV), Cantagalo virus (CTGV), Mariana virus (MARV), Guarani P2 virus (GP2V), Pelotas 2 virus (P2V) and others [22,23,24,26,28,31,32]. It is worth mentioning that the GP1V and GP2V were isolated during the same outbreak (co-circulation), while P1V and P2V were isolated from the same clinical sample (co-infection) [26,28].
Several studies have demonstrated genetic and biological diversity among isolates of the same species of a given OPV, including Variola virus (VARV), CPXV and MPXV. Diversity studies among MPXV isolates have revealed the existence of a dichotomy supported by differences in epidemiological and clinical features, while CPXV is clustered in four or more groups. Furthermore, studies on VARV variability have demonstrated the existence of at least two groups with distinct virulence profiles: variola major and variola minor. Taken together, these studies showed a great biological and genetic diversity amongst isolates from different OPV species [33,34,35,36,37,38,39,40].
Osborne and colleagues in 2007 reported high genetic variability among clones of the smallpox vaccine Dryvax [41]. Moreover, other studies reveal that old vaccines were a complex mixture of viruses with a large number of single-nucleotide polymorphisms (SNPs) besides insertions and deletions (INDELs), resulting in genetically and phenotypically distinct clones [42,43,44,45]. Besides, recombination events can occur among OPVs since it has been observed in VACV co-cultures and between VARV and CPXV [46,47,48]. Together, these studies suggest that poxviruses from the same species exist as distinct genetical groups and can recombine inside co-infected hosts, which increases the importance of studying the genetic variability among strains and clones of poxviruses circulating in the environment. Even though studies using vaccine strains clones show that VACV can recombine, there is an absence of data concerning VACV clones from field samples. This work aims to contribute to the understanding about the biological and genetic variability among field VACV clones by using genetic and biological analyses to study the diversity of VACV populations isolated during BV outbreaks.

2. Materials and Methods

2.1. Ethics Approval

The study was approved by the Committee of Ethics in Animal Use from the Universidade Federal de Minas Gerais (CEUA/UFMG). The animals were anesthetized by intraperitoneal injection of ketamine (Vetnil, São Paulo, Brazil) and xylazine (União Química, São Paulo, Brazil) before the infection and those ones that had lost more than 25% of their initial body weight were euthanized with an overdose of anesthetics.

2.2. Cells and Viruses

African green monkey kidney BSC-40 and VERO cells (American type cell culture—ATCC) were maintained in 5% CO2 atmosphere at 37 °C in Eagle’s Minimum Essential Medium (MEM) (Gibco BRL, Invitrogen, Carlsbad, CA, USA) supplemented with 5% fetal bovine serum (FBS) (Cultilab, Brazil), 25 µg/mL fungizone (Amphotericin B) (Cristália, São Paulo, Brazil), 500 U/mL penicillin (Cristália, São Paulo, São Paulo, Brazil) and 50 µg/mL gentamicin (Schering-Plough, São Paulo, Brazil). VERO cells were used for viral isolation from clinical specimens and growth. The BSC-40 cells were used for viral clones purification, plaque phenotype, comet phenotype and growth curve assays. Only typical poxviruses cytopathic effects were visualized after inoculation of clinical samples in both cell systems. The VACV Western Reserve (VACV-WR) was gently provided by Dr C. Jungwirth (Universität Würzburg, Würzburg, Germany) and was used as a virulent control in mice assays.

2.3. Clinical Samples

Clinical samples, consisting of swabs of vesicular fluids and scabs obtained from cattle’s teats and milkers’ hands, were collected during BV outbreaks that occurred between 2005 and 2011 in five Brazilian states: Minas Gerais, Bahia, Goiás, Espírito Santo and Pará. Sample collection procedures were performed by veterinarians, following institutional recommendations. Two viruses were isolated from human samples, obtained in Bahia and Minas Gerais states; and eight viruses were isolated from bovine samples, two obtained from Bahia state, two from Para state, two from Goiás state, one from Minas Gerais state and one from Espírito Santo state. The counties and states where the samples were isolated are shown in Table 1 and Figure 1.
Table 1. Identification of viral clones isolated from clinical samples.
Table 1. Identification of viral clones isolated from clinical samples.
VACV-BR isolates (abbreviation)Isolated clonesSourceIsolation stateIsolation year
VACV-BABVVACV-BABVBos taurusBahia2011
Clone 1 to Clone 5
VACV-BAB2VVACV-BAB2VBos taurusBahia2011
Clone 1 to Clone 5
VACV-BAHVVACV- BAHVHomo sapiensBahia2011
Clone 1 to Clone 5
VACV-GOBVVACV-GOBVBos taurusGoiás2011
Clone 1 to Clone 5
VACV-GOB2VVACV-GOB2VBos taurusGoiás2011
Clone 1 to Clone 5
VACV-ESBVVACV-ESBVBos taurusEspírito Santo2008
Clone 1 to Clone 5
VACV-MGBVVACV-MGBVBos taurusMinas Gerais2005
Clone 1 to Clone 4
VACV-MGHVVACV-MGHVHomo sapiensMinas Gerais2011
Clone 1 to Clone 5
VACV-PABVVACV-PABVBos taurusPará2011
Clone 1 to Clone 5
VACV-PAB2VVACV-PAB2VBos taurusPará2010
Clone 1 to Clone 4
Figure 1. Map displaying Brazilian states and counties from which the clinical samples used for Vaccinia virus clone isolation were collected. The red circles indicate the counties where the samples were collected: Miguel Calmon, Baixa Grande and Mundo Novo counties localized in Bahia state, Vila Propício county localized in Goiás state, Bom Jesus do Tocantins county localized in Pará state, Resplendor and Argirita counties localized in Minas Gerais States and Alegre county in Espírito Santo state.
Figure 1. Map displaying Brazilian states and counties from which the clinical samples used for Vaccinia virus clone isolation were collected. The red circles indicate the counties where the samples were collected: Miguel Calmon, Baixa Grande and Mundo Novo counties localized in Bahia state, Vila Propício county localized in Goiás state, Bom Jesus do Tocantins county localized in Pará state, Resplendor and Argirita counties localized in Minas Gerais States and Alegre county in Espírito Santo state.
Viruses 07 01218 g001

2.4. Sample Processing and Virus Isolation

Scabs were macerated using a homogenizer (Polytron, Lucerne, Switzerland) in phosphate buffer saline (PBS) (0.1 g scab/0.9 mL PBS) containing 200 U/mL penicillin, 4 µg/mL amphotericin B and 100 µg/mL gentamicin, and then clarified by centrifugation at 2000× g for 3 min (min). To the swabs, 200 µL of PBS was added, followed by a centrifugation at 2000× g for 3 min [32]. For virus isolation, 200 µL of supernatant was added to VERO cells monolayers which were incubated after adsorption at 37 °C for 72 h or until detection of cytopathic effect (CPE).

2.5. Viral Clone Plaque Purification

BSC-40 cells were used for viral plaque purification as described elsewhere [49]. After isolation, each virus isolate was serially diluted (log10) and inoculated into BSC-40 cells monolayers cultured in 6-well plates for clone selection and viral plaque purification. After 1 h of incubation (at 37 °C, 5%CO2), culture medium was removed, cell monolayers were washed twice with PBS and solid medium (MEM supplemented with 1% FBS + agarose 1%—Gibco, São Paulo, Brazil) was added. After 48 hours, five viral plaques from each sample were collected and purified in BSC40 cells with solid medium. The criteria used to choose the clones was the plaque size: after measurement, small and large viral plaques were collected from each sample. Individual clones were then isolated by serial plaque-purifications in three passages in BSC40 cells using 1% agarose overlay. The resulting clones were propagated in Vero cells monolayers, purified and titrated as described [49,50]. Isolated VACV strains and clones are shown in Table 1. Two clones, from two distinct clinical samples, were not able to grow after some steps of clone plaque purification (Table 1). After clone plaque purification, clones were manipulated separately to avoid cross-contamination.

2.6. Biological Assays

2.6.1. Plaque Phenotype

For plaque phenotype assays, BSC40 cells seeded in 6-well plates at 90%–95% confluence were infected with specific clones. After 1 h of adsorption (37 °C, 5%CO2), monolayers were washed twice with PBS and overlaid with solid medium prepared by mixing equal proportion of 1% agarose and 2× Eagle’s minimum essential medium (MEM) (Gibco, São Paulo, Brazil) supplemented with 2% FBS (Gibco, São Paulo, Brazil). After 48 h of incubation (37 °C, 5% CO2), cells were fixed with formaldehyde and stained with crystal violet for plaque size analysis.

2.6.2. Comet Phenotype Assay

For comet phenotype assays, BSC40 cells seeded in 6-well plates at 90%–95% confluence were infected with the large plaque clones. After 1 h of adsorption (37 °C, 5% CO2), the medium was removed, monolayers were washed twice with PBS and fresh MEM supplemented with 1% FBS was added. After 48 h of incubation (37 °C, 5% CO2), the cells were fixed with paraformaldehyde and stained with crystal violet.

2.6.3. Growth Curve Assays

BSC-40 cells were grown in 24-well plates and infected with VACV clones in biological duplicates. Infection was carried out at a MOI (multiplicity of infection) of 0.01 for 0, 12, 24 and 48 h post infection (hpi). Infected monolayers were then harvested, frozen (−70 °C) and thawed (37 °C) three times and titrated in BSC-40 cells as described. The analysis was performed with GraphPad Prism software (Version 5.00, GraphPad Software, San Diego, CA, USA, 2012).

2.7. Virulence in BALB/c Mice

To compare the virulence profile of these new clonal VACV isolates, five clones were selected for virulence assays. Mice were housed in filter-top microisolator cages and provided with commercial mouse food and water ad libitum. All the animal experiments were carried out in accordance with regulations and guidelines of the Committee of Ethics in Animal Use of the Universidade Federal de Minas Gerais/Brazil. BALB/c mice were anesthetized by intraperitoneal injection of ketamine and xylazine (3.2 mg and 0.16 mg/mice in 0.9% PBS, respectively) before the procedure (dose: 80–120 mg/kg per animal). Groups of four-week old male BALB/c mice (n = 4) were inoculated intranasally with 10 µL of viral suspensions containing 106 plaque forming units (pfu) VACV-WR was used as positive control, given its virulence for mice, and the negative control group was inoculated with 10 µL of PBS as described in [25]. Mice were weighed daily, and clinical signs were recorded for 10 days post infection (dpi) [25]. Mice that lost more than 25% of their initial weight were euthanized.

2.8. Molecular Assays

Amplification of A56R, A26L and C23L Genes and Phylogenetic Analyses

For the clones’ molecular characterization, DNA extractions were carried out using phenol-chloroform-isoamyl alcohol (PCI) [51] and used as template for Polymerase chain reaction (PCR) amplification of the viral genes A56R (hemagglutinin), C23L (chemokine binding protein) and A26L (A-type inclusion body). These genes are traditionally used for OPV phylogenetic studies, including Brazilian VACV strains [27,29,30].
Two different reactions were made for C23L: reaction 1 was made using the primer pair 5'GCGTGTCCCCAGGACAAGGT3' 5' ATGTCGCTGTCTTTCTCTTCTTCGC 3', amplifying a 124 base pairs (bp) DNA fragment found in both VACV-BR groups. Reaction 2 was made with the primer pair 5'GCGTGTCCCCAGGACAAGGT3' and 5'CTGGATGGGTCTTG3', amplifying a 138 bp DNA fragment of VACV-BR Group 2 viruses but not from VACV-BR Group 1 viruses. This is possible since the reverse primer anneals to a region present on group 2 but deleted in group 1. PCR conditions were 94 °C for 10 min, 30 cycles of 94 °C for 30 s, 50 °C for 30 s and 72 °C for 30 s, followed by 72 °C for 10 min. A semi-nested PCR was made for A26L, using the primer pair 5'ACCACGTCTACACTCGGCGA3' and 5'TGCATCGAGAGCGGAGGAGGA3' for the first reaction and the primer pair 5'ACCACGTCTACACTCGGCGA3' and 5'CGATGCCAAGTACATCGACGA3' for the second reaction. Amplicon sizes are 750 bp and 160 bp respectively, and reaction conditions were: 95 °C for 10 min, 30 cycles of 95 °C for 1 min, 60 °C for 1 min and 72 °C for 1 min, followed by 72 °C for 10 min. Reactions were carried out using 10–50 ng of DNA as template with 0.4 mM of each primer; 0.1 M tris-HCl (pH 8.4); 0.5M KCl; 2.5 mM MgCl2; 10 mM of each dNTPs and 2U Taq DNA polymerase to a final volume of 20 μL. Brazilian isolates GP1V and GP2V were used as positive controls. Reactions for the A56R gene were made as previously described [52]. PCR products were electrophoresed in 8%-PAGE and silver stained [51].
For sequencing, the A56R gene was amplified with primers and conditions described elsewhere [53], resulting in amplicons of approximately 900 bp. The A56R fragments were sequenced in both orientations and in duplicate in an automated DNA sequencer (ABI PRISM, Applied Biosystems, São Paulo, Brazil). The sequences were aligned with previously published OPV sequences from GenBank (accession numbers can be found in the figures) using MEGA software version 6.1 (Arizona State University, Phoenix, AZ, USA, 2014). The phylogenetic tree based on the A56R gene sequences was constructed by the Maximum Likelihood method based on the Tamura-Nei model implemented in MEGA6.1 (Arizona State University, Phoenix, AZ, USA, 2014). A56R sequences from clonal VACV-BR obtained in this study were deposited in GenBank. (BAB2V-Clone 1: KP282618; BAB2V-Clone 2: KP282619; BAB2V-Clone 3: KP282620; BAB2V Clone 4: KP282621; BAB2V-Clone 5: KP282622; MGHV-Clone 1: KP282623; MGHV-Clone 2: KP282624; MGHV-Clone 3: KP282625; MGHV-Clone 4: KP282626; MGHV-Clone 5: KP282627; BAHV-Clone 1: KP282628; BAHV-Clone 2: KP282629; BAHV-Clone 3: KP282630; BAHV-Clone 4: KP282631; BAHV-Clone 5: KP282632; ESBV-Clone 1: KP282633; ESBV-Clone 2: KP282634; ESBV-Clone 3: KP282635; ESBV-Clone 4: KP282636; ESBV-Clone 5: KP282637; GOB2V-Clone 1: KP282638; GOB2V-Clone 2: KP282639; GOB2V-Clone 3: KP282640; GOB2V-Clone 4: KP282641; GOB2V-Clone 5: KP282642; MGBV-Clone 1: KP282643; MGBV-Clone 2: KP282644; MGBV-Clone 3: KP282645; MGBV-Clone 4: KP282646; PAB2V-Clone 1: KP282647; PAB2V-Clone 2: KP282648; PAB2V-Clone 3: KP282649; PAB2V-Clone 4: KP282650; BABV-Clone 1: KP282651; BABV-Clone 2: KP282652; BABV-Clone 3: KP282653; BABV-Clone 4: KP282654; BABV-Clone 5: KP282655; GOBV-Clone 1: KP282656; GOBV-Clone 2: KP282657; GOBV-Clone 3: KP282658; GOBV-Clone 4: KP282659; GOBV-Clone 5: KP282660; PABV-Clone 1: KP282661; PABV-Clone 2: KP282662; PABV-Clone 3:KP282663; PABV-Clone 4: KP282664; PABV-Clone 5: KP282665).

3. Results

3.1. Viral Clone Isolation and Biological Assays

Upon inoculation onto Vero cells, the 10 clinical samples presented typical OPV CPEs, such as lytic plaques, cell agglomerates and formation of vacuoles were observed. No changes were observed in negative control cells. Then, a total of 48 viral clones from the 10 clinical samples were selected. It is important to emphasize that a higher prevalence of small plaques was observed during plaque selection. After three rounds of plaque purification, the selected clones were expanded and purified in sucrose cushion. Biological assays such as plaque phenotype, comet phenotype and growth curve assays were performed to compare the isolated clones. Plaque phenotype assays showed that the clones produced distinct plaques size (Figure 2 and Supplementary Figure S1). The VACV-BABV and VACV-GOBV revealed interesting results. Three clones isolated from the VACV-BABV (BABV-Clone 2, BABV-Clone 3 and BABV-Clone 4) showed a large-plaque profile , while two other clones (BABV-Clone 1 and BAB-Clone 5) presented a small-plaque profile; and the VACV-GOBV showed a single clone (GOBV-Clone 2) with a large-plaque profile and four clones (GOBV-Clone 1, GOBV-Clone 3, GOBV-Clone 4 and GOBV-Clone 5) with the small-plaque profile. These results demonstrate the co-circulation of VACV with at least two distinct plaque-size phenotypes in the same VACV population. The clones obtained from the other isolations showed only small plaque phenotype.
Figure 2. Viral clone plaque phenotype assays. BSC40 cells were cultured in a 6-well dish and then infected with VACV clones. Infection was carried out in the presence of 0.5% agarose for 48 h, followed by fixation and staining. All clones of the VACV-BAB2V, VACV-BAHV and VACV-GOB2V presented only small plaques, while VACV-BABV and VACV-GOBV strains presented small plaque clones as well as large plaque clones. VACV-BABV presented two small plaque clones (Clone 1 and 5) and three large plaque clones (Clone 2, 3 and 4). VACV-GOBV presented four small plaque clones (Clone 1, 3, 4 and 5) and one large plaque clone (Clone 2).
Figure 2. Viral clone plaque phenotype assays. BSC40 cells were cultured in a 6-well dish and then infected with VACV clones. Infection was carried out in the presence of 0.5% agarose for 48 h, followed by fixation and staining. All clones of the VACV-BAB2V, VACV-BAHV and VACV-GOB2V presented only small plaques, while VACV-BABV and VACV-GOBV strains presented small plaque clones as well as large plaque clones. VACV-BABV presented two small plaque clones (Clone 1 and 5) and three large plaque clones (Clone 2, 3 and 4). VACV-GOBV presented four small plaque clones (Clone 1, 3, 4 and 5) and one large plaque clone (Clone 2).
Viruses 07 01218 g002
To proceed with the biological characterization, the comet phenotype assays for clones were carried out. Prominent comets were observed only in the four large plaque clones (data not shown). In order to compare the growth of VACV clones, we performed a growth curve assay showing that small-plaque clones presented a similar replication profile, but distinct from that observed for large plaque clones, that exhibited an increase of 2–4 logs when compared to the small plaque clones (Figure 3). Thereby, plaque phenotype and growth curve assays showed that large plaque clones display higher replication capacity when compared to small plaque clones and that those different profiles can be found in clones of a same sample.
Figure 3. Growth curves assay of VACV-BR clones in BSC-40 cells. Cells were infected at an MOI of 0.01. Viral replication was determined by measuring the viral titer at 0, 12, 24 and 48 h after viral adsorption. The legend is the same for all the graphics: clone 1 (O), clone 2 (◊), clone 3 (*), clone 4 (□), and clone 5 (∆) when present. Each data point represents the average of two independent experiments. The analysis was performed using GraphPad Prism 6.0 software (Version 5.00, GraphPad Software, San Diego, CA, USA, 2012).
Figure 3. Growth curves assay of VACV-BR clones in BSC-40 cells. Cells were infected at an MOI of 0.01. Viral replication was determined by measuring the viral titer at 0, 12, 24 and 48 h after viral adsorption. The legend is the same for all the graphics: clone 1 (O), clone 2 (◊), clone 3 (*), clone 4 (□), and clone 5 (∆) when present. Each data point represents the average of two independent experiments. The analysis was performed using GraphPad Prism 6.0 software (Version 5.00, GraphPad Software, San Diego, CA, USA, 2012).
Viruses 07 01218 g003

3.2. Virulence of VACV-Clones in Balb/c Mice

For virulence assays, we selected the large plaque clones 2 and 3, and the small plaque clone 1 from VACV-BABV, and the large plaque clone 2 and small plaque clone 3 from VACV-GOBV. This selection was based on previous studies that correlated VACV plaque size and virulence in Balb/c mice [25]. Four-week-old BALB/c mice (n = 4) were intranasally infected with 1 × 106 pfu of each VACV clone and observed during 10 days for weight change, survival rate and clinical signs. The VACV-WR and PBS were used as virulence-positive and negative controls, respectively. We observed clinical signs, such as ruffling fur and arching back in mice inoculated with large plaque clones and VACV-WR within 3–7 dpi, similar to what has been previously reported for the virulent G2 VACV-BR [25]. However, some clinical signs, such as balanopostitis and periocular alopecia, were only observed in animals infected with VACV-WR (Figure 4C). Survival rates of mice inoculated with large plaque clones ranged from 0% (BABV-clone 2) to 50% (VACV-BABV-clone 3 and VACV-GOBV-clone 2), showing a variation in survival rate among virulent clones (Figure 4A). All mice inoculated with large plaque clones lost weight, while mice inoculated with small-plaque clones and with PBS gained weight over the 10-day experiment (Figure 4B). Moreover, mice inoculated with small-plaque clones, as well as the negative control group, did not present clinical signs and stayed alive during the 10 days of observation similar to what has been previously reported for the non-virulent VACV-BR group 1 [25].
Figure 4. Virulence assays in a BALB/c mice model. Groups of four-week-old mice (n = 4) were inoculated by intranasal route with 106 PFU/10 µL of VACV-BABV strain clones 2 and 3 (large-plaque) and clone 1 (small plaque), and VACV-GOBV clone 2 (large plaque) and clone 3 (small plaque). VACV-WR was used as virulent control, and the negative control group was inoculated with PBS. No mice infected with small plaque clones died while the mice infected with large plaque clones had survival rates of 0% and 50%, showing variations in survival rates between the clones (A); Uninfected mice (negative control) and animals infected with small plaque clones gained weight, while mice infected with large plaque clones lost weight similarly to mice infected with virulent control VACV-WR (B). Ruffling fur and arching back were observed in mice infected with large plaque clones on day 4 p.i. No clinical signals were observed in mice infected with small plaque clones from the same strain. Balanopostitis and periocular alopecy were observed only in mice infected with virulent control VACV-WR on day 3 p.i. No clinical signals were observed in control mice inoculated with PBS.
Figure 4. Virulence assays in a BALB/c mice model. Groups of four-week-old mice (n = 4) were inoculated by intranasal route with 106 PFU/10 µL of VACV-BABV strain clones 2 and 3 (large-plaque) and clone 1 (small plaque), and VACV-GOBV clone 2 (large plaque) and clone 3 (small plaque). VACV-WR was used as virulent control, and the negative control group was inoculated with PBS. No mice infected with small plaque clones died while the mice infected with large plaque clones had survival rates of 0% and 50%, showing variations in survival rates between the clones (A); Uninfected mice (negative control) and animals infected with small plaque clones gained weight, while mice infected with large plaque clones lost weight similarly to mice infected with virulent control VACV-WR (B). Ruffling fur and arching back were observed in mice infected with large plaque clones on day 4 p.i. No clinical signals were observed in mice infected with small plaque clones from the same strain. Balanopostitis and periocular alopecy were observed only in mice infected with virulent control VACV-WR on day 3 p.i. No clinical signals were observed in control mice inoculated with PBS.
Viruses 07 01218 g004

3.3. Phylogeny

For molecular comparison of the clones, the A56R, C23L, and A26L genes were amplified. Large plaque clones presented a typical G2 VACV-BR profile while the small plaque clones presented a profile similar to G1, which was corroborated by our biological data. However, two exceptions were found: large plaque clone 2 from VACV-BABV showed a G2 VACV-BR profile for the A56R and C23L genes and G1 profile for the A26L gene, and the small plaque clone 2 from the VACV-MGHV strain showed a G1 VACV-BR profile for the A56R and A26L genes and G2 profile to C23L gene. Table 2 summarizes these results. For phylogeny, we sequenced a partial fragment (900 bp) of the A56R gene of all VACV clones and the obtained sequences were compared to other VACV sequences available at NCBI nucleotide database. Phylogenetic analyses clustered the large plaque clones into G2 VACV-BR branch, while the small plaque clones were branched with G1 VACV-BR, further confirming the biological data (Figure 5). The branches were statistically supported by high bootstrap values, mainly due the presence of a molecular signature of 18-bp deletion in A56R gene from G1 of VACV-BR (Figure 6). Moreover, small plaque clones sequences showed specific polymorphisms that reflected the separation of some clones into different branches within this group (Figure 5 and Supplementary Figure S2).
Table 2. VACV-BR group 1 or 2 profile of the isolated clones based on A56R, A26L and C23L genes. Bold font indicates clones with mixed profile (Group 1 and 2).
Table 2. VACV-BR group 1 or 2 profile of the isolated clones based on A56R, A26L and C23L genes. Bold font indicates clones with mixed profile (Group 1 and 2).
ClonesGenesClonesGenes
A56RA26LC23LA56RA26LC23L
BABV-Clone 1Group 1Group 1Group 1GOB2V-Clone 1Group 1Group 1Group 1
BABV-Clone 2Group 2Group 1Group 2GOB2V-Clone 2Group 1Group 1Group 1
BABV-Clone 3Group 2Group 2Group 2GOB2V-Clone 3Group 1Group 1Group 1
BABV-Clone 4Group 2Group 2Group 2GOB2V-Clone 4Group 1Group 1Group 1
BABV-Clone 5Group 1Group 1Group 1GOB2V-Clone 5Group 1Group 1Group 1
BAB2V-Clone 1Group 1Group 1Group 1ESBV-Clone 1Group 1Group 1Group 1
BAB2V-Clone 2Group 1Group 1Group 1ESBV-Clone 2Group 1Group 1Group 1
BAB2V-Clone 3Group 1Group 1Group 1ESBV-Clone 3Group 1Group 1Group 1
BAB2V-Clone 4Group 1Group 1Group 1ESBV-Clone 4Group 1Group 1Group 1
BAB2V-Clone 5Group 1Group 1Group 1ESBV-Clone 5Group 1Group 1Group 1
BAHV-Clone 1Group 1Group 1Group 1MGHV-Clone 1Group 1Group 1Group 1
BAHV-Clone 2Group 1Group 1Group 1MGHV-Clone 2Group 1Group 1Group 2
BAHV-Clone 3Group 1Group 1Group 1MGHV-Clone 3Group 1Group 1Group 1
BAHV-Clone 4Group 1Group 1Group 1MGHV-Clone 4Group 1Group 1Group 1
BAHV-Clone 5Group 1Group 1Group 1MGHV-Clone 5Group 1Group 1Group 1
GOBV-Clone 1Group 1Group 1Group 1PABV-Clone 1Group 1Group 1Group 1
GOBV-Clone 2Group 2Group 2Group 2PABV-Clone 2Group 1Group 1Group 1
GOBV-Clone 3Group 1Group 1Group 1PABV-Clone 3Group 1Group 1Group 1
GOBV-Clone 4Group 1Group 1Group 1PABV-Clone 4Group 1Group 1Group 1
GOBV-Clone 5Group 1Group 1Group 1PABV-Clone 5Group 1Group 1Group 1
MGBV-Clone 1Group 1Group 1Group 1PAB2V-Clone 1Group 1Group 1Group 1
MGBV-Clone 2Group 1Group 1Group 1PAB2V-Clone 2Group 1Group 1Group 1
MGBV-Clone 3Group 1Group 1Group 1PAB2V-Clone 3Group 1Group 1Group 1
MGBV-Clone 4Group 1Group 1Group 1PAB2V-Clone 4Group 1Group 1Group 1
Figure 5. Phylogenetic analysis of viral clones based on the nucleotide sequences of Orthopoxvirus A56R gene. The tree was constructed by the Maximum Likelihood method with 1000 bootstrap replicates based in Tamura-Nei model implemented in MEGA6.1 software program (Arizona State University, Phoenix, AZ, USA, 2014) (www.megasoftware.net). Bootstrap confidence intervals are shown on branches. Nucleotide sequences were obtained from GenBank. Distinct colors indicate a strain and its clones as shown in the legend. Clones isolated from the same clinical sample were grouped in separate branches (Groups 1 and 2). Most of the clones were grouped in G1 and subgroups within the G1 were observed.
Figure 5. Phylogenetic analysis of viral clones based on the nucleotide sequences of Orthopoxvirus A56R gene. The tree was constructed by the Maximum Likelihood method with 1000 bootstrap replicates based in Tamura-Nei model implemented in MEGA6.1 software program (Arizona State University, Phoenix, AZ, USA, 2014) (www.megasoftware.net). Bootstrap confidence intervals are shown on branches. Nucleotide sequences were obtained from GenBank. Distinct colors indicate a strain and its clones as shown in the legend. Clones isolated from the same clinical sample were grouped in separate branches (Groups 1 and 2). Most of the clones were grouped in G1 and subgroups within the G1 were observed.
Viruses 07 01218 g005
Figure 6. Alignment of a fragment of viral clones A56R gene nucleotide sequences with other Orthopoxvirus sequences. The sequences were obtained from GenBank and aligned using the standard parameters of CLUSTAL W. Nucleotide positions are shown according to the VACV-MARV (G1 VACV-BR isolate). (.) indicates identity and (-) indicates deletion of nucleotides. The box highlights significant variations among clones from the same sample in which the clones 1 and 5 from the VACV-BABV showed a deletion of the 18 nucleotides observed in the A56R gene sequence of G1 VACV-BR, while clones 2, 3 and 4 from the same clinical sample did not show this deletion (G2). This alignment displays only representative clones from each sample, except for VACV-BABV and VACV-GOBV, in which all clones are represented.
Figure 6. Alignment of a fragment of viral clones A56R gene nucleotide sequences with other Orthopoxvirus sequences. The sequences were obtained from GenBank and aligned using the standard parameters of CLUSTAL W. Nucleotide positions are shown according to the VACV-MARV (G1 VACV-BR isolate). (.) indicates identity and (-) indicates deletion of nucleotides. The box highlights significant variations among clones from the same sample in which the clones 1 and 5 from the VACV-BABV showed a deletion of the 18 nucleotides observed in the A56R gene sequence of G1 VACV-BR, while clones 2, 3 and 4 from the same clinical sample did not show this deletion (G2). This alignment displays only representative clones from each sample, except for VACV-BABV and VACV-GOBV, in which all clones are represented.
Viruses 07 01218 g006

4. Discussion

Bovine vaccinia outbreaks have been increasingly reported in Brazil since 1999 causing economic, social, veterinary and public health problems [11]. Despite the current and historical importance of VACV, little information is available about its circulation, prevalence, origins, maintenance in environment and natural reservoirs.
In this study, we were able to demonstrate the genetic and biologic variability among VACV and even among clones isolated from the same sample. Previous studies have demonstrated the genetic variability present in vaccine preparations of VACV, but in this study we observed the same profile among viruses isolated during BV outbreaks. Our study reinforced the circulation of the two groups (G1 and 2) of VACV-BR, as previously reported [24,26,27,28]. In addition, our data demonstrate the co-circulation of mixed VACV populations in the same sample, as well as a great genetic diversity among VACV-BR G1 clones in the same sample. In addition, we isolated two clones (VACV-BABV-clone 2 and VACV-MGHV clone 2) that presented a mosaic profile, with both G1 and G2 features based on the molecular analysis of A56R, A26L and C23L genes. These clones may have arisen as product of a recombination event.
Previous studies have shown a high recombination rate during a poxvirus co-infection [46,47,48]. Co-circulation of the genetically different clones in the same sample as shown here may favor the recombination generating new genetic variants which can also have influence in virus evolution and pathogenesis. Moreover, genetic changes can favor adaptation to new hosts, such as described for other viruses of the same genus contributing to expand its circulation and host range [5].
Despite the existence of two VACV groups circulating in Brazil, generally G1 viruses has been more frequently isolated when compared to G2 viruses, including this study in which 92% of the isolated clones were grouped in VACV-BR-G1 while only 8% were grouped in VACV-BR-G2 were based on A56R gene. However, we believe that in addition to the higher prevalence of G1 VACV-BR there is a laboratory filter that favors the isolation and characterization of G1 VACV-BR, since the cellular systems commonly used in isolation and characterization, like embryonated chicken eggs and Vero cells, are unable to highlight differences between G1 and 2 such as plaque size variations. Probably, in a mixed sample (which has G1 and 2 viruses), the prevalence of G2 may be lower as was demonstrated in the case of VACV-GOBV sample, which would not allow the identification of the G2 viruses if the sample is directly processed and characterized by conducting PCR from VACV without previous plaque purification in an appropriate cellular system for the identification of both groups. Thus, we suggest that the viruses that have been isolated represent a small part of the diversity of VACV present in the natural environment. Our findings raise more discussions about the origin of VACV in Brazil, reinforcing the notion that VACV-BR viruses may have more than one origin [54].
We confirmed the virulence differences between large plaque and small plaque viruses as demonstrated previously [25,28]. While small plaque viral clones did not cause weight loss or death during infection in BALB/c mice, large plaque clones caused weight loss, clinical signs and death of infected BALB/C mice, showing that virulent and no virulent clones may be present in the same sample.
Here, we also showed that the replication rates varies among clones of the same sample with large plaque clones showing a higher replication rates when compared to small plaque clones isolated from the same sample in BSC40 cells. The importance of the higher replication rates for VACV evolution and maintenance in the environment needs further studies.
Taken together, data regarding G2 clones show exclusive biological and molecular features (including genes profiles, plaque phenotype and mouse clinical signs) distinct to other known large plaque reference samples, including VACV-Western Reserve.
In summary, our study demonstrated the genetic and biological variability between clones isolated from natural circulating VACV-BR in BV outbreaks, evidencing that VACV can contain viral clone subpopulations. These results raise new questions about the diversity and origin of VACV circulating in Brazilian natural and rural environments. Nevertheless, the real implications of different VACV clones’ co-circulation are unknown but we can emphasize the significance in the evolution and maintenance of these viruses in the environment.

Supplementary Files

Supplementary File 1Supplementary File 2

Acknowledgments

We thank colleagues from Laboratório de Vírus (ICB-UFMG). For assistance in collection of clinical specimens we thank Instituto Mineiro de Agropecuária(IMA), Laboratório de Virologia Animal da Escola de Medicina Veterinária of UFMG and Empresa Brasileira de Pesquisa Agropecuária.
(EMBRAPA). We thank Juliano Leal de Paula from Núcleo de Análise de Genoma e Expressão Gênica (NAGE)-Departamento de Bioquímica e Imunologia (ICB—UFMG). We thank Pró-Reitoria de Pesquisa da Universidade Federal de Minas Gerais (PRPq-UFMG) and support agencies (CNPq, CAPES and FAPEMIG) for the financial support.

Author Contributions

Conceived and designed the experiments: G.O., F.A., G.A., G.T., E.K., J.A. Performed the experiments: G.O., M.L., A.C.A., R.C. Analyzed the data: G.O., J.A., J.A. Contributed reagents/materials/analysis tools: F.A., P.C.F., E.K. Wrote the paper: G.O., J.A. All authors read and approved final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References and Notes

  1. Fenner, F.; Henderson, D.A.; Arita, I.; Jezek, A.; Ladnyi, I.D. Smallpox and Its Eradication; World Health Organization Press: Geneva, Switzerland, 1988. [Google Scholar]
  2. Chen, H.; Chuai, X.; Deng, Y.; Wen, B.; Wang, W.; Xiong, S.; Ruan, L.; Tan, W. Optimisation of prime-boost immunization in mice using novel protein-based and recombinant vaccinia (Tiantan)-based HBV vaccine. PLoS One 2012, 7, e43730. [Google Scholar] [CrossRef] [PubMed]
  3. Deng, Y.; Zhang, K.; Tan, W.; Wang, Y.; Chen, H.; Wu, X.; Ruan, L. A recombinant DNA and vaccinia virus prime-boost regimen induces potent long-term T-cell responses to HCV in BALB/c mice. Vaccine 2009, 27, 2085–2088. [Google Scholar] [CrossRef] [PubMed]
  4. Liu, C.; Du, S.; Li, C.; Wang, Y.; Wang, M.; Li, Y.; Yin, R.; Li, X.; Ren, D.; Qin, Y.; et al. Immunogenicity analysis following human immunodeficiency virus recombinant DNA and recombinant vaccinia virus Tian Tan prime-boost immunization. Sci. China Life Sci. 2013, 56, 531–540. [Google Scholar] [CrossRef] [PubMed]
  5. Reynolds, M.G.; Carroll, D.S.; Karem, K.L. Factors affecting the likelihood of monkeypox’s emergence and spread in the post-smallpox era. Curr. Opin. Virol. 2012, 2, 335–343. [Google Scholar] [CrossRef] [PubMed]
  6. Pahlitzsch, R.; Hammarin, A.L.; Widell, A. A case of facial cellulitis and necrotizing lymphadenitis due to cowpox virus infection. Clin. Infect. Dis. 2006, 43, 737–742. [Google Scholar] [CrossRef] [PubMed]
  7. Karem, K.L.; Reynolds, M.; Braden, Z.; Lou, G.; Bernard, N.; Patton, J.; Damon, I.K. Characterization of acute-phase humoral immunity to monkeypox: Use ofimmunoglobulin M enzyme-linked immunosorbent assay for detection of monkeypox infection during the 2003 North American outbreak. Clin. Diagn Lab. Immunol. 2005, 12, 867–872. [Google Scholar] [PubMed]
  8. Rimoin, A.W.; Mulembakani, P.M.; Johnston, S.C.; Lloyd Smith, J.O.; Kisalu, N.K.; Kinkela, T.L.; Blumberg, S.; Thomassen, H.A.; Pike, B.L.; Fair, J.N.; et al. Major increase in human monkeypox incidence 30 years after smallpox vaccination campaigns cease in the Democratic Republic of Congo. Proc. Natl. Acad. Sci. USA 2010, 107, 16262–16267. [Google Scholar] [CrossRef] [PubMed]
  9. Singh, R.K.; Hosamani, M.; Balamurugan, V.; Bhanuprakash, V.; Rasool, T.J.; Yadav, M.P. Buffalopox: An emerging and re-emerging zoonosis. Anim. Health Res. Rev. 2007, 8, 105–114. [Google Scholar] [CrossRef] [PubMed]
  10. Trindade, G.S.; Guedes, M.I.; Drumond, B.P.; Mota, B.E.; Abrahao, J.S.; Lobato, Z.I.; Gomes, J.A.; Corrêa-Oliveira, R.; Nogueira, M.L.; Kroon, E.G.; et al. Zoonotic vaccinia virus: Clinical and immunological characteristics in a naturally infected patient. Clin. Infect. Dis. 2009, 48, 37–40. [Google Scholar] [CrossRef]
  11. Kroon, E.G.; Mota, B.E.; Abrahão, J.S.; da Fonseca, F.G.; de Souza Trindade, G. Zoonotic Brazilian Vaccinia virus: from field to therapy. Antivir. Res. 2011, 92, 150–163. [Google Scholar] [CrossRef] [PubMed]
  12. Oliveira, D.B.; Assis, F.L.; Ferreira, P.C.; Bonjardim, C.A.; de Souza Trindade, G.; Kroon, E.G.; Abrahão, J.S. Short Report: Group 1 Vaccinia virus Zoonotic Outbreak in Maranhão State, Brazil. Am. J. Trop. Med. Hyg. 2013, 89, 1142–1145. [Google Scholar] [CrossRef] [PubMed]
  13. Kalthoff, D.; Bock, W.I.; Hühn, F.; Beer, M.; Hoffmann, B. Fatal cowpox virus infection in cotton-top tamarins (Saguinus oedipus) in Germany. Vector Borne Zoonotic Dis. 2014, 14, 303–305. [Google Scholar] [CrossRef] [PubMed]
  14. Fernandes, A.T.; Travassos, C.E.; Ferreira, J.M.; Abrahao, J.S.; Rocha, E.S.; Ferreira, F.V.; dos Santos, J.R.; Bonjardim, C.A.; Ferreira, P.C.; Kroon, E.G. Natural human infections with Vaccinia virus during bovine vaccinia outbreaks. J. Clin. Virol. 2009, 44, 308–313. [Google Scholar] [CrossRef] [PubMed]
  15. Megid, J.; Appolinário, C.M.; Langoni, H.; Pituco, E.M.; Okuda, L.H. Vaccinia virus in humans and cattle in southwest region of Sao Paulo state, Brazil. Am. J. Trop. Med. Hyg. 2008, 79, 647–665. [Google Scholar] [PubMed]
  16. Damon, I.K. Poxviruses. In Fields Virology, 6th ed.; Knipe, D.M., Howley, P.M., Eds.; Lippincott, Williams and Wilkins: Philadelphia, PA, USA, 2014; Volume 2, p. 2160. [Google Scholar]
  17. Batista, V.H.; Scremin, J.; Aguiar, L.M.; Schatzmayr, H.G. Vulvar infection and possible human-to-human transmission of bovine poxvirus disease. Virus Rev. Res. 2009, 14, 1–10. [Google Scholar]
  18. Oliveira, G.P.; Fernandes, A.T.; Assis, F.L.; Alves, P.A.; Franco Luiz, A.P.M.; Figueiredo, L.B.; Costa de Almeida, C.M.; Travassos, C.E.P.F.; Trindade, G.S.; Abrahão, J.S.; et al. Intrafamilial transmission of Vaccinia virus during a bovine Vaccinia outbreak in Brazil: A new insight in viral transmission chain. Am. J. Trop. Med. Hyg. 2014, 90, 1021–1023. [Google Scholar] [CrossRef] [PubMed]
  19. Lopes, O.S.; Lacerda, J.P.; Fonseca, I.E.; Castro, D.P.; Forattini, O.P.; Rabello, E.X. CotiaVirus: A New Agent Isolated from Sentinel Mice in Sao Paulo, Brazil. Am. J. Trop. Med. Hyg. 1965, 14, 156–157. [Google Scholar] [PubMed]
  20. Fonseca, F.G.; Lanna, M.C.; Campos, M.A.; Kitajima, E.W.; Peres, J.N.; Golgher, R.R.; Ferreira, P.C.; Kroon, E.G. Morphological and molecular characterization of the poxvirus BeAn 58058. Arch. Virol. 1998, 143, 1171–1186. [Google Scholar] [CrossRef] [PubMed]
  21. Fonseca, F.G.; Trindade, G.S.; Silva, R.L.; Bonjardim, C.A.; Ferreira, P.C.; Kroon, E.G. Characterization of a vaccinia-like virus isolated in a Brazilian forest. J. Gen. Virol. 2002, 83, 223–228. [Google Scholar] [PubMed]
  22. Damaso, C.R.; Esposito, J.J.; Condit, R.C.; Moussatche, N. An emergent poxvirus from humans and cattle in Rio de Janeiro State: Cantagalo virus may derive from Brazilian smallpox vaccine. Virology 2000, 277, 439–449. [Google Scholar] [CrossRef] [PubMed]
  23. Trindade, G.S.; da Fonseca, F.G.; Marques, J.T.; Diniz, S.; Leite, J.A.; de Bodt, S.; van der Peer, Y.; Bonjardim, C.A.; Ferreira, P.C.; Kroon, E.G. Belo Horizonte virus: a vaccinia-like virus lacking the A-type inclusion body gene isolated from infected mice. J. Gen. Virol. 2004, 85, 2015–2021. [Google Scholar] [CrossRef] [PubMed]
  24. Assis, F.L.; Borges, I.A.; Ferreira, P.C.; Bonjardim, C.A.; Trindade, G.S.; Lobato, Z.I.; Guedes, M.I.; Mesquita, V.; Kroon, E.G.; Abrahão, J.S. Group 2 vaccinia virus, Brazil. Emerg. Infect. Dis. 2012, 18, 2035–2038. [Google Scholar] [CrossRef] [PubMed]
  25. Ferreira, J.M.; Drumond, B.P.; Guedes, M.I.; Pascoal-Xavier, M.A.; Almeida-Leite, C.M.; Arantes, R.M.; Mota, B.E.; Abrahão, J.S.; Alves, P.A.; Oliveira, F.M.; et al. Virulence in Murine Model Shows the Existence of Two Distinct Populations of Brazilian Vaccinia virus Strains. PLoS One 2008, 3, e3043. [Google Scholar] [CrossRef] [PubMed]
  26. Trindade, G.S.; Lobato, Z.I.; Drumond, B.P.; Leite, J.A.; Trigueiro, R.C.; Guedes, M.I.; Fonseca, F.G.; Santos, J.R.; Bonjardim, C.A.; Ferreira, P.C.; Kroon, E.G. Shortreport: Isolation of two vaccinia virus strains froma single bovine vaccinia outbreak in rural area fromBrazil: Implications on the emergence of zoonotic orthopoxviruses. Am. J. Trop. Med. Hyg. 2006, 75, 486–490. [Google Scholar] [PubMed]
  27. Drumond, B.P.; Leite, J.A.; da Fonseca, F.G.; Bonjardim, C.A.; Ferreira, P.C.; Kroon, E.G. Brazilian Vaccinia virus strains are genetically divergent and differ from the Lister vaccine strain. Microbes Infect. 2008, 10, 185–197. [Google Scholar] [CrossRef] [PubMed]
  28. Campos, R.K.; Brum, M.C.; Nogueira, C.E.; Drumond, B.P.; Alves, P.A.; Siqueira-Lima, L.; Assis, F.L.; Trindade, G.S.; Bonjardim, C.A.; Ferreira, P.C.; et al. Assessing the variability of Brazilian Vaccinia virus isolates from a horse exanthematic lesion: Coinfection with distinct viruses. Arch. Virol. 2011, 156, 275–283. [Google Scholar] [CrossRef] [PubMed]
  29. Leite, J.A.; Drumond, B.P.; de Souza Trindade, G.; Bonjardim, C.A.; Ferreira, P.C.; Kroon, E.G. Brazilian Vaccinia virus strains show genetic polymorphism at the ati gene. Virus Genes 2007, 35, 531–539. [Google Scholar] [CrossRef] [PubMed]
  30. Assis, F.L.; Almeida, G.M.; Oliveira, D.B.; Franco-Luiz, A.P.; Campos, R.K.; Guedes, M.I.; Fonseca, F.G.; Trindade, G.S.; Drumond, B.P.; Kroon, E.G.; et al. Characterization of a new Vaccinia virus isolate reveals the C23L gene as a putative genetic marker for autochthonous Group 1 Brazilian Vaccinia virus. PLoS One 2012, 11, e50413. [Google Scholar] [CrossRef]
  31. Trindade, G.S.; da Fonseca, F.G.; Marques, J.T.; Nogueira, M.L.; Mendes, L.C.; Borges, A.S.; Peiró, J.R.; Pituco, E.M.; Bonjardim, C.A.; Ferreira, P.C.; et al. Araçatuba virus: A vaccinialike virus associated with infection in humans and cattle. Emerg Infect. Dis. 2003, 9, 155–160. [Google Scholar] [CrossRef] [PubMed]
  32. Abrahão, J.S.; Guedes, M.I.; Trindade, G.S.; Fonseca, F.G.; Campos, R.K.; Mota, B.F.; Lobato, Z.I.; Fernandes, A.T.; Rodrigues, G.O.; Lima, L.S.; et al. One more piece in the VACV ecological puzzle: Could peridomestic rodents be the link between wildlife and bovine vaccinia outbreaks in Brazil? PLoS One 2009, 19, e7428. [Google Scholar] [CrossRef]
  33. Reed, K.D.; Melski, J.W.; Graham, M.B.; Regnery, R.L.; Sotir, M.J.; Wegner, M.V.; Kazmierczak, J.J.; Stratman, E.J.; Li, Y.; Fairley, J.A.; et al. The detection of monkeypox in humans in the Western Hemisphere. N. Engl. J. Med. 2004, 350, 342–350. [Google Scholar] [CrossRef] [PubMed]
  34. Gubser, C.; Hue, S.; Kellam, P.; Smith, G.L. Poxvirus genomes: A phylogenetic analysis. J. Gen. Virol. 2004, 85, 105–117. [Google Scholar] [CrossRef] [PubMed]
  35. Chen, N.; Li, G.; Liszewski, M.K.; Atkinson, J.P.; Jahrling, P.B.; Feng, Z.; Schriewer, J.; Buck, C.; Wang, C.; Lefkowitz, E.J.; et al. Virulence differences between monkeypox virus isolates from West Africa and the Congo basin. Virology 2005, 340, 46–63. [Google Scholar] [CrossRef] [PubMed]
  36. Esposito, J.J.; Sammons, S.A.; Frace, A.M.; Osborne, J.D.; Olsen-Rasmussen, M.; Zhang, M.; Govil, D.; Damon, I.K.; Kline, R.; Laker, M.; et al. Genome sequence diversity and clues to the evolution of variola (smallpox) virus. Science 2006, 313, 807–812. [Google Scholar] [CrossRef] [PubMed]
  37. Emerson, G.L.; Li, Y.; Frace, M.A.; Olsen-Rasmussen, M.A.; Khristova, M.L.; Govil, D.; Sammons, S.A.; Regnery, R.L.; Karem, K.L.; Damon, I.K.; et al. The phylogenetics and ecology of the orthopoxviruses endemic to North America. PLoS One 2009, 4, e7666. [Google Scholar] [CrossRef] [PubMed]
  38. Gubser, C.; Hue, S.; Kellam, P.; Smith, G.L.; Hendrickson, R.C.; Wang, C.; Hatcher, E.L.; Lefkowitz, E.J. Orthopoxvirus genome evolution: The role of gene loss. Viruses 2010, 2, 1933–1967. [Google Scholar] [CrossRef] [PubMed]
  39. Kugelman, J.R.; Johnston, S.C.; Mulembakani, P.M.; Kisalu, N.; Lee, M.S.; Koroleva, G.; McCarthy, S.E.; Gestole, M.C.; Wolfe, N.D.; Fair, J.N.; et al. Genomic variability of monkeypox virus among humans, Democratic Republic of the Congo. Emerg. Infect. Dis. 2014, 20, 232–239. [Google Scholar] [CrossRef] [PubMed]
  40. Okeke, M.I.; Okoli, A.S.; Nilssen, Ø.; Moens, U.; Tryland, M.; Bøhn, T.; Traavik, T. Molecular characterization and phylogenetics of Fennoscandian cowpox virus isolates based on the p4c and atip genes. Virol. J. 2014, 11, 119–120. [Google Scholar] [CrossRef] [PubMed]
  41. Osborne, J.D.; Da Silva, M.; Frace, A.M.; Sammons, S.A.; Olsen-Rasmussen, M.; Upton, C.; Buller, R.M.; Chen, N.; Feng, Z.; Roper, R.L.; et al. Genomic differences of Vaccinia virus clones from Dryvax smallpox vaccine: The Dryvax-like ACAM2000 and the mouse neurovirulent Clone-3. Vaccine 2007, 25, 8807–8832. [Google Scholar] [CrossRef] [PubMed]
  42. Garcel, A.; Perino, J.; Crance, J.M.; Drillien, R.; Garin, D.; Favier, A.L. Phenotypic and genetic diversity ofthe traditional Lister smallpox vaccine. Vaccine 2009, 27, 708–717. [Google Scholar] [CrossRef] [PubMed]
  43. Morikawa, S.; Sakiyama, T.; Hasegawa, H.; Saijo, M.; Maeda, A.; Kurane, I.; Maeno, G.; Kimura, J.; Hirama, C.; Yoshida, T.; et al. An attenuated LC16m8 smallpox vaccine: Analysis of fullgenome sequence and induction of immune protection. J. Virol. 2005, 79, 11873–11891. [Google Scholar] [CrossRef] [PubMed]
  44. Qin, L.; Upton, C.; Hazes, B.; Evans, D.H. Genomic analysis of the vaccinia virus strain variants found in Dryvax vaccine. J. Virol. 2011, 85, 13049–13060. [Google Scholar] [CrossRef] [PubMed]
  45. Qin, L.; Liang, M.; Evans, D.H. Genomic analysis of vaccinia virus strain TianTan provides new insights into the evolution and evolutionary relationships between Orthopoxviruses. Virology 2013, 442, 59–66. [Google Scholar] [CrossRef] [PubMed]
  46. Fenner, F.; Comben, B.M. Genetic studies with mammalian poxviruses. I.Demonstration of recombination between two strains of vaccina virus. Virology 1958, 5, 530–548. [Google Scholar] [CrossRef] [PubMed]
  47. Bedson, H.S.; Dumbell, K.R. Hybrids Derived from the Viruses of Variola Major and Cowpox. J. Hyg. 1964, 62, 147–158. [Google Scholar] [CrossRef] [PubMed]
  48. Qin, L.; Evans, D.H. Genome scale patterns of recombination between coinfecting vaccinia viruses. J. Virol. 2014, 88, 5277–5286. [Google Scholar] [CrossRef] [PubMed]
  49. Campos, M.A.S.; Kroon, E.G. Critical period for reversible block of vaccinia virus replication. Rev. Braz. Microbiol. 1993, 24, 104–110. [Google Scholar]
  50. Joklik, W.K. The purification of four strains of poxvirus. Virology 1962, 18, 9–18. [Google Scholar] [CrossRef] [PubMed]
  51. Sambrook, J.; Russell, D.W. Molecular Cloning: A Laboratory Manual, 3rd ed.; Nolan, C., Ford, N., Eds.; Cold Spring Harbor Laboratory Press: New York, NY, USA, 2001. [Google Scholar]
  52. Trindade, G.S.; Li, Y.; Olson, V.A.; Emerson, G.; Regnery, R.L.; da Fonseca, F.G.; Kroon, E.G.; Damon, I.K. Real-time PCR assay to identify variants of Vaccinia virus: Implications for the diagnosis of bovine vaccinia in Brazil. J. Virol. Methods 2008, 152, 63–71. [Google Scholar] [CrossRef] [PubMed]
  53. Ropp, S.L.; Jin, Q.; Knight, J.C.; Massung, R.F.; Esposito, J.J. PCR strategy for identification and differentiation of small pox and other Orthopoxviruses. J. Clin. Microbiol. 1995, 33, 2069–2076. [Google Scholar] [PubMed]
  54. Trindade, G.S.; Emerson, G.L.; Carroll, D.S.; Kroon, E.G.; Damon, I.K. Brazilian Vaccinia Viruses and Their Origins. Emerg. Infect. Dis. 2007, 13, 965–972. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Oliveira, G.; Assis, F.; Almeida, G.; Albarnaz, J.; Lima, M.; Andrade, A.C.; Calixto, R.; Oliveira, C.; Diomedes Neto, J.; Trindade, G.; et al. From Lesions to Viral Clones: Biological and Molecular Diversity amongst Autochthonous Brazilian Vaccinia Virus. Viruses 2015, 7, 1218-1237. https://0-doi-org.brum.beds.ac.uk/10.3390/v7031218

AMA Style

Oliveira G, Assis F, Almeida G, Albarnaz J, Lima M, Andrade AC, Calixto R, Oliveira C, Diomedes Neto J, Trindade G, et al. From Lesions to Viral Clones: Biological and Molecular Diversity amongst Autochthonous Brazilian Vaccinia Virus. Viruses. 2015; 7(3):1218-1237. https://0-doi-org.brum.beds.ac.uk/10.3390/v7031218

Chicago/Turabian Style

Oliveira, Graziele, Felipe Assis, Gabriel Almeida, Jonas Albarnaz, Maurício Lima, Ana Cláudia Andrade, Rafael Calixto, Cairo Oliveira, José Diomedes Neto, Giliane Trindade, and et al. 2015. "From Lesions to Viral Clones: Biological and Molecular Diversity amongst Autochthonous Brazilian Vaccinia Virus" Viruses 7, no. 3: 1218-1237. https://0-doi-org.brum.beds.ac.uk/10.3390/v7031218

Article Metrics

Back to TopTop