Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2006, p. 5631-5636, Vol. 72, No. 8
0099-2240/06/$08.00+0 doi:10.1128/AEM.00440-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
| SHORT REPORT |
The G. Eliava Institute of Bacteriophage, Microbiology and Virology, 0160 Tbilisi, Georgia,1 Dipartimento di Scienze e Tecnologie Alimentari e Microbiologiche, Università degli Studi di Milano, 20133 Milano, Italy,2 Health Protection Agency, Centre for Emergency Preparedness and Response, Porton Down, Salisbury SP4 OJG, United Kingdom,3 Department of Communicable Disease Surveillance and Response, World Health Organization, CH-1211 Geneva 27, Switzerland,4 Istituto Nazionale per le Malattie Infettive I.R.C.C.S. Lazzaro Spallanzani, 00142 Rome, Italy5
Received 23 February 2006/ Accepted 30 May 2006
|
|
|---|
|
|
|---|
In recent years efforts to evaluate the diversity of B. anthracis isolates from different geographic areas improved our knowledge of the distribution of this species throughout the world (3, 8, 9, 14, 18, 22, 25). These studies were also undertaken to trace strains deliberately released into the environment back to their origins (10, 15). Most of these studies took advantage of multiple-locus variable-number tandem repeat (VNTR) analysis (MLVA) (14), which is recognized as the "gold standard" for molecular typing of B. anthracis strains. MLVA is based on analysis of VNTRs in eight loci, including six chromosomal loci, vrrA, vrrB1, vrrB2, vrrC1, vrrC2, and CG3, and two plasmid loci, pXO1-aat and pXO2-at. VNTR markers have high genetic stability, which makes MLVA highly reproducible (10) and thus permits direct comparison of data from different laboratories. Following the original work of Keim et al. (13, 14), a number of studies addressed the biogeography of B. anthracis, characterizing strains originating from Kruger National Park in South Africa (25), France (8, 18), Poland (9), China (3), Italy (7), and Korea (22) and strains isolated following terrorist attacks in Japan (15) and the United States (10). These studies provided a clearer picture of the global spread and regional diversity of B. anthracis.
In the Caucasus, anthrax can be considered endemic, and the frequency of outbreaks has increased dramatically in the last 10 years according to the World Anthrax Data Site (http://www.vetmed.lsu.edu/whocc/georgia.htm). Georgia has been considered one of the countries most affected by outbreaks, including human cases (11). However, only a few sporadic reports have documented the occurrence of B. anthracis in the Caucasus (11; http://www.vetmed.lsu.edu/whocc/georgia.htm; http://www.ncdc.ge/Eng/statistics/Communicable%20Diseases%202004_files/frame .html), and very little information and incomplete data on the diversity of B. anthracis in this region are available.
In the other republics of the former Soviet Union (FSU), a similar situation exists, and in general, few strains have been characterized previously using MLVA and other molecular methods. An exception is the strains involved in the incident that occurred in Sverdlovsk (now Ekaterinburg, Russia) in 1979 (12).
There have also not been careful molecular characterizations of strains used in the FSU for both human and veterinary vaccines. Molecular characterization of the vaccine strains that were used in the FSU and are still used in many countries of the FSU would provide further data for evaluating the diversity of the isolates that are used worldwide for production of anthrax vaccines (14).
In the present study MLVA and other molecular typing methods were used to investigate B. anthracis diversity in Georgia and to characterize vaccine strains used in the FSU.
|
|
|---|
From August 2002 to May 2003 samples were obtained from several sites in eight different areas of Georgia for isolation of B. anthracis (Table 1). The sites sampled covered a 300-km transect in Georgia, including valley and mountain areas. The sites selected were considered reservoirs of B. anthracis in which the organism is endemic, and in these areas several cases were recorded in the past 5 to 7 years. Three blood samples (from two humans and a goat) and a total of 239 soil samples were collected from animal burial sites, abattoirs, sites where animals died, and animal stalls. Because of the high levels of infection in the Gardabani area, 70 samples were collected there. Sampling expeditions that resulted in isolation of B. anthracis were sent to eastern Georgia in the Gachiani, Gamarjveba, and Gardabani areas; in central Georgia in the Abisi area; and in western Georgia in the Rukhi, Zugdidi, Senaki and Kutaisi areas (Table 1).
|
View this table: [in a new window] |
TABLE 1. Characteristics of B. anthracis isolates used in this study, including Georgian isolates and vaccine strains used in the FSU
|
Notably, control samples collected from the Akhaltsikhe area in southern Georgia, where cases of B. anthracis have never been recorded, did not yield B. anthracis isolates (data not shown).
|
|
|---|
Strain STI-1 was used in the FSU for preparation of human vaccine (24, 27). The other strains were considered strains for veterinary use.
|
|
|---|
|
|
|---|
The pag, cya, lef, and cap virulence genes were detected in all 18 newly isolated strains and in vaccine strain 71/12. The other four vaccine strains lacked the cap gene located on the pXO2 plasmid (data not shown). No PCR products were obtained for the remaining isolates of the initial strain collection when sequences were amplified with the primer sets for the four virulence genes. When analyzed to determine the presence of the SG749 marker, all the virulent and vaccine strains produced the typical AluI restriction profile of B. anthracis, i.e., two restriction fragments at 660 bp and 90 bp (4). All of the new isolates and the vaccine strains exhibited the same intergenic transcribed spacer homoduplex-heteroduplex polymorphism profile typical of B. anthracis, with intergenic transcribed spacer homoduplex bands at about 240 and 460 bp and a series of slowly migrating heteroduplex bands (5). The strains were characterized by repetitive element polymorphism PCR with the BOX-A1R primer, and all of them produced major bands at approximately 130, 400, 630, and 770 bp, as previously shown for B. anthracis (2).
|
|
|---|
PCR products were analyzed with an ABI Prism 310 capillary sequencer (Applied Biosystems). One microliter of each of the four PCR mixtures was mixed with 1 µl of the 1,000-bp ROX internal size standard (Applied Biosystems) labeled with the phosphoramidite dye 6-carboxyrhodamine. Fourteen microliters of deionized formamide was added, and the mixture was denatured at 95°C for 10 min and cooled in an ice bath. For separation of the PCR products, a capillary (47 cm by 50 µm) filled with 4% performance-optimized polymer (Applied Biosystems) was used. The samples were electrophoresed under standard ABI 310 denaturing electrophoresis conditions for 45 min at 60°C. The data were analyzed using the GeneScan 3.1 software program (Applied Biosystems). The output of the analysis was a series of peaks (an electropherogram), the sizes of which were estimated by comparison with the fragments of the internal size standard. A threshold of 50 fluorescent units was used, and analysis of the size standard was performed by the local Southern method and with light data smoothing. If the baseline varied inconsistently, the sample was rerun. The results were compared to B. anthracis genotypes described by Keim et al. (14).
Diversity group clusters were analyzed using the sequential agglomerative hierarchical nested cluster analysis and unweighted pair group method analysis using average linkage clustering routines in the NTSYSpc 3.0 software package (21). Although the isolates lacking pXO2 were subjected to genotyping with seven MLVA markers, the data obtained were not included in the cluster analysis due to incompatibility of the marker number tested. A simple percent distance similarity matrix of the allele values for the isolates was constructed prior to cluster analysis. Individual marker diversity (Simpson diversity index) was calculated by using 1
(allele frequency)2 (14) and was based on the entire virulent isolate collection used in the study.
The sizes of the alleles detected by MLVA in all the strains analyzed are listed for each of the eight VNTR markers in Table 1, along with the designations of the global diversity groups, based on the analysis performed by Keim et al. (14). Based on marker allele combinations eight genotypes, designated GG I to GG VII in Table 1, were found among the strains analyzed. Table 1 also shows strain genotypes (GK genotypes) defined by Keim et al. (14).
Three distinct MLVA genotypes (GG I to GG III) were found for the Georgian isolates, and all of them belonged to cluster A3.a in the clustering system of Keim et al. (14). In the Georgian isolates, the informative VNTR markers were those from pXO1 and pXO2. These two markers are considered the most discriminatory markers for describing the diversity of B. anthracis (14). The diversity indexes of the pXO1 and pXO2 markers for the Georgian strains were 0.43, and 0.21, respectively. No diversity was observed for the other six VNTR markers.
Most of the isolates (13 of 18) represented the rare genotype GK 44, previously reported for only three strains found in Turkey and Namibia (14). The genotype of three other Georgian isolates was genotype GK 35, which has been found previously in Turkey, Namibia, the United Kingdom, and the United States. The three isolates with this genotype originated from different sources in one place (Gardabani) (Table 1); two were obtained from different animal burial sites, and one was obtained from the blood of a dead goat. The two other isolates represented the GK 40 genotype previously found only in Africa (in particular, Namibia, Mozambique, and Zimbabwe).
To describe the genetic relationship between genotypes found in the Turkish-Southern Caucasian region, a mutational step model, based on the minimum number of mutations in the MLVA loci needed to create the diversity observed, was developed (Fig. 1). All the GK genotypes belonging to cluster A3.a previously found in Turkey (14) and in Georgia were included in the model. Keim et al. (14) thought that the geographical distribution of the genotypes could be a function of isolate availability rather than the intrinsic diversity in a limited geographic area. Most of the Georgian isolates used in this study were obtained from soils that had a previous history of anthrax, but the anthrax outbreaks occurred long ago, so the collection should have been less biased than other collections that included only strains isolated during outbreaks. The sampling and isolation strategies used to obtain the Turkish isolates are unknown, making any estimate of the sampling bias infeasible.
![]() View larger version (17K): [in a new window] |
FIG. 1. Mutational step model showing the genetic relationships between the B. anthracis MLVA GK genotypes found in Turkey and Georgia and belonging to the A3.a genetic cluster, as defined by Keim et al. (14). Each circle represents a different genotype, and the notations on the lines connecting the genotypes indicate the minimal mutational changes needed to convert one genotype to another. The data are data for Georgian isolates characterized in this work and for the Turkish isolates characterized by Keim et al. (14). The dark gray circles indicate genotypes found in Turkey, the open circle indicates a genotype found only in Georgia, and the light gray circles indicate the genotypes found in both the countries.
|
Genotype GK 44 was found in both Turkey and Georgia, but considering the number of isolates with this genotype from the two countries, it seems to be a more prevalent genotype in Georgia (South Caucasus) than in Turkey. Thirteen of the 18 Georgian isolates had this genotype, and it was found at seven of the eight locations. The only location where it was not found was Gardabani, a city close to the border with Azerbaijan. Isolates from Gardabani had genotype GK 35, which is a prevalent Turkish genotype. The Gardabani region is a center for commercial exchange of livestock with Azerbaijan, and it is used as a place to pasture Azerbaijani livestock. The fact that genotype GK 35 was found only in Gardabani suggested that it could have been derived from adjacent Azerbaijan and indicates that a careful description of the pattern of B. anthracis spread in the Turkish-Southern Caucasian region requires typing of isolates from Azerbaijan and Armenia.
Although in the study of Keim et al. (14) the third Georgian genotype, genotype GK 40, was found only in African strains, this genotype differs from the other Georgian and Turkish genotypes only by the VNTR marker pXO2; hence, it is related to most of the other Turkish genotypes. The fact that genotype GK 40 was found in Georgia suggests that it could be typical of the Turkish-Southern Caucasian region, resulting from mutation of other genotypes in the region (e.g., genotype GK 35) rather than from independent introduction from other places.
The MLVA-based diversity data presented in this study establish the first picture of the genetic diversity of B. anthracis in the Southern Caucasus. A comparison of the data for the Georgian and Turkish isolates suggests that there is a common regional strain pattern in the Turkish South-Caucasian region, dominated by strictly related genotypes that differ only in a few polymorphisms in the pXO-based VNTR markers. Hence, it can be hypothesized that B. anthracis strains responsible for anthrax outbreaks in the region result from recent evolution of a common ancestral strain. However, the possibility that the pXO1 and pXO2 marker differences between the observed genotypes could be due to spontaneous mutations in these rapidly evolving markers due to strain subculturing cannot be excluded. Indeed, we did not intensively subculture the Georgian isolates to minimize this possible bias. Hence, to confirm the previous conclusion and to clarify the pattern of B. anthracis spread in the east-west and north-south directions in the Caucasian-Turkish region, it will be necessary to type by MLVA clinical and virulent isolates that occur in natural outbreaks in Azerbaijan and Armenia, as well as in the northern Caucasian region.
Four of the five vaccine strains lacked pXO2, and hence, the related plasmid marker could not be included in the analysis. Based on the other seven VNTR marker profiles (Table 1), these strains could be affiliated with several GK genotypes. Strains I-17 and 55 could be affiliated with genotypes in the A1.a cluster (genotype GK 15, GK 16, or GK 17 and genotype GK 3, GK 4, or GK 5, respectively). The strain STI-1 genotype confirmed the genotype reported by Keim et al. (14), with allele sizes (excluding the pXO2 marker) of 313, 229, 162, 613, 604, 153, and 129 bp (Table 1), and although it did not match any previously defined genotype due to the lack of the pXO2 plasmid, it could be affiliated with the A1.a cluster. Also, the Sterne 34F2 strain lacked the pXO2 marker and for the remaining seven markers had the MLVA genotype described by Keim et al. (14), with allele sizes of 313, 229, 162, 583, 532, 158, and 129 bp (Table 1). According to these markers the most similar genotypes were genotypes GK 59 and GK 61 in the A3.b cluster (14).
Strain 71/12 had an MLVA profile that does not match any of the previously described genotypes (Table 1). According to the 153-bp CG3 allele typical of cluster A1.a isolates (14), the strain 71/12 genotype is a novel genotype in the A1.a cluster, which can be distinguished from the other genotypes by means of a 532-bp vrrC2 marker.
In addition to the well-known Sterne 34F2 strain and the previously described STI-1 strain, there were three other strains which represented the Russian strains widely used in the former Soviet Union for vaccine preparation. Despite their origin all of these strains were grouped in the A1.a cluster with other isolates, most of which were from North America. Hence, their relatively high polymorphism is probably the vestigial result of a rather different geographical origin.
|
|
|---|
uzewski, A. Rakin, M. Jagielski, A. Zasada, A. Jakubczak, B. Borkowska-Opacka, and W. Rastawicki. 2004. Intriguing diversity of Bacillus anthracis in eastern Poland the molecular echoes of the past outbreaks. FEMS Microbiol. Lett. 239:235-240.[CrossRef][Medline]This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»