Research article | Open | Published:
MLVA distribution characteristics of Yersinia pestisin China and the correlation analysis
BMC Microbiologyvolume 9, Article number: 205 (2009)
Yersinia pestis, the aetiological agent of plague, has been well defined genotypically on local and worldwide scales. In November 2005, five cases of severe pneumonia of unknown causes, resulting in two deaths, were reported in Yulong, Yunnan province. In this study, we compared Y. pestis isolated from the Yulong focus to strains from other areas.
Two hundred and thirteen Y. pestis strains collected from different plague foci in China and a live attenuated vaccine strain of Y. pestis (EV76) were genotyped using multiple-locus variable-number tandem repeat analysis (MLVA) on 14 loci. A total of 214 Y. pestis strains were divided into 85 MLVA types, and Nei's genetic diversity indices of the various loci ranged between 0.02 - 0.76. Minimum spanning tree analysis showed that Y. pestis in China could be divided into six complexes. It was observed that Microtus strains were different from the other three biovar strains. Each plague focus had its own unique MLVA types.
The strains isolated from Yulong, Yunnan province had a unique MLVA type, indicating a new clone group. Our results suggest that Yulong strains may have a close relationship with strains from the Qinghai-Tibet Plateau plague focus.
Plague is an infectious disease caused by Yersinia pestis, a naturally occurring bacterium found primarily in wild rodents. It is highly transmissible and brings a high mortality, leading to major public health disasters throughout the history of humanity . In the early 1990s, the incidence of human plague increased significantly , with outbreaks occurring in Africa  and India . WHO has classified plague as a reemerging infectious disease for the past 20 years, and Y. pestis has been identified as a bioterrorism agent, posing as a significant threat to human health and safety . In November 2005, a natural focus of human plague was discovered in Yulong, Yunnan province, China. In this study, we compared Y. pestis isolated from the Yulong focus to strains from other areas.
Y. pestis couldn't be separated by serotype and phage-type, but could be classified into three biovars: Antiqua, Mediaevalis and Orientalis, according to their ability to ferment glycerol and to reduce nitrate as described by Devignat in the 1950s . Recently, a new biovar Microtus was proposed based on whole genome sequencing and genetic analysis [8, 9].
Y. pestis has a broad host and vector range . These hosts and vectors have their own natural environment, resulting in the diversity of micro-ecological environments for Y. pestis. During its expansion and adaption into new niches, Y. pestis undergoes considerable genome variability in response to natural selection. This variability can partly explain the genomic diversity of strains from different plague foci . At present, natural plague foci are widespread inChina. Through systematic analysis of Y. pestis in these areas, it is possible to understand the evolution of Y. pestis and investigate the source of new plague foci.
Previous studies have revealed a large number of tandem repeat sequences (TRSs) in the Y. pestis genome, and these TRSs introduce diversity into various plague strains . These loci are called variable-number tandem repeats (VNTRs). Multiple-locus VNTR analysis (MLVA) is an individual identification method that detects VNTR loci. MLVA is widely used in Y. pestis genotyping, and is useful for performing phylogenetic analysis [12–16]. In this study, 213 Y. pestis strains collected from different plague foci in China and a live attenuated vaccine strain of Y. pestis (EV76) were genotyped by MLVA using 14 loci.
Bacterial strains and DNA preparation
A total of 214 Y. pestis strains were included in this study. 208 strains were isolated from 13 natural plague foci in China between 1952 and 2002, an additional five strains were isolated from Yulong Yunnan in 2006, and the EV76 strain was also included in this study (Table 1). The bacteria were cultivated in Hottinger's medium at 28°C for 24 - 36 h, and then the genome DNAs were extracted by using conventional SDS lysis and phenol-chloroform extraction method. The bacterial culture and extraction of DNAs were performed in biosafety level 3 (BSL-3) laboratories.
VNTR locus selection
A total of 14 VNTR loci with core sequences >9 bp were selected from previously described VNTR loci [12, 17] (Table 2). The 14 VNTR loci had shown at least two alleles in six sequenced strains of Y. pestis (CO92, KIM, 91001, Nepal516, Antiqua, Angola). In order to provide an assay that is useful and widely accessible to research and public health laboratories, the present investigation favors markers with relatively large repeat units.
MLVA was performed as previously described . The loci and corresponding PCR primers were listed in Table 2. PCR reactions were prepared in a 20 μl volume with 10× PCR buffer, 0.5 U Taq polymerase, 200 μM each of the four dNTPs, 10 μM each primer set, 1 ng template DNA and filtered sterile water. PCR conditions were as follows: initial denaturation at 94°C for 5 min, and then 35 cycles of 94°C for 20 s, 60°C for 20 s and 72°C for 45 s, followed by a final polymerase extension step at 72°C for 5 min. Strain 91001 was used as a positive control. The PCR products were separated by electrophoresis with a 3% agarose gel (15 cm) in 0.5× Tris-borate-EDTA (TBE) buffer for about 11 h (4 v/cm), and the sizes of amplified DNA fragments were visualized by UV light. Copy numbers were calculated by size and assigned according to the number of repeats for each locus. The PCR products of samples, randomly selected from each number of repeats in each locus, were sequenced and the sequences were compared.
All data were input into BioNumerics version 4.0 software (Applied Maths, Belgium), cluster analysis and minimal spanning tree (MST) analysis was performed. Nei's indices  of VNTR loci and Simpson's diversity index  of MLVA typing method were calculated. The formulas were as follows: Nei's index = 1-Σ (allele frequency)2; Simpson's diversity index = 1 - Σ [n j (n j -1)]/[N (N-1)], where N refers to the total number of experimental strains and n j is the number of strains belonging to the jth type.
Genetic polymorphisms of Yersinia pestisin China
Based on analysis of 14 VNTR loci, 214 Y. pestis strains were divided into 85 MLVA types (MTs), named as MT01-MT85 (Figure 1). The Simpson's diversity index was 0.9790. The main characteristics of the 14 loci in 214 Y. pestis strains were shown in Table 3. Nei's diversity indices of the 14 loci were between 0.02 and 0.76. Locus M54 had higher Nei's diversity index than others. The numbers of alleles of the 14 loci were between 2 and 7. Loci M54, M58 and M61 had the largest number of alleles (n = 7). Locus M58 could further distinguish strains within 15 foci or sub-foci and the other 13 loci could further distinguish strains within 5 to 13 foci or sub-foci. Nei's diversity index of locus M72 was the smallest at 0.02, and it was just able to distinguish strains within two foci (Table 4).
According to cluster analysis based on 85 MTs (Figure 1), Microtus isolates comprising MT79 to MT85 (Group A) could be correctly distinguished from Antiqua, Orientalis and Medievalis. Within Microtus, the strains from two foci (L and M) could be further separated. Most Orientalis isolates were clustered into one major branch (Group B) with the exception of the strains representing MT09 and MT35. Similarly, 80% (20/25) of Medievalis isolates were also clustered together (Group C).
As a complementary analysis, a MST analysis was performed based on the categorical data sets (Figure 2). Six complexes and 3 single MTs were obtained. Complex 1, 4 and 5 represented Antiqua isolates and complex 2, 3 and 6 represented Orientalis, Medievalis and Microtus isolates, respectively. Complex 1 contained the largest number of strains (n = 130), which could be divided into 50 MTs. 84.35% (124/147) Antiqua strains were divided into complex 1. It was interesting that the strains isolated from the Xinjiang region (Figure 2, Foci A, B2, B3 and B4) constructed a long branch in complex 1. Complex 2 contained most of the Orientalis isolates, which were all isolated from Focus F (Figure 3). Complex 3 contained 18 Medievalis strains, which was account 72.00% (18/25) of all the Medievalis strains in this study, and three Antiqua strains. Complex 4 and complex 5 were constructed by Antiqua strains. Most of strains in complex 4 were from Focus G, while most of strains in complex 5 were from Focus H. All the Microtus isolates constituted complex 6, which was a well-defined complex representing Microtus isolates.
Strains from each focus presented their own unique MTs. For example, MT39 to MT43 were only found in Focus A, MT44 to MT51 were only found in Focus B, and MT17 was only found in Focus P. A total of 72 MTs were found in the specific foci (Figure 1). However, some strains isolated from different foci could share the same MTs. There were a total of 12 MTs (MT09, 18, 19, 21, 22, 26, 27, 35, 44, 52, 63, and 76) covering strains isolated from different foci. MT09 was shared by 10 strains isolated from 4 foci (C, D, J, F), including the main strains from Focus C. MT19 was shared by 10 isolates from 3 foci (D, C, K), including the main strains from Focus D. The other 10 MTs covered strains of 2 foci.
Most strains from the same focus presented the same or similar MTs (Figure 1). For example, the five strains in Focus P had exactly the same MT (MT17), and 6 of 9 bacteria isolated from Focus J had the same MT (MT53). The 11 strains of Focus A as a subgroup had very similar MTs, including seven strains for MT42 and other four strains for 4 MTs with one to three loci different from MT42. Most of the strains in Focus F were clustered together, including 14 strains for MT76 and the other six strains presenting in 6 MTs. On the other hand, strains from the same focus were dispersed in the cluster tree. For example, strains isolated from Focus G were dispersed in complex 1, 3 and 4, and strains from Focus C were scattered in complex 1 and 4.
MLVA comparison of Yersinia pestisin Yulong and the adjacent foci
Five strains isolated from Yulong, Yunnan had the same MT (MT17: 2-2-2-4-4-7-7-6-2-4-3-3-3-5). Three MTs with a difference in only one locus from MT17 were as follows: MT18 (2-2-2-4-4-7-7-7-2-4-3-3-3-5), including the strains from Foci C and G, had one copy difference on locus M58 with MT17; MT16 (2-2-2-4-4-7-7-6-2-4-3-2-3-5), including a strain which was isolated from Focus H, had one copy difference on locus M51 with MT17; MT29 (2-2-2-4-4-7-7-6-2-4-3-3-3-4), including a strain which was isolated from Focus C, had one copy difference on locus M37 with MT17.
The geographic locations of the natural plague foci adjacent to Yulong were C, E, and F (Figure 3). All the strains from Focus F were Orientalis, and the strains from Foci C, E and Yulong (Focus P) were Antiqua. A further MT comparisons between the Yulong strains and the strains isolated from Foci C and E were as follows: compared with Focus C, It was found that the five Yulong strains and five Focus C strains (belonging to MT29 to MT 33,) were clustered into group D (Figure 1); compared with Focus E, we found one copy difference located at three loci (M66, M58, and M54) in MT35 (major MT) and one copy difference located at four loci (M66, M58, M54, and M49) in MT23 (major MT); The MST analysis (Figure 2) showed that strains from Foci P, C, and E had a close relationship, and almost all strains belonged to one group.
In 2001, Klevytska et al. performed a systematic, whole genome analysis of Y. pestis CO92, and found that TRSs were widespread and randomly distributed in the bacterial chromosomes and plasmids . Subsequent studies had shown that MLVA could distinguish Y. pestis isolated from different natural plague foci [13–15, 20]. Our results showed that the loci selected in this study can distinguish the strains from different natural plague foci and even from the same focus. 214 Y. pestis strains used in this study were divided into 85 MTs. Simpson's diversity index was 0.9790, indicating that the probability of two unrelated strains being characterized as the same type was 2.10% (1 - 0.9790), showing high resolution and the combination of these 14 loci could be used as a typing method for Y. pestis with the generally accepted probability of 5% of type I errors .
However, a small number of strains from different foci had the same or similar MTs, suggesting that more VNTR loci should be tried to find out the combination of VNTR loci which could distinguish the strains among different foci completely, or that these strains from different foci may have the same source from another point of view. For example, we observed no considerable differences in the isolation times and places between the only human isolates (N010024, MT03) and the other strains isolated from Focus M. However, we did find a marked difference in MT. In previous studies, epidemiological investigations and traditional ecological typing studies confirmed that this case was imported from Focus C [22, 23]. In this study, N010024 was significantly different from the other strains isolated from Focus M, but had very similar MT with the strains from Focus C and gathered with them in the same subgroup. These results coincided with the conclusion of epidemiological investigations and the ecological typing, which further supported MLVA as a bacterial typing method suitable for field epidemiological investigations.
There were cross-types among the MTs of strains from different foci, with MT09 and MT19 being the most prominent. Foci that contained the same MT were geographically close to each other (Figure 3). For example, Foci C, D, F, and J contained MT09, and Foci C, D, and K contained MT19, indicating that there were close relationships among the strains of adjacent foci. It is possible that these strains have the same source. Foci C, D, G, and K have locations adjacent to the border and even similar topography, climate conditions and hosts. The Marmota Himalayana plague focus of the Qinghai-Tibet Plateau  was sub-divided into four foci in recent years . Cluster analysis showed that majority of the strains in the four foci were in complex 1, indicating a close relationship. Therefore, we suggest that more accurate results will be obtained by combining the four foci in a unit when performing epidemiological and phylogenetic analysis.
Foci A, B, and K are in Xinjiang province (Figure 1). The strains from Foci A and B were in the long branch of complex 1 and obviously different from other strains isolated in China. On the contrary, most strains from Focus K were together with the strains from foci around the Qinghai-Tibet Plateau. Foci A and B are adjacent to the Central Asia foci. Due to the lack of strains outside China in this study, it is impossible to provide a detailed and integrated relationship between the strains in Xinjiang and those of the Central Asia. However, we can confirm that there is a long genetic distance between strains of Foci A, B and other domestic strains isolated in China.
To date, all the strains from Foci L and M belonged to biovar Microtus, except for one imported strain (N010024). Microtus is a newly-identified biovar that is phenotypically and genotypically different from the other three biovars . Our results showed that MLVA could not only differentiate between Microtus and the other three biovars, but also divided the Microtus strains into two subclusters containing strains from foci L and M respectively. Five loci (M61, ms01, M58, M49 and M37) showed differences between the two subclusters. These loci could be used in the subsequent studies focused on Microtus strains.
In November 2005, five cases of severe pneumonia of unknown causes were reported in Yulong, Yunnan province, resulting in two deaths. These cases were subsequently diagnosed as plague, and the natural plague focus was confirmed by field investigation. Five strains of Y. pestis were isolated from host animals and vectors. Our results revealed that these five strains had exactly the same MT, suggesting that they had the same source. Furthermore, MT17 was different from the MTs of all the other strains, suggesting that the Yulong strains were a newly-discovered clone. In the 14 selected VNTR loci, M58 was a necessary locus which distinguishes the Yulong strains from the other strains. Moreover, it is also the marker with the second strongest discriminatory ability and the largest number of alleles. Consequently, we propose that M58 is a key locus for MLVA typing of Y. pestis in China.
The Yulong focus has distinct geographical features: it is adjacent to Focus E, and both of these foci are in the longitudinal valley area of Western Yunnan, located at the southeast edge of the Qinghai-Tibet Plateau. The two foci are also near Foci C and F. The cluster analysis and MST results suggested that the Yulong strains show a closer genetic relationship with the strains from Focus C than those from Focus E, as is consistent with the results of biological character comparisons  and insertion sequence typing . Therefore, it was predicted that the Yulong strains were more likely to be a new branch that evolved from Focus C, rather than the result of expansion and spread of Focus E.
The Yulong natural plague focus is adjacent to the previously-discovered Jianchuan focus (Focus E). Their natural conditions are the same, but the VNTR characteristics of strains from the two foci are critically different, suggesting that the two foci have relatively independent properties because of the hindrance of an ecological barrier. The recent occurrence of "severe pneumonia of unknown causes" in Yulong suggests that plague in this region is a threat to the human population. Since plague has not occurred in the Jianchuan region for a long time, the public health authorities in that area should remain vigilant in monitoring potential plague outbreaks.
MLVA is a very powerful and reproducible genotyping method and it is promising to be used as a simple molecular tool for characterization and epidemiological studies of Y. pestis. It could also unravel the molecular phylogeny of Y. pestis when being applied to a larger number of isolates. The 14 loci used in this study gave a high discriminatory power and successfully separated isolates of different biovars and from different natural plague foci. The strains isolated from a new focus, Yulong, Yunnan province had a unique MLVA type, indicating a new clone group. The results of cluster analysis and MST analysis suggest that the Yulong focus strains may have a close relationship with strains from the Qinghai-Tibet Plateau plague focus.
variable-number tandem repeat
multiple-locus VNTR analysis
tandem repeat sequences
minimum spanning tree.
Perry RD, Fetherston JD: Yersinia pestis-etiologic agent of plague. Clin Microbiol Rev. 1997, 10: 35-66.
Anonymous: Human plague in 1992. Wkly Epidemiol Rec. 1994, 69: 8-10.
Ratsitorahina M, Chanteau S, Rahalison L, Ratsifasoamanana L, Boisier P: Epidemiological and diagnostic aspects of the outbreak of pneumonic plague in Madagascar. Lancet. 2000, 355: 111-113. 10.1016/S0140-6736(99)05163-6.
Centers for Disease Control and Prevention (CDC): Update: human plague-India. MMWR Morb Mortal, Wkly Rep. 1994, 43: 761-762.
Broussard LA: Biological agents: weapons of warfare and bioterrorism. Mol Diagn. 2001, 6: 323-333.
Song ZZ, Xia LX, Liang Y, Guo Y, Lu L, Wang GL, Cai WF, Zhang ZF, He YT, Zhang FX, Dong XQ, Yu GL, Wang J, Yu DZ: Confirmation and study of Plague Natural Foci for Yulong County and Guchengqu in Yunnan Province. Chin J Ctrl Endem Dis. 2008, 23: 3-7.
Devignat R: Varieties of Pasteurella pestis; new hypothesis. Bull World Health Organ. 1951, 4: 247-253.
Song Y, Tong Z, Wang J, Wang L, Guo Z, Han Y, Zhang J, Pei D, Zhou D, Qin H, Pang X, Han Y, Zhai J, Li M, Cui B, Qi Z, Jin L, Dai R, Chen F, Li S, Ye C, Du Z, Lin W, Wang J, Yu J, Yang H, Wang J, Huang P, Yang R: Complete genome sequence of Yersinia pestis strain 9 an isolate avirulent to humans. DNA Res. 1001, 11: 179-197. 10.1093/dnares/11.3.179.
Zhou D, Tong Z, Song Y, Han Y, Pei D, Pang X, Zhai J, Li M, Cui B, Qi Z, Jin L, Dai R, Du Z, Wang J, Guo Z, Wang J, Huang P, Yang R: Genetics of Metabolic Variations between Yersinia pestis Biovars and the Proposal of a New Biovar, microtus. J Bacteriol. 2004, 186: 5147-5152. 10.1128/JB.186.15.5147-5152.2004.
Anisimov AP, Lindler LE, Pier GB: Intraspecific Diversity of Yersinia pestis . Clin Microbiol Rev. 2004, 17: 434-464. 10.1128/CMR.17.2.434-464.2004.
Li Y, Dai E, Cui Y, Li M, Zhang Y, Wu M, Zhou D, Guo Z, Dai X, Cui B, Qi Z, Wang Z, Wang H, Dong X, Song Z, Zhai J, Song Y, Yang R: Different Region Analysis for Genotyping Yersinia pestis Isolates from China. PLoS ONE. 2008, 3: e2166-10.1371/journal.pone.0002166.
Klevytska AM, Price LB, Schupp JM, Worsham PL, Wong J, Keim P: Identification and characterization of variable-number tandem repeats in the Yersinia pestis genome. J Clin Microbiol. 2001, 39: 3179-3185. 10.1128/JCM.39.9.3179-3185.2001.
Adair DM, Worsham PL, Hill KK, Klevytska AM, Jackson PJ, Friedlander AM, Keim P: Diversity in a variable-number tandem repeat from Yersinia pestis. J Clin Microbiol. 2000, 38: 1516-1519.
Pourcel C, André-Mazeaud F, Neubauer H, Ramisse F, Vergnaud G: Tandem repeats analysis for the high resolution phylogenetic analysis of Yersinia pestis. BMC Microbiol. 2004, 4: 22-10.1186/1471-2180-4-22.
Girard JM, Wagner DM, Vogler AJ, Keys C, Allender CJ, Drickamer LC, Keim P: Differential plague-transmission dynamics determine Yersinia pestis population genetic structure on local, regional, and global scales. Proc Natl Acad Sci USA. 2004, 101: 8408-8413. 10.1073/pnas.0401561101.
Lowell JL, Wagner DM, Atshabar B, Antolin MF, Vogler AJ, Keim P, Chu MC, Gage KL: Identifying Sources of Human Exposure to Plague. J Clin Microbiol. 2005, 43: 650-656. 10.1128/JCM.43.2.650-656.2005.
Le Flèche P, Hauck Y, Onteniente L, Prieur A, Denoeud F, Ramisse V, Sylvestre P, Benson G, Ramisse F, Vergnaud G: A tandem repeats database for bacterial genomes: application to the genotyping of Yersinia pestis and Bacillus anthracis. BMC Microbiol. 2001, 1: 2-10.1186/1471-2180-1-2.
Malorny B, Junker E, Helmuth R: Multi-locus variable-number tandem repeat analysis for outbreak studies of Salmonella enterica serotype Enteritidis. BMC Microbiol. 2008, 8: 84-10.1186/1471-2180-8-84.
Hunter PR, Gaston MA: Numerical index of the discriminatory ability of typing systems: an application of Simpson's index of diversity. J Clin Microbiol. 1988, 26: 2465-2466.
Huang XZ, Chu MC, Engelthaler DM, Lindler LE: Genotyping of a Homogeneous Group of Yersinia pestis Strains Isolated in the United States. J Clin Microbiol. 2002, 40: 1164-1173. 10.1128/JCM.40.4.1164-1173.2002.
Struelens MJ: Consensus guidelines for appropriate use and evaluation of microbial epidemiologic typing systems. Clin Microbiol Infect. 1996, 2: 2-11.
Hai R, Wei JC, Cai H, Zhou Y, Zhang JH: Molecular Biology Characters of Yersinia pestis Strains Isolated from Shiqu County, Sichuan Province, China. Chin J Vector Bio & Control. 2002, 13: 48-52.
Wei JC, Hai R, Sun LZ: Pulse - field gel electrophoresis of Yersinia pestis that isolated from Shiqu county, Sichuan province, China. Chinese Journal of Endemiology. 2002, 21: 309-311.
Ji SL, Zhang HB, Liu YP: The pattern of Yersinia pestis and its ecology significance in China. Zhongguo Yixue Kexue Yuan Xue Bao. 1983, 1-8. Special issue for Yersina pestis,
Zhu XY, Hai R, Song ZZ, Wei JC, Xia LX, Guo Y, Zhang HJ, Yu DZ: The primary genetic analysis of Yersinia pestis from Yulong with insertion sequence IS 285. Chin J Vector Bio & Control. 2008, 19: 144-147.
We gratefully thank Lijiang Center for Disease Control and Prevention, Yunnan, and Yunnan Institute for Endemic Disease Control and Prevention, China, for epidemiological investigation.
This work was supported by grant (200802016) from Ministry of Health of the People's Republic of China, grants (2004BA718B07 and 2008zx10004-008) from Ministry of Science and Technology of the People's Republic of China and Ministry of Health of the People's Republic of China.
ZXA did most of the typing work. WJC prepared the DNA samples. CZG was in charge of the Bionumerics database, clustering analyses and minimum spanning tree analyses. ZEM and SZZ were in charge of the epidemiological investigation and collection of Yulong strains. HR and YDZ initiated and managed the project, and ZXA and HR wrote the report. All authors read and approved the final manuscript.