Search tips
Search criteria 


Logo of jcmPermissionsJournals.ASM.orgJournalJCM ArticleJournal InfoAuthorsReviewers
J Clin Microbiol. 2010 May; 48(5): 1758–1766.
Published online 2010 March 10. doi:  10.1128/JCM.02255-09
PMCID: PMC2863856

Identification and Genotyping of Mycobacterium tuberculosis Complex Species by Use of a SNaPshot Minisequencing-Based Assay[down-pointing small open triangle]


The aim of the present study was to investigate the use of the SNaPshot minisequencing method for the identification of Mycobacterium tuberculosis complex (MTBC) isolates to the species level and for further genotyping of M. tuberculosis isolates. We developed an innovative strategy based on two multiplex allele-specific minisequencing assays that allowed detection of eight species-specific and eight lineage-specific single nucleotide polymorphisms (SNPs). Each assay consisted of an eightplex PCR amplification, followed by an eightplex minisequencing reaction with the SNaPshot multiplex kit (Applied Biosystems) and, finally, analysis of the extension products by capillary electrophoresis. The whole strategy was developed with a panel of 56 MTBC strains and 15 negative controls. All MTBC strains tested except one M. africanum clinical isolate were accurately identified to the species level, and all M. tuberculosis isolates were successfully further genotyped. This two-step strategy based on SNaPshot minisequencing allows the simultaneous differentiation of closely related members of the MTBC, the distinction between principal genetic groups, and the characterization of M. tuberculosis isolates into one of the seven prominent SNP cluster groups (SCGs) and could be a useful tool for diagnostic and epidemiological purposes.

Although tuberculosis (TB) is an age-old disease, it still represents a major health problem worldwide, accounting for nearly 2 million deaths annually (World Health Organization, Tuberculosis Facts 2009 []). Besides the need for an improved therapy and for new vaccines, effective TB control requires the initiation of appropriate therapy as well as an increased understanding of its epidemiology.

The causative agents of TB in humans and animals, including Mycobacterium tuberculosis, M. africanum, M. bovis, M. canettii, M. microti, M. caprae, and M. pinnipedii, form the Mycobacterium tuberculosis complex (MTBC). Although they differ widely in terms of their host tropisms, phenotypes, and pathogenicities, all members of the MTBC are closely related genetically (51). M. tuberculosis species, the most common pathogen in humans, can be further divided into genetic groups that also show differences in their levels of virulence, immunogenicities, and geographical distributions (21). On the one hand, it is important to differentiate MTBC species to distinguish between strict human and zoonotic TB and to initiate an appropriate therapy (15). In particular, the distinction between M. tuberculosis and M. bovis is necessary, as the latter species is naturally resistant to the antituberculous drug pyrazinamide (48). On the other hand, genotyping of M. tuberculosis isolates is useful as a means of addressing evolutionary questions but also as a means of surveying the transmission dynamics of this pathogen and identifying new outbreaks.

Identification of MTBC isolates to the species level is so far routinely performed by analysis of the phenotypic and biochemical characteristics of the bacteria after culture. However, this is a time-consuming and subjective process, and that is why various molecular methods have been developed in recent years. Most molecular methods for MTBC species differentiation described in the literature are based on the analysis of genomic deletions by PCR, followed by agarose gel electrophoresis (30, 40, 56), or are based on the detection of single nucleotide polymorphisms (SNPs) by PCR-restriction fragment length polymorphism (PCR-RFLP) analysis, which consists of digestion of the PCR product with a restriction enzyme, followed by agarose gel electrophoresis (3, 13, 23, 33, 38). A commercial kit based on a DNA strip assay for the analysis of gyrB gene SNPs, 23S rRNA, and the RD1 deletion, named GenoType MTBC (Hain Lifescience GmbH, Nehren, Germany), has recently become available (37, 42, 50). This line probe assay enables identification of the presence of various members of the MTBC but does not differentiate M. canettii from M. tuberculosis and M. africanum type I from M. pinnipedii.

Various molecular methods for genotyping of M. tuberculosis have also been developed and were recently reviewed by Mathema and colleagues (36). The most commonly used methods include the “gold standard” analysis of highly conserved DNA fingerprinting patterns obtained by RFLP analysis of the insertion sequence IS6110 (IS6110-RFLP fingerprinting) (54), the study of variations within the genomic direct repeat region by spoligotyping (10, 32), and the typing of mycobacterial interspersed repetitive-unit-variable number of DNA tandem repeats (MIRU-VNTR) (52). They are powerful molecular epidemiological typing methods since they are based on the analysis of mobile DNA elements (IS6110-RFLP fingerprinting) or repetitive DNA elements (spoligotyping and MIRU-VNTR typing) that change quite rapidly and thus provide a high degree of discriminatory power, but they are less useful for definition of phylogenetic groupings (14, 21). Recently, two studies have revealed that SNPs are valuable phylogenetically informative markers (18, 26), and the use of a minimal number of sets of SNPs has been proposed as a means of resolving M. tuberculosis and M. bovis isolates into seven genetic groups called SNP cluster groups (SCGs) (1, 18).

Thus, SNPs can be useful markers for both the identification of MTBC species and the genotyping of M. tuberculosis isolates. We describe here the development of an SNP typing-based strategy that was designed to simultaneously distinguish members of the MTBC to the species level and M. tuberculosis lineages. Among the various SNP typing methodologies described in the literature, we selected the SNaPshot minisequencing-based approach because of its high multiplexing capacity, robustness, and extreme sensitivity (22, 49). Moreover, it is a relatively simple and affordable method that requires the use of a thermal cycler and a genetic analyzer, which are types of equipment commonly available in molecular biology laboratories. This approach is based on the single-base extension (SBE) of an unlabeled minisequencing primer that anneals one base upstream of the relevant SNP with a fluorochrome-labeled dideoxynucleotide (ddNTP). The allelic state is then determined after separation of the extension products and the detection of fluorescence by capillary electrophoresis. This SNP typing method is currently widely used in the fields of forensic and population genetics (8, 9, 25, 44, 47) and has recently raised interest in clinical research (2, 5, 17, 31, 39). It was also applied in a large-scale study for the elucidation of the M. tuberculosis population substructure and evolution (26, 27), and a recent study reports on its use for the identification of six common mycobacterial species (M. tuberculosis, M. avium, M. intracellulare, M. chelonae, M. kansasii, and M. gordonae) (55). We propose here its application for the rapid and simultaneous identification of MTBC species and the genotypic characterization of M. tuberculosis isolates.


Bacterial strains tested and reference tests.

The following 56 MTBC strains were used for the development of the assays: M. tuberculosis H37Rv ATCC 27294; M. bovis CIP 102426; M. bovis BCG CIP 105226; and clinical isolates of M. tuberculosis (n = 35), M. bovis (n = 6), M. bovis BCG (n = 4), M. africanum (n = 4), M. canettii (n = 1), M. caprae (n = 1), M. microti (n = 1), and M. pinnipedii (n = 1) (the M. pinnipedii isolate was recovered from a tapir, and the others were isolated from human clinical samples). Eleven clinical isolates of mycobacteria other than members of the M. tuberculosis complex (MOTT), comprising M. fortuitum (n = 2), M. kansasii (n = 2), M. abscessus (n = 1), M. avium (n = 3), M. chelonae (n = 2), and M. gordonae (n = 1), as well as DNA samples from four bacteria that do not belong to the genus Mycobacterium (Nocardia nova, Corynebacterium amycolatum, Staphylococcus aureus, and Escherichia coli) were used as negative controls.

DNA samples of control strains M. tuberculosis H37Rv ATCC 27294, M. bovis CIP 102426, and M. bovis BCG CIP 105226 were provided by the Laboratory of Bacteriology of Strasbourg, Strasbourg, France. DNA extraction was performed from solid cultures by use of a MagNA Pure LC DNA III isolation kit (Roche Diagnostics, Indianapolis, IN), according to the manufacturer's protocol. Samples of clinical isolates, which were prepared by simple thermolysis, were provided by the Laboratories of Bacteriology of Strasbourg, France, and Toulouse, France, and the Centre National de Référence des Mycobactéries in Paris, France. A loop of culture (i.e., bacteria grown on Lowenstein-Jensen or Coletsos medium) was suspended in 300 to 500 μl water, and the mixture was heated at 95°C for 15 to 20 min. Mycobacterial cells were disrupted by sonication for 15 min or heat shock treatment (1 min at 95°C, followed by 1 min on ice, repeated five times). The supernatants obtained after centrifugation at full speed for 5 min contained the mycobacterial extracts used for the PCR amplifications.

All isolates were previously identified to the species level by phenotypic and biochemical characterization methods after culture (43) or a gene probe assay, according to the manufacturer's protocol (GenoType MTBC; Hain Lifescience GmbH). Ten samples, including seven clinical samples of M. tuberculosis, one clinical sample of M. africanum, M. tuberculosis H37Rv, and M. bovis control strains, were further genotyped by 24-locus MIRU-VNTR typing (52). The person who performed the MIRU-VNTR typing was blinded to the species identification and to the results of SNP typing.

SNP selection.

A set of 16 well-characterized SNPs was chosen from the literature (1, 18, 19, 24, 29, 33, 38, 51). The following eight SNPs were selected for use for the differentiation of MTBC species: hsp65631 (C → T in M. canettii), gyrB(675) (numbers in parentheses refer to the SNP position on the gene) (C → T in M. microti), gyrB(756) (G → A in M. caprae and M. bovis), gyrB(1410) (C → T in M. bovis and M. bovis BCG), 16S rRNA1249 (T → C in M. pinnipedii), and the three SNPs that determined the so-called principal genetic groups (PGGs), which are katG203 (ACT → ACC for PGG 1a and PGG 1b differentiation), katG463 (CTG → CGG for PGG 2), and gyrA95 (ACC → AGC for PGG 3). For genotyping of the M. tuberculosis isolates, the minimal SNP set that has recently been proposed for use by Alland and coworkers was chosen (the bases at the following SNP positions in H37Rv: 1977, 3352929 [instead of 54394], 74092, 105139, 2460626 [instead of 144390], 232574, 311613, and 913274) (1). These SNPs define the seven prominent SCGs that were first described by Filliol et al. (18). The sequence variation and lineage defined by them can be obtained from previously published reports (1, 18). It is important to note that the SNP at position 2154724, which also belongs to the minimal set, corresponds to katG463, which was already selected for use in the first SNP set. In our study, the first eight SNPs were qualified as “species specific,” whereas the last eight were “lineage specific” or “genotype specific.”

Design and validation of PCR and minisequencing primers.

Primer design was performed according to the recommendations of Sanchez et al. (45) so that the species-specific and lineage-specific SNPs were typed in two different multiplex assays: multiplex PCR 1 (mPCR 1) followed by SBE 1 and mPCR 2 followed by SBE 2, respectively. The sizes of the amplicons ranged from 72 to 150 bp with mPCR 1 and 81 to 141 bp with mPCR 2. The minisequencing primers were tailed at the 5′ end with a nonhomologous sequence (34) and poly(C), if necessary, to produce extension products that ranged in length from 28 to 76 nucleotides (nt) with SBE 1 and from 31 to 73 nt with SBE 2 and that differed in length from each other by 6 nt to allow sufficient separation by capillary electrophoresis. The sequences of the PCR and minisequencing primers used in this study are shown in Tables Tables11 and and2,2, respectively. All primers were tested in singleplex reactions with 1 ng template DNA under the conditions outlined by Sanchez et al. (45). Since the M. tuberculosis genome has a high GC content, 1 M betaine (Sigma-Aldrich, St. Louis, MO) was added to the singleplex PCR mixtures in order to reduce the formation of secondary structures in GC-rich regions (28). The amplification products were then verified by classical sequencing.

PCR primers for mPCR 1 and mPCR 2 used in this study
Minisequencing primers for SBE 1 and SBE 2 used in this study

Multiplex PCR amplifications and purification of PCR products.

Each mPCR was performed with 1 ng template DNA or 1 to 10 μl mycobacterial extract for clinical samples in a 50-μl final volume composed of 1× PCR Gold buffer (Applied Biosystems [AB], Foster City, CA), 8 mM MgCl2, 400 μM each deoxynucleoside triphosphate (dNTP), 0.2 to 0.5 μM each primer, and 2 U of AmpliTaq Gold DNA polymerase (AB). The concentrations of the PCR primers used in each reaction mixture are specified in Table Table1.1. The thermal cycling consisted of a first denaturation step at 94°C for 5 min, followed by 33 cycles of denaturation at 95°C for 30 s, annealing at 60°C for 30 s, and extension at 65°C for 30 s and with a final extension at 65°C for 7 min. Excess PCR primers and dNTPs were removed by using a NucleoSpin Extract II kit (Macherey-Nagel, Düren, Germany), according to the manufacturer's recommendations. The PCR products were eluted from the filter in 50 μl sterile water.

Multiplex SBE reactions and purification of extension products.

Each SBE reaction was carried out in a 10 μl final volume containing 4 μl of SNaPshot multiplex ready reaction mixture (AB), 20 mM ammonium sulfate, 3 μl cleaned PCR products, and 0.13 to 0.58 μM each minisequencing primer. The concentrations of primers in the reaction mixture are specified in Table Table2.2. Extension was performed for 25 cycles of denaturation at 96°C for 10 s, annealing at 50°C for 5 s, and extension at 60°C for 30 s. Unincorporated ddNTPs were removed by addition of 1 U of shrimp alkaline phosphatase (GE Healthcare, Little Chalfont, United Kingdom) to the SBE mixtures and incubation at 37°C for 1 h, followed by incubation at 75°C for 15 min.

Capillary electrophoresis and data analysis.

Two microliters of the purified SBE products was mixed with 17.6 μl Hi-Di formamide (AB) and 0.4 μl of GeneScan-120 LIZ internal size standard (AB), and the mixture was further analyzed by capillary electrophoresis on an ABI Prism 3100 genetic analyzer (AB) with 36-cm capillary arrays and performance-optimized polymer 4 (POP-4) (AB). The data were analyzed and the alleles were automatically called by the use of GeneMapper ID (version 3.2.1) software (AB).


Multiplex reactions for species differentiation within MTBC.

The eight species-specific SNPs selected in this study were successfully genotyped for all 56 MTBC samples tested except M. canettii, for which the gyrA95 locus failed to be amplified. The genotype data for the 56 MTBC samples are given in Table Table3.3. For each sample, the MTBC species and PGG inferred from the genotype data were always concordant with the previous identification by use of a reference test and data from the literature (29). The reference strain, M. tuberculosis H37Rv, was identified as M. tuberculosis PGG3. The clinical M. tuberculosis samples were found to belong to PGG1b (n = 7), PGG2 (n = 23), and PGG3 (n = 5). All M. bovis and M. bovis BCG samples tested were found to belong to PGG1a, and only those samples had both the A allele for SNP position 756 and the T allele for the SNP position 1410 on the gyrB gene. M. caprae was easily distinguished from M. bovis because it had the A allele for SNP position 756 and the ancestral C allele for SNP position 1410 on the gyrB gene. M. microti and M. pinnipedii samples were assigned to PGG1a, and M. canettii was assigned to PGG1b. These three species were easily distinguishable from other species. Indeed, for all MTBC samples tested, the allele at the hsp65631 locus was always C except in M. canettii (in which it was T), the allele at 16S rRNA1249 locus was always T except in M. pinnipedii (in which it was C), and gyrB gene SNP position 675 was always C except in M. microti (in which it was T). Finally, three M. africanum samples were differentiated from the other species because the ancestral allele was detected for all SNPs, which also enabled us to determine that they belonged to PGG1a and thus probably belonged to M. africanum subtype Ia, as referred to by Huard and colleagues (29). The fourth M. africanum sample tested was assigned to PGG1b since it had a C allele at the katG203 SNP locus and thus probably belonged to M. africanum subtype Ib. However, the alleles detected from this M. africanum isolate and all PGG1b M. tuberculosis isolates tested were identical for all SNP loci analyzed in this study, showing that this assay does not allow discrimination between M. africanum subtype Ib and M. tuberculosis PGG1b. All negative controls always failed to amplify all loci except for, sometimes, hsp65631 and/or 16S rRNA1249. Examples of electropherograms obtained from each MTBC species tested are shown in Fig. Fig.11.

FIG. 1.
Electropherograms obtained from the seven MTBC species by the multiplex assay for MTBC species identification. Each plot was obtained with GeneMapper ID (version. 3.2.1) software (AB) and shows the relative fluorescence units (RFUs) versus the measured ...
Genotypes obtained for the eight species-specific SNPs typed for the 56 MTBC strains

Multiplex reactions for genotyping of M. tuberculosis isolates (SCG identification).

The laboratory strain M. tuberculosis H37Rv, the M. bovis control strain, all clinical isolates of M. tuberculosis (n = 35), and one clinical M. africanum strain were further successfully characterized into one of the seven SCGs by means of the second multiplex assay: SCG1 (n = 2), SCG2 (n = 3), SCG3a (n = 2), SCG3b (n = 9), SCG5 (n = 14), SCG6a (n = 5), SCG6b (M. tuberculosis H37Rv), and SCG7 (M. bovis). All samples assigned to SCGs 1, 2 and 3a were found to belong to PGG1b, those assigned to SCGs 3b and 5 belonged to PGG2, and finally, those assigned to SCGs 6a and 6b belonged to PGG3, which is consistent with the results reported in the literature (18). Genotype data for the eight lineage-specific SNPs and the MIRU-VNTR typing results obtained for seven clinical M. tuberculosis samples, one clinical sample of M. africanum, M. tuberculosis H37Rv, and the M. bovis control samples are shown in Table Table4.4. The SCGs inferred from SNP typing were concordant with the lineages inferred from MIRU-VNTR typing (21). For example, clinical strain M. tuberculosis 1 was assigned to the Beijing clade by MIRU-VNTR and to SCG2 by SNP typing. According to the literature, Beijing clade isolates are exclusively present in SCG2 (18, 21). Negative controls were always negative for all SNPs except some MOTT, for which position 232574 was amplified.

SNP and MIRU-VNTR genotyping results for 10 MTBC strains


The work presented here describes the successful development of two simple, rapid, and specific assays for the differentiation of MTBC members and for the further classification of M. tuberculosis isolates into previously described phylogenetic groups. These assays are based on the SNaPshot minisequencing approach, which was found to be a reliable and reproducible method for SNP typing. Indeed, the SNP genotyping results were always confirmed by at least two independent PCRs for each mycobacterial sample, and the results were always in complete accordance with the results obtained by previous reference identification tests.

As supported by our results, the first assay for the typing of species-specific SNPs allows the accurate differentiation of all MTBC species except M. africanum subtype Ib and M. tuberculosis PGG1b. The differentiation of these two species could be done by adding the gyrBafr SNP (SNP position 1450 of the gyrB gene, G → T) to the species-specific SNP panel chosen for use in this study (29). Thus, the assay described here differentiates M. canettii from M. tuberculosis and M. africanum type I from M. pinnipedii, which is not feasible by use of the commercial GenoType MTBC DNA strip assay (Hain Lifescience GmbH). It also presents many technical advantages with respect to other PCR-based methods previously developed for MTBC species differentiation. In contrast to the previous methods, which are mostly based on an approach involving agarose gel electrophoresis (e.g., analysis of genomic deletions and PCR-RFLP assays for analysis of SNPs of the gyrB gene) or DNA hybridization (e.g., commercial GenoType MTBC DNA strip assay [Hain Lifescience GmbH]), the assay developed here is based on SNP detection by a multiplex minisequencing reaction, followed by capillary electrophoresis. First, this technique does not need the skilled and subjective interpretation of bands required for the other methods. Indeed, the SNaPshot minisequencing method always produced precise, unambiguous, and objective results, since allele calling was done automatically by the GeneMapper ID analysis software (AB). Second, there is no need to interpret negative results, since the SNaPshot minisequencing approach is not based on the presence or the absence of a specific signal on the resulting electropherograms. For each SNP locus investigated, the specific SBE primers must incorporate one fluorochrome-labeled ddNTP, the one that is complementary to the nucleotide present at the position of interest. On the contrary, molecular identifications based on the analysis of genomic deletions rely on the success or failure of PCR amplification in the case of a specific deletion and are thus limited by the necessity of interpreting negative results. Lastly, high-throughput analysis is possible because the process can be automated: the detection platform (i.e., the genetic analyzer) is automated, the data generated by capillary electrophoresis are automatically stored in a portable format, and alleles can be automatically called by the ad hoc software.

Besides the differentiation of MTBC species, the first assay presented here allows the simultaneous classification of MTBC isolates into three distinct genotypic groups, termed PGGs (PGGs 1a/b, 2 and 3), which is not achieved by previously developed MTBC species identification methods. All MTBC strains can be assigned to one of the three PGGs on the basis of the comparative analysis of three SNPs (katG203, katG463, and gyrA95) (19, 29, 51). To date, M. tuberculosis strains have been found to segregate to PGG1b, PGG2, and PGG3, while other MTBC members are restricted to PGG1a and PGG1b, which was also illustrated by our results (29). In the evolutionary pathway proposed by Sreevatsan et al., M. tuberculosis PGG1 strains are ancestral to PGG2 and PGG3 strains (51). Moreover, those authors observed through their study that PGG3 organisms are rather associated with sporadic TB cases, while PGG1 and PGG2 organisms are associated with clustered TB cases, suggesting differences in transmissibility or virulence between these genotypic groups (51). Thus, the determination of these major genetic groups, which is currently widely used, can be useful for both evolutionary and epidemiological applications.

In recent years, several additional comparative genomic studies of human-adapted members of the MTBC have identified other phylogenetically informative polymorphisms that define discrete strain lineages (1, 4, 10, 18, 20, 26). Because of the highly clonal population structure of the MTBC, the phylogenetic groupings recognized in the different studies are relatively concordant, even though they are based on the analysis of different markers (21). Although molecular epidemiological informative markers like MIRU-VNTR loci have been proposed for use in phylogenetic and population genetic analyses of MTBC, a recent comparative study conducted by Comas et al. revealed that SNPs and large sequence polymorphisms (LSPs) are the most robust and appropriate phylogenetically informative markers (14). The classification of MTBC strains into these distinct lineages is relevant for evolutionary purposes but also for TB control. Indeed, these discrete lineages are associated with particular human populations or geographical regions and show differences in virulence and disease outcomes (12, 20, 21, 35, 41, 53).

The second assay presented here allows the further accurate and unambiguous classification of M. tuberculosis isolates into six more resolved genotypic groups called SCGs. Because these SCGs are congruent with the lineages defined on the basis of other phylogenetically informative markers (21), interlaboratory comparisons can be done even if different markers are studied. This second assay can be run after the first assay, if the PGG determination is not sufficient, or it can be simultaneously performed.

The main advantage of the two-step strategy presented here is that it allows the reliable simultaneous identification of MTBC species and the classification of M. tuberculosis lineages into distinct lineages. It is easy to perform, gives unambiguous results, and can be easily introduced in a bacteriology laboratory with an automated sequencer, a detection platform that is also needed for other common applications like DNA sequencing and MIRU-VNTR typing.

Comas et al. recently suggested a combined SNP typing-MIRU-VNTR typing scheme to generate more accurate data for epidemiological and evolutionary applications (14). Indeed, on the one hand, SNPs are robust phylogenetic markers but offer insufficient discriminatory power for routine molecular epidemiological identifications, and on the other hand, the standard 24 MIRU-VNTR loci are powerful epidemiological markers but are unable to detect all strain lineages and their discriminatory powers vary by strain lineage. The new approach proposed by Comas et al. (14) consists of the identification of the main strain lineages by SNP typing, followed by further molecular epidemiological discrimination by the use of lineage-specific sets of the most discriminatory MIRU-VNTR markers for the lineage of interest. The SNaPshot minisequencing assay used in the present study seems to be an ideal method for SNP typing for such an approach because it uses the same equipment as MIRU-VNTR typing (i.e., a thermal cycler and a genetic analyzer). Furthermore, our study shows that the reliable identification of MTBC species can also be achieved by the SNaPshot minisequencing assay. This suggests that the identification of MTBC species could easily be added as a first step to the scheme proposed by Comas et al. (14).

It is also important to note that this two-step strategy was initially developed for paleomicrobiological analyses, i.e., investigation of ancient human tissues showing lesions suggestive of TB that are suspected to contain low levels of degraded mycobacterial DNA. That is why PCR primers were carefully designed in order to obtain amplicons of less than 150 bp. This feature could be of special interest in a clinical context when laboratories must deal with samples containing small amounts of degraded DNA (e.g., fixed tissues), since short DNA fragments are amplified more efficiently than long fragments in such cases. We are currently testing the use of this two-step strategy with ancient human remains. However, in recent studies describing the development and validation of a similar assay for the typing of human SNPs designed for the screening of ancient human remains, we already reported on the sensitivity and robustness of the SNaPshot minisequencing approach for SNP typing (6, 7). These advantages of this SNP typing method were also reported by other authors, who successfully used it to evaluate ancient human samples or archival clinical samples (6, 7, 11, 16, 22, 31, 46).

To conclude, the present study describes the successful development of reliable SNaPshot minisequencing-based SNP assays for the identification of MTBC species, determination of the principal genetic groups (PGGs 1a/b, 2, and 3), and the further classification of M. tuberculosis isolates into more resolved phylogenetic groups called SCGs. SNaPshot minisequencing for SNP genotyping is currently used for various applications; here, we propose the use of this technique as a complement to current methods for the clinical diagnosis of TB and epidemiological investigations.


We thank Rosine Bauriaud, Cathie Barthel, and Marie Durand for preparation of the DNA from the MTBC strains and Daniel Montagnon and Gerald Millot for technical assistance.


[down-pointing small open triangle]Published ahead of print on 10 March 2010.


1. Alland, D., D. W. Lacher, M. H. Hazbon, A. S. Motiwala, W. Qi, R. D. Fleischmann, and T. S. Whittam. 2007. Role of large sequence polymorphisms (LSPs) in generating genomic diversity among clinical isolates of Mycobacterium tuberculosis and the utility of LSPs in phylogenetic analysis. J. Clin. Microbiol. 45:39-46. [PMC free article] [PubMed]
2. Alvarez-Iglesias, V., F. Barros, A. Carracedo, and A. Salas. 2008. Minisequencing mitochondrial DNA pathogenic mutations. BMC Med. Genet. 9:26. [PMC free article] [PubMed]
3. Arnold, C., L. Westland, G. Mowat, A. Underwood, J. Magee, and S. Gharbia. 2005. Single-nucleotide polymorphism-based differentiation and drug resistance detection in Mycobacterium tuberculosis from isolates or directly from sputum. Clin. Microbiol. Infect. 11:122-130. [PubMed]
4. Baker, L., T. Brown, M. C. Maiden, and F. Drobniewski. 2004. Silent nucleotide polymorphisms and a phylogeny for Mycobacterium tuberculosis. Emerg. Infect. Dis. 10:1568-1577. [PMC free article] [PubMed]
5. Bardien, S., H. Human, T. Harris, G. Hefke, R. Veikondis, H. S. Schaaf, L. van der Merwe, J. H. Greinwald, J. Fagan, and G. de Jong. 2009. A rapid method for detection of five known mutations associated with aminoglycoside-induced deafness. BMC Med. Genet. 10:2. [PMC free article] [PubMed]
6. Bouakaze, C., C. Keyser, S. Amory, E. Crubezy, and B. Ludes. 2007. First successful assay of Y-SNP typing by SNaPshot minisequencing on ancient DNA. Int. J. Legal Med. 121:493-499. [PubMed]
7. Bouakaze, C., C. Keyser, E. Crubezy, D. Montagnon, and B. Ludes. 2009. Pigment phenotype and biogeographical ancestry from ancient skeletal remains: inferences from multiplexed autosomal SNP analysis. Int. J. Legal Med. 123:315-325. [PubMed]
8. Brandstatter, A., A. Salas, H. Niederstatter, C. Gassner, A. Carracedo, and W. Parson. 2006. Dissection of mitochondrial superhaplogroup H using coding region SNPs. Electrophoresis 27:2541-2550. [PubMed]
9. Brion, M., B. Sobrino, A. Blanco-Verea, M. V. Lareu, and A. Carracedo. 2005. Hierarchical analysis of 30 Y-chromosome SNPs in European populations. Int. J. Legal Med. 119:10-15. [PubMed]
10. Brudey, K., J. R. Driscoll, L. Rigouts, W. M. Prodinger, A. Gori, S. A. Al-Hajoj, C. Allix, L. Aristimuno, J. Arora, V. Baumanis, L. Binder, P. Cafrune, A. Cataldi, S. Cheong, R. Diel, C. Ellermeier, J. T. Evans, M. Fauville-Dufaux, S. Ferdinand, D. Garcia de Viedma, C. Garzelli, L. Gazzola, H. M. Gomes, M. C. Guttierez, P. M. Hawkey, P. D. van Helden, G. V. Kadival, B. N. Kreiswirth, K. Kremer, M. Kubin, S. P. Kulkarni, B. Liens, T. Lillebaek, M. L. Ho, C. Martin, C. Martin, I. Mokrousov, O. Narvskaia, Y. F. Ngeow, L. Naumann, S. Niemann, I. Parwati, Z. Rahim, V. Rasolofo-Razanamparany, T. Rasolonavalona, M. L. Rossetti, S. Rusch-Gerdes, A. Sajduda, S. Samper, I. G. Shemyakin, U. B. Singh, A. Somoskovi, R. A. Skuce, D. van Soolingen, E. M. Streicher, P. N. Suffys, E. Tortoli, T. Tracevska, V. Vincent, T. C. Victor, R. M. Warren, S. F. Yap, K. Zaman, F. Portaels, N. Rastogi, and C. Sola. 2006. Mycobacterium tuberculosis complex genetic diversity: mining the Fourth International Spoligotyping Database (SpolDB4) for classification, population genetics and epidemiology. BMC Microbiol. 6:23. [PMC free article] [PubMed]
11. Burger, J., M. Kirchner, B. Bramanti, W. Haak, and M. G. Thomas. 2007. Absence of the lactase-persistence-associated allele in early neolithic Europeans. Proc. Natl. Acad. Sci. U. S. A. 104:3736-3741. [PubMed]
12. Caws, M., G. Thwaites, S. Dunstan, T. R. Hawn, N. T. Lan, N. T. Thuong, K. Stepniewska, M. N. Huyen, N. D. Bang, T. H. Loc, S. Gagneux, D. van Soolingen, K. Kremer, M. van der Sande, P. Small, P. T. Anh, N. T. Chinh, H. T. Quy, N. T. Duyen, D. Q. Tho, N. T. Hieu, E. Torok, T. T. Hien, N. H. Dung, N. T. Nhu, P. M. Duy, N. van Vinh Chau, and J. Farrar. 2008. The influence of host and bacterial genotype on the development of disseminated disease with Mycobacterium tuberculosis. PLoS Pathog. 4:e1000034. [PMC free article] [PubMed]
13. Chimara, E., L. Ferrazoli, and S. C. Leao. 2004. Mycobacterium tuberculosis complex differentiation using gyrB-restriction fragment length polymorphism analysis. Mem. Inst. Oswaldo Cruz 99:745-748. [PubMed]
14. Comas, I., S. Homolka, S. Niemann, and S. Gagneux. 2009. Genotyping of genetically monomorphic bacteria: DNA sequencing in Mycobacterium tuberculosis highlights the limitations of current methodologies. PLoS One 4:e7815. [PMC free article] [PubMed]
15. Djelouadji, Z., D. Raoult, M. Daffe, and M. Drancourt. 2008. A single-step sequencing method for the identification of Mycobacterium tuberculosis complex species. PLoS Negl. Trop. Dis. 2:e253. [PMC free article] [PubMed]
16. Endicott, P., M. Metspalu, C. Stringer, V. Macaulay, A. Cooper, and J. J. Sanchez. 2006. Multiplexed SNP typing of ancient DNA clarifies the origin of Andaman mtDNA haplogroups amongst south Asian tribal populations. PLoS One 1:e81. [PMC free article] [PubMed]
17. Filippini, S., A. Blanco, A. Fernandez-Marmiesse, V. Alvarez-Iglesias, C. Ruiz-Ponte, A. Carracedo, and A. Vega. 2007. Multiplex SNaPshot for detection of BRCA1/2 common mutations in Spanish and Spanish related breast/ovarian cancer families. BMC Med. Genet. 8:40. [PMC free article] [PubMed]
18. Filliol, I., A. S. Motiwala, M. Cavatore, W. Qi, M. H. Hazbon, M. Bobadilla del Valle, J. Fyfe, L. Garcia-Garcia, N. Rastogi, C. Sola, T. Zozio, M. I. Guerrero, C. I. Leon, J. Crabtree, S. Angiuoli, K. D. Eisenach, R. Durmaz, M. L. Joloba, A. Rendon, J. Sifuentes-Osornio, A. Ponce de Leon, M. D. Cave, R. Fleischmann, T. S. Whittam, and D. Alland. 2006. Global phylogeny of Mycobacterium tuberculosis based on single nucleotide polymorphism (SNP) analysis: insights into tuberculosis evolution, phylogenetic accuracy of other DNA fingerprinting systems, and recommendations for a minimal standard SNP set. J. Bacteriol. 188:759-772. [PMC free article] [PubMed]
19. Frothingham, R., P. L. Strickland, G. Bretzel, S. Ramaswamy, J. M. Musser, and D. L. Williams. 1999. Phenotypic and genotypic characterization of Mycobacterium africanum isolates from West Africa. J. Clin. Microbiol. 37:1921-1926. [PMC free article] [PubMed]
20. Gagneux, S., K. DeRiemer, T. Van, M. Kato-Maeda, B. C. de Jong, S. Narayanan, M. Nicol, S. Niemann, K. Kremer, M. C. Gutierrez, M. Hilty, P. C. Hopewell, and P. M. Small. 2006. Variable host-pathogen compatibility in Mycobacterium tuberculosis. Proc. Natl. Acad. Sci. U. S. A. 103:2869-2873. [PubMed]
21. Gagneux, S., and P. M. Small. 2007. Global phylogeography of Mycobacterium tuberculosis and implications for tuberculosis product development. Lancet Infect. Dis. 7:328-337. [PubMed]
22. Gilbert, M. T., J. J. Sanchez, T. Haselkorn, L. D. Jewell, S. B. Lucas, E. Van Marck, C. Borsting, N. Morling, and M. Worobey. 2007. Multiplex PCR with minisequencing as an effective high-throughput SNP typing method for formalin-fixed tissue. Electrophoresis 28:2361-2367. [PubMed]
23. Goh, K. S., M. Fabre, R. C. Huard, S. Schmid, C. Sola, and N. Rastogi. 2006. Study of the gyrB gene polymorphism as a tool to differentiate among Mycobacterium tuberculosis complex subspecies further underlines the older evolutionary age of ‘Mycobacterium canettii.’ Mol. Cell. Probes 20:182-190. [PubMed]
24. Goh, K. S., E. Legrand, C. Sola, and N. Rastogi. 2001. Rapid differentiation of “Mycobacterium canettii” from other Mycobacterium tuberculosis complex organisms by PCR-restriction analysis of the hsp65 gene. J. Clin. Microbiol. 39:3705-3708. [PMC free article] [PubMed]
25. Grignani, P., C. Turchi, A. Achilli, G. Peloso, M. Alu, U. Ricci, C. Robino, S. Pelotti, E. Carnevali, I. Boschi, A. Tagliabracci, and C. Previdere. 2009. Multiplex mtDNA coding region SNP assays for molecular dissection of haplogroups U/K and J/T. Forensic Sci. Int. Genet. 4:21-25. [PubMed]
26. Gutacker, M. M., B. Mathema, H. Soini, E. Shashkina, B. N. Kreiswirth, E. A. Graviss, and J. M. Musser. 2006. Single-nucleotide polymorphism-based population genetic analysis of Mycobacterium tuberculosis strains from 4 geographic sites. J. Infect. Dis. 193:121-128. [PubMed]
27. Gutacker, M. M., J. C. Smoot, C. A. Migliaccio, S. M. Ricklefs, S. Hua, D. V. Cousins, E. A. Graviss, E. Shashkina, B. N. Kreiswirth, and J. M. Musser. 2002. Genome-wide analysis of synonymous single nucleotide polymorphisms in Mycobacterium tuberculosis complex organisms: resolution of genetic relationships among closely related microbial strains. Genetics 162:1533-1543. [PubMed]
28. Henke, W., K. Herdel, K. Jung, D. Schnorr, and S. A. Loening. 1997. Betaine improves the PCR amplification of GC-rich DNA sequences. Nucleic Acids Res. 25:3957-3958. [PMC free article] [PubMed]
29. Huard, R. C., M. Fabre, P. de Haas, L. C. Lazzarini, D. van Soolingen, D. Cousins, and J. L. Ho. 2006. Novel genetic polymorphisms that further delineate the phylogeny of the Mycobacterium tuberculosis complex. J. Bacteriol. 188:4271-4287. [PMC free article] [PubMed]
30. Huard, R. C., L. C. Lazzarini, W. R. Butler, D. van Soolingen, and J. L. Ho. 2003. PCR-based method to differentiate the subspecies of the Mycobacterium tuberculosis complex on the basis of genomic deletions. J. Clin. Microbiol. 41:1637-1650. [PMC free article] [PubMed]
31. Hurst, C. D., T. C. Zuiverloon, C. Hafner, E. C. Zwarthoff, and M. A. Knowles. 2009. A SNaPshot assay for the rapid and simple detection of four common hotspot codon mutations in the PIK3CA gene. BMC Res. Notes 2:66. [PMC free article] [PubMed]
32. Kamerbeek, J., L. Schouls, A. Kolk, M. van Agterveld, D. van Soolingen, S. Kuijper, A. Bunschoten, H. Molhuizen, R. Shaw, M. Goyal, and J. van Embden. 1997. Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. J. Clin. Microbiol. 35:907-914. [PMC free article] [PubMed]
33. Kasai, H., T. Ezaki, and S. Harayama. 2000. Differentiation of phylogenetically related slowly growing mycobacteria by their gyrB sequences. J. Clin. Microbiol. 38:301-308. [PMC free article] [PubMed]
34. Lindblad-Toh, K., E. Winchester, M. J. Daly, D. G. Wang, J. N. Hirschhorn, J. P. Laviolette, K. Ardlie, D. E. Reich, E. Robinson, P. Sklar, N. Shah, D. Thomas, J. B. Fan, T. Gingeras, J. Warrington, N. Patil, T. J. Hudson, and E. S. Lander. 2000. Large-scale discovery and genotyping of single-nucleotide polymorphisms in the mouse. Nat. Genet. 24:381-386. [PubMed]
35. Malik, A. N., and P. Godfrey-Faussett. 2005. Effects of genetic variability of Mycobacterium tuberculosis strains on the presentation of disease. Lancet Infect. Dis. 5:174-183. [PubMed]
36. Mathema, B., N. E. Kurepina, P. J. Bifani, and B. N. Kreiswirth. 2006. Molecular epidemiology of tuberculosis: current insights. Clin. Microbiol. Rev. 19:658-685. [PMC free article] [PubMed]
37. Neonakis, I. K., Z. Gitti, E. Petinaki, S. Maraki, and D. A. Spandidos. 2007. Evaluation of the GenoType MTBC assay for differentiating 120 clinical Mycobacterium tuberculosis complex isolates. Eur. J. Clin. Microbiol. Infect. Dis. 26:151-152. [PubMed]
38. Niemann, S., D. Harmsen, S. Rusch-Gerdes, and E. Richter. 2000. Differentiation of clinical Mycobacterium tuberculosis complex isolates by gyrB DNA sequence polymorphism analysis. J. Clin. Microbiol. 38:3231-3234. [PMC free article] [PubMed]
39. Palacajornsuk, P., C. Halter, V. Isakova, M. Tarnawski, J. Farmar, M. E. Reid, and A. Chaudhuri. 2009. Detection of blood group genes using multiplex SNaPshot method. Transfusion 49:740-749. [PubMed]
40. Parsons, L. M., R. Brosch, S. T. Cole, A. Somoskovi, A. Loder, G. Bretzel, D. Van Soolingen, Y. M. Hale, and M. Salfinger. 2002. Rapid and simple approach for identification of Mycobacterium tuberculosis complex isolates by PCR-based genomic deletion analysis. J. Clin. Microbiol. 40:2339-2345. [PMC free article] [PubMed]
41. Reed, M. B., V. K. Pichler, F. McIntosh, A. Mattia, A. Fallow, S. Masala, P. Domenech, A. Zwerling, L. Thibert, D. Menzies, K. Schwartzman, and M. A. Behr. 2009. Major Mycobacterium tuberculosis lineages associate with patient country of origin. J. Clin. Microbiol. 47:1119-1128. [PMC free article] [PubMed]
42. Richter, E., M. Weizenegger, S. Rusch-Gerdes, and S. Niemann. 2003. Evaluation of genotype MTBC assay for differentiation of clinical Mycobacterium tuberculosis complex isolates. J. Clin. Microbiol. 41:2672-2675. [PMC free article] [PubMed]
43. Roberts, G. D., E. W. Koneman, and Y. K. Kim. 1991. Mycobacterium, p. 304-339. In A. Balows, W. J. Hausler, Jr., K. L. Herrmann, H. D. Isenberg, and H. J. Shadomy (ed.), Manual of clinical microbiology, 5th ed. American Society for Microbiology, Washington, DC.
44. Sanchez, J. J., C. Borsting, K. Balogh, B. Berger, M. Bogus, J. M. Butler, A. Carracedo, D. S. Court, L. A. Dixon, B. Filipovic, M. Fondevila, P. Gill, C. D. Harrison, C. Hohoff, R. Huel, B. Ludes, W. Parson, T. J. Parsons, E. Petkovski, C. Phillips, H. Schmitter, P. M. Schneider, P. M. Vallone, and N. Morling. 2008. Forensic typing of autosomal SNPs with a 29 SNP-multiplex—results of a collaborative EDNAP exercise. Forensic Sci. Int. Genet. 2:176-183. [PubMed]
45. Sanchez, J. J., C. Borsting, and N. Morling. 2005. Typing of Y chromosome SNPs with multiplex PCR methods. Methods Mol. Biol. 297:209-228. [PubMed]
46. Sanchez, J. J., and P. Endicott. 2006. Developing multiplexed SNP assays with special reference to degraded DNA templates. Nat. Protoc. 1:1370-1378. [PubMed]
47. Sanchez, J. J., C. Phillips, C. Borsting, K. Balogh, M. Bogus, M. Fondevila, C. D. Harrison, E. Musgrave-Brown, A. Salas, D. Syndercombe-Court, P. M. Schneider, A. Carracedo, and N. Morling. 2006. A multiplex assay with 52 single nucleotide polymorphisms for human identification. Electrophoresis 27:1713-1724. [PubMed]
48. Scorpio, A., and Y. Zhang. 1996. Mutations in pncA, a gene encoding pyrazinamidase/nicotinamidase, cause resistance to the antituberculous drug pyrazinamide in tubercle bacillus. Nat. Med. 2:662-667. [PubMed]
49. Sobrino, B., M. Brion, and A. Carracedo. 2005. SNPs in forensic genetics: a review on SNP typing methodologies. Forensic Sci. Int. 154:181-194. [PubMed]
50. Somoskovi, A., J. Dormandy, J. Rivenburg, M. Pedrosa, M. McBride, and M. Salfinger. 2008. Direct comparison of the genotype MTBC and genomic deletion assays in terms of ability to distinguish between members of the Mycobacterium tuberculosis complex in clinical isolates and in clinical specimens. J. Clin. Microbiol. 46:1854-1857. [PMC free article] [PubMed]
51. Sreevatsan, S., X. Pan, K. E. Stockbauer, N. D. Connell, B. N. Kreiswirth, T. S. Whittam, and J. M. Musser. 1997. Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination. Proc. Natl. Acad. Sci. U. S. A. 94:9869-9874. [PubMed]
52. Supply, P., C. Allix, S. Lesjean, M. Cardoso-Oelemann, S. Rusch-Gerdes, E. Willery, E. Savine, P. de Haas, H. van Deutekom, S. Roring, P. Bifani, N. Kurepina, B. Kreiswirth, C. Sola, N. Rastogi, V. Vatin, M. C. Gutierrez, M. Fauville, S. Niemann, R. Skuce, K. Kremer, C. Locht, and D. van Soolingen. 2006. Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J. Clin. Microbiol. 44:4498-4510. [PMC free article] [PubMed]
53. Thwaites, G., M. Caws, T. T. Chau, A. D'Sa, N. T. Lan, M. N. Huyen, S. Gagneux, P. T. Anh, D. Q. Tho, E. Torok, N. T. Nhu, N. T. Duyen, P. M. Duy, J. Richenberg, C. Simmons, T. T. Hien, and J. Farrar. 2008. Relationship between Mycobacterium tuberculosis genotype and the clinical phenotype of pulmonary and meningeal tuberculosis. J. Clin. Microbiol. 46:1363-1368. [PMC free article] [PubMed]
54. van Embden, J. D., M. D. Cave, J. T. Crawford, J. W. Dale, K. D. Eisenach, B. Gicquel, P. Hermans, C. Martin, R. McAdam, T. M. Shinnick, et al. 1993. Strain identification of Mycobacterium tuberculosis by DNA fingerprinting: recommendations for a standardized methodology. J. Clin. Microbiol. 31:406-409. [PMC free article] [PubMed]
55. Wang, H., J. Yue, M. Han, J. Yang, and Y. Zhao. 2010. Rapid method for identification of six common species of mycobacteria based on multiplex SNP analysis. J. Clin. Microbiol. 48:247-250. [PMC free article] [PubMed]
56. Warren, R. M., N. C. Gey van Pittius, M. Barnard, A. Hesseling, E. Engelke, M. de Kock, M. C. Gutierrez, G. K. Chege, T. C. Victor, E. G. Hoal, and P. D. van Helden. 2006. Differentiation of Mycobacterium tuberculosis complex by PCR amplification of genomic regions of difference. Int. J. Tuber. Lung Dis. 10:818-822. [PubMed]

Articles from Journal of Clinical Microbiology are provided here courtesy of American Society for Microbiology (ASM)