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

Emergence of bedaquiline-resistant tuberculosis and of multidrug-resistant and extensively drug-resistant Mycobacterium tuberculosis strains with rpoB Ile491Phe mutation not detected by Xpert MTB/RIF in Mozambique: a retrospective observational study

BACKGROUND

In 2021, an estimated 4800 people developed rifampicin-resistant tuberculosis in Mozambique, 75% of which went undiagnosed. Detailed molecular data on rifampicin-resistant and multidrug-resistant (MDR) tuberculosis are not available. Here, we aimed at gaining precise data on the determinants of rifampicin-resistant and MDR tuberculosis in Mozambique.

METHODS

In this retrospective observational study, we performed whole-genome sequencing of 704 rifampicin-resistant Mycobacterium tuberculosis complex (Mtbc) strains submitted to the National Tuberculosis Reference Laboratory (NTRL) in Maputo, Mozambique, between 2015 and 2021. Phylogenetic strain classification, genomic resistance prediction, and cluster analysis were performed.

FINDINGS

Between Jan 1, 2015, and July 31, 2021, 2606 Mtbc isolates with an isoniazid or rifampicin resistance were identified in the NTRL biobank, of which, 1483 (56·9%) were from men, 1114 (42·7%) from women, and nine (0·4%) were unknown. Genome-based drug-resistant prediction classified 704 Mtbc strains as rifampicin resistant. 628 (89%) of the 704 Mtbc strains were classified MDR; of those, 146 (23%) were pre-extensively drug resistant (pre-XDR; additional fluoroquinolone resistance), and 24 (4%) extensively drug resistant (XDR; combined fluoroquinolone and bedaquiline resistance). Overall, 61 (9%) of 704 strains revealed resistance to bedaquiline: five (7%) of 76 rifampicin resistant plus bedaquiline resistant, 32 (7%) of 458 MDR plus bedaquiline resistant, and 24 (100%) of 24 XDR. Prevalence of bedaquiline resistance increased from 3% in 2016 to 14% in 2021. The cluster rate (12 single-nucleotide polymorphism threshold) was 42% for rifampicin-resistant strains, 78% for MDR strains, 94% for pre-XDR strains, and 96% for XDR Mtbc strains. 31 (4%) of 704 Mtbc strains, belonging to a diagnostic escape outbreak strain previously described in Eswatini (group_56), had an rpoB Ile491Phe mutation which is not detected by Xpert MTB/RIF (no other rpoB mutation). Of these, 23 (74%) showed additional resistance to bedaquiline, 13 (42%) had bedaquiline and fluoroquinolone resistance, and two (6%) were bedaquiline, fluoroquinolone, and delamanid resistant.

INTERPRETATION

Pre-XDR resistance is highly prevalent among MDR Mtbc strains in Mozambique and so is bedaquiline resistance; and the frequency of bedaquiline resistance quadrupled over time and was found even in Mtbc strains without fluoroquinolone resistance. Importantly, strains with Ile491Phe mutation were frequent, accounting for 31% (n=10) of MDR plus bedaquiline-resistant strains and 54% (n=13) of XDR Mtbc strains. Given the current diagnostic algorithms and treatment regimens, both the emergence of rifampicin resistance due to Ile491Phe and bedaquiline resistance might jeopardise MDR tuberculosis prevention and care unless sequencing-based technology is rolled out. The potential cross border spread of diagnostic escape strains needs further investigation.

FUNDING

The German Ministry of Health through the Seq_MDRTB-Net project, the Deutsche Forschungsgemeinschaft under Germany's Excellence Strategy Precision Medicine in Inflammation and the Research Training Group 2501 TransEvo, the Leibniz Science Campus Evolutionary Medicine of the Lung, and the German Ministry of Education and Research via the German Center for Infection Research.

Prison as a driver of recent transmissions of multidrug-resistant tuberculosis in Callao, Peru: a cross-sectional study

Background

We sought to identify resistance patterns and key drivers of recent multidrug-resistant tuberculosis (MDR-TB) transmission in a TB-prevalent area in Peru.

Methods

Cross-sectional study including MDR Mycobacterium tuberculosis complex (Mtbc) strains identified in Callao-Peru between April 2017 and February 2019. Mtbc DNA was extracted for whole genome sequencing which was used for phylogenetic inference, clustering, and resistance mutation analyses. Clusters indicative of recent transmission were defined based on a strain-to-strain distance of ≤5 (D5) single nucleotide polymorphisms (SNPs). Epidemiologic factors linked to MDR-TB clustering were analyzed using Poisson regression.

Findings

171 unique MDR-Mtbc strains were included; 22 (13%) had additional fluoroquinolone resistance and were classified as pre-XDR. Six strains (3.5%) harboured bedaquiline (BDQ) resistance mutations and were classified as MDR + BDQ. 158 (92%) Mtbc strains belonged to lineage 4 and 13 (8%) to lineage 2. Using a cluster threshold of ≤5 SNPs, 98 (57%) strains were grouped in one of the 17 D5 clusters indicative of recent transmission, ranging in size from 2 to the largest cluster formed by 53 4.3.3 strains (group_1). Lineage 4.3.3 strains showed the overall highest cluster rate (43%). In multivariate analyses, current or previous imprisonment was independently associated with being part of any MDR-TB transmission clusters (adjusted prevalence ratio [aPR], 1.45; 95% CI, 1.09-1.92).

Interpretation

Pre-XDR-TB emerged in more than 10% of the MDR-TB strains investigated. Transmission of 4.3.3 Mtbc strains especially of the dominant group_1 clone is a major driver of the MDR-TB epidemic in Callao. Current or previous imprisonment was linked to recent MDR-TB transmissions, indicating an important role of prisons in driving the MDR-TB epidemic.

Funding

This work was supported in part by the ERANet-LAC Network of the European Union, Latin America and the Caribbean Countries on Joint Innovation and Research Activities, and FONDECYT. Additional support was received from Leibniz Science Campus Evolutionary Medicine of the Lung, the Deutsche Forschungsgemeinschaft (German Research Foundation, under Germany's Excellence Strategy-EXC 2167 Precision Medicine in Inflammation), and the Research Training Group 2501 TransEvo.

Phylogenetic relationships of Mycobacterium tuberculosis isolates in Poland: The emergence of Beijing genotype among multidrug-resistant cases

Introduction

The epidemiological situation of tuberculosis (TB) in Poland urges for its continuous and scrupulous monitoring. The objective of this study was to explore the genetic diversity of multidrug-resistant (MDR) and drug-susceptible (DS) Mycobacterium tuberculosis isolates from Poland with a combination of spoligotyping and high-resolution mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) analysis. The results were placed in the Northern and Eastern Europe context.

Methods

The study included 89 (39 MDR and 50 DS) M. tuberculosis isolates collected from as many patients between 2018 and 2021 in Poland. The analysis was done using spoligotyping, and MIRU-VNTR typing at 24 standard loci. The data were compared to those available on Poland and neighbors and global M. tuberculosis datasets.

Results

The main identified families were Beijing (28.1%) and Haarlem (16.8%) while 34.8% of isolates were in the heterogeneous L4-unclassified group. Although the Beijing family was the most prevalent (61.5%) among MDR-TB cases, it accounted for only 2% of DS isolates. Among foreign-born patients, a higher ratio of MDR isolates were observed when compared with those who Poland-born (64.3% vs. 40%). Furthermore, all patients from the Former Soviet Union (FSU) countries were infected with MDR-TB.

Discussion

Whereas DS M. tuberculosis population in Poland is dominated by L4 isolates, MDR isolates are mostly of the Beijing genotype. The rise in the prevalence of the Beijing isolates in Poland, coupled with high proportion of the Beijing genotype among foreign-born TB patients may reflect an ongoing transmission of this family, imported to Poland mainly from FSU countries.

Improving tuberculosis surveillance by detecting international transmission using publicly available whole genome sequencing data

IntroductionImproving the surveillance of tuberculosis (TB) is especially important for multidrug-resistant (MDR) and extensively drug-resistant (XDR) TB. The large amount of publicly available whole genome sequencing (WGS) data for TB gives us the chance to re-use data and to perform additional analyses at a large scale.AimWe assessed the usefulness of raw WGS data of global MDR/XDR Mycobacterium tuberculosis isolates available from public repositories to improve TB surveillance.MethodsWe extracted raw WGS data and the related metadata of M. tuberculosis isolates available from the Sequence Read Archive. We compared this public dataset with WGS data and metadata of 131 MDR- and XDR M. tuberculosis isolates from Germany in 2012 and 2013.ResultsWe aggregated a dataset that included 1,081 MDR and 250 XDR isolates among which we identified 133 molecular clusters. In 16 clusters, the isolates were from at least two different countries. For example, Cluster 2 included 56 MDR/XDR isolates from Moldova, Georgia and Germany. When comparing the WGS data from Germany with the public dataset, we found that 11 clusters contained at least one isolate from Germany and at least one isolate from another country. We could, therefore, connect TB cases despite missing epidemiological information.ConclusionWe demonstrated the added value of using WGS raw data from public repositories to contribute to TB surveillance. Comparing the German with the public dataset, we identified potential international transmission events. Thus, using this approach might support the interpretation of national surveillance results in an international context.

Molecular Epidemiology of Tuberculosis in British Columbia, Canada: A 10-Year Retrospective Study

Background

Understanding regional molecular epidemiology allows for the development of more efficient tuberculosis prevention strategies in low-incidence settings.

Methods

We analyzed 24-locus mycobacterial interspersed repetitive-unit–variable-number tandem repeat (MIRU-VNTR) genotyping for 2290 Mycobacterium tuberculosis clinical isolates collected in the province of British Columbia (BC), Canada, in 2005–2014. Laboratory data for each isolate were linked to case-level clinical and demographic data. These data were used to describe the molecular epidemiology of tuberculosis across the province.

Results

We detected >1500 distinct genotypes across the 4 major M. tuberculosis lineages, reflecting BC’s diverse population. Disease site and clustering rates varied across lineages, and MIRU-VNTR was used to group the 2290 isolates into 189 clusters (2–70 isolates per cluster), with an overall clustering rate of 42.4% and an estimated local transmission rate of 34.1%. Risk factors for clustering varied between Canadian-born and foreign-born individuals; the former had increased odds (odds ratio, 7.8; 95% confidence interval [CI], 6.2–9.6) of belonging to a genotypic cluster, although nearly one-quarter of clusters included both Canadian- and foreign-born persons. Large clusters (≥10 cases) occurred more frequently within the M. tuberculosis Euro-American lineage, and individual-level risk factors associated with belonging to a large cluster included being Canadian born (adjusted odds ratio, 3.3; 95% CI, 2.3–4.8), residing in a rural area (2.3; 1.2–4.5), and illicit drug use (2.0; 1.2–3.4).

Conclusions

Although tuberculosis in BC largely arises through reactivation of latent tuberculosis in foreign-born persons, locally transmitted infections occur in discrete populations with distinct disease and risk factor profiles, representing groups for targeted interventions.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Assessment of the use and need for an integrated molecular surveillance of tuberculosis: an online survey in Germany

Background

The implementation of an integrated molecular surveillance (IMS) of tuberculosis (TB) is of high priority for TB control. IMS is defined as the systematic inclusion of molecular typing results in the national TB surveillance system. Although not standardized, an IMS of TB is already implemented in several low TB incidence countries. Germany is in the process of implementing a nationwide IMS of TB. This requires close collaboration between national and local health authorities. We conducted an online survey to understand the current use of molecular typing results for TB surveillance among the local public health offices (PHO)s in Germany, and to collect their perception and expectations towards the implementation of a nationwide IMS of TB.

Methods

The online survey was developed using the software Voxco and included 31 questions. The survey was sent to all the 377 local PHOs in Germany in April 2017. Responses were collected until June 2017.

Results

A total of 174/377 (46.2%) local PHOs participated in our survey, and 88/377 (23.3%) used molecular typing results in their routine TB surveillance work. The PHOs used molecular typing results especially as support for epidemiological contact tracing (62/88, 70.4%). We found statistically significant differences between answers of PHOs that did not use molecular typing results (n = 86) vs. PHOs that did use molecular typing results (n = 88): the latter perceived the use of molecular typing results as more beneficial for their work compared to the former (65.9% vs. 34.9%, p < 0.05). Moreover, the PHOs using molecular typing results expect for the future more support and coordination from regional and national public health institutes, especially regarding the identification and analysis of molecular clusters.

Conclusions

Our study is a step forward in the broader goal of implementing an IMS of TB in Germany. The local PHOs currently using the molecular typing results highlighted their positive attitude towards the implementation of an IMS, but also their needs of more support. Similar assessments might serve as an example for other countries which are on the way to implement a nationwide IMS of TB.

Electronic supplementary material

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

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.

DNA markers for tuberculosis diagnosis

Tuberculosis (TB), caused by Mycobacterium tuberculosis complex (MTBC), is an infectious disease with more than 10.4 million cases and 1.7 million deaths reported worldwide in 2016. The classical methods for detection and differentiation of mycobacteria are: acid-fast microscopy (Ziehl-Neelsen staining), culture, and biochemical methods. However, the microbial phenotypic characterization is time-consuming and laborious. Thus, fast, easy, and sensitive nucleic acid amplification tests (NAATs) have been developed based on specific DNA markers, which are commercially available for TB diagnosis. Despite these developments, the disease remains uncontrollable. The identification and differentiation among MTBC members with the use of NAATs remains challenging due, among other factors, to the high degree of homology within the members and mutations, which hinders the identification of specific target sequences in the genome with potential impact in the diagnosis and treatment outcomes. In silico methods provide predictive identification of many new target genes/fragments/regions that can specifically be used to identify species/strains, which have not been fully explored. This review focused on DNA markers useful for MTBC detection, species identification and antibiotic resistance determination. The use of DNA targets with new technological approaches will help to develop NAATs applicable to all levels of the health system, mainly in low resource areas, which urgently need customized methods to their specific conditions.

MIRU-profiler: a rapid tool for determination of 24-loci MIRU-VNTR profiles from assembled genomes of Mycobacterium tuberculosis

Background

Tuberculosis (TB) resulted in an estimated 1.7 million deaths in the year 2016. The disease is caused by the members of Mycobacterium tuberculosis complex, which includes Mycobacterium tuberculosis, Mycobacterium bovis and other closely related TB causing organisms. In order to understand the epidemiological dynamics of TB, national TB control programs often conduct standardized genotyping at 24 Mycobacterial-Interspersed-Repetitive-Units (MIRU)-Variable-Number-of-Tandem-Repeats (VNTR) loci. With the advent of next generation sequencing technology, whole-genome sequencing (WGS) has been widely used for studying TB transmission. However, an open-source software that can connect WGS and MIRU-VNTR typing is currently unavailable, which hinders interlaboratory communication. In this manuscript, we introduce the MIRU-profiler program which could be used for prediction of MIRU-VNTR profile from WGS of M. tuberculosis.

Implementation

The MIRU-profiler is implemented in shell scripting language and depends on EMBOSS software. The in-silico workflow of MIRU-profiler is similar to those described in the laboratory manuals for genotyping M. tuberculosis. Given an input genome sequence, the MIRU-profiler computes alleles at the standard 24-loci based on in-silico PCR amplicon lengths. The final output is a tab-delimited text file detailing the 24-loci MIRU-VNTR pattern of the input sequence.

Validation

The MIRU-profiler was validated on four datasets: complete genomes from NCBI-GenBank (n = 11), complete genomes for locally isolated strains sequenced using PacBio (n = 4), complete genomes for BCG vaccine strains (n = 2) and draft genomes based on 250 bp paired-end Illumina reads (n = 106).

Results

The digital MIRU-VNTR results were identical to the experimental genotyping results for complete genomes of locally isolated strains, BCG vaccine strains and five out of 11 genomes from the NCBI-GenBank. For draft genomes based on short Illumina reads, 21 out of 24 loci were inferred with a high accuracy, while a number of inaccuracies were recorded for three specific loci (ETRA, QUB11b and QUB26). One of the unique features of the MIRU-profiler was its ability to process multiple genomes in a batch. This feature was tested on all complete M. tuberculosis genome (n = 157), for which results were successfully obtained in approximately 14 min.

Conclusion

The MIRU-profiler is a rapid tool for inference of digital MIRU-VNTR profile from the assembled genome sequences. The tool can accurately infer repeat numbers at the standard 24 or 21/24 MIRU-VNTR loci from the complete or draft genomes respectively. Thus, the tool is expected to bridge the communication gap between the laboratories using WGS and those using the conventional MIRU-VNTR typing.

Automated IS6110-based fingerprinting of Mycobacterium tuberculosis: Reaching unprecedented discriminatory power and versatility

BACKGROUND

Several technical hurdles and limitations have restricted the use of IS6110 restriction fragment length polymorphism (IS6110 RFLP), the most effective typing method for detecting recent tuberculosis (TB) transmission events. This has prompted us to conceive an alternative modality, IS6110-5'3'FP, a plasmid-based cloning approach coupled to a single PCR amplification of differentially labeled 5' and 3' IS6110 polymorphic ends and their automated fractionation on a capillary sequencer. The potential of IS6110-5'3'FP to be used as an alternative to IS6110 RFLP has been previously demonstrated, yet further technical improvements are still required for optimal discriminatory power and versatility.

OBJECTIVES

Here we introduced critical amendments to the original IS6110-5'3'FP protocol and compared its performance to that of 24-loci multiple interspersed repetitive unit-variable number tandem repeats (MIRU-VNTR), the current standard method for TB transmission analyses.

METHODS

IS6110-5'3'FP protocol modifications involved: (i) the generation of smaller-sized polymorphic fragments for efficient cloning and PCR amplification, (ii) omission of the plasmid amplification step in E. coli for shorter turnaround times, (iii) the use of more stable fluorophores for increased sensitivity, (iv) automated subtraction of background fluorescent signals, and (v) the automated conversion of fluorescent peaks into binary data.

RESULTS

In doing so, the overall turnaround time of IS6110-5'3'FP was reduced to 4 hours. The new protocol allowed detecting almost all 5' and 3' IS6110 polymorphic fragments of any given strain, including IS6110 high-copy number Beijing strains. IS6110-5'3'FP proved much more discriminative than 24-loci MIRU-VNTR, particularly with strains of the M. tuberculosis lineage 4.

CONCLUSIONS

The IS6110-5'3'FP protocol described herein reached the optimal discriminatory potential of IS6110 fingerprinting and proved more accurate than 24-loci MIRU-VNTR in estimating recent TB transmission. The method, which is highly cost-effective, was rendered versatile enough to prompt its evaluation as an automatized solution for a TB integrated molecular surveillance.

Comparative study of RFLP-IS6110 and MIRU-VNTR from Mycobacterium tuberculosis isolated in the state of Minas Gerais, Brazil

DNA genotyping of Mycobacterium tuberculosis has been widely applied in the understanding of disease transmission in many countries. The purpose of this study was to genotype the strains of M. tuberculosis isolated in patients with new tuberculosis (TB) cases in Minas Gerais, as well as to compare the similarity, discriminatory power, and agreement of the clusters between the IS6110 Restriction Fragment Length Polymorfism (RFLP) and 12 loci Variable Number Tandem Repeat – Mycobacterial Interspersed Repetitive Units (MIRU-VNTR) techniques. It was observed that 32% (66/204) of the isolated strains in the RFLP-IS6110 and 50.9% (104/204) of the isolated strains in the MIRU-VNTR presented a similarity of equal to or above 85%. The RFLP-IS6110 and MIRU-VNTR proved to contain a high discriminatory power. The similarity index resulting from the RFLP showed no recent transmission. Good agreement was observed between the techniques when clusters were detected; however, the best epidemiological relationship was found when using the RFLP-IS6110.

Genetic diversity & drug sensitivity profiles of Mycobacterium tuberculosis isolates from two slums of Jaipur city, Rajasthan, India

Background & objectives:

Slums are considered as hotspots of tuberculosis (TB). The study of genetic diversity and drug susceptibility profile of Mycobacterium tuberculosis (MTB) will help understand the transmission dynamics and can be used for better prevention and control of the disease. The aim of this study was to determine the drug susceptibility profiles and genetic diversity using the random amplified polymorphic DNA (RAPD) and mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU VNTR) of MTB isolates from sputum samples of pulmonary TB patients residing in the two slums of Jaipur city in Rajasthan, India.

Methods:

Sputum samples collected from pulmonary TB patients, their contacts and suspects during 2010-2012 were processed for microscopy and mycobacterial culture. Drug susceptibility testing was done by one per cent indirect proportion method on Lowenstein–Jensen medium for first-line anti-TB drugs rifampicin, isoniazid, ethambutol and streptomycin. MTB DNA was extracted by physicochemical method, and DNA fingerprinting was done by RAPD and MIRU VNTR analysis.

Results:

Among 175 sputum samples collected, 75 were positive (43.8%) for acid-fast bacilli, 83 for MTB culture and four were contaminated. Fifty two isolates (62.7%) were fully sensitive to four drugs, and five (6%) were multidrug resistant (MDR). RAPD analysis of 81 isolates revealed six clusters containing 23 (28.4%) isolates, and 58 (71.6%) were unique. MIRU VNTR analysis clustered 20 (24.7%) isolates, and 61 (75.3%) were unique.

Interpretation & conclusions:

About 62.7 per cent isolates from the sputum samples from slum areas were sensitive to four drugs; six per cent of isolates were MDR. Poly-resistance other than MDR was high (16%). About one-fourth isolates were clustered by either method. RAPD was rapid, less expensive but had low reproducibility. MIRU VNTR analysis could identify to greater extent the epidemiological link in the population studied.

Spoligotype defined lineages of Mycobacterium tuberculosis and drug resistance: Merely a casual correlation?

Drug-resistant tuberculosis (TB) is a major challenge to TB control strategy worldwide. Analysis of genetic polymorphism among drug resistant Mycobacterium tuberculosis (MTB) strains may help provide some insight into the transmission dynamics of these strains. Spoligotyping is a widely used technique to identify genetic polymorphism, based on 43 known spacers interspersed between direct repeat regions. Considerable work has been done in various parts of the world using this technique to identify and analyse the polymorphic nature of MTB. Many studies have been carried out to determine the association of drug resistance with spoligotype defined lineages, and much data has been produced over the years. New information continues to be generated. This review aims to put together the findings of relevant studies in an attempt to understand the correlation of drug resistance with spoligotype defined lineages of MTB. This would help provide a perspective of the available data that can be used as a starting point to understand the molecular epidemiology of drug resistant TB.

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.

Molecular characterization of Mycobacterium tuberculosis isolates from Tehran, Iran by restriction fragment length polymorphism analysis and spoligotyping

INTRODUCTION Characterization of Mycobacterium tuberculosis (MTB) isolates by DNA fingerprinting has contributed to tuberculosis (TB) control. The aim of this study was to determine the genetic diversity of MTB isolates from Tehran province in Iran. METHODS MTB isolates from 60 Iranian and 10 Afghan TB patients were fingerprinted by standard IS6110-restriction fragment length polymorphism (RFLP) analysis and spoligotyping. RESULTS The copy number of IS6110 ranged from 10-24 per isolate. The isolates were classified into 22 clusters showing ≥ 80% similarity by RFLP analysis. Fourteen multidrug-resistant (MDR) isolates were grouped into 4 IS6110-RFLP clusters, with 10 isolates [71% (95% CI: 45-89%)] in 1 cluster, suggesting a possible epidemiological linkage. Eighteen Iranian isolates showed ≥ 80% similarity with Afghan isolates. There were no strains with identical fingerprints. Spoligotyping of 70 isolates produced 23 distinct patterns. Sixty (85.7%) isolates were grouped into 13 clusters, while the remaining 10 isolates (14.2%) were not clustered. Ural (formerly Haarlem4) (n = 22, 31.4%) was the most common family followed by Central Asian strain (CAS) (n = 18, 25.7%) and T (n = 9, 12.8%) families. Only 1strain was characterized as having the Beijing genotype. Among 60 Iranian and 10 Afghan MTB isolates, 25% (95% CI: 16-37) and 70% (95% CI: 39-89) were categorized as Ural lineage, respectively. CONCLUSIONS A higher prevalence of Ural family MTB isolates among Afghan patients than among Iranian patients suggests the possible transmission of this lineage following the immigration of Afghans to Iran.

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.

Genomic Diversity of Mycobacterium tuberculosis Complex Strains in Cantabria (Spain), a Moderate TB Incidence Setting

Background

Tuberculosis (TB) control strategies are focused mainly on prevention, early diagnosis, compliance to treatment and contact tracing. The objectives of this study were to explore the frequency and risk factors of recent transmission of clinical isolates of Mycobacterium tuberculosis complex (MTBC) in Cantabria in Northern Spain from 2012 through 2013 and to analyze their clonal complexity for better understanding of the transmission dynamics in a moderate TB incidence setting.

Methods

DNA from 85 out of 87 isolates from bacteriologically confirmed cases of MTBC infection were extracted directly from frozen stocks and genotyped using the mycobacterial interspersed repetitive units-variable number tandem repeat (MIRU-VNTR) method. The MIRU-VNTRplus database tool was used to identify clusters and lineages and to build a neighbor joining (NJ) phylogenetic tree. In addition, data were compared to the SITVIT2 database at the Pasteur Institute of Guadeloupe.

Results

The rate of recent transmission was calculated to 24%. Clustering was associated with being Spanish-born. A high prevalence of isolates of the Euro-American lineage was found. In addition, MIRU-VNTR profiles of the studied isolates corresponded to previously found MIRU-VNTR types in other countries, including Spain, Belgium, Great Britain, USA, Croatia, South Africa and The Netherlands. Six of the strains analyzed represented clonal variants.

Conclusion

Transmission of MTBC is well controlled in Cantabria. The majority of TB patients were born in Spain. The population structure of MTBC in Cantabria has a low diversity of major clonal lineages with the Euro-American lineage predominating.

High prevalence of Mycobacterium tuberculosis mixed infection in the capital of moderate tuberculosis incidence country

OBJECTIVE

Recent studies using molecular epidemiological techniques have demonstrated mixed infection with multiple strains of Mycobacterium tuberculosis especially in countries with high tuberculosis (TB) burden. We aimed to determine the prevalence of mixed infection among patients with TB in the capital of Iran as a country with moderate incidence rate.

METHODS

Samples were collected randomly from January 2011 to December 2013 in Tehran, capital of Iran. A total of 75 M. tuberculosis isolates were genotyped by 24 loci mycobacterial interspersed repetitive unit-variable number tandem repeat typing (MIRU-VNTR) for screening the mixed infection.

RESULTS

Twenty patients (20/75) were identified with mixed infection, and the estimated rate of mixed infection was 26.6%. Thirteen out of the 24 loci were able to detect the mixed infection in our study.

CONCLUSIONS

Mixed infections occur at high prevalence among studied Iranian TB patients. Further research is inevitable to evaluate the association of mixed infection and disease progression and treatment.

SpoTyping: fast and accurate in silico Mycobacterium spoligotyping from sequence reads

SpoTyping is a fast and accurate program for in silico spoligotyping of Mycobacterium tuberculosis isolates from next-generation sequencing reads. This novel method achieves high accuracy for reads of both uniform and varying lengths, and is about 20 to 40 times faster than SpolPred. SpoTyping also integrates the function of producing a report summarizing associated epidemiological data from a global database of all isolates having the same spoligotype. SpoTyping is freely available at: https://github.com/xiaeryu/SpoTyping-v2.0.

Genetic Diversity of Mycobacterium tuberculosis Isolates from Assam, India: Dominance of Beijing Family and Discovery of Two New Clades Related to CAS1_Delhi and EAI Family Based on Spoligotyping and MIRU-VNTR Typing

Tuberculosis (TB) is one of the major public health concerns in Assam, a remote state located in the northeastern (NE) region of India. The present study was undertaken to explore the circulating genotypes of Mycobacterium tuberculosis complex (MTBC) in this region. A total of 189 MTBC strains were collected from smear positive pulmonary tuberculosis cases from different designated microscopy centres (DMC) from various localities of Assam. All MTBC isolates were cultured on Lowenstein-Jensen (LJ) media and subsequently genotyped using spoligotyping and 24-loci mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) typing. Spoligotyping of MTBC isolates revealed 89 distinct spoligo patterns. The most dominant MTBC strain belonged to Beijing lineage and was represented by 35.45% (n = 67) of total isolates, followed by MTBC strains belonging to Central Asian-Delhi (CAS/Delhi) lineage and East African Indian (EAI5) lineage. In addition, in the present study 43 unknown spoligo patterns were detected. The discriminatory power of spoligotyping was found to be 0.8637 based on Hunter Gaston Discriminatory Index (HGDI). On the other hand, 24-loci MIRU-VNTR typing revealed that out of total 189 MTBC isolates from Assam 185 (97.9%) isolates had unique MIRU-VNTR profiles and 4 isolates grouped into 2 clusters. Phylogenetic analysis of 67 Beijing isolates based on 24-loci MIRU-VNTR typing revealed that Beijing isolates from Assam represent two major groups, each comprising of several subgroups. Neighbour-Joining (NJ) phylogenetic tree analysis based on combined spoligotyping and 24-loci MIRU-VNTR data of 78 Non-Beijing isolates was carried out for strain lineage identification as implemented by MIRU-VNTRplus database. The important lineages of MTBC identified were CAS/CAS1_Delhi (41.02%, n = 78) and East-African-Indian (EAI, 33.33%). Interestingly, phylogenetic analysis of orphan (23.28%) MTBC spoligotypes revealed that majority of these orphan isolates from Assam represent two new sub-clades Assam/EAI and Assam/CAS. The prevalence of multidrug resistance (MDR) in Beijing and Non-Beijing strains was found to be 10.44% and 9.01% respectively. In conclusion, the present study has shown the predominance of Beijing isolates in Assam which is a matter of great concern because Beijing strains are considered to be ecologically more fit enabling wider dissemination of M. tuberculosis. Other interesting finding of the present study is the discovery of two new clades of MTBC isolates circulating in Assam. More elaborate longitudinal studies are required to be undertaken in this region to understand the transmission dynamics of MTBC.

Genetic diversity and molecular epidemiology of multidrug-resistant Mycobacterium tuberculosis in Minas Gerais State, Brazil

Background

We aimed to characterize the genetic diversity of drug-resistant Mycobacterium tuberculosis (MTb) clinical isolates and investigate the molecular epidemiology of multidrug-resistant (MDR) tuberculosis from Minas Gerais State, Brazil.

Methods

One hundred and four MTb clinical isolates were assessed by IS6110-RFLP, 24-locus mycobacterial interspersed repetitive units variable-number tandem repeats (MIRU-VNTR), TB-SPRINT (simultaneous spoligotyping and rifampicin-isoniazid drug-resistance mutation analysis) and 3R-SNP-typing (analysis of single-nucleotide polymorphisms in the genes involved in replication, recombination and repair functions).

Results

Fifty-seven different IS6110-RFLP patterns were found, among which 50 had unique patterns and 17 were grouped into seven clusters. The discriminatory index (Hunter and Gaston, HGDI) for RFLP was 0.9937. Ninety-nine different MIRU-VNTR patterns were found, 95 of which had unique patterns and nine isolates were grouped into four clusters. The major allelic diversity index in the MIRU-VNTR loci ranged from 0.6568 to 0.7789. The global HGDI for MIRU-VNTR was 0.9991. Thirty-two different spoligotyping profiles were found: 16 unique patterns (n = 16) and 16 clustered profiles (n = 88). The HGDI for spoligotyping was 0.9009. The spoligotyped clinical isolates were phylogenetically classified into Latin-American Mediterranean (66.34 %), T (14.42 %), Haarlem (5.76 %), X (1.92 %), S (1.92 %) and U (unknown profile; 8.65 %). Among the U isolates, 77.8 % were classified further by 3R-SNP-typing as 44.5 % Haarlem and 33.3 % LAM, while the 22.2 % remaining were not classified. Among the 104 clinical isolates, 86 were identified by TB-SPRINT as MDR, 12 were resistant to rifampicin only, one was resistant to isoniazid only, three were susceptible to both drugs, and two were not successfully amplified by PCR. A total of 42, 28 and eight isolates had mutations in rpoB positions 531, 526 and 516, respectively. Correlating the cluster analysis with the patient data did not suggest recent transmission of MDR-TB.

Conclusions

Although our results do not suggest strong transmission of MDR-TB in Minas Gerais (using a classical 100 % MDR-TB identical isolates cluster definition), use of a smoother cluster definition (>85 % similarity) does not allow us to fully eliminate this possibility; hence, around 20–30 % of the isolates we analyzed might be MDR-TB transmission cases.

Mycobacterium tuberculosis Genotypes Determined by Spoligotyping to Be Circulating in Colombia between 1999 and 2012 and Their Possible Associations with Transmission and Susceptibility to First-Line Drugs

INTRODUCTION

Tuberculosis (TB) remains a primary public health problem worldwide. The number of multidrug-resistant tuberculosis (MDR TB) cases has increased in recent years in Colombia. Knowledge of M. tuberculosis genotypes defined by spoligotyping can help determine the circulation of genotypes that must be controlled to prevent the spread of TB.

OBJECTIVE

To describe the genotypes of M. tuberculosis using spoligotyping in resistant and drug-sensitive isolates and their possible associations with susceptibility to first-line drugs.

METHODS

An analytical observational study was conducted that included 741 isolates of M. tuberculosis from patients. The isolates originated from 31 departments and were obtained by systematic surveillance between 1999 and 2012.

RESULTS

In total 61.94% of the isolates were resistant to 1 or more drugs, and 147 isolates were MDR. In total, 170 genotypes were found in the population structure of Colombian M. tuberculosis isolates. The isolates were mainly represented by four families: LAM (39.9%), Haarlem (19%), Orphan (17%) and T (9%). The SIT42 (LAM 9) was the most common genotype and contained 24.7% of the isolates, followed by the genotypes SIT62 (Haarlem1), SIT53 (T1), and SIT50 (H3). A high clustering of isolates was evident with 79.8% of the isolates classified into 32 groups. The Beijing family was associated with resistant isolates, whereas the Haarlem and T families were associated with sensitive isolates. The Haarlem family was also associated with grouped isolates (p = 0.031).

CONCLUSIONS

A high proportion (approximately 80%) of isolates was found in clusters; these clusters were not associated with resistance to first-line drugs. The Beijing family was associated with drug resistance, whereas the T and Haarlem families were associated with susceptibility in the Colombian isolates studied.

Mycobacterium tuberculosis phylogeography in the context of human migration and pathogen's pathobiology: Insights from Beijing and Ural families

Here, I review the population structure and phylogeography of the two contrasting families of Mycobacterium tuberculosis, Beijing and Ural, in the context of strain pathobiology and human history and migration. Proprietary database (12-loci MIRU-VNTR profiles of 3067 Beijing genotype isolates) was subjected to phylogenetic and statistical analysis. The highest rate (90%) and diversity (HGI 0.80-0.95) of the Beijing genotype in North China suggest it to be its area of origin. Under VNTR-based MDS analysis the interpopulation genetic distances correlated with geography over uninterrupted landmasses. In contrast, large water distances together with long time generated remarkable outliers. Weak and less expected affinities of the distant M. tuberculosis populations may reflect hidden epidemiological links due to unknown migration. Association with drug-resistance or increased virulence/transmissibility along with particular human migration flows shape global dissemination of some Beijing clones. The paucity of data on the Ural genotype prevents from high-resolution analysis that was mainly based on the available spoligotyping data. The North/East Pontic area marked with the highest prevalence of the Ural family may have been the area of its origin and primary dispersal in Eurasia. Ural strains are not marked by increased pathogenic capacities, increased transmissibility and association with drug resistance (but most recent reports describe an alarming increase of MDR Ural strains in some parts of eastern Europe and northwestern Russia). Large-scale SNP or WGS population-based studies targeting strains from indigenous populations and, eventually, analysis of ancient DNA will better test these hypotheses. Host genetics factors likely play the most prominent role in differential dissemination of particular M. tuberculosis genotypes.

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.

Prevalence and occurrence rate of Mycobacterium tuberculosis Haarlem family multi-drug resistant in the worldwide population: A systematic review and meta-analysis

BACKGROUND

Transmission of Mycobacterium tuberculosis (M. tuberculosis) can occur in different ways. Furthermore, drug resistant in M. tuberculosis family is a major problem that creates obstacles in treatment and control of tuberculosis (TB) in the world. One of the most prevalent families of M. tuberculosis is Haarlem, and it is associated with drug resistant. Our objectives of this study were to determine the prevalence and occurrence rate of M. tuberculosis Haarlem family multi-drug resistant (MDR) in the worldwide using meta-analysis based on a systematic review that performed on published articles.

MATERIALS AND METHODS

Data sources of this study were 78 original articles (2002-2012) that were published in the literatures in several databases including PubMed, Science Direct, Google Scholar, Biological abstracts, ISI web of knowledge and IranMedex. The articles were systematically reviewed for prevalence and rate of MDR. Data were analyzed using meta-analysis and random effects models with the software package Meta R, Version 2.13 (P < 0.10).

RESULTS

Final analysis included 28601 persons in 78 articles. The highest and lowest occurrence rate of Haarlem family in M. tuberculosis was in Hungary in 2006 (66.20%) with negative MDR-TB and in China in 2010 (0.8%), respectively. From 2002 to 2012, the lowest rate of prevalence was in 2010, and the highest prevalence rate was in 2012. Also 1.076% were positive for MDR and 9.22% were negative (confidence interval: 95%).0020.

CONCLUSION

Many articles and studies are performed in this field globally, and we only chose some of them. Further studies are needed to be done in this field. Our study showed that M. tuberculosis Haarlem family is prevalent in European countries. According to the presence of MDR that was seen in our results, effective control programs are needed to control the spread of drug-resistant strains, especially Haarlem family.

Consequences of genomic diversity in Mycobacterium tuberculosis

The causative agent of human tuberculosis, Mycobacterium tuberculosis complex (MTBC), comprises seven phylogenetically distinct lineages associated with different geographical regions. Here we review the latest findings on the nature and amount of genomic diversity within and between MTBC lineages. We then review recent evidence for the effect of this genomic diversity on mycobacterial phenotypes measured experimentally and in clinical settings. We conclude that overall, the most geographically widespread Lineage 2 (includes Beijing) and Lineage 4 (also known as Euro-American) are more virulent than other lineages that are more geographically restricted. This increased virulence is associated with delayed or reduced pro-inflammatory host immune responses, greater severity of disease, and enhanced transmission. Future work should focus on the interaction between MTBC and human genetic diversity, as well as on the environmental factors that modulate these interactions.

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.

Prospective genotyping of Mycobacterium tuberculosis from fresh clinical samples

Shorter time-to-result is key for improving molecular-guided epidemiological investigation of tuberculosis (TB) cases. We performed a prospective study to evaluate the use of standardized MIRU-VNTR (mycobacterial interspersed repetitive-unit-variable-number tandem-repeat) typing of Mycobacterium tuberculosis directly on 79 fresh clinical samples from 26 TB patients consecutively enrolled over a 17-month period. Overall, complete 24-locus types were obtained for 18 out of the 26 (69.2%) patients and 14 of the 16 grade 3+ and grade 2+ samples (87.5%). The degree of completion of the genotypes obtained significantly correlated with smear microscopy grade both for 26 first samples (p = 0.0003) and for 53 follow-up samples (p = 0.002). For 20 of the 26 patients for whom complete or even incomplete M. tuberculosis isolate genotypes were obtained, typing applied to the clinical samples allowed the same unambiguous conclusions regarding case clustering or uniqueness as those that could have been drawn based on the corresponding cultured isolates. Standard 24 locus MIRU-VNTR typing of M. tuberculosis can be applied directly to fresh clinical samples, with typeability depending on the bacterial load in the sample.

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.

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.

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.

Assessment of Mycobacterium tuberculosis transmission in Oxfordshire, UK, 2007-12, with whole pathogen genome sequences: an observational study

Background

Patients born outside the UK have contributed to a 20% rise in the UK’s tuberculosis incidence since 2000, but their effect on domestic transmission is not known. Here we use whole-genome sequencing to investigate the epidemiology of tuberculosis transmission in an unselected population over 6 years.

Methods

We identified all residents with Oxfordshire postcodes with a Mycobacterium tuberculosis culture or a clinical diagnosis of tuberculosis between Jan 1, 2007, and Dec 31, 2012, using local databases and checking against the national Enhanced Tuberculosis Surveillance database. We used Illumina technology to sequence all available M tuberculosis cultures from identified cases. Sequences were clustered by genetic relatedness and compared retrospectively with contact investigations. The first patient diagnosed in each cluster was defined as the index case, with links to subsequent cases assigned first by use of any epidemiological linkage, then by genetic distance, and then by timing of diagnosis.

Findings

Although we identified 384 patients with a diagnosis of tuberculosis, country of birth was known for 380 and we sequenced isolates from 247 of 269 cases with culture-confirmed disease. 39 cases were genomically linked within 13 clusters, implying 26 local transmission events. Only 11 of 26 possible transmissions had been previously identified through contact tracing. Of seven genomically confirmed household clusters, five contained additional genomic links to epidemiologically unidentified non-household members. 255 (67%) patients were born in a country with high tuberculosis incidence, conferring a local incidence of 109 cases per 100 000 population per year in Oxfordshire, compared with 3·5 cases per 100 000 per year for those born in low-incidence countries. However, patients born in the low-incidence countries, predominantly UK, were more likely to have pulmonary disease (adjusted odds ratio 1·8 [95% CI 1·2–2·9]; p=0·009), social risk factors (4·4 [2·0–9·4]; p<0·0001), and be part of a local transmission cluster (4·8 [1·6–14·8]; p=0·006).

Interpretation

Although inward migration has contributed to the overall tuberculosis incidence, our findings suggest that most patients born in high-incidence countries reactivate latent infection acquired abroad and are not involved in local onward transmission. Systematic screening of new entrants could further improve tuberculosis control, but it is important that health care remains accessible to all individuals, especially high-risk groups, if tuberculosis control is not to be jeopardised.

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.

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.

Whole genome sequencing versus traditional genotyping for investigation of a Mycobacterium tuberculosis outbreak: a longitudinal molecular epidemiological study

BACKGROUND

Understanding Mycobacterium tuberculosis (Mtb) transmission is essential to guide efficient tuberculosis control strategies. Traditional strain typing lacks sufficient discriminatory power to resolve large outbreaks. Here, we tested the potential of using next generation genome sequencing for identification of outbreak-related transmission chains.

METHODS AND FINDINGS

During long-term (1997 to 2010) prospective population-based molecular epidemiological surveillance comprising a total of 2,301 patients, we identified a large outbreak caused by an Mtb strain of the Haarlem lineage. The main performance outcome measure of whole genome sequencing (WGS) analyses was the degree of correlation of the WGS analyses with contact tracing data and the spatio-temporal distribution of the outbreak cases. WGS analyses of the 86 isolates revealed 85 single nucleotide polymorphisms (SNPs), subdividing the outbreak into seven genome clusters (two to 24 isolates each), plus 36 unique SNP profiles. WGS results showed that the first outbreak isolates detected in 1997 were falsely clustered by classical genotyping. In 1998, one clone (termed "Hamburg clone") started expanding, apparently independently from differences in the social environment of early cases. Genome-based clustering patterns were in better accordance with contact tracing data and the geographical distribution of the cases than clustering patterns based on classical genotyping. A maximum of three SNPs were identified in eight confirmed human-to-human transmission chains, involving 31 patients. We estimated the Mtb genome evolutionary rate at 0.4 mutations per genome per year. This rate suggests that Mtb grows in its natural host with a doubling time of approximately 22 h (400 generations per year). Based on the genome variation discovered, emergence of the Hamburg clone was dated back to a period between 1993 and 1997, hence shortly before the discovery of the outbreak through epidemiological surveillance.

CONCLUSIONS

Our findings suggest that WGS is superior to conventional genotyping for Mtb pathogen tracing and investigating micro-epidemics. WGS provides a measure of Mtb genome evolution over time in its natural host context.

Molecular typing of Mycobacterium tuberculosis isolates circulating in Henan, central China

The aim of this study was to characterize the genotypes of Mycobacterium tuberculosis (MT) in isolates obtained from Henan, China, and to study the distribution features of Beijing strains in Henan. A total of 443 MT strains isolated in Henan Province were tested for susceptibility to isonicotinylhydrazide (INH), rifampicin (RFP), ethambutol (EMB) and streptomycin (SM), and genotyped by spoligotyping. The clustering of genotypes revealed 4 gene clusters (Beijing and Beijing-like, T, Manu, and S and LAM3) and 24 genotypes. In total, 387 (87.4%) of the strains were Beijing strains. The frequency of multidrug-resistant (MDR) strains was significantly higher in the Beijing and Beijing-like strains than in the other strains (χ(2)=4.6564, P<0.05). However, the percentages of drug resistance and sensitivity in the Beijing strains were almost the same as those in the non-Beijing strains. The proportion of Beijing strains in the ≤60-year-old TB patients was significantly higher than in the >60-year-old TB patients (χ(2)=32.053, P<0.001). The distribution of Beijing strains deceased gradually from the east to the west in Henan Province (P<0.0001). The data demonstrate that the Beijing genotype is a major type in the area and may be related to enhanced transmissibility. The transmission of Beijing family strains has increased in Henan and its incidence is highest in the east of Henan. The MDR strains were significant in the spread of MT.