Clonal expansion of a globally disseminated lineage of Mycobacterium tuberculosis with low IS6110 copy numbers

Distribution and Clonality of drug-resistant tuberculosis in South Africa

Background

Studies have shown that drug-resistant tuberculosis (DR-TB) in South Africa (SA) is clonal and is caused mostly by transmission. Identifying transmission chains is important in controlling DR-TB. This study reports on the sentinel molecular surveillance data of Rifampicin-Resistant (RR) TB in SA, aiming to describe the RR-TB strain population and the estimated transmission of RR-TB cases.

Method

RR-TB isolates collected between 2014 and 2018 from eight provinces were genotyped using combination of spoligotyping and 24-loci mycobacterial interspersed repetitive-units-variable-number tandem repeats (MIRU-VNTR) typing.

Results

Of the 3007 isolates genotyped, 301 clusters were identified. Cluster size ranged between 2 and 270 cases. Most of the clusters (247/301; 82.0%) were small in size (< 5 cases), 12.0% (37/301) were medium sized (5–10 cases), 3.3% (10/301) were large (11–25 cases) and 2.3% (7/301) were very large with 26–270 cases. The Beijing genotype was responsible for majority of RR-TB cases in Western and Eastern Cape, while the East-African-Indian-Somalian (EAI1_SOM) genotype accounted for a third of RR-TB cases in Mpumalanga. The overall proportion of RR-TB cases estimated to be due to transmission was 42%, with the highest transmission-rate in Western Cape (64%) and the lowest in Northern Cape (9%).

Conclusion

Large clusters contribute to the burden of RR-TB in specific geographic areas such as Western Cape, Eastern Cape and Mpumalanga, highlighting the need for community-wide interventions. Most of the clusters identified in the study were small, suggesting close contact transmission events, emphasizing the importance of contact investigations and infection control as the primary interventions in SA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02232-z.

The Epidemiological Significance and Temporal Stability of Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats-Based Method Applied to Mycobacterium tuberculosis in China

This study aimed to validate the epidemiological significance and temporal stability of Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) typing in a genetically and geographically diverse set of clinical isolates from patients diagnosed with pulmonary tuberculosis in China. Between 2010 and 2013, a total of 982 Mycobacterium tuberculosis isolates were collected from four population-based investigations in China. Apart from the currently applied 24-locus MIRU-VNTR, six additional hypervariable loci were analyzed in order to validate the MIRU-VNTR combinations in terms of their epidemiological links, clustering time span, and paired geographic distance. In vitro temporal stability was analyzed for both individual MIRU-VNTR loci, and for several combinations of loci. In the present study, four MIRU-VNTR combinations, including the hypervariable loci 3820, 3232, 2163a, and 4120, were evaluated. All of these combinations obtained a Hunter-Gaston discriminatory index (HGDI) value over 0.9900 with a reduced clustering proportion (from 32.0% to 25.6%). By comparing epidemiological links, clustering time span, and paired geographic distance, we found that the performances of the four MIRU-VNTR combinations were comparable to the insertion sequence 6110 restriction fragment length polymorphism (IS6110-RFLP), and significantly better than that of 24-locus MIRU-VNTR genotyping alone. The proportion of temporally stable loci ranged from 90.5% to 92.5% within the combined MIRU-VNTR genotyping, which is higher than IS6110-RFLP (85.4%). By adding four hypervariable loci to the standard 24-locus MIRU-VNTR genotyping, we obtained a high discriminatory power, stability and epidemiological significance. This algorithm could therefore be used to improve tuberculosis transmission surveillance and outbreak investigation in China.

The Evolution of Strain Typing in the Mycobacterium tuberculosis Complex

Tuberculosis (TB) is a contagious disease with a complex epidemiology. Therefore, molecular typing (genotyping) of Mycobacterium tuberculosis complex (MTBC) strains is of primary importance to effectively guide outbreak investigations, define transmission dynamics and assist global epidemiological surveillance of the disease. Large-scale genotyping is also needed to get better insights into the biological diversity and the evolution of the pathogen. Thanks to its shorter turnaround and simple numerical nomenclature system, mycobacterial interspersed repetitive unit-variable-number tandem repeat (MIRU-VNTR) typing, based on 24 standardized plus 4 hypervariable loci, optionally combined with spoligotyping, has replaced IS6110 DNA fingerprinting over the last decade as a gold standard among classical strain typing methods for many applications. With the continuous progress and decreasing costs of next-generation sequencing (NGS) technologies, typing based on whole genome sequencing (WGS) is now increasingly performed for near complete exploitation of the available genetic information. However, some important challenges remain such as the lack of standardization of WGS analysis pipelines, the need of databases for sharing WGS data at a global level, and a better understanding of the relevant genomic distances for defining clusters of recent TB transmission in different epidemiological contexts. This chapter provides an overview of the evolution of genotyping methods over the last three decades, which culminated with the development of WGS-based methods. It addresses the relative advantages and limitations of these techniques, indicates current challenges and potential directions for facilitating standardization of WGS-based typing, and provides suggestions on what method to use depending on the specific research question.

Genetic Lineages of Mycobacterium tuberculosis Isolates in Isfahan, Iran

In this study, we aimed to identify the genetic lineages of Mycobacterium tuberculosis isolates in Isfahan via the mycobacterial interspersed repetitive-unit-variable number tandem repeat typing method based on 15 loci. Forty-nine M. tuberculosis isolates were collected between 2013 and 2015 from Tuberculosis patients in Mollahadi Sabzevari Tuberculosis Center in Isfahan. All isolates were typed by 15-locus MIRU-VNTR typing. The highest percentage of isolates, 44.89 % (22/49), belonged to the Euro-American lineage, while the frequencies of the East-African-Indian, East-Asian, and Indo-Oceanic lineages were 28.57 % (14/49), 24.4 % (12/49), and 2.04 % (1/49), respectively. Among the 22 isolates of the Euro-American lineage, those belonging to the NEW-1 sub-lineage were most prevalent (24.4 %). Approximately, the same proportion of isolates belonging to the Delhi/CAS, Beijing, and NEW-1 sub-lineages were identified in Iranian and Afghan immigrant patients. The Delhi/CAS and Beijing sub-lineage isolates were prevalent among patients who had been previously treated for TB. Results showed that all of the 49 MIRU-VNTR patterns were unique and the clustering rate of the 15-locus MIRU-VNTR was 0.0 (minimum recent transmission). The results of this study show that the lineages of M. tuberculosis isolates in Isfahan are similar to those reported in the Eastern Mediterranean region (indicative of the epidemiological relationship between the countries in the region). The low clustering rate in our results reveals that transmission of tuberculosis in Isfahan is, in most cases, a reactivation of previous tuberculosis infection and the role of recently transmitted disease is minor.

Molecular epidemiology of Mycobacterium tuberculosis among South African gold miners

RATIONALE

HIV-associated tuberculosis remains a major health problem among the gold-mining workforce in South Africa. We postulate that high levels of recent transmission, indicated by strain clustering, are fueling the tuberculosis epidemic among gold miners.

OBJECTIVES

To combine molecular and epidemiologic data to describe Mycobacterium tuberculosis genetic diversity, estimate levels of transmission, and examine risk factors for clustering.

METHODS

We conducted a cross-sectional study of culture-positive M. tuberculosis isolates in 15 gold mine shafts across three provinces in South Africa. All isolates were subject IS6110-based restriction fragment length polymorphisms, and we performed spoligotyping analysis and combined it with basic demographic and clinical information.

MEASUREMENTS AND MAIN RESULTS

Of the 1,602 M. tuberculosis patient isolates, 1,240 (78%) had genotyping data available for analysis. A highly diverse bacillary population was identified, comprising a total of 730 discrete genotypes. Four genotypic families (Latin American Mediterranean spoligotype family; W-Beijing; AH or X; and T1-T4) accounted for over 50% of all strains. Overall, 45% (560/1,240) of strains were genotypically clustered. The minimum estimate for recent transmission (n - 1 method) was 32% (range, 27-34%). There were no individual-level risk factors for clustering, apart from borderline evidence for being non-South African and having self-reported HIV infection.

CONCLUSIONS

The high M. tuberculosis genetic diversity and lack of risk factors for clustering are indicative of a universal risk for disease among gold miners and likely mixing with nonmining populations. Our results underscore the urgent need to intensify interventions to interrupt transmission across the entire gold-mining workforce in South Africa.

Molecular Typing of Mycobacterium Tuberculosis Complex by 24-Locus Based MIRU-VNTR Typing in Conjunction with Spoligotyping to Assess Genetic Diversity of Strains Circulating in Morocco

Background

Standard 24-locus Mycobacterial Interspersed Repetitive Unit Variable Number Tandem Repeat (MIRU-VNTR) typing allows to get an improved resolution power for tracing TB transmission and predicting different strain (sub) lineages in a community.

Methodology

During 2010–2012, a total of 168 Mycobacterium tuberculosis Complex (MTBC) isolates were collected by cluster sampling from 10 different Moroccan cities, and centralized by the National Reference Laboratory of Tuberculosis over the study period. All isolates were genotyped using spoligotyping, and a subset of 75 was genotyped using 24-locus based MIRU-VNTR typing, followed by first line drug susceptibility testing. Corresponding strain lineages were predicted using MIRU-VNTRplus database.

Principal Findings

Spoligotyping resulted in 137 isolates in 18 clusters (2–50 isolates per cluster: clustering rate of 81.54%) corresponding to a SIT number in the SITVIT database, while 31(18.45%) patterns were unique of which 10 were labelled as “unknown” according to the same database. The most prevalent spoligotype family was LAM; (n = 81 or 48.24% of isolates, dominated by SIT42, n = 49), followed by Haarlem (23.80%), T superfamily (15.47%), >Beijing (2.97%), > U clade (2.38%) and S clade (1.19%). Subsequent 24-Locus MIRU-VNTR typing identified 64 unique types and 11 isolates in 5 clusters (2 to 3isolates per cluster), substantially reducing clusters defined by spoligotyping only. The single cluster of three isolates corresponded to two previously treated MDR-TB cases and one new MDR-TB case known to be contact a same index case and belonging to a same family, albeit residing in 3 different administrative regions. MIRU-VNTR loci 4052, 802, 2996, 2163b, 3690, 1955, 424, 2531, 2401 and 960 were highly discriminative in our setting (HGDI >0.6).

Conclusions

24-locus MIRU-VNTR typing can substantially improve the resolution of large clusters initially defined by spoligotyping alone and predominating in Morocco, and could therefore be used to better study tuberculosis transmission in a population-based, multi-year sample context.

Diversity and evolution of Mycobacterium tuberculosis: moving to whole-genome-based approaches

Genotyping of clinical Mycobacterium tuberculosis complex (MTBC) strains has become a standard tool for epidemiological tracing and for the investigation of the local and global strain population structure. Of special importance is the analysis of the expansion of multidrug (MDR) and extensively drug-resistant (XDR) strains. Classical genotyping and, more recently, whole-genome sequencing have revealed that the strains of the MTBC are more diverse than previously anticipated. Globally, several phylogenetic lineages can be distinguished whose geographical distribution is markedly variable. Strains of particular (sub)lineages, such as Beijing, seem to be more virulent and associated with enhanced resistance levels and fitness, likely fueling their spread in certain world regions. The upcoming generalization of whole-genome sequencing approaches will expectedly provide more comprehensive insights into the molecular and epidemiological mechanisms involved and lead to better diagnostic and therapeutic tools.

Genetic diversity of Mycobacterium tuberculosis isolates from Beijing, China assessed by Spoligotyping, LSPs and VNTR profiles

Background

Tuberculosis is one of the most infectious diseases in the world. Molecular typing methods such as spoligotyping, and VNTR (variable number tandem repeats), IS6110 in the NTF region and LSP (large sequence polymorphisms) analysis are generally useful tools for the resolution of various issues related to the classical epidemiology of Mycobacterium tuberculosis (M. tuberculosis).

Methods

To determine the transmission characteristics of M. tuberculosis strains isolated in Beijing, China, and their genetic relationships, especially those among Beijing family strains, 260 M. tuberculosis strains isolated from patients presenting pulmonary tuberculosis were analyzed by spoligotyping, and by examining 22 VNTR loci and the presence/absence of IS6110 in the NTF region, RD105 and RD181.

Results

81% (211 strains) of the isolates studied were Beijing family strains, 174 (82.5%) of which were identified as modern Beijing strains based on the presence of IS6110 upstream of the NTF region. RD181 was intact in 9 of the other 37 (17.5%) ancestral Beijing strains. The percentage of Beijing family strains in this study was consistent with previous reports. There are many differences, however, in allele diversity among VNTR loci between reports on strains from different areas.

Conclusions

The Beijing family is the most prevalent genotype in Beijing city and the predominance of Beijing family strains has not altered in almost twenty years. Differences in the alleles and discrimination ability of VNTR loci between different regions is likely due to population differences in the regions where these M. tuberculosis strains were isolated or to differences in sampling times.

gyrA mutations and phenotypic susceptibility levels to ofloxacin and moxifloxacin in clinical isolates of Mycobacterium tuberculosis

OBJECTIVES

To compare mutations in the quinolone resistance-determining region of the gyrA gene and flanking sequences with the MICs of ofloxacin and moxifloxacin for Mycobacterium tuberculosis.

METHODS

The presence of mutations in 177 drug-resistant M. tuberculosis isolates was determined by DNA sequencing and the MICs quantified by MGIT 960.

RESULTS

Single nucleotide polymorphisms were detected at codons 94 (n = 30), 90 (n = 12), 91 (n = 3), 89 (n = 1), 88 (n = 1) and 80 (n = 1). Four isolates with double mutations D94G plus A90V (n = 2) and D94G plus D94N (n = 2) reflect mixed populations. Agreement between genotypic and phenotypic susceptibility was high (≥97%) for both drugs. Mutant isolates had an MIC(50) of 8.0 mg/L and an MIC(90) of >10 mg/L for ofloxacin compared with an MIC(50) and MIC(90) of 2.0 mg/L for moxifloxacin. Codons 94 and 88 were linked to higher levels of fluoroquinolone resistance compared with codons 90, 91 and 89. The MIC distributions for the wild-type isolates ranged from ≤0.5 to 2.0 mg/L for ofloxacin and from ≤0.125 to 0.25 mg/L for moxifloxacin. However, 96% of the isolates with genetic alterations had MICs ≤2.0 mg/L for moxifloxacin, which is within its achievable serum levels.

CONCLUSIONS

This study provides quantitative evidence that the addition of moxifloxacin to extensively drug-resistant tuberculosis (XDR-TB) regimens based on a clinical breakpoint of 2.0 mg/L has merit. The use of moxifloxacin in the treatment of multidrug-resistant tuberculosis may prevent the acquisition of additional mutations and development of XDR-TB.

Genetic diversity of the Mycobacterium tuberculosis Beijing family based on IS6110, SNP, LSP and VNTR profiles from Thailand

The Beijing family of Mycobacterium tuberculosis (MTB) has been shown to be highly virulent and is associated with tuberculosis outbreaks and drug resistance. To explore the evolutionary pathway of this genotype, 165 epidemiologically unlinked Beijing MTB isolates from Thailand were characterized based on IS6110, variable number of tandem repeats (VNTRs), region of differences (RDs) and single nucleotide polymorphisms (SNPs). Analysis based on VNTRs showed that 7 VNTRs were phylogenetically informative loci. Three of the seven VNTRs, 4156, 1982 and 1955, were promising markers for defining ancestral and modern Beijing sublineages. The evolutionary history of Beijing sublineages based on VNTR profiles showed dynamic changes in the repeat number, which also correlated with specific SNP STs. We constructed a refined, comprehensive phylogenetic tree based on multiple genetic markers from the Beijing strains from Thailand and combined the data with previous studies that were performed in South Africa and Japan. The phylogenetic tree showed how the modern Beijing and W-Beijing strains evolved from ancestral Beijing strains. The points at which the RD150 deletion occurred repeatedly were also demonstrated. The evolutionary scheme showed that most events were concordant with the SNP STs. Nevertheless, some SNPs used for sub-typing the Beijing family of MTB may not be irreversible, unique events.

MIRU-VNTR typing of drug-resistant tuberculosis isolates in Greece

The increasing immigration rate in Greece from countries with a high prevalence of Mycobacterium tuberculosis (MTB) and multidrug-resistant tuberculosis (MDR-TB) may have an impact οn the number of MDR-TB cases in Greece. The aim of this study was to genotypically characterize the MTB isolates from patients with pulmonary drug-resistant tuberculosis (DR-TB) in Greece, and to determine whether there is any association between the prevalent genotypes and drug resistance. Fifty-three drug-resistant MTB strains isolated from culture specimens of clinical material from native Greeks and immigrant patients with pulmonary tuberculosis were genotyped using the mycobacterial interspersed repetitive units–variable number of tandem repeats (MIRU-VNTR) method. The phylogenetically distinct groups of isolates identified were: the Beijing (34%), the LAM (11%), the Haarlem (24.5%), the Uganda I (9.4%), the Ural (3.8%), the Delhi/CAS (9.4%) and the Cameroon (3.8%) families. Greek patients were more likely to have monoresistant and polyresistant TB with the most prevalent isolates belonging to the Haarlem family. Among foreign-born patients with MDR-TB, the most prevalent genotypes belonged to the Beijing family. MIRU-VNTR rapidly obtained clinically useful genotyping data, by characterizing clonal MTB heterogeneity in the isolated strains. Our results underline the need for more effective antituberculosis control programs in order to control the expansion of DR-TB in Greece.

The forest behind the tree: phylogenetic exploration of a dominant Mycobacterium tuberculosis strain lineage from a high tuberculosis burden country

BACKGROUND

Genotyping of Mycobacterium tuberculosis isolates is a powerful tool for epidemiological control of tuberculosis (TB) and phylogenetic exploration of the pathogen. Standardized PCR-based typing, based on 15 to 24 mycobacterial interspersed repetitive unit-variable number of tandem repeat (MIRU-VNTR) loci combined with spoligotyping, has been shown to have adequate resolution power for tracing TB transmission and to be useful for predicting diverse strain lineages in European settings. Its informative value needs to be tested in high TB-burden countries, where the use of genotyping is often complicated by dominance of geographically specific, genetically homogeneous strain lineages.

METHODOLOGY/PRINCIPAL FINDINGS

We tested this genotyping system for molecular epidemiological analysis of 369 M. tuberculosis isolates from 3 regions of Brazil, a high TB-burden country. Deligotyping, targeting 43 large sequence polymorphisms (LSPs), and the MIRU-VNTRplus identification database were used to assess phylogenetic predictions. High congruence between the different typing results consistently revealed the countrywide supremacy of the Latin-American-Mediterranean (LAM) lineage, comprised of three main branches. In addition to an already known RDRio branch, at least one other branch characterized by a phylogenetically informative LAM3 spoligo-signature seems to be globally distributed beyond Brazil. Nevertheless, by distinguishing 321 genotypes in this strain population, combined MIRU-VNTR typing and spoligotyping demonstrated the presence of multiple distinct clones. The use of 15 to 24 loci discriminated 21 to 25% more strains within the LAM lineage, compared to a restricted lineage-specific locus set suggested to be used after SNP analysis. Noteworthy, 23 of the 28 molecular clusters identified were exclusively composed of patient isolates from a same region, consistent with expected patterns of mostly local TB transmission.

CONCLUSIONS/SIGNIFICANCE

Standard MIRU-VNTR typing combined with spoligotyping can reveal epidemiologically meaningful clonal diversity behind a dominant M. tuberculosis strain lineage in a high TB-burden country and is useful to explore international phylogenetical ramifications.

The non-clonality of drug resistance in Beijing-genotype isolates of Mycobacterium tuberculosis from the Western Cape of South Africa

Background

The Beijing genotype of M. tuberculosis is a virulent strain that is disseminating worldwide and has a strong association with drug resistance. In the Western Cape of South Africa, epidemiological studies have identified the R220 cluster of the Beijing genotype as a major contributor to a recent outbreak of drug-resistant tuberculosis. Although the outbreak is considered to be due to clonal transmission, the relationship among drug resistant isolates has not yet been established.

Results

To better understand the evolution of drug resistance among these strains, 14 drug-resistant clinical isolates of the Beijing genotype were sequenced by whole-genome sequencing, including eight from R220 and six from a more ancestral Beijing cluster, R86, for comparison. While each cluster shares a distinct resistance mutation for isoniazid, mapping of other drug-resistance mutations onto a phylogenetic tree constructed from single nucleotide polymorphisms shows that resistance mutations to many drugs have arisen multiple times independently within each cluster of isolates. Thus, drug resistance among these isolates appears to be acquired, not clonally derived. This observation suggests that, although the Beijing genotype as a whole might have selective advantages enabling its rapid dissemination, the XDR isolates are relatively less fit and do not propagate well. Although it has been hypothesized that the increased frequency of drug resistance in some Beijing lineages might be caused by a mutator phenotype, no significant shift in synonymous substitution patterns is observed in the genomes.

Conclusion

While MDR-TB is spreading by transmission in the Western Cape, our data suggests that further drug resistance (i.e. XDR-TB) at this stage is acquired.

Insertion sequences as highly resolutive genomic markers for sequence type 1 Legionella pneumophila Paris

The causative agent of legionellosis, Legionella pneumophila, colonizes all natural and human-made water networks, thus constituting the source of contaminated aerosols responsible for airborne human infections. Efficient control of infections, especially during epidemics, necessitates the fastest and most resolutive identification possible of the bacterial source for subsequent disinfection of reservoirs. We thus compared recognized typing approaches for Legionella with a method based on characterization of insertion sequence (IS) content. A total of 86 clinical or environmental isolates of L. pneumophila, including 84 Paris isolates, sampled from 25 clinical investigations in France between 2001 and 2007, were obtained from the Legionella National Reference Center. All strains were typed by monoclonal antibody subgrouping, sequence-based typing, pulsed-field gel electrophoresis, and restriction fragment length polymorphism based on the presence or absence of IS elements. We identified six different types of IS elements in L. pneumophila Paris and used them as genomic markers in hybridization experiments. One IS type, ISLpn11, revealed a high discriminatory power. Simpson's index of discrimination, calculated from the distribution of IS elements, was higher than that obtained with the other typing methods used for L. pneumophila Paris. Moreover, specific ISLpn11 copies were found only in strains isolated from particular cities. In more than half of the cases, each clinical isolate had an ISLpn11 profile that was recovered in at least one environmental isolate from the same geographical location, suggesting that our method could identify the infection source. Phylogenetic analysis suggests a clonal expansion for the L. pneumophila Paris strain.

Novel multiplex PCR using dual-priming oligonucleotides for detection and discrimination of the Mycobacterium tuberculosis complex and M. bovis BCG

We developed a novel multiplex PCR assay using dual-priming oligonucleotide primers targeting the RD1 gene for simultaneous identification of the Mycobacterium tuberculosis complex and M. bovis bacillus Calmette-Guérin (BCG). This assay would be useful both for detection of the M. tuberculosis complex and for differentiation of M. bovis BCG from pathogenic M. tuberculosis complex species.

Discriminatory ability of hypervariable variable number tandem repeat loci in population-based analysis of Mycobacterium tuberculosis strains, London, UK

To address conflicting results about the stability of variable number tandem repeat (VNTR) loci and their value in prospective molecular epidemiology of Mycobacterium tuberculosis, we conducted a large prospective population-based analysis of all M. tuberculosis strains in a metropolitan setting. Optimal and reproducible conditions for reliable PCR and fragment analysis, comprising enzymes, denaturing conditions, and capillary temperature, were identified for a panel of hypervariable loci, including 3232, 2163a, 1982, and 4052. A total of 2,261 individual M. tuberculosis isolates and 265 sets of serial isolates were analyzed by using a standardized 15-loci VNTR panel, then an optimized hypervariable loci panel. The discriminative ability of loci varied substantially; locus VNTR 3232 varied the most, with 19 allelic variants and Hunter-Gaston index value of 0.909. Hypervariable loci should be included in standardized panels because they can provide consistent comparable results at multiple settings, provided the proposed conditions are adhered to.

Mycobacterium tuberculosis spoligotypes in Monterrey, Mexico

Although tuberculosis is still a public health problem in Mexico, there is little information about the genetic characteristics of the isolates. In the present study, we analyzed by spoligotyping 180 Mycobacterium tuberculosis clinical isolates from the urban area of Monterrey, Mexico, including drug-susceptible and drug-resistant isolates. The spoligotype patterns were compared with those in the international SITVIT2 spoligotyping database. Four isolates presented spoligotype patterns not found in the database (orphan types); the rest were distributed among 44 spoligo international types (SITs). SIT53 (clade T1) and SIT119 (clade X1) were predominant and included 43 (23.8%) and 28 (15.5%) of the isolates, respectively. In order to determine if there was a dominant spoligotype in the group of multidrug-resistant isolates, 37 of them were analyzed by IS6110-based restriction fragment length polymorphism assays, and scarce clustering of strains with more than five bands was observed. Fourteen isolates of this multidrug-resistant group presented four bands or less and were distributed in four SITs: SIT53 (n = 8), SIT92 (n = 3), SIT70 (n = 2), and SIT3038 (n = 1). When the molecular detection of mutations in the katG and rpoB genes were analyzed in these isolates with low copy numbers of IS6110, only two isolates shared the same IS6110, spoligotyping, and mutations patterns. When the distribution of the spoligotypes was analyzed by age cohort, SIT119 was predominantly found in patients 0 to 20 years old, especially in males, accounting for up to 40% of the isolates. In contrast, SIT53 was more prevalent in older females. This analysis demonstrates the variability of M. tuberculosis isolates in Monterrey and the partial dominance of SIT53 and SIT119 in that area of Mexico.

Multiple, independent, identical IS6110 insertions in Mycobacterium tuberculosis PPE genes

IS6110 is a transposable element found in Mycobacterium tuberculosis complex members. Regions of preferential IS6110 integration occur within the M. tuberculosis genome but the element has not previously been shown to exhibit any sequence-specific integration preferences. Here we provide evidence for multiple independent IS6110 insertions into identical, or near-identical, positions within the highly homologous 5' region of 3 PPE genes.

Changing Mycobacterium tuberculosis population highlights clade-specific pathogenic characteristics

Mycobacterium tuberculosis strains can be classified into a number of major clades according to defined evolutionary markers. It is hypothesised that strains comprising these clades have evolved different properties which may influence a local strain population structure. To investigate this, we analysed the incidence of tuberculosis caused by the predominant clades (Beijing, Haarlem, LAM, Quebec and the Low-Copy Clade) found in a community within the Cape Town metropole in South Africa over a 12-year period. We found that while the incidence of cases infected with strains of the Haarlem, LAM, Quebec and the Low-Copy Clades remained relatively stable, that of cases of the Beijing clade increased exponentially over time, with a doubling time of 4.86 years (P=0.018). This growth was exclusively attributable to drug-susceptible strains. Although drug-resistant Beijing cases remained constant in number, non-Beijing drug-resistant cases declined over time (P=0.007). Drug-susceptible Beijing-infected cases had a greater proportion of smear-positive sputa than their non-Beijing counterparts (P=0.013) and were less likely to be successfully treated (retreatment cases) (P=0.026). Recent evidence suggests that these differences likely reflect enhanced pathogenicity rather than transmissibility. The rapid emergence of Beijing strains demonstrates adaptation to conditions within the study community and poses a grave challenge to future TB control.

Discordance between mycobacterial interspersed repetitive-unit-variable-number tandem-repeat typing and IS6110 restriction fragment length polymorphism genotyping for analysis of Mycobacterium tuberculosis Beijing strains in a setting of high incidence of tuberculosis

IS6110 restriction fragment length polymorphism (RFLP) genotyping is the most widely used genotyping method to study the epidemiology of Mycobacterium tuberculosis. However, due to the complexity of the IS6110 RFLP genotyping technique, and the interpretation of RFLP data, mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) genotyping has been proposed as the new genotyping standard. This study aimed to determine the discriminatory power of different MIRU-VNTR locus combinations relative to IS6110 RFLP genotyping, using a collection of Beijing genotype M. tuberculosis strains with a well-established phylogenetic history. Clustering, diversity index, clustering concordance, concordance among unique genotypes, and divergent and convergent evolution were calculated for seven combinations of 27 different MIRU-VNTR loci and compared to IS6110 RFLP results. Our results confirmed previous findings that MIRU-VNTR genotyping can be used to estimate the extent of recent or ongoing transmission. However, molecular epidemiological linking of cases varied significantly depending on the genotyping method used. We conclude that IS6110 RFLP and MIRU-VNTR loci evolve independently and at different rates, which leads to discordance between transmission chains predicted by the respective genotyping methods. Concordance between the two genotyping methods could be improved by the inclusion of genetic distance (GD) into the clustering formulae for some of the MIRU-VNTR loci combinations. In summary, our findings differ from previous reports, which may be explained by the fact that in settings of low tuberculosis incidence, the genetic distance between epidemiologically unrelated isolates was sufficient to define a strain using either marker, whereas in settings of high incidence, continuous evolution and persistence of strains revealed the weaknesses inherent to these markers.

Delayed behavioral and endocrine effects of sarin and stress exposure in mice

The organophosphorus agent sarin is a potent inhibitor of acetylcholinesterase. Experiments tested the influence of exposure to low doses of sarin along with psychological stress on delayed behavioral and endocrine changes in mice. Motor activity, acoustic startle response (ASR), pre-pulse inhibition (PPI) of ASR, activity of cholinesterase in blood and catecholamine levels in adrenals were evaluated after low dose sarin exposure (3 x 0.4 LD50 subcutaneously) combined with chronic intermittent stress in C57BL/6J mice. While sarin alone produced depression of motor activity, no interaction of the stress with sarin exposure was observed. Cholinesterase activity was significantly reduced 24 h after exposure to sarin; however, the basal activity was re-established 3 weeks later. The combination of low dose sarin exposure and stress produced delayed behavioral change manifested as excessive grooming together with endocrine alterations in adrenals 7 weeks after exposure. The size of the adrenals in the combined exposure group was increased and the concentration of catecholamines was significantly decreased. In conclusion, these findings indicate that sarin in low doses is more dangerous when combined with shaker stress inducing delayed behavioral and endocrine changes.

Assessment of population structure and major circulating phylogeographical clades of Mycobacterium tuberculosis complex in Bangladesh suggests a high prevalence of a specific subclade of ancient M. tuberculosis genotypes

Spoligotyping was performed to study the population structure of Mycobacterium tuberculosis complex strains (n = 224) from Bangladesh. Strains were split into principal genetic group 1 (PGG 1 [75.0%]) and PGG 2 and 3 (25%). Forty-nine strains with a new spoligotype signature and considered as south or southeast Asian-linked emerging clones were designated as "Matlab type."

Evaluation of variable numbers of tandem repeat as molecular epidemiological markers of Mycobacterium tuberculosis in Japan

Using 243 Mycobacterium tuberculosis isolates obtained in 2001 in Osaka City, Japan, the discriminatory power of variable numbers of tandem repeats (VNTRs) of 12 standard mycobacterial interspersed repetitive units (MIRUs) was assessed. The biggest cluster defined by MIRU-VNTRs consisted of 57 (23.5 %) isolates and they belonged to the Beijing family based on spoligotyping. When additional VNTR loci were included in the MIRU-VNTR analysis, the 57 originally clustered strains were further differentiated by the addition of Queen's University Belfast (QUB)-VNTRs, but not exact tandem repeat-VNTR. The allelic diversity of additional VNTR loci such as VNTR 3232 (QUB-3232), VNTR 2163a (QUB-11a), VNTR 2163b (QUB-11b) and VNTR 1982 (QUB-18) was high in the 57 strains. When the 243 M. tuberculosis isolates were analysed using 16-locus VNTR (the 12 standard MIRUs and the 4 QUB loci) and IS6110 RFLP, the respective Hunter-Gaston discriminatory indexes were 0.9966 and 0.9971. The discrimination power of 16-locus VNTR was equal to that of IS6110 RFLP analysis. If appropriate loci are added to the standard MIRU analysis, VNTR genotyping could be a valuable tool for strain typing and epidemiological research of M. tuberculosis in Japan.

Characterization of Mycobacterium tuberculosis Central Asian Strain 1 using mycobacterial interspersed repetitive unit genotyping

BACKGROUND

The Central Asian Strain 1 (CAS1) genogroup of Mycobacterium tuberculosis (MTB) is the most prevalent in Pakistan, India and Bangladesh. Mycobacterial interspersed repetitive units variable number tandem repeat (MIRU-VNTR) typing is a reliable and reproducible method for differentiation of MTB isolates. However, information of its utility in determining the diversity of CAS1 strain is limited. We performed standard 12 loci based MIRU-VNTR typing on previously spoligotyped CAS1 strains and 'unique' strains in order to evaluate its discriminatory power for these isolates.

METHODS

Twelve loci based MIRU- VNTR typing was used to type 178 CAS1 and 189 'unique' MTB strains. The discriminatory index for each of the loci was calculated using the Hunter Gaston Discriminatory Index (HGDI). A subset of these strains (n = 78) were typed using IS6110 restriction fragment length polymorphism (RFLP). MIRU-VNTR profiles were studied together with their drug susceptibility patterns.

RESULTS

A total of 349 MIRU patterns were obtained for the 367 strains tested. The CAS1 strains were subdivided into 160 distinct patterns; 15 clusters of 2 strains each, 1 cluster of four strains and 144 unique patterns. Using HGDI, seven MIRU loci, (numbers 26, 31, 27, 16, 10, 39, and 40) were found to be "highly discriminatory" (DI: >or=0.6), four MIRU loci (numbers 20, 24, 23, and 4) were "moderately discriminatory" (DI: 0.3-0.59), and one locus (number 2) was "poorly discriminatory" (DI< 0.3). Loci 26 and 31 were the most discriminatory for the CAS1 isolates. Amongst 'unique' strains in addition to loci 26, 31, 27, 16, 10, 39, and 40, locus 23 was highly discriminatory, while no locus was poorly discriminating. DI values for loci 4, 10 and 26 were significantly lower (P-value < .01) in CAS1 strains than in 'unique' strains. The association between CAS1 strains and MDR was not found to be significant (p value = 0.21).

CONCLUSION

We propose that MIRU typing could be used to estimate the phylogenetic relatedness amongst prevalent CAS1 strains, for which MIRU loci 26, 31, 16, 10, 27, 39 and 40 were found to be the most discriminatory.

Predominance of ancestral lineages of Mycobacterium tuberculosis in India

Molecular epidemiologic findings suggest an ancient focus of TB.

Although India has the highest prevalence of tuberculosis (TB) worldwide, the genetic diversity of Mycobacterium tuberculosis in India is largely unknown. A collection of 91 isolates originating from 12 different regions spread across the country were analyzed by genotyping using 21 loci with variable-number tandem repeats (VNTRs), by spoligotyping, by principal genetic grouping (PGG), and by deletion analysis of M. tuberculosis–specific deletion region 1. The isolates showed highly diverse VNTR genotypes. Nevertheless, highly congruent groupings identified by using the 4 independent sets of markers permitted a clear definition of 3 prevalent PGG1 lineages, which corresponded to the "ancestral" East African–Indian, the Delhi, and the Beijing/W genogroups. A few isolates from PGG2 lineages and a single representative of the presumably most recent PGG3 were identified. These observations suggest a predominance of ancestral M. tuberculosis genotypes in the Indian subcontinent, which supports the hypothesis that India is an ancient endemic focus of TB.

Spoligotype signatures in the Mycobacterium tuberculosis complex

Evolution of the direct repeat region in Mycobacterium tuberculosis has created unique spoligotype signatures specifically associated with IS6110-defined strain families. Spoligotyping signatures may enable the analysis of the strain population structure in different settings and will enable the rapid identification of strain families that acquire drug resistance or escape protective immunity in drug and vaccine trials.

Molecular epidemiology of tuberculosis: current insights

Molecular epidemiologic studies of tuberculosis (TB) have focused largely on utilizing molecular techniques to address short- and long-term epidemiologic questions, such as in outbreak investigations and in assessing the global dissemination of strains, respectively. This is done primarily by examining the extent of genetic diversity of clinical strains of Mycobacterium tuberculosis. When molecular methods are used in conjunction with classical epidemiology, their utility for TB control has been realized. For instance, molecular epidemiologic studies have added much-needed accuracy and precision in describing transmission dynamics, and they have facilitated investigation of previously unresolved issues, such as estimates of recent-versus-reactive disease and the extent of exogenous reinfection. In addition, there is mounting evidence to suggest that specific strains of M. tuberculosis belonging to discrete phylogenetic clusters (lineages) may differ in virulence, pathogenesis, and epidemiologic characteristics, all of which may significantly impact TB control and vaccine development strategies. Here, we review the current methods, concepts, and applications of molecular approaches used to better understand the epidemiology of TB.

Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis

Molecular typing based on 12 loci containing variable numbers of tandem repeats of mycobacterial interspersed repetitive units (MIRU-VNTRs) has been adopted in combination with spoligotyping as the basis for large-scale, high-throughput genotyping of Mycobacterium tuberculosis. However, even the combination of these two methods is still less discriminatory than IS6110 fingerprinting. Here, we define an optimized set of MIRU-VNTR loci with a significantly higher discriminatory power. The resolution and the stability/robustness of 29 loci were analyzed, using a total of 824 tubercle bacillus isolates, including representatives of the main lineages identified worldwide so far. Five loci were excluded for lack of robustness and/or stability in serial isolates or isolates from epidemiologically linked patients. The use of the 24 remaining loci increased the number of types by 40%--and by 23% in combination with spoligotyping--among isolates from cosmopolitan origins, compared to those obtained with the original set of 12 loci. Consequently, the clustering rate was decreased by fourfold--by threefold in combination with spoligotyping--under the same conditions. A discriminatory subset of 15 loci with the highest evolutionary rates was then defined that concentrated 96% of the total resolution obtained with the full 24-locus set. Its predictive value for evaluating M. tuberculosis transmission was found to be equal to that of IS6110 restriction fragment length polymorphism typing, as shown in a companion population-based study. This 15-locus system is therefore proposed as the new standard for routine epidemiological discrimination of M. tuberculosis isolates and the 24-locus system as a high-resolution tool for phylogenetic studies.

Evolution of short sequence repeats in Mycobacterium tuberculosis

Whole genome comparison has revealed the presence of short sequence repeats (also called mycobacterial interspersed repeat units and variable number tandem repeat units) used for genotyping schemes. In this study, we have used deletion analysis, single nucleotide polymorphism data and spoligotype taken from published data from others to investigate the evolution of selected repeats that form the common denominators of the majority of established schemes. Analysis of the number of repeats per locus from over 400 isolates revealed that the general trend globally appears to be loss of repeats in modern strains compared with ancestral strains.

Mixed infection and clonal representativeness of a single sputum sample in tuberculosis patients from a penitentiary hospital in Georgia

Background

Studies on recurrent tuberculosis (TB), TB molecular epidemiology and drug susceptibility testing rely on the analysis of one Mycobacterium tuberculosis isolate from a single sputum sample collected at different disease episodes. This scheme rests on the postulate that a culture of one sputum sample is homogeneous and representative of the total bacillary population in a patient.

Methods

We systematically analysed several pre-treatment isolates from each of 199 smear-positive male adult inmates admitted to a prison TB hospital by standard IS6110 DNA fingerprinting, followed by PCR typing based on multiple loci containing variable number of tandem repeats (VNTRs) on a subset of isolates. Drug susceptibility testing (DST) was performed on all isolates for isoniazid, rifampicin, streptomycin and ethambutol.

Results

We found mixed infection in 26 (13.1%) cases. In contrast, analysis of a single pre-treatment isolate per patient would have led to missed mixed infections in all or 14 of these 26 cases by using only standard DNA fingerprinting or the PCR multilocus-based method, respectively. Differences in DST among isolates from the same patient were observed in 10 cases, of which 6 were from patients with mixed infection.

Conclusion

These results suggest that the actual heterogeneity of the bacillary population in patients, especially in high TB incidence settings, may be frequently underestimated using current analytical schemes. These findings have therefore important implications for correct interpretation and evaluation of molecular epidemiology data and in treatment evaluations.

Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology

Background

The Direct Repeat locus of the Mycobacterium tuberculosis complex (MTC) is a member of the CRISPR (Clustered regularly interspaced short palindromic repeats) sequences family. Spoligotyping is the widely used PCR-based reverse-hybridization blotting technique that assays the genetic diversity of this locus and is useful both for clinical laboratory, molecular epidemiology, evolutionary and population genetics. It is easy, robust, cheap, and produces highly diverse portable numerical results, as the result of the combination of (1) Unique Events Polymorphism (UEP) (2) Insertion-Sequence-mediated genetic recombination. Genetic convergence, although rare, was also previously demonstrated. Three previous international spoligotype databases had partly revealed the global and local geographical structures of MTC bacilli populations, however, there was a need for the release of a new, more representative and extended, international spoligotyping database.

Results

The fourth international spoligotyping database, SpolDB4, describes 1939 shared-types (STs) representative of a total of 39,295 strains from 122 countries, which are tentatively classified into 62 clades/lineages using a mixed expert-based and bioinformatical approach. The SpolDB4 update adds 26 new potentially phylogeographically-specific MTC genotype families. It provides a clearer picture of the current MTC genomes diversity as well as on the relationships between the genetic attributes investigated (spoligotypes) and the infra-species classification and evolutionary history of the species. Indeed, an independent Naïve-Bayes mixture-model analysis has validated main of the previous supervised SpolDB3 classification results, confirming the usefulness of both supervised and unsupervised models as an approach to understand MTC population structure. Updated results on the epidemiological status of spoligotypes, as well as genetic prevalence maps on six main lineages are also shown. Our results suggests the existence of fine geographical genetic clines within MTC populations, that could mirror the passed and present Homo sapiens sapiens demographical and mycobacterial co-evolutionary history whose structure could be further reconstructed and modelled, thereby providing a large-scale conceptual framework of the global TB Epidemiologic Network.

Conclusion

Our results broaden the knowledge of the global phylogeography of the MTC complex. SpolDB4 should be a very useful tool to better define the identity of a given MTC clinical isolate, and to better analyze the links between its current spreading and previous evolutionary history. The building and mining of extended MTC polymorphic genetic databases is in progress.

Molecular typing of Mycobacterium tuberculosis by mycobacterial interspersed repetitive unit-variable-number tandem repeat analysis, a more accurate method for identifying epidemiological links between patients with tuberculosis

IS6110 fingerprinting of Mycobacterium tuberculosis is the standard identification method in studies on transmission of tuberculosis. However, intensive epidemiological investigation may fail to confirm transmission links between patients clustered by IS6110-restriction fragment length polymorphism (RFLP) typing. We applied typing based on variable numbers of tandem repeats (VNTRs) of mycobacterial interspersed repetitive units (MIRUs) to isolates from 125 patients in 42 IS6110 clusters, for which thorough epidemiological data were available, to investigate the potential of this method in distinguishing epidemiologically linked from nonlinked patients. Of seven IS6110 clusters without epidemiological links, five were split by MIRU-VNTR typing, while nearly all IS6110 clusters with proven or likely links displayed conserved MIRU-VNTR types. These results provide molecular evidence that not all clusters determined on the basis of multibanded IS6110 RFLP patterns necessarily reflect transmission of tuberculosis. They support the use of MIRU-VNTR typing as a more reliable and faster method for transmission analysis.

Identification of moaA3 gene in patient isolates of Mycobacterium tuberculosis in Kerala, which is absent in M. tuberculosis H37Rv and H37Ra

Background

Tuberculosis is endemic to developing countries like India. Though the whole genome sequences of the type strain M. tuberculosis H37Rv and the clinical strain M. tuberculosis CDC1551 are available, the clinical isolates from India have not been studied extensively at the genome level. This study was carried out in order to have a better understanding of isolates from Kerala, a state in southern India.

Results

A PCR based strategy was followed making use of the deletion region primers to understand the genome level differences between the type strain H37Rv and the clinical isolates of M. tuberculosis from Kerala. PCR analysis of patient isolates using RD1 region primers revealed the amplification of a 386 bp region, in addition to the expected 652 bp amplicon. Southern hybridization of genomic DNA with the 386 bp amplicon confirmed the presence of this new region in a majority of the patient isolates from Kerala. Sequence comparison of this amplicon showed close homology with the moaA3 gene of M. bovis. In M. bovis this gene is present in the RvD5 region, an IS6110 mediated deletion that is absent in M. tuberculosis H37Rv.

Conclusion

This study demonstrates the presence of moaA3 gene, that is absent in M. tuberculosis H37Rv and H37Ra, in a large number of local isolates. Whether the moaA3 gene provides any specific advantage to the field isolates of the pathogen is unclear. Field strains from Kerala have fewer IS6110 sequences and therefore are likely to have fewer IS6110 dependent rearrangements. But as deletions and insertions account for much of the genomic diversity of M. tuberculosis, the mechanisms of formation of sequence polymorphisms in the local isolates should be further examined. These results suggest that studies should focus on strains from endemic areas to understand the complexities of this pathogen.

Multiple-locus variable number tandem repeats analysis for genetic fingerprinting of pathogenic bacteria

DNA fingerprinting has attracted considerable interest as means for identifying, tracing and preventing the dissemination of infectious agents. Various methods have been developed for typing of pathogenic bacteria, which differ in discriminative power, reproducibility and ease of interpretation. During recent years a typing method, which uses the information provided by whole genome sequencing of bacterial species, has gained increased attention. Short sequence repeat (SSR) motifs are known to undergo frequent variation in the number of repeated units through cellular mechanisms most commonly active during chromosome replication. A class of SSRs, named variable number of tandem repeats (VNTRs), has proven to be a suitable target for assessing genetic polymorphisms within bacterial species. This review attempts to give an overview of bacterial agents where VNTR-based typing, or multiple-locus variant-repeat analysis (MLVA) has been developed for typing purposes, together with addressing advantages and drawbacks associated with the use of tandem repeated DNA motifs as targets for bacterial typing and identification.

Genotypic and phenotypic heterogeneity among Mycobacterium tuberculosis isolates from pulmonary tuberculosis patients

Although the heterogeneity of Mycobacterium tuberculosis populations and the existence of mixed infections are now generally accepted, systematic studies on their relative importance are rare. In the present study, 10 individual colonies of each M. tuberculosis isolate (primary isolate) from 97 tuberculosis patients in a primarily human immunodeficiency virus-negative population were screened for heterogeneity and detectable mixed infections by spoligotyping, IS6110-based restriction fragment length polymorphism analysis, and mycobacterial interspersed repetitive unit-variable number of tandem repeat typing. The MICs of antituberculosis drugs for colonies with divergent fingerprints were determined. Infections with different bacterial subpopulations were detected in the samples from eight patients (8.2%), and the frequency of detectable mixed infections in the study population was estimated to be 2.1%. Genotypic variations were found to be independent of the drug susceptibilities, and the various molecular markers evolved independently in most cases. The predominant strains and the primary isolates always had concordant drug susceptibility and MIC testing results. These findings have implications on the interpretation of molecular epidemiology results for patient follow-up and in transmission studies.