Molecular typing of Mycobacterium tuberculosis based on variable number of tandem DNA repeats used alone and in association with spoligotyping

Spoligotyping of the Mycobacterium tuberculosis complex using on-Chip PCR

AIMS

The aim of this study was to develop a rapid PCR-based method for spoligotyping of Mycobacteria in the microarray format and to compare it to conventional spoligotyping by hybridization.

METHODS AND RESULTS

The method employs the on-Chip PCR technique with primers specific for 43 spacers that separate direct repeats (DRs) in the DR region of mycobacterial DNA. The primers were immobilized on gel-based microarrays, and PCR was performed directly on the chips. The PCR fluorescence images were acquired and processed using a portable fluorescence analyzer equipped with dedicated software. Analysis takes 1.5-2 hours and can be carried out on clinical samples without additional handling. The analytical sensitivity of the method was 103 copies of target DNA. The spoligotyping results of 51 samples produced by the proposed method and by conventional reverse hybridization approach were in full concordance.

CONCLUSIONS

High throughput capacity, computerized data analysis, compact equipment, and reliable results make the on-Chip PCR an attractive alternative to intra- and interspecific spoligotyping of Mycobacterium tuberculosis complex bacteria.

SIGNIFICANCE AND IMPACT OF STUDY

Fast microarray-based spoligotyping technique using on-Chip PCR was developed.

Mycobacterium bovis: From Genotyping to Genome Sequencing

Mycobacterium bovis is the main pathogen of bovine, zoonotic, and wildlife tuberculosis. Despite the existence of programs for bovine tuberculosis (bTB) control in many regions, the disease remains a challenge for the veterinary and public health sectors, especially in developing countries and in high-income nations with wildlife reservoirs. Current bTB control programs are mostly based on test-and-slaughter, movement restrictions, and post-mortem inspection measures. In certain settings, contact tracing and surveillance has benefited from M. bovis genotyping techniques. More recently, whole-genome sequencing (WGS) has become the preferential technique to inform outbreak response through contact tracing and source identification for many infectious diseases. As the cost per genome decreases, the application of WGS to bTB control programs is inevitable moving forward. However, there are technical challenges in data analyses and interpretation that hinder the implementation of M. bovis WGS as a molecular epidemiology tool. Therefore, the aim of this review is to describe M. bovis genotyping techniques and discuss current standards and challenges of the use of M. bovis WGS for transmission investigation, surveillance, and global lineages distribution. We compiled a series of associated research gaps to be explored with the ultimate goal of implementing M. bovis WGS in a standardized manner in bTB control programs.

Genotyping Mycobacterium tuberculosis isolates with few copies of IS6110: Value of additional genetic markers

PURPOSE

IS6110 restriction fragment length polymorphism (RFLP) analysis is widely used for molecular epidemiological studies of tuberculosis. Role of spoligotyping and Fluorescent Amplified Fragment Length Polymorphism (FAFLP) was studied in low-copy number IS6110 strains of Mycobacterium tuberculosis complex (Mtbc).

METHODS

The study isolates included 70 strains of Mtbc collected from different regions of India. IS6110 restriction fragment, spoligotyping and FAFLP were performed for genotypic analysis.

RESULTS

A single copy of IS6110 was found in 30% of isolates with 90.5% of them harboring characteristic 1.5-Kb IS6110 restriction fragment.IS6110RFLP identified 51 different types, FAFLP 41 types, and spoligotyping 31 types. Combination of all three techniques identified 67 different types.IS6110 RFLP analysis was found sensitive for genotyping isolates with more than one copy of IS6110 (Hunter Gaston Discriminatory Index (HGDI-1) while, neither spoligotyping (HGI-0.89) nor FAFLP (HGDI-0.92) or their combinations were as good. The discriminatory power of spoligotyping (HGDI- 0.89) in isolates with a single copy of IS6110 was higher than IS6110-RFLP.Clustering was reduced to 67% using spoligotyping and to 38% with FAFLP.

CONCLUSION

Combination of FAFLP and Spoligotyping may prove to be valuable in studying the epidemiology of M. tuberculosis strains harboring few copies of IS6110 element.

Exploring the "Latin American Mediterranean" family and the RDRio lineage in Mycobacterium tuberculosis isolates from Paraguay, Argentina and Venezuela

Background

The Latin American & Mediterranean (LAM) spoligotype family is one of the most successful genotype of Mycobacterium tuberculosis worldwide and particularly prevalent in South-America. Within this family, a sublineage named Region of Difference Rio (RD^Rio) was reported initially in Brazil and is characterized by a genomic deletion of about 26.3 kb. This lineage seems to show a specific adaptation to the Euro-Latin American population. In this context, we sought to evaluate the LAM family and the presence of the RD^Rio genotype in samples from three Latin American countries including Paraguay, Venezuela and Argentina. To detect LAM strains reliably we applied a typing scheme using spoligotyping, 12 loci MIRU-VNTR, the Ag85C¹⁰³ SNP and the regions of difference RD^Rio and RD174. IS6110-RFLP results were also used when available.

Results

Genotyping of 413 M. tuberculosis isolates from three Latin-American countries detected LAM (46%) and the ill-defined T clade (16%) as the most frequent families. The highest clustering rate was detected in the sample population from the city of Caracas in Venezuela. We observed considerable differences in the presence of the RD^Rio lineage, with high frequency in Caracas-Venezuela (55%) and low frequency in Buenos Aires-Argentina (11%) and Paraguay (10%). The molecular markers (RD174, Ag85C¹⁰³, MIRU02-MIRU40 signature) of the RD^Rio lineage were essentially confirmed. For the LAM family, the most polymorphic loci were MIRU40, MIRU31, MIRU10, MIRU26, MIRU16 and the least polymorphic MIRU24, MIRU20, MIRU04, MIRU23.

Conclusions

Our results suggest a differential adaptation of LAM-sublineages in neighboring populations and that RD^Rio strains spread regionally with different rates of distribution. The Ag85C SNP and RDs (RD174, RD^Rio) tested in this study can in fact facilitate molecular epidemiological studies of LAM strains in endemic settings and low-income countries.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1479-6) contains supplementary material, which is available to authorized users.

The Evolutionary History, Demography, and Spread of the Mycobacterium tuberculosis Complex

With the advent of next-generation sequencing technology, the genotyping of clinical Mycobacterium tuberculosis strains went through a major breakup that dramatically improved the field of molecular epidemiology but also revolutionized our deep understanding of the M. tuberculosis complex evolutionary history. The intricate paths of the pathogen and its human host are reflected by a common geographical origin in Africa and strong biogeographical associations that largely reflect the past migration waves out of Africa. This long coevolutionary history is cardinal for our understanding of the host-pathogen dynamic, including past and ongoing demographic components, strains’ genetic background, as well as the immune system genetic architecture of the host. Coalescent- and Bayesian-based analyses allowed us to reconstruct population size changes of M. tuberculosis through time, to date the most recent common ancestor and the several phylogenetic lineages. This information will ultimately help us to understand the spread of the Beijing lineage, the rise of multidrug-resistant sublineages, or the fall of others in the light of socioeconomic events, antibiotic programs, or host population densities. If we leave the present and go through the looking glass, thanks to our ability to handle small degraded molecules combined with targeted capture, paleomicrobiology covering the Pleistocene era will possibly unravel lineage replacements, dig out extinct ones, and eventually ask for major revisions of the current model.

Molecular Epidemiology of Mycobacterium bovis in Humans and Cattle

Bovine tuberculosis (bTB), caused by Mycobacterium bovis (M. bovis), is a serious re-emerging disease in both animals and humans. The evolution of the Multi- and Extensively drug-resistant M. bovis strains (MDR-TB and XDR-TB) represents a global threat to public health. Worldwide, the disease is responsible for great economic losses in the veterinary field, serious threat to the ecosystem, and about 3.1% of human TB cases, up to 16% in Tanzania. Only thorough investigation to understand the pathogen's epidemiology can help in controlling the disease and minimizing its threat. For this purpose, various tools have been developed for use in advanced molecular epidemiological studies of bTB, either alone or in combination with standard conventional epidemiological approaches. These techniques enable the analysis of the intra- and inter-species transmission dynamics of bTB. The delivered data can reveal detailed insights into the source of infection, correlations among human and bovine isolates, strain diversity and evolution, spread, geographical localization, host preference, tracing of certain virulence factors such as antibiotic resistance genes, and finally the risk factors for the maintenance and spread of M. bovis. They also allow for the determination of epidemic and endemic strains. This, in turn, has a significant diagnostic impact and helps in vaccine development for bTB eradication programs. The present review discusses many topics including the aetiology, epidemiology and importance of M. bovis, the prevalence of bTB in humans and animals in various countries, the molecular epidemiology of M. bovis, and finally applied molecular epidemiological techniques.

Phylogenetic Distribution of CRISPR-Cas Systems in Antibiotic-Resistant Pseudomonas aeruginosa

Pseudomonas aeruginosa is an antibiotic-refractory pathogen with a large genome and extensive genotypic diversity. Historically, P. aeruginosa has been a major model system for understanding the molecular mechanisms underlying type I clustered regularly interspaced short palindromic repeat (CRISPR) and CRISPR-associated protein (CRISPR-Cas)-based bacterial immune system function. However, little information on the phylogenetic distribution and potential role of these CRISPR-Cas systems in molding the P. aeruginosa accessory genome and antibiotic resistance elements is known. Computational approaches were used to identify and characterize CRISPR-Cas systems within 672 genomes, and in the process, we identified a previously unreported and putatively mobile type I-C P. aeruginosa CRISPR-Cas system. Furthermore, genomes harboring noninhibited type I-F and I-E CRISPR-Cas systems were on average ~300 kb smaller than those without a CRISPR-Cas system. In silico analysis demonstrated that the accessory genome (n = 22,036 genes) harbored the majority of identified CRISPR-Cas targets. We also assembled a global spacer library that aided the identification of difficult-to-characterize mobile genetic elements within next-generation sequencing (NGS) data and allowed CRISPR typing of a majority of P. aeruginosa strains. In summary, our analysis demonstrated that CRISPR-Cas systems play an important role in shaping the accessory genomes of globally distributed P. aeruginosa isolates.

IMPORTANCE

P. aeruginosa is both an antibiotic-refractory pathogen and an important model system for type I CRISPR-Cas bacterial immune systems. By combining the genome sequences of 672 newly and previously sequenced genomes, we were able to provide a global view of the phylogenetic distribution, conservation, and potential targets of these systems. This analysis identified a new and putatively mobile P. aeruginosa CRISPR-Cas subtype, characterized the diverse distribution of known CRISPR-inhibiting genes, and provided a potential new use for CRISPR spacer libraries in accessory genome analysis. Our data demonstrated the importance of CRISPR-Cas systems in modulating the accessory genomes of globally distributed strains while also providing substantial data for subsequent genomic and experimental studies in multiple fields. Understanding why certain genotypes of P. aeruginosa are clinically prevalent and adept at horizontally acquiring virulence and antibiotic resistance elements is of major clinical and economic importance.

Genomics and Machine Learning for Taxonomy Consensus: The Mycobacterium tuberculosis Complex Paradigm

Infra-species taxonomy is a prerequisite to compare features such as virulence in different pathogen lineages. Mycobacterium tuberculosis complex taxonomy has rapidly evolved in the last 20 years through intensive clinical isolation, advances in sequencing and in the description of fast-evolving loci (CRISPR and MIRU-VNTR). On-line tools to describe new isolates have been set up based on known diversity either on CRISPRs (also known as spoligotypes) or on MIRU-VNTR profiles. The underlying taxonomies are largely concordant but use different names and offer different depths. The objectives of this study were 1) to explicit the consensus that exists between the alternative taxonomies, and 2) to provide an on-line tool to ease classification of new isolates. Genotyping (24-VNTR, 43-spacers spoligotypes, IS6110-RFLP) was undertaken for 3,454 clinical isolates from the Netherlands (2004-2008). The resulting database was enlarged with African isolates to include most human tuberculosis diversity. Assignations were obtained using TB-Lineage, MIRU-VNTRPlus, SITVITWEB and an algorithm from Borile et al. By identifying the recurrent concordances between the alternative taxonomies, we proposed a consensus including 22 sublineages. Original and consensus assignations of the all isolates from the database were subsequently implemented into an ensemble learning approach based on Machine Learning tool Weka to derive a classification scheme. All assignations were reproduced with very good sensibilities and specificities. When applied to independent datasets, it was able to suggest new sublineages such as pseudo-Beijing. This Lineage Prediction tool, efficient on 15-MIRU, 24-VNTR and spoligotype data is available on the web interface “TBminer.” Another section of this website helps summarizing key molecular epidemiological data, easing tuberculosis surveillance. Altogether, we successfully used Machine Learning on a large dataset to set up and make available the first consensual taxonomy for human Mycobacterium tuberculosis complex. Additional developments using SNPs will help stabilizing it.

Current methods in the molecular typing of Mycobacterium tuberculosis and other mycobacteria

In the epidemiology of tuberculosis (TB) and nontuberculous mycobacterial (NTM) diseases, as in all infectious diseases, the key issue is to define the source of infection and to disclose its routes of transmission and dissemination in the environment. For this to be accomplished, the ability of discerning and tracking individual Mycobacterium strains is of critical importance. Molecular typing methods have greatly improved our understanding of the biology of mycobacteria and provide powerful tools to combat the diseases caused by these pathogens. The utility of various typing methods depends on the Mycobacterium species under investigation as well as on the research question. For tuberculosis, different methods have different roles in phylogenetic analyses and person-to-person transmission studies. In NTM diseases, most investigations involve the search for environmental sources or phylogenetic relationships. Here, too, the type of setting determines which methodology is most suitable. Within this review, we summarize currently available molecular methods for strain typing of M. tuberculosis and some NTM species, most commonly associated with human disease. For the various methods, technical practicalities as well as discriminatory power and accomplishments are reviewed.

Use of genotyping based clustering to quantify recent tuberculosis transmission in Guadeloupe during a seven years period: analysis of risk factors and access to health care

BACKGROUND

The present study aimed to characterize Mycobacterium tuberculosis population structure and to identify transmission chains and risk factors by prospective molecular typing in conjunction with conventional epidemiological investigations in the French overseas department of Guadeloupe.

METHODS

The study included all the culture-positive TB cases (1 clinical isolate per patient; n = 129) diagnosed between a seven year period (April 4th, 1999 to December 31st, 2005). Prospective molecular typing was performed using spoligotyping and VNTRs, and a subset of 44 M. tuberculosis isolates found to be clustered was retrospectively typed using 12-loci MIRUs. Data were compared using the SITVIT2 database, followed by analysis of risk factors in function of clustering of the isolates and available demographic and socioeconomic data.

RESULTS

The study sample was characterized by a majority of new cases (87.4%); a moderate proportion of drug-resistance (7.8%); a high level of immigration (51.2% foreign-born) originating from high TB/HIV incidence neighboring islands such as Haiti or Dominican Republic; lower socioeconomic conditions (70.7% of jobless, average income 824 EUR/month); and a significantly higher proportion of TB/HIV co-infected cases (38.2% vs. 8.5%; p < 0.001), and extrapulmonary disease (18.2% vs. 4.8%; p < 0.02) among migrants as compared to French patients. The study revealed an important delay in access to healthcare with a median delay of 74.5 days between the 1st symptoms and clinical suspicion of TB. Prospective molecular typing based on spoligotyping and 5-loci VNTRs showed that evolutionary recent Euro-American lineages predominated in Guadeloupe (91.5% of isolates). In conjunction with epidemiological data, it allowed to estimate a recent transmission rate of 18.6%, which was close to the rate of 16.7% estimated using retrospective 12-loci MIRU typing. Although a higher proportion of cases in older age-group were apparently linked to reactivation; univariate analysis of risk factors did not allow pinpointing specific risk factors for a patient to belong to a TB transmission group.

CONCLUSIONS

Ongoing TB transmission in the insular, low TB-incidence setting of Guadeloupe can be defined as follows: (i) a significant proportion of imported cases of the disease from neighboring islands; (ii) significantly higher TB/HIV coinfection among foreign-born cases; and, (iii) a higher proportion of cases affecting older age-group among French patients due to reactivation. This study emphasizes the need for universal typing using spoligotyping and 15-loci MIRUs in prospective studies.

Investigation on Mycobacterium tuberculosis diversity in China and the origin of the Beijing clade

Background

Investigation of the genetic diversity of Mycobacterium tuberculosis in China has shown that Beijing genotype strains play a dominant role in the tuberculosis (TB) epidemic. In order to examine the strain diversity in the whole country, and to study the evolutionary development of Beijing strains, we sought to genotype a large collection of isolates using different methods.

Methodology/Principal Findings

We applied a 15-loci VNTR typing analysis on 1,586 isolates from the Beijing municipality and 12 Chinese provinces or autonomous regions. The data was compared to that of 900 isolates from various other worldwide geographic regions outside of China. A total of 1,162/1,586 (73.2%) of the isolates, distributed into 472 VNTR types, were found to belong to the Beijing genotype family and this represented 56 to 94% of the isolates in each of the localizations. VNTR typing revealed that the majority of the non-Beijing isolates fall into two genotype families, which represented 17% of the total number of isolates, and seem largely restricted to China. A small number of East African Indian genotype strains was also observed in this collection. Ancient Beijing strains with an intact region of difference (RD) 181, as well as strains presumably resembling ancestors of the whole Beijing genotype family, were mainly found in the Guangxi autonomous region.

Conclusions/Significance

This is the largest M. tuberculosis VNTR-based genotyping study performed in China to date. The high percentage of Beijing isolates in the whole country and the presence in the South of strains representing early branching points may be an indication that the Beijing lineage originated from China, probably in the Guangxi region. Two modern lineages are shown here to represent the majority of non-Beijing Chinese isolates. The observed geographic distribution of the different lineages within China suggests that natural frontiers are major factors in their diffusion.

Rapid deletion-based subtyping system for the Manila family of Mycobacterium tuberculosis

The Manila family of Mycobacterium tuberculosis is a group of clonal isolates seen throughout the Pacific Basin. Commonly used rapid molecular typing methods often leave large groups of Manila family isolates clustered together. Here we describe a simple deletion-based PCR method that improves the discrimination for Manila family isolates, with or without the use of mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) classification, and that is both rapid and affordable. Twenty-eight Manila family isolates, classified by spoligotyping, were collected from around the Pacific Basin from 1995 to 2003 and were tested for known genomic deletions. Nine of 15 regions of difference tested were identified as potentially discriminatory, with 18 distinct patterns; of these 9, 5 were selected for optimal discrimination using 61 Manila family isolates collected from California in 2009. For this geographically limited sample, the single large cluster was reduced to 14 distinct patterns. When the isolates were tested by spoligotyping and MIRU-VNTR, the addition of deletion analysis increased the number of distinct patterns from 43 to 56. In summary, the two study groups, which together form a single group of 89 isolates by spoligotyping, were segregated into 17 subgroups by our deletion-based subtyping system.

Three-year longitudinal study of genotypes of Mycobacterium tuberculosis in a low prevalence population

OBJECTIVE

To investigate the molecular epidemiology of tuberculosis, temporal and spatial distribution of Mycobacterium tuberculosis isolates and associations between genotypes and clinical characteristics, in a low prevalence population.

METHODS

A total of 930 M. tuberculosis isolates referred to the New South Wales (NSW, Australia) Mycobacterium Reference Laboratory in 2004-2006 were characterised by mycobacterial interspersed repetitive unit (MIRU) and spacer oligonucleotide (spoligo) typing. Associations between genotypes, patient age, disease site and drug resistance were explored and the predictive power of molecular typing was analysed using Bayesian Belief Networks.

RESULTS

Among isolates from 855 NSW residents, there were 287 spoligotypes, 494 MIRU types and 643 unique spoligotype-MIRU type combinations. They formed 73 spoligotype, 104 MIRU type and 76 spoligo-MIRU clusters, most of which contained only two isolates. The majority (87.7%) of spoligotype clusters contained several MIRU profiles and 64.4% of MIRU clusters contained several spoligotypes. The three most common M. tuberculosis clades were Beijing (24.1%), East African Indian (11.8%) and Central Asian (6.5%); 6.9% and 0.7% isolates were resistant to isoniazid and rifampicin, respectively. There was no proof of association between genotype and drug resistance but isoniazid resistance increased independently over time. Given the low rates of genotype clustering, statistical analysis of genotype-phenotype associations was limited. Potential associations were not confirmed by Bayesian classifiers.

CONCLUSIONS

Spoligo and MIRU typing demonstrated low levels of M. tuberculosis clustering in NSW; temporal and spatial changes in M. tuberculosis genotypes reflected migration patterns to Australia. No analytically significant associations between M. tuberculosis genotypes and clinical phenotypes were detected.

[Strategical use of genotyping of Mycobacterium tuberculosis in tuberculosis control]

The tuberculosis situation in Portugal justifies the use of a strategy for the genotyping of Mycobacterium tuberculosis, particularly as Portugal is part of the global backdrop of human mobility, something which has a knock-on effect on the pandemic. Several international studies have placed spoligotyping and MIRU- VNTR typing as first line techniques for the molecular epidemiology of Mycobacterium tuberculosis as these techniques rely on simple technologies (PCR) and produce patterns which are easily translated into a direct interpretation numerical code. Spoligotyping has been accordingly proposed for all the isolates, while MIRU-VNTR typing should be applied to isolates with a common spoliotype. Other techniques, including IS6110-RFLP, should be reserved for use ill accordance with selected criteria. Previous studies in Portugal using spoligotyping have underlined the advantages of a strategy based on sampling consecutive patient isolates with no prior selection criteria. This allows characterisation of the M. tuberculosis population structure through monitoring the distribution of the genotypes geographically over time and within the various risk groups. On the other hand, the association of spoligotyping, MIRU-VNTF (typing and, possibly, other techniques, needs evaluating as part of bigger pictures, including identifying recent transmission situations, distinguishing between reinfection and relapse episodes and mapping the size and dynamics of disease transmission. The solution to the tuberculosis problem in Portugal implies structuring genotyping's role in tuberculosis prevention and control and its evaluation through concrete examples and results.

Molecular characterization of Mycobacterium tuberculosis isolated in the State of Parana in southern Brazil

Sequence IS6110 has been successfully used throughout the world for characterizing the Mycobacterium tuberculosis lineages. The aim of this study was to obtain data about circulating strains of M. tuberculosis in patients from the State of Parana in southern Brazil. Sixty-two clinical specimens obtained from sputum, bronchial aspirate, biopsy and urine from 62 patients clinically diagnosed with tuberculosis and admitted to the SUS-Brazil - The Brazilian Centralized Health Service System - were genotyped by the mixed-linker PCR DNA fingerprinting technique. The analysis demonstrated that the number of copies of the IS6110 sequence per isolates varied from four to 13 bands, with an average number of 8.5. From this, 93% of the isolates presented multiple copies. Isolates with no copies of the IS6110 element were not observed. The genetic analysis by UPGMA grouped the 62 isolates by similarity into three different groups: the first group contained two strains, the second was composed of 23, and the third, a more heterogeneous group, contained 37 isolates. Only two isolates (3.2%) formed a cluster; in other words, they presented a pattern of polymorphism with similarity above 95%. Such findings suggest that in the State of Parana, illness predominantly develops through reactivation of the latent infection as opposed to exogenous transmission. The methodology used (mixed-linker PCR DNA fingerprinting) allowed for 93.5% differentiation of the isolates tested, and proved to be a powerful tool for differentiation in the molecular genotyping of M. tuberculosis.

Use of cluster-graphs from spoligotyping data to study genotype similarities and a comparison of three indices to quantify recent tuberculosis transmission among culture positive cases in French Guiana during a eight year period

Background

French Guiana has the highest tuberculosis (TB) burden among all French departments, with a strong increase in the TB incidence over the last few years. It is now uncertain how best to explain this incidence. The objective of this study was to compare three different methods evaluating the extent of recent TB transmission in French Guiana.

Methods

We conducted a population-based molecular epidemiology study of tuberculosis in French Guiana based on culture-positive TB strains (1996 to 2003, n = 344) to define molecular relatedness between isolates, i.e. potential transmission events. Phylogenetic relationships were inferred by comparing two methods: a "cluster-graph" method based on spoligotyping results, and a minimum spanning tree method based on both spoligotyping and variable number of tandem DNA repeats (VNTR). Furthermore, three indices attempting to reflect the extent of recent TB transmission (RTIn, RTIn-1 and TMI) were compared.

Results

Molecular analyses showed a total amount of 120 different spoligotyping patterns and 273 clinical isolates (79.4%) that were grouped in 49 clusters. The comparison of spoligotypes from French Guiana with an international spoligotype database (SpolDB4) showed that the majority of isolates belonged to major clades of M. tuberculosis (Haarlem, 22.6%; Latin American-Mediterranean, 23.3%; and T, 32.6%). Indices designed to quantify transmission of tuberculosis gave the following values: RTIn = 0.794, RTIn-1 = 0.651, and TMI = 0.146.

Conclusion

Our data showed a high number of Mycobacterium tuberculosis clusters, suggesting a high level of recent TB transmission, nonetheless an estimation of transmission rate taking into account cluster size and mutation rate of genetic markers showed a low ongoing transmission rate (14.6%). Our results indicate an endemic mode of TB transmission in French Guiana, with both resurgence of old spatially restricted genotypes, and a significant importation of new TB genotypes by migration of TB infected persons from neighgouring high-incidence countries.

Contribution of spoligotyping to the characterization of the population structure of Mycobacterium tuberculosis isolates in Portugal

Tuberculosis is a major health problem in Portugal. To begin characterizing the population structure of Mycobacterium tuberculosis, spoligotyping was used for the systematic typing, through consecutive sampling, of patient isolates from the Amadora-Sintra area of Greater Lisbon. Distribution amongst major spoligotype families, including the Latin American Mediterranean (LAM), T, Haarlem and Beijing, was compared to that of the international spoligotype database SpolDB4 and to the European countries of traditional Portuguese immigration represented in SpolDB4. Spoligotypes from 665 isolates were analyzed and 97 shared international types (SITs) identified. In SpolDB4 Portugal is represented by part of the spoligotypes from this study explaining the reduced number of unidentified patterns. The importance of the LAM family, and especially of LAM1 and LAM9 sub-families that alone represented 38% of all the isolates in this study as compared to 8% relative to the European sub group, led us to believe that at least in this respect the population structure was closer to that of Africa and South America than to Europe. Spoligotypes characteristic of Portugal or Portuguese related settings were identified. These included SIT244 a T1 sub-family predominant in Portugal and Bangladesh, SIT64 a LAM 6 sub-family common to Portugal and Brazil, and SIT1106 a LAM 9 sub-family. These studies were the first in Portugal stressing the importance of monitoring the population structure of M. tuberculosis isolates, an important step towards gaining an understanding of tuberculosis and the dynamics of this disease.

Progression toward an improved DNA amplification-based typing technique in the study of Mycobacterium tuberculosis epidemiology

While high-copy-number IS6110-based restriction fragment length polymorphism (HCN-RFLP) is the gold standard for typing most Mycobacterium tuberculosis strains, the time taken for culturing and low throughput make it impractical for large-scale prospective typing of large numbers of isolates. The development of a new method, mycobacterial interspersed repetitive units (MIRU), a variation of the original variable-number tandem repeat (VNTR) technique, may provide a viable alternative. Panels based on the original 12-loci MIRU (12MIRU), a combination of 12MIRU and remaining ETR loci (15MIRU-VNTR), and an extended panel with an additional 10 novel regions (25VNTR) were used to study three populations with varying degrees of epidemiological data. MIRU discrimination increased with panel size and the addition of spoligotyping. Combining these two techniques enabled a reduction in the panel size from 25 to 14 loci without a significant loss in discrimination. However, 25VNTR alone or in combination with spoligotyping still possessed weaker discrimination than RFLP for high-copy-number isolates.

Efficient differentiation of Mycobacterium tuberculosis strains of the W-Beijing family from Russia using highly polymorphic VNTR loci

The W-Beijing family is a widespread Mycobacterium tuberculosis clonal lineage that frequently causes epidemic outbreaks. This family is genetically homogeneous and conserved, so ETR-VNTR (exact tandem repeat-variable number of tandem repeats) typing is insufficient for strain differentiation, due to a common ETR-A to E profile (42435). This leads to the false clustering in molecular epidemiological studies, especially in the regions of predominance of the W-Beijing family. In this study, we searched for VNTR loci with a high evolutionary rate of polymorphism in the W-Beijing genome. Here we further evaluated VNTR typing on a set of 99 Mycobacterium tuberculosis clinical isolates and reference strains. These isolates were characterized and classified into several genotype families based on three ETR loci (A, C, E) and eight additional loci [previously described as QUB (Queen's University Belfast) or MIRU (Mycobacterial Interspersed Repetitive Units) or Mtubs]. Ninety-nine strains were divided into 74 VNTR-types, 51 isolates of the W-Beijing family identified by IS6110 RFLP-typing (the restriction fragment length polymorphism-typing) and/or spoligotyping were subdivided into 30 VNTR-types. HGDI (the Hunter-Gaston discriminatory index) for all studied loci was close to that of IS6110 RFLP typing, a "gold standard" method for subtyping M. tuberculosis complex strains. The QUB 26 and QUB 18 loci located in the PPE genes were highly polymorphic and more discriminative than other loci (HGDI is 0.8). Statistically significant increase of tandem repeats number in loci ETR-A, -E, QUB 26, QUB 18, QUB 11B, Mtub21 was revealed in the W-Beijing group compared to genetically divergent non-W-Beijing strains. Thirty-six isolates were subjected to IS6110 RFLP typing. The congruence between results of the IS6110 RFLP typing and 11-loci VNTR typing was estimated on 23 isolates of the W-Beijing family. These isolates were subdivided into 9 IS6110-RFLP types and 13 VNTR types. The poor profiles correlation (0.767) reflects the differences in the rate and type of evolution between genome regions targeted by IS6110-RFLP and VNTR typing. VNTR typing in proposed format is powerful tool for discrimination of M. tuberculosis strains with different level of genetic relationship.

Épidémiologie moléculaire de la tuberculose en Guadeloupe de 1994 à 2000

In Guadeloupe, the incidence of tuberculosis decreased between 1994 and 2000. The rate of resistance to at least one antibiotic remained constant at 11%, whereas the rate of multiple-drug resistance increased from 0.9 to 2.4% in 2000. The proportion of patients of foreign origin (mainly from Haiti and the Dominican Republic) increased whereas the number of French patients decreased. These results show that the epidemiology of tuberculosis in Guadeloupe is similar to industrialized countries as older people, foreigners from countries where TB is endemic, and HIV+ patients are at a higher risk to declare tuberculosis disease. Molecular typing realized by spoligotyping showed the importance of previous successive colonizations and migrations as characterized by the presence of major phylogenetic families originating essentially from Northern Europe (Haarlem), Latin America and Mediterranean (LAM) and from Anglo-Saxon countries (X). The sub-typing of clustered strains by IS6110-RFLP and by a PCR method based on the variable number of tandem DNA repeats (VNTR), highlighted 29 clusters, corresponding to 44.8% of clustered strains, and allowed to estimate the rate of recent transmission at 32.2%. The epidemiologic data associated with fingerprinting results underlined the importance of reactivation cases among older people, a significant number of imported TB cases without evident links, and casual contacts.

Molecular diagnostics in tuberculosis

Molecular diagnostics in tuberculosis has enabled rapid detection of Mycobacterium tuberculosis complex in clinical specimens, identification of mycobacterial species, detection of drug resistance, and typing for epidemiological investigation. In the laboratory diagnosis of tuberculosis, the nucleic acid amplification (NAA) test is rapid and specific but not as sensitive as culture of mycobacteria. The primary determinant of successful NAA testing for tuberculosis depends on the shedding of mycobacterial DNA in secretions from caseating granulomas and its dissemination into sterile body fluids or tissue biopsies. In multibacillary diseases with a high mycobacterial load, a positive Ziehl-Neelsen smear with a positive NAA test is diagnostic of active tuberculosis, whereas a positive Ziehl-Neelsen smear with a negative NAA test in the absence of inhibitors would indicate nontuberculous mycobacterial disease. The role of the NAA test is more important in paucibacillary diseases with low mycobacterial loads. The presence of polymerase chain reaction (PCR) inhibitors, however, especially in extrapulmonary specimens, may produce false-negative results. Although this problem can be overcome to some extent by extra extraction steps, the additional processing invariably leads to the loss of mycobacterial DNA. To circumvent this problem, a brief culture augmentation step is carried out before the NAA test is performed, which can enhance the mycobacterial load while concomitantly diluting inhibitors, thereby maintaining the sensitivity of the test without excessively increasing turnaround time.

Molecular epidemiology of bovine tuberculosis in wild animals in Spain: A first approach to risk factor analysis

In human tuberculosis (Mycobacterium tuberculosis), molecular epidemiology has accurately indicated the risk factors involved in active transmission of the disease, by comparing individuals whose isolates belong to a cluster with patients whose strains are considered unique. Nevertheless, this application has not been used in bovine tuberculosis (Mycobacterium bovis). Our study describes the integration of epidemiological data into molecular classification data on M. bovis isolates. These were isolated from wild ungulates in Extremadura (western Spain) with the objective of detecting the risk factors linked to the association of strains in clades, which are indicators of the active spread of the disease. The molecular markers used were spoligotyping + VNTR typing (loci: VNTR 2165, VNTR 2461, VNTR 0577, VNTR 0580, VNTR 3192 VNTR 2163a and VNTR 2163b) on a population of 59 M. bovis strains isolated from deer (Cervus elaphus), 112 from wild boar (Sus scrofa), six from bovines, 28 from pigs and 2 from goats (n=207). Epidemiological variables included the animal species from which the strain was isolated, pathological condition of the host (incipient lesion, early and late generalisation), date of sampling (during or after the reproductive period) and hunting season. Bivariant analysis was used to establish the risk factors connected to the association of strains and later, the variables were evaluated by means of logistic regression. Molecular typing grouped a total of 131 strains (64.21%) in 28 clusters and 76 isolates shows unique profiles. The association of strains was connected to the appearance of macroscopic lesions during the reproductive period (O.R. 4.80; 95% CI 1.09-22.99, P<0.005), showing a possible higher transmission during the courting period. This happened mainly during the last hunting season analysed (2002-2003, O.R. 3.69; 95% CI 1.27-11.9, P<0.05), clashing with the time of higher prevalence of the disease in wild ungulates. Active spread was not connected to any species in particular, or to any concrete pathological condition.

Comparison between spoligotyping and IS6110 restriction fragment length polymorphisms in molecular genotyping analysis of Mycobacterium tuberculosis strains

Spoligotyping was compared with RFLP fingerprinting analysis in the identification of Mycobacterium tuberculosis strains. Spoligotyping sensitivity was 97.6% with a specificity of 47%. The global probability for two strains clustered with spoligotyping to be clustered also with RFLP analysis was 33%; the probability for two strains clustered with RFLP analysis to be clustered also with spoligotyping analysis was 95%. However, comparing the two methods in five outbreak episodes, full concordance was evidenced between spoligotyping and RFLP. Moreover, we evaluated the presence of our 17 largest spoligotyping clusters in spoligotyping databases from Caribbean countries, London and Cuba. Only five out of 17 patterns were present in all the cohorts. The conditional probability comparing spoligotyping and RFLP methods related to these patterns resulted in very low concordance (range from 2 to 38%). In conclusion, we confirm that spoligotyping when used alone overestimates the number of recent transmission and does not represent a suitable method for wide clinical practice application. However, it allows to get a first good picture of strain identity in a new setting and in more localized or confined settings, the probability of reaching the same result compared to RFLP was 100% confirming the usefulness of spoligotyping in the management of epidemic events, especially in hospitals, prisons and close communities.

Molecular characterization of Mycobacterium tuberculosis isolates in a region of Brazil with a high incidence of tuberculosis

One hundred and seventy Mycobacterium tuberculosis clinical isolates were characterized by spoligotyping to evaluate the biodiversity of tubercle bacilli in a region of Brazil with a high incidence of tuberculosis (Pelotas and Rio Grande cities - Rio Grande do Sul State). The spoligotyping results were compared to the World Spoligotyping Database (Institut Pasteur de Guadeloupe), which contains data from >14,000 worldwide isolates of M. tuberculosis. The isolates clustered by spoligotyping were further characterized by IS6110-RFLP to confirm the clonal relationship. Sixty-six different spoligotypes were identified, grouping 125 of the isolates (74%). Approximately half of the isolates belonged to seven of the most frequently occurring spoligotypes in the database. Three shared types (with two or more isolates) not previously identified were given the type numbers 826, 827 and 863. An additional 45 spoligotypes were identified that did not match any existing database pattern. RFLP characterization reduced the number of isolates in most of the clusters, thereby showing a higher differentiation capacity than spoligotyping. These results highlight the importance of molecular epidemiology studies of tuberculosis in insufficiently studied regions with a high TB burden, in order to uncover the true extent of genetic diversity of the pathogen.

Spoligotyping of Mycobacterium tuberculosis isolates from patients with pulmonary tuberculosis in Mumbai, India

Tuberculosis remains a major health problem in India, with 2 million new cases and 421,000 deaths each year. In this paper, we describe the spoligotyping results of 216 Mycobacterium tuberculosis culture isolates from patients with pulmonary tuberculosis in Mumbai, India. As spoligotyping data from India have rarely been described until now, and as there is limited information on the major circulating clades of M. tuberculosis, the data obtained were also compared to an international spoligotype database (SpolDB4) that contained patterns from 22,546 isolates from more than 100 countries. Eighty-four (39%) of the isolates were definitively marked as orphan strains, indicating the paucity of such data from India. The remaining 132 isolates clustered among 59 shared types; among these, 42 shared types were already present in the database, 17 were newly created, and 5 of them were specifically reported from Mumbai. A total of 9 major types in this study clustered 32% of the isolates. At the phylogenetic level, 30% of the isolates belonged to the Central Asian families CAS1 and CAS2, of the major genetic group (MGG) 1, 29% to MGG 2 and 3 families (spacers 33-36 missing) and 17% to the ancestral East African Indian (EAI) family. Finally, nearly 10% of the isolates belonged to the W-Beijing family in a broad sense, also in the MGG 1 group. In conclusion, historic clones of the MGG 1 group of M. tuberculosis are responsible for roughly 60% of all tuberculosis cases in Mumbai. Together with the fact that organisms presumably of European descent (such as the Haarlem family) were only rarely found, our observations suggest that tuberculosis in Mumbai, India is essentially caused by historical clones of tubercle bacilli undergoing active circulation due to uncontrolled demography, high prevalence of the disease, and a paucity of resources.

Mycobacterium tuberculosis molecular evolution in western Mediterranean Islands of Sicily and Sardinia

In this study, a total of 204 Mycobacterium tuberculosis DNAs from Sicily (n = 144) and Sardinia (n = 60) were studied by three genotyping methods. Results were analyzed both within and across islands, to define the phylogeographical specificities of the genotypes, look for their diversity and infer a molecular evolutionary scenario. A strong link between geography and tuberculosis genotypes was observed in Sardinia. The results were also matched against a world-wide genetic diversity database to compare the population structure of the tubercle bacilli in the islands. Eight common genotypes between Sicily, Sardinia and continental Italy were found which underlines the influences of the Italian mainland on the population structure on the islands and vice versa. A unified evolutionary scenario of M. tuberculosis evolution was built using numerical taxonomy and maximum parsimony (MP) methods. The finding of multiple families of M. tuberculosis strains (S, T, LAM, Haarlem), their presumed links with the major genetic groups (MGG) of M. tuberculosis complex, supports the view of independent introduction of several ancestral genotypes in Sicily and in Sardinia. We conclude that the two PCR-based genotyping combination (spoligotyping-VNTR) is an excellent tool to reconstruct M. tuberculosis phylogeny, that may be used to construct global and local evolutionary scenarios of the M. tuberculosis complex. The results obtained are paradigmatic of the complex interplay that exists between epidemic dynamics and evolutionary genetics of M. tuberculosis.

Characterization of ancestral Mycobacterium tuberculosis by multiple genetic markers and proposal of genotyping strategy

Sixty-eight ancestral Mycobacterium tuberculosis isolates were previously identified by using the tuberculosis-specific deletion 1 (TbD1) PCR and mycobacterial interspersed-repetitive-unit-variable-number tandem repeat (MIRU-VNTR) typing (Y. J. Sun, R. Bellamy, A. S. G. Lee, S. T. Ng, S. Ravindran, S.-Y. Wong, C. Locht, P. Supply, and N. I. Paton, J. Clin. Microbiol. 42:1986-1993, 2004). These TbD1(+) ancestral isolates were further characterized and typed in this study by IS6110 restriction fragment length polymorphism (RFLP) typing, VNTR typing using exact tandem repeats (VNTR-ETR), and spoligotyping of the direct-repeat region. To our knowledge, this is the first characterization of this genogroup by multiple genetic markers based on a fairly large sample size. In this genogroup, all spoligotypes were characterized by the absence of spacers 29 to 32 and 34. In addition, VNTR-ETR typing could add further resolution to the clustered isolates identified by MIRU-VNTR, and the combination of MIRU-VNTR and VNTR-ETR, called MIRU-ETR, showed the highest discriminatory power for these strains compared to IS6110 RFLP typing and spoligotyping alone. However, MIRU-ETR appeared to still cluster some probably epidemiologically unrelated strains, as judged by IS6110 RFLP divergence. Therefore, a typing strategy based on stepwise combination of MIRU-ETR and IS6110 RFLP is proposed to achieve maximal discrimination for unrelated TbD1(+) strains. This typing strategy may be useful in areas where TbD1(+) ancestral strains are prevalent.

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.

Transmission of mycobacterium tuberculosis via lung transplantation

Organ donors are not routinely screened for tuberculosis (TB) in the United States. We investigated a case of pulmonary TB in a double-lung transplant recipient. We reviewed the donor's and recipient's records, and used molecular methods to compare the lung recipient's isolate with others from three sources: her hospital, the California state health department's genotyping database, and the donor's resident-nation of Guatemala. A respiratory specimen obtained from the lung recipient 1 day after transplantation grew Mycobacterium tuberculosis. Donor chest radiograph had a previously unnoticed pulmonary opacity that was present on post-transplant recipient chest radiographs and computed tomographs. The recipient's isolate was molecularly distinct from others at her hospital and in the state database, but was identical to two isolates from Guatemala. Tuberculosis was transmitted from lung donor to recipient. As organ transplantation becomes more common worldwide, similar cases could occur. Screening for TB in potential organ donors should be considered.

Molecular epidemiology of Mycobacterium tuberculosis in western Sweden

The genetic diversity of Mycobacterium tuberculosis isolates among patients from Sweden was determined by a combination of two PCR-based techniques (spoligotyping and variable number of tandem repeats analysis). It resulted in a clustering of 23.6% of the isolates and a rate of recent transmission of 14.1%. The clustered isolates mainly belonged to the Haarlem family (23.2%), followed by the Beijing (9.8%), Latin American and Mediterranean (LAM; 8%), and East African-Indian (EAI; 6.2%) families. A comparison of the spoligotypes with those in the international spoligotyping database showed that 62.5% of the clustered isolates and 36.6% of all isolates typed were grouped into six major shared types. A comparison of the spoligotypes with those in databases for Scandinavian countries showed that 33% of the isolates belonged to an ill-defined T family, followed by the EAI (22%), Haarlem (20%), LAM (11%), Central Asian (5%), X (5%), and Beijing (4%) families. Both the highest number of cases and the proportion of clustered cases were observed in patients ages 15 to 39 years. Nearly 10% of the isolates were resistant to one or more drugs (essentially limited to isoniazid monoresistance). However, none of the strains were multidrug resistant. Data on the geographic origins of the patients showed that more than two-thirds of the clustered patients with tuberculosis were foreign-born individuals or refugees. These results are explained on the basis of both the historical links within specific countries and recently imported cases of tuberculosis into Sweden.

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.

Epidemiology of Mycobacterium bovis infections of pigs and wild boars using a molecular approach

A molecular epidemiological approach was applied to establishing a possible role for the wild boar as a natural reservoir of Mycobacterium bovis in Sierra de Villuercas, Western Spain; an area free of farmed cattle and wild deer populations. Spoligo and VNTR typing were used over a three year period to study the epidemiological relationship between the occurrence of bovine tuberculosis (TB) in extensively bred Iberian pigs and indigenous wild boar. The 37 sampled wild boar showed different degree of calcified granulomatous lesions in retropharyngeal, mediastinal and pulmonary lymph nodes. The 25 sampled Iberian pigs showed calcified lesions, mainly in the respiratory tract. Lesions located in the mesenteric lymph nodes appeared secondarily. M. bovis was isolated from all affected animals. Twenty-five and 37 isolates of M. bovis were obtained from domestic pigs and wild boar, respectively. Our findings provide evidence that supports the possibility of cross infection between wild boar and domestic pig populations. This is contrary to the generally held belief that swine represent an epidemiological dead end host and play no role in the epidemiology of M. bovis.

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.

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.

Evaluation and utilization as a public health tool of a national molecular epidemiological tuberculosis outbreak database within the United Kingdom from 1997 to 2001

The aim of this study was to develop a national model and analyze the value of a molecular epidemiological Mycobacterium tuberculosis DNA fingerprint-outbreak database. Incidents were investigated by the United Kingdom PHLS Mycobacterium Reference Unit (MRU) from June 1997 to December 2001, inclusive. A total of 124 incidents involving 972 tuberculosis cases, including 520 patient cultures from referred incidents and 452 patient cultures related to two population studies, were examined by using restriction fragment length polymorphism IS6110 fingerprinting and rapid epidemiological typing. Investigations were divided into the following three categories, reflecting different operational strategies: retrospective passive analysis, retrospective active analysis, and retrospective prospective analysis. The majority of incidents were in the retrospective passive analysis category, i.e., the individual submitting isolates has a suspicion they may be linked. Outbreaks were examined in schools, hospitals, farms, prisons, and public houses, and laboratory cross-contamination events and unusual clinical presentations were investigated. Retrospective active analysis involved a major outbreak centered on a high school. Contact tracing of a teenager with smear-positive pulmonary tuberculosis matched 14 individuals, including members of his class, and another 60 cases were identified in schools clinically and radiologically and by skin testing. Retrospective prospective analysis involved an outbreak of 94 isoniazid-resistant tuberculosis cases in London, United Kingdom, that began after cases were identified at one hospital in January 2000. Contact tracing and comparison with MRU databases indicated that the earliest matched case had occurred in 1995. Subsequently, the MRU changed to an active prospective analysis targeting linked isoniazid-monoresistant isolates for follow up. The patients were multiethnic, born mainly in the United Kingdom, and included professionals, individuals from the music industry, intravenous drug abusers, and prisoners.

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.

A variable number of tandem repeats result in polymorphic alpha -isopropylmalate synthase in Mycobacterium tuberculosis

A locus of variable number of the tandem repeat, VNTR4155, resides in the putative leuA gene, encoding for alpha -isopropylmalate synthase (alpha -IPMS) of Mycobacterium tuberculosis, a repeat that is unique to the bacterium. The objective was to determine whether the polymorphic VNTR4155 was translated and resulted in a polymorphic protein. The putative leuA gene of the M. tuberculosis H37Rv strain was cloned by PCR and expressed in a His-tagged form in Escherichia coli. The enzymatic properties of the purified protein were studied. The protein was used as an antigen to immunize rabbits. Soluble proteins of several strains of M. tuberculosis were examined by Western blot analysis. The polymorphism of VNTR4155 was due to the presence of different copy number of the 57-bp tandem repeat. The putative alpha -IPMS of various strains of M. tuberculosis had different sizes, varying directly with the length of their VNTR4155.

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.

Molecular characterization of multiple-drug-resistant Mycobacterium tuberculosis isolates from northwestern Russia and analysis of rifampin resistance using RNA/RNA mismatch analysis as compared to the line probe assay and sequencing of the rpoB gene

This investigation evaluated the potential of RNA/RNA mismatch analysis for the detection of rifampin resistance among 38 multiple-drug-resistant (MDR) isolates of Mycobacterium tuberculosis from northwestern Russia. The results obtained were compared with a commercialized line probe assay and rpoB sequencing, and the genetic diversity of the isolates was also investigated in parallel using spoligotyping and variable number of tandem DNA repeats (VNTR). The mismatch analysis revealed 3 distinct RNA cleavage profiles permitting the subdivision of the strains into mutation groups 1 to 3, the most common being group 1 (28 of 38 isolates) that contained a majority of strains with a TCG531>TTG (Ser->Leu) mutation, followed by group 2 (6 of 38 isolates) characterized by different mutations in the codon CAC526 (His), and group 3 (4 of 38 isolates), all characterized by a GAC516(Asp) mutation. Spoligotyping revealed the Beijing type to be the most prevalent among mismatch group 1 (24 out of 28 strains), suggesting that the most frequent rpoB mutation among the Beijing family in our setting was TCG531 >TTG (Ser->Leu). All the Beijing type isolates were also characterized by a unique VNTR pattern made up of exact tandem repeats (ETR)-A to E of 42435. We conclude that the Beijing genotype constitutes the major family of MDR-TB isolates currently circulating in northwestern Russia, and that the in-house RNA/RNA mismatch analysis may be successfully used for rapid and reliable diagnosis of rifampin-resistant tuberculosis in this setting.

Discrimination of Mycobacterium tuberculosis complex bacteria using novel VNTR-PCR targets

The lack of a convenient high-resolution strain-typing method has hampered the application of molecular epidemiology to the surveillance of bacteria of the Mycobacterium tuberculosis complex, particularly the monitoring of strains of Mycobacterium bovis. With the recent availability of genome sequences for strains of the M. tuberculosis complex, novel PCR-based M. tuberculosis-typing methods have been developed, which target the variable-number tandem repeats (VNTRs) of minisatellite-like mycobacterial interspersed repetitive units (MIRUs), or exact tandem repeats (ETRs). This paper describes the identification of seven VNTR loci in M. tuberculosis H37Rv, the copy number of which varies in other strains of the M. tuberculosis complex. Six of these VNTRs were applied to a panel of 100 different M. bovis isolates, and their discrimination and correlation with spoligotyping and an established set of ETRs were assessed. The number of alleles varied from three to seven at the novel VNTR loci, which differed markedly in their discrimination index. There was positive correlation between spoligotyping, ETR- and VNTR-typing. VNTR-PCR discriminates well between M. bovis strains. Thirty-three allele profiles were identified by the novel VNTRs, 22 for the ETRs and 29 for spoligotyping. When VNTR- and ETR-typing results were combined, a total of 51 different profiles were identified. Digital nomenclature and databasing were intuitive. VNTRs were located both in intergenic regions and annotated ORFs, including PPE (novel glycine-asparigine-rich) proteins, a proposed source of antigenic variation, where VNTRs potentially code repeating amino acid motifs. VNTR-PCR is a valuable tool for strain typing and for the study of the global molecular epidemiology of the M. tuberculosis complex. The novel VNTR targets identified in this study should additionally increase the power of this approach.

Spoligotype diversity of Mycobacterium bovis strains isolated in France from 1979 to 2000

The molecular fingerprints of 1,349 isolates of Mycobacterium bovis received between 1979 and August 2000 at Agence Française de Sécurité Sanitaire des Aliments (Afssa) have been obtained by spoligotyping. The majority of the isolates (1,266) were obtained from cattle living in France. An apparently high level of heterogeneity was observed between isolates. One hundred sixty-one spoligotypes were observed in total, of which 153 were from French isolates. The two predominant spoligotypes, designated BCG-like and GB54, accounted for 26 and 12% of the isolates, respectively. In addition, 84% of the spoligotypes were found fewer than 10 times. Analysis of the results by clustering and parsimony-based algorithms revealed that the majority of the spoligotypes were closely related. The predominant spoligotype was identical to that of the vaccine strain Mycobacterium bovis BCG, which was isolated in France at the end of the 19th century. Some spoligotypes were closely associated with restricted geographical areas. Interestingly, some spoligotypes, which were frequently observed in France, were also observed in neighboring countries. Conversely, few spoligotypes were common to France and England, and those that were shared were observed at very different frequencies. This last point illustrates the potential role for an international data bank, which could help trace the spread of M. bovis across national borders.

Molecular epidemiology of tuberculosis: recent developments and applications

The standard method for the typing of Mycobacterium tuberculosis is still IS6110 restriction fragment length polymorphism (RFLP). This method has been widely used and has provided information on the variety and distribution of tuberculosis strain types across the globe. Recently, IS6110 RFLP has been used to investigate the question of reinfection versus reactivation, examine the existence of multiple infection, and track the spread of multidrug-resistant tuberculosis. There have also been efforts to increase our understanding of the biologic characteristics of IS6110. These studies have resulted in a clearer understanding of fingerprinting data and increased our understanding of the evolution and pathogenicity of this organism.

Comparison of variable number tandem repeat and IS6110-restriction fragment length polymorphism analyses for discrimination of high- and low-copy-number IS6110 Mycobacterium tuberculosis isolates

The present study was designed to evaluate the use of variable number tandem repeat (VNTR) and IS6110-restriction fragment length polymorphism (RFLP) analyses in combination as a two-step strategy for discrimination (as measured by the Hunter-Gaston Discrimination Index [HGDI]) of both high- and low-copy-number IS6110 Mycobacterium tuberculosis isolates compared to IS6110-RFLP alone with an unselected collection of isolates. Individually, IS6110-RFLP fingerprinting produced six clusters that accounted for 69% of the low-copy-number IS6110 isolates (five clusters) and 5% of the high-copy-number IS6110 isolates (one cluster). A total of 39% of all the isolates were clustered (HGDI = 0.97). VNTR analysis generated a total of 35 different VNTR allele profile sets from 93 isolates (HGDI = 0.938). Combining IS6110-RFLP analysis with VNTR analysis reduced the overall percentage of clustered isolates to 29% (HGDI = 0.988) and discriminated a further 27% of low-copy-number isolates that would have been clustered by IS6110-RFLP alone. The use of VNTR analysis as an initial typing strategy facilitates further analysis by IS6110-RFLP, and more importantly, VNTR analysis subdivides some IS6110-RFLP-defined clusters containing low- and single-copy IS6110 isolates.

Genetic diversity of Mycobacterium tuberculosis in Sicily based on spoligotyping and variable number of tandem DNA repeats and comparison with a spoligotyping database for population-based analysis

In a previous study, we proposed to associate spoligotyping and typing with the variable number of tandem DNA repeats (VNTR) as an alternative strategy to IS6110-restriction fragment length polymorphism (RFLP) for molecular epidemiological studies on tuberculosis. The aim of the present study was to further evaluate this PCR-based typing strategy and to describe the population structure of Mycobacterium tuberculosis in another insular setting, Sicily. A collection of 106 DNA samples from M. tuberculosis patient isolates was characterized by spoligotyping and VNTR typing. All isolates were independently genotyped by the standard IS6110-RFLP method, and clustering results between the three methods were compared. The totals for the clustered isolates were, respectively, 15, 60, and 82% by IS6110-RFLP, spoligotyping, and VNTR typing. The most frequent spoligotype included type 42 that missed spacers 21 to 24 and spacers 33 to 36 and derived types 33, 213, and 273 that, together represented as much as 26% of all isolates, whereas the Haarlem clade of strains (types 47 and 50, VNTR allele 32333) accounted for 9% of the total strains. The combination of spoligotyping and VNTR typing results reduced the number of clusters to 43% but remained superior to the level of IS6110-RFLP clustering (ca. 15%). All but one IS6110-defined cluster were identified by the combination of spoligotyping and VNTR clustering results, whereas 9 of 15 spoligotyping-defined clusters could be further subdivided by IS6110-RFLP. Reinterpretation of previous IS6110-RFLP results in the light of spoligotyping-VNTR typing results allowed us to detect an additional cluster that was previously missed. Although less discriminative than IS6110-RFLP, our results suggest that the use of the combination of spoligotyping and VNTR typing is a good screening strategy for detecting epidemiological links for the study of tuberculosis epidemiology at the molecular level.