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

Implementation of national whole-genome sequencing of Mycobacterium tuberculosis, National Public Health Laboratory, Singapore, 2019-2022

The National Tuberculosis Programme (NTBP) monitors the occurrence and spread of tuberculosis (TB) and multidrug-resistant TB (MDR-TB) in Singapore. Since 2020, whole-genome sequencing (WGS) of Mycobacterium tuberculosis isolates has been performed at the National Public Health Laboratory (NPHL) for genomic surveillance, replacing spoligotyping and mycobacterial interspersed repetitive unit-variable number tandem repeats analysis (MIRU-VNTR). Four thousand three hundred and seven samples were sequenced from 2014 to January 2023, initially as research projects and later developed into a comprehensive public health surveillance programme. Currently, all newly diagnosed culture-positive cases of TB in Singapore are prospectively sent for WGS, which is used to perform lineage classification, predict drug resistance profiles and infer genetic relationships between TB isolates. This paper describes NPHL's operational and technical experiences with implementing the WGS service in an urban TB-endemic setting, focusing on cluster detection and genomic drug susceptibility testing (DST). Cluster detection: WGS has been used to guide contact tracing by detecting clusters and discovering unknown transmission networks. Examples have been clusters in a daycare centre, housing apartment blocks and a horse-racing betting centre. Genomic DST: genomic DST prediction (gDST) identifies mutations in core genes known to be associated with TB drug resistance catalogued in the TBProfiler drug resistance mutation database. Mutations are reported with confidence scores according to a standardized approach referencing NPHL's internal gDST confidence database, which is adapted from the World Health Organization (WHO) TB drug mutation catalogue. Phenotypic-genomic concordance was observed for the first-line drugs ranging from 2959/2998 (98.7 %) (ethambutol) to 2983/2996 (99.6 %) (rifampicin). Aspects of internal database management, reporting standards and caveats in results interpretation are discussed.

Developing customized stepwise MIRU-VNTR typing for tuberculosis surveillance in Georgia

Introduction

Mycobacterial Interspersed Repetitive Units–Variable Tandem Repeats (MIRU-VNTR) typing has been widely used for molecular epidemiological studies of tuberculosis (TB). However, genotyping tools for Mycobacterium tuberculosis (Mtb) may be limiting in some settings due to high cost and workload. In this study developed a customized stepwise MIRU-VNTR typing that prioritizes high discriminatory loci and validated this method using penitentiary system cohort in the country of Georgia.

Methods

We used a previously generated MIRU-VNTR dataset from recurrent TB cases (32 cases) in Georgia and a new dataset of TB cases from the penitentiary system (102 cases) recruited from 2014 to 2015. A Hunter-Gaston Discriminatory Index (HGDI) was calculated utilizing a 24 standard loci panel, to select high discriminatory power loci, subsequently defined as the customized Georgia-specific set of loci for initial typing. The remaining loci were scored and hierarchically grouped for second and third step typing of the cohort. We then compared the processing time and costs of the customized stepwise method to the standard 24-loci method.

Results

For the customized Georgia-specific set that was used for initial typing, 10 loci were selected with a minimum value of 0.32 to the highest HGDI score locus. Customized 10 loci (step 1) typing of 102 Mtb patient isolates revealed 35.7% clustered cases. This proportion was reduced to 19.5% after hierarchical application of 2nd and 3rd step typing with the corresponding groups of loci. Our customized stepwise MIRU-VNTR genotyping approach reduced the quantity of samples to be typed and therefore overall processing time and costs by 42.6% each.

Conclusion

Our study shows that our customized stepwise MIRU-VNTR typing approach is a valid alternative of standard MIRI-VNTR typing panels for molecular epidemiological investigation in Georgia that saves time, workload and costs. Similar approaches could be developed for other settings.

A comparative analysis of molecular genotypes of Mycobacterium tuberculosis isolates from HIV-positive and HIV-negative patients

SETTING

Tuberculosis Research Laboratory, Division of Clinical Microbiology and Molecular Medicine, Department of Laboratory Medicine, All India Institute of Medical Sciences, and the National Institute of Tuberculosis and Respiratory Diseases (NITRD), both situated in New Delhi.

OBJECTIVES

We aimed to identify the distribution of various genotypes of M. tuberculosis among HIV-positive and HIV-negative patients suspected of having Tuberculosis, seen at the National Institute of Tuberculosis and Respiratory Diseases, New Delhi, which is a tertiary care dedicated TB hospital.

PATIENTS AND METHODS

Genotyping by Spoligotyping and 24 loci MIRU-VNTR was performed and analyzed using SITVITWEB and MIRU-VNTRplus. Drug susceptibility patterns were also analyzed.

RESULTS

A total of 503 subjects who were PTB/EPTB suspected were recruited and 287 were culture positive. Among them, 276 had growth of Mycobacterium tuberculosis (MTB) and in 11 patients non-tuberculous mycobacteria (NTM) were grown. The isolation rate of NTM was predominantly from HIV positive [10 of 130 (7.6%)] patients. Of the total isolates of MTB, 156 (56.5%) were from HIV negative patients and 120 (43.5%) were from HIV positive patients. All 276 M. tuberculosis isolates were genotyped and tested for drug susceptibility patterns. The CAS genotype was most predominant [153 (55.4%)], followed by Beijing lineage [44 (15.9%)], East African India [25 (9.1%)] and others [54 (19.6%)]. Beijing genotype was significantly more common in HIV positive patients (22.5%) than in HIV negative patients (10.9%). In MIRU-VNTR analysis, clustering was found to be more frequent in CAS strains irrespective of HIV status. In the HIV positive group, spoligotyping could differentiate various genotypes in 90% of isolates and MIRU-VNTR analysis in 84.2% of isolates. The clustering of various MTB strains was more associated with drug resistance.

CONCLUSION

The Beijing lineage was predominant in HIV-TB coinfected cases, even though the Central Asian Strain (CAS) was overall more predominant in the region.

The geno-spatio analysis of Mycobacterium tuberculosis complex in hot and cold spots of Guangxi, China

BACKGROUND

At present, there are few studies on polymorphism of Mycobacterium tuberculosis (Mtb) gene and how it affects the TB epidemic. This study aimed to document the differences of polymorphisms between tuberculosis hot and cold spot areas of Guangxi Zhuang Autonomous Region, China.

METHODS

The cold and hot spot areas, each with 3 counties, had been pre-identified by TB incidence for 5 years from the surveillance database. Whole genome sequencing analysis was performed on all sputum Mtb isolates from the detected cases during January and June 2018. Single nucleotide polymorphism (SNP) of each isolate compared to the H37Rv strain were called and used for lineage and sub-lineage identification. Pairwise SNP differences between every pair of isolates were computed. Analyses of Molecular Variance (AMOVA) across counties of the same hot or cold spot area and between the two areas were performed.

RESULTS

As a whole, 59.8% (57.7% sub-lineage 2.2 and 2.1% sub-lineage 2.1) and 39.8% (17.8% sub-lineage 4.4, 6.5% sub-lineage 4.2 and 15.5% sub-lineage 4.5) of the Mtb strains were Lineage 2 and Lineage 4 respectively. The percentages of sub-lineage 2.2 (Beijing family strains) are significantly higher in hot spots. Through the MDS dimension reduction, the genomic population structure in the three hot spot counties is significantly different from those three cold spot counties (T-test p = 0.05). The median of SNPs distances among Mtb isolates in cold spots was greater than that in hot spots (897 vs 746, Rank-sum test p < 0.001). Three genomic clusters, each with genomic distance ≤12 SNPs, were identified with 2, 3 and 4 consanguineous strains. Two clusters were from hot spots and one was from cold spots.

CONCLUSION

Narrower genotype diversity in the hot area may indicate higher transmissibility of the Mtb strains in the area compared to those in the cold spot area.

The Evolution of Genotyping Strategies To Detect, Analyze, and Control Transmission of Tuberculosis

The introduction of genotypic tools to analyze Mycobacterium tuberculosis isolates has transformed our knowledge of the transmission dynamics of this pathogen. We discuss the development of the laboratory methods that have been applied in recent years to study the epidemiology of M. tuberculosis. This review integrates two approaches: on the one hand, it considers how genotyping techniques have evolved over the years; and on the other, it looks at how the way we think these techniques should be applied has changed. We begin by examining the application of fingerprinting tools to suspected outbreaks only, before moving on to universal genotyping schemes, and finally we describe the latest real-time strategies used in molecular epidemiology. We also analyze refined approaches to obtaining epidemiological data from patients and to increasing the discriminatory power of genotyping by techniques based on genomic characterization. Finally, we review the development of integrative solutions to reconcile the speed of PCR-based methods with the high discriminatory power of whole-genome sequencing in easily implemented formats adapted to low-resource settings. Our analysis of future considerations highlights the need to bring together the three key elements of high-quality surveillance of transmission in tuberculosis, namely, speed, precision, and ease of implementation.

Genomic sequencing is required for identification of tuberculosis transmission in Hawaii

Background

Tuberculosis (TB) caused an estimated 1.4 million deaths and 10.4 million new cases globally in 2015. TB rates in the United States continue to steadily decline, yet rates in the State of Hawaii are perennially among the highest in the nation due to a continuous influx of immigrants from the Western Pacific and Asia. TB in Hawaii is composed of a unique distribution of genetic lineages, with the Beijing and Manila families of Mycobacterium tuberculosis (Mtb) comprising over two-thirds of TB cases. Standard fingerprinting methods (spoligotyping plus 24-loci Mycobacterial Interspersed Repetitive Units-Variable Number Tandem Repeats [MIRU-VNTR] fingerprinting) perform poorly when used to identify actual transmission clusters composed of isolates from these two families. Those typing methods typically group isolates from these families into large clusters of non-linked isolates with identical fingerprints. Next-generation whole-genome sequencing (WGS) provides a new tool for molecular epidemiology that can resolve clusters of isolates with identical spoligotyping and MIRU-VNTR fingerprints.

Methods

We performed WGS and SNP analysis and evaluated epidemiological data to investigate 19 apparent TB transmission clusters in Hawaii from 2003 to 2017 in order to assess WGS’ ability to resolve putative Mtb clusters from the Beijing and Manila families. This project additionally investigated MIRU-VNTR allele prevalence to determine why standard Mtb fingerprinting fails to usefully distinguish actual transmission clusters from these two Mtb families.

Results

WGS excluded transmission events in seven of these putative clusters, confirmed transmission in eight, and identified both transmission-linked and non-linked isolates in four. For epidemiologically identified clusters, while the sensitivity of MIRU-VNTR fingerprinting for identifying actual transmission clusters was found to be 100%, its specificity was only 28.6% relative to WGS. We identified that the Beijing and Manila families’ significantly lower Shannon evenness of MIRU-VNTR allele distributions than lineage 4 was the cause of standard fingerprinting’s poor performance when identifying transmission in Beijing and Manila family clusters.

Conclusions

This study demonstrated that WGS is necessary for epidemiological investigation of TB in Hawaii and the Pacific.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3502-1) contains supplementary material, which is available to authorized users.

Multi-clonal evolution of multi-drug-resistant/extensively drug-resistant Mycobacterium tuberculosis in a high-prevalence setting of Papua New Guinea for over three decades

An outbreak of multi-drug resistant (MDR) tuberculosis (TB) has been reported on Daru Island, Papua New Guinea. Mycobacterium tuberculosis strains driving this outbreak and the temporal accrual of drug resistance mutations have not been described. Whole genome sequencing of 100 of 165 clinical isolates referred from Daru General Hospital to the Supranational reference laboratory, Brisbane, during 2012–2015 revealed that 95 belonged to a single modern Beijing sub-lineage strain. Molecular dating suggested acquisition of streptomycin and isoniazid resistance in the 1960s, with potentially enhanced virulence mediated by an mycP1 mutation. The Beijing sub-lineage strain demonstrated a high degree of co-resistance between isoniazid and ethionamide (80/95; 84.2 %) attributed to an inhA promoter mutation combined with inhA and ndh coding mutations. Multi-drug resistance, observed in 78/95 samples, emerged with the acquisition of a typical rpoB mutation together with a compensatory rpoC mutation in the 1980s. There was independent acquisition of fluoroquinolone and aminoglycoside resistance, and evidence of local transmission of extensively drug resistant (XDR) strains from 2009. These findings underline the importance of whole genome sequencing in informing an effective public health response to MDR/XDR TB.

Methodological and Clinical Aspects of the Molecular Epidemiology of Mycobacterium tuberculosis and Other Mycobacteria

Molecular typing has revolutionized epidemiological studies of infectious diseases, including those of a mycobacterial etiology. With the advent of fingerprinting techniques, many traditional concepts regarding transmission, infectivity, or pathogenicity of mycobacterial bacilli have been revisited, and their conventional interpretations have been challenged. Since the mid-1990s, when the first typing methods were introduced, a plethora of other modalities have been proposed. So-called molecular epidemiology has become an essential subdiscipline of modern mycobacteriology. It serves as a resource for understanding the key issues in the epidemiology of tuberculosis and other mycobacterial diseases. Among these issues are disclosing sources of infection, quantifying recent transmission, identifying transmission links, discerning reinfection from relapse, tracking the geographic distribution and clonal expansion of specific strains, and exploring the genetic mechanisms underlying specific phenotypic traits, including virulence, organ tropism, transmissibility, or drug resistance. Since genotyping continues to unravel the biology of mycobacteria, it offers enormous promise in the fight against and prevention of the diseases caused by these pathogens. In this review, molecular typing methods for Mycobacterium tuberculosis and nontuberculous mycobacteria elaborated over the last 2 decades are summarized. The relevance of these methods to the epidemiological investigation, diagnosis, evolution, and control of mycobacterial diseases is discussed.

Whole Genome Sequencing Demonstrates Limited Transmission within Identified Mycobacterium tuberculosis Clusters in New South Wales, Australia

Australia has a low tuberculosis incidence rate with most cases occurring among recent immigrants. Given suboptimal cluster resolution achieved with 24-locus mycobacterium interspersed repetitive unit (MIRU-24) genotyping, the added value of whole genome sequencing was explored. MIRU-24 profiles of all Mycobacterium tuberculosis culture-confirmed tuberculosis cases diagnosed between 2009 and 2013 in New South Wales (NSW), Australia, were examined and clusters identified. The relatedness of cases within the largest MIRU-24 clusters was assessed using whole genome sequencing and phylogenetic analyses. Of 1841 culture-confirmed TB cases, 91.9% (1692/1841) had complete demographic and genotyping data. East-African Indian (474; 28.0%) and Beijing (470; 27.8%) lineage strains predominated. The overall rate of MIRU-24 clustering was 20.1% (340/1692) and was highest among Beijing lineage strains (35.7%; 168/470). One Beijing and three East-African Indian (EAI) clonal complexes were responsible for the majority of observed clusters. Whole genome sequencing of the 4 largest clusters (30 isolates) demonstrated diverse single nucleotide polymorphisms (SNPs) within identified clusters. All sequenced EAI strains and 70% of Beijing lineage strains clustered by MIRU-24 typing demonstrated distinct SNP profiles. The superior resolution provided by whole genome sequencing demonstrated limited M. tuberculosis transmission within NSW, even within identified MIRU-24 clusters. Routine whole genome sequencing could provide valuable public health guidance in low burden settings.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

Genetic diversity of Mycobacterium tuberculosis from Guadalajara, Mexico and identification of a rare multidrug resistant Beijing genotype

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

Identifying locations of recent TB transmission in rural Uganda: a multidisciplinary approach

OBJECTIVES

Targeting high Tuberculosis (TB) transmission sites may offer a novel approach to TB prevention in sub-Saharan Africa. We sought to characterise TB transmission sites in a rural Ugandan township.

METHODS

We recruited adults starting TB treatment in Tororo, Uganda, over 1 year. Fifty four TB cases provided names of frequent contacts, sites of residence, health care, work and social activities, and two sputum samples. Mycobacterium tuberculosis (MTB) culture-positive specimens underwent spoligotyping to identify strains with shared genotypes. We visualised TB case social networks, and obtained, mapped and geo-coded global positioning system measures for every location that cases reported frequenting 1 month before treatment. Locations of spatial overlap among genotype-clustered cases were considered potential transmission sites.

RESULTS

Six distinct genotypic clusters were identified involving 21 of 33 (64%) MTB culture-positive, genotyped cases; none shared a home. Although 18 of 54 (33%) TB cases shared social network ties, none of the genotype-clustered cases shared social ties. Using spatial analysis, we identified potential sites of within-cluster TB transmission for five of six genotypic clusters. All sites but one were healthcare and social venues, including sites of drinking, worship and marketplaces. Cases reported spending the largest proportion of pre-treatment person-time (22.4%) at drinking venues.

CONCLUSIONS

Using molecular epidemiology, geospatial and social network data from adult TB cases identified at clinics, we quantified person-time spent at high-risk locations across a rural Ugandan community and determined the most likely sites of recent TB transmission to be healthcare and social venues. These sites may not have been identified using contact investigation alone.

Effect of study design and setting on tuberculosis clustering estimates using Mycobacterial Interspersed Repetitive Units-Variable Number Tandem Repeats (MIRU-VNTR): a systematic review

OBJECTIVES

To systematically review the evidence for the impact of study design and setting on the interpretation of tuberculosis (TB) transmission using clustering derived from Mycobacterial Interspersed Repetitive Units-Variable Number Tandem Repeats (MIRU-VNTR) strain typing.

DATA SOURCES

MEDLINE, EMBASE, CINHAL, Web of Science and Scopus were searched for articles published before 21st October 2014.

REVIEW METHODS

Studies in humans that reported the proportion of clustering of TB isolates by MIRU-VNTR were included in the analysis. Univariable meta-regression analyses were conducted to assess the influence of study design and setting on the proportion of clustering.

RESULTS

The search identified 27 eligible articles reporting clustering between 0% and 63%. The number of MIRU-VNTR loci typed, requiring consent to type patient isolates (as a proxy for sampling fraction), the TB incidence and the maximum cluster size explained 14%, 14%, 27% and 48% of between-study variation, respectively, and had a significant association with the proportion of clustering.

CONCLUSIONS

Although MIRU-VNTR typing is being adopted worldwide there is a paucity of data on how study design and setting may influence estimates of clustering. We have highlighted study design variables for consideration in the design and interpretation of future studies.

IS6110-5'3'FP: an automated typing approach for Mycobacterium tuberculosis complex strains simultaneously targeting and resolving IS6110 5' and 3' polymorphisms

OBJECTIVES

Fingerprinting of Mycobacterium tuberculosis complex strains based on the IS6110 insertion sequence would considerably gain in terms of discriminatory power and versatility if both 5' and 3' polymorphisms were simultaneously targeted, and if it benefited from automated capillary electrophoresis. In response to these requirements, we developed IS6110-5'3'FP (IS6110 5' and 3' fluorescent polymorphisms).

METHODS

IS6110-5'3'FP involves the construction of an M. tuberculosis genomic library in a plasmid vector using HincII endonuclease, which cuts within the IS6110 sequence. After amplification in Escherichia coli, the library is subjected to selective and simultaneous PCR amplification of IS6110 5' and 3' polymorphic fragments, using differentially labeled fluorescent primers. The resulting amplicons are then fractionated on a capillary sequencer and the signal peaks analyzed as digital data.

RESULTS

IS6110-5'3'FP consistently detected and resolved both 5' and 3' IS6110 polymorphic fragments (35% and 65%, respectively) with a high level of reproducibility. The method differentiated all M. tuberculosis strains, as did IS6110 restriction fragment length polymorphism (RFLP), the gold standard of IS6110-based typing. Strikingly, the potential of IS6110-5'3'FP to resolve more polymorphic fragments than IS6110 RFLP was demonstrated.

CONCLUSIONS

IS6110-5'3'FP demonstrated sufficient potential to be a promising automated alternative to IS6110 RFLP, amenable to high throughput analysis and inter-laboratory comparison.

Temporal dynamics of Mycobacterium tuberculosis genotypes in New South Wales, Australia

BACKGROUND

Molecular epidemiology of Mycobacterium tuberculosis, its transmission dynamics and population structure have become important determinants of targeted tuberculosis control programs. Here we describe recent changes in the distribution of M. tuberculosis genotypes in New South Wales (NSW), Australia and compared strain types with drug resistance, site of disease and demographic data.

METHODS

We evaluated all culture-confirmed newly identified tuberculosis cases in NSW, Australia, from 2010-2012. M. tuberculosis population structure and clustering rates were assessed using 24-loci Mycobacterial interspersed repetitive unit (MIRU) analysis and compared to MIRU data from 2006-2008.

RESULTS

Of 1177 tuberculosis cases, 1128 (95.8%) were successfully typed. Beijing and East African Indian (EAI) lineage strains were most common (27.6% and 28.5%, respectively) with EAI strains increasing in relative abundance from 11.8% in 2006-2008 to 28.5% in 2010-2012. Few cases of multi-drug resistant tuberculosis were identified (18; 1.7%). Compared to 12-loci, 24-loci MIRU provided improved cluster resolution with 695 (61.6%) and 227 (20.1%) clustered cases identified, respectively. Detailed analysis of the largest cluster identified (an 11 member Beijing cluster) revealed wide geographic diversity in the absence of documented social contact.

CONCLUSIONS

EAI strains of M. tuberculosis recently overtook Beijing family as a prevalent cause of tuberculosis in NSW, Australia. This lineage appeared to be less commonly related to multi-drug resistant tuberculosis as compared to Beijing strain lineage. The resolution provided by 24-loci MIRU typing was insufficient for reliable assessment of transmissions, especially of Beijing family strains.

Whole-genome sequencing to detect recent transmission of Mycobacterium tuberculosis in settings with a high burden of tuberculosis

Whole genome sequencing (WGS) of Mycobacterium tuberculosis has been used to trace the transmission of M. tuberculosis, the causative agent of tuberculosis (TB). Previously published studies using WGS were conducted in developed countries with a low TB burden. We sought to evaluate the relative usefulness of traditional VNTR and SNP typing methods, WGS and epidemiological investigations to study the recent transmission of M. tuberculosis in a high TB burden country. We conducted epidemiological investigations of 42 TB patients whose M. tuberculosis isolates were classified into three clusters based on variable-number tandem repeat (VNTR) typing. We applied WGS to 32 (76.2%) of the 42 strains and calculated the pairwise genomic distances between strains within each cluster. Eighteen (56.3%) of the 32 strains had genomic differences ≥100 SNPs with every other strain, suggesting that direct transmission did not likely occurred. Ten strains were grouped into four WGS-based clusters with genomic distances ≤5 SNPs within each cluster, and confirmed epidemiological links were identified in two of these clusters. Our results indicate that WGS provides reliable resolution for tracing the transmission of M. tuberculosis in high TB burden settings. The high resolution of WGS is particularly useful to confirm or exclude the possibility of direct transmission events defined by traditional typing methods.

Comparison between RFLP and MIRU-VNTR genotyping of Mycobacterium tuberculosis strains isolated in Stockholm 2009 to 2011

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

Development of a hierarchical variable-number tandem repeat typing scheme for Mycobacterium tuberculosis in China

Molecular typing based on variable-number tandem repeats (VNTR) analysis is a promising tool for identifying transmission of Mycobacterium tuberculosis. However, the currently proposed 15- and 24-locus VNTR sets (VNTR-15/24) only have limited resolution and contain too many loci for large-scale typing in high burden countries. To develop an optimal typing scheme in China, we evaluated the resolution and robustness of 25 VNTR loci, using population-based collections of 1362 clinical isolates from six provinces across the country. The resolution of most loci showed considerable variations among regions. By calculating the average resolution of all possible combinations of 20 robust loci, we identified an optimal locus set with a minimum of 9 loci (VNTR-9) that could achieve comparable resolution of the standard VNTR-15. The VNTR-9 had consistently high resolutions in all six regions, and it was highly concordant with VNTR-15 for defining both clustered and unique genotypes. Furthermore, VNTR-9 was phylogenetically informative for classifying lineages/sublineages of M. tuberculosis. Three hypervariable loci (HV-3), VNTR 3232, VNTR 3820 and VNTR 4120, were proved important for further differentiating unrelated clustered strains based on VNTR-9. We propose the optimized VNTR-9 as first-line method and the HV-3 as second-line method for molecular typing of M. tuberculosis in China and surrounding countries. The development of hierarchical VNTR typing methods that can achieve high resolution with a small number of loci could be suitable for molecular epidemiology study in other high burden countries.

Molecular epidemiology of Mycobacterium tuberculosis complex in Singapore, 2006-2012

Background

Tuberculosis remains common in Singapore, increasing in incidence since 2008. We attempted to determine the molecular epidemiology of Mycobacterium tuberculosis complex (MTC) isolates locally, identifying major circulating genotypes and obtaining a glimpse of transmission dynamics.

Methodology

Non-duplicate MTC isolates archived between 2006 and 2012 at the larger clinical tuberculosis laboratory in Singapore were sampled for spoligotyping and MIRU-VNTR typing, with case data obtained from the Singapore Tuberculosis Elimination Program registry database. Isolates between 2008 and 2012 were selected because of either multidrug-resistance or potential epidemiological linkage, whereas earlier isolates were randomly selected. Separate analyses were performed for the early (2006-2007) and later (2008-2012) study phases in view of potential selection bias.

Principal Findings

A total of 1,612 MTC isolates were typed, constituting 13.1% of all culture-positive tuberculosis cases during this period. Multidrug-resistance was present in 91 (5.6%) isolates – higher than the national prevalence in view of selection bias. The majority of isolates belonged to the Beijing (45.8%) and EAI (22.8%) lineages. There were 347 (30.7%) and 133 (27.5%) cases clustered by combined spoligotyping and MIRU-VNTR typing from the earlier and later phases respectively. Patients within these clusters tended to be of Chinese ethnicity, Singapore resident, and have isolates belonging to the Beijing lineage. A review of prior contact investigation results for all patients with clustered isolates failed to reveal epidemiological links for the majority, suggesting either unknown transmission networks or inadequate specificity of the molecular typing methods in a country with a moderate incidence of tuberculosis.

Conclusion

Our work demonstrates that Singapore has a large and heterogeneous distribution of MTC strains, and with possible cross-transmission over the past few years based on our molecular typing results. A universal MTC typing program coupled with enhanced contact investigations may be useful in further understanding the transmission dynamics of tuberculosis locally.

Proposal of a consensus set of hypervariable mycobacterial interspersed repetitive-unit-variable-number tandem-repeat loci for subtyping of Mycobacterium tuberculosis Beijing isolates

Mycobacterium tuberculosis Beijing strains represent targets of special importance for molecular surveillance of tuberculosis (TB), especially because they are associated with spread of multidrug resistance in some world regions. Standard 24-locus mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing lacks resolution power for accurately discriminating closely related clones that often compose Beijing strain populations. Therefore, we evaluated a set of 7 additional, hypervariable MIRU-VNTR loci for better resolution and tracing of such strains, using a collection of 535 Beijing isolates from six world regions where these strains are known to be prevalent. The typeability and interlaboratory reproducibility of these hypervariable loci were lower than those of the 24 standard loci. Three loci (2163a, 3155, and 3336) were excluded because of their redundant variability and/or more frequent noninterpretable results compared to the 4 other markers. The use of the remaining 4-locus set (1982, 3232, 3820, and 4120) increased the number of types by 52% (from 223 to 340) and reduced the clustering rate from 58.3 to 36.6%, when combined with the use of the standard 24-locus set. Known major clonal complexes/24-locus-based clusters were all subdivided, although the degree of subdivision varied depending on the complex. Only five single-locus variations were detected among the hypervariable loci of an additional panel of 92 isolates, representing 15 years of clonal spread of a single Beijing strain in a geographically restricted setting. On this calibrated basis, we propose this 4-locus set as a consensus for subtyping Beijing clonal complexes and clusters, after standard typing.

Mycobacterium tuberculosis Beijing genotype is associated with HIV infection in Mozambique

The Beijing genotype is a lineage of Mycobacterium tuberculosis that is distributed worldwide and responsible for large epidemics, associated with multidrug-resistance. However, its distribution in Africa is less understood due to the lack of data. Our aim was to investigate the prevalence and possible transmission of Beijing strains in Mozambique by a multivariate analysis of genotypic, geographic and demographic data. A total of 543 M. tuberculosis isolates from Mozambique were spoligotyped. Of these, 33 were of the Beijing lineage. The genetic relationship between the Beijing isolates were studied by identification of genomic deletions within some Regions of Difference (RD), Restriction Fragment Length Polymorphism (RFLP) and Mycobacterial Interspersed Repetivie Unit – variable number tandem repeat (MIRU-VNTR). Beijing strains from South Africa, representing different sublineages were included as reference strains. The association between Beijing genotype, Human Immunodeficiency Virus (HIV) serology and baseline demographic data was investigated. HIV positive serostatus was significantly (p=0.023) more common in patients with Beijing strains than in patients with non-Beijing strains in a multivariable analysis adjusted for age, sex and province (14 (10.9%) of the 129 HIV positive patients had Beijing strains while 6/141 (4.3%) of HIV negative patients had Beijing strains). The majority of Beijing strains were found in the Southern region of Mozambique, particularly in Maputo City (17%). Only one Beijing strain was drug resistant (multi-drug resistant). By combined use of RD and spoligotyping, three genetic sublineages could be tentatively identified where a distinct group of four isolates had deletion of RD150, a signature of the “sublineage 7” recently emerging in South Africa. The same group was very similar to South African “sublineage 7” by RFLP and MIRU-VNTR, suggesting that this sublineage could have been recently introduced in Mozambique from South Africa, in association with HIV infection.

First molecular epidemiology study of Mycobacterium tuberculosis in Kiribati

Tuberculosis incidence rates in Kiribati are among the highest in the Western Pacific Region, however the genetic diversity of circulating Mycobacterium tuberculosis complex strains (MTBC) and transmission dynamics are unknown. Here, we analysed MTBC strains isolated from culture positive pulmonary tuberculosis (TB) cases from the main TB referral centre between November 2007 and October 2009. Strain genotyping (IS6110 typing, spoligotyping, 24-loci MIRU-VNTR and SNP typing) was performed and demographic information collected. Among 73 MTBC strains analysed, we identified seven phylogenetic lineages, dominated by Beijing strains (49%). Beijing strains were further differentiated in two main branches, Beijing-A (n = 8) and -B (n = 28), that show distinct genotyping patterns and are characterized by specific deletion profiles (Beijing A: only RD105, RD207 deleted; Beijing B: RD150 and RD181 additionally deleted). Many Kiribati strains (59% based on IS6110 typing of all strains) occurred in clusters, suggesting ongoing local transmission. Beijing-B strains and over-crowded living conditions were associated with strain clustering (likely recent transmission), however little evidence of anti-tuberculous drug resistance was observed. We suggest enhanced case finding amongst close contacts and continued supervised treatment of all identified cases using standard first-line drugs to reduce TB burden in Kiribati. Beijing strains can be subdivided in different principle branches that might be associated with differential spreading patterns in the population.

Comparative study of IS6110 restriction fragment length polymorphism and variable-number tandem-repeat typing of Mycobacterium tuberculosis isolates in the Netherlands, based on a 5-year nationwide survey

In order to switch from IS6110 and polymorphic GC-rich repetitive sequence (PGRS) restriction fragment length polymorphism (RFLP) to 24-locus variable-number tandem-repeat (VNTR) typing of Mycobacterium tuberculosis complex isolates in the national tuberculosis control program in The Netherlands, a detailed evaluation on discriminatory power and agreement with findings in a cluster investigation was performed on 3,975 tuberculosis cases during the period of 2004 to 2008. The level of discrimination of the two typing methods did not differ substantially: RFLP typing yielded 2,733 distinct patterns compared to 2,607 in VNTR typing. The global concordance, defined as isolates labeled unique or identically distributed in clusters by both methods, amounted to 78.5% (n = 3,123). Of the remaining 855 cases, 12% (n = 479) of the cases were clustered only by VNTR, 7.7% (n = 305) only by RFLP typing, and 1.8% (n = 71) revealed different cluster compositions in the two approaches. A cluster investigation was performed for 87% (n = 1,462) of the cases clustered by RFLP. For the 740 cases with confirmed or presumed epidemiological links, 92% were concordant with VNTR typing. In contrast, only 64% of the 722 cases without an epidemiological link but clustered by RFLP typing were also clustered by VNTR typing. We conclude that VNTR typing has a discriminatory power equal to IS6110 RFLP typing but is in better agreement with findings in a cluster investigation performed on an RFLP-clustering-based cluster investigation. Both aspects make VNTR typing a suitable method for tuberculosis surveillance systems.

Molecular characterization of Mycobacterium tuberculosis isolates from North Indian patients with extrapulmonary tuberculosis

Genotypic studies are important to understand the molecular epidemiology and transmission routes of Mycobacterium tuberculosis. In the first and largest study from India, spoligotyping and 24 loci mycobacterial interspersed repetitive units (MIRU) were performed to find genetic profiles of 125 M. tuberculosis strains isolated from patients with extrapulmonary tuberculosis (EPTB) and their drug susceptibility test was performed using BACTEC-MGIT 960. Spoligotyping results were compared with the world Spoligotyping Database of Institute Pasteur de Guadeloupe (SpolDB4). The spoligotyping results showed that 110 (88%) displayed known patterns while 15 (12%) isolates had no matching database. Predominant spoligotypes belonged to CAS family (57.27%). The largest clade comprised of 38 isolates belonging to the CAS1_DEL lineage. Though there was no significant association between specific mycobacterial lineage and extrapulmonary site, a significantly high (p < 0.001) number of Beijing type isolates (28.6%) were isolated from bone and joint samples as compared to cerebrospinal fluid (5%). There was a significant association between Beijing family isolates and multi-drug-resistance, while all MANU genotypes were pan-drug sensitive. The CAS family lineage was most prevalent genotype in the EPTB cases in our population.

Combination of single nucleotide polymorphism and variable-number tandem repeats for genotyping a homogenous population of Mycobacterium tuberculosis Beijing strains in China

The standard 15- and 24-locus variable-number tandem repeat (VNTR) genotyping methods have demonstrated adequate discriminatory power and a small homoplasy effect for tracing tuberculosis (TB) transmission and predicting Mycobacterium tuberculosis lineages in European and North American countries. However, its validity for the definition of transmission in homogenous M. tuberculosis populations in settings with high TB burdens has been questioned. Here, we genotyped a population-based collection of 191 Beijing strains based on standard 15-locus VNTR (VNTR-15) and 8 single nucleotide polymorphisms (SNPs) in Shanghai, China. Limited discriminatory power and high rates of VNTR homoplasy were observed in the homogenous population of evolutionarily "modern" Beijing strains. Additional typing of three hypervariable loci (VNTR3820, VNTR4120, and VNTR3232) was performed for VNTR-15-based clusters. High variations of hypervariable alleles were observed in clusters with inconsistent SNP sublineages. We concluded that SNPs and hypervariable VNTR loci are helpful to enhance the discriminatory power and decrease the VNTR homoplasy effect for defining clusters. We recommend the combination of standard VNTR-15 and SNPs as first-line typing methods and the hypervariable loci for second-line typing of clustered strains for molecular epidemiology studies of homogenous M. tuberculosis populations.

DNA fingerprinting of Mycobacterium tuberculosis: from phage typing to whole-genome sequencing

Current typing methods for Mycobacterium tuberculosis complex evolved from simple phenotypic approaches like phage typing and drug susceptibility profiling to DNA-based strain typing methods, such as IS6110-restriction fragment length polymorphisms (RFLP) and variable number of tandem repeats (VNTR) typing. Examples of the usefulness of molecular typing are source case finding and epidemiological linkage of tuberculosis (TB) cases, international transmission of MDR/XDR-TB, the discrimination between endogenous reactivation and exogenous re-infection as a cause of relapses after curative treatment of tuberculosis, the evidence of multiple M. tuberculosis infections, and the disclosure of laboratory cross-contaminations. Simultaneously, phylogenetic analyses were developed based on single nucleotide polymorphisms (SNPs), genomic deletions usually referred to as regions of difference (RDs) and spoligotyping which served both strain typing and phylogenetic analysis. National and international initiatives that rely on the application of these typing methods have brought significant insight into the molecular epidemiology of tuberculosis. However, current DNA fingerprinting methods have important limitations. They can often not distinguish between genetically closely related strains and the turn-over of these markers is variable. Moreover, the suitability of most DNA typing methods for phylogenetic reconstruction is limited as they show a high propensity of convergent evolution or misinfer genetic distances. In order to fully explore the possibilities of genotyping in the molecular epidemiology of tuberculosis and to study the phylogeny of the causative bacteria reliably, the application of whole-genome sequencing (WGS) analysis for all M. tuberculosis isolates is the optimal, although currently still a costly solution. In the last years WGS for typing of pathogens has been explored and yielded important additional information on strain diversity in comparison to the classical DNA typing methods. With the ongoing cost reduction of DNA sequencing it is possible that WGS will become the sole diagnostic tool in the secondary laboratory diagnosis of tuberculosis for identification, drug susceptibility testing and genetic characterization.

Evaluation of Mycobacterium tuberculosis typing methods in a 4-year study in Schleswig-Holstein, Northern Germany

In order to evaluate the discriminatory power of different methods for genotyping of Mycobacterium tuberculosis complex (MTBC) isolates, we compared the performance of (i) IS6110 DNA fingerprint typing, (ii) spoligotyping, and (iii) 24-loci mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) typing in a long-term study on the epidemiology of tuberculosis (TB) in Schleswig-Holstein, the northernmost federal state of Germany. In total, we analyzed 277 MTBC isolates collected from patients between the years 2006 and 2010. The collection comprised a broad spectrum of 13 different genotypes, among which strains of the Haarlem genotype (31%) were most prominent, followed by strains belonging to the Delhi and Beijing lineages (7% and 6%, respectively). On the basis of IS6110 restriction fragment length polymorphism (RFLP) and spoligotyping analyses, 211 isolates had unique patterns (76%) and 66 isolates (24%) were in 20 clusters. MIRU-VNTR combined with spoligotyping analyses revealed 202 isolates with unique patterns (73%) and 75 isolates in 18 clusters (27%). Overall, there was 93.1% concordance between the typing results obtained; 198 strains were identified as unique, and 60 isolates were clustered by both typing combinations (including all 31 isolates with confirmed epidemiological links). Of the remaining 19 isolates with discrepant results, 15 were falsely clustered by MIRU-VNTR (six Beijing genotype strains) and four were clustered by IS6110 RFLP (low IS6110 copy number) only. In conclusion, in the study population investigated, a minority of isolates, especially of the Beijing genotype, clustered by standard 24-loci MIRU-VNTR and without an obvious epidemiological link may require second-line typing by IS6110 RFLP or hypervariable MIRU-VNTR loci.

An appraisal of PCR-based technology in the detection of Mycobacterium tuberculosis

Tuberculosis is an under-recognized yet catastrophic health problem, particularly in developing countries. The HIV pandemic has served to increase the number of susceptible individuals, and multidrug-resistance and poor socioeconomic conditions also augment the prevalence and the consequences of the disease. To control the disease and its spread, it is vital that tuberculosis diagnostics are accurate and rapid. Whereas microscopy and culture have several limitations (low sensitivity is a problem for the former, while the latter has a delayed turnaround time), PCR-based techniques targeting regions of the Mycobacterium tuberculosis genome such as IS6110 have proved to be useful. The purpose of this review is to assess the use of PCR-RFLP, nested PCR and real-time PCR protocols and the choice of target regions for the detection of M. tuberculosis. Real-time PCR for the detection of M. tuberculosis target genes in clinical specimens has contributed to improving diagnosis and epidemiologic surveillance in the past decade. However, targeting one genome sequence such as IS6110 may not by itself be sufficiently sensitive to reach 100% diagnosis, especially in the case of pulmonary tuberculosis. Additional testing for target genome sequences such as hsp65 seems encouraging. An interesting approach would be a multiplex real-time PCR targeting both IS6110 and hsp65 to achieve comprehensive and specific molecular diagnosis. This technology needs development and adequate field testing before it becomes the acceptable gold standard for diagnosis.

Application of single-nucleotide polymorphism and mycobacterial interspersed repetitive units-variable number of tandem repeats analyses to clinical Mycobacterium tuberculosis isolates from Korea

Background

Single-nucleotide polymorphism (SNP) analysis is a powerful strategy for large-scale molecular population studies examining phylogenetic relationships among bacterial strains. Mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) can be easily digitized to share data among laboratories. This study applied SNP and MIRU-VNTR analyses for molecular strain typing of Mycobacterium tuberculosis isolates collected throughout Korea.

Methods

We studied 102 clinical M. tuberculosis isolates, including 6 paired strains, collected from 11 university hospitals in Korea in 2008 and 2009. SNPs were detected using hairpin primer assays, and then, MIRU-VNTR analysis was performed.

Results

Thirty-five SNPs contained polymorphisms that helped differentiate the 96 tested isolates. The isolates were classified into 15 clusters. The Beijing family strains were distributed within closely related clusters in the SNP dendrogram. For MIRU-VNTR analysis, the 96 isolates were divided into 12 groups. The discriminatory index in 8 of these groups (MIRU-10, -23, -26, and -31; ETR-A, -B, -C, and -F) was high (Hunter-Gaston diversity index > 0.6). Unlike the SNP method, MIRU-VNTR analysis did not identify any notable localizations of Beijing or non-Beijing family isolates in specific clusters.

Conclusions

SNP and MIRU-VNTR analyses are surrogate molecular strain-typing methods for M. tuberculosis in Korea where Beijing family isolates are predominant.

Genotyping of Mycobacterium tuberculosis: application in epidemiologic studies

Genotyping is used to track specific isolates of Mycobacterium tuberculosis in a community. It has been successfully used in epidemiologic research (termed 'molecular epidemiology') to study the transmission dynamics of TB. In this article, we review the genetic markers used in molecular epidemiologic studies including the use of whole-genome sequencing technology. We also review the public health application of molecular epidemiologic tools.

A novel method for the rapid and prospective identification of Beijing Mycobacterium tuberculosis strains by high-resolution melting analysis

Genotypic analysis of Mycobacterium tuberculosis (MTB) has enabled the definition of several lineages. The Beijing family, which is considered highly virulent and transmissible, has been associated with resistance in certain settings and involved in severe outbreaks, making it one of the most closely-monitored lineages. Therefore, rapid prospective identification of Beijing MTB strains could be relevant. In the present study, we evaluate a real-time PCR followed by high-resolution melting (HRM) based on the identification of a single nucleotide polymorphism (SNP) in the Rv2629 gene which defines Beijing lineage (A191C for Beijing genotype and A191A for non-Beijing genotype). This combined methodology efficiently differentiated Beijing and non-Beijing strains in 100% of the isolates from a collection of reference strains without requiring specific DNA probes. Additionally, HRM was able to assign a Beijing/non-Beijing genotype in 90.9% of the respiratory specimens assayed. Its applicability was tested on a Peruvian sample of circulating MTB strains, in which it identified 10.7% as belonging to the Beijing genotype; this proportion reached 20% in the North Lima area. HRM analysis of the A191C SNP is a rapid, reliable, and sensitive method for the efficient prospective survey of high-risk Beijing MTB strains, even in developing settings where MTB culture is often not available.

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.

Mechanisms behind variation in the Clostridium difficile 16S-23S rRNA intergenic spacer region

Clostridium difficile infection is an increasing problem in hospitals worldwide, mainly due to the recent emergence of a hypervirulent C. difficile strain. C. difficile PCR ribotyping, based on size variation of the 16S-23S rRNA intergenic spacer region (16S-23S ISR), is widely used in Europe for molecular epidemiological investigation. The mechanism underlying the 16S-23S ISR size variations in the genome of C. difficile is currently not completely understood. To elucidate this mechanism, isolates of six different PCR ribotypes were analysed by cloning and sequencing the 16S-23S ISR. A direct repeat, IB, of 9&emsp14;bp was detected up to five times in the 16S-23S ISR in all 47 clones investigated. Thirty-five clones displayed differences either by ribotype or by nucleotide sequence. The sequences of the 16S-23S ISR of C. difficile showed a uniformly organized structure, composed of a tRNA(Ala) gene and spacers of 33 and 53&emsp14;bp separated by the 9&emsp14;bp direct repeat IB. The results of the study support the hypothesis that this composition is responsible for the length variations seen in the 16S-23S ISR, and indicate that these length variations result from slipped-strand mispairing and intra- and possibly interchromosomal homologous recombination.

Microevolution of Mycobacterium tuberculosis in a tuberculosis patient

Five Mycobacterium tuberculosis isolates were obtained from three body sites from a Dutch patient. The isolates displayed a single genotype by 24-locus MIRU-VNTR typing (except for a single locus not amplified from one isolate) but were differentiated by small variations in IS6110 fingerprints, spoligotypes, 6 hypervariable MIRU-VNTR loci, and/or DiversiLab profiles, revealing patterns of microevolution in a clonal infection.

Drug resistant Mycobacterium tuberculosis of the Beijing genotype does not spread in Sweden

BACKGROUND

Drug resistant (DR) and multi-drug resistant (MDR) tuberculosis (TB) is increasing worldwide. In some parts of the world 10% or more of new TB cases are MDR. The Beijing genotype is a distinct genetic lineage of Mycobacterium tuberculosis, which is distributed worldwide, and has caused large outbreaks of MDR-TB. It has been proposed that certain lineages of M. tuberculosis, such as the Beijing lineage, may have specific adaptive advantages. We have investigated the presence and transmission of DR Beijing strains in the Swedish population.

METHODOLOGY/PRINCIPAL FINDINGS

All DR M. tuberculosis complex isolates between 1994 and 2008 were studied. Isolates that were of Beijing genotype were investigated for specific resistance mutations and phylogenetic markers. Seventy (13%) of 536 DR strains were of Beijing genotype. The majority of the patients with Beijing strains were foreign born, and their country of origin reflects the countries where the Beijing genotype is most prevalent. Multidrug-resistance was significantly more common in Beijing strains than in non-Beijing strains. There was a correlation between the Beijing genotype and specific resistance mutations in the katG gene, the mabA-inhA-promotor and the rpoB gene. By a combined use of RD deletions, spoligotyping, IS1547, mutT gene polymorphism and Rv3135 gene analysis the Beijing strains could be divided into 11 genomic sublineages. Of the patients with Beijing strains 28 (41%) were found in altogether 10 clusters (2-5 per cluster), as defined by RFLP IS6110, while 52% of the patients with non-Beijing strains were in clusters. By 24 loci MIRU-VNTR 31 (45%) of the patients with Beijing strains were found in altogether 7 clusters (2-11 per cluster). Contact tracing established possible epidemiological linkage between only two patients with Beijing strains.

CONCLUSIONS/SIGNIFICANCE

Although extensive outbreaks with non-Beijing TB strains have occurred in Sweden, Beijing strains have not taken hold, in spite of the proximity to high prevalence countries such as Russia and the Baltic countries. The Beijing sublineages so far introduced in Sweden may not be adapted to spread in the Scandinavian population.

Evaluation of 24 locus MIRU-VNTR genotyping of Mycobacterium tuberculosis isolates in Canada

The current gold standard for Mycobacterium tuberculosis complex (MTBC) genotyping is insertion sequence (IS) 6110 restriction fragment length polymorphism (RFLP) as it provides the highest discriminatory power of all available MTBC genotyping methods. However, RFLP is labour intensive and the interpretation of data from this method can be susceptible to errors. In 2001 a rapid, reproducible variable number of tandem repeat (VNTR) based typing method using 12 mycobacterial interspersed repetitive units (MIRU) was developed. Despite this advancement, this method lacked the discriminatory power of IS6110-RFLP. More recently a set of 24 MIRU-VNTR loci was reported to have greater discriminatory power than the original 12 locus system and may exceed that of RFLP when combined with spoligotyping. We compared the 24 locus method to the 12 locus method in order to improve surveillance of tuberculosis in Canada. A random sample of 650 MTBC isolates from British Columbia, Saskatchewan, Manitoba and Quebec Canada was genotyped using the 24 MIRU loci. Comparison of the data for the 12 and 24 MIRU loci showed an increase of the Hunter-Gaston discriminatory index (HGDI) from 0.895 (12 loci) to 0.920 (24 loci). The implementation of the 24 locus MIRU-VNTR methods offers improvement in discriminatory power over the traditional 12 locus method. For long-term surveillance of MTBC within Canada, the use of 24 MIRU-VNTR loci will provide rapid, highly discriminatory molecular epidemiology information.

Genetic diversity of Mycobacterium tuberculosis isolates from a tertiary care tuberculosis hospital in South Korea

Tuberculosis (TB) remains an immense public health problem in the Republic of Korea despite a more than fivefold decrease in the prevalence of the disease over the last 3 decades. The rise in drug-resistant TB has compounded the situation. We analyzed 208 clinical isolates of M. tuberculosis from the National Masan Tuberculosis Hospital by spoligotyping, IS6110 restriction fragment length polymorphism (RFLP), and 24-locus-based mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) typing to assess the diversity and transmission dynamics of the tubercle bacilli in the Republic of Korea. The majority of the isolates (97.1%) belonged to the Beijing genotype. Cluster analysis by MIRU-VNTR yielded a low clustering rate of 22.3%, with most of the clusters comprising isolates with diverse drug resistance patterns. The discriminatory capacity of the typing methods was high for RFLP and MIRU-VNTR (allelic diversity [h] = 0.99) but low for spoligotyping (h = 0.31). Although analysis of 19 MIRU-VNTR loci was needed to achieve maximum discrimination, an informative set of 8 loci (960, 1955, 2163b, 2165, 2996, 3192, 4052, and 4348) (h = 0.98) that was able to differentiate most of the closely related strains was identified. These findings suggest that 24-locus-based MIRU-VNTR typing is a likely suitable alternative to RFLP to differentiate clinical isolates in this setting, which is dominated by M. tuberculosis Beijing strains. Within the study limits, our results also suggest that the problem of drug-resistant TB in the Republic of Korea may be largely due to acquired resistance as opposed to transmission.

Mycobacterium tuberculosis complex CRISPR genotyping: improving efficiency, throughput and discriminative power of 'spoligotyping' with new spacers and a microbead-based hybridization assay

The aims of the present study were to implement a microbead-based 'spoligotyping' technique and to evaluate improvements by the addition of a panel of 25 extra spacers that we expected to provide an increased resolution on principal genetic group 1 (PGG 1) strains. We confirmed the high sensitivity and reproducibility of the classical technique using the 43 spacer panel and we obtained perfect agreement between the membrane-based and the microbead-based techniques. We further demonstrated an increase in the discriminative power of an extended 68 spacer format for differentiation of PGG 1 clinical isolates, in particular for the East African-Indian clade. Finally, we define a limited yet highly informative reduced 10 spacer panel set which could offer a more cost-effective option for implementation in resource-limited countries and that could decrease the need for additional VNTR (variable number of tandem repeats) genotyping work in molecular epidemiological studies. We also present an economic analysis comparing membrane-based and microbead-based techniques.

Polymorphic exact tandem repeat A (PETRA): a newly defined lineage of mycobacterium tuberculosis in israel originating predominantly in Sub-Saharan Africa

As part of the Israel National Program for Prevention and Control of Tuberculosis, the molecular epidemiology of new tuberculosis cases is monitored. Prospective screening showed that about 20% of all new cases of culture-positive tuberculosis (43 of 222) in Israel in the year 2008 were caused by certain Mycobacterium tuberculosis strains of the central Asian (CAS) spoligotype lineage. The identity and similarity of these strains by mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing form a lineage we call PETRA for polymorphic at locus ETR A. The name PETRA was given to 79 strains we have found since the year 2000, because the largest number of strains with MIRU-VNTR profiles identical other than at locus A formed three groups, including 5 of 10 strains that had deleted the ETR A region from their genomes. No PETRA strain was found to be multiple drug resistant (resistant to both isoniazid and rifampin [rifampicin]). Most patients (75% [58 of 77 patients of known origin]) infected with PETRA were of sub-Saharan African origins. The genotypes associated with the 79 PETRA lineage strains presented in this paper suggest that the PETRA lineage is a large, major contributor to new tuberculosis cases in Israel.

High diversity of Mycobacterium tuberculosis genotypes in South Africa and preponderance of mixed infections among ST53 isolates

The reemergence of tuberculosis (TB) has become a major health problem worldwide, especially in Asia and Africa. Failure to combat this disease due to nonadherence or inappropriate drug regimens has selected for the emergence of multiple-drug-resistant (MDR) TB. The development of new molecular genotyping techniques has revealed the presence of mixed Mycobacterium tuberculosis infections, which may accelerate the emergence of drug-resistant strains. There are some studies describing the local distribution of circulating strains in South Africa, but to date, reports describing the frequency and distribution of M. tuberculosis genotypes, and specifically MDR genotypes, across the different provinces are limited. Thus, 252 isolates (of which 109 were MDR) from eight of the nine provinces of South Africa were analyzed by spoligotyping. Spoligotyping showed 10 different lineages, and ST53 (11.1%) and ST1 (10.3%) were the most frequent genotypes. Of the 75 different spoligopatterns observed, 20 (7.9%) were previously unreported. Analysis of the mycobacterial interspersed repetitive units of variable-number tandem repeats of the ST53 and ST1 isolates revealed that approximately 54% of the ST53 isolates were of mixed M. tuberculosis subpopulations. Drug resistance (defined as resistance to at least isoniazid and/or rifampin) could only be linked to a history of previous anti-TB treatment (adjusted odds ratio, 4.0; 95% confidence interval, 2.27 to 7.10; P = <0.0001). This study describes a high diversity of circulating genotypes in South Africa in addition to a high frequency of mixed M. tuberculosis subpopulations among the ST53 isolates. MDR TB in South Africa could not be attributed to the spread of any single lineage.

In-depth molecular characterization of Mycobacterium tuberculosis from New Delhi--predominance of drug resistant isolates of the 'modern' (TbD1) type

Background

India has the highest estimated burden of tuberculosis in the world, accounting for 21% of all tuberculosis cases world-wide. However, due to lack of systematic analysis using multiple markers the available information on the genomic diversity of Mycobacterium tuberculosis in India is limited.

Methodology/Principal Findings

Thus, 65 M. tuberculosis isolates from New Delhi, India were analyzed by spoligotyping, MIRU-VNTR, large deletion PCR typing and single nucleotide polymorphism analysis (SNP). The Central Asian (CAS) 1 _DELHI sub-lineage was the most prevalent sub-lineage comprising 46.2% (n = 30) of all isolates, with shared-type (ST) 26 being the most dominant genotype comprising 24.6% (n = 16) of all isolates. Other sub-lineages observed were: East-African Indian (EAI)-5 (9.2%, n = 6), EAI6_BGD1 (6.2%, n = 4), EAI3_IND, CAS and T1 with 6.2% each (n = 4 each), Beijing (4.6%, n = 3), CAS2 (3.1%, n = 2), and X1 and X2 with 1 isolate each. Genotyping results from five isolates (7.7%) did not match any existing spoligopatterns, and one isolate, ST124, belonged to an undefined lineage. Twenty-six percent of the isolates belonged to the TbD1+ PGG1 genogroup. SNP analysis of the pncA gene revealed a CAS-lineage specific silent mutation, S65S, which was observed for all CAS-lineage isolates (except two ST26 isolates) and in 1 orphan. Mutations in the pncA gene, conferring resistance to pyrazinamide, were observed in 15.4% of all isolates. Collectively, mutations in the rpoB gene, the katG gene and in both rpoB and katG genes, conferring resistance to rifampicin and isoniazid, respectively, were more frequent in CAS1_DELHI isolates compared to non-CAS_DELHI isolates (OR: 3.1, CI95% [1.11, 8.70], P = 0.045). The increased frequency of drug-resistance could not be linked to the patients' history of previous anti-tuberculosis treatment (OR: 1.156, CI95% [0.40, 3.36], P = 0.79). Fifty-six percent of all new tuberculosis patients had mutations in either the katG gene or the rpoB gene, or in both katG and rpoB genes.

Conclusion

CAS1_DELHI isolates circulating in New Delhi, India have a high frequency of mutations in the rpoB and katG genes. A silent mutation (S65S) in the pncA gene can be used as a putative genetic marker for CAS-lineage isolates.

Genetic diversity of the Mycobacterium tuberculosis Beijing family in East Asia revealed through refined population structure analysis

The Beijing/W family is the endemic lineage of Mycobacterium tuberculosis in East Asia: it has disseminated worldwide. To elucidate its genetic diversity in Japan, phylogenetic reconstruction was performed using 403 M. tuberculosis Beijing family clinical isolates. Variable number of tandem repeats analysis revealed the strains from Japan to be dispersed mainly among five subgroups in a phylogenetic tree. Interestingly, the genotypes of the strains from China and Mongolia were restricted mainly to a single branch; they exhibited high clonality. IS6110 insertion in the NTF region was also analyzed. The majority (78.6%) of Japanese isolates belonged to the ancient sublineage. The modern Beijing strains were observed to correspond to the branch containing the foreign strains, although the ancient Beijing strains were dispersed among the tree's other branches. Our results reflect the singular genetic diversity and the epidemiological pattern of Beijing M. tuberculosis in Japan.