|
|
||||||||

1 Division of Vector-Borne Infectious Diseases, Centers for Disease Control and Prevention, 3500 Rampart Rd, Fort Collins, CO 80522, USA
2 M. Aikimbayev's Kazakh Scientific Centre for Quarantine and Zoonotic Diseases, Almaty, Kazakhstan
3 Department of Biology, Colorado State University, Fort Collins, CO, USA
Correspondence
Jennifer L. Lowell
rzl9{at}cdc.gov
| ABSTRACT |
|---|
|
|
|---|
Present address: World Health Organization, 5120 Geneva Place, Dulles, VA 20189-5120, USA.
| INTRODUCTION |
|---|
|
|
|---|
Defining plague foci on the basis of phenotypic, biochemical and ecological properties of bacterial isolates has proven useful. Anisimov et al. (2004)
summarized the characteristics of Y. pestis isolates from central Asia, one of the hypothesized regions from which Y. pestis emerged as a flea-borne pathogen of rodents and other mammals (Achtman et al., 2004
). This summary described plague focus designations and geographical extent for each of the 43 purportedly distinct plague foci found in the southern and southeastern regions of the former Soviet Union (FSU), based on epizootic activity, main rodent hosts and main flea vectors. Subgroups of Y. pestis were also characterized using biochemical properties, rodent host and geographical region of isolation.
While these phenotypic properties have historically represented Y. pestis subtypes and helped define the foci from which they originated, relatively new molecular techniques have provided bacteriologists with tools capable of revealing genetic diversity in Y. pestis. Molecular typing methods such as multi-locus variable number tandem repeat analysis (MLVA) have helped to define Y. pestis population structure on several spatial scales (Achtman et al., 2004
; Girard et al., 2004
; Lowell et al., 2005
). On local scales, where isolates were collected within a few kilometres of each other, MLVA genotyping has demonstrated matches between environmental isolates and associated human plague infections in the southwestern United States (Lowell et al., 2005
), and genotype relationships of Y. pestis isolates among 19 prairie dog colonies in Arizona (Girard et al., 2004
). At broader scales, MLVA has inferred relationships between biovars based on clustering algorithms, but has not typically been used to show hierarchical relationships between the biovars (Achtman et al., 2004
; Klevytska et al., 2001
). Recent single-nucleotide polymorphism (SNP) analyses, however, have indicated that biovar Antiqua is polyphyletic (Achtman et al., 2004
; Chain et al., 2006
), and that the atypical Antiqua isolate Angola and the Y. pestis Pestoides' isolates are ancestral to all three biovars (Achtman et al., 2004
).
While phenotypic and molecular characterizations have provided valuable insight into worldwide Y. pestis collections and their relationships, studies examining collections of Y. pestis isolates from an intermediate scale across central Asia are lacking. Combining ecological, phenotypic and molecular information into single datasets could be especially useful for central-Asian isolates, as they potentially contain diversity not seen in parts of the world where Y. pestis is relatively new. The goals of this study were twofold. The first was to phenotypically characterize these isolates according to the tests commonly used by plague diagnostic and reference laboratories around the world (Chu, 2000
). The second was to examine intermediate geographical scale MLVA on Y. pestis isolates that have potentially circulated in nature for thousands of years, generating MLVA relationships different from those in North American (NA) isolates.
| METHODS |
|---|
|
|
|---|
|
|
Haemin-binding capabilities were determined by plating each isolate on Congo red agar and incubating for 48 h at 25 °C (Burrows & Jackson, 1956
; Surgalla & Beesley, 1969
). The presence or absence of the pgm locus was verified using irp2-directed PCR. Each 50 µl PCR reaction contained 1x PCR buffer with 1.5 mM MgCl2, 200 µM dNTPs, 0.5 U Taq polymerase (Promega), 0.5 ng Y. pestis DNA and 10 µM each of primers IRP2F (5'-TTTCCGGGGGCGTCTCAAC-3') and IRP2R (5'-CTCGACCTTTTGCAGTTTGATGTG-3'). Primers were designed using Primer Select (DNASTAR v 7.0) and the irp2 sequences of Y. pestis CO92 and KIM (GenBank accession nos NC_003143 and NC_004088, respectively). A negative control containing master mix only and a Y. pestis CO92 positive control were included in the set of reactions. All reactions were amplified in a PTC-200 thermal cycler (MJ Research) with an initial denaturation step of 95 °C for 5 min, followed by 35 cycles of 95 °C for 1 min, 55 °C for 1 min and 72 °C for 1 min, and a final extension step of 72 °C for 10 min. Products were 234 bp long and visualized by separating 5 µl reaction mixture on 2 % agarose gels containing 0.1 % ethidium bromide.
Sugar requirements and carbohydrate fermentation capabilities of all isolates were determined at KSCQZD according to the method of Aparin & Golubinskii (1989)
.
Phylogenetic analyses.
MLVA was used to compare isolates from Kazakhstan and from southwestern USA (NA isolates) on three geographical scales. A geographically distant scale was represented by comparing Orientalis, Medievalis, and Antiqua isolates from Kazakhstan to Orientalis isolates from NA. An intermediate geographical scale was represented by comparing the Kazakh isolates from 10 autonomous plague foci, while the NA Y. pestis set represented a close geographical scale analysis with isolates from human cases in Arizona, New Mexico and Colorado (Four Corners region), compared to the animal or flea isolates obtained during the associated epidemiological investigations. These isolates were previously described by Lowell et al. (2005)
, but for this analysis an additional 25 variable number of tandem repeat (VNTR) markers were added to the dataset to complete the set described by Klevytska et al. (2001)
, plus marker M19 from the work of Girard et al. (2004)
, for a total of 42 markers.
DNA was extracted and VNTR marker genotypes were determined according to the method of Lowell et al. (2005)
. Each 20 µl PCR multiplex reaction for the additional 25 markers contained 1x PCR buffer with 1.5 mM MgCl2, 200 µM dNTPs and 0.5 U Taq polymerase (Promega), 1.0 µl DNA template (approx. 0.5 ng DNA) and one of the following multiplex phosphoramidite linkage dye-labelled primer sets: 0.1 µM primer M26 and 0.25 µM M36; 0.1 µM M15 and 0.25 µM M37; 0.1 µM M42 and 0.2 µM M49; 0.1 µM M51 and 0.2 µM M52; 0.1 µM M43 and 0.25 µM M54; 0.1 µM M66 and 0.2 µM M55; 0.2 µM M68 and 0.25 µM M65; 0.1 µM M69 and 0.2 µM M56; 0.2 µM M61 and 0.25 µM M70; 0.1 µM each of M73 and M74; 0.1 µM each of M72 and M76; 0.1 µM each of M75 and M79; and 0.1 µM M71. Primer sets were run in pairs, as more consistent amplifications were achieved with fewer primers per multiplex. Pairs were chosen based on the primer dye label and on the size of the PCR fragment. All 42 markers were also analysed for the Kazakh isolates according to Lowell et al. (2005)
and the above primer combinations.
Of the 42 markers, three were variable and 34 were parsimony-informative. The phylogeny was inferred using equally weighted parsimony and 1000 tree-bisection-reconnection (TBR) searches with a maximum of 20 trees held per search. TBR branch swapping was then performed on all of the most parsimonious trees found with a maximum of 100 000 trees held, from which a strict consensus tree was calculated (Maddison, 1991
). Jackknife (JK) support (Farris et al., 1996
) was inferred using 1000 replicates, each consisting of 10 TBR searches and a maximum of 20 trees held (Fig. 2
). The atypical Antiqua Angola (http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=genomeprj&cmd=Retrieve&dopt=Overview&list_uids=16067) and Y. pestis strain 15-70 (Pestoides F) isolates (https://maple.lsd.ornl.gov/microbial/ypes_1570/) were chosen as outgroups (Farris et al., 1996
). Negative PCR results for Angola and Pestoides F were scored as missing data. All phylogenetic analyses were performed using PAUP 4.0b10 (Swofford, 2002
).
|
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
It has been established that the Y. pestis plasmids may spontaneously delete during laboratory storage (Perry & Fetherston, 1997
), but recent studies have shown that isolates undergoing several passages under laboratory conditions were highly variable in their plasmid stability and that one subculture lost both the pPMT1 and pPCP1 plasmids after only four passages, while others lost none after 32 passages (Leal-Balbino et al., 2004
). Atypical isolates lacking plasmids have also been recovered in nature from around the world (Cavalcanti et al., 2002
), including FSU foci (Filippov et al., 1990
). Because plasmid profiles were not performed on our Kazakh isolates at the time of collection, it is not possible to determine if plasmid loss occurred naturally or during storage. These isolates may have circulated in nature for thousands of years, undergoing high passage during epizootic activity, and isolates lacking pPCP1 or pMT1 may have persisted in mammals if the route of transmission was oral or pulmonary (Hinnebusch et al., 1998
, 2000
).
Thirty-two of the 48 Kazakh isolates (66 %) did not bind Congo red (Table 1
) and were also negative for irp2 PCR. Absence of the pgm locus including the yersiniabactin siderophore iron-acquisition system implies avirulence in these isolates (Hinnebusch et al., 1996
; Perry et al., 2001
). This result is most likely a laboratory artefact caused by spontaneous deletion because the isolates were passed many times during storage, or biased selection of pgm colonies during original bacterial isolation. This locus is known to spontaneously delete en bloc at a high rate (Perry & Fetherston, 1997
) and mixed populations of pgm+ and pgm isolates, with high frequency of pgm isolates, have been observed from the field (KSCQZD unpublished data; Cavalcanti et al., 2002
). If the majority of the bacterial colonies from the original post-epizootic isolations were pgm, selection of this phenotype for laboratory storage would be a likely scenario.
Nine leucine- and four tryptophan-dependent auxotrophic phenotypes were found in the isolate collection (Table 1
). Seven of the nine leucine-dependent auxotrophs were from the Pre-Ustyurt autonomous focus and one was from the Muyun-Kum autonomous focus, both located in the central-Asian natural region, where leucine dependence is typically seen. The remaining leucine-dependent isolate (3828) was collected from the Sarydzhaz autonomous focus, located in the Tien Shan natural focus, an area that does not typically yield this phenotype. Tryptophan-dependent auxotrophs are not usually identified in Kazakh isolate subtypes or autonomous plague foci; however, considering the variety of amino acid auxotrophs noted in the subtypes (Martinevsky, 1973
), finding these phenotypes in this collection was not a surprise. No additional amino acid auxotrophs were present.
Sugar fermentation characteristics for these 48 Kazakh isolates were typical for Y. pestis found in central-Asian plague foci. Isolate 3803 was phage resistant (Table 1
).
Kazakh and NA MLVA comparison
NA (Orientalis) and Kazakh (Medievalis) isolates formed distinct clades. Kazakh isolates 3831 and 3790 fell in the NA clade because they are biovar Orientalis. The NA clade was supported by 87 % of JK replicates while the remaining Kazakh isolates formed a monophyletic group with 50 % JK support (Fig. 2
). Markers M58, M59 and M52 were unique to biovar Orientalis.
Kazakh predefined foci versus Kazakh MLVA
The Kazakh phylogeny generally did not group isolates according to focus, with most isolates unresolved in a polytomy, although there were some exceptions (Fig. 2
). Two isolates from Volga-Ural sandy (3786, 3787) clustered with 87 % JK support, two from Pre-Ustyurt (3789, 3794) with 57 % support, two from Pre-Balkhash (3800, 3801) with 68 % support, two from the Volga-Ural steppe (3823, 3825) with 67 % support, and several isolates from three different foci [3802, 3804, 3805, (3820, 3821), 3822, 3824, 3832] formed a clade with 76 % support. Within this clade, two isolates from the Tau-Kum focus (3820, 3821) were supported in 87 % of JK replicates. An additional group of isolates from the North Pre-Aral (3808) and Pre-Aral Kara-Kum (3815) were supported with 91 % JK (Fig. 2
). These foci border one another, so it is possible that the isolates supported in this clade were collected in close geographical proximity.
Isolates that formed clades supported by at least 70 % JK differed by as many as four VNTR markers, and clades that were supported by between 50 and 70 % differed by as many as seven markers. Clades without support differed by at least eight markers.
This result is expected, given the local variation generated by relatively fast mutation rates in VNTR markers (Girard et al., 2004
) and by sampling on an intermediate geographical scale. Our sample set from Kazakhstan consisted of only a few isolates from each focus but the overall collection spanned most of the country. Isolates from Kazakh plague foci that have potentially existed for thousands of years may have undergone adaptation to local hosts, vectors or micro-environments, creating a high level of genetic diversity throughout the larger geographical area. Even the relatively young NA isolates that clustered in geographical space (collected less than
3 km apart) were strongly supported statistically (Girard et al., 2004
; Lowell et al., 2005
), providing evidence that individual Y. pestis clones causing plague epizootics tend to remain fairly isolated. While epizootics may eventually spread to adjacent rodent colonies after rapid animal die-offs like those observed within some rodent species, outbreaks typically remain isolated within a burrow system or rodent community for some time before the epizootic spreads to another location via dispersing juveniles or other mammals carrying infected fleas (Barnes, 1982
; Girard et al., 2004
). This pattern of isolated epizootics with little dispersal would maintain local bacterial populations, and supports the hypothesis that Y. pestis populations are maintained in relatively small areas rather than on spatial scales defined by plague foci spanning several thousand square kilometres (Fig. 1
). The low resolution of the Kazakh phylogeny using MLVA likely resulted from a combination of the geographical scale at which the Kazakh samples were collected, the relative age of Medievalis and Antiqua isolates, and the molecular marker mutation rates.
In contrast, MLVA distinguished between biovars and also provided information about their relative diversity. Biotypes Orientalis and Medievalis formed distinct clades (Fig. 2
). Antiqua did not form a monophyletic clade, indicating that these isolates were more diverse in their MLVA profiles. Orientalis showed the least amount of marker variability with 87 % JK support for the NA clade, Medievalis showed only 50 % support, and Antiqua was not supported. This result showed signs of genetic diversity that reflect previous analyses that described Antiqua as ancestral to Orientalis and Medievalis (Achtman et al., 1999
; Devignat, 1951
). Recent SNP analyses have demonstrated that one Antiqua isolate from Asia (Nepal 516) and several from Africa (Kenya) are polyphyletic and not ancestral to Asian Medievalis and worldwide Orientalis isolates, respectively (Achtman et al., 2004
; Chain et al., 2006
). However, Asian Antiqua isolates from Japan and the USSR are still shown ancestral to Asian Medievalis, and the atypical Antiqua isolate Angola is ancestral to all other biovars (Achtman et al., 2004
). When our MLVA tree was rooted with Angola and Y. pestis 15-70, the overall Y. pestis phylogeny was highly supported and one Kazakh Antiqua isolate (3828) formed a branch leading to the other Y. pestis biovars (Fig. 2
).
These isolates provided insight into the phenotypic and molecular diversity that exists in Kazakh Y. pestis. If these isolates naturally lacked plasmids important for blocked flea transmission, this may suggest that direct mammal-to-mammal transmission or mechanical transmission by fleas may be important in epizootic cycles, but evidence for this is sparse. If plasmids were lost because of storage practices, perhaps researchers should consider other possible mutations that may occur in vitro. It was possible that laboratory storage also affected VNTR marker diversity in the faster-mutating markers (Girard et al., 2004
). Although MLVA did not group these isolates by focus, we are confident that additional collections of fresh isolates on local scales would yield similar results to those seen in the USA. While current molecular typing and whole-genome sequencing of Y. pestis have provided invaluable insight into the evolution and diversity of Y. pestis, much of this work has been performed on old isolates that have probably been stored in a similar manner to this set. Future studies should include phenotypic and genotypic analyses conducted on fresh field isolates and stored isolates that have undergone minimal passages to differentiate between naturally occurring diversity and laboratory artefacts.
| ACKNOWLEDGEMENTS |
|---|
Edited by: S. MacIntyre
| REFERENCES |
|---|
|
|
|---|
Achtman, M., Morelli, G., Zhu, P. & 14 other authors (2004). Microevolution and history of the plague bacillus, Yersinia pestis. Proc Natl Acad Sci U S A 101, 1783717842.
Anisimov, A. P., Lindler, L. E. & Pier, G. B. (2004). Intraspecific diversity of Yersinia pestis. Clin Microbiol Rev 17, 434464.
Anonymous (1979). The instruction on epidemiological surveillance in the natural plague foci of Soviet Union. Proceedings of the Microbe Central Research Antiplague Institute, Saratov, 1979 (in Russian).
Aparin, G. P. & Golubinskii, E. P. (1989). Plague Microbiology Manual. Irkutsk, USSR: Irkutsk State University.
Barnes, A. M. (1982). Surveillance and control of bubonic plague in the United States. Symp Zool Soc Lond 50, 237270.
Begier, E. M., Asiki, G., Anywaine, Z., Yockey, B., Schriefer, M. E., Aleti, P., Ogden-Odoi, A., Staples, J. E., Sexton, C. & other authors (2006). Pneumonic plague cluster, Uganda, 2004. Emerg Infect Dis 12, 460467.[Medline]
Burrows, T. W. & Jackson, S. (1956). The pigmentation of Pasteurella pestis on a defined medium containing haemin. Br J Exp Pathol 37, 570576.[Medline]
Cavalcanti, Y. V., Leal, N. C. & De Almeida, A. M. (2002). Typing of Yersinia pestis isolates from the state of Ceara, Brazil. Lett Appl Microbiol 35, 543547.[CrossRef][Medline]
Cavanaugh, D. C. & Quan, S. F. (1953). Rapid identification of Pasteurella pestis using specific bacteriophage lyophilized on strips of filter paper; a preliminary report. Am J Clin Pathol 23, 619620.[Medline]
Chain, P. S., Carniel, E., Larimer, F. W., Lamardin, J., Stoutland, P. O., Regala, W. M., Georgescu, A. M., Verges, L. M., Land, M. L. & other authors (2004). Insights into the evolution of Yersinia pestis through whole-genome comparison with Yersinia pseudotuberculosis. Proc Natl Acad Sci U S A 101, 1382613831.
Chain, P. S., Hu, P., Malfatti, S. A., Radnedge, L., Larimer, F., Vergez, L. M., Worsham, P., Chu, M. C. & Andersen, G. L. (2006). Complete genome sequence of Yersinia pestis strains Antiqua and Nepal516: evidence of gene reduction in an emerging pathogen. J Bacteriol 188, 44534463.
Chu, M. (2000). Laboratory Manual of Plague Diagnostic Tests. CDC, Department of Health and Human Services.
Devignat, R. (1951). Varieties of Pasteurella pestis; new hypothesis. Bull World Health Organ 4, 247263.[Medline]
Farris, J. S., Albert, V. A., Kallersjo, M., Lipscomb, D. & Kluge, A. (1996). Parsimony jackknifing outperforms neighbor-joining. Cladistics 12, 99181.[CrossRef]
Filippov, A. A., Solodovnikov, N. S., Kookleva, L. M. & Protsenko, O. A. (1990). Plasmid content in Yersinia pestis strains of different origin. FEMS Microbiol Lett 55, 4548.[Medline]
Gage, K. L. & Kosoy, M. Y. (2005). Natural history of plague: perspectives from more than a century of research. Annu Rev Entomol 50, 505528.[CrossRef][Medline]
Girard, J. M., Wagner, D. M., Vogler, A. J., Keys, C., Allender, C. J., Drickamer, L. C. & Keim, P. (2004). Differential plague-transmission dynamics determine Yersinia pestis population genetic structure on local, regional, and global scales. Proc Natl Acad Sci U S A 101, 84088413.
Hinnebusch, B. J., Perry, R. D. & Schwan, T. G. (1996). Role of the Yersinia pestis hemin storage (hms) locus in the transmission of plague by fleas. Science 273, 367370.[Abstract]
Hinnebusch, B. J., Fischer, E. R. & Schwan, T. G. (1998). Evaluation of the role of the Yersinia pestis plasminogen activator and other plasmid-encoded factors in temperature-dependent blockage of the flea. J Infect Dis 178, 14061415.[CrossRef][Medline]
Hinnebusch, J., Cherepanov, P., Du, Y., Rudolph, A., Dixon, J. D., Schwan, T. & Forsberg, A. (2000). Murine toxin of Yersinia pestis shows phospholipase D activity but is not required for virulence in mice. Int J Med Microbiol 290, 483487.[Medline]
Huang, X. Z., Chu, M. C., Engelthaler, D. M. & Lindler, L. E. (2002). Genotyping of a homogeneous group of Yersinia pestis strains isolated in the United States. J Clin Microbiol 40, 11641173.
Kado, C. I. & Liu, S. T. (1981). Rapid procedure for detection and isolation of large and small plasmids. J Bacteriol 145, 13651373.
Klevytska, A. M., Price, L. B., Schupp, J. M., Worsham, P. L., Wong, J. & Keim, P. (2001). Identification and characterization of variable-number tandem repeats in the Yersinia pestis genome. J Clin Microbiol 39, 31793185.
Leal-Balbino, T. C., Leal, N. C., Lopes, C. V. & Almeida, A. M. (2004). Differences in the stability of the plasmids of Yersinia pestis cultures in vitro: impact on virulence. Mem Inst Oswaldo Cruz 99, 727732.[Medline]
Link, V. B. (1955). A History of Plague in the USA, pp. 12. Washington, DC: United States Department of Health, Education, and Welfare.
Lowell, J. L., Wagner, D. M., Atshabar, B., Antolin, M. F., Vogler, A. J., Keim, P., Chu, M. C. & Gage, K. L. (2005). Identifying sources of human exposure to plague. J Clin Microbiol 43, 650656.
Maddison, D. R. (1991). The discovery and importance of multiple islands of most-parsimonious trees. Syst Zool 40, 315328.[CrossRef]
Martinevsky, I. L. (1973). Materials for typing natural foci of plague according to the genetic properties of plague bacillus strains. J Hyg Epidemiol Microbiol Immunol 17, 272278.[Medline]
Motin, V. L., Georgescu, A. M., Elliott, J. M. & 8 other authors (2002). Genetic variability of Yersinia pestis isolates as predicted by PCR-based IS100 genotyping and analysis of structural genes encoding glycerol-3-phosphate dehydrogenase (glpD). J Bacteriol 184, 10191027.
Parkhill, J., Dougan, G., James, K. D. & 38 other authors (2001). Complete genome sequence of a multiple drug resistant Salmonella enterica serovar Typhi CT18. Nature 413, 848852.[CrossRef][Medline]
Perry, R. D. & Fetherston, J. D. (1997). Yersinia pestis etiologic agent of plague. Clin Microbiol Rev 10, 3566.[Abstract]
Perry, R. D., Bearden, S. W. & Fetherston, J. D. (2001). Iron and heme acquisition and storage systems of Yersinia pestis. Rec Res Dev Microbiol 5, 1327.
Protsenko, O. A., Filippov, A. A. & Kutyrev, V. V. (1991). Integration of the plasmid encoding the synthesis of capsular antigen and murine toxin into Yersinia pestis chromosome. Microb Pathog 11, 123128.[CrossRef][Medline]
Rahalison, L., Vololonirina, E., Ratsitorahina, M. & Chanteau, S. (2000). Diagnosis of bubonic plague by PCR in Madagascar under field conditions. J Clin Microbiol 38, 260263.
Surgalla, M. J. & Beesley, E. D. (1969). Congo red-agar plating medium for detecting pigmentation in Pasteurella pestis. Appl Microbiol 18, 834837.[Medline]
Swofford, D. L. (2002). PAUP* Phylogenetic Analysis Using Parsimony* and Other Methods. Sunderland, MA: Sinauer Associates.
Received 11 April 2006;
revised 1 July 2006;
accepted 19 September 2006.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |