Genetic diversity of Mycobacterium tuberculosis isolates from Beijing, China assessed by Spoligotyping, LSPs and VNTR profiles
BMC Infectious Diseases volume 12, Article number: 372 (2012)
Tuberculosis is one of the most infectious diseases in the world. Molecular typing methods such as spoligotyping, and VNTR (variable number tandem repeats), IS6110 in the NTF region and LSP (large sequence polymorphisms) analysis are generally useful tools for the resolution of various issues related to the classical epidemiology of Mycobacterium tuberculosis (M. tuberculosis).
To determine the transmission characteristics of M. tuberculosis strains isolated in Beijing, China, and their genetic relationships, especially those among Beijing family strains, 260 M. tuberculosis strains isolated from patients presenting pulmonary tuberculosis were analyzed by spoligotyping, and by examining 22 VNTR loci and the presence/absence of IS6110 in the NTF region, RD105 and RD181.
81% (211 strains) of the isolates studied were Beijing family strains, 174 (82.5%) of which were identified as modern Beijing strains based on the presence of IS6110 upstream of the NTF region. RD181 was intact in 9 of the other 37 (17.5%) ancestral Beijing strains. The percentage of Beijing family strains in this study was consistent with previous reports. There are many differences, however, in allele diversity among VNTR loci between reports on strains from different areas.
The Beijing family is the most prevalent genotype in Beijing city and the predominance of Beijing family strains has not altered in almost twenty years. Differences in the alleles and discrimination ability of VNTR loci between different regions is likely due to population differences in the regions where these M. tuberculosis strains were isolated or to differences in sampling times.
Tuberculosis (TB) remains one of the most serious infectious diseases in the world. The Mycobacterium tuberculosis (M. tuberculosis) lineage of greatest significance in China is the Beijing family. This lineage, which was first described in 1995 , has been reported to cause major outbreaks worldwide [2, 3].
M. tuberculosis Beijing family strains can be recognized on the basis of highly conserved spoligotyping (spacer oligonucletide typing) patterns (spacers 1–34 are absent) [4–6], and characteristic IS6110-RFLP patterns containing a high number of bands . Strains lacking one or more of the last nine spacers in the spoligotyping pattern are called Beijing-like strains .
Beijing family strains have been divided into two major groups, modern and ancient, based on a specific IS6110 in the NTF region [6, 9]. The group that has no IS6110 in the NTF region is thought to be the most ancient .
As with other PCR-based genotyping techniques, VNTR (variable number tandem repeats) typing is a promising, quick and easy genotyping method and can be used for epidemiological studies of Beijing strains . This method is based on size analysis of PCR-amplified VNTR loci, and it requires only basic PCR and agarose electrophoresis equipment . The number of repeat copies per locus may vary among strains, and the use of several such loci has sufficient discriminating power . VNTR profiles are presented as multiple digit numerical codes, each digit representing the copy number of a locus. However these VNTR loci present multiple independent genetic markers and are therefore ideally suited for evolutionary analyses .
LSPs (Large sequence polymorphisms), described in numerous studies of the M. tuberculosis complex, are reported to give important insights into its evolution and biology . One LSP (RD105) has been observed in all Beijing family strains and ancestral spoligotyping pattern strains . Additional LSPs (RD181) further divide the Beijing family into subgroups .
The aims of the study were to analyze the distribution of the tuberculosis lineage in Beijing city and to compare our results with those from other areas in China. In this study, we collected 260 M. tuberculosis clinical isolates from Beijing city in 2009, and genotyped them by spoligotyping, VNTR, IS6110 in the NTF and LSP (RD105 and RD180) to provide some information of potential transmission and evolution.
M. tuberculosisclinical isolates
260 M. tuberculosis strains were collected from TB patients in Beijing, China in 2009. The patients were all diagnosed to have pulmonary tuberculosis according to the national guidelines of China, and were treated in local lung disease hospitals or tuberculosis hospitals. There was no direct transmission link between patients selected in this study. All patient samples were taken as part of standard care. Demographic data, including the identification number, sex, age, occupation, registered residence, current address, as well as results of mycobacterial smears, symptoms and treatment history were obtained from medical records provided by Chaoyang Centre for Disease Control and Prevention, Beijing.
The study was approved by the Ethics Committee of the National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention. All patients in the study signed informed consent forms.
Culture of clinical isolates and DNA extraction
M. tuberculosis genomic DNA was extracted from mycobacterial colonies grown on Löwenstein-Jensen (LJ) medium by resuspending one loopful of mycobacterial colonies in 100 μl TE buffer (10 mM Tris-Cl, 1 mMEDTA) which was then incubated at 95°C for 15 min. Supernatants containing the DNA were collected by centrifugation at 12,000 g for 3 min and stored at −20°C for further use.
Spoligotyping which relies on the highly polymorphic DR (direct repeat) locus in the genome of M. tuberculosis, is a novel method for defining the Beijing family. Spoligotyping was performed according to the method described by Kamerbeek et al. .
22 VNTR loci, including 5 exact tandem repeat (ETR A,B,C,D and E) loci, 7 mycobacterial interspersed repetitive-units (MIRU10,16,23,26,27,39 and 40), 4 Mtub loci (Mtub4, 21, 30 and 39), 5 Qub (Queens’s University of Belfast) loci (Qub11a, 11b, 18, 26, and 4156), and VNTR3820 locus, were analyzed in this study. PCR products were analyzed on 12 cm 2% agarose gels with 11 lanes. Each gel included two 100-bp DNA Markers (one on either side of the gel), an H37Rv sample as a positive control, and 8 samples. The copy number at each locus was calculated by BioNumerics software. Determination of the discriminatory power of the VNTR loci was calculated using the Hunter-Gaston discriminatory Index (HGDI) . The HGDI was calculated using the following formula: , where HGDI is the numerical index of discrimination, N is the total number of strains in the typing scheme, s is the total number of different strain types, and nj is the number of strains belonging to the jth type [17, 18].
IS6110in the NTF region
The IS6110 in the NTF region was analyzed according to methods used in previous studies . All Beijing family strains identified by the spoligotyping method were amplified by PCR to detect the presence or absence of IS6110 in the NTF region. PCR analysis was performed using the following oligonucleotide primers: MDR-6, 5′-CCAGATATCGGGTGTGTCGAC-3′; MDR-6r, 5′-TGCCGTTGTCGAAATCTAAACCC-3′. Strains with the insert yielded an amplified product of ≈1,800-bp, while those without the insert yielded an ≈700-bp PCR product [19, 20].
Large sequence polymorphisms (LSPs)
RD105, and RD181 were analyzed in all the Beijing family strains from our study population, using a method described previously . PCR analysis was performed using the following oligonucleotide primers: RD105L, 5′-GGAGTCGTTGAGGGTGTTCATCAGCTCAGTC-3′; RD105R, 5′-CGCCAAGGCCGCATAGTCACGGTCG-3′; RD181L, 5′-CGCAACGGCCGCGGTGAACTCT-3′; RD181R, 5′-CGGGCGGCTGCGGGAACCTT-3′.
Spoligotypes in binary format were entered in an Excel spreadsheet and compared with the spoligotyping database SpolDB4 (http://www.pasteur-guadeloupe.fr:8081/SITVITDemo/index.jsp) . BioNumerics software version 5.0 (Applied Maths, Belgium) was used for data analysis and generating a cluster map. Clustering was based on VNTR results, and was performed using the categorical co-efficient and UPGMA in BioNumerics 5.0.
Of the 260 isolates, 246 (94.62%) were classified into one of the 17 shared international types (SITs) using SpolDB4.0 (Table 1 and Additional file 1: Table S1). The remaining 14 isolates generated 12 different spoligotypes that have not previously been described in the database. The most frequent strains were SIT1 (77.3%), followed by SIT53 (6.92%) and SIT190 (1.92%). 211 strains (81%) belonged to either the Beijing family (including SIT1, SIT190, SIT1364 and two new types), or the Beijing-like family (including SIT269, SIT632, and SIT255).
Allelic diversity of the VNTR loci
We analyzed all 260 isolates by screening 22 VNTR loci to investigate the genotypic diversity of the population in detail. This genotyping method is used as a standard discrimination tool for M. tuberculosis because of its high powers of resolution among isolates from cosmopolitan origins . VNTR results are shown in Additional file 1: Table S1 and VNTR profiles analyzed in this study are listed in Table 2. The cluster map in Figure 1 of the 260 M. tuberculosis strains was generated by BioNumerics 5.0. The 260 strains were classified as G1 to G4 strains using the cluster cutoff value calculated by the BioNumerics software according to VNTR results. G1 strains were further divided into three subgroups, G1-1, G1-2, and G1-3 (Figure 1). VNTR loci differed significantly in their allelic diversity and power of discrimination (Table 2). The VNTR3820 locus showed the highest allelic diversity among the 260 strains (0.851), while the ETRC locus showed the lowest allelic diversity (0.03).
The VNTR profiles of Beijing family strains analyzed in this and previous studies were combined and are listed in Table 2. The VNTR3820 locus showed the highest allelic diversity among the 211 Beijing family strains (0.789), while the ETRC locus showed the lowest allelic diversity (0.028). Comparing the allelic diversity observed here with that of previous studies carried out in regions such as Jiangsu province or Shanghai city, both of which are located in the south of China [22, 23], we found more differences in the allelic diversity of some loci (for example, ETRE and MIRU26). Furthermore, in a previous study on Beijing city, published by Jiao et al. four years ago , the alleles and discrimination power of MIRU16 and MIRU23 were very different to that observed in this study (Table 2). For instance, we obtained an HGDI for MIRU16 of 0.265, while Jiao et al., obtained a HGDI of 0.068; we observed five MIRU23 alleles (2, 3, 4, 5, 6), while only 2 alleles (5, 6) were observed by Jiao et al.
NTF and LSP
The RD105 deletion was found in 211 Beijing family strains. Of the 211 epidemiologically unlinked Beijing MTB isolates, 174 (82.5%) were identified as modern Beijing strains based on the presence of IS6110 upstream of the NTF region. In the remaining 37 (17.5%) ancestral Beijing strains, 9 strains had an intact RD181.
The definition of the phylogenetic sublineages of the 260 M. tuberculosis isolates from Beijing is listed in Table 3. G1 strains are Beijing family strains and have an RD105 deletion. G1-1 subgroup strains have an intact RD 181. The G1-1 and G1-2 subgroups are ancient Beijing family strains (Table 3).
Molecular typing methods are generally useful tools for resolving various issues related to the classical epidemiology of human pathogens, including M. tuberculosis. Spoligotyping is definitely a “gold standard” method for detecting a Beijing family strain of M. tuberculosis. The Beijing genotype was identified for the first time in strains isolated in the Beijing area of China in 1995, giving rise to the name of the genotype .
Beijing family strains having a characteristic spoligotype pattern (loci 1–34 are absent; loci 35–43 are present) have been identified in many areas. Moreover, in an M. tuberculosis molecular epidemiological investigation in China, Beijing family strains were the most prevalent lineage in China (74.08%), based on the SpolDB4.0 spoligotype database and SpotClust results . In this study 211 strains were Beijing family strains, accounting for 81% of all isolates examined.
Interestingly, in the first molecular epidemiology study in Beijing city, 89.4% of the M. tuberculosis strains collected from 1992 to 1994 were identified as Beijing family strains , and 85.1% of strains collected in Beijing from 2002 to 2005 were also identified as Beijing family strains . Is this study of isolates collected in 2009, 81% of the strains were Beijing family strains. Since there is no statistically significant difference between these three percentages, as determined using the X 2 test (p > 0.05) we conclude that the predominance of Beijing family strains has not altered for almost twenty years.
There is some debate over the use of RD181 for classifying isolates. An early study showed that RD181 can be used to divide Beijing isolates into ancestral (RD181 intact) and modern (RD181 deleted) strains . Faksri et al., however, suggested that RD181 alone does not provide sufficient discrimination to define ancestral or modern Beijing lineages, but with RD181 bing intact or deleted ancestral Beijing family strains can be divided into two subgroups . Results from this study are consistent with this report in that ancestral Beijing family strains could be divided into RD181 intact strains (9 strains) and RD181 deletion strains (28 strains).
During the process of evolution in Beijing family strains, numbers of VNTR repeats can increase or decrease. Changes in the numbers of repeats in VNTRs can differ with locus. It has been shown that expansion and contraction in the number of repeats of VNTRs can occur during the evolutionary process . Many previous studies have shown that VNTR loci can be good markers for phylogenetic estimation . In this study, we found that the number of Mtub21 repeats was different between ancient and modern Beijing family strains, with the number of repeats being much lower in non-Beijing family strains. In agreement with other reports, we conclude that Mtub21 is an evolutionarily informative VNTR locus .
Many previous studies have suggested that VNTR loci vary in their ability to discriminate Beijing genotype strains from geographically distant areas [11, 22–24] (Table 2). Some highly polymorphic VNTR loci have different powers of discrimination among strains from different regions; for example the discriminatory power of ETRE among strains from Jiangsu  and MIRU26 among strains from Shanghai  is much greater than among strains from Beijing. We compared the discrimination power of VNTR loci in this study with those in a previous study on Beijing strains published four years ago . Jiao et al. typed 72 Beijing family strains isolated in Beijing city from 2002 to 2005 using 24 VNTR loci. It is interesting that the discrimination power and alleles of some loci were different from our results. The differences are likely due to the different populations in the distinct geographic areas where the M. tuberculosis strains were isolated or to differences in sampling times.
This study suffers from one main limitation. Our study based on a comparison of strains from Beijing, Shanghai and Jiangsu provinces suggests that VNTR results for M. tuberculosis strains differ according to geographic region. However, more data from other provinces is needed to precisely identify the nature of differences between different regions of China.
Beijing family strains are the most prevalent genotypes in Beijing city. The predominance of Beijing family strains has not altered for almost twenty years. Many of the differences in allelic diversity between VNTR locus studies from different areas are likely due to the population differences in the regions from which the M. tuberculosis strains were isolated or to different sampling times.
van Soolingen D, Qian L, de Haas PE, Douglas JT, Traore H, Portaels F, Qing HZ, Enkhsaikan D, Nymadawa P, van Embden JD: Predominance of a single genotype of Mycobacterium tuberculosis in countries of east Asia. J Clin Microbiol. 1995, 33 (12): 3234-3238.
Brosch R, Gordon SV, Marmiesse M, Brodin P, Buchrieser C, Eiglmeier K, Garnier T, Gutierrez C, Hewinson G, Kremer K, et al: A new evolutionary scenario for the Mycobacterium tuberculosis complex. Proc Natl Acad Sci USA. 2002, 99 (6): 3684-3689. 10.1073/pnas.052548299.
Sola C, Filliol I, Legrand E, Mokrousov I, Rastogi N: Mycobacterium tuberculosis phylogeny reconstruction based on combined numerical analysis with IS1081, IS6110, VNTR, and DR-based spoligotyping suggests the existence of two new phylogeographical clades. J Mol Evol. 2001, 53 (6): 680-689. 10.1007/s002390010255.
Parwati I, van Crevel R, van Soolingen D: Possible underlying mechanisms for successful emergence of the Mycobacterium tuberculosis Beijing genotype strains. Lancet Infect Dis. 2010, 10 (2): 103-111. 10.1016/S1473-3099(09)70330-5.
Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, Kuijper S, Bunschoten A, Molhuizen H, Shaw R, Goyal M, et al: Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. J Clin Microbiol. 1997, 35 (4): 907-914.
Kremer K, Glynn JR, Lillebaek T, Niemann S, Kurepina NE, Kreiswirth BN, Bifani PJ, van Soolingen D: Definition of the Beijing/W lineage of Mycobacterium tuberculosis on the basis of genetic markers. J Clin Microbiol. 2004, 42 (9): 4040-4049. 10.1128/JCM.42.9.4040-4049.2004.
van Embden JD, Cave MD, Crawford JT, Dale JW, Eisenach KD, Gicquel B, Hermans P, Martin C, McAdam R, Shinnick TM, et al: Strain identification of Mycobacterium tuberculosis by DNA fingerprinting: recommendations for a standardized methodology. J Clin Microbiol. 1993, 31 (2): 406-409.
Chan MY, Borgdorff M, Yip CW, de Haas PE, Wong WS, Kam KM, Van Soolingen D: Seventy percent of the Mycobacterium tuberculosis isolates in Hong Kong represent the Beijing genotype. Epidemiol Infect. 2001, 127 (1): 169-171.
Kremer K, Van-der-Werf MJ, Au BK, Anh DD, Kam KM, Van-Doorn HR, Borgdorff MW, Van-Soolingen D: Vaccine-induced immunity circumvented by typical Mycobacterium tuberculosis Beijing strains. Emerg Infect Dis. 2009, 15 (2): 335-339. 10.3201/eid1502.080795.
Mokrousov I, Ly HM, Otten T, Lan NN, Vyshnevskyi B, Hoffner S, Narvskaya O: Origin and primary dispersal of the Mycobacterium tuberculosis Beijing genotype: clues from human phylogeography. Genome Res. 2005, 15 (10): 1357-1364. 10.1101/gr.3840605.
Kremer K, Au BK, Yip PC, Skuce R, Supply P, Kam KM, van Soolingen D: Use of variable-number tandem-repeat typing to differentiate Mycobacterium tuberculosis Beijing family isolates from Hong Kong and comparison with IS6110 restriction fragment length polymorphism typing and spoligotyping. J Clin Microbiol. 2005, 43 (1): 314-320. 10.1128/JCM.43.1.314-320.2005.
Supply P, Mazars E, Lesjean S, Vincent V, Gicquel B, Locht C: Variable human minisatellite-like regions in the Mycobacterium tuberculosis genome. Mol Microbiol. 2000, 36 (3): 762-771.
Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rusch-Gerdes S, Willery E, Savine E, de Haas P, van Deutekom H, Roring S, et al: Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J Clin Microbiol. 2006, 44 (12): 4498-4510. 10.1128/JCM.01392-06.
Flores L, Van T, Narayanan S, DeRiemer K, Kato-Maeda M, Gagneux S: Large sequence polymorphisms classify Mycobacterium tuberculosis strains with ancestral spoligotyping patterns. J Clin Microbiol. 2007, 45 (10): 3393-3395. 10.1128/JCM.00828-07.
Tsolaki AG, Gagneux S, Pym AS, Goguet de la Salmoniere YO, Kreiswirth BN, Van Soolingen D, Small PM: Genomic deletions classify the Beijing/W strains as a distinct genetic lineage of Mycobacterium tuberculosis. J Clin Microbiol. 2005, 43 (7): 3185-3191. 10.1128/JCM.43.7.3185-3191.2005.
Hunter PR: Reproducibility and indices of discriminatory power of microbial typing methods. J Clin Microbiol. 1990, 28 (9): 1903-1905.
Gaston MA, Hunter PR: Efficient selection of tests for bacteriological typing schemes. J Clin Pathol. 1989, 42 (7): 763-766. 10.1136/jcp.42.7.763.
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 (11): 2465-2466.
Wada T, Iwamoto T, Maeda S: Genetic diversity of the Mycobacterium tuberculosis Beijing family in East Asia revealed through refined population structure analysis. FEMS Microbiol Lett. 2009, 291 (1): 35-43. 10.1111/j.1574-6968.2008.01431.x.
Plikaytis BB, Marden JL, Crawford JT, Woodley CL, Butler WR, Shinnick TM: Multiplex PCR assay specific for the multidrug-resistant strain W of Mycobacterium tuberculosis. J Clin Microbiol. 1994, 32 (6): 1542-1546.
Brudey K, Driscoll JR, Rigouts L, Prodinger WM, Gori A, Al-Hajoj SA, Allix C, Aristimuno L, Arora J, Baumanis V, et al: Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology. BMC Microbiol. 2006, 6: 23-10.1186/1471-2180-6-23.
Zhang L, Chen J, Shen X, Gui X, Mei J, Deriemer K, Gao Q: Highly polymorphic variable-number tandem repeats loci for differentiating Beijing genotype strains of Mycobacterium tuberculosis in Shanghai, China. FEMS Microbiol Lett. 2008, 282 (1): 22-31. 10.1111/j.1574-6968.2008.01081.x.
Liu Q, Yang D, Xu W, Wang J, Lv B, Shao Y, Song H, Li G, Dong H, Wan K, et al: Molecular typing of Mycobacterium tuberculosis isolates circulating in Jiangsu province, China. BMC Infect Dis. 2011, 11: 288-10.1186/1471-2334-11-288.
Jiao WW, Mokrousov I, Sun GZ, Guo YJ, Vyazovaya A, Narvskaya O, Shen AD: Evaluation of new variable-number tandem-repeat systems for typing Mycobacterium tuberculosis with Beijing genotype isolates from Beijing, China. J Clin Microbiol. 2008, 46 (3): 1045-1049. 10.1128/JCM.01869-07.
Dong H, Liu Z, Lv B, Zhang Y, Liu J, Zhao X, Wan K: Spoligotypes of Mycobacterium tuberculosis from different Provinces of China. J Clin Microbiol. 2010, 48 (11): 4102-4106. 10.1128/JCM.00549-10.
Dou HY, Tseng FC, Lin CW, Chang JR, Sun JR, Tsai WS, Lee SY, Su IJ, Lu JJ: Molecular epidemiology and evolutionary genetics of Mycobacterium tuberculosis in Taipei. BMC Infect Dis. 2008, 8: 170-10.1186/1471-2334-8-170.
Faksri K, Drobniewski F, Nikolayevskyy V, Brown T, Prammananan T, Palittapongarnpim P, Prayoonwiwat N, Chaiprasert A: Genetic diversity of the Mycobacterium tuberculosis Beijing family based on IS6110, SNP, LSP and VNTR profiles from Thailand. Infect Genet Evol. 2011, 11 (5): 1142-1149. 10.1016/j.meegid.2011.04.007.
Warren RM, Victor TC, Streicher EM, Richardson M, van der Spuy GD, Johnson R, Chihota VN, Locht C, Supply P, van Helden PD: Clonal expansion of a globally disseminated lineage of Mycobacterium tuberculosis with low IS6110 copy numbers. J Clin Microbiol. 2004, 42 (12): 5774-5782. 10.1128/JCM.42.12.5774-5782.2004.
Wada T, Iwamoto T: Allelic diversity of variable number of tandem repeats provides phylogenetic clues regarding the Mycobacterium tuberculosis Beijing family. Infect Genet Evol. 2009, 9 (5): 921-926. 10.1016/j.meegid.2009.06.012.
Millet J, Miyagi-Shiohira C, Yamane N, Mokrousov I, Rastogi N: High-resolution MIRU-VNTRs typing reveals the unique nature of Mycobacterium tuberculosis Beijing genotype in Okinawa, Japan. Infect Genet Evol. 2012, 12 (4): 637-641. 10.1016/j.meegid.2011.06.013.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/12/372/prepub
This research was supported by the State Key Laboratory for Infectious Diseases Prevention and Control (2011SKLID208), and the National Key Program of Mega Infectious Diseases (2008ZX100/03-010).
The authors declare that they have no competing interests.
BL participated in data analysis and drafted the manuscript; PZ and XZ participated in sample collection; BL and QY carried out the molecular genetic studies; HD and KW participated in the design of the study. All authors read and approved the final manuscript.
Bing Lu, Ping Zhao, Binbin Liu, Haiyan Dong contributed equally to this work.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Lu, B., Zhao, P., Liu, B. et al. Genetic diversity of Mycobacterium tuberculosis isolates from Beijing, China assessed by Spoligotyping, LSPs and VNTR profiles. BMC Infect Dis 12, 372 (2012). https://doi.org/10.1186/1471-2334-12-372
- M. tuberculosis
- Beijing family