Automated high-throughput genotyping for study of global epidemiology of Mycobacterium tuberculosis based on mycobacterial interspersed repetitive units

High rates of clustering of strains causing tuberculosis in Harare, Zimbabwe: a molecular epidemiological study

We examined the pattern of tuberculosis (TB) transmission (i.e., reactivation versus recent transmission) and the impact of human immunodeficiency virus (HIV) infection in Harare, Zimbabwe. Consecutive adult smear-positive pulmonary TB patients presenting to an urban hospital in Harare were enrolled. A detailed epidemiological questionnaire was completed, and tests for HIV type 1 and CD4 cell counts were performed for each patient. Molecular fingerprinting of the genomic DNA recovered from cultures of sputum was performed by two molecular typing methods: spacer oligonucleotide typing (spoligotyping) and analysis of variable number of tandem DNA repeats (VNTRs). A cluster was defined as isolates from two or more patients that shared the same spoligotype pattern or the same VNTR pattern, or both. DNA suitable for typing was recovered from 224 patients. The prevalence of HIV infection was 79%. Of 187 patient isolates (78.6%) typed by both spoligotyping and analysis of VNTRs, 147 were identified as part of a cluster by both methods. By spoligotyping alone, 84.1% of patient isolates were grouped into 20 clusters. The cluster size was generally <8 patient isolates, although three large clusters comprised 68, 25, and 23 patient isolates. A total of 89.4% of the patient isolates grouped into 12 clusters defined by analysis of VNTRs, with 2 large clusters consisting of 127 and 13 patient isolates, respectively. Thirty-six percent of patient isolates with a shared spoligotype and 17% with a shared VNTR pattern were geographically linked within Harare, but they were not linked on the basis of the patient's home district. In a multivariate analysis, there were no independent predictors of clustering, including HIV infection status. Comparison with the International Spoligotype database (Pasteur Institute, Pointe a Pitre, Guadeloupe) demonstrated that our three largest spoligotype clusters are well recognized and ubiquitous in Africa. In this epidemiologically well characterized urban population with a high prevalence of HIV infection, we identified a very high level of strain clustering, indicating substantial ongoing recent TB transmission. Geographic linkage could be detected in a proportion of these clusters. A small group of actively circulating strains accounted for most of the cases of TB transmission.

Automated high-throughput mycobacterial interspersed repetitive unit typing of Mycobacterium tuberculosis strains by a combination of PCR and nondenaturing high-performance liquid chromatography

Mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) typing of Mycobacterium tuberculosis complex isolates is portable, 100% reproducible, and highly discriminatory. Nondenaturing high-performance liquid chromatography (non-dHPLC) with use of a WAVE microbial analysis system is a promising method of PCR amplicon analysis as it is low cost and requires no preanalysis processing. The aims of this study were to validate the application of WAVE microbial analysis system technology to MIRU-VNTR typing. A collection of 70 strains were cultivated in liquid culture and extracted using the QIAamp DNA minikit. Novel primers were designed to target the 12 MIRU-VNTR loci (P. Supply et al., J. Clin. Microbiol. 39:3563-3571, 2001). After amplification, each PCR product was analyzed on a WAVE microbial analysis system. The fragment size was calculated from the chromatogram, and the number of tandem repeats at each locus was determined. For the collection of 70 strains 100% concordance was achieved when comparing MIRU-VNTR profiles obtained from agarose gel electrophoresis and PCRs analyzed on a WAVE microbial analysis system. The calculated fragment sizes, obtained from the WAVE microbial analysis system, were sufficiently accurate to ensure 100% confidence when assigning the number of tandem repeats to a MIRU-VNTR locus. This study is the first to report the successful use of non-dHPLC for screening for variations in the number of MIRU-VNTRs in mycobacterial DNA. Non-dHPLC analysis was demonstrated to be a rapid, low-labor input method for the detection and analysis of MIRU-VNTR amplicons. The combination with non-dHPLC further enhances the utility of MIRU-VNTR typing.

Data mining of Mycobacterium tuberculosis complex genotyping results using mycobacterial interspersed repetitive units validates the clonal structure of spoligotyping-defined families

Recently, a combination of spoligotyping and bioinformatics was proposed as a potential tool for defining major circulating clades of tuberculosis bacilli. In the present study, we attempted to validate the above mentioned classification using a new high-throughput marker, named mycobacterial interspersed repetitive units (MIRUs). Using 12 MIRU loci and spoligotyping, we performed data mining of results on clinical isolates of the Mycobacterium tuberculosis complex representative of global mycobacterial allelic diversity. Knowledge rules permitting automatic labeling of major M. tuberculosis families were defined. Using this strategy, MIRU 24 appeared to be most appropriate for classifying our dataset. The Bovis family was shown to be perfectly classified by a maximum of 3 MIRUs, followed by Africanum and East African Indian (EAI) families by 4 MIRUs, the Beijing family by 6 MIRUs, Haarlem and X families by 8 MIRUs, the T family by 9, and the Latin-American and Mediterranean (LAM) family by 10 MIRUs. Considering the hierarchy of family divergence, our results corroborate a recent suggestion that EAI is the ancestral family followed by Africanum and Bovis. On the other hand, T, X, LAM and Haarlem families appear to be of more recent evolution. These results indicate that data mining of MIRUs is a valuable new tool for analyzing the evolutionary dynamics of the M. tuberculosis complex, and for monitoring an infectious disease such as tuberculosis.

Characterization of a trinucleotide repeat sequence (CGG)5 and potential use in restriction fragment length polymorphism typing of Mycobacterium tuberculosis

The genomes of 28 bacterial strains, including mycobacterial species Mycobacterium tuberculosis and Mycobacterium bovis, were analyzed for the presence of a special class of microsatellite, that of trinucleotide repeat sequences (TRS). Results of a search of all 10 possible TRS motifs (i.e., CCT, CGG, CTG, GAA, GAT, GTA, GTC, GTG, GTT, and TAT) with five or more repeating units showed that (CGG)(5) was highly represented within the genomic DNA of M. tuberculosis and M. bovis. Most of the (CGG)(5) repeats in the genome were within the open reading frames of two large gene families encoding PE_PGRS and PPE proteins that have the motifs Pro-Glu (PE) and Pro-Pro-Glu (PPE). (CGG)(5)-probed Southern hybridization showed that some mycobacterial species, such as Mycobacterium marinum, Mycobacterium kansasii, and Mycobacterium szulgai, possess many copies of (CGG)(5) in their genomes. Analysis of clinical isolates obtained from Tokyo and Warsaw with both IS6110 and (CGG)(5) probes showed that there is an association between the fingerprinting patterns and the geographic origin of the isolates and that (CGG)(5) fingerprinting patterns were relatively more stable than IS6110 patterns. The (CGG)(5) repeat is a unique sequence for some mycobacterial species, and (CGG)(5) fingerprinting can be used as an epidemiologic method for these species as well as IS6110 fingerprinting can. If these two fingerprinting methods are used together, the precise analysis of M. tuberculosis isolates will be accomplished. (CGG)(5)-based fingerprinting is particularly useful for M. tuberculosis isolates with few or no insertion elements and for the identification of other mycobacterial species when informative probes are lacking.

Use of mycobacterial interspersed repetitive unit-variable-number tandem repeat typing to examine genetic diversity of Mycobacterium tuberculosis in Singapore

Strain typing using variable-number tandem repeats of mycobacterial interspersed repetitive units (MIRU-VNTR) is a powerful tool for studying the epidemiology and genetic relationships of Mycobacterium tuberculosis isolates. For this study, isolates from 291 patients in Singapore were genotyped by this method. One hundred sixty-six distinct MIRU-VNTR patterns were detected. One hundred sixty-two strains were grouped into 1 of 35 different MIRU-VNTR clusters and 131 isolates were unique. In this sample collection, 9 of the 12 MIRU-VNTR loci were moderately or highly discriminative according to their allelic diversities. The Hunter-Gaston discriminatory index was 0.975, indicating the high power of discrimination of MIRU-VNTR typing. By direct comparisons with previously typed MIRU-VNTR patterns and by genetic relationship analyses, we could identify and clearly define four epidemic groups of M. tuberculosis in our sample, corresponding to the W/Beijing, East-Africa-Indian, Haarlem, and Delhi genotype families. Furthermore, MIRU-VNTR typing was able to clearly distinguish ancestral and modern M. tuberculosis strains as defined by TbD1 genomic deletion analysis. These results indicate that MIRU-VNTR typing can be a useful first-line tool for studying the genetic diversity of M. tuberculosis isolates in a large urban setting such as Singapore.

Utility of fast mycobacterial interspersed repetitive unit-variable number tandem repeat genotyping in clinical mycobacteriological analysis

BACKGROUND

Analysis of variable numbers of tandem repeats (VNTR) of genetic elements called mycobacterial interspersed repetitive units (MIRUs) is a recently described, polymerase chain reaction (PCR)-based method used to genotype Mycobacterium tuberculosis. It is much faster, requires a smaller amount of DNA, and has approximately the same discriminatory power as the standard IS6110 restriction fragment-length polymorphism (RFLP) method. We report the adaptation and optimization of MIRU-VNTR genotyping on a capillary electrophoresis system. We describe its application to 3 typical clinical situations encountered in our laboratory (Institut Pasteur de Bruxelles, Laboratoire Tuberculose et Mycobacteries; Brussels, Belgium).

METHODS

MIRU-VNTR genotyping was performed on heat-inactivated M. tuberculosis cultures obtained from clinical specimens on Lowenstein solid medium or in mycobacteria growth indicator liquid tubes (Becton Dickinson). After amplification of 12 genomic loci using 4 different multiplex PCRs, DNA fragments were separated by capillary electrophoresis using the ABI Prism 3100-Avant Genetic Analyzer (Applied Biosystems). Sizing of the PCR fragments and assignment of the various MIRU-VNTR alleles were done using the GeneScan and customized Genotyper software packages (PE Applied Biosystem).

RESULTS

Clustering on the basis of IS6110 fingerprinting of isolates from 3 different patients attending the same hospital was confirmed by MIRU-VNTR typing. This concordance between 2 independent, highly discriminatory techniques was decisive in triggering an epidemiological inquiry that led to identification of a bronchoscopy-related tuberculosis nosocomial infection. A mixed tuberculosis infection in a patient whose infection was initially suspected as a result of the IS6110 RFLP method was clearly identified by MIRU-VNTR typing. Finally, automated MIRU-VNTR analysis permitted the identification of laboratory contamination in 6 liquid cultures of M. tuberculosis within several hours.

CONCLUSION

These examples illustrate the utility of this genotyping technique for quick and accurate resolution of problems commonly encountered in clinical mycobacteriology.

Analysis of the allelic diversity of the mycobacterial interspersed repetitive units in Mycobacterium tuberculosis strains of the Beijing family: practical implications and evolutionary considerations

A study set comprised 44 Mycobacterium tuberculosis strains of the Beijing family selected for their representativeness among those previously characterized by IS6110-RFLP and spoligotyping (Northwest Russia, 1997 to 2003). In the present study, these strains were subjected to mycobacterial interspersed repetitive units (MIRU) typing to assess a discriminatory power of the 12-MIRU-loci scheme (P. Supply et al., J. Clin. Microbiol. 39:3563-3571, 2001). The 44 Russian Beijing strains were subdivided into 12 MIRU types with identical profiles: 10 unique strains and two major types shared by 10 and 24 strains. Thus, basically, two distinct sublineages appear to shape the evolution of the Beijing strains in Russia. Most of the MIRU loci were found to be (almost) monomorphic in the Russian Beijing strains; the Hunter-Gaston discriminatory index (HGDI) for all 12 loci taken together was 0.65, whereas MIRU26 (the most variable in our study) showed a moderate level of discrimination (0.49). The results were compared against all available published MIRU profiles of Beijing strains from Russia (3 strains) and other geographic areas (51 strains in total), including South Africa (38 strains), East Asia (7 strains), and the United States (4 strains). A UPGMA (unweighted pair-group method with arithmetic averages)-based tree was constructed. Interestingly, no MIRU types were shared by Russian and South African strains (the two largest samples in this analysis), whereas both major Russian types included also isolates from other locations (United States and/or East Asia). This implies the evolution of the Beijing genotype to be generally strictly clonal, although a possibility of a convergent evolution of the MIRU loci cannot be excluded. We propose a dissemination of the prevailing local Beijing clones to have started earlier in South Africa rather than in Russia since more monomorphic loci were identified in Russian samples than in South African samples (mean HGDI scores, 0.08 versus 0.17). To conclude, we suggest to use a limited number of MIRUs for preliminary subdivision of Beijing strains in Russian (loci 26 + 31), South African (10 + 26 + 39), and global settings (10 + 26 + 39).

Genotypic variation and stability of four variable-number tandem repeats and their suitability for discriminating strains of Mycobacterium leprae

It has not been possible to distinguish different strains of Mycobacterium leprae according to their genetic sequence. However, the genome contains several variable-number tandem repeats (VNTR), which have been used effectively in strain typing of other bacteria. To determine their suitability for differentiating M. leprae, we developed PCR systems to amplify 5 different VNTR loci and examined a battery of 12 M. leprae strains derived from patients in different regions of the United States, Brazil, Mexico, and the Philippines, as well as from wild armadillos and a sooty mangabey monkey. We found diversity at four VNTR (D = 0.74), but one system (C(16)G(8)) failed to yield reproducible results. Alleles for the GAA VNTR varied in length from 10 to 16 copies, those for AT(17) varied in length from 10 to 15 copies, those for GTA varied in length from 9 to 12 copies, and those for TA(18) varied in length from 13 to 20 copies. Relatively little variation was seen with interspecies transfer of bacilli or during short-term passage of strains in nude mice or armadillos. The TA(18) locus was more polymorphic than other VNTR, and genotypic variation was more common after long-term expansion in armadillos. Most strain genotypes remained fairly stable in passage, but strain Thai-53 showed remarkable variability. Statistical cluster analysis segregated strains and passage samples appropriately but did not reveal any particular genotype associable with different regions or hosts of origin. VNTR polymorphisms can be used effectively to discriminate M. leprae strains. Inclusion of additional loci and other elements will likely lead to a robust typing system that can be used in community-based epidemiological studies and select clinical applications.

Multiple polymorphic loci for molecular typing of strains of Mycobacterium leprae

The need for molecular tools for the differentiation of isolates of Mycobacterium leprae, the organism that causes leprosy, is urgent in view of the continuing high levels of new case detection, despite years of aggressive chemotherapy and the consequent reduction in the prevalence of leprosy. The slow onset of leprosy and the reliance on physical examination for detection of disease have restricted the epidemiological tracking necessary to understand and control transmission. Two genetic loci in several isolates of M. leprae have previously been demonstrated to contain variable-number tandem repeats (VNTRs). On the basis of these reports and the availability of the full genome sequence, multiple-locus VNTR analysis for strain typing has been undertaken. A panel of 11 short tandem repeat (STR) loci with repeat units of 1, 2, 3, 6, 12, 18, 21, and 27 bp from four clinical isolates of M. leprae propagated in armadillo hosts were screened by PCR. Fragment length polymorphisms were detected at 9 of the 11 loci by agarose gel electrophoresis. Sequencing of representative DNA products confirmed the presence of VNTRs between isolates. The application of nine new polymorphic STRs in conjunction with automated methods for electrophoresis and size determination allows greater discrimination between isolates of M. leprae and enhances the potential of this technique to track the transmission of leprosy.

Genotypic analysis of Mycobacterium tuberculosis in Bangladesh and prevalence of the Beijing strain

Genotypic analysis was performed on 48 Mycobacterium tuberculosis complex strains collected from a hospital in Dhaka city. Deletion analysis showed that the isolates were all M. tuberculosis; 13 of them were found to be of the "ancestral" type, while 35 were of the "modern" type, indicating that both endemic (ancestral type) and epidemic (modern type) strains cause tuberculosis in Bangladesh. Genotyping based on the spoligotype and variable-number tandem repeats (VNTR) of mycobacterial interspersed repetitive units (MIRU) was also done. A total of 34 strains (71%) were grouped by spoligotyping into nine different clusters; the largest comprised 15 isolates of the Beijing genotype, whereas the remaining eight clusters consisted of two to five isolates. MIRU-VNTR typing detected 32 different patterns among 44 tested strains, and the 15 Beijing strains were further discriminated by MIRU-VNTR typing (7 distinct patterns for the 15 isolates). These results indicate that MIRU-VNTR typing, along with spoligotyping and deletion analysis, can be used effectively for molecular epidemiological studies to determine ongoing transmission clusters; to our knowledge, this is the first report about the type of strains prevailing in Bangladesh.

A new high-throughput AFLP approach for identification of new genetic polymorphism in the genome of the clonal microorganism Mycobacterium tuberculosis

We have here applied high-throughput amplified fragment length polymorphism (htAFLP) analysis to strains belonging to the five classical species of the Mycobacterium tuberculosis complex. Using 20 strains, three enzyme combinations and eight selective amplification primer pairs, 24 AFLP reactions were performed per strain. Overall, this resulted in 480 DNA fingerprints and more than 1200 htAFLP-amplified PCR fragments were visualised per strain. The cumulative dendrogram correctly clustered strains from the various species, albeit within a distance of 6.5% for most of them. The single isolate of Mycobacterium canettii presented separately at 19% distance. All over, 169 fragments (14%) appeared to be polymorphic. Sixty-eight were specific for M. canetti and forty-five for Mycobacterium bovis. For the 10 different M. tuberculosis strains included in the present analysis, 56 polymorphic markers were identified. Upon sequencing 20 of these marker regions and comparisons with the H37Rv genome sequence, 25% appeared to share homology to members of the antigenically variable PE/PPE surface protein encoding gene family confirming previous findings on the genetic heterogeneity within these genes. In addition, homologues for phage genes and insertion element-encoded genes were detected. Forty-five percent of the sequences derived from ORFs with a currently unknown function, which was corroborated by genome sequence comparison for the clinical M. tuberculosis CD 1551 isolate. Sequence variation in M. tuberculosis was assessed in more detail for a subset of these loci by newly designed PCR restriction fragment length polymorphism (RFLP) tests and direct sequencing. Fourteen novel PCR RFLP tests were developed and twelve novel single nucleotide polymorphisms (SNPs) were identified, all suited for epidemiological analysis of M. tuberculosis. The tests allowed for identification of the major Mycobacterium species and M. tuberculosis variants and clones.

Genetic diversity and population structure of Mycobacterium tuberculosis in Casablanca, a Moroccan city with high incidence of tuberculosis

Although lower-resource countries have by far the highest burden of tuberculosis, knowledge of Mycobacterium tuberculosis population structure and genetic diversity in these regions remains almost nonexistent. In this paper, 150 Moroccan M. tuberculosis isolates circulating in Casablanca were genotyped by random amplified polymorphic DNA analysis using 10 different primers and by mycobacterial interspersed repetitive units-variable number of tandem repeats typing at 12 loci. The population genetic tests revealed a basically clonal structure for this population, without excluding rare genetic exchanges. Genetic analysis also showed a notable genetic polymorphism for the species M. tuberculosis, a weak cluster individualization, and an unexpected genetic diversity for a population in such a high-incidence community. Phylogenetic analyses of this Moroccan sample also supported that these isolates are genetically heterogeneous.

Molecular differentiation of Mycobacterium bovis isolates. Review of main techniques and applications

Until recently, none of the Mycobacterium bovis typing techniques permitted a satisfactory differentiation of isolates. During the last 10 years, the genome of pathogenic mycobacteria has been extensively studied, and phylogenetic analyses have shown that all (except Mycobacterium avium) belong to a single genetic species: the Mycobacterium tuberculosis complex. This increase in knowledge about the genome of these bacteria has lead to the discovery of molecular markers that allow us to differentiate isolates. Because of the phylogenetic proximity of the strains, even if most of these markers have been discovered in M. tuberculosis, they could be successfully adapted to the other bacteria of the M. tuberculosis complex, especially M. bovis. The most common markers in use today are the IS6110 insertion sequence, the direct repeat (DR) region, the poly(GC) rich (PGRS) sequences and the variable number tandem repeats (VNTR) sequences. The corresponding typing techniques are briefly described, and current knowledge of polymorphism and marker stability is detailed. If molecular markers are to offer wide perspectives for field studies, these two characteristics (polymorphism and stability) must be taken into account when choosing the marker(s) used in a study. In this context, examples of the application of molecular typing techniques for M. bovis are reviewed, on the one hand with epidemiological studies for which the major problem is the comparison between isolates and, on the other, with more general studies about the population genetics of M. bovis in a given country, and about its history and its phylogeny.

Identification of polymorphic tandem repeats by direct comparison of genome sequence from different bacterial strains: a web-based resource

Background

Polymorphic tandem repeat typing is a new generic technology which has been proved to be very efficient for bacterial pathogens such as B. anthracis, M. tuberculosis, P. aeruginosa, L. pneumophila, Y. pestis. The previously developed tandem repeats database takes advantage of the release of genome sequence data for a growing number of bacteria to facilitate the identification of tandem repeats. The development of an assay then requires the evaluation of tandem repeat polymorphism on well-selected sets of isolates. In the case of major human pathogens, such as S. aureus, more than one strain is being sequenced, so that tandem repeats most likely to be polymorphic can now be selected in silico based on genome sequence comparison.

Results

In addition to the previously described general Tandem Repeats Database, we have developed a tool to automatically identify tandem repeats of a different length in the genome sequence of two (or more) closely related bacterial strains. Genome comparisons are pre-computed. The results of the comparisons are parsed in a database, which can be conveniently queried over the internet according to criteria of practical value, including repeat unit length, predicted size difference, etc. Comparisons are available for 16 bacterial species, and the orthopox viruses, including the variola virus and three of its close neighbors.

Conclusions

We are presenting an internet-based resource to help develop and perform tandem repeats based bacterial strain typing. The tools accessible at now comprise four parts. The Tandem Repeats Database enables the identification of tandem repeats across entire genomes. The Strain Comparison Page identifies tandem repeats differing between different genome sequences from the same species. The "Blast in the Tandem Repeats Database" facilitates the search for a known tandem repeat and the prediction of amplification product sizes. The "Bacterial Genotyping Page" is a service for strain identification at the subspecies level.

Recent advances in laboratory procedures for pathogenic mycobacteria

Just as tuberculosis has persisted for many centuries as one of most serious and deadly infectious diseases in many parts of the world, so has the motivation to develop improved laboratory methods for characterizing M. tuberculosis isolates. Modern technology has lead to great improvements in mycobacteriology laboratory procedures, particularly in detection, identification, epidemiologic strain typing, and drug susceptibility testing. Although the usefulness of some of these newer methods is under evaluation, many already are showing potential as adjuncts to clinical diagnostic procedures.

The population structure of Mycobacterium bovis in Great Britain: clonal expansion

We have analyzed 11,500 isolates of Mycobacterium bovis (the cause of tuberculosis in cattle and other mammals) isolated in Great Britain (England, Wales and Scotland)] and characterized by spoligotype. Genetic exchange between cells is rare or absent in strains of the Mycobacterium tuberculosis complex so that, by using spoligotypes, it is possible to recognize "clones" with a recent common ancestor. The distribution of variable numbers of tandem repeats types in the most common clone in the data set is incompatible with random mutation and drift. The most plausible explanation is a series of "clonal expansions," and this interpretation is supported by the geographical distribution of different genotypes. We suggest that the clonal expansion of a genotype is caused either by the spread of a favorable mutation, together with all other genes present in the ancestral cell in which the mutation occurred, or by the invasion of a novel geographical region by a limited number of genotypes. A similar pattern is observed in M. tuberculosis (the main cause of tuberculosis in humans). The significance of clonal expansion in other bacteria that have recombination is discussed.

Tuberculosis: amplification-based clinical diagnostic techniques

Tuberculosis (TB) is one of the major infectious causes of morbidity and mortality worldwide. TB is difficult to control due to the time taken for the microbiological diagnosis; typically culture on solid media takes 6-8 weeks. There are number of rapid molecular methods that have been developed to diagnose new cases of tuberculosis, detect drug resistance and identify the type of mycobacteria. These assays are based on recognition of mycobacterial DNA sequences and the subsequent amplification of nucleic acid sequences to facilitate detection. This review will describe some of the molecular assays that are in use for TB diagnosis and the considerations in designing and performing such assays. Early diagnosis of tuberculosis is critical for the successful management of patients allowing informed use of chemotherapy ensuring that the right patients are treated with the right antimicrobials.

Molecular typing of Mycobacterium tuberculosis by using nine novel variable-number tandem repeats across the Beijing family and low-copy-number IS6110 isolates

Molecular epidemiological tools for genotyping clinical isolates of Mycobacterium tuberculosis have been developed and used to help track and contain transmission of tuberculosis. We identified 87 short sequence repeat loci within the genome of the M. tuberculosis H37Rv strain. Nine tandem repeats were found to be variable (variable-number tandem repeats [VNTRs]) in a set of 91 isolates. Fifty-seven of the isolates had only four IS6110 bands. The other 34 isolates were members of the Beijing strain family. The number of alleles of each these nine VNTRs was determined by examining each isolate. Six of the loci (Mtb-v1, -v4, -v10, -v15, -v18, and -v20) were able to differentiate the Beijing spoligotype identical isolates into seven distinct genotypes. Five of the loci (Mtb-v3, -v5, -v6, -v10, and -v15) were informative in discriminating the four-band IS6110 restriction fragment length polymorphism isolates from each other. The Nei's diversity values of each marker ranged from 0.02 to 0.59, with the number of alleles ranging from two to eight across the entire strain set. These nine loci provide a useful, discriminatory extension of VNTR typing methods for application to molecular epidemiologic studies of M. tuberculosis.

Mycobacterial interspersed repetitive unit typing of Mycobacterium tuberculosis compared to IS6110-based restriction fragment length polymorphism analysis for investigation of apparently clustered cases of tuberculosis

An evaluation of the utility of IS6110-based restriction fragment length polymorphism (RFLP) typing compared to a combination of variable number tandem repeat (VNTR) typing and mycobacterial interspersed repetitive unit (MIRU) typing was undertaken. A total of 53 patient isolates of Mycobacterium tuberculosis from four presumed episodes of cross-infection were examined. Genomic DNA was extracted from the isolates by a cetyl trimethylammonium bromide method. The number of copies of tandem repeats of the five loci ETR(A) to ETR(E) and 12 MIRU loci was determined by PCR amplification and agarose gel electrophoresis of the amplicons. VNTR typing identified the major clusters of strains in the three investigations in which they occurred (each representing a different evolutionary clade: 32333, 42235, and 32433). The majority of unrelated isolates (by epidemiology and RFLP typing) were also identified by VNTR typing. The concordance between the RFLP and MIRU typing was complete, with the exception of two isolates with RFLP patterns that differed by one band each from the rest of the major epidemiologically linked groups of isolates in investigation A. All of these isolates had identical MIRU and VNTR types. A further pair of isolates differed in the number of tandem repeat copies at two MIRU alleles but had identical RFLP patterns. The speed of the combined VNTR and MIRU typing approach enabled results for some of the investigations to be supplied in "real time," influencing choices in contact tracing. The ease of comparison of results of MIRU and VNTR typing, which are recorded as single multidigit numbers, was also found to greatly facilitate investigation management and the communication of results to health care professionals.

Genomic analysis of Mycobacterium tuberculosis complex strains used for production of purified protein derivative

The genomes of the tuberculin production strains Mycobacterium bovis AN5 and Mycobacterium tuberculosis DT were compared to genome-sequenced tubercle bacilli by using DNA microarrays. Neither the AN5 nor DT strain suffered extensive gene deletions during in vitro passage. This suggests that bovine tuberculin made from M. bovis AN5 is suitable to detect infection with presently prevalent M. bovis strains.

Molecular methods for Mycobacterium tuberculosis strain typing: a users guide

There are now a wide range of techniques available to type Mycobacterium tuberculosis, the problem is to choose the correct technique. For large scale epidemiological studies the portability and standardization of IS6110 restriction fragment length polymorphism (RFLP) means that this remains the gold standard technique. In the next few years the internationally standard mycobacterial interspersed repetitive unit (MIRU) may come to challenge this primacy. Low copy number stains remain a problem and these can be typed by either polymorphic Guanine cytosine-rich repetitive sequence (PGRS) or MIRU-variable numbers of tandem repeat (VNTR). To confirm whether strains are part of a true cluster PGRS remains the method of choice. For local outbreaks and investigations of laboratory cross contamination where speed is of greatest importance suspect strains should be initially investigated using a PCR-based method. The superior reproducibility and discrimination of MIRU-VNTR means that these methods should be favoured. If matches are found, then further confirmation of identity can be achieved using IS6110 RFLP or PGRS if the strains prove to have a low IS6110 copy number.

Genotyping of the Mycobacterium tuberculosis complex using MIRUs: association with VNTR and spoligotyping for molecular epidemiology and evolutionary genetics

The recent introduction of molecular methods has gained increased acceptance as a powerful tool for epidemiology and phylogeny of tuberculosis (TB). In this investigation, the efficiency of molecular typing using mycobacterial interspersed repetitive units (MIRUs) was assessed on a set of 116 Mycobacterium tuberculosis complex clinical isolates from 11 different geographic origins. The results obtained were compared with spoligotyping and variable number of tandem DNA repeats (VNTRs) typing data. Eighty-nine different MIRU profiles were obtained on the sample studied. Spoligotyping- or VNTR-defined clusters were split into subclusters by MIRU typing. Conversely, almost all of the clinical isolates clustered by MIRUs were shown to belong to spoligotyping-based defined clusters. The calculation of the discriminative power by the Hunter-Gaston index (HGI) for VNTR, spoligotyping and MIRU typing gave the values of, respectively, 0.959, 0.965 and 0.988, showing the high discriminative power of the MIRUs. The allelic diversity of the sample was calculated for each of the MIRU-VNTR loci; five MIRU loci (MIRU nos. 10, 23, 26, 31 and 40) were "highly discriminant", four (MIRU nos. 4, 16, 24 and 39) were "moderately discriminant", and three (MIRU nos. 2, 20 and 27) were "poorly discriminant". Among the three complementary VNTRs (exact tandem repeats ETR-A, ETR-B and ETR-C), ETR-A was the most discriminant locus. A combined numerical analysis of spoligotyping, VNTR and MIRU typing results partly corroborated a recently hypothesized evolutionary scenario for the M. tuberculosis complex. M. canettii would be the first branch to have diverged from a common M. tuberculosis complex ancestor. The East-African Indian (EAI) clade could be the first family to have diverged thereafter. A third branching separated a M. africanum-M. bovis clade, followed by a node separating Beijing versus non-Beijing M. tuberculosis. The Beijing clade was distinct from the Central Asian 1 (CAS1) family. Among non-Beijing strains, branches such as the Latin-American and Mediterranean (LAM), X and Haarlem clades diverged later. In conclusion, the results obtained show the congruence between clades defined by spoligotyping, and MIRU-VNTR, and underline the potential of these methods for M. tuberculosis phylogeny reconstruction. We also conclude that MIRU typing is a very promising method that may be used in a "two PCR-based" genotyping strategy, in conjunction to conventional epidemiological investigations.

Mutations in putative mutator genes of Mycobacterium tuberculosis strains of the W-Beijing family

Alterations in genes involved in the repair of DNA mutations (mut genes) result in an increased mutation frequency and better adaptability of the bacterium to stressful conditions. W-Beijing genotype strains displayed unique missense alterations in three putative mut genes, including two of the mutT type (Rv3908 and mutT2) and ogt. These polymorphisms were found to be characteristic and unique to W-Beijing phylogenetic lineage. Analysis of the mut genes in 55 representative W-Beijing isolates suggests a sequential acquisition of the mutations, elucidating a plausible pathway of the molecular evolution of this clonal family. The acquisition of mut genes may explain in part the ability of the isolates of W-Beijing type to rapidly adapt to their environment.

Modeling bacterial evolution with comparative-genome-based marker systems: application to Mycobacterium tuberculosis evolution and pathogenesis

The comparative-genomic sequencing of two Mycobacterium tuberculosis strains enabled us to identify single nucleotide polymorphism (SNP) markers for studies of evolution, pathogenesis, and epidemiology in clinical M. tuberculosis. Phylogenetic analysis using these "comparative-genome markers" (CGMs) produced a highly unusual phylogeny with a complete absence of secondary branches. To investigate CGM-based phylogenies, we devised computer models to simulate sequence evolution and calculate new phylogenies based on an SNP format. We found that CGMs represent a distinct class of phylogenetic markers that depend critically on the genetic distances between compared "reference strains." Properly distanced reference strains generate CGMs that accurately depict evolutionary relationships, distorted only by branch collapse. Improperly distanced reference strains generate CGMs that distort and reroot outgroups. Applying this understanding to the CGM-based phylogeny of M. tuberculosis, we found evidence to suggest that this species is highly clonal without detectable lateral gene exchange. We noted indications of evolutionary bottlenecks, including one at the level of the PHRI "C" strain previously associated with particular virulence characteristics. Our evidence also suggests that loss of IS6110 to fewer than seven elements per genome is uncommon. Finally, we present population-based evidence that KasA, an important component of mycolic acid biosynthesis, develops G312S polymorphisms under selective pressure.

DNA fingerprinting of Salmonella enterica subsp. enterica serovar typhimurium with emphasis on phage type DT104 based on variable number of tandem repeat loci

Seventy-eight human and environmental strains of Salmonella enterica subsp. enterica serovar Typhimurium, as well as 18 isolates of other Salmonella serovars and 6 isolates of Escherichia coli, were subjected to a novel variable number of tandem repeats (VNTR)-based fingerprinting method that showed high discrimination and reproducibility for typing serovar Typhimurium isolates. The method is based on capillary separation of PCR products from fluorescence-labeled VNTR in the serovar Typhimurium genome. The serovar Typhimurium isolates displayed 54 VNTR patterns, and the VNTR assay correctly identified strains from a well-characterized outbreak. Among 37 serovar Typhimurium phage type DT104 isolates, 28 distinct VNTR patterns were found. This VNTR-based method is fast and suitable for complete automation. Our VNTR-based method was capable of high discrimination within the homogeneous serovar Typhimurium DT104 phage type and can be used to trace outbreaks and to monitor DT104 as well as other phage types. The VNTR assay was compared to XbaI pulsed-field gel electrophoresis, amplified fragment length polymorphism analysis, integron-cassette profiles and gene PCR of intI1, qacEDelta1, sulI1, and floR. The VNTR assay showed greatly improved resolution compared to all other tested methods in this study.

Evolutionary relationships among strains of Mycobacterium tuberculosis with few copies of IS6110

Molecular typing of Mycobacterium tuberculosis by using IS6110 shows low discrimination when there are fewer than five copies of the insertion sequence. Using a collection of such isolates from a study of the epidemiology of tuberculosis in London, we have shown a substantial degree of congruence between IS6110 patterns and both spoligotype and PGRS type. This indicates that the IS6110 types mainly represent distinct families of strains rather than arising through the convergent insertion of IS6110 into favored positions. This is supported by identification of the genomic sites of the insertion of IS6110 in these strains. The combined data enable identification of the putative evolutionary relationships of these strains, comprising three lineages broadly associated with patients born in South Asia (India and Pakistan), Africa, and Europe, respectively. These lineages appear to be quite distinct from M. tuberculosis isolates with multiple copies of IS6110.

[Epidemiological study of tuberculosis in the area of Angers, France, as studied by 3 PCR-based fingerprinting methods]

The new genotyping methods efficiently complement classical epidemiological investigation in order to attempt a global approach to TB control. In the present work, we have studied the genomic diversity of Mycobacterium tuberculosis isolated during the year 1998 within the district of Angers, France (260,000 inhabitants distributed in 29 districts), in order to identify recent transmission events and any related risk factors. The methods used included "spacer oligonucleotide typing" or spoligotyping, "variable number of DNA tandem repeats" or VNTR, and "double repetitive element PCR" or DRE-PCR. The resulting spoligotyping and VNTR results were also feeded to international databases and compared with >10,000 isolates for spoligotyping and 500 isolates for VNTR, representative of about 60 countries. The results obtained underlined that most of the TB cases in our setting probably reflected reactivation cases, as clustered cases indicative of potential events of recent transmission were rare. Furthermore, interrogation of international databases showed that most of the isolates from the Angers region belonged to major conserved families of TB isolates representative of Europe, with only rare cases of Asian origin, or those previously reported in specific epidemies reported from elsewhere.

High-throughput epidemiologic typing in clinical microbiology

Mapping, and ultimately preventing, the dissemination of infectious agents is an important topic in public health. Newly developed molecular-microbiological methods have contributed significantly to recent advances in the efficient tracking of the nosocomial and environmental spread of microbial pathogens. Not only has the application of novel technologies led to improved understanding of microbial epidemiology, but the concepts of population structure and dynamics of many of the medically significant microorganisms have advanced significantly also. Currently, genetic identification of microbes is also within the reach of clinical microbiology laboratory professionals including those without specialized technology research interests. This review summarizes the possibilities for high-throughput molecular-microbiological typing in adequately equipped medical microbiology laboratories from both clinical and fundamental research perspectives. First, the development and application of methods for large-scale comparative typing of serially isolated microbial strains are discussed. The outcome of studies employing these methods allows for long-term epidemiologic surveillance of infectious diseases. Second, recent methods enable an almost nucleotide-by-nucleotide genetic comparison of smaller numbers of strains, thereby facilitating the identification of the genetic basis of, for instance, medically relevant microbiological traits. Whereas the first approach provides insights into the dynamic spread of infectious agents, the second provides insights into intragenomic dynamics and genetic functionality. The current state of technology is summarized, and future perspectives are sketched.

Characterization of Mycobacterium tuberculosis complex DNAs from Egyptian mummies by spoligotyping

Bone and soft tissue samples from 85 ancient Egyptian mummies were analyzed for the presence of ancient Mycobacterium tuberculosis complex DNA (aDNA) and further characterized by spoligotyping. The specimens were obtained from individuals from different tomb complexes in Thebes West, Upper Egypt, which were used for upper social class burials between the Middle Kingdom (since ca. 2050 BC) and the Late Period (until ca. 500 BC). A total of 25 samples provided a specific positive signal for the amplification of a 123-bp fragment of the repetitive element IS6110, indicating the presence of M. tuberculosis DNA. Further PCR-based tests for the identification of subspecies failed due to lack of specific amplification products in the historic tissue samples. Of these 25 positive specimens, 12 could be successfully characterized by spoligotyping. The spoligotyping signatures were compared to those in an international database. They all show either an M. tuberculosis or an M. africanum pattern, but none revealed an M. bovis-specific pattern. The results from a Middle Kingdom tomb (used exclusively between ca. 2050 and 1650 BC) suggest that these samples bear an M. africanum-type specific spoligotyping signature. The samples from later periods provided patterns typical for M. tuberculosis. This study clearly demonstrates that spoligotyping can be applied to historic tissue samples. In addition, our results do not support the theory that M. tuberculosis originated from the M. bovis type but, rather, suggest that human M. tuberculosis may have originated from a precursor complex probably related to M. africanum.

Linkage disequilibrium between minisatellite loci supports clonal evolution of Mycobacterium tuberculosis in a high tuberculosis incidence area

Deciphering the structure of pathogen populations is instrumental for the understanding of the epidemiology and history of infectious diseases and for their control. Although Mycobacterium tuberculosis is the most widespread infectious agent in humans, its actual population structure has remained hypothetical until now because: (i) its structural genes are poorly polymorphic; (ii) adequate samples and appropriate statistics for population genetic analysis have not been considered. To investigate this structure, we analysed the statistical associations (linkage disequilibrium) between 12 independent M. tuberculosis minisatellite-like loci by high-throughput genotyping within a model population of 209 isolates representative of the genetic diversity in an area with a very high incidence of tuberculosis. These loci contain variable number tandem repeats (VNTRs) of genetic elements named mycobacterial interspersed repetitive units (MIRUs). Highly significant linkage disequilibrium was detected among the MIRU-VNTR loci in this model. This linkage disequilibrium was also evident when the MIRU-VNTR types were compared with the IS6110 restriction fragment length polymorphism types. These results support a predominant clonal evolution of M. tuberculosis.

Genome-wide analysis of synonymous single nucleotide polymorphisms in Mycobacterium tuberculosis complex organisms: resolution of genetic relationships among closely related microbial strains

Several human pathogens (e.g., Bacillus anthracis, Yersinia pestis, Bordetella pertussis, Plasmodium falciparum, and Mycobacterium tuberculosis) have very restricted unselected allelic variation in structural genes, which hinders study of the genetic relationships among strains and strain-trait correlations. To address this problem in a representative pathogen, 432 M. tuberculosis complex strains from global sources were genotyped on the basis of 230 synonymous (silent) single nucleotide polymorphisms (sSNPs) identified by comparison of four genome sequences. Eight major clusters of related genotypes were identified in M. tuberculosis sensu stricto, including a single cluster representing organisms responsible for several large outbreaks in the United States and Asia. All M. tuberculosis sensu stricto isolates of previously unknown phylogenetic position could be rapidly and unambiguously assigned to one of the eight major clusters, thus providing a facile strategy for identifying organisms that are clonally related by descent. Common clones of M. tuberculosis sensu stricto and M. bovis are distinct, deeply branching genotypic complexes whose extant members did not emerge directly from one another in the recent past. sSNP genotyping rapidly delineates relationships among closely related strains of pathogenic microbes and allows construction of genetic frameworks for examining the distribution of biomedically relevant traits such as virulence, transmissibility, and host range.

Stability of variable-number tandem repeats of mycobacterial interspersed repetitive units from 12 loci in serial isolates of Mycobacterium tuberculosis

Variable number tandem repeats (VNTRs) of elements named mycobacterial interspersed repetitive units (MIRUs) have previously been identified in 12 minisatellite loci of the Mycobacterium tuberculosis genome. These markers allow reliable high-throughput genotyping of M. tuberculosis and represent a portable approach to global molecular epidemiology of M. tuberculosis. To assess their temporal stability, we genotyped 123 serial isolates, separated by up to 6 years and belonging to a variety of distinct IS6110 restriction fragment length polymorphism (RFLP) families, from 56 patients who had positive sputum cultures. All 12 MIRU VNTR loci were completely identical within the groups of serial isolates in 55 out of 56 groups (98.2%), although 11 pairs of isolates from the same patients with conserved MIRU VNTRs displayed slightly different IS6110 RFLP profiles. In a single case, serial isolates with an unchanged IS6110 RFLP profile showed a change in 1 out of 12 MIRU VNTR loci. These results indicate that MIRU VNTRs are stable over time and therefore are suitable for reliable follow-up of patients chronically infected with tuberculosis over long periods. Moreover, they support MIRU VNTR genotyping as a powerful first-line method followed by subtyping by IS6110 RFLP to define ongoing transmission clusters.

Improvement of differentiation and interpretability of spoligotyping for Mycobacterium tuberculosis complex isolates by introduction of new spacer oligonucleotides

The direct repeat (DR) region in Mycobacterium tuberculosis complex strains is composed of multiple well-conserved 36-bp DRs interspersed with nonrepetitive DNA spacer sequences of similar size. Clinical isolates show extensive polymorphism in this DR region, and this has led to the development of a 43-spacer reversed line blot methodology: spoligotyping. Although this method has contributed significantly to the molecular epidemiology of tuberculosis in the last decade, the discriminatory power and the readability of this method were not found to be optimal. In order to improve the discriminatory power, the usefulness of 43 redesigned oligonucleotides and the usefulness of 51 new spacer oligonucleotides were evaluated. For 314 M. tuberculosis complex strains isolated in the central part of The Netherlands over a 5-year period, 264 different IS6110 RFLP types could be distinguished, and 160 different spoligotype patterns were identified by traditional spoligotyping. After the introduction of 51 new spacer oligonucleotides, 14 additional spoligotypes were recognized. This enabled us to split 11 clusters of isolates identified by the traditional spoligotyping. Furthermore, on the basis of the new spacer oligonucleotides a dichotomy was found among the Beijing genotype isolates. Among 76 Mycobacterium bovis strains, 20 patterns were found by traditional spoligotyping and 30 patterns were found by novel probe spoligotyping, respectively. Nine M. bovis subsp. caprae isolates yielded six patterns by traditional spoligotyping and eight patterns by novel probe spoligotyping. A part of the redesigned oligonucleotides slightly improved the reading of spoligotype patterns. The reproducibility of spoligotyping, based on internal control probes, invariably yielded a high score; only 4 (1%) of the 314 patient isolates gave discrepant results. Analysis of a set of 31 duplicate M. tuberculosis complex strains demonstrated a 10% error rate for the identification of blinded duplicate samples. In a redundancy analysis, 40 essential spacer oligonucleotides of the 94-spacer sequences were selected, yielding the same number of spoligotype patterns. We propose to leave the traditional commercialized first-generation membrane for spoligotyping unchanged for current applications and to introduce a second-generation spoligotyping membrane whenever extended discrimination is required, e.g., for low-copy-number IS6110 strains or for phylogenetic studies of Beijing genotype strains.

High resolution, on-line identification of strains from the Mycobacterium tuberculosis complex based on tandem repeat typing

BACKGROUND

Currently available reference methods for the molecular epidemiology of the Mycobacterium tuberculosis complex either lack sensitivity or are still too tedious and slow for routine application. Recently, tandem repeat typing has emerged as a potential alternative. This report contributes to the development of tandem repeat typing for M. tuberculosis by summarising the existing data, developing additional markers, and setting up a freely accessible, fast, and easy to use, internet-based service for strain identification.

RESULTS

A collection of 21 VNTRs incorporating 13 previously described loci and 8 newly evaluated markers was used to genotype 90 strains from the M. tuberculosis complex (M. tuberculosis (64 strains), M. bovis (9 strains including 4 BCG representatives), M. africanum (17 strains)). Eighty-four different genotypes are defined. Clustering analysis shows that the M. africanum strains fall into three main groups, one of which is closer to the M. tuberculosis strains, and an other one is closer to the M. bovis strains. The resulting data has been made freely accessible over the internet http://bacterial-genotyping.igmors.u-psud.fr/bnserver to allow direct strain identification queries.

CONCLUSIONS

Tandem-repeat typing is a PCR-based assay which may prove to be a powerful complement to the existing epidemiological tools for the M. tuberculosis complex. The number of markers to type depends on the identification precision which is required, so that identification can be achieved quickly at low cost in terms of consumables, technical expertise and equipment.

Spoligologos: a bioinformatic approach to displaying and analyzing Mycobacterium tuberculosis data

Spacer oligonucleotide (spoligotyping) analysis is a rapid polymerase chain reaction-based method of DNA fingerprinting the Mycobacterium tuberculosis complex. We examined spoligotype data using a bioinformatic tool (sequence logo analysis) to elucidate undisclosed phylogenetic relationships and gain insights into the global dissemination of strains of tuberculosis. Logo analysis of spoligotyping data provides a simple way to describe a fingerprint signature and may be useful in categorizing unique spoligotypes patterns as they are discovered. Large databases of DNA fingerprint information, such as those from the U.S. National Tuberculosis Genotyping and Surveillance Network and the European Concerted Action on Tuberculosis, contain information on thousands of strains from diverse regions. The description of related spoligotypes has depended on exhaustive listings of the individual spoligotyping patterns. Logo analysis may become another useful graphic method of visualizing and presenting spoligotyping clusters from these databases.

Global distribution of Mycobacterium tuberculosis spoligotypes

We present a short summary of recent observations on the global distribution of the major clades of the Mycobacterium tuberculosis complex, the causative agent of tuberculosis. This global distribution was defined by data-mining of an international spoligotyping database, SpolDB3. This database contains 11708 patterns from as many clinical isolates originating from more than 90 countries. The 11708 spoligotypes were clustered into 813 shared types. A total of 1300 orphan patterns (clinical isolates showing a unique spoligotype) were also detected.

Quality assessment of Mycobacterium tuberculosis genotyping in a large laboratory network

Quality assessment exercises were conducted to evaluate the reproducibility of IS6110 DNA fingerprinting performed by eight laboratories in the National Tuberculosis Genotyping and Surveillance Network. Three panels, each with 8 to 16 isolates, were typed at all laboratories, resulting in 280 images. When the pattern obtained by the majority for each isolate was used as the standard, exact matches were obtained for 73% of patterns; 90% and 97% of patterns matched within one- and two-band differences, respectively. A second approach involved retyping of randomly selected isolates at the Centers for Disease Control and Prevention. Retyping was done for 8-19 isolates per laboratory (76 total). Paired images matched exactly for 54% of isolates and within one and two band differences, 78% and 93%, respectively. We evaluated reasons for mismatching. We also evaluated the reproducibility of spoligotyping using a test panel of 13 isolates; a discrepancy of 1 in 91 results was noted.

Genotyping analyses of tuberculosis cases in U.S.- and foreign-born Massachusetts residents

We used molecular genotyping to further understand the epidemiology and transmission patterns of tuberculosis (TB) in Massachusetts. The study population included 983 TB patients whose cases were verified by the Massachusetts Department of Public Health between July 1, 1996, and December 31, 2000, and for whom genotyping results and information on country of origin were available. Two hundred seventy-two (28%) of TB patients were in genetic clusters, and isolates from U.S-born were twice as likely to cluster as those of foreign-born (odds ratio [OR] 2.29, 95% confidence interval [CI] 1.69 to 3.12). Our results suggest that restriction fragment length polymorphism analysis has limited capacity to differentiate TB strains when the isolate contains six or fewer copies of IS6110, even with spoligotyping. Clusters of TB patients with more than six copies of IS6110 were more likely to have epidemiologic connections than were clusters of TB patients with isolates with few copies of IS6110 (OR 8.01, 95%; CI 3.45 to 18.93).

Development of variable-number tandem repeat typing of Mycobacterium bovis: comparison of results with those obtained by using existing exact tandem repeats and spoligotyping

Various genetic markers have been exploited for fingerprinting the Mycobacterium tuberculosis complex (MTBC) in molecular epidemiological studies, mainly through identifying restriction fragment length polymorphisms (RFLP). In large-scale studies, RFLP typing has practical processing and analysis limitations; therefore, attempts have been made to move towards PCR-based typing techniques. Spoligotyping (spacer oligotyping) and, more recently, variable-number tandem repeat (VNTR) typing have provided PCR-derived typing techniques. This study describes the identification and characterization of novel VNTR loci, consisting of tandem repeats in the size range of 53 to 59 bp in the MTBC, and their assessment as typing tools in 47 Mycobacterium bovis field isolates and nine MTBC strains. Spoligotyping and the previously described set of exact tandem repeats (ETRs) (R. Frothingham and W. A. Meeker-O'Connell, Microbiology 144:1189-1196, 1998) were also applied to the same panel of isolates. The allelic diversity of the individual VNTR loci was calculated, and a comparison of the novel VNTRs was made against the results obtained by spoligotyping and the existing set of ETRs. Eleven unique spoligotypes were discriminated in the panel of 47 M. bovis isolates. Greater resolution was obtained through the combination of the most-discriminating VNTRs from both sets. Considerable discrimination was achieved, with the 47 M. bovis isolates resolved into 14 unique profiles, while all nine MTBC isolates were uniquely differentiated. The novel VNTR markers described increased the discrimination possible in strain typing of M. bovis, with the added benefit of an intuitive digital nomenclature, with the allele copy number of the individual VNTRs providing a profile. VNTR typing was shown to be a valuable technique with great potential for further development and application to epidemiological tracing of tuberculosis transmissions.

Variable-number tandem repeat typing of Mycobacterium tuberculosis isolates with low copy numbers of IS6110 by using mycobacterial interspersed repetitive units

A study set of 180 Mycobacterium tuberculosis and Mycobacterium bovis isolates having low copy numbers of IS6110 were genotyped using the recently introduced method based on the variable-number tandem repeats of mycobacterial interspersed repetitive units (MIRU-VNTR). The results were compared with results of the more commonly used methods, IS6110 restriction fragment length polymorphism (RFLP) and spoligotyping. The isolates were collected in Michigan from 1996 to 1999 as part of a project to genotype all isolates from new cases of tuberculosis in the state. Twelve MIRU loci were amplified, and the amplicons were analyzed by agarose gel electrophoresis to determine the copy number at each MIRU locus. MIRU-VNTR produced more distinct patterns (80 patterns) than did IS6110 RFLP (58 patterns), as would be expected in this study set. Spoligotyping identified 59 patterns. No single method defined all unique isolates, and the combination of all three typing methods generated 112 distinct patterns identifying 90 unique isolates and 90 isolates in 22 clusters. The results confirm the potential utility of MIRU-VNTR typing and show that typing with multiple methods is required to attain maximum specificity.

A new evolutionary scenario for the Mycobacterium tuberculosis complex

The distribution of 20 variable regions resulting from insertion-deletion events in the genomes of the tubercle bacilli has been evaluated in a total of 100 strains of Mycobacterium tuberculosis, Mycobacterium africanum, Mycobacterium canettii, Mycobacterium microti, and Mycobacterium bovis. This approach showed that the majority of these polymorphisms did not occur independently in the different strains of the M. tuberculosis complex but, rather, resulted from ancient, irreversible genetic events in common progenitor strains. Based on the presence or absence of an M. tuberculosis specific deletion (TbD1), M. tuberculosis strains can be divided into ancestral and "modern" strains, the latter comprising representatives of major epidemics like the Beijing, Haarlem, and African M. tuberculosis clusters. Furthermore, successive loss of DNA, reflected by region of difference 9 and other subsequent deletions, was identified for an evolutionary lineage represented by M. africanum, M. microti, and M. bovis that diverged from the progenitor of the present M. tuberculosis strains before TbD1 occurred. These findings contradict the often-presented hypothesis that M. tuberculosis, the etiological agent of human tuberculosis evolved from M. bovis, the agent of bovine disease. M. canettii and ancestral M. tuberculosis strains lack none of these deleted regions, and, therefore, seem to be direct descendants of tubercle bacilli that existed before the M. africanum-->M. bovis lineage separated from the M. tuberculosis lineage. This observation suggests that the common ancestor of the tubercle bacilli resembled M. tuberculosis or M. canettii and could well have been a human pathogen already.