Drug resistance and IS6110‐RFLP patterns of Mycobacterium tuberculosis in patients with recurrent tuberculosis in northern Thailand

The emergence of drug resistant Mycobacterium tuberculosis has become a global threat to tuberculosis (TB) prevention and control efforts. This study aimed to determine the drug resistance profiles and DNA fingerprints of M. tuberculosis strains isolated from patients with relapsed or retreatment pulmonary TB in Chiang Rai province in northern Thailand. Significant differences in multidrug resistance (MDR) (P = 0.025) and resistance to isoniazid (P = 0.025) and rifampin (P = 0.046) between first and second registrations of patients with retreatment TB were found. However, there were no significant differences in resistance to any drugs in patients with relapsed TB. The rate of MDR‐TB strains was 12.2% among new patients at first registration, 22.5% among patients with recurrence who had previously undergone treatment at second registration and 12.5% at third registration. Two retreatment patients whose initial treatment had failed had developed MDR‐TB with resistance to all TB drugs tested, including rifampin, isoniazid, streptomycin and ethambutol. IS6110‐RFLP analysis revealed that 66.7% (10/15 isolates) of MDR‐TB belonged to the Beijing family. In most cases, IS6110‐RFLP patterns of isolates from the same patients were identical in relapse and retreatment groups. However, some pairs of isolates from retreatment patients after treatment failure had non‐identical IS6110‐RFLP patterns. These results suggest that, after failure and default treatment, patients with retreatment tuberculosis have a significantly greater risk of MDR‐TB, isoniazid and rifampin resistance than do other patients.     

Studies of Polymorphic DNA Fingerprinting and Lipid Pattern of Mycobacterium tuberculosis Patient Isolates in Japan

Strain differentiation by DNA restriction fragment length polymorphism (RFLP) has been used mainly for the epidemiological purpose of Mycobacterium tuberculosis infection. In this study, we tried to connect the molecular and phenotypic characteristics of M. tuberculosis patient isolates by comparing the DNA fingerprints obtained by RFLP using IS6110 and lipid patterns using two-dimensional thin-layer chromatography (2-D TLC) with silica gel, since M. tuberculosis has a lipid-rich cell envelope which contributes to the virulence and immunomodulatory properties. We found that 66 isolates of M. tuberculosis from tuberculosis patients showed that the occurrence of IS6110 varied from 1 to 24 copies. The IS6110 patterns were highly variable among isolates. Fifty different RFLP patterns were observed, and 12 RFLP patterns were shared by two or more strains. By computerized analysis of the RFLP patterns of M. tuberculosis patient isolates, we found that 95% of the isolates fell into seven clusters, from A to G, with at least two isolates in each (> 30% similarity). Among the cellular lipids, the phospholipid composition did not differ by strain, whereas the glycolipid pattern differed markedly. Especially, the relative concentration of cord factor and sulfolipid, both of which were known as virulent factors, varied by strain. The fingerprints of some strains showed an association between the DNA and glycolipid patterns, even though some of the same DNA fingerprint strains showed differences in lipid patterns. Among the patient isolates, M. tuberculosis strain 249 possessed a specific glycolipid with 2-O-methyl-L -rhamnose and L-rhamnose, which is rarely found in other strains. This glycolipid showed serological activity against the sera of tuberculosis patients, even if the reactivity was not as strong as trehalose dimycolate. It also showed the inhibition of phagosome-lysosome fusion in macrophages, suggesting involvement with virulence. These results suggest that RFLP analysis using IS6110 is useful for clustering the human isolates of M. tuberculosis, however, for further strain differentiation on virulence, a lipid analysis provides more information.     

Genetic Diversity of Mycobacterium tuberculosis from Guadalajara,Mexico and Identification of a Rare Multidrug Resistant Beijing Genotype

Determining the genetic diversity of M. tuberculosis strains allows identification of the distinct Mycobacterium tuberculosis genotypes responsible for tuberculosis in different regions. Several studies have reported the genetic diversity of M. tuberculosis strains in Mexico, but little information is available from the state of Jalisco. Therefore, the aim of this study was to determine the genetic diversity of Mycobacterium tuberculosis clinical isolates from Western Mexico. Sixty-eight M. tuberculosis isolates were tested for susceptibility to first-line drugs using manual Mycobacteria Growth Indicator Tube method and genotyped using spoligotyping and IS6110-restriction fragment length polymorphism (RFLP) pattern analyses. Forty-seven (69.1%) isolates were grouped into 10 clusters and 21 isolates displayed single patterns by spoligotyping. Three of the 21 single patterns corresponded to orphan patterns in the SITVITWEB database, and 1 new type that contained 2 isolates was created. The most prevalent lineages were T (38.2%), Haarlem (17.7%), LAM (17.7%), X (7.4%), S (5.9%), EAI (1.5%) and Beijing (1.5%). Six (12.8%) of the clustered isolates were MDR, and type 406 of the Beijing family was among the MDR isolates. Seventeen (26.2%) isolates were grouped into 8 clusters and 48 isolates displayed single patterns by IS6110-RFLP. Combination of IS6110-RFLP and spoligotyping reduced the clustering rate to 20.0%. The results show that T, Haarlem, and LAM are predominant lineages among clinical isolates of M. tuberculosis in Guadalajara, Mexico. Clustering rates indicated low transmission of MDR strains. We detected a rare Beijing genotype, SIT406, which was a highly resistant strain. This is the first report of this Beijing genotype in Latin America.     

IS6110 Transposition and Evolutionary Scenario of the Direct Repeat Locus in a Group of Closely Related Mycobacterium tuberculosis Strains

In recent years, various polymorphic loci and multicopy insertion elements have been discovered in the Mycobacterium tuberculosis genome, such as the direct repeat (DR) locus, the major polymorphic tandem repeats, the polymorphic GC-rich repetitive sequence, IS6110, and IS1081. These, especially IS6110 and the DR locus, have been widely used as genetic markers to differentiate M. tuberculosis isolates and will continue to be so used, due to the conserved nature of the genome of M. tuberculosis. However, little is known about the processes involved in generating these or of their relative rates of change. Without an understanding of the biological characteristics of these genetic markers, it is difficult to use them to their full extent for understanding the population genetics and epidemiology of M. tuberculosis. To address these points, we identified a cluster of 7 isolates in a collection of 101 clinical isolates and investigated them with various polymorphic genetic markers, which indicated that they were highly related to each other. This cluster provided a model system for the study of IS6110 transposition, evolution at the DR locus, and the effects of these on the determination of evolutionary relationships among M. tuberculosis strains. Our results suggest that IS6110 restriction fragment length polymorphism patterns are useful in grouping closely related isolates together; however, they can be misleading if used for making inferences about the evolutionary relationships between closely related isolates. DNA sequence analysis of the DR loci of these isolates revealed an evolutionary scenario, which, complemented with the information from IS6110, allowed a reconstruction of the evolutionary steps and relationships among these closely related isolates. Loss of the IS6110 copy in the DR locus was noted, and the mechanisms of this loss are discussed.     

Organization and copy number of initiator tRNA genes in slow- and fast-growing mycobacteria

We have previously reported the isolation and characterization of a functional initiator tRNA gene,metA, and a second initiator tRNA-like sequence,metB, fromMycobacterium tuberculosis. Here we describe the fine mapping of the initiator tRNA gene locus of the avirulent (H37Ra) and virulent (H37Rv) strains ofM. tuberculosis. The genomic blot analyses show that the 1.7 kb (harbouringmetE) and the 6.0 kb BamHI (harbouringmetA) fragments are linked. Further, sequencing of a portion of the 6.0kb fragment, in conjunction with the sequence of the 1.7 kb fragment confirmed the presence of an IS6110 element in the vicinity ofmetB. The IS element is flanked by inverted (28 bp, with 3 contiguous mismatches in the middle) and direct (3 bp) repeats considered to be the hallmarks of IS6110 integration sites. The organization of the initiator tRNA gene locus is identical in both the H37Ra and H37Rv strains and they carry a single copy of the functional initiator tRNA gene. Interestingly, the fast growingMycobacterium smegmatis also bears a single initiator tRNA gene. This finding is significant in view of the qualitative differences in total tRNA pools and the copy number of rRNA genes in the fast- and slow-growing mycobacteria. Finally, we discuss hypotheses related to the origin ofmetB inM. tuberculosis.     

A Novel Molecular Strategy for Surveillance of Multidrug Resistant Tuberculosis in High Burden Settings

Background

In South Africa and other high prevalence countries, transmission is a significant contributor to rising rates of multidrug resistant tuberculosis (MDR-TB). Thus, there is a need to develop an early detection system for transmission clusters suitable for high burden settings. We have evaluated the discriminatory power and clustering concordance of a novel and simple genotyping approach, combining spoligotyping with pncA sequencing (SpoNC), against two well-established methods: IS6110-RFLP and 24-loci MIRU-VNTR.

Methods

A total of 216 MDR-TB isolates collected from January to June 2010 from the NHLS Central TB referral laboratory in Braamfontein, Johannesburg, representing a diversity of strains from South Africa, were included. The isolates were submitted for genotyping, pncA sequencing and analysis to the Centre for Tuberculosis in South Africa and the Public Health Research Institute Tuberculosis Center at Rutgers University in the United States. Clustering rates, Hunter-Gaston Discriminatory Indexes (HGI) and Wallace coefficients were compared between the methods.

Results

Overall clustering rates were high by both IS6110-RFLP (52.8%) and MIRU-VNTR (45.8%), indicative of on-going transmission. Both 24-loci MIRU-VNTR and IS6110-RFLP had similar HGI (0.972 and 0.973, respectively), with close numbers of unique profiles (87 vs. 70), clustered isolates (129 vs. 146), and cluster sizes (2 to 26 vs. 2 to 25 isolates). Spoligotyping alone was the least discriminatory (80.1% clustering, HGI 0.903), with 28 unique types. However, the discriminatory power of spoligotyping was improved when combined with pncA sequencing using the SpoNC approach (61.8% clustering, HGI 0.958). A high proportion of MDR-TB isolates had mutations in pncA (68%, n = 145), and pncA mutations were significantly associated with clustering (p = 0.007 and p = 0.0013 by 24-loci MIRU-VNTR and IS6110-RFLP, respectively), suggesting high rates of resistance to pyrazinamide among all MDR-TB cases and particularly among clustered cases.

Conclusion

We conclude that SpoNC provides good discrimination for MDR-TB surveillance and early identification of outbreaks in South Africa, with 24-loci MIRU-VNTR applied for pncA wild-type strains as needed.     

IS6110-Mediated Deletions of Wild-Type Chromosomes of Mycobacterium tuberculosis

The ipl locus is a site for the preferential insertion of IS6110 and has been identified as an insertion sequence, IS1547, in its own right. Various deletions around the ipl locus of clinical isolates of Mycobacterium tuberculosis were identified, and these deletions ranged in length from several hundred base pairs up to several kilobase pairs. The most obvious feature shared by these deletions was the presence of an IS6110 copy at the deletion sites, which suggested two possible mechanisms for their occurrence, IS6110 transposition and homologous recombination. To clarify the mechanism, an investigation was conducted; the results suggest that although deletion transpositionally mediated by IS6110 was a possibility, homologous recombination was a more likely one. The implications of such chromosomal rearrangements for the evolution of M. tuberculosis, for IS6110-mediated mutagenesis, and for the development of genetic tools are discussed. The deletion of genomic DNA in isolates of M. tuberculosis has previously been noted at only a few sites. This study examined the deletional loss of genetic material at a new site and suggests that such losses may occur elsewhere too and may be more prevalent than was previously thought. Distinct from the study of laboratory-induced mutations, the detailed analysis of clinical isolates, in combination with knowledge of their evolutionary relationships to each other, gives us the opportunity to study mutational diversity in isolates that have survived in the human host and therefore offers a different perspective on the importance of particular genetic markers in pathogenesis.     

Comparison between RFLP and MIRU-VNTR Genotyping of Mycobacterium tuberculosis Strains Isolated in Stockholm 2009 to 2011

Our aim was to analyze the difference between methods for genotyping of Mycobacterium tuberculosis complex isolates. We collected genotyping results from Restriction Fragment Length Polymorphism (RFLP) and Mycobacterial Interspersed Repetitive Units - Variable Numbers of Tandem Repeat (MIRU-VNTR) in a geographically limited area (Stockholm) during a period of three years. The number and proportion of isolates belonging to clusters was reduced by 45 and 35% respectively when combining the two methods compared with using RFLP or MIRU-VNTR only. The mean size of the clusters was smaller when combining methods and smaller with RFLP compared to MIRU-VNTR. In clusters with confirmed epidemiological links RFLP coincided slightly better than MIRU-VNTR but where there was a difference, the variation in MIRU-VNTR pattern was only in a single locus. In isolates with few IS6110 bands in RFLP, MIRU-VNTR differentiated the isolates more, dividing the RFLP clusters. Since MIRU-VNTR is faster and less labour-intensive it is the method of choice for routine genotyping. In most cases it will be sufficient for epidemiological purposes but true clustering might still be considered if there are epidemiological links and the MIRU-VNTR results differ in only one of its 24 loci.     

Genotyping of <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis> using whole genome sequencing

Tuberculosis (TB) is considered one of the most serious infectious diseases worldwide. Effective control of tuberculosis infection involves multiple steps, such as reliable detection, treatment, an epidemiological control as a part of case management, and further surveillance and monitoring of TB spread in the human population. Due to the accelerating advances in molecular biology, especially in DNA sequencing, in the past decade, the application of these methods has become crucial for TB evolution studies, differentiation of Mycobacterium tuberculosis genotypes, and their distribution. Currently, several molecular genetic methods are available. The oldest typing methods (e.g., IS6110-RFLP, spoligotyping, and MIRU-VNTR) can discover the chain of transmission to the patient. Currently, whole genome sequencing facilitates is furthermore able to identify the source of infection, the transmission trays among individuals sharing the same isolate, as well as determination of the TB evolution and its resistance to antituberculotic agents. It is obvious that this technique will become a new gold standard in genotyping methods in tuberculosis molecular epidemiological studies. In this article, molecular genetic typing methods with a special focus on whole genome sequencing and data management are reviewed.     

Distribution of spoligotyping defined genotypic lineages among drug-resistant Mycobacterium tuberculosis complex clinical isolates in Ankara, Turkey

Background

Investigation of genetic heterogeneity and spoligotype-defined lineages of drug-resistant Mycobacterium tuberculosis clinical isolates collected during a three-year period in two university hospitals and National Tuberculosis Reference and Research Laboratory in Ankara, Turkey.

Methods and Findings

A total of 95 drug-resistant M. tuberculosis isolates collected from three different centers were included in this study. Susceptibility testing of the isolates to four major antituberculous drugs was performed using proportion method on Löwenstein–Jensen medium and BACTEC 460-TB system. All clinical isolates were typed by using spoligotyping and IS6110-restriction fragment length polymorphism (RFLP) methods. Seventy-three of the 95 (76.8%) drug resistant M. tuberculosis isolates were isoniazid-resistant, 45 (47.4%) were rifampicin-resistant, 32 (33.7%) were streptomycin-resistant and 31 (32.6%) were ethambutol-resistant. The proportion of multidrug-resistant isolates (MDR) was 42.1%. By using spoligotyping, 35 distinct patterns were observed; 75 clinical isolates were grouped in 15 clusters (clustering rate of 79%) and 20 isolates displayed unique patterns. Five of these 20 unique patterns corresponded to orphan patterns in the SITVIT2 database, while 4 shared types containing 8 isolates were newly created. The most prevalent M. tuberculosis lineages were: Haarlem (23/95, 24.2%), ill-defined T superfamily (22/95, 23.2%), the Turkey family (19/95, 20%; previously designated as LAM7-TUR), Beijing (6/95, 6.3%), and Latin-America & Mediterranean (LAM, 5/95 or 5.3%), followed by Manu (3/95, 3.2%) and S (1/95, 1%) lineages. Four of the six Beijing family isolates (66.7%) were MDR. A combination of IS6110-RFLP and spoligotyping reduced the clustering rate from 79% to 11.5% among the drug resistant isolates.

Conclusions

The results obtained showed that ill-defined T, Haarlem, the Turkey family (previously designated as LAM7-TUR family with high phylogeographical specifity for Turkey), Beijing and LAM were predominant lineages observed in almost 80% of the drug-Resistant M. tuberculosis complex clinical isolates in Ankara, Turkey.     

Clinical evaluation of mtp40 polymerase chain reaction for the diagnosis of extra pulmonary tuberculosis

Rapid and sensitive detection of Mycobacterium tuberculosis from patient samples is vital for clinical diagnosis and treatment. The emergence of M. tuberculosis strains with either no copies or only a single copy of IS6110 in Asian countries makes the standard PCR based diagnosis of M. tuberculosis using IS6110 not reliable. We studied the diagnostic efficacy of the in-house PCR amplification of the candidate gene mtp40 as an alternative to IS6110 element based diagnosis. Clinical samples included pulmonary and extra-pulmonary specimens from TB suspected patients residing in Puducherry, South India and were analyzed using in-house PCR procedures targeting IS6110 element and mtp40 genes. Out of 317 clinical specimens analyzed, 132 (41.6 %) and 114 (36 %) were found positive for mtp40 PCR and IS6110 PCR, respectively. However, 18 specimens that were found to negative for IS6110 PCR were found positive for mtp40 PCR, which was further confirmed by DNA sequencing method. PCR amplification of mtp40 gene for the diagnosis of M. tuberculosis in clinical samples is fast, sensitive, and further identified clinical strains that lack IS6110 element in this region. It is clearly demonstrated that there is a significant difference between the two PCR procedures and the sensitivity and specificity levels of mtp40 PCR were found to be higher when compared with DNA sequencing method. Thus, mtp40 based PCR technique will be beneficial in diagnosis of TB where M. tuberculosis strains lack of IS6110 element is predominant.     

Restriction Fragment Length Polymorphism Analysis of Epidemiologically Related Mycobacterium tuberculosis Isolates

Restriction fragment length polymorphism (RFLP) analysis of a large number of Japanese isolates of Mycobacterium tuberculosis, containing isolates from small outbreaks of M. tuberculosis infection, and clinical isolates of M. bovis BCG, was carried out using a DNA probe derived from the insertion sequence IS986. Clinical isolates of M. tuberculosis had a high degree of RFLP. The occurrences of the IS element varied from 1 to 19, the majority of isolates having 8 to 15 copies. Very similar fingerprints, however, were seen among strains isolated in the Kanto district. In particular, 3 strains were of the same pattern with or without an additional band. Similarity of the banding patterns of strains islated in the same district was observed in other areas. Six groups of strains, each group arising from a suspected common source of infection, were analyzed. Of these, 5 showed identical fingerprints within each group, but one showed different fingerprints. RFLP patterns of three strains isolated from individuals with lymphadenitis developed about two months after BCG vaccination, and one strain isolated from a bladder cancer patient with BCG instillation therapy were identical to those of BCG-Tokyo which had been used for the vaccination and therapy. These results confirm that RFLP analysis using IS986 is a suitable tool for epidemiology of tuberculosis.     

Rapid detection of Mycobacterium tuberculosis in clinical samples by multiplex polymerase chain reaction (mPCR)

Although the multi-copy and specific element IS6110 provides a good target for the detection of Mycobacterium tuberculosis complex by PCR techniques, the emergence of IS6110-negative strains suggested that false negative may occur if IS6110 alone is used as the target for detection. In this report, a multiplex polymerase chain reaction (mPCR) system was developed using primers derived from the insertion sequence IS6110 and an IS-like elements designated as B9 (GenBank accession no. U78639.1) to overcome the problem of detecting negative or low copy IS6110 containing strains of M. tuberculosis. The mPCR was evaluated using 346 clinical samples which included 283 sputum, 19 bronchial wash, 18 pleural fluid, 9 urine, 7 CSF, 6 pus, and 4 gastric lavage samples. Our results showed that the sensitivity (93.1 %) and specificity (89.6 %) of the mPCR system exceeds that of the conventional method of microscopy and culture. The mPCR assay provides an efficient strategy to detect and identify M. tuberculosis from clinical samples and enables prompt diagnosis when rapid identification of infecting mycobacteria is necessary.     

Isolation, characterization and transposition of an (IS2)2 intermediate

We have isolated and characterized a dimer derivative of the extensively studiedEscherichia coli insertion sequence IS2. The dimer structure — called (IS2)₂ — consists of two IS2 elements arranged as a direct repeat, separated by 1 bp. The junction between the (IS2)₂ dimer and target sequences is located at various positions in independent isolates; however, one position was preferred. The transposition of (IS2)₂ into a target plasmid resulted in cointegrate-type structures. The transposition frequency of the (IS2)₂ dimer itself was significantly higher than that of the isogenic monomer IS2 insertion. The poor stability and high activity of (IS2)₂ indicates that this is an active transposition intermediate. The mode of transposition of (IS2)₂ is analogous to the joined dimer model described in the case of (IS21)₂ and (IS30)₂.     

Optimized molecular resolution of cross-contamination alerts in clinical mycobacteriology laboratories

Background

The phenomenon of misdiagnosing tuberculosis (TB) by laboratory cross-contamination when culturing Mycobacterium tuberculosis (MTB) has been widely reported and it has an obvious clinical, therapeutic and social impact. The final confirmation of a cross-contamination event requires the molecular identification of the same MTB strain cultured from both the potential source of the contamination and from the false-positive candidate. The molecular tool usually applied in this context is IS6110-RFLP which takes a long time to provide an answer, usually longer than is acceptable for microbiologists and clinicians to make decisions. Our purpose in this study is to evaluate a novel PCR-based method, MIRU-VNTR as an alternative to assure a rapid and optimized analysis of cross-contamination alerts.

Results

MIRU-VNTR was prospectively compared with IS6110-RFLP for clarifying 19 alerts of false positivity from other laboratories. MIRU-VNTR highly correlated with IS6110-RFLP, reduced the response time by 27 days and clarified six alerts unresolved by RFLP. Additionally, MIRU-VNTR revealed complex situations such as contamination events involving polyclonal isolates and a false-positive case due to the simultaneous cross-contamination from two independent sources.

Conclusion

Unlike standard RFLP-based genotyping, MIRU-VNTR i) could help reduce the impact of a false positive diagnosis of TB, ii) increased the number of events that could be solved and iii) revealed the complexity of some cross-contamination events that could not be dissected by IS6110-RFLP.     

Detection of streptomycin resistance in Mycobacterium tuberculosis clinical isolates from China as determined by denaturing HPLC analysis and DNA sequencing

China is regarded by the World Health Organization as a major hot-spot region for Mycobacterium tuberculosis infection. Streptomycin has been deployed in China for over 50 years and is still widely used for tuberculosis treatment. We have developed a denaturing HPLC (DHPLC) method for detecting various gene mutations conferring drug resistance in M. tuberculosis. The present study focused on rpsL and rrs mutation analysis. Two hundred and fifteen M. tuberculosis clinical isolates (115 proved to be streptomycin-resistant and 100 susceptible by a routine proportional method) from China were tested to determine the streptomycin minimal inhibitory concentration (MIC), and subjected to DHPLC and concurrent DNA sequencing to determine rpsL and rrs mutations. The results showed that 85.2% (98/115) of streptomycin-resistant isolates harbored rpsL or rrs mutation, while rpsL mutation (76.5%, 88/115) dominated. MIC of 98 mutated isolates revealed no close correlation between mutation types and levels of streptomycin resistance. No mutation was found in any of the susceptible isolates. The DHPLC results were completely consistent with those of sequencing. The DHPLC method devised in this study can be regarded as a useful and powerful tool for detection of streptomycin resistance. This is the first report to describe DHPLC analysis of mutations in the rpsL and rrs genes of M. tuberculosis in a large number of clinical isolates.     

Models of deletion for visualizing bacterial variation: an application to tuberculosis spoligotypes

Background  

Molecular typing methods are commonly used to study genetic relationships among bacterial isolates. Many of these methods have become standardized and produce portable data. A popular approach for analyzing such data is to construct graphs, including phylogenies. Inferences from graph representations of data assist in understanding the patterns of transmission of bacterial pathogens, and basing these graph constructs on biological models of evolution of the molecular marker helps make these inferences. Spoligotyping is a widely used method for genotyping isolates of Mycobacterium tuberculosis that exploits polymorphism in the direct repeat region. Our goal was to examine a range of models describing the evolution of spoligotypes in order to develop a visualization method to represent likely relationships among M. tuberculosis isolates.     

Highly polymorphic variable-number tandem repeats loci for differentiating Beijing genotype strains of Mycobacterium tuberculosis in Shanghai, China

The PCR-based variable-number tandem repeats (VNTR) typing method is a very promising tool for the molecular epidemiological study of Mycobacterium tuberculosis. The discriminatory power of the VNTR loci that were optimized in many previous studies has not been evaluated in Shanghai, an area where Beijing genotype strains dominate. In the present study, we first performed a literature search to identify VNTR loci that were at least 45 bp in length. Second, we determined the Hunter-Gaston discriminatory index (HGI) values of each of the 45 VNTR loci that we identified, using Beijing genotype strains from a 'test set' of isolates from a population with low migration in Chongming Island, Shanghai, China. Third, we optimized two sets of VNTR loci, which we named VNTR-7 and VNTR-16. The HGI value of VNTR-7 was slightly lower than that of IS6110 restriction fragment length polymorphisms (RFLP), and the HGI values of VNTR-16 and IS6110 RFLP were comparable. Our results suggest that VNTR-7, followed by VNTR-16 and IS6110 RFLP, can be used routinely as a tool to discriminate between M. tuberculosis isolates in population-based epidemiologic studies of M. tuberculosis transmission in Shanghai, China.